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

Increased seawater temperature triggers thermal, oxidative and metabolic response of Ostrea edulis, leading to anaerobiosis

Bivalves are among the marine organisms most influenced by climate change. Despite the flat oyster's Ostrea edulis high economic value, its culture is developed on a very small scale, since this species possesses a strong susceptibility to abiotic stressors. Due to climate change, temperature is one of the most critical environmental parameters for the welfare of the Mediterranean basin's marine inhabitants. The present study's purpose was to investigate the physiological performance of the Mediterranean's native O. edulis as it faces exposure to different temperatures. Since juveniles are more susceptible to abiotic stressors, this experimental procedure was focused on young individuals. The seawater temperatures studied included a standard control temperature of 21 °C (often observed in several marine areas throughout the Mediterranean), as well as increased seawater temperatures of 25 °C and 28 °C, occasionally occurring in shallow Mediterranean waters inhabited by bivalve spat. These were selected since the tissues of O. edulis becomes partly anaerobic in temperatures exceeding 26 °C, while cardiac dysfunction (arrhythmia) emerges at 28 °C. The results demonstrate that temperatures above 25 °C trigger both the transcriptional upregulation of hsp70 and hsp90, and the antioxidant genes Cu/Zn sod and catalase. Enhancement of thermal tolerance and increased defense against increased ROS production during thermal stress, were observed. As the intensity and duration of thermal stress increases, apoptotic damage may also occur. The increased oxidative and thermal stress incurred at the highest temperature of 28 °C, seemed to trigger the switch from aerobic to anaerobic metabolism, reflected by higher pepck mRNA expressions and lower ETS activity.

Genome-Wide Identification and Characterization of MITF Genes in Ruditapes philippinarum and Their Involvement in the Immune Response to Vibrio anguillarum Infection

Studies have shown that the shellfish have innate immune system, which is a very important immune form of shellfish, and they rely on the innate immune system to resist diseases. As a transcription factor, Microphthalmia-associated transcription factor (MITF) plays a regulatory role in immune response and the shell color is also an important index for the breeding of excellent varieties of R. philippinarum. The research on immune response mechanism of RPMITFs can provide important reference data for the breeding of excellent clam varieties. In the genome of R. philippinarum, the RPMITF genes family of shell color-related gene family was selected as the target gene of this experiment. There are 12 RpMITF genes named RpMITF1, RpMITF2, RpMITF3, RpMITF4, RpMITF5, RpMITF6, RpMITF7, RpMITF8, RpMITF9, RpMITF10, RpMITF11, and RpMITF12. The open reading frame length is 639, 1233, 996, 1239, 675, 624, 816, 1365, 612, 1614, 1122, and 486 bp, encoding 212, 410, 331, 412, 224, 207, 271, 454, 203, 537, 373, and 161 aa, respectively. The predicted molecular weight range of amino acids is 18.85-62.61 kda, and the isoelectric point range is 5.26-9.44. Real-time quantitative PCR was used to detect the gene expression of RpMITF gene family in hepatopancreas tissues of two populations of Manila clam at 6 time points (0, 3, 6, 12, 24, and 48 h) after Vibrio anguillarum stress. The results show that RpMITF gene family was significantly expressed in hepatopancreas of two clam populations after V. anguillarum stress (P < 0.05).

Multidisciplinary long-term survey of Manila clam grown in farming sites subjected to different environmental conditions

In recent years recurrent bivalve mass mortalities considerably increased around the world, causing the collapse of natural and farmed populations. Venice Lagoon has historically represented one of the major production areas of the Manila clam Ruditapes philippinarum in Europe. However, in the last 20 years a 75 % decrease in the annual production has been experienced. While climate change and anthropogenic interventions may have played a key role in natural and farmed stocks reductions, no studies investigated at multiple levels the environmental stressors affecting farmed Manila clam to date. In this work we carried out a long-term monitoring campaign on Manila clam reared in four farming sites located at different distances from the southern Venice Lagoon inlet, integrating (meta)genomic approaches (i.e. RNA-seq; microbiota characterization), biometric measurements and chemical-physical parameters. Our study allowed to characterize the molecular mechanisms adopted by this species to cope with the different environmental conditions characterizing farming sites and to propose hypotheses to explain mortality events observed in recent years. Among the most important findings, the disruption of clam's immune response, the spread of Vibrio spp., and the up-regulation of molecular pathways involved in xenobiotic metabolism suggested major environmental stressors affecting clams farmed in sites placed close to Chioggia's inlet, where highest mortality was also observed. Overall, our study provides knowledge-based tools for managing Manila clam farming on-growing areas. In addition, the collected data is a snapshot of the time immediately before the commissioning of MoSE, a system of mobile barriers aimed at protecting Venice from high tides, and will represent a baseline for future studies on the effects of MoSE on clams farming and more in general on the ecology of the Venice Lagoon.

Gene Expression Profiles in Two Razor Clam Populations: Discerning Drivers of Population Status

With rapidly changing marine ecosystems, shifts in abundance and distribution are being documented for a variety of intertidal species. We examined two adjacent populations of Pacific razor clams (Siliqua patula) in lower Cook Inlet, Alaska. One population (east) supported a sport and personal use fishery, but this has been closed since 2015 due to declines in abundance, and the second population (west) continues to support commercial and sport fisheries. We used gene expression to investigate potential causes of the east side decline, comparing razor clam physiological responses between east and west Cook Inlet. The target gene profile used was developed for razor clam populations in Alaska based on physiological responses to environmental stressors. In this study, we identified no differences of gene expression between east and west populations, leading to two potential conclusions: (1) differences in factors capable of influencing physiology exist between the east and west and are sufficient to influence razor clam populations but are not detected by the genes in our panel, or (2) physiological processes do not account for the differences in abundance, and other factors such as predation or changes in habitat may be impacting the east Cook Inlet population.

Invasive clams (Ruditapes philippinarum) are better equipped to deal with harmful algal blooms toxins than native species (R. decussatus): evidence of species-specific toxicokinetics and DNA vulnerability

This study aims to assess and compare the kinetics (accumulation/elimination) of the marine biotoxins okadaic acid (OA) and dinophysistoxin-1 (DTX1), between native (Ruditapes decussatus) and invasive (Ruditapes philippinarum) clam species, and their genotoxic effects and DNA recover capacity after, exposure to toxic dinoflagellates Prorocentrum lima. Clams were fed with P. lima for 5 days and then to non-toxic algae (post-exposure) during other 5 days. Toxin concentrations determined in clams by LC-MS/MS were related with DNA damage and repair assessment through the comet and base excision repair (BER) assays, respectively. Differential accumulation patterns were observed between the invasive and native species. The invasive species consistently and progressively accumulated the toxins during the first 24 h of exposure, while the native clams showed drastic variations in the toxin accumulation. Nevertheless, at the end of a 5 days of exposure period, the native clams presented higher toxin concentrations, nearly reaching the legal regulatory limit for human consumption. In addition, native clams were vastly affected by OA and DTX1, presenting an increment in the DNA damage since the first day, with a correspondent increase in the repair activity. On the other hand, invasive clams were not affected by the dinoflagellate toxins, exhibiting only some signs of the challenge, namely an increase in the DNA repair mechanisms in the post-exposure period. Invasive clams R. philippinarum are better adapted to cope with harmful algal blooms and OA-group toxins than native species. These results may increase farming interest and may lead to new introductions of the invasive clams. In sympatry sites, exposure to OA-group toxins may unbalance clams species biomass and distribution as exposure to toxic dinoflagellates affects the native clams from cellular to a population level, representing a significant threat to development and maintenance of R. decussatus populations.

New insights into the Manila clam and PAMPs interaction based on RNA-seq analysis of clam through in vitro challenges with LPS, PGN, and poly(I:C)

Background

Manila clam (Ruditapes philippinarum) is a worldwide commercially important marine bivalve species. In recent years, however, microbial diseases caused high economic losses and have received increasing attention. To understand the molecular basis of the immune response to pathogen-associated molecular patterns (PAMPs) in R. philippinarum, transcriptome libraries of clam hepatopancreas were constructed at 24 h post-injection with Lipopolysaccharide (LPS), peptidoglycan (PGN), and polyinosinic-polycytidylic acid (poly(I:C)) and phosphate-buffered saline (PBS) control by using RNA sequencing technology (RNA-seq).

Results

A total of 832, 839, and 188 differentially expressed genes (DEGs) were found in LPS, PGN, and poly(I:C) challenge group compared with PBS control, respectively. Several immune-related genes and pathways were activated in response to the different PAMPs, suggesting these genes and pathways might specifically participate in the immune response to pathogens. Besides, the analyses provided useful complementary data to compare different PAMPs challenges in vivo. Functional enrichment analysis of DEGs demonstrated that PAMPs responsive signal pathways were related to apoptosis, signal transduction, immune system, and signaling molecules and interaction. Several shared or specific DEGs response to different PAMPs were revealed in R. philippinarum, including pattern recognition receptors (PRRs), antimicrobial peptides (AMPs), interferon-induced proteins (IFI), and some other immune-related genes were found in the present work.

Conclusions

This is the first study employing high throughput transcriptomic sequencing to provide valuable genomic resources and investigate Manila clam response to different PAMPs through in vivo challenges with LPS, PGN, and poly(I:C). The results obtained here provide new insights to understanding the immune characteristics of R. philippinarum response to different PAMPs. This information is critical to elucidate the molecular basis of R. philippinarum response to different pathogens invasion, which potentially can be used to develop effective control strategies for different pathogens.

Clinical Applications of Antimicrobial Peptides (AMPs): Where do we Stand Now?

In this era of multi-drug resistance (MDR), antimicrobial peptides (AMPs) are one of the most promising classes of potential drug candidates to combat communicable as well as noncommunicable diseases such as cancers and diabetes. AMPs show a wide spectrum of biological activities which include antiviral, antifungal, anti-mitogenic, anticancer, and anti-inflammatory properties. Apart from these prospective therapeutic potentials, the AMPs can act as food preservatives and immune modulators. Therefore, AMPs have the potential to replace conventional drugs and may gain a significant global drug market share. Although several AMPs have shown therapeutic potential in vitro or in vivo, in most cases they have failed the clinical trial owing to various issues. In this review, we discuss in brief (i) molecular mechanisms of AMPs in various diseases, (ii) importance of AMPs in pharmaceutical industries, (iii) the challenges in using AMPs as therapeutics and how to overcome, (iv) available AMP therapeutics in market, and (v) AMPs under clinical trials. Here, we specifically focus on the therapeutic AMPs in the areas of dermatology, surgery, oncology and metabolic diseases.

A proteomic study of resistance to Brown Ring disease in the Manila clam, Ruditapes philippinarum

Marine mollusk aquaculture has more than doubled over the past twenty years, accounting for over 15% of total aquaculture production in 2016. Infectious disease is one of the main limiting factors to the development of mollusk aquaculture, and the difficulties inherent to combating pathogens through antibiotic therapies or disinfection have led to extensive research on host defense mechanisms and host-pathogen relationships. It has become increasingly clear that characterizing the functional profiles of response to a disease is an essential step in understanding resistance mechanisms and moving towards more effective disease control. The Manila clam, Ruditapes philippinarum, is a main cultured bivalve species of economic importance which is affected by Brown Ring disease (BRD), an infection induced by the bacterium Vibrio tapetis. In this study, juvenile Manila clams were subjected to a 28-day controlled challenge with Vibrio tapetis, and visual and molecular diagnoses were carried out to distinguish two extreme phenotypes within the experimental clams: uninfected ("RES", resistant) and infected ("DIS", diseased) post-challenge. Total protein extractions were carried out for resistant and diseased clams, and proteins were identified using LC-MS/MS. Protein sequences were matched against a reference transcriptome of the Manila clam, and protein intensities based on label-free quantification were compared to reveal 49 significantly accumulated proteins in resistant and diseased clams. Proteins with known roles in pathogen recognition, lysosome trafficking, and various aspects of the energy metabolism were more abundant in diseased clams, whereas those with roles in redox homeostasis and protein recycling were more abundant in resistant clams. Overall, the comparison of the proteomic profiles of resistant and diseased clams after a month-long controlled challenge to induce the onset of Brown Ring disease suggests that redox homeostasis and maintenance of protein structure by chaperone proteins may play important and interrelated roles in resistance to infection by Vibrio tapetis in the Manila clam.

Monitoring nearshore ecosystem health using Pacific razor clams (Siliqua patula) as an indicator species

An emerging approach to ecosystem monitoring involves the use of physiological biomarker analyses in combination with gene transcription assays. For the first time, we employed these tools to evaluate the Pacific razor clam (Siliqua patula), which is important both economically and ecologically, as a bioindicator species in the northeast Pacific. Our objectives were to (1) develop biomarker and gene transcription assays with which to monitor the health of the Pacific razor clam, (2) acquire baseline biomarker and gene transcription reference ranges for razor clams, (3) assess the relationship between physiological and gene transcription assays and (4) determine if site-level differences were present. Pacific razor clams were collected in July 2015 and 2016 at three sites within each of two national parks in southcentral Alaska. In addition to determining reference ranges, we found differences in biomarker assay and gene transcription results between parks and sites which indicate variation in both large-scale and local environmental conditions. Our intent is to employ these methods to evaluate Pacific razor clams as a bioindicator of nearshore ecosystem health. Links between the results of the biomarker and gene transcription assays were observed that support the applicability of both assays in ecosystem monitoring. However, we recognize the need for controlled studies to examine the range of responses in physiology and gene transcripts to different stressors.

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.

Activity Dependence of a Novel Lectin Family on Structure and Carbohydrate-Binding Properties

A GalNAc/Gal-specific lectins named CGL and MTL were isolated and characterized from the edible mussels Crenomytilus grayanus and Mytilus trossulus. Amino acid sequence analysis of these lectins showed that they, together with another lectin MytiLec-1, formed a novel lectin family, adopting β-trefoil fold. In this mini review we discuss the structure, oligomerization, and carbohydrate-binding properties of a novel lectin family. We describe also the antibacterial, antifungal, and antiproliferative activities of these lectins and report about dependence of activities on molecular properties. Summarizing, CGL, MTL, and MytiLec-1 could be involved in the immunity in mollusks and may become a basis for the elaboration of new diagnostic tools or treatments for a variety of cancers.

Development and validation of a specific real-time PCR assay for the detection of the parasite Perkinsus olseni

Perkinsus olseni is a protozoan parasite that infects a wide variety of molluscs worldwide, causing economic losses in the aquaculture sector. In the present study, a quantitative PCR (qPCR) assay was developed for the detection and quantification of P. olseni in clam gill tissue and hemolymph (Ruditapes philippinarum and R. decussatus), and the results were compared with those of the standard diagnostic methods recommended by the O.I.E. (World Organisation for Animal Health): Ray's fluid thioglycollate culture method (RFTM), a polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) assay and histopathology. The efficiency, sensitivity and reproducibility of the newly described qPCR assay were also determined. The highest prevalence was detected using the qPCR assay, and the strongest linear correlation was obtained between the RFTM infection levels and the threshold cycle (Ct) number from the gill tissue. Although better results were obtained from gill than from the hemolymph in the qPCR assays, especially with lower infection levels of the parasite, a significant linear correlation was observed between Ct values from the gill and hemolymph. The qPCR assay that was developed in this study showed high sensitivity, specificity and reproducibility for the detection and quantification of P. olseni.

High individual variability in the transcriptomic response of Mediterranean mussels to Vibrio reveals the involvement of myticins in tissue injury

Mediterranean mussels (Mytilus galloprovincialis) are sessile filter feeders that live in close contact with numerous marine microorganisms. As all invertebrates, they lack an adaptive immune response and how these animals are able to respond to a bacterial infection and discriminate it from their normal microbiome is difficult to understand. In this work, we conducted Illumina sequencing of the transcriptome of individual mussels before and after being infected with Vibrio splendidus. The control mussels were injected with filtered seawater. We demonstrate that a great variability exists among individual transcriptomes and that each animal showed an exclusive repertoire of genes not shared with other individuals. The regulated genes in both the control and infected mussels were also analyzed and, unexpectedly, the sampling before the injection was considered a stress stimulus strong enough to trigger and modulate the response in hemocytes, promoting cell migration and proliferation. We found a clear response against the injection of filtered seawater, suggesting a reaction against a tissue injury in which the myticins, the most expressed antimicrobial peptides in mussel, appeared significantly up regulated. Functional experiments with flow cytometry confirmed the transcriptomic results since a significant alteration of hemocyte structures and a decrease in the number of hemocytes positive for myticin C were found only after a Vibrio infection and not observed when mussels were bled before, generating a tissue injury. Therefore, we report the involvement of myticins in the response to a danger signal such as a simple injection in the adductor muscle.

Gene expression and genotoxicity in Manila clam (Ruditapes philippinarum) modulated by sediment contamination and lagoon dynamics in the Po river delta

The lagoons of the Po River delta are potentially exposed to complex mixtures of contaminants, nevertheless, there is a substantial lack of information about the biological effects of these contaminants in the Po delta lagoons. These environments are highly dynamic and the interactions between chemical and environmental stressors could prevent the proper identification of biological effects and their causes. In this study, we aimed to disentangle such interactions focusing on Manila clams, previously exposed to six lagoons of the Po delta, adopting three complementary tools: a) the detailed description via modelling techniques of lagoon dynamics for salinity and water temperature; b) the response sensitivity of a number of target genes (ahr, cyp4, ρ-gst, σ-gst, hsp22, hsp70, hsp90, ikb, dbh, ach, cat, Mn-sod, Cu/Zn-sod, cyp-a, flp, grx, TrxP) investigated in clam digestive glands by Real Time PCR; and c) the relevance of DNA adducts determined in clams as markers of exposure to genotoxic chemicals. The lagoons showed specific dynamics, and two of them (Marinetta and Canarin) could induce osmotic stress. A group of genes (ahr, cyp4, Mn-sod, σ-gst, hsp-22, cyp-a, TrxP) seemed to be associated with overall lagoon characteristics as may be described by salinity and its variations. Lagoon modelling and a second group of genes (hsp70, hsp90, cat, ikb, ach, grx, Cu/Zn-sod) also suggested that moderate increases of river discharge may imply worse exposure conditions. Oxidative stress seemed to be associated with such events but it was slightly evident also under normal exposure conditions. DNA adduct formation was mainly associated with overwhelmed antioxidant defences (e.g. low Cu/Zn-sod) or seemingly with their lack of response in due time. In Po delta lagoons, Manila clam can be affected by chemical and environmental factors which can contribute to induce oxidative stress, DNA adduct formation and, ultimately, to affect clam condition and health.

Effects of carbamazepine and cetirizine under an ocean acidification scenario on the biochemical and transcriptome responses of the clam Ruditapes philippinarum

Several works evaluated the toxicity of pharmaceutical drugs and climate related changes in invertebrates but few explored the combined effects of both stressors, namely considering their mode of action (MoA). Carbamazepine (CBZ) and cetirizine (CTZ) are pharmaceutical drugs detected in the environment and the toxicity derived from the combined effects of these drugs with ocean acidification (OA) is poorly explored. Thus, the present study investigated the biochemical parameters related to an oxidative stress response and the transcription of genes related to the MoA of CBZ (1.0 μg/L) and CTZ (0.6 μg/L) in the clam Ruditapes philippinarum chronically exposed (28 days) to control (7.8) and low (7.5) pH conditions. The results obtained showed that despite the clams accumulated both drugs, at low pH the clams exposed to CTZ decreased drug concentration and BCF values (CTZ uptake: 2.0 ± 0.5 ng/g fresh weight; BCF: 3.8 ± 0.9) in comparison with clams exposed to control pH (CTZ uptake: 2.9 ± 0.3 ng/g fresh weight; BCF: 5.5 ± 0.6). No oxidative stress was induced by the exposure to CBZ or CTZ at each pH level, but the transcription of several genes related with the MoA (neurotransmission, immunity and biomineralization) was altered by low pH, drug exposure and the combination of both stressors. At both pH conditions, CBZ increased the transcription of GABA receptor gene (neurotransmission) and CTZ led to a decrease of Perlucin gene (biomineralization) transcription. The transcription of MyD88 gene (immunity) decreased at low pH (7.5) combined with drug exposure (CBZ or CTZ). Thus, it was highlighted that the interaction of drug exposure and low pH conditions can change bivalves' sensitivity to drugs or alter drugs toxicity.

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.

Ancient DNA analysis identifies marine mollusc shells as new metagenomic archives of the past

Marine mollusc shells enclose a wealth of information on coastal organisms and their environment. Their life history traits as well as (palaeo-) environmental conditions, including temperature, food availability, salinity and pollution, can be traced through the analysis of their shell (micro-) structure and biogeochemical composition. Adding to this list, the DNA entrapped in shell carbonate biominerals potentially offers a novel and complementary proxy both for reconstructing palaeoenvironments and tracking mollusc evolutionary trajectories. Here, we assess this potential by applying DNA extraction, high-throughput shotgun DNA sequencing and metagenomic analyses to marine mollusc shells spanning the last ~7,000 years. We report successful DNA extraction from shells, including a variety of ancient specimens, and find that DNA recovery is highly dependent on their biomineral structure, carbonate layer preservation and disease state. We demonstrate positive taxonomic identification of mollusc species using a combination of mitochondrial DNA genomes, barcodes, genome-scale data and metagenomic approaches. We also find shell biominerals to contain a diversity of microbial DNA from the marine environment. Finally, we reconstruct genomic sequences of organisms closely related to the Vibrio tapetis bacteria from Manila clam shells previously diagnosed with Brown Ring Disease. Our results reveal marine mollusc shells as novel genetic archives of the past, which opens new perspectives in ancient DNA research, with the potential to reconstruct the evolutionary history of molluscs, microbial communities and pathogens in the face of environmental changes. Other future applications include conservation of endangered mollusc species and aquaculture management.

Native (Ruditapes decussatus) and non-indigenous (R. philippinarum) shellfish species living in sympatry: Comparison of regulated and non-regulated biotoxins accumulation

The native Ruditapes decussatus and the non-indigenous Ruditapes philippinarum are an important target of shellfish industries. The aim of this study was to compare an invader with a native species living in sympatry in the view of marine biotoxins accumulation. Samples were analysed for regulated and non-regulated biotoxins. The consistently occurrence of okadaic acid-group toxins and BMAA, may cause human health problems and economical losses. A strong positive relationship was observed between species, with significantly higher DSP toxicity in R. decussatus. Similar toxin profiles dominated by DTX3 in both species suggests similar metabolic pathways. Lower DSP toxicity in R. philippinarum may favour their cultivation, but a tendency for higher levels of the non-regulated BMAA was observed, indicating risks for consumers that are not monitored. This study highlights the need to better understand the physiological responses and adaptations allowing similar species exposed to the same conditions to present different toxicity levels.

Validation of reference genes for RT-qPCR in marine bivalve ecotoxicology: Systematic review and case study using copper treated primary Ruditapes philippinarum hemocytes

The appropriate selection of reference genes for the normalization of non-biological variance in reverse transcription real-time quantitative PCR (RT-qPCR) is essential for the accurate interpretation of the collected data. The use of multiple validated reference genes has been shown to substantially increase the robustness of the normalization. It is therefore considered good practice to validate putative genes under specific conditions, determine the optimal number of genes to be employed, and report the method or methods used. Under this premise, we assessed the current state of reference gene based normalization in RT-qPCR bivalve ecotoxicology studies (post 2011), employing a systematic quantitative literature review. A total of 52 papers met our criteria and were analysed for genes used, the use of multiple reference genes, as well as the validation method employed. We further critically discuss methods for reference gene validation based on a case study using copper exposed primary hemocytes from the marine bivalve Ruditapes philippinarum; including the established algorithms geNorm, NormFinder and BestKeeper, as well as the popular online tool RefFinder. We identified that RT-qPCR normalization is largely performed using single reference genes, while less than 40% of the studies attempted to experimentally validate the expression stability of the genes used. 18s rRNA and β-Actin were the most popular genes, yet their un-validated use did introduce artefactual variance that altered the interpretation of the resulting data. Our findings further suggest that combining the results from multiple individual algorithms and calculating the overall best-ranked gene, as computed by the RefFinder tool, does not by default lead to the identification of the most suitable reference genes.

Effects of seawater temperature increase on economically relevant native and introduced clam species

As a consequence of climate change, global warming is expected to increase during the 21^st century. Taking this into account, the impact of rising temperatures on the native Ruditapes decussatus and introduced R. philippinarum bivalve species was assessed, through biochemical and mRNA transcription analyses. Our findings showed that at 21 °C the electron transport system and antioxidant enzyme activity, as well as the expression of Hsp70 gene were induced in R. decussatus when compared with 17 °C. On the other hand, at 25 °C results suggested that R. decussatus closed their valves during short periods, as a behavioral strategy, down-regulating the expression of genes associated with mitochondrial metabolism (Cox-1 and 16S) and chaperone function (Hsp70) compared with organisms at 17 °C. In addition, the introduced species (R. philippinarum) increased the electron transport system and antioxidant activities, as well as gene expression of antioxidant enzymes and molecular chaperone (Hsp70) at 21 °C. However, antioxidant mechanisms were not enough to prevent lipid membrane damages at 21 °C. At 25 °C R. philippinarum presented increased electron transport system and antioxidant activity, as well as the expression of genes associated with apoptosis regulation and molecular chaperone. Overall, the present findings indicate that in a global warming scenario both species are able to induce different mechanisms to mitigate the impacts of temperature increase.

Moving from Histopathology to Molecular Tools in the Diagnosis of Molluscs Diseases of Concern under EU Legislation

One of the main factors limiting molluscs production is the presence of pathogens and diseases. Disease agent transfer via transfers of live molluscs has been a major cause of disease outbreaks and epizootics. Because of that, the European Union has adopted several decisions and directives, the last in 2006 (2006/88/EC) to control movements of marine organisms over the European countries. Once the disease is established in a determined area its eradication is a complicated task because life cycle of pathogens are not completely known and only a good and early diagnosis of the disease could be the most appropriate way to deal with it. Besides, molluscs do not have an adaptive immune response and vaccination strategies are not possible. Molluscs listed diseases under EU legislation are mainly protozoan parasites, that's why histological techniques are recognized for their diagnosis. However, molecular techniques are being increasingly used primarily as confirmatory techniques of the presence of the pathogens but also in disease monitoring programs. Research perspectives are mainly focussed in the optimization, of the already described techniques to gain in sensitivity and sensibility and in the development of new molecular biology techniques (quantitative real time PCRs), that are faster and easier to apply and that allow a positive diagnosis even in early stages of infection. However, molecular tools detect DNA sequences of the pathogen which does not imply that pathogen is viable in the cell host and the infection is established. Consequently, it needs to be validated against other techniques, such as histology or in situ hybridization, so that its reliability can be determined.

А new Gal/GalNAc-specific lectin from the mussel Mytilus trossulus: Structure, tissue specificity, antimicrobial and antifungal activity

In the present study, a new Gal/GalNAc specific lectin from the mussel Mytilus trossulus (designated as MTL) was identified, and its expression levels, both in tissues and toward pathogen stimulation, were then characterized. The MTL primary structure was determined via cDNA sequencing. Deduced sequence of 150 amino acid residues showed 89% similarity to lectins from the mussels Crenomytilus grayanus and Mytilus galloprovincialis that were the first members of a new family of zoolectins. The results indicated that the MTL might be involved in immune response toward pathogen infection, and it might perform different recognition specificity toward bacteria or fungi.

Clams sensitivity towards As and Hg: A comprehensive assessment of native and exotic species

To assess the environmental impact of As and Hg, bioindicator organisms such as bivalves have been used. Nevertheless, few studies have assessed the impacts of As and Hg in Ruditapes decussatus and Ruditapes philippinarum, which are native and exotic species in Europe, respectively. The main goal of the present study was to assess elements' partitioning and detoxification strategies of R. decussatus and R. philippinarum. Both clams showed a higher capacity to bioconcentrate Hg (BCF 2.29-7.49), when compared to As (0.59-1.09). Furthermore, As accumulation in both species was similar in the soluble and insoluble fractions, while in both species the majority of Hg was found in the insoluble fraction. Clams exposed to As showed different detoxification strategies, since R. decussatus had higher ability to enhance antioxidant enzymes and metallothioneins in order to reduce toxicity, and R.philippinarum increased glutathione S-transferase Ω activity, that catalyzes monomethyl arsenate reduction, the rate-limiting reaction in arsenic biotransformation. When exposed to Hg, R. decussatus presented, higher synthesis of antioxidant enzymes and lower LPO, being able to better tolerate Hg than the exotic species R. philippinarum. Thus under relevant levels of As and Hg contamination our work evidenced the higher ability of R. decussatus to survive and inhabit coastal environments not heavily contaminated by Hg and As.

Clam focal and systemic immune responses to QPX infection revealed by RNA-seq technology

Background

The hard clam Mercenaria mercenaria is an important seafood species widely exploited along the eastern coasts of the United States and play a crucial role in coastal ecology and economy. Severe hard clam mortalities have been associated with the protistan parasite QPX (Quahog Parasite Unknown). QPX infection establishes in pallial organs with the lesions typically characterized as nodules, which represent inflammatory masses formed by hemocyte infiltration and encapsulation of parasites. QPX infection is known to induce host changes on both the whole-organism level and at specific lesion areas, which imply systemic and focal defense responses, respectively. However, little is known about the molecular mechanisms underlying these alterations.

Results

RNA-seq was performed using Illumina Hiseq 2000 (641 Million 100 bp reads) to characterize M. mercenaria focal and systemic immune responses to QPX. Transcripts were assembled and the expression levels were compared between nodule and healthy tissues from infected clams, and between these and tissues from healthy clams. De novo assembly reconstructed a consensus transcriptome of 62,980 sequences that was functionally-annotated. A total of 3,131 transcripts were identified as differentially expressed in different tissues. Results allowed the identification of host immune factors implicated in the systemic and focal responses against QPX and unraveled the pathways involved in parasite neutralization. Among transcripts significantly modulated upon host-pathogen interactions, those involved in non-self recognition, signal transduction and defense response were over-represented. Alterations in pathways regulating hemocyte focal adhesion, migration and apoptosis were also demonstrated.

Conclusions

Our study is the first attempt to thoroughly characterize M. mercenaria transcriptome and identify molecular features associated with QPX infection. It is also one of the first studies contrasting focal and systemic responses to infections in invertebrates using high-throughput sequencing. Results identified the molecular signatures of clam systemic and focal defense responses, to collectively mediate immune processes such as hemocyte recruitment and local inflammation. These investigations improve our understanding of bivalve immunity and provide molecular targets for probing the biological bases of clam resistance towards QPX.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2493-9) contains supplementary material, which is available to authorized users.

The complexity of apoptotic cell death in mollusks: An update

Apoptosis is a type of programmed cell death that produces changes in cell morphology and in biochemical intracellular processes without inflammatory reactions. The components of the apoptotic pathways are conserved throughout evolution. Caspases are key molecules involved in the transduction of the death signal and are responsible for many of the biochemical and morphological changes associated with apoptosis. Nowadays, It is known that caspases are activated through two major apoptotic pathways (the extrinsic or death receptor pathway and the intrinsic or mitochondrial pathway), but there are also evidences of at least other alternative pathway (the perforin/granzyme pathway). Apoptosis in mollusks seems to be similar in complexity to apoptosis in vertebrates but also has unique features maybe related to their recurrent exposure to environmental changes, pollutants, pathogens and also related to the sedentary nature of some stages in the life cycle of mollusks bivalves and gastropods. As in other animals, apoptotic process is involved in the maintenance of tissue homeostasis and also constitutes an important immune response that can be triggered by a variety of stimuli, including cytokines, hormones, toxic insults, viruses, and protozoan parasites. The main goal of this work is to present the current knowledge of the molecular mechanisms of apoptosis in mollusks and to highlight those steps that need further study.

2-DE Mapping of the Blue Mussel Gill Proteome: The Usual Suspects Revisited

The Blue Mussel (Mytilus edulis, L. 1758) is an ecologically important and commercially relevant bivalve. Because of its ability to bioconcentrate xenobiotics, it is also a widespread sentinel species for environmental pollution, which has been used in ecotoxicological studies for biomarker assessment. Consequently, numerous proteomics studies have been carried out in various research contexts using mussels of the genus Mytilus, which intended to improve our understanding of complex physiological processes related to reproduction, adaptation to physical stressors or shell formation and for biomarker discovery. Differential-display 2-DE proteomics relies on an extensive knowledge of the proteome with as many proteoforms identified as possible. To this end, extensive characterization of proteins was performed in order to increase our knowledge of the Mytilus gill proteome. On average, 700 spots were detected on 2-DE gels by colloidal blue staining, of which 122 different, non-redundant proteins comprising 203 proteoforms could be identified by tandem mass spectrometry. These proteins could be attributed to four major categories: (i) “metabolism”, including antioxidant defence and degradation of xenobiotics; (ii) “genetic information processing”, comprising transcription and translation as well as folding, sorting, repair and degradation; (iii) “cellular processes”, such as cell motility, transport and catabolism; (iv) “environmental information processing”, including signal transduction and signalling molecules and interaction. The role of cytoskeleton proteins, energetic metabolism, chaperones/stress proteins, protein trafficking and the proteasome are discussed in the light of the exigencies of the intertidal environment, leading to an enhanced stress response, as well as the structural and physiological particularities of the bivalve gill tissue.

Identification, expression, and responses to bacterial challenge of the cathepsin C gene from the razor clam Sinonovacula constricta

Cathepsin C (dipeptidyl-peptidaseI, DPPI) is a lysosomal cysteine proteinase that belongs to the papain superfamily, and it is involved in protein degradation and proenzyme activation. However, very little is known about the function of cathepsin C in bivalves. In the present study, we identified the cathepsin C gene in the razor clam Sinonovacula constricta (Sc-CTSC). The full-length Sc-CTSC cDNA contained a complete open reading frame (ORF) of 1371 nt encoding 456 amino acids, a 98 bp 5' UTR, and a 1043 bp 3' UTR. The ORF of Sc-CTSC consisted of a putative signal peptide of 22 aa, a propeptide of 229 aa, and a mature peptide of 205 aa containing the active site triad of Cys, His, and Asn. The Sc-CTSC transcript was expressed in a wide range of tissues but exhibited the greatest level of expression in the digestive gland. During the early developmental stages, the transcript was detected widely. Upon injection with Vibrio anguillarum, the Sc-CTSC transcript was significantly up-regulated in digestive gland, mantle, and gill tissues. The results provided important information for further exploring the roles of cathepsin C in the innate immune responses.

A galectin with quadruple-domain from red abalone Haliotis rufescens involved in the immune innate response against to Vibrio anguillarum

Galectins are proteins that recognize and bind specifically β-galactosidase residues, playing important roles in the innate immune response of vertebrates and invertebrates. The cDNA of a tandem repeat galectin from the red abalone Haliotis rufescens cDNA (HrGal) was cloned and characterized using rapid amplification of cDNA end technique. The full-length cDNA of HrGal was 2471 bp, with a 5' terminal untranslated region (UTR) of 131 bp, a 3' UTR of 672 pb, and an open reading frame (ORF) of 1668 bp encoding a polypeptide of 556 amino acid. The ORF contains four domains carbohydrate recognition (CRD) with typical conserved motifs, which are important for carbohydrate recognition, and it appear to posses neither a signal peptide nor a transmembrane domain. The deduced amino acid sequence and the multi-domain organization of HrGal were highly similar to those described for other tandem repeat galectins of invertebrate organisms. Quantitative real time PCR analyses indicated that HrGal mRNA was highly expressed in hemocytes and gills tissues. The temporal expression of HrGal mRNA in hemocytes challenged to Vibrio anguillarum was time-dependent, showing u-regulation at 32 h post challenge. The results suggest that HrGal may be involved in the immune innate response against bacterial infection.

Galatrox is a C-type lectin in Bothrops atrox snake venom that selectively binds LacNAc-terminated glycans and can induce acute inflammation

Previous studies indicate that snake venom contains glycan-binding proteins (GBPs), although the binding specificity and biological activities of many of these GBPs is unclear. Here we report our studies on the glycan binding specificity and activities of galatrox, a Bothrops atrox snake venom-derived GBP. Glycan microarray analysis indicates that galatrox binds most strongly to glycans expressing N-acetyllactosamine (LacNAc), with a significant preference for Galβ1-4GlcNAcβ over Galβ1-3GlcNAcβ compounds. Galatrox also bound immobilized laminin, a LacNAc-dense extracellular matrix component, suggesting that this GBP can bind LacNAc-bearing glycoproteins. As several endogenous mammalian GBPs utilize a similar binding LacNAc binding preference to regulate neutrophil and monocyte activity, we hypothesized that galatrox may mediate B. atrox toxicity through regulation of leukocyte activity. Indeed, galatrox bound neutrophils and promoted leukocyte chemotaxis in a carbohydrate-dependent manner. Similarly, galatrox administration into the mouse peritoneal cavity induced significant neutrophil migration and the release of pro-inflammatory cytokines IL-1α and IL-6. Exposure of bone marrow-derived macrophages to galatrox induced generation of pro-inflammatory mediators IL-6, TNF-α, and keratinocyte-derived chemokine. This signaling by galatrox was mediated via its carbohydrate recognition domain by activation of the TLR4-mediated MyD88-dependent signaling pathway. These results indicate that galatrox has pro-inflammatory activity through its interaction with LacNAc-bearing glycans on neutrophils, macrophages and extracellular matrix proteins and induce the release of pro-inflammatory mediators.

Gene expression profile analysis of Manila clam (Ruditapes philippinarum) hemocytes after a Vibrio alginolyticus challenge using an immune-enriched oligo-microarray

Background

The Manila clam (Ruditapes philippinarum) is a cultured bivalve with worldwide commercial importance, and diseases cause high economic losses. For this reason, interest in the immune genes in this species has recently increased. The present work describes the construction of the first R. philippinarum microarray containing immune-related hemocyte sequences and its application to study the gene transcription profiles of hemocytes from clams infected with V. alginolyticus through a time course.

Results

The complete set of sequences from R. philippinarum available in the public databases and the hemocyte sequences enriched in immune transcripts were assembled successfully. A total of 12,156 annotated sequences were used to construct the 8 ×15 k oligo-microarray. The microarray experiments yielded a total of 579 differentially expressed transcripts. Using the gene expression results, the associated Gene Ontology terms and the enrichment analysis, we found different response mechanisms throughout the experiment. Genes related to signaling, transcription and apoptosis, such as IL-17D, NF-κB or calmodulin, were typically expressed as early as 3 hours post-challenge (hpc), while characteristic immune genes, such as PGRPs, FREPs and defense proteins appeared later at 8 hpc. This immune-triggering response could have affected a high number of processes that seemed to be activated 24 hpc to overcome the Vibrio challenge, including the expression of many cytoskeleton molecules, which is indicative of the active movement of hemocytes. In fact functional studies showed an increment in apoptosis, necrosis or cell migration after the infection. Finally, 72 hpc, activity returned to normal levels, and more than 50% of the genes were downregulated in a negative feedback of all of the previously active processes.

Conclusions

Using a new version of the R. philippinarum oligo-microarray, a putative timing for the response against a Vibrio infection was established. The key point to overcome the challenge seemed to be 8 hours after the challenge, when we detected immune functions that could lead to the destruction of the pathogen and the activation of a wide variety of processes related to homeostasis and defense. These results highlight the importance of a fast response in bivalves and the effectiveness of their innate immune system.

A microarray-based analysis of gametogenesis in two Portuguese populations of the European clam Ruditapes decussatus

The European clam, Ruditapes decussatus is a species with a high commercial importance in Portugal and other Southern European countries. Its production is almost exclusively based on natural recruitment, which is subject to high annual fluctuations. Increased knowledge of the natural reproductive cycle of R. decussatus and its molecular mechanisms would be particularly important in providing new highly valuable genomic information for better understanding the regulation of reproduction in this economically important aquaculture species. In this study, the transcriptomic bases of R. decussatus reproduction have been analysed using a custom oligonucleotide microarray representing 51,678 assembled contigs. Microarray analyses were performed in four gonadal maturation stages from two different Portuguese wild populations, characterized by different responses to spawning induction when used as progenitors in hatchery. A comparison between the two populations elucidated a specific pathway involved in the recognition signals and binding between the oocyte and components of the sperm plasma membrane. We suggest that this pathway can explain part of the differences in terms of spawning induction success between the two populations. In addition, sexes and reproductive stages were compared and a correlation between mRNA levels and gonadal area was investigated. The lists of differentially expressed genes revealed that sex explains most of the variance in gonadal gene expression. Additionally, genes like Foxl2, vitellogenin, condensing 2, mitotic apparatus protein p62, Cep57, sperm associated antigens 6, 16 and 17, motile sperm domain containing protein 2, sperm surface protein Sp17, sperm flagellar proteins 1 and 2 and dpy-30, were identified as being correlated with the gonad area and therefore supposedly with the number and/or the size of the gametes produced.

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.

Genes of the mitochondrial apoptotic pathway in Mytilus galloprovincialis

Bivalves play vital roles in marine, brackish, freshwater and terrestrial habitats. In recent years, these ecosystems have become affected through anthropogenic activities. The ecological success of marine bivalves is based on the ability to modify their physiological functions in response to environmental changes. One of the most important mechanisms involved in adaptive responses to environmental and biological stresses is apoptosis, which has been scarcely studied in mollusks, although the final consequence of this process, DNA fragmentation, has been frequently used for pollution monitoring. Environmental stressors induce apoptosis in molluscan cells via an intrinsic pathway. Many of the proteins involved in vertebrate apoptosis have been recognized in model invertebrates; however, this process might not be universally conserved. Mytilus galloprovincialis is presented here as a new model to study the linkage between molecular mechanisms that mediate apoptosis and marine bivalve ecological adaptations. Therefore, it is strictly necessary to identify the key elements involved in bivalve apoptosis. In the present study, six mitochondrial apoptotic-related genes were characterized, and their gene expression profiles following UV irradiation were evaluated. This is the first step for the development of potential biomarkers to assess the biological responses of marine organisms to stress. The results confirmed that apoptosis and, more specifically, the expression of the genes involved in this process can be used to assess the biological responses of marine organisms to stress.

Ruditapes decussatus and Ruditapes philippinarum exposed to cadmium: toxicological effects and bioaccumulation patterns

Since differences in metal accumulation may exist between bivalve species, the aim of this study was to assess the impact of cadmium (Cd) on Ruditapes decussatus and Ruditapes philippinarum. For this, the Cd accumulation, mortality rates and biochemical responses were analysed in the two species after 5days of exposure, under laboratory-controlled conditions. The concentration of Cd that caused 50% of mortality on clams was two-times higher in R. decussatus than in R. philippinarum. For both species, higher percentage (84.5-98.2%) of the Cd was in the insoluble fraction, but the Cd concentration in solution was 3 to 8 times higher in R. decussatus. Nevertheless, R. philippinarum presented higher oxidative stress and higher CAT activity. The paradox observed between the two clams can be explained by the higher capacity of R. decussatus to increase the expression of MTs when exposed to Cd.

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.

First molluscan theta-class Glutathione S-Transferase: identification, cloning, characterization and transcriptional analysis post immune challenges

Glutathione S-Transferases (GSTs) are multifunctional cytosolic isoenzymes, distinctly known as phase II detoxification enzymes. GSTs play a significant role in cellular defense against toxicity and have been identified in nearly all organisms studied to date, from bacteria to mammals. In this study, we have identified a full-length cDNA of the theta class GST from Ruditapes philippinarum (RpGSTθ), an important commercial edible molluscan species. RpGSTθ was cloned and the recombinant protein expressed, in order to study its biochemical characteristics and determine its physiological activities. The cDNA comprised an ORF of 693 bp, encoding 231 amino acids with a predicted molecular mass of 27 kDa and an isoelectric point of 8.2. Sequence analysis revealed that RpGSTθ possessed characteristic conserved domains of the GST_N family, Class Theta subfamily (PSSM: cd03050) and GST_C_family Super family (PSSM: cl02776). Phylogenetic analysis showed that RpGSTθ evolutionarily linked with other theta class homologues. The recombinant protein was expressed in Escherichia coli BL21(DE3) cells and the purified enzyme showed high activity with GST substrates like CDNB and 4-NBC. Glutathione dependent peroxidase activity of GST, investigated with cumene hydroperoxide as substrate affirmed the antioxidant property of rRpGSTθ. By quantitative PCR, RpGSTθ was found to be ubiquitously expressed in all tissues examined, with the highest levels occurring in gills, mantle, and hemocytes. Since GSTs may act as detoxification enzymes to mediate immune defense, the effects of pathogen associated molecular pattern, lipopolysaccharide and intact Vibrio tapetis bacteria challenge on RpGSTθ gene transcription were studied. Furthermore, the RpGSTθ expression changes induced by immune challenges were similar to those of the antioxidant defense enzyme manganese superoxide dismutase (RpMnSOD). To our knowledge, RpGSTθ is the first molluscan theta class GST reported, and its immune-related role in Manila clam may provide insights into potential therapeutic targets for protecting this important aquaculture species.