MytiBase: a knowledgebase of mussel (M. galloprovincialis) transcribed sequences

Hydrocarbon exposure effect on energetic metabolism and immune response in Crassostrea virginica

Crassostrea virginica was exposed to different light crude oil levels to assess the effect on transcriptomic response and metabolic rate. The exposure time was 21 days, and levels of 100 and 200 μg/L were used, including a control. The most significant difference among treatments was the overexpression of several genes associated with energy production, reactive oxygen species (ROS) regulation, immune system response, and inflammatory response. Also, a hydrocarbon concentration-related pattern was identified in ROS regulation, with a gene expression ratio near 1.8:1 between 200 and 100 μg/L treatments. Statistical analysis showed no interaction effect for metabolic rate; however, significant differences were found for oil concentration and time factors, with a higher oxygen consumption at 200 μg/L. Our findings provide novel information about the metabolic response of C. virginica during hydrocarbons exposure. In addition, our results point out which biological processes should be investigated as targets for searching bioindicators.

Transcriptomic analysis of shell repair and biomineralization in the blue mussel, Mytilus edulis

BACKGROUND

Biomineralization by molluscs involves regulated deposition of calcium carbonate crystals within a protein framework to produce complex biocomposite structures. Effective biomineralization is a key trait for aquaculture, and animal resilience under future climate change. While many enzymes and structural proteins have been identified from the shell and in mantle tissue, understanding biomieralization is impeded by a lack of fundamental knowledge of the genes and pathways involved. In adult bivalves, shells are secreted by the mantle tissue during growth, maintenance and repair, with the repair process, in particular, amenable to experimental dissection at the transcriptomic level in individual animals.

RESULTS

Gene expression dynamics were explored in the adult blue mussel, Mytilus edulis, during experimentally induced shell repair, using the two valves of each animal as a matched treatment-control pair. Gene expression was assessed using high-resolution RNA-Seq against a de novo assembled database of functionally annotated transcripts. A large number of differentially expressed transcripts were identified in the repair process. Analysis focused on genes encoding proteins and domains identified in shell biology, using a new database of proteins and domains previously implicated in biomineralization in mussels and other molluscs. The genes implicated in repair included many otherwise novel transcripts that encoded proteins with domains found in other shell matrix proteins, as well as genes previously associated with primary shell formation in larvae. Genes with roles in intracellular signalling and maintenance of membrane resting potential were among the loci implicated in the repair process. While haemocytes have been proposed to be actively involved in repair, no evidence was found for this in the M. edulis data.

CONCLUSIONS

The shell repair experimental model and a newly developed shell protein domain database efficiently identified transcripts involved in M. edulis shell production. In particular, the matched pair analysis allowed factoring out of much of the inherent high level of variability between individual mussels. This snapshot of the damage repair process identified a large number of genes putatively involved in biomineralization from initial signalling, through calcium mobilization to shell construction, providing many novel transcripts for future in-depth functional analyses.

A bivalve biomineralization toolbox

Mollusc shells are a result of the deposition of crystalline and amorphous calcite catalysed by enzymes and shell matrix proteins. Developing a detailed understanding of bivalve mollusc biomineralization pathways is complicated not only by the multiplicity of shell forms and microstructures in this class, but also by the evolution of associated proteins by domain co-option and domain shuffling. In spite of this, a minimal biomineralization toolbox comprising proteins and protein domains critical for shell production across species has been identified. Using a matched pair design to reduce experimental noise from inter-individual variation, combined with damage-repair experiments and a database of biomineralization shell matrix proteins (SMP) derived from published works, proteins were identified that are likely to be involved in shell calcification. Eighteen new, shared proteins likely to be involved in the processes related to the calcification of shells were identified by analysis of genes expressed during repair in Crassostrea gigas, Mytilus edulis and Pecten maximus. Genes involved in ion transport were also identified as potentially involved in calcification either via the maintenance of cell acid-base balance or transport of critical ions to the extrapallial space, the site of shell assembly. These data expand the number of candidate biomineralization proteins in bivalve molluscs for future functional studies and define a minimal functional protein domain set required to produce solid microstructures from soluble calcium carbonate. This is important for understanding molluscan shell evolution, the likely impacts of environmental change on biomineralization processes, materials science, and biomimicry research.

Omics-based molecular analyses of adhesion by aquatic invertebrates

Many aquatic invertebrates are associated with surfaces, using adhesives to attach to the substratum for locomotion, prey capture, reproduction, building or defence. Their intriguing and sophisticated biological glues have been the focus of study for decades. In all but a couple of specific taxa, however, the precise mechanisms by which the bioadhesives stick to surfaces underwater and (in many cases) harden have proved to be elusive. Since the bulk components are known to be based on proteins in most organisms, the opportunities provided by advancing 'omics technologies have revolutionised bioadhesion research. Time-consuming isolation and analysis of single molecules has been either replaced or augmented by the generation of massive data sets that describe the organism's translated genes and proteins. While these new approaches have provided resources and opportunities that have enabled physiological insights and taxonomic comparisons that were not previously possible, they do not provide the complete picture and continued multi-disciplinarity is essential. This review covers the various ways in which 'omics have contributed to our understanding of adhesion by aquatic invertebrates, with new data to illustrate key points. The associated challenges are highlighted and priorities are suggested for future research.

The structure of SeviL, a GM1b/asialo-GM1 binding R-type lectin from the mussel Mytilisepta virgata

SeviL is a recently isolated lectin found to bind to the linear saccharides of the ganglioside GM1b (Neu5Ac\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document}α(2-3)Gal\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta$$\end{document}β(1-3)GalNAc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta$$\end{document}β(1-4)Gal\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta$$\end{document}β(1-4)Glc) and its precursor, asialo-GM1 (Gal\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta$$\end{document}β(1-3)GalNAc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta$$\end{document}β(1-4)Gal\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta$$\end{document}β(1-4)Glc). The crystal structures of recombinant SeviL have been determined in the presence and absence of ligand. The protein belongs to the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta$$\end{document}β-trefoil family, but shows only weak sequence similarity to known structures. SeviL forms a dimer in solution, with one binding site per subunit, close to the subunit interface. Molecular details of glycan recognition by SeviL in solution were analysed by ligand- and protein-based NMR techniques as well as ligand binding assays. SeviL shows no interaction with GM1 due to steric hindrance with the sialic acid branch that is absent from GM1b. This unusual specificity makes SeviL of great interest for the detection and control of certain cancer cells, and cells of the immune system, that display asialo-GM1.

Cryptic diversity in smooth-shelled mussels on Southern Ocean islands: connectivity, hybridisation and a marine invasion

Background

Large numbers of endemic species inhabit subantarctic continental coasts and islands that are characterised by highly variable environmental conditions. Southern hemisphere populations of taxa that are morphologically similar to northern counterparts have traditionally been considered to be extensions of such Northern hemisphere taxa, and may not exhibit differentiation amongst geographically isolated populations in the Southern Ocean. Smooth-shelled blue mussels of the genus Mytilus that exhibit an anti-tropical distribution are a model group to study phylogeography, speciation and hybridisation in the sea, and contribute to the theory and practice of marine biosecurity.

Methods

We used a single nucleotide polymorphism (SNPs) panel that has the ability to accurately identify reference Northern and Southern hemisphere Mytilus taxa to test for evolutionary differentiation amongst native Southern Ocean island populations.

Results

Native mussels from the Falkland Islands and the Kerguelen Islands exhibited greatest affinity to native M. platensis d’Orbigny 1846 from the Atlantic coast of South America. The major Southern Ocean current flow from west to east is likely to explain the spreading of M. platensis to remote offshore islands, as adults via the process of rafting or perhaps directly as larvae. SNPs variation revealed that mussels from Tasmania were native and clearly differentiated from all other blue mussel groups in the Southern and Northern hemispheres. The native mussels M. planulatus from Tasmania and from mainland New Zealand (NZ), and tentatively M. aoteanus from the two NZ Southern Ocean offshore island groups (the Auckland Islands and Campbell Island), formed a distinct M. galloprovincialis–like Southern hemisphere group with closest affinity to Northern hemisphere M. galloprovincialis from the Mediterranean Sea. In all cases, the SNPs revealed evidence of hybridisation between two or more distinct taxa. The invasive Northern hemisphere M. galloprovincialis was identified only in Tasmania, amongst native mussels of a distinct Australian M. planulatus lineage.

Conclusion

Overall, our results reveal that Southern hemisphere island mussels have mixed genome ancestry and are native, not introduced by human activities. The preservation of distinct evolutionary lineages of Southern hemisphere species needs to be an ongoing focus of conservation efforts, given that population sizes on some of the remote offshore oceanic islands will be small and may be more easily adversely affected by invasion and subsequent hybridisation and introgression than larger populations elsewhere.

Electronic supplementary material

The online version of this article (10.1186/s12983-019-0332-y) contains supplementary material, which is available to authorized users.

Season influences the transcriptomic effects of dietary exposure to PVP/PEI coated Ag nanoparticles on mussels Mytilus galloprovincialis

Toxicity of AgNPs has been widely studied in waterborne exposed aquatic organisms. However, toxic effects caused by AgNPs ingested through the diet and depending on the season are still unexplored. The first cell response after exposure to xenobiotics occurs at gene transcription level. Thus, the aim of this study was to assess transcription level effects in the digestive gland of female mussels after dietary exposure to AgNPs both in autumn and in spring. Mussels were fed daily for 21 days with Isochrysis galbana microalgae previously exposed for 24 h to a dose close to environmentally relevant concentrations of 1 μg Ag/L PVP/PEI coated 5 nm AgNPs (in spring) and to a higher dose of 10 μg Ag/L of the same AgNPs both in autumn and in spring. After 1 and 21 days, mussels RNA was hybridized in a custom microarray containing 7806 annotated genes. Mussels were more responsive to the high dose compared to the low dose of AgNPs and a higher number of probes were altered in autumn than in spring. In both seasons, significantly regulated genes were involved in the cytoskeleton and lipid transport and metabolism COG categories, among others, while genes involved in carbohydrate transport and metabolism were specifically altered in autumn. Overall, transcription patterns were differently altered depending on the exposure time and season, indicating that season should be considered in ecotoxicological studies of metal nanoparticles in mussels.

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.

Transcriptomic approach: A promising tool for rapid screening nanomaterial-mediated toxicity in the marine bivalve Mytilus edulis-Application to copper oxide nanoparticles

The extensive development of nanotechnologies will inevitably lead to the release of nanomaterials (NMs) in the environment. As the aquatic environments represent the ultimate sink for various contaminants, it is highly probable that they also constitute a reservoir for NMs and hence aquatic animals represent potential targets. In a regulatory perspective, it is necessary to develop tools to rapidly screen the impact of NMs on model organisms, given that the number of NMs on the market will be increasing. In this context High Throughput Screening approaches represent relevant tools for the investigation of NM-mediated toxicity. The objective of this work was to study the effects of copper oxide nanoparticles (CuONPs) in the marine bivalve Mytilus edulis, using a transcriptomic approach. Mussels were exposed in vivo to CuONPs (10 μg·L^-1CuO NPs) for 24 h and analysis of mRNA expression levels of genes implicated in immune response, antioxidant activities, cell metabolism, cell transport and cytoskeleton was investigated by qPCR on hemocytes and gills. Results showed common effects of CuONPs and its ionic counterpart. However, greater effects of CuONPs on GST, SOD, MT, Actin, ATP synthase gene expressions were observed compared to ionic form indicating that toxicity of CuONPs is not solely due to the release of Cu²⁺. Even though M. edulis genome is not fully characterized, this study provides additional knowledge on the signaling pathways implicated in CuONP-mediated toxicity and demonstrates the reliability of using a qPCR approach to go further in the cellular aspects implicated in response to NPs in marine bivalves.

Omics tools: New challenges in aquatic nanotoxicology?

In recent years, the implication of genomics into ecotoxicological studies has been studied closely to allow a better understanding of organism's responses to environmental contaminants including engineering nanomaterials (ENMs). ENMs are increasingly produced for various applications including cosmetics, electronics, sports equipment, biomedicine and agriculture. Because of their small size, ENMs possess chemical or physical characteristics improved compared to the corresponding macro-sized material. As their application expend, the release of manufactured ENMs into the environment is likely to increase and concern over impacts for the aquatic ecosystem is growing. Several studies reported deleterious effect of ENMs to aquatic organisms, but there is little information about the molecular mechanisms of toxicity. The development of ecotoxicogenomic approaches will improve the characterization of cellular and molecular modes of action of ENMs to aquatic organisms and allow a better prediction of contaminants toxicity. This paper presents an overview of transciptomic/proteomic studies in freshwater and marine organisms exposed to ENMs. Overall, induction of gene expression in relations to defense mechanisms, immune responses, growth and reproduction were measured after ENMs exposures of organisms, but with different patterns depending on exposure duration and concentrations used. In addition, some studies reported a positive correlation between gene expression and cellular modifications, but not at the individual level.

The purplish bifurcate mussel Mytilisepta virgata gene expression atlas reveals a remarkable tissue functional specialization

BACKGROUND

Mytilisepta virgata is a marine mussel commonly found along the coasts of Japan. Although this species has been the subject of occasional studies concerning its ecological role, growth and reproduction, it has been so far almost completely neglected from a genetic and molecular point of view. In the present study we present a high quality de novo assembled transcriptome of the Japanese purplish mussel, which represents the first publicly available collection of expressed sequences for this species.

RESULTS

The assembled transcriptome comprises almost 50,000 contigs, with a N50 statistics of ~1 kilobase and a high estimated completeness based on the rate of BUSCOs identified, standing as one of the most exhaustive sequence resources available for mytiloid bivalves to date. Overall this data, accompanied by gene expression profiles from gills, digestive gland, mantle rim, foot and posterior adductor muscle, presents an accurate snapshot of the great functional specialization of these five tissues in adult mussels.

CONCLUSIONS

We highlight that one of the most striking features of the M. virgata transcriptome is the high abundance and diversification of lectin-like transcripts, which pertain to different gene families and appear to be expressed in particular in the digestive gland and in the gills. Therefore, these two tissues might be selected as preferential targets for the isolation of molecules with interesting carbohydrate-binding properties. In addition, by molecular phylogenomics, we provide solid evidence in support of the classification of M. virgata within the Brachidontinae subfamily. This result is in agreement with the previously proposed hypothesis that the morphological features traditionally used to group Mytilisepta spp. and Septifer spp. within the same clade are inappropriate due to homoplasy.

Computational design of a symmetrical β-trefoil lectin with cancer cell binding activity

Computational protein design has advanced very rapidly over the last decade, but there remain few examples of artificial proteins with direct medical applications. This study describes a new artificial β-trefoil lectin that recognises Burkitt’s lymphoma cells, and which was designed with the intention of finding a basis for novel cancer treatments or diagnostics. The new protein, called “Mitsuba”, is based on the structure of the natural shellfish lectin MytiLec-1, a member of a small lectin family that uses unique sequence motifs to bind α-D-galactose. The three subdomains of MytiLec-1 each carry one galactose binding site, and the 149-residue protein forms a tight dimer in solution. Mitsuba (meaning “three-leaf” in Japanese) was created by symmetry constraining the structure of a MytiLec-1 subunit, resulting in a 150-residue sequence that contains three identical tandem repeats. Mitsuba-1 was expressed and crystallised to confirm the X-ray structure matches the predicted model. Mitsuba-1 recognises cancer cells that express globotriose (Galα(1,4)Galβ(1,4)Glc) on the surface, but the cytotoxicity is abolished.

Crystal structure of MytiLec, a galactose-binding lectin from the mussel Mytilus galloprovincialis with cytotoxicity against certain cancer cell types

MytiLec is a lectin, isolated from bivalves, with cytotoxic activity against cancer cell lines that express globotriaosyl ceramide, Galα(1,4)Galβ(1,4)Glcα1-Cer, on the cell surface. Functional analysis shows that the protein binds to the disaccharide melibiose, Galα(1,6)Glc, and the trisaccharide globotriose, Galα(1,4)Galβ(1,4)Glc. Recombinant MytiLec expressed in bacteria showed the same haemagglutinating and cytotoxic activity against Burkitt's lymphoma (Raji) cells as the native form. The crystal structure has been determined to atomic resolution, in the presence and absence of ligands, showing the protein to be a member of the β-trefoil family, but with a mode of ligand binding unique to a small group of related trefoil lectins. Each of the three pseudo-equivalent binding sites within the monomer shows ligand binding, and the protein forms a tight dimer in solution. An engineered monomer mutant lost all cytotoxic activity against Raji cells, but retained some haemagglutination activity, showing that the quaternary structure of the protein is important for its cellular effects.

Transcriptome Analysis of the Sydney Rock Oyster, Saccostrea glomerata: Insights into Molluscan Immunity

BACKGROUND

Oysters have important ecological functions in their natural environment, acting as global carbon sinks and improving water quality by removing excess nutrients from the water column. During their life-time oysters are exposed to a variety of pathogens that can cause severe mortality in a range of oyster species. Environmental stressors encountered in their habitat can increase the susceptibility of oysters to these pathogens and in general have been shown to impact on oyster immunity, making immune parameters expressed in these marine animals an important research topic.

RESULTS

Paired-end Illumina high throughput sequencing of six S. glomerata tissues exposed to different environmental stressors resulted in a total of 484,121,702 paired-end reads. When reads and assembled transcripts were compared to the C. gigas genome, an overall low level of similarity at the nucleotide level, but a relatively high similarity at the protein level was observed. Examination of the tissue expression pattern showed that some transcripts coding for cathepsins, heat shock proteins and antioxidant proteins were exclusively expressed in the haemolymph of S. glomerata, suggesting a role in innate immunity. Furthermore, analysis of the S. glomerata ORFs showed a wide range of genes potentially involved in innate immunity, from pattern recognition receptors, components of the Toll-like signalling and apoptosis pathways to a complex antioxidant defence mechanism.

CONCLUSIONS

This is the first large scale RNA-Seq study carried out in S. glomerata, showing the complex network of innate immune components that exist in this species. The results confirmed that many of the innate immune system components observed in mammals are also conserved in oysters; however, some, such as the TLR adaptors MAL, TRIF and TRAM are either missing or have been modified significantly. The components identified in this study could help explain the oysters' natural resilience against pathogenic microorganisms encountered in their natural environment.

Mitogenomics of southern hemisphere blue mussels (Bivalvia: Pteriomorphia): Insights into the evolutionary characteristics of the Mytilus edulis complex

Marine blue mussels (Mytilus spp.) are widespread species that exhibit an antitropical distribution with five species occurring in the Northern Hemisphere (M. trossulus, M. edulis, M. galloprovincialis, M. californianus and M. coruscus) and three in the Southern Hemisphere (M. galloprovincialis, M. chilensis and M. platensis). Species limits in this group remain controversial, in particular for those forms that live in South America. Here we investigated structural characteristics of marine mussels mitogenomes, based on published F mtDNA sequences of Northern Hemisphere species and two newly sequenced South American genomes, one from the Atlantic M. platensis and another from the Pacific M. chilensis. These mitogenomes exhibited similar architecture to those of other genomes of Mytilus, including the presence of the Atp8 gene, which is missing in most of the other bivalves. Our evolutionary analysis of mitochondrial genes indicates that purifying selection is the predominant force shaping the evolution of the coding genes. Results of our phylogenetic analyses supported the monophyly of Pteriomorphia and fully resolved the phylogenetic relationships among its five orders. Finally, the low genetic divergence of specimens assigned to M. chilensis and M. platensis suggests that these South American marine mussels represent conspecific variants rather than distinct species.

cDNA and Gene Structure of MytiLec-1, A Bacteriostatic R-Type Lectin from the Mediterranean Mussel (Mytilus galloprovincialis)

MytiLec is an α-d-galactose-binding lectin with a unique primary structure isolated from the Mediterranean mussel (Mytilus galloprovincialis). The lectin adopts a β-trefoil fold that is also found in the B-sub-unit of ricin and other ricin-type (R-type) lectins. We are introducing MytiLec(-1) and its two variants (MytiLec-2 and -3), which both possess an additional pore-forming aerolysin-like domain, as members of a novel multi-genic “mytilectin family” in bivalve mollusks. Based on the full length mRNA sequence (911 bps), it was possible to elucidate the coding sequence of MytiLec-1, which displays an extended open reading frame (ORF) at the 5′ end of the sequence, confirmed both at the mRNA and at the genomic DNA sequence level. While this extension could potentially produce a polypeptide significantly longer than previously reported, this has not been confirmed yet at the protein level. MytiLec-1 was revealed to be encoded by a gene consisting of two exons and a single intron. The first exon comprised the 5′UTR and the initial ATG codon and it was possible to detect a putative promoter region immediately ahead of the transcription start site in the MytiLec-1 genomic locus. The remaining part of the MytiLec-1 coding sequence (including the three sub-domains, the 3′UTR and the poly-A signal) was included in the second exon. The bacteriostatic activity of MytiLec-1 was determined by the agglutination of both Gram-positive and Gram-negative bacteria, which was reversed by the co-presence of α-galactoside. Altogether, these data support the classification of MytiLec-1 as a member of the novel mytilectin family and suggest that this lectin may play an important role as a pattern recognition receptor in the innate immunity of mussels.

Transcriptomics provides insight into Mytilus galloprovincialis (Mollusca: Bivalvia) mantle function and its role in biomineralisation

The mantle is an organ common to all molluscs and is at the forefront of the biomineralisation process. The present study used the Mediterranean mussel (Mytilus galloprovincialis) as a model species to investigate the structural and functional role of the mantle in shell formation. The transcriptomes of three regions of the mantle edge (umbo to posterior edge) were sequenced using Illumina technology which yielded a total of 61,674,325 reads after adapter trimming and filtering. The raw reads assembled into 179,879 transcripts with an N50 value of 1086bp. A total of 1363 transcripts (321, 223 and 816 in regions 1, 2 and 3, respectively) that differed in abundance in the three mantle regions were identified and putative function was assigned to 54% using BLAST sequence similarity searches (cut-off less than 1e(-10)). Morphological differences detected by histology of the three mantle regions was linked to functional heterogeneity by selecting the top five most abundant Pfam domains in the annotated 1363 differentially abundant transcripts across the three mantle regions. Calcium binding domains dominated region two (middle segment of the mantle edge). Candidate biomineralisation genes were mined and tested by qPCR. This revealed that Flp-like, a penicillin binding protein potentially involved in shell matrix maintenance of the Pacific oyster (Crassostrea gigas), had significantly higher expression in the posterior end of the mantle edge (region one). Our findings are intriguing as they indicate that the mantle edge appears to be a heterogeneous tissue, displaying structural and functional bias.

A First Insight into the Genome of the Filter-Feeder Mussel Mytilus galloprovincialis

Mussels belong to the phylum Mollusca, one of the largest and most diverse taxa in the animal kingdom. Despite their importance in aquaculture and in biology in general, genomic resources from mussels are still scarce. To broaden and increase the genomic knowledge in this family, we carried out a whole-genome sequencing study of the cosmopolitan Mediterranean mussel (Mytilus galloprovincialis). We sequenced its genome (32X depth of coverage) on the Illumina platform using three pair-end libraries with different insert sizes. The large number of contigs obtained pointed out a highly complex genome of 1.6 Gb where repeated elements seem to be widespread (~30% of the genome), a feature that is also shared with other marine molluscs. Notwithstanding the limitations of our genome sequencing, we were able to reconstruct two mitochondrial genomes and predict 10,891 putative genes. A comparative analysis with other molluscs revealed a gene enrichment of gene ontology categories related to multixenobiotic resistance, glutamate biosynthetic process, and the maintenance of ciliary structures.

The miRNA biogenesis in marine bivalves

Small non-coding RNAs include powerful regulators of gene expression, transposon mobility and virus activity. Among the various categories, mature microRNAs (miRNAs) guide the translational repression and decay of several targeted mRNAs. The biogenesis of miRNAs depends on few gene products, essentially conserved from basal to higher metazoans, whose protein domains allow specific interactions with dsRNA. Here, we report the identification of key genes responsible of the miRNA biogenesis in 32 bivalves, with particular attention to the aquaculture species Mytilus galloprovincialis and Crassostrea gigas. In detail, we have identified and phylogenetically compared eight evolutionary conserved proteins: DROSHA, DGCR8, EXP5, RAN, DICER TARBP2, AGO and PIWI. In mussels, we recognized several other proteins participating in the miRNA biogenesis or in the subsequent RNA silencing. According to digital expression analysis, these genes display low and not inducible expression levels in adult mussels and oysters whereas they are considerably expressed during development. As miRNAs play an important role also in the antiviral responses, knowledge on their production and regulative effects can shed light on essential molecular processes and provide new hints for disease prevention in bivalves.

Sequencing, De Novo Assembly, and Annotation of the Transcriptome of the Endangered Freshwater Pearl Bivalve, Cristaria plicata, Provides Novel Insights into Functional Genes and Marker Discovery

BACKGROUND

The freshwater mussel Cristaria plicata (Bivalvia: Eulamellibranchia: Unionidae), is an economically important species in molluscan aquaculture due to its use in pearl farming. The species have been listed as endangered in South Korea due to the loss of natural habitats caused by anthropogenic activities. The decreasing population and a lack of genomic information on the species is concerning for environmentalists and conservationists. In this study, we conducted a de novo transcriptome sequencing and annotation analysis of C. plicata using Illumina HiSeq 2500 next-generation sequencing (NGS) technology, the Trinity assembler, and bioinformatics databases to prepare a sustainable resource for the identification of candidate genes involved in immunity, defense, and reproduction.

RESULTS

The C. plicata transcriptome analysis included a total of 286,152,584 raw reads and 281,322,837 clean reads. The de novo assembly identified a total of 453,931 contigs and 374,794 non-redundant unigenes with average lengths of 731.2 and 737.1 bp, respectively. Furthermore, 100% coverage of C. plicata mitochondrial genes within two unigenes supported the quality of the assembler. In total, 84,274 unigenes showed homology to entries in at least one database, and 23,246 unigenes were allocated to one or more Gene Ontology (GO) terms. The most prominent GO biological process, cellular component, and molecular function categories (level 2) were cellular process, membrane, and binding, respectively. A total of 4,776 unigenes were mapped to 123 biological pathways in the KEGG database. Based on the GO terms and KEGG annotation, the unigenes were suggested to be involved in immunity, stress responses, sex-determination, and reproduction. A total of 17,251 cDNA simple sequence repeats (cSSRs) were identified from 61,141 unigenes (size of >1 kb) with the most abundant being dinucleotide repeats.

CONCLUSIONS

This dataset represents the first transcriptome analysis of the endangered mollusc, C. plicata. The transcriptome provides a comprehensive sequence resource for the conservation of genetic information in this species and enrichment of the genetic database. The development of molecular markers will assist in the genetic improvement of C. plicata.

MytiLec, a Mussel R-Type Lectin, Interacts with Surface Glycan Gb3 on Burkitt's Lymphoma Cells to Trigger Apoptosis through Multiple Pathways

MytiLec; a novel lectin isolated from the Mediterranean mussel (Mytilus galloprovincialis); shows strong binding affinity to globotriose (Gb3: Galα1-4Galβ1-4Glc). MytiLec revealed β-trefoil folding as also found in the ricin B-subunit type (R-type) lectin family, although the amino acid sequences were quite different. Classification of R-type lectin family members therefore needs to be based on conformation as well as on primary structure. MytiLec specifically killed Burkitt's lymphoma Ramos cells, which express Gb3. Fluorescein-labeling assay revealed that MytiLec was incorporated inside the cells. MytiLec treatment of Ramos cells resulted in activation of both classical MAPK/ extracellular signal-regulated kinase and extracellular signal-regulated kinase (MEK-ERK) and stress-activated (p38 kinase and JNK) Mitogen-activated protein kinases (MAPK) pathways. In the cells, MytiLec treatment triggered expression of tumor necrosis factor (TNF)-α (a ligand of death receptor-dependent apoptosis) and activation of mitochondria-controlling caspase-9 (initiator caspase) and caspase-3 (activator caspase). Experiments using the specific MEK inhibitor U0126 showed that MytiLec-induced phosphorylation of the MEK-ERK pathway up-regulated expression of the cyclin-dependent kinase inhibitor p21, leading to cell cycle arrest and TNF-α production. Activation of caspase-3 by MytiLec appeared to be regulated by multiple different pathways. Our findings, taken together, indicate that the novel R-type lectin MytiLec initiates programmed cell death of Burkitt's lymphoma cells through multiple pathways (MAPK cascade, death receptor signaling; caspase activation) based on interaction of the lectin with Gb3-containing glycosphingolipid-enriched microdomains on the cell surface.

Structural and Antimicrobial Features of Peptides Related to Myticin C, a Special Defense Molecule from the Mediterranean Mussel Mytilus galloprovincialis

Mussels (Mytilus spp.) have a large repertoire of cysteine-stabilized α,β peptides, and myticin C (MytC) was identified in some hundreds of transcript variants after in vivo immunostimulation. Using a sequence expressed in Italian mussels, we computed the MytC structure and synthesized the mature MytC and related peptide fragments (some of them also prepared in oxidized form) to accurately assess their antibacterial and antifungal activity. Only when tested at pH 5 was the reduced MytC as well as reduced and oxidized fragments including structural β-elements able to inhibit Gram-positive and -negative bacteria (MIC ranges of 4-32 and 8-32 μM, respectively). Such fragments caused selective Escherichia coli killing (MBC of 8-32 μM) but scarcely inhibited two fungal strains. In detail, the antimicrobial β-hairpin MytC[19-40]SOX caused membrane-disrupting effects in E. coli despite its partially ordered conformation in membrane-mimetic environments. In perspective, MytC-derived peptides could be employed to protect acidic mucosal tissues, in cosmetic and food products, and, possibly, as adjuvants in aquaculture.

IL-17 signaling components in bivalves: Comparative sequence analysis and involvement in the immune responses

The recent discovery of soluble immune-regulatory molecules in invertebrates takes advantage of the rapid growth of next generation sequencing datasets. Following protein domain searches in the transcriptomes of 31 bivalve spp. and in few available mollusk genomes, we retrieved 59 domains uniquely identifying interleukin 17 (IL-17) and 96 SEFIR domains typical of IL-17 receptors and CIKS/ACT1 proteins acting downstream in the IL-17 signaling pathway. Compared to the Chordata IL-17 family members, we confirm a separate clustering of the bivalve domain sequences and a consistent conservation pattern of amino acid residues. Analysis performed at transcript and genome level allowed us to propose an updated view of the components outlining the IL-17 signaling pathway in Mytilus galloprovincialis and Crassostrea gigas (in both species, homology modeling reduced the variety of IL-17 domains to only two 3D structures). Digital expression analysis indicated more heterogeneous expression levels for the mussel and oyster IL-17 ligands than for IL-17 receptors and CIKS/CIKSL proteins. Besides, new qPCR analyses confirmed such gene expression trends in hemocytes and gills of mussels challenged with heat-killed bacteria. These results uphold the involvement of an ancient IL-17 signaling pathway in the bivalve immune responses and, likewise in humans, suggest the possibility of distinctive modulatory roles of individual IL-17s/IL-17 receptors. Overall, the common evidence of pro-inflammatory cytokines and inter-related intracellular signaling pathways in bivalves definitely adds complexity to the invertebrate immunity.

RNA-Seq in Mytilus galloprovincialis: comparative transcriptomics and expression profiles among different tissues

Background

The Mediterranean mussel (Mytilus galloprovincialis) is a cosmopolitan, cultured bivalve with worldwide commercial and ecological importance. However, there is a qualitative and quantitative lack of knowledge of the molecular mechanisms involved in the physiology and immune response of this mollusc. In order to start filling this gap, we have studied the transcriptome of mantle, muscle and gills from naïve Mediterranean mussels and hemocytes exposed to distinct stimuli.

Results

A total of 393,316 million raw RNA-Seq reads were obtained and assembled into 151,320 non-redundant transcripts with an average length of 570 bp. Only 55 % of the transcripts were shared across all tissues. Hemocyte and gill transcriptomes shared 60 % of the transcripts while mantle and muscle transcriptomes were most similar, with 77 % shared transcripts. Stimulated hemocytes showed abundant defense and immune-related proteins, in particular, an extremely high amount of antimicrobial peptides. Gills expressed many transcripts assigned to both structure and recognition of non-self patterns, while in mantle many transcripts were related to reproduction and shell formation. Moreover, this tissue presented additional and interesting hematopoietic, antifungal and sensorial functions. Finally, muscle expressed many myofibril and calcium-related proteins and was found to be unexpectedly associated with defense functions. In addition, many metabolic routes related to cancer were represented.

Conclusions

Our analyses indicate that whereas the transcriptomes of these four tissues have characteristic expression profiles in agreement with their biological structures and expected functions, tissue-specific transcriptomes reveal a complex and specialized functions.

Electronic supplementary material

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

Transcriptome analysis of the pearl oyster (Pinctada fucata) hemocytes in response to Vibrio alginolyticus infection

The pearl oyster Pinctada fucata is cultured widely for production of marine pearls in China, while mass mortalities, likely related to pathogenic infections, have occurred frequently in juvenile, mother and operated oysters. To address this issue, understanding host defense mechanisms of P. fucata against pathogenic challenge is extremely important. In the present study, a comparative analysis of hemocyte transcriptomes of P. fucata before and after Vibrio alginolyticus infection was conducted using the Illumina/Hiseq-2000 RNA-Seq technology. A total of 56,345,139 clean reads were generated and then assembled into 74,007 unigenes with an average length of 680 bp and an N50 of 1197 bp. Unigenes were annotated by comparing against non-redundant protein sequence (nr), non-redundant nucleotide (nt), Swiss-Prot, Pfam, Gene Ontology database (GO), Clusters of Orthologous Groups (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases, and 29,615 unigenes (40.01%) were annotated in at least one database. There were 636 genes (518 up-regulated and 118 down-regulated) that were significantly differentially expressed after bacterial challenge, and among which 369 were associated with 122 pathways, including classical immune-related pathways, such as 'MAPK signaling pathway', 'Chemokine signaling pathway', 'Apoptosis' and 'Wnt signaling pathway'. These findings provide information on the pearl oyster innate immunity and may contribute to developing strategies for management of diseases and long-term sustainability of P. fucata culture.

An updated molecular basis for mussel immunity

Non-self recognition with the consequent tolerance or immune reaction is a crucial process to succeed as living organisms. At the same time the interactions between host species and their microbiome, including potential pathogens and parasites, significantly contribute to animal life diversity. Marine filter-feeding bivalves, mussels in particular, can survive also in heavily anthropized coastal waters despite being constantly surrounded by microorganisms. Based on the first outline of the Mytilus galloprovincialis immunome dated 2011, the continuously growing transcript data and the recent release of a draft mussel genome, we explored the available sequence data and scientific literature to reinforce our knowledge on the main gene-encoded elements of the mussel immune responses, from the pathogen recognition to its clearance. We carefully investigated molecules specialized in the sensing and targeting of potential aggressors, expected to show greater molecular diversification, and outlined, whenever relevant, the interconnected cascades of the intracellular signal transduction. Aiming to explore the diversity of extracellular, membrane-bound and intracellular pattern recognition receptors in mussel, we updated a highly complex immune system, comprising molecules which are described here in detail for the first time (e.g. NOD-like receptors) or which had only been partially characterized in bivalves (e.g. RIG-like receptors). Overall, our comparative sequence analysis supported the identification of over 70 novel full-length immunity-related transcripts in M. galloprovincialis. Nevertheless, the multiplicity of gene functions relevant to immunity, the involvement of part of them in other vital processes, and also the lack of a refined mussel genome make this work still not-exhaustive and support the development of more specific studies.

Identification and Characterization of a Novel Family of Cysteine-Rich Peptides (MgCRP-I) from Mytilus galloprovincialis

We report the identification of a novel gene family (named MgCRP-I) encoding short secreted cysteine-rich peptides in the Mediterranean mussel Mytilus galloprovincialis. These peptides display a highly conserved pre-pro region and a hypervariable mature peptide comprising six invariant cysteine residues arranged in three intramolecular disulfide bridges. Although their cysteine pattern is similar to cysteines-rich neurotoxic peptides of distantly related protostomes such as cone snails and arachnids, the different organization of the disulfide bridges observed in synthetic peptides and phylogenetic analyses revealed MgCRP-I as a novel protein family. Genome- and transcriptome-wide searches for orthologous sequences in other bivalve species indicated the unique presence of this gene family in Mytilus spp. Like many antimicrobial peptides and neurotoxins, MgCRP-I peptides are produced as pre-propeptides, usually have a net positive charge and likely derive from similar evolutionary mechanisms, that is, gene duplication and positive selection within the mature peptide region; however, synthetic MgCRP-I peptides did not display significant toxicity in cultured mammalian cells, insecticidal, antimicrobial, or antifungal activities. The functional role of MgCRP-I peptides in mussel physiology still remains puzzling.

Proteomic study of the brackish water mussel Mytilopsis leucophaeata

BACKGROUND

We encountered the opportunity to study proteochemically a brackish water invertebrate animal, Mytilopsis leucophaeata, belonging to the bivalves which stem from the second half of the Cambrian Period (about 510 million years ago). This way, we were able to compare it with the vertebrate animal, the frilled shark (Chlamydoselachus anguineus) that stems from a much later period of geologic time (Permian: 245-286 MYA).

RESULTS

The mussel contains a well-adapted system of protein synthesis on the ER, protein folding on the ER, protein trafficking via COPI or clathrin-coated vesicles from endoplasmic reticulum (ER) to Golgi and plasmalemma, an equally well-developed system of actin filaments that with myosin forms the transport system for vesicular proteins and tubulin, which is also involved in ATP-driven vesicular protein transport via microtubules or transport of chromosomes in mitosis and meiosis. A few of the systems that we could not detect in M. leucophaeata in comparison with C. anguineus are the synaptic vesicle cycle components as synaptobrevin, cellubrevin (v-snare) and synaptosomal associated protein 25-A (t-snare), although one component: Ras-related protein (O-Rab1) could be involved in synaptic vesicle traffic. Another component that we did not find in M. leucophaeata was Rab11 that is involved in the tubulovesicular recycling process of H+/K+-ATPase in C. anguineus. We have not been able to trace the H+/K+-ATPase of M. leucophaeata, but Na+/K+-ATPase was present. Furthermore, we have studied the increase of percent protein expression between 1,070 MYA (the generation of the Amoeba Dictyostelium discoideum) and present (the generation of the mammal Sus scrofa = wild boar). In this time span, three proteomic uprises did occur: 600 to 500 MYA, 47.5 to 4.75 MYA, and 1.4 to 0 MYA. The first uprise covers the generation of bivalves, the second covers gold fish, chicken, brine shrimp, house mouse, rabbit, Japanese medaka and Rattus norvegicus, and the third covers cow, chimpanzee, Homo sapiens, dog, goat, Puccinia graminis and wild boar. We hypothesise that the latter two uprises are related to geological and climate changes and their compensation in protein function expression.

CONCLUSIONS

The proteomic and evolutionary data demonstrate that M. leucophaeata is a highly educatioanal animal to 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.

RNA sequencing and de novo assembly of the digestive gland transcriptome in Mytilus galloprovincialis fed with toxinogenic and non-toxic strains of Alexandrium minutum

Background

The Mediterranean mussel Mytilus galloprovincialis is marine bivalve with a relevant commercial importance as well as a key sentinel organism for the biomonitoring of environmental pollution. Here we report the RNA sequencing of the mussel digestive gland, performed with the aim: a) to produce a high quality de novo transcriptome assembly, thus improving the genetic and molecular knowledge of this organism b) to provide an initial assessment of the response to paralytic shellfish poisoning (PSP) on a molecular level, in order to identify possible molecular markers of toxin accumulation.

Results

The comprehensive de novo assembly and annotation of the transcriptome yielded a collection of 12,079 non-redundant consensus sequences with an average length of 958 bp, with a high percentage of full-length transcripts. The whole-transcriptome gene expression study indicated that the accumulation of paralytic toxins produced by the dinoflagellate Alexandrium minutum over a time span of 5 days scarcely affected gene expression, but the results need further validation with a greater number of biological samples and naturally contaminated specimens.

Conclusion

The digestive gland reference transcriptome we produced significantly improves the data collected from previous sequencing efforts and provides a basic resource for expanding functional genomics investigations in M. galloprovincialis. Although not conclusive, the results of the RNA-seq gene expression analysis support the classification of mussels as bivalves refractory to paralytic shellfish poisoning and point out that the identification molecular biomarkers of PSP in the digestive gland of this organism is problematic.

Electronic supplementary material

The online version of this article (doi:10.1186/1756-0500-7-722) contains supplementary material, which is available to authorized users.

Comparative study of immune responses in the deep-sea hydrothermal vent mussel Bathymodiolus azoricus and the shallow-water mussel Mytilus galloprovincialis challenged with Vibrio bacteria

The deep-sea hydrothermal vent mussel Bathymodiolus azoricus and the continental European coast Mytilus galloprovincialis are two bivalves species living in highly distinct marine habitats. Mussels are filter-feeding animals that may accumulate rapidly bacteria from the environment. Contact with microorganism is thus inevitable during feeding processes where gill tissues assume a strategic importance at the interface between the external milieu and the internal body cavities promoting interactions with potential pathogens during normal filtration and a constant challenge to their immune system. In the present study B. azoricus and M. galloprovincialis were exposed to Vibrio alginolyticus, Vibrio anguillarum and Vibrio splendidus suspensions and to a mixture of these Vibrio suspensions, in order to ascertain the expression level of immune genes in gill samples, from both mussel species. The immune gene expressions were analyzed by means of quantitative-Polymerase Chain Reaction (qPCR). The gene expression results revealed that these bivalve species exhibit significant expression differences between 12 h and 24 h post-challenge times, and between the Vibrio strains used. V. splendidus induced the strongest gene expression level in the two bivalve species whereas the NF-κB and Aggrecan were the most significantly differentially expressed between the two mussel species. When comparing exposure times, both B. azoricus and M. galloprovincialis showed similar percentage of up-regulated genes at 12 h while a marked increased of gene expression was observed at 24 h for the majority of the immune genes in M. galloprovincialis. This contrasts with B. azoricus where the majority of the immune genes were down-regulated at 24 h. The 24 h post-challenge gene expression results clearly bring new evidence supporting time-dependent transcriptional activities resembling acute phase-like responses and different immune responses build-up in these two mussel species when challenged with Vibrio bacteria. High Pressure Liquid Chromatography (HPLC)-Electrospray ionization mass spectrometry (ESI-MS/MS) analyses resulted in different peptide sequences from B. azoricus and M. galloprovincialis gill tissues suggesting that naïve animals present differences, at the protein synthesis level, in their natural environment. B. azoricus proteins sequences, mostly of endosymbiont origin, were related to metabolic, energy production, protein synthesis processes and nutritional demands whereas in M. galloprovincialis putative protein functions were assumed to be related to structural and cellular integrity and signaling functions.

Immune response of apoptosis-related cysteine peptidases from the red abalone Haliotis rufescens (HrCas8 and HrCas3): molecular characterization and transcription expression

Caspases play an important role in the different stages of programmed cell death, or apoptosis, which has been related to the immune response in multicellular organisms. The present study characterized an initiator caspase (HrCas8) and an effector caspase (HrCas3) from the red abalone Haliotis rufescens using the RACE method and qPCR analysis. HrCas8 showed a complete sequence of 2529 base pairs (bp) with an open-reading frame (ORF) of 1911 bp, a 5'UTR of 201 bp, and a 3'UTR of 417 bp. The estimated molecular mass for the 636 amino acids from HrCas8 was 71.5 kDa with an isoelectric point of 6.2. The HrCas8 sequence had two death-effector domains (DEDs) and the subunits p20 and p10, in addition to an active site characteristic of cysteine proteins. Meanwhile, the effector caspase HrCas3 showed a complete sequence of 1404 bp, a 5'UTR of 82 bp, and a 3'UTR of 574 bp. The ORF of this caspase had 747 bp that coded for 248 residues. Moreover, the predicted molecular mass of HrCas3 was 29.4 kDa; the theoretical isoelectric point was 5.70, and the sequence evidenced a conserved caspase recruitment domain (CARD). The distribution of the caspases in distinct tissues revealed that HrCas8 was principally expressed in the hemolymph, while HrCas3 had a higher expression in the gills. A basal level of expression was found for both caspases in muscle tissue. The immune response of caspases in H. rufescens was evaluated through an injection of Vibrio anguillarum. The results showed an increase in the transcription of HrCas8 post-challenge, as well as an activation of HrCas3, which together suggest the initiation of apoptosis as a response to bacterial infection in H. rufescens.

Target capture and massive sequencing of genes transcribed in Mytilus galloprovincialis

Next generation sequencing (NGS) allows fast and massive production of both genome and transcriptome sequence datasets. As the genome of the Mediterranean mussel Mytilus galloprovincialis is not available at present, we have explored the possibility of reducing the whole genome sequencing efforts by using capture probes coupled with PCR amplification and high-throughput 454-sequencing to enrich selected genomic regions. The enrichment of DNA target sequences was validated by real-time PCR, whereas the efficacy of the applied strategy was evaluated by mapping the 454-output reads against reference transcript data already available for M. galloprovincialis and by measuring coverage, SNPs, number of de novo sequenced introns, and complete gene sequences. Focusing on a target size of nearly 1.5 Mbp, we obtained a target coverage which allowed the identification of more than 250 complete introns, 10,741 SNPs, and also complete gene sequences. This study confirms the transcriptome-based enrichment of gDNA regions as a good strategy to expand knowledge on specific subsets of genes also in nonmodel organisms.

Deep transcriptome sequencing of Pecten maximus hemocytes: a genomic resource for bivalve immunology

Pecten maximus, the king scallop, is a bivalve species with important commercial value for both fisheries and aquaculture, traditionally consumed in several European countries. Major problems in larval rearing, however, still limit hatchery-based seed production. High mortalities during early larval stages, likely related to bacterial pathogens, represent the most relevant bottleneck. To address this issue, understanding host defense mechanisms against microbes is extremely important. In this study next-generation RNA-sequencing was carried on scallop hemocytes. To enrich for immune-related transcripts, cDNA libraries from hemocytes challenged in vivo with inactivated-Vibrio anguillarum and in vitro with pathogen-associated molecular patterns, as well as unchallenged controls, were sequenced yielding 216,444,674 sequence reads. De novo assembly of the scallop hemocyte transcriptome consisted of 73,732 contigs (31% annotated). A total of 934 contigs encoded proteins with a known immune function, grouped into several functional categories. Particular attention was reserved to Toll-like receptors (TLRs), a family of pattern recognition receptors (PRRs) involved in non-self recognition. Through mining the scallop hemocyte transcriptome, at least four TLRs could be identified. The organization of canonical TLR domains demonstrated that single cysteine cluster and multiple cysteine cluster TLRs co-exist in this species. In addition, preliminary data concerning their mRNA level following bacterial challenge suggested that different members of this family could exhibit opposite responses to pathogenic stimuli. Finally, a global analysis of differential expression comparing gene-expression levels in in vitro and in vivo stimulated hemocytes against controls provided evidence on a large set of transcripts involved in the great scallop immune response.

Identification of immune-related SNPs in the transcriptome of Mytilus chilensis through high-throughput sequencing

Single nucleotide polymorphisms (SNPs) identified in coding regions represent a useful tool for understanding the immune response against pathogens and stressful environmental conditions. In this study, a SNPs database was generated from transcripts involved in the innate immune response of the mussel Mytilus chilensis. The SNPs were identified through hemocytes transcriptome sequencing from 18 individuals, and SNPs mining was performed in 225,336 contigs, yielding 20,306 polymorphisms associated to immune-related genes. Classification of identified SNPs was based on different pathways of the immune response for Mytilus sp. A total of 28 SNPs were identified in the Toll-like receptor pathway and included 5 non-synonymous polymorphisms; 19 SNPs were identified in the apoptosis pathway and included 3 non-synonymous polymorphisms; 35 SNPs were identified in the Ubiquitin-mediated proteolysis pathway and included 4 non-synonymous variants; and 54 SNPs involved in other molecular functions related to the immune response, such as molecular chaperones, antimicrobial peptides, and genes that interacts with marine toxins were also identified. The molecular markers identified in this work could be useful for novel studies, such as those related to associations between high-resolution molecular markers and functional response to pathogen agents.

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

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

Gene discovery in the tropical scallop Nodipecten nodosus: construction and sequencing of a normalized cDNA library

We report the construction and characterization of a normalized cDNA library from the digestive gland of the marine bivalve Nodipecten nodosus, a commercially valuable tropical scallop. A total of 288 clones were sequenced, and 250 unique sequences were obtained. The cDNA library showed a small sequence redundancy (2.3%) and high numbers of recombinant (99.9%) and independent clones (2.0 × 10(6) cfu), indicating that the cDNA library generated in this study is a profitable resource for efficient gene discovery for N. nodosus. EST functional annotation by Gene Ontology term assignment revealed the identification of sequences potentially involved in aquaculture and ecotoxicology relevant processes such as apoptosis, growth, lipid metabolism, reproduction, development, response to stress and immunity.

cDNA cloning and structural characterization of a lectin from the mussel Crenomytilus grayanus with a unique amino acid sequence and antibacterial activity

An amino acid sequence of GalNAc/Gal-specific lectin from the mussel Crenomytilus grayanus (CGL) was determined by cDNA sequencing. CGL consists of 150 amino acid residues, contains three tandem repeats with high sequence similarities to each other (up to 73%) and does not belong to any known lectins family. According to circular dichroism results CGL is a β/α-protein with the predominance of β-structure. CGL was predicted to adopt a ß-trefoil fold. The lectin exhibits antibacterial activity and might be involved in the recognition and clearance of bacterial pathogens in the shellfish.

Epigenetic modification in neurons of the mollusc Pomacea canaliculata after immune challenge

In human and rodents, the transcriptional response of neurons to stress is related to epigenetic modifications of both DNA and histone proteins. To assess the suitability of simple invertebrate models in studying the basic mechanisms of stress-related epigenetic modifications, we analyzed epigenetic modifications in neurons of the freshwater snail Pomacea canaliculata after the injection of Escherichia coli-derived lipopolysaccharide (LPS). The phospho-acetylation of histone H3, together with the induction of stress-related factors, c-Fos and HSP70, were evaluated in large and small neurons of the pedal ganglia of sham- and LPS-injected snails. Immunocytochemical investigations showed that after LPS injection, the immunopositivity towards phospho (Ser10)-acetyl (Lys14)-histone H3 and c-Fos increases in the nuclei of small gangliar neurons. Western blot analysis confirmed a significant increase of phospho (Ser10)-acetyl (Lys14)-histone H3 in nuclear extracts from 2h LPS-injected animals. c-Fos protein levels were significantly augmented 6h after LPS injection. Immunocytochemistry and western blot indicated that no changes occurred in HSP70 distribution and protein levels. To our knowledge this is the first demonstration of epigenetic changes in molluscan neurons after an immune challenge and indicate the gastropod P. canaliculata as a suitable model for evolutionary and translational studies on stress-related epigenetic modifications.

Construction and characterization of a forward subtracted library of blue mussels Mytilus edulis for the identification of gene transcription signatures and biomarkers of styrene exposure

Transcriptional profiling can elucidate adaptive/toxicity pathways participating in achieving homeostasis or leading to pathogenesis in marine biota exposed to chemical substances. With the aim of analyzing transcriptional responses in the mussel Mytilus edulis exposed to the corrosive and putatively carcinogenic hydrocarbon styrene (3-5 ppm, 3days), a forward subtracted (SSH) cDNA library was produced. Female mussels were selected and digestive gland mRNA was isolated. A library with 1440 clones was produced and a total of 287 clones were sequenced, 53% being identified through BlastN analysis against Mytibase and DeepSeaVent databases. Those genes included GO terms such as 'response to drugs', 'immune defense' and 'cell proliferation'. Furthermore, sequences related to chitin and beta-1-3-glucan metabolism were also up-regulated by styrene. Many of the obtained sequences could not be annotated constituting new mussel sequences. In conclusion, this SSH study reveals novel sequences useful to generate molecular biomarkers of styrene exposure in mussels.

Transcriptional and catalytic responses of antioxidant and biotransformation pathways in mussels, Mytilus galloprovincialis, exposed to chemical mixtures

Antioxidant and biotransformation pathways are widely studied in marine organisms exposed to environmental stressors. However, mechanisms of responses and links between different intracellular levels are not always easy to elucidate and conflicting results are frequently observed between molecular and enzymatic data. In this study, transcriptional and catalytic responses of antioxidant and biotransformation parameters were analyzed after a 4-week exposure of a marine invertebrate, Mytilus galloprovincialis, to chemical mixtures from low polluted and highly polluted sediments. A significant, dose-dependent bioaccumulation was observed for polycyclic aromatic hydrocarbons, especially low molecular weight compounds. Among antioxidant defences, catalase and glutathione peroxidases did not exhibit variations in enzymatic activity, while the corresponding gene transcriptions were up- and down-regulated, respectively; unchanged mRNA levels of superoxide dismutase confirmed the non-synchronous pathways of variations for such antioxidants. Biotransformation responses also revealed inconsistent trends between transcriptional and catalytic variations of glutathione S-transferases, and a significant increase in mRNA levels for cytochrome P450 3A1. The overall results indicated that transcriptional responses might be sensitive but do not necessarily correspond to functional changes, being more useful as "exposure" rather than "effect" biomarkers. Data on gene transcription and catalytic activities should be carefully interpreted when assessing the impact of chemical pollutants and additional studies are needed on modulation of post-transcriptional mechanisms by environmental stressors.

Toll-like receptors and MyD88 adaptors in Mytilus: complete cds and gene expression levels

TLR- and MyD88-related sequences have been previously investigated in Mytibase and then in new transcript reads obtained by Illumina technology from the mussel, Mytilus galloprovincialis. Based on full cds and domain organizations of virtual translations, we identified 23 Toll-like receptors (TLRs) and 3 MyD88 adaptors. MgTLRs can be arranged in 4 clusters according to extra-cellular LRR domain content. MgTLR-b, -i and -k were the only ones containing a multiple cysteine cluster (mccTLR), a domain composition also found in Drosophila Toll-1 and 18-wheeler. The 3 MyD88 we identified in M. galloprovincialis were also retrieved from Mytilus edulis, as well as MgTLR-b and -i. All MgTLRs were constitutively expressed in digestive gland whereas only 4 of them were also present in hemocytes. On the opposite, the 3 MgMyD88s were constitutively expressed in all the tissues. In vivo challenge of M. galloprovincialis with bacteria caused the up regulation of only MgTLR-i, but of all the 3 MgMyD88s. Highest response was induced by Gram-negative Vibrio anguillarum at 9h p.i. Injection of filamentous fungus, Fusarium oxysporum, resulted in up regulation of MgTLR-i and MgMyD88-c at 9h p.i. Such similar pattern of responses suggested MgMyD88-c represents the intra cytoplasm partner of MgTLR-i. Their interaction constituted the first cellular event revealing the existence of a Toll-signaling pathway in Lophotrochozoa.

Gigabase-scale transcriptome analysis on four species of pearl oysters

Pearl oysters have been found to secrete nacre and form pearls with good quality and significant commercial interest. However, the transcriptomic and genomic resources for pearl oysters are still limited. To improve this situation, transcriptome sequencing was conducted from four species of pearl oysters with Illumina HiSeq™ 2000. There were four gigabase-scale transcriptomes for four species of pearl oysters, ∼26.3 million reads with ∼2.37 gigabase base pairs (Gbp) in Pinctada fucata, ∼26.5 million reads with ∼2.39 Gbp in Pinctada margaritifera, ∼27.0 million reads with ∼2.43 Gbp in Pinctada maxima, and ∼25.9 million reads with ∼2.33 Gbp in Pteria penguin, respectively. After sequence assembly and blastx alignment, the numbers of annotated unigenes ≥200 bp were 33,882 in P. fucata, 30,666 in P. margaritifera, 26,420 in P. maxima, and 29,928 in P. penguin. Based on these annotated unigenes among four species of pearl oysters, CDSs were extracted and predicted and furthermore, analyses of GO and KEGG assignments were performed. In addition, 60 putative genes of growth factors and their receptors from four species of pearl oysters were predicted. This study established an excellent resource for gene discovery and expression in pearl oysters, but also offered a significant platform for functional genomics and comparative genomic studies for mollusks.