Deep sequencing-based transcriptome profiling analysis of Chlamys farreri exposed to benzo[a]pyrene

Influence of temperature changes on oxidative stress and antioxidant defense system in the bay scallop, Argopecten irradians

In this study, we aimed to understand the effects of changes in temperature on biochemical and molecular responses associated with the antioxidant defense system in the bay scallop, Argopecten irradians. We measured the contents of H2O2 and malondialdehyde (MDA), as well as the activities of antioxidant enzymes (e.g., glutathione S-transferase [GST], superoxide dismutase [SOD], and catalase [CAT]), and the regulation of stress-related genes (e.g., GST, SOD, CAT, and heat shock protein 70 [HSP70]). In addition, total antioxidant capacity (TAC) was examined in scallops exposed to different temperatures. A. irradians showed high levels of H2O2 and MDA in response to acute thermal stress (48 and 72 h of exposure). Temperature changes also led to a significant increase in antioxidant enzyme activity and mRNA expression levels in A. irradians. Interestingly, the TAC increased in response to acute thermal stress (28 °C) for up to 12 h and decreased thereafter. The oxidative stress induced by high temperatures could not be alleviated by an increase in levels of antioxidant enzymes, such as GST, SOD, and CAT, resulting in high levels of H2O2 and MDA and low levels of TAC. In addition, significant changes (P < 0.05) in HSP70 levels were observed in response to changes in temperature, suggesting that HSP70 played an important role in the heat tolerance of A. irradians. In conclusion, A. irradians exhibits a greater degree of oxidative stress responses in high-temperature environments than that in low-temperature environments. Overall, these findings indicate that temperature changes lead to oxidative stress, resulting in cellular damage and activation of the antioxidant defense system in bay scallops. Further experiments are required to elucidate other antioxidants and fully understand the redox system in A. irradians.

Comparative toxicity study of waterborne two booster biocides (CuPT and ZnPT) on embryonic flounder (Paralichthys olivaceus)

A new generation of booster biocides that include metal pyrithiones (PTs) such as copper pyrithione (CuPT) and zinc pyrithione (ZnPT) are being used as tributyltin alternatives. In the marine environment, ZnPT can easily transchelate Cu to form CuPT, and the environmental fate and persistence of these two metal pyrithiones are closely related. Although some data on the toxicity of biocides on marine fish are available, little is known about their toxicity and toxic pathway. We thus compared the toxic effects of CuPT and ZnPT on embryonic olive flounder (Paralichthys olivaceus) by investigating their adverse effects based on developmental morphogenesis and transcriptional variation. In our study, the toxic potency of CuPT was greater with respect to developmental malformation and mortality than ZnPT. Consistent with the developmental effects, the expression of genes related to tail fin malformation (including plod2, furin, and wnt3a) was higher in embryonic flounder exposed to CuPT than in those exposed to ZnPT. Genes related to muscle and nervous system development exhibited significant changes on differential gene expression profiles using RNA sequencing (cutoff value P < 0.05). Gene ontology analysis of embryos exposed to CuPT revealed affected cellular respiration and kidney development, whereas genes associated with cell development, nervous system development and heart development showed significant variation in embryonic flounder exposed to ZnPT. Overall, our study clarifies the common and unique developmental toxic effects of CuPT and ZnPT through transcriptomic analyses in embryonic flounder.

Is hull cleaning wastewater a potential source of developmental toxicity on coastal non-target organisms?

Chemical contaminants can be discharged by vessel hull cleaning processes, such as scraping, jet spraying, and painting, all of which produce readily transportable contaminants into the marine environment, where they are referred to as 'hotspots' of contamination in coastal areas. However, many countries have not yet established effective evaluation methods for disposal of waste mixtures or management guidelines for areas of hull cleaning. To define the toxic effects of wastewater from vessel hull cleaning in dry docks on resident non-target organisms, we investigated the chemical concentrations and developmental toxicity on embryonic flounder, which is an organism sensitive to chemical contamination. In this study, the dominant inorganic metal discharged was zinc when cleaning Ship A (300 tons) and copper for Ship B (5,000 tons). The wastewater from high-pressure water blasting (WHPB) of Ship A (300 tons) and Ship B (5,000 tons) produced a largely overlapping suite of developmental malformations including pericardial edema, spinal curvature, and tail fin defects. Forty-eight hours after exposure, the frequency percentage of malformation began to increase in embryos exposed to a 500-fold dilution of WHPB from Ships A and B. We performed transcriptome sequencing to characterize the toxicological developmental effects of WHPB exposure at the molecular level. The results of the analysis revealed significantly altered expression of genes associated with muscle cell differentiation, actin-mediated cell contraction, and nervous system development (cutoff P < 0.01) in embryonic flounder exposed to high-pressure cleaning effluent from Ship A. Genes associated with chromatin remodeling, cell cycling, and insulin receptor signaling pathways were significantly altered in embryonic flounder exposed to WHPB of Ship B (cutoff P < 0.01). These findings provide a greater understanding of the developmental toxicity and potential effects of WHPB effluent on coastal embryonic fish. Furthermore, our results could inform WHPB effluent management practices to reduce impacts on non-target coastal organisms.

Transcriptomic changes in western mosquitofish (Gambusia affinis) liver following benzo[a]pyrene exposure

Widely distributed western mosquitofish (Gambusia affinis) has been used as a new model species for hazard assessment of environmental stressors such as polycyclic aromatic hydrocarbons (PAHs). However, most of the PAH studies using G. affinis rely on targeted biomarker-based analysis, and thus may not adequately address the complexity of the toxic mechanisms of the stressors. In the present study, the whole transcriptional sequencing of G. affinis liver after exposure to a PAH model, benzo[a]pyrene (BaP) (100 μg/L), for 20 days was performed by using the HiSeq XTen sequencers. In total, 58,156,233 and 51,825,467 clean nucleotide reads were obtained in the control and BaP-exposed libraries, respectively, with average N50 lengths of 1419 bp. In addition, after G. affinis was exposed for 20 days, 169 genes were upregulated, and 176 genes were downregulated in liver. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were applied to all the genes to determine the genes' biological functions and processes. The results clearly showed that the differentially expressed genes were mainly related to immune pathways and metabolic correlation pathways. Interestingly, almost all the pathways related with the immunity were upregulated, while the metabolism pathways were downregulated. Lastly, quantitative real-time PCR (qRT-PCR) was performed to measure expressional levels of twelve genes confirmed through the DGE analysis. These results demonstrate that BaP damages immunity and enhances the consumption of all available energy storage to activate mechanisms of the detoxification in G. affinis. Up until now, the present study is the first time that a whole transcriptome sequencing analysis in the liver of G. affinis exposed to BaP has been reported.

Differential responses in the biotransformation systems of the oyster Crassostrea gasar (Adanson, 1757) elicited by pyrene and fluorene: molecular, biochemical and histological approach - Part I

Polycyclic aromatic hydrocarbons (PAHs) are among the main contaminants in aquatic environments. PAHs can affect organisms due to their carcinogenic, mutagenic and/or teratogenic characteristics. Depending on the PAHs, concentration, and period of exposure, biological damage can occur leading to histopathologic alterations. This study aimed to evaluate the molecular, biochemical and histological responses of the oyster Crassostrea gasar exposed to pyrene (0.25 and 0.5 μM) and fluorene (0.6 and 1.2 μM), after exposure for 24 and 96 h. Concentrations of both PAHs were quantified in the water and in oyster tissues. Transcript levels of phase I (CYP3475C1, CYP2-like, CYP2AU1 and CYP356A) and phase II (GSTO-like, MGST-like and SULT-like) biotransformation-related genes and the activities of ethoxyresorufin-O-deethylase (EROD), total and microsomal glutathione S-transferase (GST and MGST) were evaluated in the gills. Also, histological changes and localization of mRNA transcripts CYP2AU1 in gills, mantle, and digestive diverticula were evaluated. Both PAHs accumulated in oyster tissues. Pyrene half-life in water was significantly lower than fluorene. Transcript levels of all genes were higher in oysters exposed to of pyrene 0.5 μM (24 h). Only CYP2AU1 gene was up-regulated by fluorene exposure. EROD and MGST activities were higher in oysters exposed to pyrene. Tubular atrophy in the digestive diverticula and an increased number of mucous cells in the mantle were observed in oysters exposed to pyrene. CYP2AU1 transcripts were observed in different tissues of pyrene-exposed oysters. A significant correlation was observed between tubular atrophy and the CYP2AU1 hybridization signal in oysters exposed to pyrene, suggesting the sensibility of the species to this PAH. These results suggest an important role of biotransformation-related genes and enzymes and tissue alterations associated to pyrene metabolism but not fluorene. In addition, it reinforces the role of CYP2AU1 gene in the biotransformation process of PAHs in the gills of C. gasar.

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.

Overlapping and unique toxic effects of three alternative antifouling biocides (Diuron, Irgarol 1051®, Sea-Nine 211®) on non-target marine fish

The use of alternative biocides has increased due to their economic and ecological relevance. Although data regarding the toxicity of commercial alternative biocides in marine organisms are accumulating, little is known about their toxic pathways or mechanisms. To compare the toxic effects of commercial alternative biocides on non-target pelagic fish (flounder) embryos, we investigated the adverse effects of developmental malformation and transcriptional changes. Three biocides including Diuron, Irgarol 1051^® and Sea-Nine 211^® produced a largely overlapping suite of developmental malformations, including tail-fin fold defects and dorsal body axis curvature. In our test, the potencies of these biocides were ranked in the following order with respect to malformation and mortalities: Sea-Nine 211^® > Irgarol 1051^® > Diuron. Consistent with the toxicity rankings, the expression of genes related to heart formation was greater in embryonic flounder exposed to Sea-Nine 211^® than in those exposed to Irgarol 1051^® or Diuron, while expression of genes related to fin malformation was greater in the Irgarol 1051^® exposure group. In analyses of differential gene expression (DEG) profiles (fold change of genes with a cutoff P < 0.05) by high-throughput sequencing (RNA-seq), genes associated with nervous system development, transmembrane transport activity, and muscle cell development were significantly changed commonly. Embryos exposed to Diuron showed changes related to cellular protein localization, whereas genes associated with immune system processes were up-regulated significantly in embryos exposed to Irgarol 1051^®. Genes related to actin filament organization and embryonic morphogenesis were up-regulated in embryos exposed to Sea-Nine 211^®. Overall, our study provides a better understanding of the overlapping and unique developmental toxic effects of three commercial booster biocides through transcriptomic analyses in a non-target species, embryonic flounder.

Detoxification, Apoptosis, and Immune Transcriptomic Responses of the Gill Tissue of Bay Scallop Following Exposure to the Algicide Thiazolidinedione 49

Thiazolidinedione 49 (TD₄₉), a newly synthesized algicide, shows strong toxicity at low concentrations of 0.1-2.0 μM. However, its potential effects on non-target species at the transcript level were not investigated. Differentially expressed genes (DEGs) in the gills of the bay scallop, Argopecten irradians, were accessed after treatment with 0.68 μM TD₄₉ for up to 48 h. Following exposure, it was observed that 5214 genes were upregulated and 3497 were downregulated. Functional enrichment analysis revealed that the apoptosis pathway was activated. The extrinsic apoptosis pathway was activated and the survival factors related pathway was suppressed. Furthermore, gene expressions related to ATP-binding cassette, nuclear factor erythroid 2-related factor, B cell lymphoma-2 family protein, glutathione reductase, glutathione peroxidase, catalase, NADPH2:quinone reductase, and superoxide dismutase were decreased. Conversely, gene expressions related to FAS-associated death domain protein, glutathione S-transferase, caspase 6, 8, cytochrome P450 1A1, and 2C8 were increased. These results comprehensively demonstrated the toxicity of the novel algicide TD₄₉, and should draw the attention of researchers to the importance of analyzing the potential impact of chemical compounds as algicides to control the proliferation of harmful algae, due to the secondary pollution caused by their application.

Detoxification and Immune Transcriptomic Response of the Gill Tissue of Bay Scallop (Argopecten irradians) Following Exposure to the Algicide Palmitoleic Acid

Palmitoleic acid (PA) is an effective algicide against Alexandrium tamarense. However, the toxicological mechanism of PA exposure is unclear. The transcript abundance and differentially expressed genes (DEGs) in gills of bay scallop were investigated following 80 mg/L PA exposure up to 48 h using the Illumina HiSeq 4000 deep-sequencing platform with the recommended read length of 100 bp. De novo assembly of paired-end reads yielded 62,099 unigenes; 5414 genes were identified as being significantly increased, and 4452 were decreased. Based on gene ontology classification and enrichment analysis, the ‘cellular process’, ‘metabolic process’, ‘response to stimulus’, and ‘catalytic process’ with particularly high functional enrichment were revealed. The DEGs, which are related to detoxification and immune responses, revealed that acid phosphatase, fibrinogen C domain-containing protein, cyclic AMP-responsive element-binding protein, glutathione reductase, ATP-binding cassette, nuclear factor erythroid 2-related factor, NADPH2:quinone reductase, and cytochrome P450 4F22, 4B1, and 2C8-related gene expression decreased. In contrast, some genes related to glutathione S-transferase, C-type lectin, superoxide dismutase, toll-like receptors, and cytochrome P450 2C14, 2U1, 3A24 and 4A2 increased. The results of current research will be a valuable resource for the investigation of gene expression stimulated by PA, and will help understanding of the molecular mechanisms underlying the scallops’ response to PA exposure.

Transcriptomic response of the Crassostrea virginica gonad after exposure to a water-accommodation fraction of hydrocarbons and the potential implications in reproduction

The Crassostrea virginica oyster has biological and economic importance in the Gulf of Mexico, an area with a high extraction and production of hydrocarbons. Exposure to hydrocarbons affects the reproductive processes in bivalves. In C. virginica, the effect of hydrocarbons on the gonad of the undifferentiated organism has not been evaluated to determine the possible damage during the maturation process. To evaluate this effect, RNA-seq data was generated from C. virginica gonads exposed to a 200 μg/L of hydrocarbons at different exposure times (7, 14 and 21 days) and a control treatment (without hydrocarbons). The analysis of the gonad transcriptome showed the negative effect of hydrocarbons on maturation, with a sub-expression of 22 genes involved in different stages of this process. Additionally, genes in the immune system were down-regulated, which may indicate that exposure to hydrocarbons causes immunosuppression in bivalves. A group of oxidative stress genes was also reduced. These data contribute to a better understanding of the effect of hydrocarbons on the reproductive process in bivalves and, at the same time, allow us to identify possible biomarkers associated with hydrocarbon contamination in the gonad of C. virginica.

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

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

Transcriptional responses of Acropora hyacinthus embryo under the benzo(a)pyrene stress by deep sequencing

Coral embryos are a critical and sensitive period for the early growth and development of coral. Benzo(a)pyrene (BaP) is widely distributed in the ocean and has strong toxicity, but there is little information on the toxic effects to coral embryos exposed to this widespread environmental contaminant. Thus, in this study, we utilized the Illumina Hiseq™ 4000 platform to explore the gene response of Acropora hyacinthus embryos under the BaP stress. A total of 130,042 Unigenes were obtained and analyzed, and approximately 37.67% of those matched with sequences from four different species. In total, 2606 Unigenes were up-regulated, and 3872 Unigenes were down-regulated. After Gene Ontology (GO) annotation, the results show that the "cellular process" and "metabolic process" were leading in the category of biological processes, which the "binding" and "catalytic activity" were the most abundant subcategories in molecular function. Based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, the most differentially expressed genes (DEGs) were enriched, as well as down-regulated in the pathways of oxidative phosphorylation, metabolism of xenobiotics, immune-related genes, apoptosis and human disease genes. At the same time, 388,197 of Single-nucleotide Polymorphisms (SNPs) and 6164 of Simple Sequence Repeats (SSRs) were obtained, which can be served as the richer and more valuable SSRs molecular markers in the future. The results of this study can help to better understand the toxicological mechanism of coral embryo exposed to BaP, and it is also essential for the protection and restoration of coral reef ecosystem in the future.

Detoxification- and Immune-Related Transcriptomic Analysis of Gills from Bay Scallops (Argopectenirradians) in Response to Algal Toxin Okadaic Acid

To reveal the molecular mechanisms triggered by okadaic acid (OA)-exposure in the detoxification and immune system of bay scallops, we studied differentially-expressed genes (DEGs) and the transcriptomic profile in bay scallop gill tissue after 48 h exposure to 500 nM of OA using the Illumina HiSeq 4000 deep-sequencing platform. De novo assembly of paired-end reads yielded 55,876 unigenes, of which 3204 and 2620 genes were found to be significantly up- or down-regulated, respectively. Gene ontology classification and enrichment analysis of the DEGs detected in bay scallops exposed to OA revealed four ontologies with particularly high functional enrichment, which were ‘cellular process’ (cellular component), ‘metabolic process’ (biological process), ‘immune system process’ (biological process), and ‘catalytic process’ (molecular function). The DEGs revealed that cyclic AMP-responsive element-binding proteins, acid phosphatase, toll-like receptors, nuclear erythroid 2-related factor, and the NADPH2 quinone reductase-related gene were upregulated. In contrast, the expression of some genes related to glutathione S-transferase 1, C-type lectin, complement C1q tumor necrosis factor-related protein, Superoxide dismutase 2 and fibrinogen C domain-containing protein, decreased. The outcomes of this study will be a valuable resource for the study of gene expression induced by marine toxins, and will help understanding of the molecular mechanisms underlying the scallops’ response to OA exposure.

The various components implied the diversified Toll-like receptor (TLR) signaling pathway in mollusk Chlamys farreri

Toll-like receptor (TLR) signaling pathway, composed of various components, plays pivotal roles in host innate immune defense mechanism. In the present study, twenty-nine TLR signaling pathway components, including receptors, adaptors, transduction molecules and immune effectors, were identified in Zhikong scallop Chlamys farreri via assembling and screening public available transcriptomic data and expression sequence tags (ESTs). These identified TLR signaling pathway components were constitutively expressed and detectable in various tissues, and almost all of them were highly expressed in gill and hepatopancreas. These results indicated the presence of TLR signaling pathways in both MyD88-dependent and MyD88-independent forms in scallop, and implied the diversified TLR signaling pathway in mollusk C. farreri.

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

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

RNA seq- and DEG-based comparison of developmental toxicity in fish embryos of two species exposed to Iranian heavy crude oil

To determine and compare the toxic effects of Iranian heavy crude oil (IHCO) on the embryonic development of two fish species, we examined transcriptome profiles using RNA-seq. The assembled contigs were 66,070 unigenes in olive flounder embryos and 76,498 unigenes in spotted seabass embryos. In the differential gene expression (DEG) profiles, olive flounder embryos showed different up- and down-regulated patterns than spotted seabass embryos in response to fresh IHCO (FIHCO) and weathered IHCO (WIHCO). In this work, we categorized DEG profiles into six pathways: ribosome, oxidative phosphorylation, Parkinson's disease, Alzheimer's disease, Huntington's disease, and cardiac muscle contraction, validating the expression patterns of 13 DEGs using real-time quantitative RT-PCR. The expression of the CYP1A, CYP1B1, and CYP1C1 genes in spotted seabass embryos was higher than in olive flounder embryos, whereas genes related to cell processing, development, and the immune system showed the opposite trend. Orthologous gene cluster analysis showed that olive flounder embryos were sensitive (fold change of genes with cutoff P<0.05) to both FIHCO and WIHCO, but spotted seabass embryos exhibited higher sensitivity to WIHCO than FIHCO, indicating that species-specific differences are likely to be reflected in population levels after oil spills. Overall, our study provides new insight on the different embryonic susceptibilities of two marine fish species to FIHCO and WIHCO and a better understanding of the underlying molecular mechanisms via RNA-seq and DEGs.

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

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

Combined exposure to pyrene and fluoranthene and their molecular effects on the Sydney rock oyster, Saccostrea glomerata

Polycyclic aromatic hydrocarbons (PAHs) are ubiquitously detected in the water column, associated with particulate matter or in the tissue of marine organisms such as molluscs. PAH exposure and their resultant bioaccumulation in molluscs can cause a range of serious physiological effects in the affected animals. To examine the molecular response of these xenobiotics in bivalves, Sydney rock oysters (Saccostrea glomerata) were exposed to pyrene and fluoranthene for seven days. Chemical analysis of the soft-tissue of PAH stressed S. glomerata confirmed that pyrene and fluoranthene could be bioaccumulated by these oysters. RNA-Seq analysis of PAH-exposed S. glomerata showed a total of 765 transcripts differentially expressed between control and PAH-stressed oysters. Closer examination of the transcripts revealed a range genes encoding enzymes involved in PAH detoxification (e.g. cytochrome P450), innate immune responses (e.g. pathogen recognition, phagocytosis) and protein synthesis. Overall, pyrene and fluoranthene exposure appears to have resulted in a suppression of pathogen recognition and some protein synthesis processes, whereas transcripts of genes encoding proteins involved in clearance of cell debris and some transcripts of genes involved in PAH detoxification were induced in response to the stressors. Pyrene and fluoranthene exposure thus invoked a complex molecular response in S. glomerata, with results suggesting that oysters focus on removing the stressors from their system and dealing with the downstream effects of PAH exposure, potentially at the exclusion of other, less immediate concerns (e.g. protection from infection).

Comparative Analysis and Distribution of Omega-3 lcPUFA Biosynthesis Genes in Marine Molluscs

Recent research has identified marine molluscs as an excellent source of omega-3 long-chain polyunsaturated fatty acids (lcPUFAs), based on their potential for endogenous synthesis of lcPUFAs. In this study we generated a representative list of fatty acyl desaturase (Fad) and elongation of very long-chain fatty acid (Elovl) genes from major orders of Phylum Mollusca, through the interrogation of transcriptome and genome sequences, and various publicly available databases. We have identified novel and uncharacterised Fad and Elovl sequences in the following species: Anadara trapezia, Nerita albicilla, Nerita melanotragus, Crassostrea gigas, Lottia gigantea, Aplysia californica, Loligo pealeii and Chlamys farreri. Based on alignments of translated protein sequences of Fad and Elovl genes, the haeme binding motif and histidine boxes of Fad proteins, and the histidine box and seventeen important amino acids in Elovl proteins, were highly conserved. Phylogenetic analysis of aligned reference sequences was used to reconstruct the evolutionary relationships for Fad and Elovl genes separately. Multiple, well resolved clades for both the Fad and Elovl sequences were observed, suggesting that repeated rounds of gene duplication best explain the distribution of Fad and Elovl proteins across the major orders of molluscs. For Elovl sequences, one clade contained the functionally characterised Elovl5 proteins, while another clade contained proteins hypothesised to have Elovl4 function. Additional well resolved clades consisted only of uncharacterised Elovl sequences. One clade from the Fad phylogeny contained only uncharacterised proteins, while the other clade contained functionally characterised delta-5 desaturase proteins. The discovery of an uncharacterised Fad clade is particularly interesting as these divergent proteins may have novel functions. Overall, this paper presents a number of novel Fad and Elovl genes suggesting that many mollusc groups possess most of the required enzymes for the synthesis of lcPUFAs.

The use of -omic tools in the study of disease processes in marine bivalve mollusks

Our understanding of disease processes and host-pathogen interactions in model species has benefited greatly from the application of medium and high-throughput genomic, metagenomic, epigenomic, transcriptomic, and proteomic analyses. The rate at which new, low-cost, high-throughput -omic technologies are being developed has also led to an expansion in the number of studies aimed at gaining a better understanding of disease processes in bivalves. This review provides a catalogue of the genetic and -omic tools available for bivalve species and examples of how -omics has contributed to the advancement of marine bivalve disease research, with a special focus in the areas of immunity, bivalve-pathogen interactions, mechanisms of disease resistance and pathogen virulence, and disease diagnosis. The analysis of bivalve genomes and transcriptomes has revealed that many immune and stress-related gene families are expanded in the bivalve taxa examined thus far. In addition, the analysis of proteomes confirms that responses to infection are influenced by epigenetic, post-transcriptional, and post-translational modifications. The few studies performed in bivalves show that epigenetic modifications are non-random, suggesting a role for epigenetics in regulating the interactions between bivalves and their environments. Despite the progress -omic tools have enabled in the field of marine bivalve disease processes, there is much more work to be done. To date, only three bivalve genomes have been sequenced completely, with assembly status at different levels of completion. Transcriptome datasets are relatively easy and inexpensive to generate, but their interpretation will benefit greatly from high quality genome assemblies and improved data analysis pipelines. Finally, metagenomic, epigenomic, proteomic, and metabolomic studies focused on bivalve disease processes are currently limited but their expansion should be facilitated as more transcriptome datasets and complete genome sequences become available for marine bivalve species.

RNA-seq based on transcriptome reveals differ genetic expressing in Chlamys farreri exposed to carcinogen PAHs

The effects of a mixture of carcinogen polycyclic aromatic hydrocarbons (BaP, BaA, BbF and CHR) on transcriptional responses in the digestive gland of scallop, Chlamys farreri, were investigated by high-throughput RNA-sequencing (RNA-seq) technologies. In total, 209 and 260 genes were identified as significantly up- or down-regulated, respectively. Functional analysis based on gene ontology (GO) classification system and the Kyoto encyclopedia of genes and genomes (KEGG) database revealed that PAHs significantly altered the expression of genes involved in stress response, detoxication, antioxidation which were extensively discussed. In particular, CYP450 2P1 and QO mRNA expression were found to be up-regulated by exposure to PAHs mixture, suggesting that CYP450 2P1 and QO mRNA expression can be a potential marker for prediction of the biological effects of a mixture of PAHs on scallops.