Pyrosequencing of Mytilus galloprovincialis cDNAs: tissue-specific expression patterns

Cell type and gene regulatory network approaches in the evolution of spiralian biomineralisation

Biomineralisation is the process by which living organisms produce hard structures such as shells and bone. There are multiple independent origins of biomineralised skeletons across the tree of life. This review gives a glimpse into the diversity of spiralian biominerals and what they can teach us about the evolution of novelty. It discusses different levels of biological organisation that may be informative to understand the evolution of biomineralisation and considers the relationship between skeletal and non-skeletal biominerals. More specifically, this review explores if cell type and gene regulatory network approaches could enhance our understanding of the evolutionary origins of biomineralisation.

Long-read RNA sequencing of Pacific abalone Haliotis discus hannai reveals innate immune system responses to environmental stress

Haliotis discus hannai is a commercially important mollusk species, and the abalone aquaculture sector has been jeopardized by deteriorating environmental circumstances such as bacterial infection and thermal stress during the hot summers. However, due to a paucity of genetic information, such as transcriptome resources, our understanding of their stress adaptation is restricted. In this research, using single-molecule long-read (SMRT) sequencing technology, a library composed of ten tissues (i.e., haemocytes, gills, muscle, hepatopancreas, digestive tract, mantle, mucous gland, ovary, testis and head) was constructed and sequenced. In all, 41,855 high-quality unique transcripts, among which 24,778 were successfully annotated. Additionally, 13,463 SSRs, 1,169 transcription factors, and 18,124 lncRNAs were identified in H. discus hannai transcriptome. Furthermore, multiple immune-related transcripts were identified according to KEGG annotation, and a portion of these transcripts were mapped into several classical immune-related pathways, including the PI3K-AKT signaling pathway and Toll-like receptor signaling pathway. Additionally, 24 typical sequences related to the immunity pathway were detected by RT-PCR; the results showed that most of the immune-related genes showed significantly high expression at 72 h after bacterial challenges and thermal stress, especially the expression level of genes in gills was significantly higher than that in haemocytes under V. parahaemolyticus stress at 24 h. At the same time. The analysis of alternative splicing identified several innate immunity-related functions genes, including CD109 and caspase 2. These results suggest that the complex immune system, particularly the powerful innate immunity system, was crucial for H. discus hannai response to numerous environmental challenges.

RNA-seq analysis of the mantle transcriptome from Mytilus edulis during a seasonal spawning event in deep and shallow water culture sites on the northeast coast of Newfoundland, Canada

The blue mussel (Mytilus edulis) has global commercial and ecological importance both in wild and cultured conditions. However there is a qualitative and quantitative lack of knowledge of the molecular mechanisms associated with its reproductive physiology, especially with reference to environmental interactions. Here we initiated a transcriptomic analysis (RNA-sequencing (RNA-seq)) of the mantle from both sexes sampled during a seasonal spawning event and from two culture depths (shallow-5 m; deep- 15 m). Mantle libraries were produced from 3 males and 3 females sampled from each of two shallow sites and two deep sites for a total of 12 replicate male and 12 replicate female libraries (24 total libraries). Overall a total of 2.3 billion raw 100 base reads with an average of 96.5 million reads/library were obtained and assembled into 296,118 transcripts with an average length of 568 bp. Overall, 315 transcripts from male libraries and 25 from female libraries were found to be upregulated in deep water as compared to shallow (edgeR adjusted p value ≤ 0.05). Conversely, 126 transcripts from male libraries and 135 from female libraries were found to be significantly downregulated at the same depth. Thirteen transcripts were selected for qPCR validation based on importance in reproduction, antimicrobial defense and metabolism. Of these, 9 RNA-seq identified transcripts were shown by qPCR to be differentially expressed between groups: 2 were upregulated in deep compared with shallow water (dhx38, mt-co1), 2 were upregulated for female compared with male mantle (pias2, mapkap1) and 6 genes (fndc3a, acbd3, klhl10, ccnb3, armc4, mt-co1) showed to be upregulated in males compared to females. The majority of qPCR studied transcripts were identified as involved in gamete development based on the UniProt database. This study further characterizes the importance of the mantle transcriptome during reproductive activities of M. edulis.

Deciphering mollusc shell production: the roles of genetic mechanisms through to ecology, aquaculture and biomimetics

Most molluscs possess shells, constructed from a vast array of microstructures and architectures. The fully formed shell is composed of calcite or aragonite. These CaCO₃ crystals form complex biocomposites with proteins, which although typically less than 5% of total shell mass, play significant roles in determining shell microstructure. Despite much research effort, large knowledge gaps remain in how molluscs construct and maintain their shells, and how they produce such a great diversity of forms. Here we synthesize results on how shell shape, microstructure, composition and organic content vary among, and within, species in response to numerous biotic and abiotic factors. At the local level, temperature, food supply and predation cues significantly affect shell morphology, whilst salinity has a much stronger influence across latitudes. Moreover, we emphasize how advances in genomic technologies [e.g. restriction site-associated DNA sequencing (RAD-Seq) and epigenetics] allow detailed examinations of whether morphological changes result from phenotypic plasticity or genetic adaptation, or a combination of these. RAD-Seq has already identified single nucleotide polymorphisms associated with temperature and aquaculture practices, whilst epigenetic processes have been shown significantly to modify shell construction to local conditions in, for example, Antarctica and New Zealand. We also synthesize results on the costs of shell construction and explore how these affect energetic trade-offs in animal metabolism. The cellular costs are still debated, with CaCO₃ precipitation estimates ranging from 1-2 J/mg to 17-55 J/mg depending on experimental and environmental conditions. However, organic components are more expensive (~29 J/mg) and recent data indicate transmembrane calcium ion transporters can involve considerable costs. This review emphasizes the role that molecular analyses have played in demonstrating multiple evolutionary origins of biomineralization genes. Although these are characterized by lineage-specific proteins and unique combinations of co-opted genes, a small set of protein domains have been identified as a conserved biomineralization tool box. We further highlight the use of sequence data sets in providing candidate genes for in situ localization and protein function studies. The former has elucidated gene expression modularity in mantle tissue, improving understanding of the diversity of shell morphology synthesis. RNA interference (RNAi) and clustered regularly interspersed short palindromic repeats - CRISPR-associated protein 9 (CRISPR-Cas9) experiments have provided proof of concept for use in the functional investigation of mollusc gene sequences, showing for example that Pif (aragonite-binding) protein plays a significant role in structured nacre crystal growth and that the Lsdia1 gene sets shell chirality in Lymnaea stagnalis. Much research has focused on the impacts of ocean acidification on molluscs. Initial studies were predominantly pessimistic for future molluscan biodiversity. However, more sophisticated experiments incorporating selective breeding and multiple generations are identifying subtle effects and that variability within mollusc genomes has potential for adaption to future conditions. Furthermore, we highlight recent historical studies based on museum collections that demonstrate a greater resilience of molluscs to climate change compared with experimental data. The future of mollusc research lies not solely with ecological investigations into biodiversity, and this review synthesizes knowledge across disciplines to understand biomineralization. It spans research ranging from evolution and development, through predictions of biodiversity prospects and future-proofing of aquaculture to identifying new biomimetic opportunities and societal benefits from recycling shell products.

Influence of light and temperature cycles on the expression of circadian clock genes in the mussel Mytilus edulis

Clock genes and environmental cues regulate essential biological rhythms. The blue mussel, Mytilus edulis, is an ecologically and economically important intertidal bivalve undergoing seasonal reproductive rhythms. We previously identified seasonal expression differences in M. edulis clock genes. Herein, the effects of light/dark cycles, constant darkness, and daily temperature cycles on the circadian expression patterns of such genes are characterised. Clock genes Clk, Cry1, ROR/HR3, Per and Rev-erb/NR1D1, and Timeout-like, show significant mRNA expression variation, persisting in darkness indicating endogenous control. Rhythmic expression was apparent under diurnal temperature cycles in darkness for all except Rev-erb. Temperature cycles induced a significant expression difference in the non-circadian clock-associated gene aaNAT. Furthermore, Suppression Subtractive Hybridisation (SSH) was used to identify seasonal genes with potential links to molecular clock function and revealed numerous genes meriting further investigation. Understanding the relationship between environmental cues and molecular clocks is crucial in predicting the outcomes of environmental change on fundamental rhythmic processes.

Molecular mechanisms of biomineralization in marine invertebrates

Much recent marine research has been directed towards understanding the effects of anthropogenic-induced environmental change on marine biodiversity, particularly for those animals with heavily calcified exoskeletons, such as corals, molluscs and urchins. This is because life in our oceans is becoming more challenging for these animals with changes in temperature, pH and salinity. In the future, it will be more energetically expensive to make marine skeletons and the increasingly corrosive conditions in seawater are expected to result in the dissolution of these external skeletons. However, initial predictions of wide-scale sensitivity are changing as we understand more about the mechanisms underpinning skeletal production (biomineralization). These studies demonstrate the complexity of calcification pathways and the cellular responses of animals to these altered conditions. Factors including parental conditioning, phenotypic plasticity and epigenetics can significantly impact the production of skeletons and thus future population success. This understanding is paralleled by an increase in our knowledge of the genes and proteins involved in biomineralization, particularly in some phyla, such as urchins, molluscs and corals. This Review will provide a broad overview of our current understanding of the factors affecting skeletal production in marine invertebrates. It will focus on the molecular mechanisms underpinning biomineralization and how knowledge of these processes affects experimental design and our ability to predict responses to climate change. Understanding marine biomineralization has many tangible benefits in our changing world, including improvements in conservation and aquaculture and exploitation of natural calcified structure design using biomimicry approaches that are aimed at producing novel biocomposites.

The calcitonin-like system is an ancient regulatory system of biomineralization

Biomineralization is the process by which living organisms acquired the capacity to accumulate minerals in tissues. Shells are the biomineralized exoskeleton of marine molluscs produced by the mantle but factors that regulate mantle shell building are still enigmatic. This study sought to identify candidate regulatory factors of molluscan shell mineralization and targeted family B G-protein coupled receptors (GPCRs) and ligands that include calcium regulatory factors in vertebrates, such as calcitonin (CALC). In molluscs, CALC receptor (CALCR) number was variable and arose through lineage and species-specific duplications. The Mediterranean mussel (Mytilus galloprovincialis) mantle transcriptome expresses six CALCR-like and two CALC-precursors encoding four putative mature peptides. Mussel CALCR-like are activated in vitro by vertebrate CALC but only receptor CALCRIIc is activated by the mussel CALCIIa peptide (EC50 = 2.6 ×10-5 M). Ex-vivo incubations of mantle edge tissue and mantle cells with CALCIIa revealed they accumulated significantly more calcium than untreated tissue and cells. Mussel CALCIIa also significantly decreased mantle acid phosphatase activity, which is associated with shell remodelling. Our data indicate the CALC-like system as candidate regulatory factors of shell mineralization. The identification of the CALC system from molluscs to vertebrates suggests it is an ancient and conserved calcium regulatory system of mineralization.

Bacterial metagenome analysis of Mytilus galloprovincialis collected from Istanbul and Izmir coastal stations of Turkey

Mytilus galloprovincialis is a marine mollusk belonging to the Bivalvia class. It has been distributed largely in Turkish shores and worldwide aquatic environments. Besides being known as an environmental pollution indicator, it is highly consumed as a food and has a high economic value. Due to their nutritional mechanisms by filtering water, they are affected by pollution in seawater and mussels can host-microbial diversity of environmental origin as well as pathogenic bacteria. Therefore, in this study, bacterial species found in Mediterranean mussels collected from the coastal stations of Istanbul [Rumeli Kavagi (RK), Kucukcekmece (KC)], and Izmir [(Foca (MF), Urla (MU)] were investigated and compared with microbiological and metagenomic analyses. According to microbiological analysis results, 34 mussel-associated Enterobacteriaceae and Vibrionaceae family members were identified. As a result of the culture-independent metagenomic analysis, taxonomic groups for each station were identified and compared based on Operational Taxonomic Unit data. For all stations, the most abundant bacterial genera were the unclassified bacterial genera. The total number of mussel-related total richness identified in all groups was 4889 (RK = 1605; KC = 1930; MF = 1508; and MU = 1125). According to the metagenomic data obtained in this study, different relative amounts of Lachnospiraceae and Bacteroidetes taxa groups were reported for all stations. The pathogenic bacterial genera identified by metagenomic analyses which may be significant for the public health are Arcobacter, Clostridium, Aeromonas, Vibrio, Escherichia_Shigella, Klebsiella, Campylobacter, Helicobacter, Pseudomonas, Morganella, Serratia, Corynebacterium, Enterococcus, Staphylococcus, Yersinia, Mycoplasma, Brucellaceae_unclassified, Pantoea, and Proteus.

Immunity-related genes and signaling pathways under hypoxic stresses in Haliotis diversicolor: a transcriptome analysis

Due to increased temperatures and aquaculture density, thermal and hypoxia stresses have become serious problems for the aquaculture of abalone Haliotis diversicolor. Stresses lead to immunosuppression, which can cause severe negative impacts on aquaculture farms. To study the mechanism of immunosuppression after hypoxia stress and bacterial challenge, transcriptomes of H. diversicolor hemocytes involved in immunity were profiled. A total of 307,395,572 clean reads were generated and assembled into 99,774 unigenes. KEGG analysis indicated that 225 unigenes with immunologic function were mapped into immune-related pathways. Expression of 41 unigenes measured by quantitative real-time PCR (qRT-PCR) showed consistent results with that of transcriptome analysis. When exposure challenge of Vibrio parahaemolyticus, it is indicated that the PI3K-AKT, MAPK, NF-κB and P53 signal pathways were involved in the hypoxia-induced immunosuppression of H. diversicolor. Furthermore, when the AKT gene (HdAKT) was inhibited by double-stranded RNA (dsRNA), expression levels of HdAKT was lower than the blank and control group in hemocytes at 4 h, 12 h and 24 h (p < 0.05).

Evolution and diversity of alpha-carbonic anhydrases in the mantle of the Mediterranean mussel (Mytilus galloprovincialis)

The α-carbonic anhydrases (α-CAs) are a large and ancient group of metazoan-specific enzymes. They generate bicarbonate from metabolic carbon dioxide and through calcium carbonate crystal formation play a key role in the regulation of mineralized structures. To better understand how α-CAs contribute to shell mineralization in the marine Mediterranean mussel (Mytilus galloprovincialis) we characterized them in the mantle. Phylogenetic analysis revealed that mollusc α-CA evolution was affected by lineage and species-specific events. Ten α-CAs were found in the Mediterranean mussel mantle and the most abundant form was named, MgNACR, as it grouped with oyster nacreins (NACR). Exposure of the Mediterranean mussel to reduced water salinity (18 vs 37 ppt), caused a significant reduction (p < 0.05) in mantle esterase activity and MgNACR transcript abundance (p < 0.05). Protonograms revealed multiple proteins in the mantle with α-CA hydratase activity and mapped to a protein with a similar size to that deduced for monomeric MgNACR. Our data indicate that MgNACR is a major α-CA enzyme in mantle and that by homology with oyster nacreins likely regulates mussel shell production. We propose that species-dependent α-CA evolution may contribute to explain the diversity of bivalve shell structures and their vulnerability to environmental changes.

Genomics and immunity of the Mediterranean mussel Mytilus galloprovincialis in a changing environment

The Mediterranean mussel (Mytilus galloprovincialis) is a marine invasive species cultured all over the world. Mussels are an appreciated resource in local aquaculture enterprises because of their robust production and resilience that translates into a reliable economic value. So far, no massive mortalities have been reported in natural or cultured populations of this species. In the last years, the knowledge about its immune system has greatly improved but there are still many questions to be answered. One of them is why mussels, with their high filtering activity, are able to be exposed to a high number of potential pathogens without getting infected and without developing an elevated inflammatory response. The sequencing of the mussel genome has revealed a very complex organization with high heterozygosity, abundance of repetitive sequences and extreme intraspecific sequence diversity among individuals, mainly in immune related genes. Among those genes, antimicrobial peptides are the most expressed gene families in mussels, highly polymorphic and with antimicrobial effect against molluscs pathogens, but also against pathogens of lower vertebrates and humans. The combination of a complex genome with the adaptation of mussel immune system to a changing environment could explain this high variability, not only in immune-related genes, but also in the functional response among individuals sampled in the same location and date.

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.

Combined genotyping, microbial diversity and metabolite profiling studies on farmed Mytilus spp. from Kiel Fjord

The blue mussel Mytilus is a popular food source with high economical value. Species of the M. edulis complex (M. edulis, M. galloprovincialis and M. trossulus) hybridise whenever their geographic ranges overlap posing difficulties to species discrimination, which is important for blue mussel aquaculture. The aim of this study was to determine the genetic structure of farmed blue mussels in Kiel Fjord. Microbial and metabolic profile patterns were studied to investigate a possible dependency on the genotype of the bivalves. Genotyping confirmed the complex genetic structure of the Baltic Sea hybrid zone and revealed an unexpected dominance of M. trossulus alleles being in contrast to the predominance of M. edulis alleles described for wild Baltic blue mussels. Culture-dependent and -independent microbial community analyses indicated the presence of a diverse Mytilus-associated microbiota, while an LC-MS/MS-based metabolome study identified 76 major compounds dominated by pigments, alkaloids and polyketides in the whole tissue extracts. Analysis of mussel microbiota and metabolome did not indicate genotypic dependence, but demonstrated high intraspecific variability of farmed mussel individuals. We hypothesise that individual differences in microbial and metabolite patterns may be caused by high individual plasticity and might be enhanced by e.g. nutritional condition, age and gender.

Transcriptome profiling reveals candidate cleft palate-related genes in cultured Chinese sturgeons (Acipenser sinensis)

The Chinese sturgeon (Acipenser sinensis) is an anadromous fish distributed in the Yangtze River and the East China Sea. In this study, we report the novel finding of cleft palates in Chinese sturgeons derived from artificial fertilization. To explore the genetic basis of palate malformation in A. sinensis, Illumina RNA-seq technology was used to analyze the transcriptome data of farmed Chinese sturgeons with normal palates and cleft-palates. Raw reads were obtained and assembled into 808,612 unigenes, with an average length of 509.33 bp and an N50 of 574 bp. Sequence similarity analyses against four public databases (Nr, UniProt, KEGG, and COGs) found 158,642 unigenes that could be annotated. GABAergic synapses and TGF-β signal pathways were the two most enriched pathways with high Rich Factors in the analyses of differentially expressed genes. In these two signal pathways, six genes (GABRA4, GS, GNS, S6K, PITX2, and BMP8) were found as candidate cleft-palate genes in Chinese sturgeon. These findings contribute to our understanding of cleft palate genetics in sturgeon, while simultaneously adding to our knowledge about craniofacial development.

Comparative 16SrDNA Gene-Based Microbiota Profiles of the Pacific Oyster (Crassostrea gigas) and the Mediterranean Mussel (Mytilus galloprovincialis) from a Shellfish Farm (Ligurian Sea, Italy)

The pacific oyster Crassostrea gigas and the Mediterranean mussel Mytilus galloprovincialis are two widely farmed bivalve species which show contrasting behaviour in relation to microbial diseases, with C. gigas being more susceptible and M. galloprovincialis being generally resistant. In a recent study, we showed that different susceptibility to infection exhibited by these two bivalve species may depend on their different capability to kill invading pathogens (e.g., Vibrio spp.) through the action of haemolymph components. Specific microbial-host interactions may also impact bivalve microbiome structure and further influence susceptibility/resistance to microbial diseases. To further investigate this concept, a comparative study of haemolymph and digestive gland 16SrDNA gene-based bacterial microbiota profiles in C. gigas and M. galloprovincialis co-cultivated at the same aquaculture site was carried out using pyrosequencing. Bacterial communities associated with bivalve tissues (hemolymph and digestive gland) were significantly different from those of seawater, and were dominated by relatively few genera such as Vibrio and Pseudoalteromonas. In general, Vibrio accounted for a larger fraction of the microbiota in C. gigas (on average 1.7-fold in the haemolymph) compared to M. galloprovincialis, suggesting that C. gigas may provide better conditions for survival for these bacteria, including potential pathogenic species such as V. aestuarianus. Vibrios appeared to be important members of C. gigas and M. galloprovincialis microbiota and might play a contrasting role in health and disease of bivalve species. Accordingly, microbiome analyses performed on bivalve specimens subjected to commercial depuration highlighted the ineffectiveness of such practice in removing Vibrio species from bivalve tissues.

The impact of depuration on mussel hepatopancreas bacteriome composition and predicted metagenome

Due to the rapid elimination of bacteria through normal behaviour of filter feeding and excretion, the decontamination of hazardous contaminating bacteria from shellfish is performed by depuration. This process, under conditions that maximize shellfish filtering activity, is a useful method to eliminate microorganisms from bivalves. The microbiota composition in bivalves reflects that of the environment of harvesting waters, so quite different bacteriomes would be expected in shellfish collected in different locations. Bacterial accumulation within molluscan shellfish occurs primarily in the hepatopancreas. In order to assess the effect of the depuration process on these different bacteriomes, in this work we used 16S RNA pyrosequencing and metagenome prediction to assess the impact of 15 h of depuration on the whole hepatopancreas bacteriome of mussels collected in three different locations.

Transcriptional response after exposure to domoic acid-producing Pseudo-nitzschia in the digestive gland of the mussel Mytilus galloprovincialis

Bivalve molluscs are filter feeding species that can accumulate biotoxins in their body tissues during harmful algal blooms. Amnesic Shellfish Poisoning (ASP) is caused by species of the diatom genus Pseudo-nitzschia, which produces the toxin domoic acid. The Mytilus galloprovincialis digestive gland transcriptome was de novo assembled based on the sequencing of 12 cDNA libraries, six obtained from control mussels and six from mussels naturally exposed to domoic acid-producing diatom Pseudo-nitzschia australis. After de novo assembly 94,727 transcripts were obtained, with an average length of 1015 bp and a N50 length of 761 bp. The assembled transcripts were clustered (homology > 90%) into 69,294 unigenes. Differential gene expression analysis was performed (DESeq2 algorithm) in the digestive gland following exposure to the toxic algae. A total of 1158 differentially expressed unigenes (absolute fold change > 1.5 and p-value < 0.05) were detected: 686 up-regulated and 472 down-regulated. Several membrane transporters belonging to the family of the SLC (solute carriers) were over-expressed in exposed mussels. Functional enrichment was performed using Pfam annotations obtained from the genes differentially expressed, 37 Pfam families were found to be significantly (FDR adjusted p-value < 0.1) enriched. Some of these families (sulfotransferases, aldo/keto reductases, carboxylesterases, C1q domain and fibrinogen C-terminal globular domain) could be putatively involved in detoxification processes, in the response against of the oxidative stress and in immunological processes. Protein network analysis with STRING algorithm found alteration of the Notch signaling pathway under the action of domoic acid-producing Pseudo-nitzschia. In conclusion, this study provides a high quality reference transcriptome of M. galloprovincialis digestive gland and identifies potential genes involved in the response to domoic acid.

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.

Hepatocellular adenoma in a European flatfish (Limanda limanda): Genetic alterations in laser-capture micro-dissected tissue and global transcriptomic approach

Liver tumours in flatfish have been diagnosed using histopathology for decades to monitor the impacts of marine pollution. Here we describe the application of specific gene (retinoblastoma, Rb) profiling in laser capture micro-dissected samples, and a suppression subtractive hybridization (SSH) approach to isolate differentially expressed genes in hepatocellular adenoma (HCA) samples from dab, Limanda limanda. The Rb profiles from apparently normal and HCA micro-dissected samples of fish from the North Sea showed no significant difference, and genotypic heterogeneity within defined histological phenotypes was observed. In the SSH, sequences associated with cell signalling, cell cycle, gene expression regulation, protein transport and protein degradation were isolated. These included up-regulation of arrestin domain containing 3 (arrdc3), Rac-1 and tribbles, and down-regulation of ankyrin repeat/sterile alpha-motif domain-containing protein 1B-like (ANKS1B-like), c-fos, CDKN1B and RhoA-like sequences, previously implicated in mammalian HCA. This study offers new candidates involved in fish liver tumour development.

Transcriptomic Analysis of Differentially Expressed Genes During Larval Development of Rapana venosa by Digital Gene Expression Profiling

During the life cycle of shellfish, larval development, especially metamorphosis, has a vital influence on the dynamics, distribution, and recruitment of natural populations, as well as seed breeding. Rapana venosa, a carnivorous gastropod, is an important commercial shellfish in China, and is an ecological invader in the United States, Argentina, and France. However, information about the mechanism of its early development is still limited, because research in this area has long suffered from a lack of genomic resources. In this study, 15 digital gene expression (DGE) libraries from five developmental stages of R. venosa were constructed and sequenced on the IIIumina Hi-Sequation 2500 platform. Bioinformaticsanalysis identified numerous differentially and specifically expressed genes, which revealed that genes associated with growth, nervous system, digestive system, immune system, and apoptosis participate in important developmental processes. The functional analysis of differentially expressed genes was further implemented by gene ontology, and Kyoto encyclopedia of genes and genomes enrichment. DGE profiling provided a general picture of the transcriptomic activities during the early development of R. venosa, which may provide interesting hints for further study. Our data represent the first comparative transcriptomic information available for the early development of R. venosa, which is a prerequisite for a better understanding of the physiological traits controlling development.

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.

Exposure to low pH induces molecular level changes in the marine worm, Platynereis dumerilii

Fossil fuel emissions and changes in net land use lead to an increase in atmospheric CO2 concentration and a subsequent decrease of ocean pH. Noticeable effects on organisms' calcification rate, shell structure and energy metabolism have been reported in the literature. To date, little is known about the molecular mechanisms altered under low pH exposure, especially in non-calcifying organisms. We used a suppression subtractive hybridisation (SSH) approach to characterise differentially expressed genes isolated from Platynereis dumerilii, a non-calcifying marine polychaeta species, kept at normal and low pH conditions. Several gene sequences have been identified as differentially regulated. These are involved in processes previously considered as indicators of environment change, such as energy metabolism (NADH dehydrogenase, 2-oxoglutarate dehydrogenase, cytochrome c oxidase and ATP synthase subunit F), while others are involved in cytoskeleton function (paramyosin and calponin) and immune defence (fucolectin-1 and paneth cell-specific alpha-defensin) processes. This is the first study of differential gene expression in a non-calcifying, marine polychaete exposed to low pH seawater conditions and suggests that mechanisms of impact may include additional pathways not previously identified as impacted by low pH in other species.

De novo transcriptome sequencing and analysis of Rapana venosa from six different developmental stages using Hi-seq 2500

The carnivorous whelk Rapana venosa is regarded as a biological invader with strong ecological fitness in the United States, Argentina, France and other countries. R. venosa may seriously damage bivalve resources. Nonetheless, in China, R. venosa is an important commercial species. Larval development, especially metamorphosis, influences the natural population and industrial breeding. However, there are few studies on the early development of R. venosa, and our understanding is further limited by a lack of genomic information. In this study, de novo sequencing was performed to obtain a comprehensive transcriptome profile during early development. A Hi-seq 2500 sequencing run produced 148,737,902 raw reads that were assembled into 1,137,556 unigenes (average length of 619 nucleotides, of which 49,673 could be annotated). The unigenes were assigned to biological processes and functions after annotation in Gene Ontology, eukaryotic Ortholog Groups and Kyoto Encyclopedia of Genes and Genomes. We also identified 93,196 simple sequence repeats among the unigenes. Six unique sequences associated with neuroendocrine function were analyzed by quantitative real-time PCR. Our data represent the first comprehensive transcriptomic resource for R. venosa. Functional annotation of the unigenes involved in various biological processes could stimulate research on the mechanisms of early development in this species. Understanding the mechanism of early development and metamorphosis would benefit antifouling research and aquaculture of R. venosa.

Comparative physiological, biochemical, and molecular thermal stress response profiles for two Unionid freshwater mussel species

Freshwater mussels, aquatic keystone species, are in global decline. Long life spans, sedentary lifestyles, and unique reproductive strategies involving obligate parasitic stages make unionid freshwater mussels particularly sensitive to environmental perturbations resulting from global climate change. A greater understanding of the mechanisms by which closely related species differ in their response to thermal challenge is critical for successful conservation and management practices. As such, both an acute heat shock and a chronic warming simulation were conducted in order to evaluate responses between hypothesized thermally tolerant, Villosa lienosa, and thermally sensitive, Villosa nebulosa, freshwater mussels in response to predicted thermal warming. Multiple biological responses were quantified, including mortality, condition index, growth rates, glycogen and triglyceride content, and candidate gene expression. During acute heat shock, both species upregulated HSP90 and HSP70, though V. lienosa showed consistently greater transcript levels during upregulation. This pattern was consistent during the chronic warming simulation, with V. nebulosa showing greater induction of HSP60. Chronic warming stimulated increases in condition index for V. nebulosa, however declines in growth rates during a recovery period were observed with no concurrent tissue glycogen levels changes. This contrasts with V. lienosa, where tissue glycogen significantly increased during chronic warming, though no response was observed for condition index or growth rates. These differences might indicate disparate thermal stress response mechanisms correlated with metabolic demands and resource utilization. These biological differences could thus be a factor influencing current ranges and these two species‘ future ability to cope with persistent warming in their native habitats.

Allatostatin-type A, kisspeptin and galanin GPCRs and putative ligands as candidate regulatory factors of mantle function

Allatostatin-type A (AST-A), kisspeptin (KISS) and galanin (GAL) G-protein coupled receptor (GPCR) systems share a common ancestral origin in arthropods and the vertebrates where they regulate metabolism and reproduction. The molluscs are the second most diverse phylum in the animal kingdom, they occupy an important phylogenetic position, and their genome is more similar to deuterostomes than the arthropods and nematodes and thus they are good models for studies of gene family evolution and function. This mini-review intends to extend the current knowledge about AST-A, KISS and GAL GPCR system evolution and their putative function in the mollusc mantle. Comparative evolutionary analysis of the target GPCR systems was established by identifying homologues in genomes and tissue transcriptome datasets available for molluscs and comparing them to those of other metazoan systems. Studies in arthropods have revealed the existence of the AST-A system but the loss of homologues of the KISS and GAL systems. Homologues of the insect AST-AR and vertebrate KISSR genes were found in molluscs but putative GALR genes were absent. Receptor gene number suggested that members of this family have suffered lineage specific evolution during the molluscan radiation. In molluscs, orthologues of the insect AST-A peptides were not identified but buccalin peptides that are structurally related were identified and are putative receptor agonists. The identification of AST-AR and KISSR genes in molluscs strengthens the hypotheses that in metazoans members of the AST-AR subfamily share evolutionary proximity with KISSRs. The variable number of receptors and large repertoire of buccalin peptides may be indicative of the functional diversity of the AST-AR/KISSR systems in molluscs. The identification of AST-A and KISS receptors and ligands in the mantle transcriptome indicates that in molluscs they may have acquired a novel function and may play a role in shell development or sensory detection in the mantle.

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.

Key metabolic pathways involved in xenobiotic biotransformation and stress responses revealed by transcriptomics of the mangrove oyster Crassostrea brasiliana

The Brazilian oyster Crassostrea brasiliana was challenged to three common environmental contaminants: phenanthrene, diesel fuel water-accommodated fraction (WAF) and domestic sewage. Total RNA was extracted from the gill and digestive gland, and cDNA libraries were sequenced using the 454 FLX platform. The assembled transcriptome resulted in ̃20,000 contigs, which were annotated to produce the first de novo transcriptome for C. brasiliana. Sequences were screened to identify genes potentially involved in the biotransformation of xenobiotics and associated antioxidant defence mechanisms. These gene families included those of the cytochrome P450 (CYP450), 70kDa heat shock, antioxidants, such as glutathione S-transferase, superoxide dismutase, catalase and also multi-drug resistance proteins. Analysis showed that the massive expansion of the CYP450 and HSP70 family due to gene duplication identified in the Crassostrea gigas genome also occurred in C. brasiliana, suggesting these processes form the base of the Crassostrea lineage. Preliminary expression analyses revealed several candidates biomarker genes that were up-regulated during each of the three treatments, suggesting the potential for environmental monitoring.

Predicting growth and mortality of bivalve larvae using gene expression and supervised machine learning

It is commonly known that the nature of the diet has diverse consequences on larval performance and longevity, however it is still unclear which genes have critical impacts on bivalve development and which pathways are of particular importance in their vulnerability or resistance. First we show that a diet deficient in essential fatty acid (EFA) produces higher larval mortality rates, a reduced shell growth, and lower postlarval performance, all of which are positively correlated with a decline in arachidonic and eicosapentaenoic acids levels, two EFAs known as eicosanoid precursors. Eicosanoids affect the cell inflammatory reactions and are synthesized from long-chain EFAs. Second, we show for the first time that a deficiency in eicosanoid precursors is associated with a network of 29 genes. Their differential regulation can lead to slower growth and higher mortality of Mytilus edulis larvae. Some of these genes are specific to bivalves and others are implicated at the same time in lipid metabolism and defense. Several genes are expressed only during pre-metamorphosis where they are essential for muscle or neurone development and biomineralization, but only in stress-induced larvae. Finally, we discuss how our networks of differentially expressed genes might dynamically alter the development of marine bivalves, especially under dietary influence.

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.

Intersex related gene expression profiles in clams Scrobicularia plana: Molecular markers and environmental application

Intersex, the appearance of female characteristics in male gonads, has been identified in several aquatic species. It is a widespread phenomenon in populations of the bivalve, Scrobicularia plana, from the southwest coast of the U.K. Genes previously identified as differentially expressed (ferritin, testicular haploid expressed gene, THEG, proliferating cell nuclear antigen, PCNA; receptor activated protein kinase C, RACK; cytochrome B, CYB; and cytochrome c oxidase 1, COX1) in intersex clams relative to normal male clams, were selected for characterisation and an environmental survey of the Channel region. Transcripts were significantly differentially expressed at sites with varying intersex incidence and contaminant burdens. Significant correlations between specific gene expressions, key contaminants and sampling locations have been identified, though no single gene was associated with intersex incidence. The results highlight the difficulty in understanding the intersex phenomenon in molluscs where there is still a lack of knowledge on the control of normal reproduction.

Transcriptome-wide analysis of the response of the thecosome pteropod Clio pyramidata to short-term CO2 exposure

Thecosome pteropods, a group of calcifying holoplanktonic mollusks, have recently become a research focus due to their potential sensitivity to increased levels of anthropogenic dissolved CO2 in seawater and the accompanying ocean acidification. Some populations, however, already experience high CO2 in their natural distribution during diel vertical migrations. To achieve a better understanding of the mechanisms of pteropod calcification and physiological response to this sort of short duration CO2 exposure, we characterized the gene complement of Clio pyramidata, a cosmopolitan diel migratory thecosome, and investigated its transcriptomic response to experimentally manipulated CO2 conditions. Individuals were sampled from the Northwest Atlantic in the fall of 2011 and exposed to ambient conditions (~380ppm) and elevated CO2 (~800ppm, similar to levels experienced during a diel vertical migration) for ~10h. Following this exposure the respiration rate of the individuals was measured. We then performed RNA-seq analysis, assembled the C. pyramidata transcriptome de novo, annotated the genes, and assessed the differential gene expression patterns in response to exposure to elevated CO2. Within the transcriptome, we identified homologs of genes with known roles in biomineralization in other mollusks, including perlucin, calmodulin, dermatopontin, calponin, and chitin synthases. Respiration rate was not affected by short-term exposure to CO2. Gene expression varied greatly among individuals, and comparison between treatments indicated that C. pyramidata down-regulated a small number of genes associated with aerobic metabolism and up-regulated genes that may be associated with biomineralization, particularly collagens and C-type lectins. These results provide initial insight into the effects of short term CO2 exposure on these important planktonic open-ocean calcifiers, pairing respiration rate and the gene expression level of response, and reveal candidate genes for future ecophysiological, biomaterial and phylogenetic studies.

Transcriptome profiles of wild and cultured South African abalone, Haliotis midae

This report describes the use of pyrosequencing technologies to generate the first comparative analysis of de novo assembled transcriptome data from cultured and wild specimens of the South African abalone. The transcriptome data and database described here provide a significant genomic resource for abalone research. The data set annotated 11,240 genes, which matched genes with known functions in other species. A large number of transmembrane protein domains (4087) that may indicate a high portion of undiscovered gene receptors were identified. Further, we detected an interesting set of transcription factors (516) that are valuable candidates for participating in regulatory events in developmental (such as cell proliferation and differentiation) and reproductive processes.

Composition and predicted functional ecology of mussel-associated bacteria in Indonesian marine lakes

In the present study, we sampled bacterial communities associated with mussels inhabiting two distinct coastal marine ecosystems in Kalimantan, Indonesia, namely, marine lakes and coastal mangroves. We used 16S rRNA gene pyrosequencing and predicted metagenomic analysis to compare microbial composition and function. Marine lakes are small landlocked bodies of seawater isolated to varying degrees from the open sea environment. They contain numerous endemic taxa and represent natural laboratories of speciation. Our primary goals were to (1) use BLAST search to identify closely related organisms to dominant bacterial OTUs in our mussel dataset and (2) to compare bacterial communities and enrichment in the predicted bacterial metagenome among lakes. Our sequencing effort yielded 3553 OTUs belonging to 44 phyla, 99 classes and 121 orders. Mussels in the largest marine lake (Kakaban) and the coastal mangrove habitat were dominated by bacteria belonging to the phylum Proteobacteria whereas smaller lakes, located on the island of Maratua, were dominated by bacteria belonging to the phyla Firmicutes and Tenericutes. The single most abundant OTU overall was assigned to the genus Mycoplasma. There were several significant differences among locations with respect to metabolic pathways. These included enrichment of xenobiotic biodegradation pathways in the largest marine lake and coastal mangrove. These locations were also the most enriched with respect to nitrogen metabolism. The presence of genes related to isoquinoline alkaloids, polyketides, hydrolases, mono and dioxygenases in the predicted analysis of functional pathways is an indication that the bacterial communities of Brachidontes mussels may be potentially important sources of new marine medicines and enzymes of industrial interest. Future work should focus on measuring how mussel microbial communities influence nutrient dynamics within the marine lake environment and isolating microbes with potential biotechnological applications.

Lectin-like molecules in transcriptome of Littorina littorea hemocytes

The common periwinkle Littorina littorea was introduced in the list of models for comparative immunobiology as a representative of phylogenetically important taxon Caenogastropoda. Using Illumina sequencing technology, we de novo assembled the transcriptome of Littorina littorea hemocytes from 182 million mRNA-Seq pair-end 100 bp reads into a total of 15,526 contigs clustered in 4472 unigenes. The transcriptome profile was analyzed for presence of carbohydrate-binding molecules in a variety of architectural contexts. Hemocytes' repertoire of lectin-like proteins bearing conserved carbohydrate-recognition domains (CRDs) is highly diversified, including 11 of 15 lectin families earlier described in animals, as well as the novel members of lectin family found for the first time in mollusc species. The new molluscan lineage-specific domain combinations were confirmed by cloning and sequencing, including the fuco-lectin related molecules (FLReMs) composed of N-terminal region with no sequence homology to any known protein, a middle Fucolectin Tachylectin-4 Pentaxrin (FTP) domain, and a C-terminal epidermal growth factor (EGF) repeat region. The repertoire of lectin-like molecules is discussed in terms of their potential participation in the receptor phase of immune response. In total, immune-associated functions may be attributed to 70 transcripts belonging to 6 lectin families. These lectin-like genes show low overlap between species of invertebrates, suggesting relatively rapid evolution of immune-associated genes in the group. The repertoire provides valuable candidates for further characterization of the gene functions in mollusc immunity.

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.

De novo transcriptome analysis of Perna viridis highlights tissue-specific patterns for environmental studies

BACKGROUND

The tropical green-lipped mussel Perna viridis is a common biomonitor throughout the Indo-Pacific region that is used for environmental monitoring and ecotoxicological investigations. However, there is limited molecular data available regarding this species. We sought to establish a global transcriptome database from the tissues of adductor muscle, gills and the hepatopancreas of P. viridis in an effort to advance our understanding of the molecular aspects involved during specific toxicity responses in this sentinel species.

RESULTS

Illumina sequencing results yielded 544,272,542 high-quality filtered reads. After de novo assembly using Trinity, 233,257 contigs were generated with an average length of 1,264 bp and an N50 length of 2,868 bp; 192,879 assembled transcripts and 150,111 assembled unigenes were obtained after clustering. A total of 93,668 assembled transcripts (66,692 assembled genes) with putative functions for protein domains were predicted based on InterProScan analysis. Based on similarity searches, 44,713 assembled transcripts and 25,319 assembled unigenes were annotated with at least one BLAST hit. A total of 21,262 assembled transcripts (11,947 assembled genes) were annotated with at least one well-defined Gene Ontology (GO) and 5,131 assembled transcripts (3,181 assembled unigenes) were assigned to 329 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. The quantity of assembled unigenes and transcripts obtained from male and female mussels were similar but varied among the three studied tissues, with the highest numbers recorded in the gills, followed by the hepatopancreas, and then the adductor muscle. Multivariate analyses revealed strong tissue-specific patterns among the three different tissues, but not between sexes in terms of expression profiles for annotated genes in various GO terms, and genes associated with stress responses and degradation of xenobiotics. The expression profiles of certain selected genes in each tissue type were further validated using real-time quantitative polymerase chain reaction assays and a similar tissue-specific trend was seen.

CONCLUSIONS

The extensive sequence data generated from this study will provide a valuable molecular resource for facilitating environmental studies with P. viridis, and highlight the importance of tissue-specific approaches in the future.

454 pyrosequencing-based analysis of gene expression profiles in the amphipod Melita plumulosa: transcriptome assembly and toxicant induced changes

Next generation sequencing using Roche's 454 pyrosequencing platform can be used to generate genomic information for non-model organisms, although there are bioinformatic challenges associated with these studies. These challenges are compounded by a lack of a standardized protocol to either assemble data or to evaluate the quality of a de novo transcriptome. This study presents an assembly of the control and toxicant responsive transcriptome of Melita plumulosa, an Australian amphipod commonly used in ecotoxicological studies. RNA was harvested from control amphipods, juvenile amphipods, and from amphipods exposed to either metal or diesel contaminated sediments. This RNA was used as the basis for a 454 based transcriptome sequencing effort. Sequencing generated 1.3 million reads from control, juvenile, metal-exposed and diesel-exposed amphipods. Different read filtering and assembly protocols were evaluated to generate an assembly that (i) had an optimal number of contigs; (ii) had long contigs; (iii) contained a suitable representation of conserved genes; and (iv) had long ortholog alignment lengths relative to the length of each contig. A final assembly, generated using fixed-length trimming based on the sequence quality scores, followed by assembly using the MIRA algorithm, produced the best results. The 26,625 contigs generated via this approach were annotated using Blast2GO, and the differential expression between treatments and control was determined by mapping with BWA followed by DESeq. Although the mapping generated low coverage, many differentially expressed contigs, including some with known developmental or toxicological function, were identified. This study demonstrated that 454 pyrosequencing is an effective means of generating reference transcriptome information for organisms, such as the amphipod M. plumulosa, that have no genomic information available in databases or in closely related sequenced species. It also demonstrated how optimization of read filtering protocols and assembly approaches changes the utility of results obtained from next generation sequencing studies, and establishes criteria to determine the quality of a de novo assembly in species lacking a reference genome. This new transcriptomic knowledge provides the genomic foundation for the creation of microarray and qPCR assays, serving as a reference transcriptome in future RNAseq studies, and allowing both the biology and ecotoxicology of this organism to be better understood. This approach will allow genomics-based methodology to be applied to a wider range of environmentally relevant species.

Deep sequencing of the mantle transcriptome of the great scallop Pecten maximus

RNA-Seq transcriptome data were generated from mantle tissue of the great scallop, Pecten maximus. The consensus data were produced from a time course series of animals subjected to a 56-day thermal challenge at 3 different temperatures. A total of 26,064 contigs were assembled de novo, providing a useful resource for both the aquaculture community and researchers with an interest in mollusc shell production.

Combined effects of thermal stress and Cd on lysosomal biomarkers and transcription of genes encoding lysosomal enzymes and HSP70 in mussels, Mytilus galloprovincialis

In estuaries and coastal areas, intertidal organisms may be subject to thermal stress resulting from global warming, together with pollution. In the present study, the combined effects of thermal stress and exposure to Cd were investigated in the endo-lysosomal system of digestive cells in mussels, Mytilus galloprovincialis. Mussels were maintained for 24h at 18°C and 26°C seawater temperature in absence and presence of 50 μg Cd/L seawater. Cadmium accumulation in digestive gland tissue, lysosomal structural changes and membrane stability were determined. Semi-quantitative PCR was applied to reveal the changes elicited by the different experimental conditions in hexosaminidase (hex), β-glucuronidase (gusb), cathepsin L (ctsl) and heat shock protein 70 (hsp70) gene transcription levels. Thermal stress provoked lysosomal enlargement whilst Cd-exposure led to fusion of lysosomes. Both thermal stress and Cd-exposure caused lysosomal membrane destabilisation. hex, gusb and ctsl genes but not hsp70 gene were transcriptionally up-regulated as a result of thermal stress. In contrast, all the studied genes were transcriptionally down-regulated in response to Cd-exposure. Cd bioaccumulation was comparable at 18°C and 26°C seawater temperatures but interactions between thermal stress and Cd-exposure were remarkable both in lysosomal biomarkers and in gene transcription. hex, gusb and ctsl genes, reacted to elevated temperature in absence of Cd but not in Cd-exposed mussels. Therefore, thermal stress resulting from global warming might influence the use and interpretation of lysosomal biomarkers in marine pollution monitoring programmes and, vice versa, the presence of pollutants may condition the capacity of mussels to respond against thermal stress in a climate change scenario.

High-throughput transcriptome sequencing and preliminary functional analysis in four Neotropical tree species

Background

The Amazonian rainforest is predicted to suffer from ongoing environmental changes. Despite the need to evaluate the impact of such changes on tree genetic diversity, we almost entirely lack genomic resources.

Results

In this study, we analysed the transcriptome of four tropical tree species (Carapa guianensis, Eperua falcata, Symphonia globulifera and Virola michelii) with contrasting ecological features, belonging to four widespread botanical families (respectively Meliaceae, Fabaceae, Clusiaceae and Myristicaceae). We sequenced cDNA libraries from three organs (leaves, stems, and roots) using 454 pyrosequencing. We have developed an R and bioperl-based bioinformatic procedure for de novo assembly, gene functional annotation and marker discovery. Mismatch identification takes into account single-base quality values as well as the likelihood of false variants as a function of contig depth and number of sequenced chromosomes. Between 17103 (for Symphonia globulifera) and 23390 (for Eperua falcata) contigs were assembled. Organs varied in the numbers of unigenes they apparently express, with higher number in roots. Patterns of gene expression were similar across species, with metabolism of aromatic compounds standing out as an overrepresented gene function. Transcripts corresponding to several gene functions were found to be over- or underrepresented in each organ. We identified between 4434 (for Symphonia globulifera) and 9076 (for Virola surinamensis) well-supported mismatches. The resulting overall mismatch density was comprised between 0.89 (S. globulifera) and 1.05 (V. surinamensis) mismatches/100 bp in variation-containing contigs.

Conclusion

The relative representation of gene functions in the four transcriptomes suggests that secondary metabolism may be particularly important in tropical trees. The differential representation of transcripts among tissues suggests differential gene expression, which opens the way to functional studies in these non-model, ecologically important species. We found substantial amounts of mismatches in the four species. These newly identified putative variants are a first step towards acquiring much needed genomic resources for tropical tree species.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-238) contains supplementary material, which is available to authorized users.

Co-expressed mitochondrial genomes: recently masculinized, recombinant mitochondrial genome is co-expressed with the female-transmitted mtDNA genome in a male Mytilus trossulus mussel from the Baltic Sea

BACKGROUND

Few exceptions have been described from strict maternal inheritance of mitochondrial DNA in animals, including sea mussels (Mytilidae), clams (Donacidae, Veneridae and Solenidae) and freshwater mussels (Unionoidae) order. In these bivalves mitochondria and their DNA are transferred through two separate routes. The females inherit only the maternal mitochondrial DNA whereas the males inherit maternal as well as paternal mitochondrial DNA, which is usually present only in gonads and sperm. The mechanism controlling this phenomenon is unclear but leads to the existence of two separate mitochondrial DNA lineages in a single species. The lineages are usually well differentiated: up to 20-50% divergence in nucleotide sequence. Occasionally, a maternal mitochondrial DNA can invade the paternal transmission route, eventually replacing the diverged M-type and lowering the divergence. Such role reversal (masculinization) event has happened recently in the Mytilus population of the Baltic Sea which consists of M. edulis × M. trossulus hybrids, but the functional status of the resulting mitochondrial genome was unknown.

RESULTS

In this paper we sequenced transcripts from one specimen that was identified as male carrying both the female mitochondrial genome and a recently masculinized mitochondrial genome. Additionally, the analysis of the control region has showed that the recently masculinized, recombinant genome, not only has an M-type control region and all coding regions derived from the F-type, but also is transcriptionally active along side the maternally inherited F-type genome. In the comparative analysis, the two genomes exhibit different substitution patterns, typical for the M vs. F genome comparisons. The genetic distances and ratios of non-synonymous substitutions also suggest that one of the genomes is transitioning from the maternal to the paternal inheritance mode, consistent with its recent masculinization.

CONCLUSION

We have shown, for the first time, that the recently masculinized mitochondrial genome is active and that it accumulates excess of non-synonymous substitutions across its coding sequence. This suggests, that, under certain cytonuclear incompatibility conditions, masculinization may serve to restore the endangered functionality of the paternally inherited genome. This is also another example of a mitochondrial genome in which the recombination in the control region predated its transition from paternal to maternal transmission route.

Biomineral proteins from Mytilus edulis mantle tissue transcriptome

The common blue mussel, Mytilus edulis, has a bimineralic shell composed of approximately equal proportions of the two major polymorphs of calcium carbonate: calcite and aragonite. The exquisite biological control of polymorph production is the focus of research interest in terms of understanding the details of biomineralisation and the proteins involved in the process of complex shell formation. Recent advances in ease and availability of pyrosequencing and assembly have resulted in a sharp increase in transcriptome data for invertebrate biominerals. We have applied Roche 454 pyrosequencing technology to profile the transcriptome for the mantle tissue of the bivalve M. edulis. A comparison was made between the results of several assembly programs: Roche Newbler assembler versions 2.3, 2.5.2 and 2.6 and MIRA 3.2.1 and 3.4.0. The Newbler and MIRA assemblies were subsequently merged using the CAP3 assembler to give a higher consensus in alignments and a more accurate estimate of the true size of the M. edulis transcriptome. Comparison sequence searches show that the mantle transcripts for M. edulis encode putative proteins exhibiting sequence similarities with previously characterised shell proteins of other species of Mytilus, the Bivalvia Pinctada and haliotid gastropods. Importantly, this enhanced transcriptome has detected several transcripts that encode proteins with sequence similarity with previously described shell biomineral proteins including Shematrins and lysine-rich matrix proteins (KRMPs) not previously found in Mytilus.

Reproduction-related genes in the pearl oyster genome

Molluscan reproduction has been a target of biological research because of the various reproductive strategies that have evolved in this phylum. It has also been studied for the development of fisheries technologies, particularly aquaculture. Although fundamental processes of reproduction in other phyla, such as vertebrates and arthropods, have been well studied, information on the molecular mechanisms of molluscan reproduction remains limited. The recently released draft genome of the pearl oyster Pinctada fucata provides a novel and powerful platform for obtaining structural information on the genes and proteins involved in bivalve reproduction. In the present study, we analyzed the pearl oyster draft genome to screen reproduction-related genes. Analysis was mainly conducted for genes reported from other molluscs for encoding orthologs of reproduction-related proteins in other phyla. The gene search in the P. fucata gene models (version 1.1) and genome assembly (version 1.0) were performed using Genome Browser and BLAST software. The obtained gene models were then BLASTP searched against a public database to confirm the best-hit sequences. As a result, more than 40 gene models were identified with high accuracy to encode reproduction-related genes reported for P. fucata and other molluscs. These include vasa, nanos, doublesex- and mab-3-related transcription factor, 5-hydroxytryptamine (5-HT) receptors, vitellogenin, estrogen receptor, and others. The set of reproduction-related genes of P. fucata identified in the present study constitute a new tool for research on bivalve reproduction at the molecular level.

A large-scale gene discovery for the red palm weevil Rhynchophorus ferrugineus (Coleoptera: Curculionidae)

The red palm weevil (RPW; Rhynchophorus ferrugineus) is a devastating pest of palms, prevalent in the Middle East as well as many other regions of the world. Here, we report a large-scale de novo complementary DNA (cDNA) sequencing effort that acquired ∼5 million reads and assembled them into 26 765 contigs from 12 libraries made from samples of different RPW developmental stages based on the Roche/454 GS FLX platform. We annotated these contigs based on the publically available known insect genes and the Tribolium castaneum genome assembly. We find that over 80% of coding sequences (CDS) from the RPW contigs have high-identity homologs to known proteins with complete CDS. Gene expression analysis shows that the pupa and larval stages have the highest and lowest expression levels, respectively. In addition, we also identified more than 60 000 single nucleotide polymorphisms and 1 200 simple sequence repeat markers. This study provides the first large-scale cDNA dataset for RPW, a much-needed resource for future molecular studies.