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

Unveiling the impact of bisphenol A on date mussels: Insights into oxidative stress, hormonal imbalance, gonadal atresia, and immune resilience

Sedentary organisms, such as mussels, may be susceptible to environmental estrogenic compounds, including bisphenol A (BPA). This study aimed to evaluate the interplay between BPA exposure and the immune response, hormonal imbalance, tissue damage (specifically in the digestive glands, labial palps, and male gonads), gonadal atresia, and antioxidant mechanisms in the marine mussel, Lithophaga lithophaga. Over a period of 28 days, mussels were exposed to BPA concentrations of 0, 0.25, 1, 2, and 5 μg/L. The exposure resulted in notable morphological alterations in the hemocytes of L. lithophaga, characterized by irregularities in the outer cell membranes of granulocytes and hyalinocytes, with some cells exhibiting filopodia formation. Granulocytes displayed an increased number of granules and vacuoles, while the nuclei of hyalinocytes appeared shrunken. The condition index, along with levels of testosterone and 17β-estradiol, significantly decreased with increasing BPA concentration, except for the 1 and 2 μg/L treatments. BPA exposure led to a marked increase in malondialdehyde (MDA) levels and a reduction in reduced glutathione (GSH) across all tissues at every concentration tested. The activity of antioxidant enzymes varied among the gonads, digestive glands, and labial palps. Notably, there was a significant increase in superoxide dismutase (SOD) activity in the gonads of mussels exposed to 2 μg/L of BPA, as well as in the digestive glands and labial palps of those exposed to 1 μg/L, suggesting a potential alteration in redox homeostasis. Additionally, structural changes in the digestive tubules of BPA-exposed mussels were observed. The observed pathological symptoms were characteristic of an inflammatory response, including hemocyte diapedesis and infiltration, the formation of syncytia, and the sloughing of epithelial tissue, indicated by an increased ratio of mean luminal radius to mean epithelial thickness in a dose-dependent manner. In the BPA-exposed group, testicular follicles exhibited atrophy, deformation, and a reduction in both size and number per area, appearing nearly empty and lacking spermatids and spermatozoa, alongside hypertrophy and hyperplasia of auxiliary cells. Scanning electron microscopy further revealed structural abnormalities in the heads and flagella of spermatids from the BPA-exposed group. Thus, this study demonstrates the risk of long-term exposure to BPA in immune response, tissue, and biochemical responses of date mussel L. lithophaga. The gonad was the most affected tissues followed by the digestive gland and labial palps.

2-Ethylhexyl-4-methoxycinnamate on marine and coastal environments: A comprehensive review of its environmental significance and biological impact

Marine and coastal environments are constantly subjected to increasing pressures associated with population growth, industrialization development, pollution and higher demand feeding society's consumerism. Among these pressures, there has been an increasing concern towards UV filters occurrence in aquatic ecosystems due to a greater use of personal care products (PCPs). 2-ethylhexyl-4-methoxycinnamate (EHMC) is one of the most used UV filters in sunscreen formulations, yet few reports address its effects in biota. This literature review intends to collect the available information concerning the environmental presence of EHMC in marine and coastal ecosystems and their effects in biota. The EHMC effects have been reported for the taxonomic groups: Actinomycetes, Alphaproteobacteria, Bacilli, Cytophagia, Flavobacteriia, Gammaproteobacteria, Actinopterygii, Anthozoa, Bacillariophyceae, Bivalvia, Branchiopoda, Coccolithophyceae, Echinoidea, Gastropoda, Malacostraca, Annelida and Thecostraca. The reported literature evaluated endpoints mainly related to development, viability, mortality, estrogenicity, gene transcription disruptions, biochemical alterations and morphophysiological changes. Based on the available information, there is still a clear need for further investigations related to EHMC and its toxicological effects on marine and coastal organisms.

Effect of salinity on the physiological response and transcriptome of spotted seabass (Lateolabrax maculatus)

Salinity fluctuations significantly impact the reproduction, growth, development, as well as physiological and metabolic activities of fish. To explore the osmoregulation mechanism of aquatic organisms acclimating to salinity stress, the physiological and transcriptomic characteristics of spotted seabass (Lateolabrax maculatus) in response to varying salinity gradients were investigated. In this study, different salinity stress exerted inhibitory effects on lipase activity, while the impact on amylase activity was not statistically significant. Notably, a moderate increase in salinity (24 psu) demonstrated the potential to enhance the efficient utilization of proteins by spotted seabass. Both Na+/K+-ATPase and malondialdehyde showed a fluctuating trend of increasing and then decreasing, peaking at 72 h. Transcriptomic analysis revealed that most differentially expressed genes were involved in energy metabolism, signal transduction, the immune response, and osmoregulation. These results will provide insights into the molecular mechanisms of salinity adaptation and contribute to sustainable development of the global aquaculture industry.

Expansion of the HSP70 gene family in Tegillarca granosa and expression profiles in response to zinc toxicity

Zinc (Zn) is an essential micronutrient in organisms and an abundant element in the Earth's crust. Trace amounts of Zn released from natural sources can enter aquatic ecosystems through weathering and erosion. Zn accumulates in organisms, and when its intracellular concentration exceeds a certain level, it can induce oxidative stress and trigger oxidative stress-mediated heat shock protein (HSP) modulation. HSP70 is the most evolutionarily conserved among the HSP families. Despite extensive research on HSP70 genes in bivalves, the HSP70 gene family of Tegillarca granosa is still poorly characterized. We identified 65 HSP70 genes belonging to 6 families in the T. granosa genome, with 50 HSPa12 and 11 HSPa B2 genes highly expanded. On chromosome 11, 39 HSP70 (60%) genes were identified, and the HSPa12A genes were highly duplicated. A total of 527 and 538 differentially expressed genes were identified in the gills and mantle based on Zn exposure, respectively. The Gene Ontology of cellular anatomical entities was significantly enriched with upregulated differentially expressed genes in the gills and mantle. Eight of the 11 HSPa B2 genes were upregulated in both tissues. Most of the genes identified in both tissues were involved in "protein homeostasis" and "inhibition of apoptosis," which are associated with the HSP70 family's resistance to extrinsic and intrinsic stress. Hence, this study identified that the HSP70 gene family plays a vital role in the adaptation of aquatic organisms to heavy metal (e.g., Zn) stress in contaminated environments by compiling the different physiological responses to preserve homeostasis.

A widespread picornavirus affects the hemocytes of the noble pen shell (Pinna nobilis), leading to its immunosuppression

Introduction

The widespread mass mortality of the noble pen shell (Pinna nobilis) has occurred in several Mediterranean countries in the past 7 years. Single-stranded RNA viruses affecting immune cells and leading to immune dysfunction have been widely reported in human and animal species. Here, we present data linking P. nobilis mass mortality events (MMEs) to hemocyte picornavirus (PV) infection. This study was performed on specimens from wild and captive populations.

Methods

We sampled P. nobilis from two regions of Spain [Catalonia (24 animals) and Murcia (four animals)] and one region in Italy [Venice (6 animals)]. Each of them were analyzed using transmission electron microscopy (TEM) to describe the morphology and self-assembly of virions. Illumina sequencing coupled to qPCR was performed to describe the identified virus and part of its genome.

Results and discussion

In 100% of our samples, ultrastructure revealed the presence of a virus (20 nm diameter) capable of replicating within granulocytes and hyalinocytes, leading to the accumulation of complex vesicles of different dimensions within the cytoplasm. As the PV infection progressed, dead hemocytes, infectious exosomes, and budding of extracellular vesicles were visible, along with endocytic vesicles entering other cells. The THC (total hemocyte count) values observed in both captive (eight animals) (3.5 × 104-1.60 × 105 ml-1 cells) and wild animals (14 samples) (1.90-2.42 × 105 ml-1 cells) were lower than those reported before MMEs. Sequencing of P. nobilis (six animals) hemocyte cDNA libraries revealed the presence of two main sequences of Picornavirales, family Marnaviridae. The highest number of reads belonged to animals that exhibited active replication phases and abundant viral particles from transmission electron microscopy (TEM) observations. These sequences correspond to the genus Sogarnavirus-a picornavirus identified in the marine diatom Chaetoceros tenuissimus (named C. tenuissimus RNA virus type II). Real-time PCR performed on the two most abundant RNA viruses previously identified by in silico analysis revealed positive results only for sequences similar to the C. tenuissimus RNA virus. These results may not conclusively identify picornavirus in noble pen shell hemocytes; therefore, further study is required. Our findings suggest that picornavirus infection likely causes immunosuppression, making individuals prone to opportunistic infections, which is a potential cause for the MMEs observed in the Mediterranean.

Invertebrate C1q Domain-Containing Proteins: Molecular Structure, Functional Properties and Biomedical Potential

C1q domain-containing proteins (C1qDC proteins) unexpectedly turned out to be widespread molecules among a variety of invertebrates, despite their lack of an integral complement system. Despite the wide distribution in the genomes of various invertebrates, data on the structure and properties of the isolated and characterized C1qDC proteins, which belong to the C1q/TNF superfamily, are sporadic, although they hold great practical potential for the creation of new biotechnologies. This review not only summarizes the current data on the properties of already-isolated or bioengineered C1qDC proteins but also projects further strategies for their study and biomedical application. It has been shown that further broad study of the carbohydrate specificity of the proteins can provide great opportunities, since for many of them only interactions with pathogen-associated molecular patterns (PAMPs) was evaluated and their antimicrobial, antiviral, and fungicidal activities were studied. However, data on the properties of C1qDC proteins, which researchers originally discovered as lectins and therefore studied their fine carbohydrate specificity and antitumor activity, intriguingly show the great potential of this family of proteins for the creation of targeted drug delivery systems, vaccines, and clinical assays for the differential diagnosis of cancer. The ability of invertebrate C1qDC proteins to recognize patterns of aberrant glycosylation of human cell surfaces and interact with mammalian immunoglobulins indicates the great biomedical potential of these molecules.

Molecular Features Associated with Resilience to Ocean Acidification in the Northern Quahog, Mercenaria mercenaria

The increasing concentration of CO₂ in the atmosphere and resulting flux into the oceans will further exacerbate acidification already threatening coastal marine ecosystems. The subsequent alterations in carbonate chemistry can have deleterious impacts on many economically and ecologically important species including the northern quahog (Mercenaria mercenaria). The accelerated pace of these changes requires an understanding of how or if species and populations will be able to acclimate or adapt to such swift environmental alterations. Thus far, studies have primarily focused on the physiological effects of ocean acidification (OA) on M. mercenaria, including reductions in growth and survival. However, the molecular mechanisms of resilience to OA in this species remains unclear. Clam gametes were fertilized under normal pCO₂ and reared under acidified (pH ~ 7.5, pCO₂ ~ 1200 ppm) or control (pH ~ 7.9, pCO₂ ~ 600 ppm) conditions before sampled at 2 days (larvae), 32 days (postsets), 5 and 10 months (juveniles) and submitted to RNA and DNA sequencing to evaluate alterations in gene expression and genetic variations. Results showed significant shift in gene expression profiles among clams reared in acidified conditions as compared to their respective controls. At 10 months of exposure, significant shifts in allele frequency of single nucleotide polymorphisms (SNPs) were identified. Both approaches highlighted genes coding for proteins related to shell formation, bicarbonate transport, cytoskeleton, immunity/stress, and metabolism, illustrating the role these pathways play in resilience to OA.

RNA-Seq comparative study reveals molecular effectors linked to the resistance of Pinna nobilis to Haplosporidium pinnae parasite

With the intensification of maritime traffic, recently emerged infectious diseases have become major drivers in the decline and extinction of species. Since 2016, mass mortality events have decimated the endemic Mediterranean Sea bivalve Pinna nobilis, affecting ca. 100% of individuals. These events have largely been driven by Haplosporidium pinnae's infection, an invasive species which was likely introduced by shipping. While monitoring wild populations of P. nobilis, we observed individuals that survived such a mass mortality event during the summer of 2018 (France). We considered these individuals resistant, as they did not show any symptoms of the disease, while the rest of the population in the area was devastated. Furthermore, the parasite was not detected when we conducted a PCR amplification of a species-specific fragment of the small subunit ribosomal DNA. In parallel, the transcriptomic analysis showed evidence of some parasite RNA indicating that the resistant individuals had been exposed to the parasite without proliferating. To understand the underlying mechanisms of resistance in these individuals, we compared their gene expression with that of susceptible individuals. We performed de novo transcriptome assembly and annotated the expressed genes. A comparison of the transcriptomes in resistant and susceptible individuals highlighted a gene expression signature of the resistant phenotype. We found significant differential expressions of genes involved in immunity and cell architecture. This data provides the first insights into how individuals escape the pathogenicity associated with infection.

Comparative transcriptome analysis reveals immunotoxicology induced by three organic UV filters in Manila clam (Ruditapes philippinarum)

Benzophenone-3 (BP-3), 4-methyl-benzylidene camphor (4-MBC) and 2-ethyl-hexyl-4-trimethoxycinnamate (EHMC) are commonly used organic ultraviolet (UV) filters and are frequently detected in water environments. In the present study, we studied the potential adverse impacts of UV filter exposures in Ruditapes philippinarum by investigating transcriptomic profiles and non-specific immune enzyme activities. Transcriptome analysis showed that more genes were differentially regulated in EHMC-treated group, and down-regulated genes (2009) were significantly more than up-regulated ones (410) at day 7. Function annotation revealed that pathways "immune system", "cell growth and death" and "infectious diseases" were significantly enriched. Generally, combined qPCR and biochemical analyses demonstrated that short-term exposure to low dose of UV filters could activate immune responses, whereas the immune system would be restrained after prolonged exposure. Taken together, the present study firstly demonstrated the immunotoxicology induced by BP-3, 4-MBC and EHMC on R. philippinarum, indicating their potential threats to the survival of marine bivalves.

Effects of low salinity on hemolymph osmolality and transcriptome of the Iwagaki oyster, Crassostrea nippona

Crassostrea nippona is a kind of oysters with great development value as it can be edible in summer for its late reproductive period. Salinity is one of the important limiting abiotic factors to the survival and distribution of this stenohaline species. To better understand the physiological and immunological response of C. nippona to varying environmental salinities, the effects of low salinity on the hemolymph osmolality and gill transcriptome were investigated in this study. The osmolality of hemolymph in vivo and surrounding water were assessed regularly over one week at five test salinities ranging from 5 psμ to 30 psμ. They reached osmotic equilibrium within hours above 15 psμ but remained hyperosmotic at 10 and 5 psμ for the whole sampling period. Through comparative transcriptome analysis, there were less differentially expressed genes (DEGs) in pairwise comparison of S1 (10 psμ) vs S3 (30 psμ) than in S2 (20 psμ) vs S3. KEGG enrichment analysis identified ubiquitin-mediated proteolysis and mitochondrial apoptosis pathway specifically enriched at 10 psμ. This study gained comprehensive insights on the low salinity response of C. nippona at the molecular level, which provide a theoretical basis for understanding the immune mechanism under low salinity stress.

Apextrin from Ruditapes philippinarum functions as pattern recognition receptor and modulates NF-κB pathway

Apextrin belongs to ApeC-containing proteins (ACPs) and features a signal-peptide, an N-terminal membrane attack complex component/perforin (MACPF) domain, and a C-terminal ApeC domain. Recently, apextrin-like proteins were identified as pattern recognition receptor (PRR), which recognize the bacterial cell wall component and participate in innate immunity. Here, an apextrin (Rpape) was identified and characterized in Ruditapes philippinarum. Our results showed that Rpape mRNA was significantly induced under bacterial challenges. The Rpape recombinant protein exhibited a significant inhibitory effect on gram-positive bacteria (Bacillus subtilis and Staphylococcus aureus) and bound with Vibrio anguillarum, S. aureus and B. subtilis. We found Rpape protein positively activated the NF-κB signaling cascade and increased the activity of Nitric oxide (NO). This study revealed the immunity role of apextrin in R. philippinarum and provided a reference for further study on the role of apextrin in bivalves.

RNA Interference-Mediated Knockdown of Bombyx mori Haemocyte-Specific Cathepsin L (Cat L)-Like Cysteine Protease Gene Increases Bacillus thuringiensis kurstaki Toxicity and Reproduction in Insect Cadavers

The silkworm’s Cat L-like gene, which encodes a lysosomal cathepsin L-like cysteine protease, is thought to be part of the insect’s innate immunity via an as-yet-undetermined mechanism. Assuming that the primary function of Cat L-like is microbial degradation in mature phagosomes, we hypothesise that the suppression of the Cat L-like gene expression would increase Bacillus thuringiensis (Bt) bacteraemia and toxicity in knockdown insects. Here, we performed a functional analysis of Cat L-like in larvae that were fed mulberry leaves contaminated with a commercial biopesticide formulation based on Bt kurstaki (Btk) (i.e., Dipel) to investigate its role in insect defence against a known entomopathogen. Exposure to sublethal doses of Dipel resulted in overexpression of the Cat L-like gene in insect haemolymph 24 and 48 h after exposure. RNA interference (RNAi)-mediated suppression of Cat L-like expression significantly increased the toxicity of Dipel to exposed larvae. Moreover, Btk replication was higher in RNAi insects, suggesting that Cat L-like cathepsin may be involved in a bacterial killing mechanism of haemocytes. Finally, our results confirm that Cat L-like protease is part of the antimicrobial defence of insects and suggest that it could be used as a target to increase the insecticidal efficacy of Bt-based biopesticides.

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

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

Cracking the Challenge of Antimicrobial Drug Resistance with CRISPR/Cas9, Nanotechnology and Other Strategies in ESKAPE Pathogens

Antimicrobial resistance is mushrooming as a silent pandemic. It is considered among the most common priority areas identified by both national and international agencies. The global development of multidrug-resistant strains now threatens public health care improvement by introducing antibiotics against infectious agents. These strains are the product of both continuous evolution and unchecked antimicrobial usage (AMU). The ESKAPE pathogens (Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa, and Enterobacter species) are the leading cause of nosocomial infections throughout the world. Most of them are now multidrug-resistant, which pose significant challenges in clinical practice. Understanding these bacteria’s resistance mechanisms is crucial for developing novel antimicrobial agents or other alternative tools to fight against these pathogens. A mechanistic understanding of resistance in these pathogens would also help predict underlying or even unknown mechanisms of resistance of other emerging multidrug-resistant pathogens. Research and development to find better antibacterial drugs and research on tools like CRISPER-Cas9, vaccines, and nanoparticles for treatment of infections that can be further explored in the clinical practice health sector have recognized these alternatives as essential and highly effective tools to mitigate antimicrobial resistance. This review summarizes the known antimicrobial resistance mechanisms of ESKAPE pathogens and strategies for overcoming this resistance with an extensive overview of efforts made in this research area.

In vitro Evaluation of Programmed Cell Death in the Immune System of Pacific Oyster Crassostrea gigas by the Effect of Marine Toxins

Programmed cell death (PCD) is an essential process for the immune system's development and homeostasis, enabling the remotion of infected or unnecessary cells. There are several PCD's types, depending on the molecular mechanisms, such as non-inflammatory or pro-inflammatory. Hemocytes are the main component of cellular immunity in bivalve mollusks. Numerous infectious microorganisms produce toxins that impair hemocytes functions, but there is little knowledge on the role of PCD in these cells. This study aims to evaluate in vitro whether marine toxins induce a particular type of PCD in hemocytes of the bivalve mollusk Crassostrea gigas during 4 h at 25°C. Hemocytes were incubated with two types of marine toxins: non-proteinaceous toxins from microalgae (saxitoxin, STX; gonyautoxins 2 and 3, GTX2/3; okadaic acid/dynophysistoxin-1, OA/DTX-1; brevetoxins 2 and 3, PbTx-2,-3; brevetoxin 2, PbTx-2), and proteinaceous extracts from bacteria (Vibrio parahaemolyticus, Vp; V. campbellii, Vc). Also, we used the apoptosis inducers, staurosporine (STP), and camptothecin (CPT). STP, CPT, STX, and GTX 2/3, provoked high hemocyte mortality characterized by apoptosis hallmarks such as phosphatidylserine translocation into the outer leaflet of the cell membrane, exacerbated chromatin condensation, DNA oligonucleosomal fragments, and variation in gene expression levels of apoptotic caspases 2, 3, 7, and 8. The mixture of PbTx-2,-3 also showed many apoptosis features; however, they did not show apoptotic DNA oligonucleosomal fragments. Likewise, PbTx-2, OA/DTX-1, and proteinaceous extracts from bacteria Vp, and Vc, induced a minor degree of cell death with high gene expression of the pro-inflammatory initiator caspase-1, which could indicate a process of pyroptosis-like PCD. Hemocytes could carry out both PCD types simultaneously. Therefore, marine toxins trigger PCD's signaling pathways in C. gigas hemocytes, depending on the toxin's nature, which appears to be highly conserved both structurally and functionally.

Phylogenetic analysis of the caspase family in bivalves: implications for programmed cell death, immune response and development

BACKGROUND

Apoptosis is an important process for an organism's innate immune system to respond to pathogens, while also allowing for cell differentiation and other essential life functions. Caspases are one of the key protease enzymes involved in the apoptotic process, however there is currently a very limited understanding of bivalve caspase diversity and function.

RESULTS

In this work, we investigated the presence of caspase homologues using a combination of bioinformatics and phylogenetic analyses. We blasted the Crassostrea gigas genome for caspase homologues and identified 35 potential homologues in the addition to the already cloned 23 bivalve caspases. As such, we present information about the phylogenetic relationship of all identified bivalve caspases in relation to their homology to well-established vertebrate and invertebrate caspases. Our results reveal unexpected novelty and complexity in the bivalve caspase family. Notably, we were unable to identify direct homologues to the initiator caspase-9, a key-caspase in the vertebrate apoptotic pathway, inflammatory caspases (caspase-1, - 4 or - 5) or executioner caspases-3, - 6, - 7. We also explored the fact that bivalves appear to possess several unique homologues to the initiator caspase groups - 2 and - 8. Large expansions of caspase-3 like homologues (caspase-3A-C), caspase-3/7 group and caspase-3/7-like homologues were also identified, suggesting unusual roles of caspases with direct implications for our understanding of immune response in relation to common bivalve diseases. Furthermore, we assessed the gene expression of two initiator (Cg2A, Cg8B) and four executioner caspases (Cg3A, Cg3B, Cg3C, Cg3/7) in C. gigas late-larval development and during metamorphosis, indicating that caspase expression varies across the different developmental stages.

CONCLUSION

Our analysis provides the first overview of caspases across different bivalve species with essential new insights into caspase diversity, knowledge that can be used for further investigations into immune response to pathogens or regulation of developmental processes.

Molecular cloning and expression analysis of mc5r like genes (mc5rl) in Ruditapes philippinarum (Manila clam) after aerial exposure and low-temperature stress

The melanocortin-5 receptor (mc5r) plays an important role in exocrine function, lipid metabolism, obesity, and stress response in the vertebrate. However, the functions of the mc5r in mollusks have been rarely investigated. We cloned the full length of Ruditapes philippinarum mc5r like gene (mc5rl) and the sequence structure and phylogenetic relationship of mc5rl were analyzed. Besides, we detected the tissue distribution and the expression pattern of R. philippinarum mc5r like (mc5rl) genes after aerial exposure and low-temperature stress. The full-length cDNA of the mc5rl-1 was 2143 bp, consisting of a 1224 bp open reading frame encoding (ORF) 408 amino acids. Sequence and phylogenetic analyses revealed that the nucleotide and amino acid sequences of Manila clam mc5rl were highly homologous with mc5r of Crassostrea virginica, Crassostrea gigas, Mizuhopecten yessoensis, and Pecten maximus (32%-36%) and low homologous with vertebrates. The results of the distribution of mc5rl genes showed that mc5rl genes were dominant in the mantle, gonad, and hepatopancreas in R. philippinarum. The expression of mc5rl genes was significantly increased after aerial exposure and low-temperature stress in R. philippinarum in hepatopancreas. Aerial exposure and low-temperature stress could induce mc5rl expressed. Mc5rl might serve as a sensor and promote stress response in R. philippinarum. The cloning and expression characteristics of mc5rl will facilitate the investigation of its function in stress response and other physiological processes in R. philippinarum.

Transcriptomic analysis of Manila clam Ruditapes philippinarum under lipopolysaccharide challenge provides molecular insights into immune response

The Manila clam, Ruditapes philippinarum, is an economically important shellfish in marine aquaculture. A better understanding of the immune system in R. philippinarum will provide the basis for the development of strategies to mitigate the impact of infectious diseases affecting this species but can also be of relevance for other bivalves of commercial interest. In this study, the transcriptional response of the Manila clam under lipopolysaccharide (LPS) challenge was characterized using RNA sequencing. The transcriptomes of LPS challenged group of clams (LH1, LH2 and LH3), and the PBS control group (CH1, CH2 and CH3), were sequenced with the Illumina HiSeq platform. Compared with the unigene expression profile of the control group, 223 unigenes were up-regulated and 389 unigenes were down-regulated in the LPS challenged group. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses revealed that signal transduction, defense response, and immune-related pathways such as Chemokine signaling pathway, Complement and coagulation cascades, NOD-like receptor signaling pathway, and Inflammatory mediator regulation of TRP channels in sensory system were the most highly enriched pathways among the genes that were differentially expressed under LPS challenge. This study present understanding of the molecular basis underpinning response to LPS challenge and provides useful information for future work on the molecular mechanism of pathogen resistance and immunity in Manila clam.

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.

Transcriptome analysis and histopathological observations of Geloina erosa gills upon Cr(VI) exposure

The heavy metal contamination like Cr(VI) has been increased by human activities and that threats the ecosystem health of mangrove areas. Bioindicator is an emerging tool in the environmental contamination assessment. The objective of this study was to investigate the Geloina erosa response mechanisms and sensitivities of several biomarkers in the Cr(VI) exposure and identify the G. erosa capability of being used as heavy metals bioindicator. In this study, G. erosa was exposed to 100 μmol·L^-1 Cr(VI) for 48 h. After transcriptome sequencing, a total of 134,817 unigenes were obtained, including 12,555 up-regulated and 18,829 down-regulated differentially expressed genes and were validated through quantitative real-time PCR. In addition, a total of 12,185 SSRs and 1,428,214 candidate SNPs were identified from all the G. erosa transcriptome libraries. Histopathology of the gill indicated the Cr(VI) exposure induced damage of the organ leading to its immunization, detoxification or apoptosis reactions. Among eight genes of the selected biomarkers, Calm, HSP70, CYP450, ATG5, TLR2, MYD88 and CASP8 were up-regulated, while TLR4 was down-regulated in response to the Cr(VI) exposure.

Transcriptomic analysis of clam extrapallial fluids reveals immunity and cytoskeleton alterations in the first week of Brown Ring Disease development

The Brown Ring Disease is an infection caused by the bacterium Vibrio tapetis on the Manila clam Ruditapes philippinarum. The process of infection, in the extrapallial fluids (EPFs) of clams, involves alteration of immune functions, in particular on hemocytes which are the cells responsible of phagocytosis. Disorganization of the actin-cytoskeleton in infected clams is a part of what leads to this alteration. This study is the first transcriptomic approach based on collection of extrapallial fluids on living animals experimentally infected by V. tapetis. We performed differential gene expression analysis of EPFs in two experimental treatments (healthy-against infected-clams by V. tapetis), and showed the deregulation of 135 genes. In infected clams, a downregulation of transcripts implied in immune functions (lysosomal activity and complement- and lectin-dependent PRR pathways) was observed during infection. We also showed a deregulation of transcripts encoding proteins involved in the actin cytoskeleton organization such as an overexpression of β12-Thymosin (which is an actin sequestration protein) or a downregulation of proteins that closely interact with capping proteins such as Coactosin, that counteract action of capping proteins, or Profilin. We validated these transcriptomic results by cellular physiological analyses that showed a decrease of the lysosome amounts and the disorganization of actin cytoskeleton in infected hemocytes.

Host-microbiota interactions shed light on mortality events in the striped venus clam Chamelea gallina

Mass mortalities due to disease outbreaks have recently affected a number of major taxa in marine ecosystems. Climate- and pollution-induced stress may compromise host immune defenses, increasing the risk of opportunistic diseases. Despite growing evidence that mass mortality events affecting marine species worldwide are strongly influenced by the interplay of numerous environmental factors, the reductionist approaches most frequently used to investigate these factors hindered the interpretation of these multifactorial pathologies. In this study, we propose a broader approach based on the combination of RNA-sequencing and 16S microbiota analyses to decipher the factors underlying mass mortality in the striped venus clam, Chamelea gallina, along the Adriatic coast. On one hand, gene expression profiling and functional analyses of microbial communities showed the over-expression of several genes and molecular pathways involved in xenobiotic metabolism, suggesting potential chemical contamination in mortality sites. On the other hand, the down-regulation of several genes involved in immune and stress response, and the over-representation of opportunistic pathogens such as Vibrio and Photobacterium spp. indicates that these microbial species may take advantage of compromised host immune pathways and defense mechanisms that are potentially affected by chemical exposure, resulting in periodic mortality events. We propose the application of our approach to interpret and anticipate the risks inherent in the combined effects of pollutants and microbes on marine animals in today's rapidly changing environment.

Immunity in mussels: An overview of molecular components and mechanisms with a focus on the functional defenses

Bivalves' immunity has received much more attention in the last decade, which resulted to a valuable growth in the availability of its molecular components. Such data availability coupled with the economical importance of these organisms aimed to shift the increase in the number of immunological and stress-related studies. Unfortunately, the crowd of generated data deciphering the involved physiological processes, investigators' differential conceptualization and the aimed objectives, has complicated the sensu stricto outlining of immune-related mechanisms. Overall, this review tried to compiles a summary about the molecular components of the mussels' immune response, surveying an overview of the mussels' functional immunity through gathering the most recent-related topics of bivalves' immunity as apoptosis and autophagy which deserves a great attention as stress-related mechanisms, the disseminated neoplasia as outbreak transmissible disease, not only within the same specie but also among different species, the hematopoiesis as topic that still generating interesting debate in the scientific community, the mucosal immunity described as the interface where host-pathogen interactions would occurs and determinate the late immune response, and innate immune memory and transgenerational priming, which described as very recent research topic with extensive applications in shellfish farming industry.

Understanding the transition from water to land: Insights from multi-omic analyses of the perivitelline fluid of apple snail eggs

Unlike most of the freshwater gastropod families, the family Ampullariidae includes members that exhibit both underwater and aerial oviposition, making it an ideal model for understanding mechanisms underlying the evolutionary transition from water to land. We applied SDS-PAGE and LC-MS/MS to analyse the proteome of the egg perivitelline fluid (PVF) of Marisa cornuarietis - an aquatic ovipositing ampullariid. Comparison with the reported PVF proteomes of two aerial ovipositing ampullariids (Pomacea canaliculata and P. maculata) showed that the three species all contain several major perivitellins that nourish the embryos. However, M. cornuarietis invests more heavily on immune-related proteins, which might be due to exposure to aquatic pathogens. Interestingly, only the PVF of out-of-water egg laying species have PV2 - a neurotoxin lethal to mice, and a calcium-binding protein which might be involved in the formation of calcareous eggshell. Integrated phylogenetic, evolutionary and gene expressional analyses detected the involvement of gene duplication, positive selection and neofunctionalisation in the formation of several major PVF proteins. Overall, our study provides multiple lines of evidence of adaptive evolution in the PVF proteins, and contributes to a better understanding of how aquatic gastropod ancestors invaded terrestrial habitats. SIGNIFICANCE: Aerial egg deposition has evolved in several groups of animals, but except for Vertebrata little is known about the mechanisms underlying this critical evolution process. We compared aquatic and aerial egg laying apple snails to understand the molecular mechanisms enabling such a transition in egg laying habitat. We found that the composition of perivitelline fluid proteomes of underwater and aerial egg depositors was remarkably different, and then gene duplication and positive selection were responsible for the formation of such novel proteins than enabled the evolutionary transition.

Comparative transcriptome analysis between the short-term stress and long-term adaptation of the Ruditapes philippinarum in response to benzo[a]pyrene

In order to monitor the pollution of polycyclic aromatic hydrocarbons (PAHs) in the seawater environment, screening biomarkers capable of monitoring PAHs is the focus of many studies. The transcriptomic profiles of the digestive gland tissue from the R. philippinarum groups after the exposure to BaP (4 μg/L) at four time points (0, 0.5, 6 and 15 days) were investigated to globally screen the key genes and pathways involved in the responses to short-term stress and long-term adaptation of BaP resistance. By comparative transcriptome analysis, 233, 282 and 58 differentially expressed genes (DEGs) were identified at 0.5 day, 6 day and 15 day (vs 0 day). The differential expression genes were related to stress response, detoxification metabolic process and innate immunity. DEGs of each group at different stages were clustered in six profiles based on gene expression pattern. Gene Ontology (GO) and Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathway analysis were used on all genes to determine the biological functions and processes. We selected Multidrug resistance protein 3 (MRP3), transcriptional regulator ATRX-like isoform X2 (ATRX) as biomarker indicator genes for short-term pollution monitoring and NADH dehydrogenase [ubiquinone] 1 (NQO1), Complement C1q-like protein 4 (C1q), Glutathione-S-transferase theta (GST), E3 ubiquitin-protein ligase (E3) for long-term pollution monitoring based on the different expression patterns and the function in detoxification and antioxidant defense system. Besides, the expression of seven genes was measured through Quantitative real-time PCR (qPCR) according to their gene expression patterns which was confirmed by the DGE analysis. Taken together, adoption of transcriptomic analysis to explore the bivalves' mRNA abundance changes and detoxification metabolic mechanism under the BaP stress at different time points can aid the development of sensitive and informed molecular endpoints for application towards ecotoxicogenomic monitoring of bivalves.

Molecular and cellular characterization of apoptosis in flat oyster a key mechanisms at the heart of host-parasite interactions

Bonamia ostreae has been associated with the decline of flat oyster Ostrea edulis populations in some European countries. This obligatory intracellular parasite persists and multiplies into hemocytes. Previous in vitro experiments showed that apoptosis is activated in hemocytes between 1 h and 4 h of contact with the parasite. The flat oyster uses the apoptosis pathway to defend against B. ostreae. However, the parasite might be also able to modulate this response in order to survive in its host. In order to investigate this hypothesis the apoptotic response of the host was evaluated using flow cytometry, transmission electron microscopy and by measuring the response of genes involved in the apoptotic pathway after 4 h. In parallel, the parasite response was investigated by measuring the expression of B. ostreae genes involved in different biological functions including cell cycle and cell death. Obtained results allow describing molecular apoptotic pathways in O. edulis and confirm that apoptosis is early activated in hemocytes after a contact with B. ostreae. Interestingly, at cellular and molecular levels this process appeared downregulated after 44 h of contact. Concurrently, parasite gene expression appeared reduced suggesting that the parasite could inhibit its own metabolism to escape the immune response.

Metabolomics Study of Immune Responses of New Zealand Greenshell™ Mussels (Perna canaliculus) Infected with Pathogenic Vibrio sp

Vibrio coralliilyticus is a bacterial pathogen which can affect a range of marine organisms, such as corals, fish and shellfish, with sometimes devastating consequences. However, little is known about the mechanisms involved in the host-pathogen interaction, especially within molluscan models. We applied gas chromatography-mass spectrometry (GC-MS)-based metabolomics to characterize the physiological responses in haemolymph of New Zealand Greenshell™ mussels (Perna canaliculus) injected with Vibrio sp. DO1 (V. coralliilyticus/neptunius-like isolate). Univariate data analyses of metabolite profiles in Vibrio-exposed mussels revealed significant changes in 22 metabolites at 6 h post-infection, compared to non-exposed mussels. Among them, 10 metabolites were up-regulated, while 12 metabolites were down-regulated in infected mussels. Multivariate analyses showed a clear distinction between infected and non-infected mussels. In addition, secondary pathway analyses indicated perturbations of the host innate immune system following infection, including oxidative stress, inflammation and disruption of the TCA cycle, change in amino acid metabolism and protein synthesis. These findings provide new insights into the pathogenic mechanisms of Vibrio infection of mussels and demonstrate our ability to detect detailed and rapid host responses from haemolymph samples using a metabolomics approach.

The immune response of Mytilus edulis hemocytes exposed to Vibrio splendidus LGP32 strain: A transcriptomic attempt at identifying molecular actors

The marine mussel Mytilus edulis, tolerant to a wide range of environmental changes, combines a key role as a sentinel species for environmental monitoring programs and a significant economic importance. Mortality events caused by infective agents and parasites have not been described in mussels, which suggests an efficient immune system. This study aims at identifying the molecular mechanisms involved in the early immune responses M. edulis' hemocytes challenged with Vibrio splendidus LGP32 strain during 2, 4 and 6 h. A total of 149,296 assembled sequences has been annotated and compared to KEGG reference pathways. Several immune related sequences were identified such as Toll-Like receptors (TLRs), transcription factors, cytokines, protease inhibitors, stress proteins and sequences encoding for proteins involved in cell adhesion, phagocytosis, oxidative stress, apoptosis and autophagy. Differential gene expression clustered 10 different groups of transcripts according to kinetics of transcript occurrence. Sequences were assigned to biological process gene ontology categories. Sequences encoding for galectins, fibrinogen-related proteins, TLRs, MyD88, some antimicrobial peptides, lysosomal hydrolases, heat shock proteins and protease inhibitors, as well as proteins of oxidative stress and apoptosis were identified as differently regulated during the exposure to V. splendidus LGP32. The levels of candidate transcripts were quantified in M. edulis' hemocytes exposed to V. splendidus LGP32 and 7SHRW by using branched DNA technology. Transcripts encoding for inhibitor kappa B, inhibitor of apoptosis proteins, tumor protein D54, serine/threonine-proteine kinase SIK2 were identified as up-regulated in hemocytes exposed to both strains.

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.

Molecular characterization and expression analysis of a novel r-spondin member (rspo2l) in Chinese tongue sole (Cynoglossus semilaevis)

Numerous studies suggest R-spondins (Rspos) plays a role in mammalian sex development and differentiation by activating WNT signaling pathways. However, Rspos are frequently less reported in teleosts. In this study, a molecular characterization and expression analysis was conducted with a new rspondin member in the Chinese tongue sole, rspondin2-like (rspo2l). The length of rspo2l cDNA is 1251 bp with 732 bp of coding sequence. A qRT-PCR analysis revealed that the transcription of rspo2l was distributed in various tissues, with high transcription levels in the liver, skin, and gills which might indicate a possible role in immunity. We next examined a time-course of transcription levels in four immune tissues (gill, liver, spleen, and kidney) after Vibrio harveyi challenge. It was found that rspo2l was up-regulated in the gills, spleen, and kidney and down-regulated in the liver, and the greatest responses occurred at 24 and 48 h after bacterial challenge. An assessment of β-catenin, the key regulator of the canonical WNT signaling pathway, at different time points in four immune organs revealed that its transcription profile was similar to that of rspo2l after bacterial challenge. The results suggest that tongue sole rspo2l might play a role in immune responses after bacterial challenge, while the potential link with the WNT signaling pathway still requires further investigation. This is the first report about the involvement of rspondins in fish immune responses.

The Whole-Genome and Transcriptome of the Manila Clam (Ruditapes philippinarum)

The manila clam, Ruditapes philippinarum, is an important bivalve species in worldwide aquaculture including Korea. The aquaculture production of R. philippinarum is under threat from diverse environmental factors including viruses, microorganisms, parasites, and water conditions with subsequently declining production. In spite of its importance as a marine resource, the reference genome of R. philippinarum for comprehensive genetic studies is largely unexplored. Here, we report the de novo whole-genome and transcriptome assembly of R. philippinarum across three different tissues (foot, gill, and adductor muscle), and provide the basic data for advanced studies in selective breeding and disease control in order to obtain successful aquaculture systems. An approximately 2.56 Gb high quality whole-genome was assembled with various library construction methods. A total of 108,034 protein coding gene models were predicted and repetitive elements including simple sequence repeats and noncoding RNAs were identified to further understanding of the genetic background of R. philippinarum for genomics-assisted breeding. Comparative analysis with the bivalve marine invertebrates uncover that the gene family related to complement C1q was enriched. Furthermore, we performed transcriptome analysis with three different tissues in order to support genome annotation and then identified 41,275 transcripts which were annotated. The R. philippinarum genome resource will markedly advance a wide range of potential genetic studies, a reference genome for comparative analysis of bivalve species and unraveling mechanisms of biological processes in molluscs. We believe that the R. philippinarum genome will serve as an initial platform for breeding better-quality clams using a genomic approach.

Transcriptome analysis of Mytilus coruscus hemocytes in response to Vibrio alginnolyficus infection

As an economically important bivalve, the Mytilus coruscus is cultured widely in the eastern coast of China. In recent years, this bivalve has been seriously affected by the pathogenic infections. To elucidate the host defense mechanisms of M. coruscus against pathogenic challenge, the hemocyte transcriptomes of M. coruscus before and after Vibrio alginnolyficus infection were analyzed using the deep-sequencing platform Illumina/HiSeq-2500, meanwhile the differentially expressed genes (DEGs) were investigated. In total, 130,031,083 clean reads were obtained and then assembled into 63,942 unigenes with an average length of 810 bp and an N50 of 1056 bp. Unigenes were annotated by comparing against nr, Swiss-Prot, KEGG, COG, KOG, GO, and Pfam databases, and 27,345 unigenes (42.77%) were annotated in at least one database. After bacterial challenge, 1270 and 265 genes were identified as remarkably up- or down-regulated, respectively, amongst 1154 were associated with 122 pathways, including classical immune-related pathways, such as 'Toll-like receptor signaling', 'the complement cascades', 'MAPK signaling pathway', 'Apoptosis' and 'Wnt signaling pathway'. Besides, nine genes which were differently-expressed immuno-related were confirmed by using quantitative real-time PCR. These findings would provide new insights on the M. coruscus innate immunity, based on which, some novel strategies for management of diseases and long-term sustainability of M. coruscus culture could be developed.

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.

Conservation and divergence of mitochondrial apoptosis pathway in the Pacific oyster, Crassostrea gigas

Apoptosis is considered a crucial part of the host defense system in oysters according to previous reports; however, the exact process by which this occurs remains unclear. Besides, mitochondrial apoptosis is the primary method of apoptosis in vertebrate cells, but has been poorly studied in invertebrates and is quite controversial. In this study, we investigated the molecular mechanism of mitochondrial apoptosis in the Pacific oyster Crassostrea gigas. Notably, we show that most key elements involved in the vertebrate mitochondrial apoptosis pathway – including mitochondrial outer membrane permeabilization, cytochrome c release, and caspase activation – are also present in C. gigas. In contrast, the lack of Bcl-2 homology 3-only subfamily members and apoptotic protease activating factor-1 (APAF-1) protein revealed evolutionary diversity from other phyla. Our results support that mitochondrial apoptosis in animals predates the emergence of vertebrates, but suggest that an unexpectedly diverse mitochondrial apoptosis pathway may exist in invertebrates. In addition, our work provided new clues for an improved understanding of how bivalve acclimate themselves to an inconstant environment.

Differential gene responses in the embryo of the green mussel Perna viridis exposed to dichlorodiphenyltrichloroethane (DDT)

The green-lipped mussel, Perna viridis, is considered to be an ideal indicator for marine environmental pollution. Dichlorodiphenyltrichloroethane (DDT), a typical persistent organic pollutant, is extensively distributed in marine environments. However, little is known about the toxic effects of DDT on the embryo of marine animals, particularly in marine bivalves. Using next-generation sequencing technology, we studied P. viridis embryo after DDT stress at the transcriptome level. A total of 99 202 unigenes were obtained based on the 2383 bp of unigene N50. These differentially expressed genes (DEGs) participated in the various molecular pathways of biological effects, including oxidative stress, detoxification, innate immunity and neurobehavioral disease. Quantitative real-time PCR was performed to verify the mRNA expression of several genes identified by differential gene expression (DGE) analysis. The results indicated that DDT was in induced a dose-dependent manner in the embryo of P. viridis, and most genes involved in oxidative stress and detoxification were up-regulated by DDT exposure; however, the immunity-related genes were down-regulated, except the genes involved in phagocytosis. Gene expression changes in embryo from P. viridis provide a preliminary basis to better understand the molecular toxic response mechanisms of embryo to DDT.

Microbial Diseases of Bivalve Mollusks: Infections, Immunology and Antimicrobial Defense

A variety of bivalve mollusks (phylum Mollusca, class Bivalvia) constitute a prominent commodity in fisheries and aquacultures, but are also crucial in order to preserve our ecosystem’s complexity and function. Bivalve mollusks, such as clams, mussels, oysters and scallops, are relevant bred species, and their global farming maintains a high incremental annual growth rate, representing a considerable proportion of the overall fishery activities. Bivalve mollusks are filter feeders; therefore by filtering a great quantity of water, they may bioaccumulate in their tissues a high number of microorganisms that can be considered infectious for humans and higher vertebrates. Moreover, since some pathogens are also able to infect bivalve mollusks, they are a threat for the entire mollusk farming industry. In consideration of the leading role in aquaculture and the growing financial importance of bivalve farming, much interest has been recently devoted to investigate the pathogenesis of infectious diseases of these mollusks in order to be prepared for public health emergencies and to avoid dreadful income losses. Several bacterial and viral pathogens will be described herein. Despite the minor complexity of the organization of the immune system of bivalves, compared to mammalian immune systems, a precise description of the different mechanisms that induce its activation and functioning is still missing. In the present review, a substantial consideration will be devoted in outlining the immune responses of bivalves and their repertoire of immune cells. Finally, we will focus on the description of antimicrobial peptides that have been identified and characterized in bivalve mollusks. Their structural and antimicrobial features are also of great interest for the biotechnology sector as antimicrobial templates to combat the increasing antibiotic-resistance of different pathogenic bacteria that plague the human population all over the world.

Transcriptome analysis of Ruditapes philippinarum hepatopancreas provides insights into immune signaling pathways under Vibrio anguillarum infection

The Manila clam, Ruditapes philippinarum, is one of the most economically important aquatic clams that are harvested on a large scale by the mariculture industry in China. However, increasing reports of bacterial pathogenic diseases have had a negative effect on the aquaculture industry of R. philippinarum. In the present study, the two transcriptome libraries of untreated (termed H) and challenged Vibrio anguillarum (termed HV) hepatopancreas were constructed and sequenced from Manila clam using an Illumina-based paired-end sequencing platform. In total, 75,302,886 and 66,578,976 high-quality clean reads were assembled from 101,080,746 and 99,673,538 raw data points from the two transcriptome libraries described above, respectively. Furthermore, 156,116 unigenes were generated from 210,685 transcripts, with an N50 length of 1125 bp, and from the annotated SwissProt, NR, NT, KO, GO, KOG and KEGG databases. Moreover, a total of 4071 differentially expressed unigenes (HV vs H) were detected, including 903 up-regulated and 3168 down-regulated genes. Among these differentially expressed unigenes, 226 unigenes were annotated using KEGG annotation in 16 immune-related signaling pathways, including Toll-like receptor, NF-kappa B, MAPK, NOD-like receptor, RIG-I-like receptor, and the TNF and chemokine signaling pathways. Finally, 20,341 simple sequence repeats (SSRs) and 214,430 potential single nucleotide polymorphisms (SNPs) were detected from the H and HV transcriptome libraries. In conclusion, these studies identified many candidate immune-related genes and signaling pathways and conducted a comparative analysis of the differentially expressed unigenes from Manila clam hepatopancreas in response to V. anguillarum stimulation. These data laid the foundation for studying the innate immune systems and defense mechanisms in R. philippinarum.

Evidence of horizontal transmission of the cancer-associated Steamer retrotransposon among ecological cohort bivalve species

Bivalve specimens from legacy frozen tissue collections, and others freshly obtained, were surveyed for the presence of the Steamer long terminal repeat (LTR)-retrotransposon associated with disseminated hemic neoplasia of the soft-shelled clam Mya areneria. Of 22 species investigated using primers for the pol region, only Atlantic M. arenaria, Atlantic and North Sea razor clams Ensis directus, and Baltic clams Macoma balthica from the North Sea were found to possess copies of Steamer in their genomes. Notably, close relatives like Mya truncata and Siliqua patula did not exhibit evidence of Steamer. Amplified Steamer sequences were uniformly identical in all M. areneria specimens, and were highly variable across specimens of E. directus. Variation in the latter included nucleotide polymorphisms among and within individuals as well as length variation in 2 specimens corresponding to the deletion of a predicted stable hairpin structure. Results implicate Atlantic razor clams as the proximal source for horizontal transmission of Steamer among ecologically similar yet markedly distantly related bivalves. The consequences of cross-species transmission of the Steamer retrotransposon are unknown, and the finding of Steamer in 3 bivalve species suggests that further spread is possible.

Construction of a High-Density Genetic Map and Quantitative Trait Locus Mapping in the Manila clam Ruditapes philippinarum

Genetic linkage maps are indispensable tools in a wide range of genetic and genomic research. With the advancement of genotyping-by-sequencing (GBS) methods, the construction of a high-density linkage maps has become achievable in marine organisms lacking sufficient genomic resources, such as mollusks. In this study, high-density linkage map was constructed for an ecologically and commercially important clam species, Ruditapes philippinarum. For the consensus linkage map, a total of 9658 markers spanning 1926.98 cM were mapped to 18 sex-averaged linkage groups, with an average marker distance of 0.42 cM. Based on the high-density linkage map, ten QTLs for growth-related traits and shell color were detected. The coverage and density of the current map are sufficient for us to effectively detect QTL for segregating traits, and two QTL positions were all coincident with the closest markers. This high-density genetic linkage map reveals basic genomic architecture and will be useful for comparative genomics research, genome assembly and genetic improvement of R. philippinarum and other bivalve molluscan species.

Infectious diseases of marine molluscs and host responses as revealed by genomic tools

More and more infectious diseases affect marine molluscs. Some diseases have impacted commercial species including MSX and Dermo of the eastern oyster, QPX of hard clams, withering syndrome of abalone and ostreid herpesvirus 1 (OsHV-1) infections of many molluscs. Although the exact transmission mechanisms are not well understood, human activities and associated environmental changes often correlate with increased disease prevalence. For instance, hatcheries and large-scale aquaculture create high host densities, which, along with increasing ocean temperature, might have contributed to OsHV-1 epizootics in scallops and oysters. A key to understanding linkages between the environment and disease is to understand how the environment affects the host immune system. Although we might be tempted to downplay the role of immunity in invertebrates, recent advances in genomics have provided insights into host and parasite genomes and revealed surprisingly sophisticated innate immune systems in molluscs. All major innate immune pathways are found in molluscs with many immune receptors, regulators and effectors expanded. The expanded gene families provide great diversity and complexity in innate immune response, which may be key to mollusc's defence against diverse pathogens in the absence of adaptive immunity. Further advances in host and parasite genomics should improve our understanding of genetic variation in parasite virulence and host disease resistance.

Comparative study of two single CRD C-type lectins, CgCLec-4 and CgCLec-5, from pacific oyster Crassostrea gigas

C-type lectins (CTLs), a superfamily of Ca²⁺-dependent carbohydrate-recognition proteins, are involved in nonself-recognition and pathogen elimination, and play crucial roles in the innate immunity. In the present study, two single CRD C-type lectins, CgCLec-4 and CgCLec-5, were identified from oyster Crassostrea gigas. The open reading frame (ORF) of CgCLec-4 and CgCLec-5 encoded polypeptides of 152 and 150 amino acids, respectively. Both CgCLec-4 and CgCLec-5 contained one CRD with six conserved cysteines to form three disulfide bridges. The motif in Ca²⁺-binding site 2 of CgCLec-4 was QPE, while it was QYE, a non-a typical motif in CgCLec-5. CgCLec-4 was a secreted lectin with a signal peptide which was highly expressed in hepatopancreas, mantle and hemocytes. CgCLec-5 was an intracellular lectin which was mostly expressed in hemocytes. The lipopolysaccharide stimulation could induce the expressions of CgCLec-4 and CgCLec-5. The recombinant proteins of CgCLec-4 and CgCLec-5 (rCgCLec-4 and rCgCLec-5) could bind to various PAMPs including LPS, PGN, GLU and mannan, while the binding affinity of rCgCLec-5 was stronger than that of rCgCLec-4. Meanwhile, rCgCLec-4 and rCgCLec-5 could bind to different kinds of microorganisms, including Staphylococcus aureus, Escherichia coli and Vibro anguillarum and Yarrowia lipolytica, and the microbial agglutinating ability of rCgCLec-4 was stronger than that of CgCLec-5. Moreover, rCgCLec-4 exhibited anti-microbial activity against bacteria and fungi, but anti-microbial activity of CgCLec-5 was not obvious. All these results suggested that CgCLec-4 and CgCLec-5 could function as an important PRR involved in immune defense against invading pathogen in oyster, and the diversity and complexity of motifs in Ca²⁺ binding site 2 in CRDs determined their comprehensive recognition spectrum and multiple immune functions.

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.

Transcriptomic responses of Perna viridis embryo to Benzo(a)pyrene exposure elucidated by RNA sequencing

The green mussel Perna viridis is an ideal biomonitor to evaluate marine environmental pollution. Benzo(a)pyrene (BaP) is a typical polycyclic aromatic hydrocarbon (PAH), which is well known for the mutagenic and carcinogenic characteristics. However, the toxicological effects of BaP on Perna viridis embryo are still unclear. In this study, we investigated the embryo transcriptomic profile of Perna viridis treated with BaP via digital gene expression analysis. A total of 92,362,742 reads were produced from two groups (control and BaP exposure) by whole transcriptome sequencing (RNA-Seq). Gene Ontology (GO) and Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathway analysis were used on all genes to determine the biological functions and processes. Genes involved in various molecular pathways of toxicological effects were enriched further. The differential expression genes (DEGs) were related to stress response, infectious disease and innate immunity. Quantitative real-time PCR (qRT-PCR) measured expressional levels of six genes confirmed through the DGE analysis. This study reveals that RNA-seq for transcriptome profiling of P. viridis embryo can better understand the embryo toxic effects of BaP. Furthermore, it also suggests that RNA-seq is a superior tool for generating novel and valuable information for revealing the toxic effects caused by BaP at transcriptional level.

Impact of ocean acidification on antimicrobial activity in gills of the blue mussel (Mytilus edulis)

Here, we aimed to investigate potential effects of ocean acidification on antimicrobial peptide (AMP) activity in the gills of Mytilus edulis, as gills are directly facing seawater and the changing pH (predicted to be reduced from ∼8.1 to ∼7.7 by 2100). The AMP activity of gill and haemocyte extracts was compared at pH 6.0, 7.7 and 8.1, with a radial diffusion assay against Escherichia coli. The activity of the gill extracts was not affected by pH, while it was significantly reduced with increasing pH in the haemocyte extracts. Gill extracts were also tested against different species of Vibrio (V. parahaemolyticus, V. tubiashii, V. splendidus, V. alginolyticus) at pH 7.7 and 8.1. The metabolic activity of the bacteria decreased by ∼65-90%, depending on species of bacteria, but was, as in the radial diffusion assay, not affected by pH. The results indicated that AMPs from gills are efficient in a broad pH-range. However, when mussels were pre-exposed for pH 7.7 for four month the gill extracts presented significantly lower inhibit of bacterial growth. A full in-depth proteome investigation of gill extracts, using LC-Orbitrap MS/MS technique, showed that among previously described AMPs from haemocytes of Mytilus, myticin A was found up-regulated in response to lipopolysaccharide, 3 h post injection. Sporadic occurrence of other immune related peptides/proteins also pointed to a rapid response (0.5-3 h p.i.). Altogether, our results indicate that the gills of blue mussels constitute an important first line defence adapted to act at the pH of seawater. The antimicrobial activity of the gills is however modulated when mussels are under the pressure of ocean acidification, which may give future advantages for invading pathogens.

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.

Prediction of Leymus arenarius (L.) antimicrobial peptides based on de novo transcriptome assembly

Leymus arenarius is a unique wild growing Poaceae plant exhibiting extreme tolerance to environmental conditions. In this study we for the first time performed whole-transcriptome sequencing of lymegrass seedlings using Illumina platform followed by de novo transcriptome assembly and functional annotation. Our goal was to identify transcripts encoding antimicrobial peptides (AMPs), one of the key components of plant innate immunity. Using the custom software developed for this study that predicted AMPs and classified them into families, we revealed more than 160 putative AMPs in lymegrass seedlings. We classified them into 7 families based on their cysteine motifs and sequence similarity. The families included defensins, thionins, hevein-like peptides, snakins, cyclotide, alfa-hairpinins and LTPs. This is the first communication about the presence of almost all known AMP families in trascriptomic data of a single plant species. Additionally, cysteine-rich peptides that potentially represent novel families of AMPs were revealed. We have confirmed by RT-PCR validation the presence of 30 transcripts encoding selected AMPs in lymegrass seedlings. In summary, the presented method of pAMP prediction developed by us can be applied for relatively fast and simple screening of novel components of plant immunity system and is well suited for whole-transcriptome or genome analysis of uncharacterized plants.