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

Comparative study of integrated bio-responses in deep-sea and nearshore mussels upon abiotic condition changes: Insight into distinct regulation and adaptation

Deep-sea mussels, one of the dominant species in most deep-sea ecosystems, have long been used as model organisms to investigate the adaptations and symbiotic relationships of deep-sea macrofauna under laboratory conditions due to their ability to survive under atmospheric pressure. However, the impact of additional abiotic conditions beyond pressure, such as temperature and light, on their physiological characteristics remains unknown. In this study, deep-sea mussels (Gigantidas platifrons) from cold seep of the South China Sea, along with nearshore mussels (Mytilus coruscus) from the East China Sea, were reared in unfavorable abiotic conditions for up to 8 days. Integrated biochemical indexes including antioxidant defense, immune ability and energy metabolism were investigated in the gill and digestive gland, while cytotoxicity was determined in hemocytes of both types of mussels. The results revealed mild bio-responses in two types of mussels in the laboratory, represented by the effective antioxidant defense with constant total antioxidant capability level and malondialdehyde content. There were also disparate adaptations in deep-sea and nearshore mussels. In deep-sea mussels, significantly increased immune response and energy reservation were observed in gills, together with the elevated cytotoxicity in hemocytes, implying the more severe biological adaptation was required, mainly due to the symbiotic bacteria loss under laboratory conditions. On the contrary, insignificant biological responses were exhibited in nearshore mussels except for the increased energy consumption, indicating the trade-off strategy to use more energy to deal with potential stress. Overall, this comparative study highlights the basal bio-responses of deep-sea and nearshore mussels out of their native environments, providing evidence that short-term culture of both mussels under easily achievable laboratory conditions would not dramatically alter their biological status. This finding will assist in broadening the application of deep-sea mussels as model organism in future research regardless of the specialized research equipment.

Variation in gene expression along an elevation gradient of Rhododendron sanguineum var. haemaleum assessed in a comparative transcriptomic analysis

Selection along environmental gradients may play a vital role in driving adaptive evolution. Nevertheless, genomic variation and genetic adaptation along environmental clines remains largely unknown in plants in alpine ecosystems. To close this knowledge gap, we assayed transcriptomic profiles of late flower bud and early leaf bud of Rhododendron sanguineum var. haemaleum from four different elevational belts between 3,000 m and 3,800 m in the Gaoligong Mountains. By comparing differences in gene expression of these samples, a gene co-expression network (WGCNA) was constructed to identify candidate genes related to elevation. We found that the overall gene expression patterns are organ-specific for the flower and leaf. Differentially expressed unigenes were identified in these organs. In flowers, these were mainly related to terpenoid metabolism (RsHMGR, RsTPS), while in leaves mainly related to anthocyanin biosynthesis (RsCHS, RsF3'5'H). Terpenoids are the main components of flower scent (fragrance) likely attracting insects for pollination. In response to fewer pollinators at higher elevation zone, it seems relatively less scent is produced in flower organs to reduce energy consumption. Secondary metabolites in leaves such as anthocyanins determine the plants' alternative adaptive strategy to extreme environments, such as selective pressures of insect herbivory from environmental changes and substrate competition in biosynthesis pathways at high elevations. Our findings indicated that the gene expression profiles generated from flower and leaf organs showed parallel expression shifts but with different functionality, suggesting the existence of flexibility in response strategies of plants exposed to heterogeneous environments across elevational gradients. The genes identified here are likely to be involved in the adaptation of the plants to these varying mountainous environments. This study thus contributes to our understanding of the molecular mechanisms of adaptation in response to environmental change.

Transcriptome analysis of halophyte Nitraria tangutorum reveals multiple mechanisms to enhance salt resistance

As a typical halophyte, Nitraria tangutorum Bobr. has attracted the interest of many researchers with the excellent salt tolerance. Elucidation of the mechanism of N. tangutorum salinity tolerance will facilitate the genetic improvement of productive plants faced with salinity. To reveal the molecular response to gradually accumulated salt stress in N. tangutorum, RNA-sequencing and analysis of gradually accumulated NaCl treated samples and control samples were performed, and a total of 1419 differentially expressed genes were identified, including 949 down-regulated genes and 470 up-regulated genes. Detailed analysis uncovered that the catabolism of organic compounds mainly based on oxidative phosphorylation genes was up-regulated. Additionally, various antioxidant genes, especially anthocyanin-related genes, were found to help N. tangutorum remove reactive oxygen species. Moreover, the Mitogen activated protein kinase signaling pathway and other signaling pathways co-regulated various salt tolerance activities. Additionally, intracellular ion homeostasis was maintained via regulation of osmotic regulator-related genes, cutin-related genes, and cell elongation-related genes to retain cellular water and reduce ion concentration. In particularly, simultaneous up-regulation in cytoskeleton-related genes, cell wall-related genes, and auxin-related genes, provided evidence of important role of cell expansion in plant salt tolerance. In conclusion, complex regulatory mechanisms modulated by multiple genes might contribute to the salt tolerance by N. tangutorum.

Differential gene expression and chemical patterns of an intertidal crab inhabiting a polluted port and an adjacent marine protected area

The acquisition of data to safeguard marine protected areas located close to ports is important in order to develop plans that allow effective protection from pollution as well as sustainable development of the port. The area Secche della Meloria is a Marine Protected Area (MPA-MEL) three miles from Livorno Harbour (LH), which is characterized by a long history of pollution. Here we studied the bioaccumulation and transcriptomic patterns of the marbled crab, Pachygrapsus marmoratus (Fabricius, 1787) (Crustacea; Brachyura, Grapsidae), inhabiting the two selected sites. Results showed that the two crab populations are significantly different in their chemical composition of trace elements and Polyciclic Aromatic Hydrocarbons (PAHs), and gene expression patterns (1280 DEGs). Enrichment analysis indicated that crabs at LH had the highest stress response genes, and they were associated with higher levels of bioaccumulation detected in body tissues. We are confident that the significant differential gene expression profiles observed between crabs, characterized by significant chemical differences, is associated with responses to contaminant exposure.

New Insight on Vitality Differences for the Penaeid Shrimp, Fenneropenaeus chinensis, in Low Salinity Environment Through Transcriptomics

Excessive rainfall changes salinity in shrimp farming ponds in short period and exerts low salinity stress on the outdoor breeding shrimp under global warming. Fenneropenaeus chinensis can have different performance on vitality in low salinity environments. To reveal mechanisms of vitality difference in shrimp living in low saline environments. This study based on the normal and moribund F. chinensis in 10 ppt salinity environment using high-throughput sequencing identifies 1,429 differentially expressed genes (DEGs), 586 of which are upregulated, while 843 of which are downregulated in the normal group (FCN10) as compared to the moribund group (FCM10). Meanwhile, another transcriptomic analysis is conducted on the normal and moribund shrimp from 25 ppt (FCN25 vs. FCM25) salinity environment as the control, in which 1,311 DEGs (upregulated: 327 genes, downregulated: 984 genes) are identified. In this study, intersective pathways, GO (Gene Ontology) categories and DEGs from the two groups of comparative transcriptome are investigated. The two intersective pathways (Metabolism of xenobiotics by cytochrome P450, Pentose, and glucuronate interconversions) significantly enriched by DEGs are related to detoxification. In these two pathways, there is one vitality regulation-related gene (VRRG), the Dhdh (dihydrodiol dehydrogenase), which is upregulated in both the groups of FCN10 and FCN25 as compared to the groups of FCM10 and FCM25, respectively. Similarly, in the 25 top intersective GO categories, four VRRGs are revealed. Three of them are upregulated (Itgbl, kielin/chordin-like protein, Slc2a8, solute carrier family 2, facilitated glucose transporter member 8-like protein and Cyp3a30, cytochrome P450 3A30-like protein); one of them is downregulated (Slc6a9, sodium-dependent nutrient amino acid transporter 1-like protein isoform X2). These GO categories are related to transmembrane transporter activity of substance, enzyme inhibitor activity, monooxygenase activity. RT-qPCR analysis further verifies the VRRGs. The study gives new insight into understanding the vitality differences for F. chinensis, in low salinity environment. The pathways and DEGs in response to low salinity stress in modulating the vitality of F. chinensis that could serve as tools in future genetic studies and molecular breeding.

Gene expression and functional analysis of different heat shock protein (HSPs) in Ruditapes philippinarum under BaP stress

Heat shock proteins (HSPs) are a class of highly conserved proteins which can protect cells against various types of stress. However, little information on the mechanism involved in the organic contaminants stress response of HSPs is available, especially in marine invertebrates. The present study was conducted to evaluate the responses of HSPs in clams (Ruditapes philippinarum) under Benzo[a] pyrene (BaP) exposure. The clams were exposed to BaP (concentrations: 0, 0.1, 1, 10 μg/L) for 15 days. 6 HSPs mRNA were classified, and the results of tissue distribution indicated that 4 HSPs gene expressed most in the digestive glands. The transcription level of 6 HSPs (HSP22-1, HSP22-2, HSP40A, HSP60, HSP70, HSP90) genes and the aryl hydrocarbon receptor signaling pathway-related genes, and detoxification system-related enzymes activities were analyzed at 0, 1, 3, 6, 10 and 15 days. The activities of phase II detoxification metabolic enzymes and signaling pathway related genes in clams were severely affected by BaP stress and presented significant difference. Our result suggested that HSPs were produced in the presence of BaP and participated in the process of detoxification metabolism to a certain extent. Additionally, the transcription of HSP40A gene may be used as a potential biomarker of BaP exposure due to its evident concentration- and time-dependent expression pattern. Overall, the study investigated the classification of HSPs in R. philippinarum, provided information about the expression profiles of various HSPs after BaP exposure and broadened the understanding mechanism of HSPs in detoxification defense system under PAHs stress in mollusks.

Transcriptomics reveal triphenyltin-induced molecular toxicity in the marine mussel Perna viridis

Triphenyltin (TPT) is widely used as an active ingredient in antifouling paints and fungicides, and continuous release of this highly toxic endocrine disruptor has caused serious pollution to coastal marine ecosystems and organisms worldwide. Using bioassays and transcriptome sequencing, this study comprehensively investigated the molecular toxicity of TPT chloride (TPTCl) to the marine mussel Perna viridis which is a commercially important species and a common biomonitor for marine pollution in Southeast Asia. Our results indicated that TPTCl was highly toxic to adult P. viridis, with a 96-h LC₁₀ and a 96-h EC₁₀ at 18.7 μg/L and 2.7 μg/L, respectively. A 21-day chronic exposure to 2.7 μg/L TPTCl revealed a strong bioaccumulation of TPT in gills (up to 36.48 μg/g dry weight) and hepatopancreas (71.19 μg/g dry weight) of P. viridis. Transcriptome analysis indicated a time course dependent gene expression pattern in both gills and hepatopancreas. Higher numbers of differentially expressed genes were detected at Day 21 (gills: 1686 genes; hepatopancreas: 1450 genes) and at Day 28 (gills: 628 genes; hepatopancreas: 238 genes) when compared with that at Day 7 (gills: 104 genes, hepatopancreas: 112 genes). Exposure to TPT strongly impaired the endocrine system through targeting on nuclear receptors and putative steroid metabolic genes. Moreover, TPT widely disrupted cellular functions, including lipid metabolism, xenobiotic detoxification, immune response and endoplasmic-reticulum-associated degradation expression, which might have caused the bioaccumulation of TPT in the tissues and aggregation of peptides and proteins in cells that further activated the apoptosis process in P. viridis. Overall, this study has advanced our understanding on both ecotoxicity and molecular toxic mechanisms of TPT to marine mussels, and contributed empirical toxicity data for risk assessment and management of TPT contamination.

Comparative Genomics Reveals 13 Different Isoforms of Mytimycins (A-M) in Mytilus galloprovincialis

Mytimycins are cysteine-rich antimicrobial peptides that show antifungal properties. These peptides are part of the immune network that constitutes the defense system of the Mediterranean mussel (Mytilus galloprovincialis). The immune system of mussels has been increasingly studied in the last decade due to its great efficiency, since these molluscs, particularly resistant to adverse conditions and pathogens, are present all over the world, being considered as an invasive species. The recent sequencing of the mussel genome has greatly simplified the genetic study of some of its immune genes. In the present work, we describe a total of 106 different mytimycin variants in 16 individual mussel genomes. The 13 highly supported mytimycin clusters (A–M) identified with phylogenetic inference were found to be subject to the presence/absence variation, a widespread phenomenon in mussels. We also identified a block of conserved residues evolving under purifying selection, which may indicate the “functional core” of the mature peptide, and a conserved set of 10 invariable plus 6 accessory cysteines which constitute a plastic disulfide array. Finally, we extended the taxonomic range of distribution of mytimycins among Mytilida, identifying novel sequences in M. coruscus, M. californianus, P. viridis, L. fortunei, M. philippinarum, M. modiolus, and P. purpuratus.

De novo transcriptome analysis of the mussel Perna viridis after exposure to the toxic dinoflagellate Prorocentrum lima

Diarrheic shellfish poisoning (DSP) toxins are produced by harmful microalgae and accumulate in bivalve mollusks, causing various toxicity. These toxic effects appear to abate with increasing DSP concentration and longer exposure time, however, the underlying mechanisms remain unclear. To explore the underlying molecular mechanisms, de novo transcriptome analysis of the digestive gland of Perna viridis was performed after Prorocentrum lima exposure. RNA-seq analysis showed that 1886 and 237 genes were up- and down-regulated, respectively after 6 h exposure to P. lima, while 265 genes were up-regulated and 217 genes were down-regulated after 96 h compared to the control. These differentially expressed genes mainly involved in Nrf2 signing pathways, immune stress, apoptosis and cytoskeleton, etc. Combined with qPCR results, we speculated that the mussel P. viridis might mainly rely on glutathione S-transferase (GST) and ABC transporters to counteract DSP toxins during short-term exposure. However, longer exposure of P. lima could activate the Nrf2 signaling pathway and inhibitors of apoptosis protein (IAP), which in turn reduced the damage of DSP toxins to the mussel. DSP toxins could induce cytoskeleton destabilization and had some negative impact on the immune system of bivalves. Collectively, our findings uncovered the crucial molecular mechanisms and the regulatory metabolic nodes that underpin the defense mechanism of bivalves against DSP toxins and also advanced our current understanding of bivalve defense mechanisms.

Transcriptomic response and hydrocarbon accumulation in the eastern oyster (Crassostrea virginica) exposed to crude oil

The adverse effect of crude oil on marine invertebrates is well known. To have a better understanding of its effects on marine invertebrates, Crassostrea virginica was exposed to different concentrations (50, 100 and 200 μg/L) of a mixture of super-light and light crude oil for two weeks, evaluating the transcriptomic response of the digestive gland using RNA-Seq and their accumulation in soft tissues. A total of 33,469,374 reads were assembled, which resulted in 61,356 genome assemblies ('Genes'). Trinotate was used for transcript annotation. At the end of this process, 86,409 transcripts were maintained, comprising a broad set of enzymes from xenobiotics metabolism, oxidative stress, stress and immune responses, and energetic metabolism. The enrichment analysis revealed a change in biological processes and molecular functions, finding from 100 to 200 μg/L. Moreover, the differential gene expression analysis showed a dose-dependent transcriptional response, generally up to 100 μg/L and in some cases up to 200 μg/L, which suggested that oysters' response decreased after 100 μg/L; the analysis of crude oil presence in soft tissues indicated that C. virginica is a suitable candidate for ecotoxicology. Finally, these results should contribute to expanding current genomic resources for C. virginica. Furthermore, they will help to develop new studies in aquatic toxicology focused on knowledge in depth of metabolic pathways, jointly with other approaches (such as proteomics) to allow obtaining a complete idea about the eastern oyster response to crude oil.

High-throughput identification of heavy metal binding proteins from the byssus of chinese green mussel (Perna viridis) by combination of transcriptome and proteome sequencing

The Byssus, which is derived from the foot gland of mussels, has been proved to bind heavy metals effectively, but few studies have focused on the molecular mechanisms behind the accumulation of heavy metals by the byssus. In this study, we integrated high-throughput transcriptome and proteome sequencing to construct a comprehensive protein database for the byssus of Chinese green mussel (Perna viridis), aiming at providing novel insights into the molecular mechanisms by which the byssus binds to heavy metals. Illumina transcriptome sequencing generated a total of 55,670,668 reads. After filtration, we obtained 53,047,718 clean reads and subjected them to de novo assembly using Trinity software. Finally, we annotated 73,264 unigenes and predicted a total of 34,298 protein coding sequences. Moreover, byssal samples were analyzed by proteome sequencing, with the translated protein database from the foot transcriptome as the reference for further prediction of byssal proteins. We eventually determined 187 protein sequences in the byssus, of which 181 proteins are reported for the first time. Interestingly, we observed that many of these byssal proteins are rich in histidine or cysteine residues, which may contribute to the byssal accumulation of heavy metals. Finally, we picked one representative protein, Pvfp-5-1, for recombinant protein synthesis and experimental verification of its efficient binding to cadmium (Cd²⁺) ions.

Systematic identification and validation of the reference genes from 60 RNA-Seq libraries in the scallop Mizuhopecten yessoensis

Background

Reverse transcription quantitative PCR (RT-qPCR) is widely used for gene expression analysis in various organisms. Its accuracy largely relies on the stability of reference genes, making reference gene selection a vital step in RT-qPCR experiments. However, previous studies in mollusks only focused on the reference genes widely used in vertebrates.

Results

In this study, we conducted the transcriptome-wide identification of reference genes in the bivalve mollusk Mizuhopecten yessoensis based on 60 transcriptomes covering early development, adult tissues and gonadal development. A total of 964, 1210 and 2097 candidate reference genes were identified, respectively, resulting in a core set of 568 genes. Functional enrichment analysis showed that these genes are significantly overrepresented in Gene Ontology (GO) terms or Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways related to ribosomes, energy production, etc. Six genes (RS23, EF1A, NDUS4, SELR1, EIF3F, and OLA1) were selected from the candidate genes for RT-qPCR validation, together with 6 commonly used reference genes (ACT, CYTC, HEL, EF1B, GAPDH and RPL16). Stability analyses using geNorm, NormFinder and the comparative delta-Ct method revealed that the new candidate reference genes are more stable than the traditionally used genes, and ACT and CYTC are not recommended under either of the three circumstances. There was a significant correlation between the Ct of RT-qPCR and the log₂(TPM) of RNA-Seq data (Ct = − 0.94 log₂(TPM) + 29.67, R² = 0.73), making it easy to estimate the Ct values from transcriptome data prior to RT-qPCR experiments.

Conclusion

Our study represents the first transcriptome-wide identification of reference genes for early development, adult tissues, and gonadal development in the Yesso scallop and will benefit gene expression studies in other bivalve mollusks.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5661-x) contains supplementary material, which is available to authorized users.

De novo male gonad transcriptome draft for the marine mussel Perumytilus purpuratus with a focus on its reproductive-related proteins

Perumytilus purpuratus is a marine mussel considered a bioengineer species with a broad distribution in the Pacific and Atlantic coast of South America. Studies have shown two geographically and genetically differentiated subpopulations at molecular level and in sperm morphological traits. To open avenues for molecular research on P. purpuratus, a global de novo transcriptome from gonadal tissue of mature males was sequenced using the Illumina platform. From a total of 126.38 million reads, 37,765 transcripts were successfully annotated. BUSCO analysis determined a level of 89% completeness for the assembled transcriptome. The functional gene ontology (GO) annotation indicated that, in terms of abundance, the transcripts related with molecular function were the most represented, followed by those related with biological process and cellular components. Additionally, a subset of GO annotations generated using the "sperm" term resulted in a total of 1,294 sequences where the biological process category was the more represented, with transcripts strongly associated to sperm-processes required for fertilization, and with processes where the sperm-egg interaction could be implicated. Our work will contribute to the evolutionary understanding of the molecular mechanisms related to tissue-specific functions. This work reports the first male gonad transcriptome for the mussel P. purpuratus, generating a useful transcriptomic resource for this species and other closely related mytilids.

Multimarker study of the effects of antifouling biocide on benthic organisms: results using Perna viridis as candidate species

Toxic effects of continuous low dose application of the antifouling biocide chlorine on marine benthic organisms were monitored using transplanted green mussels (Perna viridis) and a suite of biomarkers. Caged mussels were deployed in chlorinated and non-chlorinated sections of the cooling system of an operating electric power plant. Biomarkers indicative of general stress, oxidative stress (superoxide dismutase and catalase), and DNA integrity, along with expression of stress proteins, were studied to assess the effects. Deterioration in condition index with corresponding increase in DNA strand breaks was indicative of chlorine stress. Superoxide dismutase enzyme did not show any particular trend, but catalase activity was high during the initial days of exposure at the chlorinated site; later, it became almost equal to that at the control site. Similarly, expressions of stress proteins (HSP60, HSP70, HSP22, GSTS1, and CYP4) showed bell-shaped pattern during the period of study. Positive correlation among the endpoints indicated the utility of the multimarker approach to monitor the effects of continuous low dose chlorination on mussels.

Analysis of digital gene expression profiling in the gonad of male silkworms (Bombyx mori) under fluoride stress

Fluorine is an essential element, but excessive fluoride can cause serious effects on the respiratory, digestive, and reproductive systems. Fluorine has been suggested to cause reproductive toxicity in vertebrates, but its potential to reproductively affect invertebrates remains unknown. In the present study, the lepidopteran model insect Bombyx mori was used to assess the reproductive toxicity of NaF. The underlying molecular mechanisms were explored by RNA sequencing, and we investigated the testes transcriptomic profile of B. mori treated with NaF via a digital gene expression (DGE) analysis. Among 520 candidate genes, 297 and 223 were identified as significantly upregulated or downregulated, respectively. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were carried out on all genes to determine their biological functions and associated processes. The results indicated that numerous differentially expressed genes are involved in the stress response, detoxification, antibacterial, transport, oxidative phosphorylation, and ribosome. The reliability of the data was confirmed by a quantitative real-time polymerase chain reaction (qRT-PCR) analysis. The changed Glutathione S-transferase (GST) activity and glutathione (GSH) content in the NaF-treated groups were increased and reduced respectively. This study reveals that using RNA-sequencing for the transcriptome profiling of B. mori testes can lead to better comprehension of the male reproductive toxicity effects of NaF. Furthermore, we expect that these results will aid future molecular studies on the reproductive toxicity of NaF in other species.

Gender-specific metabolic responses in gonad of mussel Perna viridis to triazophos

Triazophos, as a lipophilic organophosphate pesticide, displays higher bioaccumulation in the gonads of shellfish. To study the reproductive toxicity of triazophos, we applied metabolomics to characterize the gender-specific metabolic responses in mussel Perna viridis exposed to triazophos. Metabolites were differently altered by triazophos in ovaries of mussel at different concentrations and time intervals, while basically similar metabolic response patterns were observed in male mussels at the two tested concentrations after exposure for 24 and 48h. The significant changes of metabolites in ovaries of mussel exhibited the disturbances in energy metabolism and osmotic regulation, while in male samples triazophos only affected the energy metabolism. Moreover, glycine, sn-glycero-3-phosphocholine, ethanol, aspartate, etc. exhibited consistent variation tendency in both male and female individuals. While the changes of homarine, betaine, taurine, hypotaurine, malonate, β-alanine, succinate, and choline showed obviously gender-specific responses. Overall, this study confirmed the gender-specific responses in gonad of P. viridis to triazophos exposure.

1H NMR-based metabolomics investigation on the effects of petrochemical contamination in posterior adductor muscles of caged mussel Mytilus galloprovincialis

Environmental metabolomics is a high-throughout approach that provides a snapshot of the metabolic status of an organism. In order to elucidate the biological effects of petrochemical contamination on aquatic invertebrates, mussels Mytilus galloprovincialis were caged at the "Augusta-Melilli-Priolo" petrochemical area and Brucoli (Sicily, south Italy), chosen as the reference site. After confirming the elevated concentrations of polycyclic aromatic hydrocarbons (PAHs) and mercury (Hg) in Augusta sediments in our previous work (Maisano et al., 2016a), herein an environmental metabolomics approach based on protonic nuclear magnetic resonance (¹H NMR), coupled with chemometrics, was applied on the mussel posterior adductor muscle (PAM), the main muscular system in bivalve molluscs. Amino acids, osmolytes, energy storage compounds, tricarboxylic acid cycle intermediates, and nucleotides, were found in PAM NMR spectra. Principal Component Analysis (PCA) indicated that mussels caged at the polluted site clustered separately from mussels from the control area, suggesting a clear differentiation between their metabolic profiles. Specifically, disorders in energy metabolism, alterations in amino acids metabolism, and disturbance in the osmoregulatory processes were observed in mussel PAM. Overall, findings from this work demonstrated the usefulness of applying an active biomonitoring strategy for environmental risk assessment, and the effectiveness of metabolomics in elucidating changes in metabolic pathways of aquatic organisms caged at sites differentially contaminated, and thus its suitability to be applied in ecotoxicological studies.

De novo assembly and comparative transcriptome analysis of the foot from Chinese green mussel (Perna viridis) in response to cadmium stimulation

The Chinese green mussel, Perna viridis, is a marine bivalve with important economic values as well as biomonitoring roles for aquatic pollution. Byssus, secreted by the foot gland, has been proved to bind heavy metals effectively. In this study, using the RNA sequencing technology, we performed comparative transcriptomic analysis on the mussel feet with or without inducing by cadmium (Cd). Our current work is aiming at providing insights into the molecular mechanisms of byssus binding to heavy metal ions. The transcriptome sequencing generated a total of 26.13-Gb raw data. After a careful assembly of clean data, we obtained a primary set of 105,127 unigenes, in which 32,268 unigenes were annotated. Based on the expression profiles, we identified 9,048 differentially expressed genes (DEGs) between Cd treatment (50 or 100 μg/L) at 48 h and the control, suggesting an extensive transcriptome response of the mussels during the Cd stimulation. Moreover, we observed that the expression levels of 54 byssus protein coding genes increased significantly after the 48-h Cd stimulation. In addition, 16 critical byssus protein coding genes were picked for profiling by quantitative real-time PCR (qRT-PCR). Finally, we reached a primary conclusion that high content of tyrosine (Tyr), cysteine (Cys), histidine (His) residues or the special motif plays an important role in the accumulation of heavy metals in byssus. We also proposed an interesting model for the confirmed byssal Cd accumulation, in which biosynthesis of byssus proteins may play simultaneously critical roles since their transcription levels were significantly elevated.

The response of Isidorella newcombi to copper exposure: Using an integrated biological framework to interpret transcriptomic responses from RNA-seq analysis

This study describes the transcriptomic response of the Australian endemic freshwater gastropod Isidorella newcombi exposed to 80±1μg/L of copper for 3days. Analysis of copper tissue concentration, lysosomal membrane destabilisation and RNA-seq were conducted. Copper tissue concentrations confirmed that copper was bioaccumulated by the snails. Increased lysosomal membrane destabilisation in the copper-exposed snails indicated that the snails were stressed as a result of the exposure. Both copper tissue concentrations and lysosomal destabilisation were significantly greater in snails exposed to copper. In order to interpret the RNA-seq data from an ecotoxicological perspective an integrated biological response model was developed that grouped transcriptomic responses into those associated with copper transport and storage, survival mechanisms and cell death. A conceptual model of expected transcriptomic changes resulting from the copper exposure was developed as a basis to assess transcriptomic responses. Transcriptomic changes were evident at all the three levels of the integrated biological response model. Despite lacking statistical significance, increased expression of the gene encoding copper transporting ATPase provided an indication of increased internal transport of copper. Increased expression of genes associated with endocytosis are associated with increased transport of copper to the lysosome for storage in a detoxified form. Survival mechanisms included metabolic depression and processes associated with cellular repair and recycling. There was transcriptomic evidence of increased cell death by apoptosis in the copper-exposed organisms. Increased apoptosis is supported by the increase in lysosomal membrane destabilisation in the copper-exposed snails. Transcriptomic changes relating to apoptosis, phagocytosis, protein degradation and the lysosome were evident and these processes can be linked to the degradation of post-apoptotic debris. The study identified contaminant specific transcriptomic markers as well as markers of general stress. From an ecotoxicological perspective, the use of a framework to group transcriptomic responses into those associated with copper transport, survival and cell death assisted with the complex process of interpretation of RNA-seq data. The broad adoption of such a framework in ecotoxicology studies would assist in comparison between studies and the identification of reliable transcriptomic markers of contaminant exposure and response.

Transcriptome profiles of the clam Meretrix petechialis hepatopancreas in response to Vibrio infection

Microbial diseases have received much attention due to their enormous destruction of aquaculture, and Vibrio parahaemolyticus is one of the main pathogens that cause bacterial disease in the clam Meretrix petechialis. To better understand the molecular mechanisms of the immune response to Vibrio in M. petechialis, RNA-Seq was applied to explore global expression changes of hepatopancreas from this clam after Vibrio challenge. There were 199,318,966 clean reads obtained by Illumina sequencing, which were further assembled into 214,577 transcripts, and then 147,255 unigenes with an N50 of 1393 bp were identified. Gene ontology (GO) analysis revealed 21 biological process subcategories, 15 cellular component subcategories and 12 molecular function subcategories. A total of 8358 unigenes were mapped onto 267 biological signaling pathways by KEGG, among which there were 16 pathways related to the immune system. In total, 206 differentially expressed genes (DEGs) were identified, including 113 up-regulated unigenes and 93 down-regulated unigenes. In these DEGs, 96 DEGs were annotated in at least one database, accounting for 46.60% of all significant DEGs. To validate the transcriptome dataset, 15 DEGs were selected for real-time qPCR confirmation and the results showed that expression patterns of 13 genes (86.7%) agreed well with the RNA-Seq analysis. Fourteen of the 206 DEGs were annotated to be immune-related genes, and we examined the expression patterns of four immune-related DEGs using clams post immersion challenge. This study enriched the M. petechialis transcriptome database and provided insight into the immune response of M. petechialis against Vibrio infection.

Sub-lethal effects of cadmium and copper on RNA/DNA ratio and energy reserves in the green-lipped mussel Perna viridis

This study aims to test if RNA/DNA ratio and various energy reserve parameters (i.e., glycogen, lipid, protein content and total energy reserves) are sensitive biomarkers for indicating stresses induced by metal contaminants in the green-lipped mussel Perna viridis, a common organism for biomonitoring in Southeast Asia. This study was, therefore, designed to examine the effects of cadmium (Cd) and copper (Cu) on these potential biomarkers in two major energy storage tissues, adductor muscle (AM) and hepatopancreas (HP), of P. viridis after sub-lethal exposure to either metal for 10 days. The results showed that neither Cd nor Cu treatments affected the RNA/DNA ratio, glycogen and protein contents in AM and HP. As the most sensitive biomarker in P. viridis, the total lipid content in both AM and HP was significantly decreased in the treatment of 5μg Cu/L and 0.01-0.1μgCd/L, respectively. However, soft-tissue body burdens of Cu or Cd did not significantly correlate with each of the four biomarkers regardless of the tissue type. In addition, AM generally stored more glycogen than HP, whereas HP stored more lipids than AM. We proposed that multiple biomarkers may be employed as an integrated diagnostic tool for monitoring the health condition of the mussels.

De novo transcriptome assembly of the marine gastropod Reishia clavigera for supporting toxic mechanism studies

The intertidal whelk Reishia clavigera is commonly used as a biomonitor of chemical contamination in the marine environment along Western Pacific region, and as a model for mechanistic studies of organotin-mediated imposex development. However, limited genomic resources of R. clavigera have restricted its role for the investigation of molecular mechanisms of such endocrine disruptions. This study, therefore, aimed to establish tissue-specific transcriptomes of the digestive gland, gonad, head ganglia, penis and the remaining body part of the male and female R. clavigera. By combining the results, a global transcriptome was obtained. A total of 578,134,720 high-quality filtered reads were obtained using Illumina sequencing. The R. clavigera transcriptome comprised of 38,466 transcripts and 32,798 unigenes with predicted open reading frames. The average length of transcripts was 1,709bp with N50 of 2,236bp. Based on sequence similarity searches against public databases, 28,657 transcripts and 24,403 unigenes had at least one BLAST hit. There were 17,530 transcripts and 14,897 unigenes annotated with at least one Gene Ontology (GO) term. Moreover, 5,776 transcripts and 5,137 unigenes were associated with 333 Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathways. The numbers of unigenes were similar among the five target tissues and between sexes, but tissue-specific expression profiles were revealed by multivariate analyses. Based on the functional annotation, putative steroid hormone-associated unigenes were identified. In particular, we highlighted the presence of steroid hormone receptor homologues that could be the targets for mechanistic studies of the organotin-mediated imposex development in marine gastropods. This newly generated transcriptome assembly of R. clavigera provides a valuable molecular resource for ecotoxicological and environmental genomic studies.

The Shell of the Invasive Bivalve Species Dreissena polymorpha: Biochemical, Elemental and Textural Investigations

The zebra mussel Dreissena polymorpha is a well-established invasive model organism. Although extensively used in environmental sciences, virtually nothing is known of the molecular process of its shell calcification. By describing the microstructure, geochemistry and biochemistry/proteomics of the shell, the present study aims at promoting this species as a model organism in biomineralization studies, in order to establish a bridge with ecotoxicology, while sketching evolutionary conclusions. The shell of D. polymorpha exhibits the classical crossed-lamellar/complex crossed lamellar combination found in several heterodont bivalves, in addition to an external thin layer, the characteristics of which differ from what was described in earlier publication. We show that the shell selectively concentrates some heavy metals, in particular uranium, which predisposes D. polymorpha to local bioremediation of this pollutant. We establish the biochemical signature of the shell matrix, demonstrating that it interacts with the in vitro precipitation of calcium carbonate and inhibits calcium carbonate crystal formation, but these two properties are not strongly expressed. This matrix, although overall weakly glycosylated, contains a set of putatively calcium-binding proteins and a set of acidic sulphated proteins. 2D-gels reveal more than fifty proteins, twenty of which we identify by MS-MS analysis. We tentatively link the shell protein profile of D. polymorpha and the peculiar recent evolution of this invasive species of Ponto-Caspian origin, which has spread all across Europe in the last three centuries.

Biodegradation of DDT by Stenotrophomonas sp. DDT-1: Characterization and genome functional analysis

A novel bacterium capable of utilizing 1,1,1-trichloro-2,2-bis(p-chlorophenyl)ethane (DDT) as the sole carbon and energy source was isolated from a contaminated soil which was identified as Stenotrophomonas sp. DDT-1 based on morphological characteristics, BIOLOG GN2 microplate profile, and 16S rDNA phylogeny. Genome sequencing and functional annotation of the isolate DDT-1 showed a 4,514,569 bp genome size, 66.92% GC content, 4,033 protein-coding genes, and 76 RNA genes including 8 rRNA genes. Totally, 2,807 protein-coding genes were assigned to Clusters of Orthologous Groups (COGs), and 1,601 protein-coding genes were mapped to Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway. The degradation half-lives of DDT increased with substrate concentration from 0.1 to 10.0 mg/l, whereas decreased with temperature from 15 °C to 35 °C. Neutral condition was the most favorable for DDT biodegradation. Based on genome annotation of DDT degradation genes and the metabolites detected by GC-MS, a mineralization pathway was proposed for DDT biodegradation in which it was orderly converted into DDE/DDD, DDMU, DDOH, and DDA via dechlorination, hydroxylation, and carboxylation, and ultimately mineralized to carbon dioxide. The results indicate that the isolate DDT-1 is a promising bacterial resource for the removal or detoxification of DDT residues in the environment.

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.

Potential Biotechnological Strategies for the Cleanup of Heavy Metals and Metalloids

Global mechanization, urbanization, and various natural processes have led to the increased release of toxic compounds into the biosphere. These hazardous toxic pollutants include a variety of organic and inorganic compounds, which pose a serious threat to the ecosystem. The contamination of soil and water are the major environmental concerns in the present scenario. This leads to a greater need for remediation of contaminated soils and water with suitable approaches and mechanisms. The conventional remediation of contaminated sites commonly involves the physical removal of contaminants, and their disposition. Physical remediation strategies are expensive, non-specific and often make the soil unsuitable for agriculture and other uses by disturbing the microenvironment. Owing to these concerns, there has been increased interest in eco-friendly and sustainable approaches such as bioremediation, phytoremediation and rhizoremediation for the cleanup of contaminated sites. This review lays particular emphasis on biotechnological approaches and strategies for heavy metal and metalloid containment removal from the environment, highlighting the advances and implications of bioremediation and phytoremediation as well as their utilization in cleaning-up toxic pollutants from contaminated environments.

High-throughput transcriptome sequencing of the cold seep mussel Bathymodiolus platifrons

Bathymodiolid mussels dominate hydrothermal vents, cold methane/sulfide-hydrocarbon seeps, and other sites of organic enrichment. Here, we aimed to explore the innate immune system and detoxification mechanism of the deep sea mussel Bathymodiolus platifrons collected from a methane seep in the South China Sea. We sequenced the transcriptome of the mussels’ gill, foot and mantle tissues and generated a transcriptomic database containing 96,683 transcript sequences. Based on GO and KEGG annotations, we reported transcripts that were related to the innate immune system, heavy metal detoxification and sulfide metabolic genes. Our in-depth analysis on the isoforms of peptidoglycan recognition protein (PGRP) that have different cellular location and potentially differential selectivity towards peptidoglycan (PGN) from gram-positive and gram-negative bacteria were differentially expressed in different tissues. We also reported a potentially novel form of metallothionein and the production of phytochelatin in B. platifrons, which has not been reported in any of its coastal relative Mytilus mussel species. Overall, the present study provided new insights into heavy metal and sulfide metabolism in B. platifrons and can be served as the basis for future molecular studies on host-symbiont interactions in cold seep mussels.

Differential gene expression analysis of benzo(a)pyrene toxicity in the clam, Ruditapes philippinarum

Polycyclic aromatic hydrocarbons (PAHs) are known for their carcinogenic, teratogenic and mutagenic properties. Benzo(a)pyrene (BaP) possesses the greatest carcinogenic potential among the various PAHs. In this study, digital gene expression (DGE) was performed to investigate the digestive gland transcriptome profile of the clam Ruditapes philippinarum exposed to BaP. A total of 10,508,312 and 11,414,297 clean reads were generated respectively, from control and BaP exposure DGE libraries. One hundred and forty-five differentially expressed genes were detected after comparing two libraries with 58 up-regulated and 87 down-regulated genes. GO annotation and KEGG pathway analyses were performed on all genes to understand their biological functions and processes. The results showed that numerous enriched differentially expressed genes are related to growth and development, antioxidant metabolism, apoptosis and detoxification metabolism. Quantitative real-time PCR was performed to verify the expressed genes of DGE. Our results provide evidences that RNA-seq is a powerful tool for toxicology and capable of generating novel and valuable information at the transcriptome level for characterizing deleterious effects caused by BaP.