Differentially expressed genes of the carpet shell clam Ruditapes decussatus against Perkinsus olseni

A saposin domain-containing protein of tongue sole Cynoglossus semilaevis: Antimicrobial activity and mechanism

Prosaposin is a precursor that can be processed into four different saposins, designated as A, B, C, and D, which have multiple functions in mammals, including neuroprotection and immune modulation. The immune function of saposin in teleost remains largely unknown. In the present study, a saposin (SAP) domain-containing protein was identified in half-smooth tongue sole Cynoglossus semilaevis and named CsSDP. CsSDP harbors one SAP A domain and two SAP B domains. When expressed in HEK293T cells, CsSDP was specifically localized in the lysosome. When overexpressed in Escherichia coli, CsSDP markedly inhibited bacterial growth, and the inhibitory effect depended on two specific regions in the SAP A and SAP B domains. Two polypeptides (P32 and P30) derived from the above SAP A and B domains could bind to and inhibit the growth of both Gram-negative and Gram-positive bacteria. The ultrastructural analysis revealed that P32 and P30 killed target bacteria by disrupting the bacterial cell wall and inducing substantial release of cytoplasmic contents. These results shed new lights on the immune function of saposin domain-containing protein in teleost.

Search for new biomarkers of tolerance to Perkinsus olseni parasite infection in Ruditapes decussatus clams

The grooved carpet shell (Ruditapes decussatus) is a clam species with high economic and social importance in several European and Mediterranean countries. Production of this species suffered a decline caused by biotic (parasite infection) and abiotic factors (environmental factors, stress, poor management methods and intensive culture of the introduced species Ruditapes philippinarum). The protozoan parasite Perkinsus olseni is also responsible for the decline of production, being nowadays one of the major issues for clam culture. Molecular biomarkers that might represent tolerance of R. decussatus to P. olseni have already been uncovered, shedding light in a possible production improvement by selecting those clams with a strongest immune response. In the present study, new tolerance biomarkers to P. olseni infection in R. decussatus were identified. The haemolymph proteomic profiles of naturally non/low-infected (tolerant) and highly-infected (susceptible) clams by the parasite across several heavy affected areas of Europe were characterized through a shotgun proteomics approach. Also, the mechanisms that might be involved in the responses against the disease in chronic infections were explored. Proteins related to energy restoration and balance, metabolic regulation, energy accumulation, ROS production, lysosomal activity, amino acid synthesis, proteolytic activity, iron regulation, iron withholding, and immune response modulation were significantly regulated in susceptible clams. In the tolerant group, proteins related to phagocytosis regulation, control of cell growth and proliferation, gonadal maturation, regulation of apoptosis, growth modulation, response to oxidative stress, iron regulation, shell development and metabolic regulation were significantly expressed. In summary, the protein expression profile of tolerant individuals suggests that an efficient pathogen elimination mechanism coupled to a better metabolic regulation leads to a tolerance to the parasite infection by limiting the spread through the tissues.

MicroRNA-mediated stress response in bivalve species

Bivalve mollusks are important aquatic organisms, which are used for biological monitoring because of their abundance, ubiquitous nature, and abilities to adapt to different environments. MicroRNAs (miRNAs) are small noncoding RNAs, which typically silence the expression of target genes; however, certain miRNAs directly or indirectly upregulate their target genes. They are rapidly modulated and play an essential role in shaping the response of organisms to stresses. Based on the regulatory function and rapid alteration of miRNAs, they could act as biomarkers for biotic and abiotic stress, including environmental stresses and contaminations. Moreover, mollusk, particularly hemocytes, rapidly respond to environmental changes, such as pollution, salinity changes, and desiccation, which makes them an attractive model for this purpose. Thus, bivalve mollusks could be considered a good animal model to examine a system's response to different environmental conditions and stressors. miRNAs have been reported to adjust the adaptation and physiological functions of bivalves during endogenous and environmental stressors. In this review, we aimed to discuss the potential mechanisms underlying the response of bivalves to stressors and how miRNAs orchestrate this process; however, if necessary, other organisms' response is included to explain specific processes.

Transcriptional Alteration of Gene Biomarkers in Hemocytes of Wild Ostrea edulis with Molecular Evidence of Infections with Bonamia spp. and/or Marteilia refringens Parasites

The European flat Ostrea edulis is highly susceptible to intracellular parasitic infections, particularly bonamiosis and marteiliosis. The defensive response of oyster to both bonamiosis and marteiliosis is typically mediated by hemocytes, which play a pivotal role in immune system homeostasis. In the present study, we first used a DNA-based tool in order to rapidly and specifically detect the presence of parasites in oysters from natural banks in the middle Adriatic Sea. In a second step, we used qRT-PCR to analyze the mRNA levels of a set of genes (i.e., superoxide dismutase (SOD), glutathione S-transferase (GST), metallothionein (MT), heat shock protein (HSP) 70 and 90, inhibitor of apoptosis (IAP), fas ligand (FAS), galectin (GAL) and extracellular superoxide dismutase (Ec-SOD)) expressed by hemocytes of flat oysters infected by the parasites, present singularly or in combination, compared to hemocytes from non-infected specimens. The results indicate that the presence of parasite DNA may be associated to a general upregulation of host genes related to apoptosis, detoxification and oxidative stress protection, with the exception of Ec-SOD, whose trend to a downregulation might reflect a mechanism for parasite escape before internalization.

New insights into the Manila clam - Perkinsus olseni interaction based on gene expression analysis of clam hemocytes and parasite trophozoites through in vitro challenges

The Manila clam (Ruditapes philippinarum) is the bivalve species with the highest global production from both fisheries and aquaculture, but its production is seriously threatened by perkinsosis, a disease caused by the protozoan parasite Perkinsus olseni. To understand the molecular mechanisms underlying R. philippinarum-P. olseni interactions, we analysed the gene expression profiles of in vitro challenged clam hemocytes and P. olseni trophozoites, using two oligo-microarray platforms, one previously validated for R. philippinarum hemocytes and a new one developed and validated in this study for P. olseni. Manila clam hemocytes were in vitro challenged with trophozoites, zoospores, and extracellular products from P. olseni in vitro cultures, while P. olseni trophozoites were in vitro challenged with Manila clam plasma along the same time-series (1 h, 8 h, and 24 h). The hemocytes showed a fast activation of the innate immune response, particularly associated with hemocyte recruitment, in the three types of challenges. Nevertheless, different immune-related pathways were activated in response to the different parasite stages, suggesting specific recognition mechanisms. Furthermore, the analyses provided useful complementary data to previous in vivo challenges, and confirmed the potential of some proposed biomarkers. The combined analysis of gene expression in host and parasite identified several processes in both the clam and P. olseni, such as redox and glucose metabolism, protease activity, apoptosis and iron metabolism, whose modulation suggests cross-talk between parasite and host. This information might be critical to determine the outcome of the infection, thus highlighting potential therapeutic targets. Altogether, the results of this study aid understanding the response and interaction between R. philippinarum and P. olseni, and will contribute to developing effective control strategies for this threatening parasitosis.

Proteomic profile of Ostrea edulis haemolymph in response to bonamiosis and identification of candidate proteins as resistance markers

European flat oyster Ostrea edulis populations have suffered extensive mortalities caused by bonamiosis. The protozoan parasite Bonamia ostreae is largely responsible for this disease in Europe, while its congener B. exitiosa has been detected more recently in various European countries. Both of these intracellular parasites are able to survive and proliferate within haemocytes, the main cellular effectors of the immune system in molluscs. Two-dimensional electrophoresis was used to compare the haemolymph protein profile between Bonamia spp.-infected and non-infected oysters within 3 different stocks, a Galician stock of oysters selected for resistance against bonamiosis, a non-selected Galician stock and a selected Irish stock. Thirty-four proteins with a presumably relevant role in the oyster-Bonamia spp. interaction were identified; they were involved in major metabolic pathways, such as energy production, respiratory chain, oxidative stress, signal transduction, transcription, translation, protein degradation and cell defence. Furthermore, the haemolymph proteomic profiles of the non-infected oysters of the 2 Galician stocks were compared. As a result, 7 proteins representative of the non-infected Galician oysters selected for resistance against bonamiosis were identified; these 7 proteins could be considered as candidate markers of resistance to bonamiosis, which should be further assessed.

Gene expression analysis of Ruditapes philippinarum haemocytes after experimental Perkinsus olseni zoospore challenge and infection in the wild

The production of Manila clam (Ruditapes philippinarum) is seriously threatened by the protistan parasite Perkinsus olseni. We characterized and compared gene expression of Manila clam haemocytes in response to P. olseni in a time-course (10 h, 24 h, 8 d) controlled laboratory challenge (LC), representing the first step of infection, and in a more complex infection in the wild (WI), using a validated oligo-microarray containing 11,232 transcripts, mostly annotated. Several immune-genes involved in NIK/NF-kappaB signalling, Toll-like receptor signalling and apoptosis were activated at LC-10 h. However, down-regulation of genes encoding lysozyme, histones, cathepsins and heat shock proteins indicated signals of immunodepression, which persisted at LC-24 h, when only down-regulated genes were detected. A rebound of haemocyte activity occurred at LC-8 d as shown by up-regulation of genes involved in cytoskeleton organization and cell survival. The WI study showed a more complex picture, and several immune-relevant processes including cytoskeleton organization, cell survival, apoptosis, encapsulation, cell redox- and lipid-homeostasis were activated, illustrating the main mechanism of host response. Our results provide useful information, including potential biomarkers, to develop strategies for controlling Manila clam perkinsosis.

Protein expression profiling in haemocytes and plasma of the Manila clam Ruditapes philippinarum in response to infection with Perkinsus olseni

The protein expression profiling in clam haemocytes and plasma in response to Perkinsus olseni was addressed. Adult Manila clams from a P. olseni-free bed were experimentally challenged with parasite zoospores to analyse immune response. In another experiment, the effects of longer term infection were assessed in adult clams collected from a P. olseni-affected bed, by comparing moderate to very heavily infected clams with non-infected ones. Haemocyte and plasma proteins were separated by two-dimensional electrophoresis; spot patterns were qualitatively compared between treatments within each experiment and the spots indicating differential protein expression associated with P. olseni challenge or with field infection were processed for protein identification. Fifteen clam proteins (four in haemocytes and eleven in plasma) of which expression was markedly affected by P. olseni were identified. Some of the identified proteins have a well-known role in clam immune response against the parasite, such as lysozyme and lectins. Rho GTPase-activating protein 6 could be a marker of resistance against P. olseni, which should be further studied.

Transcriptional response of four C1q domain containing protein (C1qDC) genes from Venerupis philippinarum exposed to the water soluble fraction of No.0 diesel oil

As pattern recognitionreceptors, the C1q-domain-containing (C1qDC) proteins play an important role in the pathogen recognition and complement pathway activation. In the present study, four novel C1q domain containing proteins (designated as VpC1qDC1, VpC1qDC2, VpC1qDC3 and VpC1qDC4) were cloned and characterized from clam Venerupis philippinarum. The four VpC1qDCs all possessed the conserved features critical for the fundamental structure and function of the C1q family. The four VpC1qDCs genes showed differential response profiles after exposure to the water soluble fraction of No.0 diesel oil (WSFD). More notably, VpC1qDC1 and VpC1qDC3 were more sensitive to low concentration of WSFD, as their mRNA level changed by higher magnitudes. In addition, VpC1qDC2 and VpC1qDC4 displayed notable increases with larger amplitude to high concentration of WSFD. All these results suggested that the transcriptional response of VpC1qDCs genes were probably a protective mechanism of the cell to oils pollution. The diverse expression patterns of VpC1qDCs demonstrated that VpC1qDC1 and VpC1qDC3 were sensitive responders to environmental stress in V. philippinarum.

Identification and characterization of a ferritin gene involved in the immune defense response of scallop Chlamys farreri

Scallop Chlamys farreri is an important aquaculture species in northern China. However, its mass mortality caused by several pathogens can result in great economic loss and negative impacts to the sustainable development of the scallop industry. Thus, improving the overall understanding of immune response mechanisms involved in host-pathogen interactions is necessary. Ferritins are conserved molecules in organisms that are involved in diverse biological processes, such as mediating host-pathogen responses. In this study, we report a novel ferritin gene from C. farreri (denoted as CfFER). The full length of CfFER is 848 bp and contains a 5'-UTR of 113 bp, a 3'-UTR of 219 bp, and a complete open reading frame (ORF) of 516 bp. The ORF encodes a polypeptide of 171 amino acid residues with a molecular weight of approximately 19.95 kDa and an isoelectric point of 5.07. The CfFER protein exhibited typical ferritin structures, namely, a ferroxidase diiron center, a ferrihydrite nucleation center, and an iron-binding response signature. Phylogenetic analysis revealed that CfFER was closely related to other mollusk ferritin proteins. Expression of CfFER in different tissues was analyzed by quantitative real-time PCR, and results showed that CfFER was ubiquitously expressed in all examined tissues. The highest and lowest expression levels of CfFER were measured in the muscle and hemocyte, respectively. The relative mRNA expression of CfFER in response to bacterial (Vibrio anguillarum) and viral (acute viral necrobiotic virus) challenges sharply increased by ca. 5-fold about12 h post-infection (hpi) and then normalized at 48 hpi. Western blot analysis with polyclonal antibodies generated from the recombinant product of CfFER also demonstrated the presence of ferritin protein in hemocytes. These findings strongly suggest that CfFER is involved in the immune response of C. farreri and protection against pathogen challenge.

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.

Oral immunostimulation of the oyster Ostrea edulis: Impacts on the parasite Bonamia ostreae

Bioactive compounds were orally administered to the native European oyster Ostrea edulis to evaluate the immune response and the progression of infection of the protozoan parasite Bonamia ostreae. The immunostimulants lipopolysaccharide and zymosan directly administrated to the water column induced an increase in lysozyme activity and the percentage of granulocytes in naïve oysters over a period of 7 days. In another set of experiments, zymosan and curdlan were microencapsulated in alginate and also administered to the water column to naïve and B. ostreae infected O. edulis. Oyster mortality, prevalence and intensity of infection and several immune parameters were evaluated up to 28 days post-administration. Lysozyme activity, nitric oxide production and the expression of galectin, lysozyme and superoxide dismutase increased after 24 h in both infected and uninfected oysters. Zymosan immunostimulated oysters displayed a decrease in the prevalence of B. ostreae infection not attributed to mortalities but which could be associated to the enhanced ability of immunostimulants to evoke an enhanced immune response in the oysters and reduce infection.

An immune-enriched oligo-microarray analysis of gene expression in Manila clam (Venerupis philippinarum) haemocytes after a Perkinsus olseni challenge

Parasites of the genus Perkinsus cause high mortality and economic losses in bivalves commonly produced in global aquaculture. Although the immune responses of oysters and clams naturally infected with Perkinsus marinus or Perkinsus olseni have been extensively studied, there is not much information on host response at the early stages of infection. In this study, we analysed how P. olseni influences the gene expression profiles of haemocytes from the Manila clam (Venerupis philippinarum) using temporal experimental infections and an immune-enriched microarray. We identified an early phase of infection that was characterised by no mortality and by the increased expression of genes associated with pathogen recognition, production of nitrogen radicals and antimicrobial activity. Cellular processes such as inhibition of serine proteases and proliferation were also involved in this early response. This phase was followed by an intermediate stage, when the pathogen was most likely multiplying and infecting new areas of the body, and animals began to die. In this stage, many genes related to cell movement were over-expressed. Thirty days after infection metabolic pathway genes were the most affected. Apoptosis appears to be important during pathogenesis. Our results provide novel observations of the broader innate immune response triggered by P. olseni at different infection stages.

The genome of the Pacific oyster Crassostrea gigas brings new insights on the massive expansion of the C1q gene family in Bivalvia

C1q domain-containing (C1qDC) proteins are regarded as important players in the innate immunity of bivalve mollusks and other invertebrates and their highly adaptive binding properties indicate them as efficient pathogen recognition molecules. Although experimental studies support this view, the molecular data available at the present time are not sufficient to fully explain the great molecular diversification of this family, present in bivalves with hundreds of C1q coding genes. Taking advantage of the fully sequenced genome of the Pacific oyster Crassostrea gigas and more than 100 transcriptomic datasets, we: (i) re-annotated the oyster C1qDC loci, thus identifying the correct genomic organization of 337 C1qDC genes, (ii) explored the expression pattern of oyster C1qDC genes in diverse developmental stages and adult tissues of unchallenged and experimentally treated animals; (iii) investigated the expansion of the C1qDC gene family in all major bivalve subclasses. Overall, we provide a broad description of the functionally relevant features of oyster C1qDC genes, their comparative expression levels and new evidence confirming that a gene family expansion event has occurred during the course of Bivalve evolution, leading to the diversification of hundreds of different C1qDC genes in both the Pteriomorphia and Heterodonta subclasses.

A novel C1q-domain-containing (C1qDC) protein from Mytilus coruscus with the transcriptional analysis against marine pathogens and heavy metals

The C1q-domain-containing (C1qDC) proteins, which are involved in various processes of vertebrates, are important pattern recognition receptors in innate immunity of invertebrates. In present study, a novel C1qDC was identified from Mytilus coruscus (designated as McC1qDC), which was 917 bp in length encoding 236 amino acids with a typical signal peptide of 19 amino acid residues in N-terminus. Based on its conserved C1q domain and molecular architecture of 10 β-strand jelly-roll folding topology structure, McC1qDC might be classified as a member of the C1q family. The mRNA transcript of McC1qDC was predominantly detectable in the hemocytes, and a less degree in gill, gonad and mantle, but trace in foot, adductor and digestive gland. Upon induction by Vibrio harveyi and Vibrio alginolyticus, McC1qDC expression was significantly up-regulated. Time-dependent mRNA expression of McC1qDC was found during copper and cadmium exposure for its heavy metal-binding domain. These results indicated that McC1qDC was a novel member of the C1qDC protein family as a pattern recognition receptor against pathogens, and might be developed as a potential indicator for monitoring heavy metals pollution.

Proteomic characterization of the hemolymph of Octopus vulgaris infected by the protozoan parasite Aggregata octopiana

The immune system of cephalopods is poorly known to date. The lack of genomic information makes difficult to understand vital processes like immune defense mechanisms and their interaction with pathogens at molecular level. The common octopus Octopus vulgaris has a high economic relevance and potential for aquaculture. However, disease outbreaks provoke serious reductions in production with potentially severe economic losses. In this study, a proteomic approach is used to analyze the immune response of O. vulgaris against the coccidia Aggregata octopiana, a gastrointestinal parasite which impairs the cephalopod nutritional status. The hemocytes and plasma proteomes were compared by 2-DE between sick and healthy octopus. The identities of 12 differentially expressed spots and other 27 spots without significant alteration from hemocytes, and 5 spots from plasma, were determined by mass spectrometry analysis aided by a six reading-frame translation of an octopus hemocyte RNA-seq database and also public databases. Principal component analysis pointed to 7 proteins from hemocytes as the major contributors to the overall difference between levels of infection and so could be considered as potential biomarkers. Particularly, filamin, fascin and peroxiredoxin are highlighted because of their implication in octopus immune defense activity. From the octopus plasma, hemocyanin was identified. This work represents a first step forward in order to characterize the protein profile of O. vulgaris hemolymph, providing important information for subsequent studies of the octopus immune system at molecular level and also to the understanding of the basis of octopus tolerance-resistance to A. octopiana.

BIOLOGICAL SIGNIFICANCE

The immune system of cephalopods is poorly known to date. The lack of genomic information makes difficult to understand vital processes like immune defense mechanisms and their interaction with pathogens at molecular level. The study herein presented is focused to the comprehension of the octopus immune defense against a parasite infection. Particularly, it is centered in the host-parasite relationship developed between the octopus and the protozoan A. octopiana, which induces severe gastrointestinal injuries in octopus that produce a malabsorption syndrome. The common octopus is a commercially important species with a high potential for aquaculture in semi-open systems, and this pathology reduces the condition of the octopus populations on-growing in open-water systems resulting in important economical loses. This is the first proteomic approach developed on this host-parasite relationship, and therefore, the contribution of this work goes from i) ecological, since this particular relationship is tending to be established as a model of host-parasite interaction in natural populations; ii) evolutionary, due to the characterization of immune molecules that could contribute to understand the functioning of the immune defense in these highly evolved mollusks; and iii) to economical view. The results of this study provide an overview of the octopus hemolymph proteome. Furthermore, proteins influenced by the level of infection and implicated in the octopus cellular response are also showed. Consequently, a set of biomarkers for disease resistance is suggested for further research that could be valuable for the improvement of the octopus culture, taken into account their high economical value, the declining of landings and the need for the diversification of reared species in order to ensure the growth of the aquaculture activity. Although cephalopods are model species for biomedical studies and possess potential in aquaculture, their genomes have not been sequenced yet, which limits the application of genomic data to research important biological processes. Similarly, the octopus proteome, like other non-model organisms, is poorly represented in public databases. Most of the proteins were identified from an octopus' hemocyte RNA-seq database that we have performed, which will be the object of another manuscript in preparation. Therefore, the need to increase molecular data from non-model organisms is herein highlighted. Particularly, here is encouraged to expand the knowledge of the genomic of cephalopods in order to increase successful protein identifications. This article is part of a Special Issue entitled: Proteomics of non-model organisms.

mRNA-Seq and microarray development for the Grooved Carpet shell clam, Ruditapes decussatus: a functional approach to unravel host-parasite interaction

Background

The Grooved Carpet shell clam Ruditapes decussatus is the autochthonous European clam and the most appreciated from a gastronomic and economic point of view. The production is in decline due to several factors such as Perkinsiosis and habitat invasion and competition by the introduced exotic species, the manila clam Ruditapes philippinarum. After we sequenced R. decussatus transcriptome we have designed an oligo microarray capable of contributing to provide some clues on molecular response of the clam to Perkinsiosis.

Results

A database consisting of 41,119 unique transcripts was constructed, of which 12,479 (30.3%) were annotated by similarity. An oligo-DNA microarray platform was then designed and applied to profile gene expression in R. decussatus heavily infected by Perkinsus olseni. Functional annotation of differentially expressed genes between those two conditionswas performed by gene set enrichment analysis. As expected, microarrays unveil genes related with stress/infectious agents such as hydrolases, proteases and others. The extensive role of innate immune system was also analyzed and effect of parasitosis upon expression of important molecules such as lectins reviewed.

Conclusions

This study represents a first attempt to characterize Ruditapes decussatus transcriptome, an important marine resource for the European aquaculture. The trancriptome sequencing and consequent annotation will increase the available tools and resources for this specie, introducing the possibility of high throughput experiments such as microarrays analysis. In this specific case microarray approach was used to unveil some important aspects of host-parasite interaction between the Carpet shell clam and Perkinsus, two non-model species, highlighting some genes associated with this interaction. Ample information was obtained to identify biological processes significantly enriched among differentially expressed genes in Perkinsus infected versus non-infected gills. An overview on the genes related with the immune system on R. decussatus transcriptome is also reported.

A preliminary study of differentially expressed genes of the scallop Chlamys farreri against acute viral necrobiotic virus (AVNV)

The scallop Chlamys farreri is one of the most important aquaculture species in northern coastal provinces. However, the sustainable development of scallop industry is currently threatened by a notorious pathogen named as acute viral necrobiotic virus (AVNV), which often causes mass mortality of the animals. Despite that great attention has been focused on this novel pathogen, little knowledge about the host-virus interactions is available. In this study, suppression subtractive hybridization (SSH) was employed to identify the up-regulated differentially expressed genes in the hemocytes of C. farreri challenged by AVNV. A forward subtracted cDNA library was finally constructed and 288 positive colonies representing differentially genes were screened to perform sequencing. A total of 275 ESTs were used for further analysis using bioinformatics tools after vector screening, among which 167 ESTs could be finally identified, with significant match (E values <1 × 10(-3)) to the deposited genes (proteins) in the corresponding databases. These genes could be classified into ten categories according to their Gene Ontology annotations of biological processes and molecular functions, i.e. cell defense and homeostasis (13.82%), cellular protein metabolic process (14.90), cellular metabolism (13.09%), cytoskeletal or cellular component (5.82%), transcription regulation or RNA processing (2.18%), cell division (meiosis)/apoptosis (2.18%), DNA metabolic process and repair (1.45%), cell adhesion/signaling (1.09%), microsatellite (0.73%), and ungrouped or unknown functions (6.88). The possible biological significance of some novel genes (mainly immune and homeostasis related genes) in the host response to AVNV were discussed. This study is the first global analysis of differentially expressed genes in hemocytes from AVNV-infected C. farreri, and in addition to increasing our understanding of the molecular pathogenesis of this virus-associated scallop disease, the results presented here should provide new insights into the molecular basis of host-pathogen interactions in C. farreri.

A C1q domain containing protein from scallop Chlamys farreri serving as pattern recognition receptor with heat-aggregated IgG binding activity

Background

The C1q domain containing (C1qDC) proteins refer to a family of all proteins that contain the globular C1q (gC1q) domain, and participate in a series of immune responses depending on their gC1q domains to bind a variety of self and non-self binding ligands.

Methodology

In the present study, the mRNA expression patterns, localization, and activities of a C1qDC protein from scallop Chlamys farreri (CfC1qDC) were investigated to understand its possible functions in innate immunity. The relative expression levels of CfC1qDC mRNA in hemocytes were all significantly up-regulated after four typical PAMPs (LPS, PGN, β-glucan and polyI:C) stimulation. During the embryonic development of scallop, the mRNA transcripts of CfC1qDC were detected in all the stages, and the expression level was up-regulated from D-hinged larva and reached the highest at eye-spot larva. The endogenous CfC1qDC was dominantly located in the hepatopancreas, gill, kidney and gonad of adult scallop through immunofluorescence. The recombinant protein of CfC1qDC (rCfC1qDC) could not only bind various PAMPs, such as LPS, PGN, β-glucan as well as polyI:C, but also enhance the phagocytic activity of scallop hemocytes towards Escherichia coli. Meanwhile, rCfC1qDC could interact with human heat-aggregated IgG, and this interaction could be inhibited by LPS.

Conclusions

All these results indicated that CfC1qDC in C. farreri not only served as a PRR involved in the PAMPs recognition, but also an opsonin participating in the clearance of invaders in innate immunity. Moreover, the ability of CfC1qDC to interact with immunoglobulins provided a clue to understand the evolution of classical pathway in complement system.

A novel C1qDC protein acting as pattern recognition receptor in scallop Argopecten irradians

The C1q domain containing (C1qDC) proteins refer to a family of proteins containing the versatile charge pattern recognition globular C1q domain in the C-terminus, which could bind various ligands including PAMPs and trigger a serial of immune response. In this study, a novel C1qDC protein was identified from Argopecten irradians (designated as AiC1qDC-2). Its full-length cDNA was of 1062 bp with an open reading frame of 720 bp encoding a polypeptide of 240 amino acids containing a typical gC1q domain. This gC1q domain possessed the typical 10-stranded β-sandwich fold with a jelly-roll topology common to all C1q family members, and shared high homology with most of the other identified gC1q domains. The mRNA transcripts of AiC1qDC-2 were mainly detected in hepatopancreas, and also marginally detectable in mantle, gonad, adductor, gill and hemocytes. Its relative expression level in hemocytes was significantly up-regulated after challenges of fungi Pichia pastoris GS115 (P < 0.05), Gram-positive bacteria Micrococcus luteus (P < 0.05) and Gram-negative bacteria Vibrio anguillarum (P < 0.05). The recombinant protein of AiC1qDC-2 (rAiC1qDC-2) could bind various PAMPs, including LPS, PGN, polyI:C, mannan, β-1,3-glucan as well as Yeast-glucan, and displayed agglutinating activity to fungi P. pastoris GS115, Gram-positive bacteria Bacillus subtilis and Gram-negative bacteria Escherichia coli TOP10F' as well as V. anguillarum. All these results indicated that AiC1qDC-2 could function as a pattern recognition receptor to recognize various PAMPs on different pathogens in the innate immune responses of scallop, and provided new clues to understand the role of invertebrate C1qDC proteins in the ancient complement system.

New insights in flat oyster Ostrea edulis resistance against the parasite Bonamia ostreae

Bonamiosis due to the parasite Bonamia ostreae has been associated with massive mortality in flat oyster stocks in Europe. Control of the disease currently relies on disease management practices and transfer restriction. Previously, massal selections based on survival to challenge to infection with B. ostreae have been applied to produce flat oyster families with resistant progeny. In an attempt to understand the molecular mechanisms involved in disease resistance, differentially expressed sequence tags between resistant and wild Ostrea edulis haemocytes, both infected with the parasite, were identified using suppression subtractive hybridisation. Expression of seven ESTs has been studied using quantitative reverse-transcriptase PCR. The base-line expression of an extracellular superoxide dismutase, inhibitor of apoptosis (OeIAP), Fas ligand (OeFas-ligand) and Cathepsin B was significantly increased, whilst cyclophilin B appeared significantly decreased in resistant oysters. Considering their great interest for further studies, the open reading frames of the OeFas-ligand and OeIAP were completely sequenced.

Transcriptomics of in vitro immune-stimulated hemocytes from the Manila clam Ruditapes philippinarum using high-throughput sequencing

Background

The Manila clam (Ruditapes philippinarum) is a worldwide cultured bivalve species with important commercial value. Diseases affecting this species can result in large economic losses. Because knowledge of the molecular mechanisms of the immune response in bivalves, especially clams, is scarce and fragmentary, we sequenced RNA from immune-stimulated R. philippinarum hemocytes by 454-pyrosequencing to identify genes involved in their immune defense against infectious diseases.

Methodology and Principal Findings

High-throughput deep sequencing of R. philippinarum using 454 pyrosequencing technology yielded 974,976 high-quality reads with an average read length of 250 bp. The reads were assembled into 51,265 contigs and the 44.7% of the translated nucleotide sequences into protein were annotated successfully. The 35 most frequently found contigs included a large number of immune-related genes, and a more detailed analysis showed the presence of putative members of several immune pathways and processes like the apoptosis, the toll like signaling pathway and the complement cascade. We have found sequences from molecules never described in bivalves before, especially in the complement pathway where almost all the components are present.

Conclusions

This study represents the first transcriptome analysis using 454-pyrosequencing conducted on R. philippinarum focused on its immune system. Our results will provide a rich source of data to discover and identify new genes, which will serve as a basis for microarray construction and the study of gene expression as well as for the identification of genetic markers. The discovery of new immune sequences was very productive and resulted in a large variety of contigs that may play a role in the defense mechanisms of Ruditapes philippinarum.

Massively parallel RNA sequencing identifies a complex immune gene repertoire in the lophotrochozoan Mytilus edulis

The marine mussel Mytilus edulis and its closely related sister species are distributed world-wide and play an important role in coastal ecology and economy. The diversification in different species and their hybrids, broad ecological distribution, as well as the filter feeding mode of life has made this genus an attractive model to investigate physiological and molecular adaptations and responses to various biotic and abiotic environmental factors. In the present study we investigated the immune system of Mytilus, which may contribute to the ecological plasticity of this species. We generated a large Mytilus transcriptome database from different tissues of immune challenged and stress treated individuals from the Baltic Sea using 454 pyrosequencing. Phylogenetic comparison of orthologous groups of 23 species demonstrated the basal position of lophotrochozoans within protostomes. The investigation of immune related transcripts revealed a complex repertoire of innate recognition receptors and downstream pathway members including transcripts for 27 toll-like receptors and 524 C1q domain containing transcripts. NOD-like receptors on the other hand were absent. We also found evidence for sophisticated TNF, autophagy and apoptosis systems as well as for cytokines. Gill tissue and hemocytes showed highest expression of putative immune related contigs and are promising tissues for further functional studies. Our results partly contrast with findings of a less complex immune repertoire in ecdysozoan and other lophotrochozoan protostomes. We show that bivalves are interesting candidates to investigate the evolution of the immune system from basal metazoans to deuterostomes and protostomes and provide a basis for future molecular work directed to immune system functioning in Mytilus.

Specific versus non-specific immune responses in an invertebrate species evidenced by a comparative de novo sequencing study

Our present understanding of the functioning and evolutionary history of invertebrate innate immunity derives mostly from studies on a few model species belonging to ecdysozoa. In particular, the characterization of signaling pathways dedicated to specific responses towards fungi and Gram-positive or Gram-negative bacteria in Drosophila melanogaster challenged our original view of a non-specific immunity in invertebrates. However, much remains to be elucidated from lophotrochozoan species. To investigate the global specificity of the immune response in the fresh-water snail Biomphalaria glabrata, we used massive Illumina sequencing of 5′-end cDNAs to compare expression profiles after challenge by Gram-positive or Gram-negative bacteria or after a yeast challenge. 5′-end cDNA sequencing of the libraries yielded over 12 millions high quality reads. To link these short reads to expressed genes, we prepared a reference transcriptomic database through automatic assembly and annotation of the 758,510 redundant sequences (ESTs, mRNAs) of B. glabrata available in public databases. Computational analysis of Illumina reads followed by multivariate analyses allowed identification of 1685 candidate transcripts differentially expressed after an immune challenge, with a two fold ratio between transcripts showing a challenge-specific expression versus a lower or non-specific differential expression. Differential expression has been validated using quantitative PCR for a subset of randomly selected candidates. Predicted functions of annotated candidates (approx. 700 unisequences) belonged to a large extend to similar functional categories or protein types. This work significantly expands upon previous gene discovery and expression studies on B. glabrata and suggests that responses to various pathogens may involve similar immune processes or signaling pathways but different genes belonging to multigenic families. These results raise the question of the importance of gene duplication and acquisition of paralog functional diversity in the evolution of specific invertebrate immune responses.

Morphological characterization and functional immune response of the carpet shell clam (Ruditapes decussatus) haemocytes after bacterial stimulation

The morphology and functionality of Ruditapes decussatus haemocytes have been characterized by light microscopy and flow cytometry, leading to the identification of three different cellular subpopulations. Granulocytes were the largest cells, the hyalinocytes were smaller and contained fewer granules and the intermediate cells showed a size similar to hyalinocytes and a higher number of granules. The phagocytosis of different particles and the associated production of oxygen radicals were measured by flow cytometric methods. Granulocytes were the most active cells, followed by the intermediate cells and hyalinocytes. The effect of stimulation of haemocytes with lipopolysaccharide (LPS), with a heat inactivated bacterial mixture or with the infection of Vibrio splendidus on the cell viability and the expression of selected immune-related genes were studied. While significant low levels of damaged cells were registered in LPS-stimulated cells, the treatment with dead bacteria or V. splendidus reduced cell viability 1 h, 3 h and 6 h after treatment. The stimulation of haemocytes with LPS and dead bacteria induced changes in the expression of defender against cell death (DAD-1), thrombin, prosaposin, inhibitor of apoptosis (IAP), factor B and C3 complement component.

Gene expression analysis of clams Ruditapes philippinarum and Ruditapes decussatus following bacterial infection yields molecular insights into pathogen resistance and immunity

The carpet shell clam (Ruditapes decussatus) and Manila clam (Ruditapes philippinarum), which are cultured bivalve species with important commercial value, are affected by diseases that result in large economic losses. Because the molecular mechanism of the immune response of bivalves, especially clams, is scarce and fragmentary, we have examined all Expressed Sequence Tags (EST) resources available in public databases for these two species in order to increase our knowledge on genes related with the immune function in these animals. After automatic annotation and classification of the 3784 not-annotated ESTs of R. decussatus and 4607 of R. philippinarum found in GenBank, 424 ESTs of R. decussatus and 464 of R. philippinarum were found to be putatively involved in immune response. These were carefully reviewed and reannotated. As a result, 13 immune-related ESTs were selected and studied to compare the immune response of R. decussatus and R. philippinarum following a Vibrio alginolyticus challenge. Quantitative PCR was performed, and the expression of each EST was determined. The results showed that, in R. philippinarum, the immune response seems to be faster than that in R. decussatus. Additionally, expression of NF-κB activating genes in R. decussatus did not seem to be sufficient to promote an immune response after Vibrio infection. R. philippinarum, however, was able to trigger and efficiently regulate the transcriptional activity of NF-κB, even when low expression values were reported.

Diversity of coding sequences and gene structures of the antifungal peptide mytimycin (MytM) from the Mediterranean mussel, Mytilus galloprovincialis

Knowledge on antifungal biomolecules is limited compared to antibacterial peptides. A strictly antifungal peptide from the blue mussel, Mytilus edulis named mytimycin (MytM) was reported in 1996 as partial NH(2) 33 amino acid sequence. Using back-translations of the previous sequence, MytM-related nucleotide sequences were identified from a normalized Mytilus galloprovincialis expressed sequence tag library. Primers designed from a consensus sequence have been used to obtain a fragment of 560 nucleotides, including the complete coding sequence of 456 nucleotides. Precursor is constituted by a signal peptide of 23 amino acids, followed by MytM of 54 amino acids (6.2-6.3 kDa, 12 cysteines) and C-terminal extension of 75 amino acids. Only two major amino acid precursor sequences emerged, one shared by M. galloprovincialis from Venice and Vigo, the other belonging to M. galloprovincialis from Palavas, with nine amino acid differences between the two MytM. Predicted disulfide bonds suggested the presence of two constrained domains joined by amino acidic NIFG track. Intriguing was the presence of conserved canonical EF hand-motif located in the C-terminus extension of the precursor. The MytM gene was found interrupted by two introns. Intron 2 existed in two forms, a long (1,112 nucleotides) and a short (716 nucleotides) one resulting from the removal of the central part of the long one. Both the short (GenBank FJ804479) and the long (GenBank FJ804478) genes are simultaneously present in the mussel genome.

Genomic organization, molecular diversification, and evolution of antimicrobial peptide myticin-C genes in the mussel (Mytilus galloprovincialis)

Myticin-C is a highly variable antimicrobial peptide associated to immune response in Mediterranean mussel (Mytilus galloprovincialis). In this study, we tried to ascertain the genetic organization and the mechanisms underlying myticin-C variation and evolution of this gene family. We took advantage of the large intron size variation to find out the number of myticin-C genes. Using fragment analysis a maximum of four alleles was detected per individual at both introns in a large mussel sample suggesting a minimum of two myticin-C genes. The transmission pattern of size variants in two full-sib families was also used to ascertain the number of myticin-C genes underlying the variability observed. Results in both families were in accordance with two myticin-C genes organized in tandem. A more detailed analysis of myticin-C variation was carried out by sequencing a large sample of complementary (cDNA) and genomic DNA (gDNA) in 10 individuals. Two basic sequences were detected at most individuals and several sequences were constituted by combination of two different basic sequences, strongly suggesting somatic recombination or gene conversion. Slight within-basic sequence variation detected in all individuals was attributed to somatic mutation. Such mutations were more frequently at the C-terminal domain and mostly determined non-synonymous substitutions. The mature peptide domain showed the highest variation both in the whole cDNA and in the basic-sequence samples, which is in accordance with the pathogen recognition function associated to this domain. Although most tests suggested neutrality for myticin-C variation, evidence indicated positive selection in the mature peptide and C-terminal region. Three main highly supported clusters were observed when reconstructing phylogeny on basic sequences, meiotic recombination playing a relevant role on myticin-C evolution. This study demonstrates that mechanisms to generate molecular variation similar to that observed in vertebrates are also operating in molluscs.

Molecular responses of Ostrea edulis haemocytes to an in vitro infection with Bonamia ostreae

Bonamiosis due to the parasite Bonamia ostreae is a disease affecting the flat oyster Ostrea edulis. B. ostreae is a protozoan, affiliated to the order of haplosporidia and to the cercozoan phylum. This parasite is mainly intracellular, infecting haemocytes, cells notably involved in oyster defence mechanisms. Suppression subtractive hybridisation (SSH) was carried out in order to identify oyster genes differentially expressed during an infection of haemocytes with B. ostreae. Forward and reverse banks allowed obtaining 1104 and 1344 clones respectively, among which 391 and 480 clones showed a differential expression between both tested conditions (haemocytes alone versus haemocytes in contact with parasites). ESTs of interest including genes involved in cytoskeleton, respiratory chain, detoxification membrane receptors, and immune system were identified. The open reading frames of two selected genes (galectin and IRF-like) were completely sequenced and characterized. Real time PCR assays were developed to study the relative expression of candidate ESTs during an in vitro infection of haemocytes by live and dead parasites. Haemocyte infection with B. ostreae induced an increased expression of omega glutathione S-transferase (OGST), superoxide dismutase (SOD), tissue inhibitor of metalloproteinase (TIMP), galectin, interferon regulatory factor (IRF-like) and filamin genes.

Gene expression patterns of abalone, Haliotis tuberculata, during successive infections by the pathogen Vibrio harveyi

Since 1998, episodic mass mortality of the abalone Haliotistuberculata has been observed along the northern Brittany coast of France caused by a complex interaction among the host, pathogen and environmental factors. In the present study, abalone were submitted to two successive infections with the pathogen Vibrioharveyi under controlled conditions. During the first challenge, infection by V.harveyi resulted in 64% mortality of mature abalone. After a second infection of those surviving the first challenge, only 44% mortality was observed. Physiological variability in the host response appears to be a major determinant in susceptibility to V.harveyi. In order to isolate differentially expressed genes in H.tuberculata challenged with this bacterium, suppression subtractive hybridization (SSH) cDNA libraries were constructed from muscle of moribund abalone (susceptibles), surviving individuals (apparently resistant to the bacterium) and control (unexposed) animals. Of the 1152 clones sequenced, 218 different partial cDNA sequences were obtained and represented 69 known genes. Of these, 65 were identified for the first time in H.tuberculata. Using real-time PCR, a time-course study was conducted on 19 of the genes identified by SSH. A majority of differentially expressed transcripts were down-regulated in susceptible individuals as compared to their resistant counterparts. Bacterial challenge of abalone resulted in the up-regulation of three transcripts (encoding ferritin, heat shock protein HSP84 and fatty acid binding protein FABP) in those that survived exposure to V.harveyi. This study has identified potential candidates for further investigation into the functional basis of resistance and susceptibility to summer vibriosis outbreaks in abalone.

MgC1q, a novel C1q-domain-containing protein involved in the immune response of Mytilus galloprovincialis

In this study, we present the characterization of a newly identified gene, MgC1q, revealed in suppression subtractive hybridization and cDNA libraries from immunostimulated mussels. Based on sequence homology, molecular architecture and domain similarity, this new C1q-domain-containing gene may be classified as a member of the C1q family and, therefore, part of the C1q-TNF superfamily. The expression of MgC1q was detected all along the mussel ontogeny, being detectable within 2h post-fertilization, with a notable increase after 1 month and continuing to increase until 3 months. Measurable transcript levels were also evident in all analyzed tissues of naïve adult mussels, and the hemocytes showed the highest expression levels. Experimental infection of adult mussels with Gram positive or Gram negative bacteria significantly modulated the MgC1q expression, and confirmed it as an immune-related gene. Intra- and inter-individual sequence analyses revealed extraordinary diversity of MgC1q at both the DNA and cDNA levels. While further research is needed to define its function, our data indicate that MgC1q is a pattern recognition molecule able to recognize pathogens during innate immune responses in Myitilus galloprovincialis. The high sequence variability suggests that somatic diversification of these nonself recognition molecules could have occurred.

Expression of Mytilus immune genes in response to experimental challenges varied according to the site of collection

Mussels live in diverse coastal environments experience various physical, chemical and biological conditions, which they counteract with functional adjustments and heritable adaptive changes. In order to investigate possible differences in immune system capabilities, we analyzed by qPCR the expression levels of 4 immune genes (defensin, mytilin B, myticin B, lysozyme) and HSP70 in the Mediterranean mussel, Mytilus galloprovincialis collected in 3 European farming areas {Atlantic Ocean-Ría de Vigo-Spain (RV), French Mediterranean Gulf of Lion-Palavas-Prévost lagoon (PP) and Northern Adriatic Sea-Venice-Italy (VI)} in response to one injection of one of the 3 bacterial species (Vibrio splendidus LGP32, Vibrio anguillarum, Micrococcus lysodeikticus), and to heat shock or cold stress. We confirmed that the 5 genes are constitutively expressed in hemocytes, defensin being the less expressed, myticin B the highest. As suspected, the same gene resulted differently expressed according to mussel group, with the biggest difference being for HSP70 and lysozyme and lowest expression of all the 5 genes in mussels from RV. In addition, gene expression levels varied according to the challenge. Most frequent effect of bacterial injections was down-regulation, especially for mytilin B and myticin B. Heat shock enhanced transcript levels, particularly in mussels from RV, whereas cold stress had no effect. In situ hybridization of labelled probes on mussel hemocytes indicated that bacterial injections did not change the mRNA patterns of defensin and myticin B whereas mytilin B mRNA almost disappeared. In conclusion, these results demonstrated that constitutive level, nature and intensity of immune gene expression regulations strongly depended from mussel group, and support the concept of gene-environment interactions.

Identification and expression of differentially expressed genes in the hard clam, Mercenaria mercenaria, in response to quahog parasite unknown (QPX)

Background

The hard clam, Mercenaria mercenaria, has been affected by severe mortality episodes associated with the protistan parasite QPX (Quahog Parasite Unknown) for several years. Despite the commercial importance of hard clams in the United States, molecular bases of defense mechanisms in M. mercenaria, especially during QPX infection, remain unknown.

Results

Our study used suppression subtractive hybridization (SSH), as well as the construction of cDNA libraries from hemocytes to identify genes related to the defense of the hard clam against its parasite. Hard clams were experimentally infected with QPX and SSH was performed on mRNA samples extracted from mantle and gill tissues at different times post-challenge. A total of 298 clones from SSH libraries and 1352 clones from cDNA libraries were sequenced. Among these sequences, homologies with genes involved in different physiological processes related to signal transduction, stress response, immunity and protein synthesis were identified. Quantitative PCR revealed significant changes in the expression of several of these genes in response to QPX challenge and demonstrated significant correlations in terms of levels of gene expression between intermediates of signalling pathways and humoral defense factors, such as big defensin and lysozyme.

Conclusion

Results of this study allowed the detection of modifications caused by QPX at the transcriptional level providing insight into clam immune response to the infection. These investigations permitted the identification of candidate genes and pathways for further analyses of biological bases of clam resistance to QPX allowing for a better understanding of bivalve immunity in general.