Solute carriers in scallop genome: Gene expansion and expression regulation after exposure to toxic dinoflagellate

Comprehensive Multi-omics Approaches Provide Insights to Summer Mortality in the Clam Meretrix petechialis

Bivalve mass mortalities have been reported worldwide, which not only can be explained as a result of pathogen infection, but may reflect changes in environments. Although these episodes were often reported, there was limited information concerning the molecular responses to various stressors leading to summer mortality. In the present work, RNA sequencing (RNA-seq), tandem mass tagging (TMT)-based quantitative proteomics, and 16S rRNA sequencing were used to explore the natural outbreak of summer mortality in the clam Meretrix petechialis. We identified a total of 172 differentially expressed genes (DEGs) and 222 differentially expressed proteins (DEPs) in the diseased group compared to the normal group. The inconsistent expression profiles of immune DEGs/DEPs may be due to the immune dysregulation of the diseased clams. Notably, 11 solute carrier family genes were found among the top 20 down-regulated genes in the diseased group, indicating that weakened transmembrane transport ability might occur in the diseased clams. Integration analysis of transcriptomic and proteomic results showed that many metabolic processes such as "arginine and proline metabolism" and "tyrosine metabolism" were inhibited in the diseased group, suggesting metabolic inhibition. Moreover, 16S rRNA sequencing revealed that the microbial composition of clam hepatopancreas was disordered in the diseased group. The comparison of DEGs expression between the natural summer mortality event and an artificial challenge experiment involving both Vibrio infection and heat stress revealed 9/15 genes showing similar expression trends between the two conditions, suggesting that the summer mortality might be caused by a combination of high temperature and Vibrio infection. These results would deepen our understanding of summer mortality and provide candidate resistance markers for clam resistance breeding.

Tissue-specific antioxidative response and metabolism of paralytic shellfish toxins in scallop (Chlamys farreri) mantle with Alexandrium dinoflagellate exposure

Bivalves show remarkable capacity to acclimate paralytic shellfish toxins (PSTs) produced by dinoflagellates, severely affecting fishery industry and public health. Here, transcriptomic response to PSTs-producing dinoflagellate (Alexandrium minutum) was investigated in Zhikong scallop (Chlamys farreri) mantle. The PSTs accumulated in C. farreri mantle continually increased during the 15 days exposure, with "oxidation-reduction" genes induced compared to the control group at the 1st and 15th day. Through gene co-expression network analysis, 16 PSTs-responsive modules were enriched with up- or down-regulated genes. The concentration of GTXs, major PSTs in A. minutum and accumulated in scallops, was correlated with the up-regulated magenta module, enriching peroxisome genes as the potential mantle-specific PSTs biomarker. Moreover, Hsp70B2s were inhibited throughout the exposure, which together with the expanded neurotransmitter transporter SLC6As, may play essential roles on neurotransmitter homeostasis in scallop mantle. These results paved the way for a comprehensive understanding of defensive mechanism and homeostatic response in scallop mantle against PSTs.

Peripheral and central employment of acid-sensing ion channels during early bilaterian evolution

Nervous systems are endowed with rapid chemosensation and intercellular signaling by ligand-gated ion channels (LGICs). While a complex, bilaterally symmetrical nervous system is a major innovation of bilaterian animals, the employment of specific LGICs during early bilaterian evolution is poorly understood. We therefore questioned bilaterian animals' employment of acid-sensing ion channels (ASICs), LGICs that mediate fast excitatory responses to decreases in extracellular pH in vertebrate neurons. Our phylogenetic analysis identified an earlier emergence of ASICs from the overarching DEG/ENaC (degenerin/epithelial sodium channel) superfamily than previously thought and suggests that ASICs were a bilaterian innovation. Our broad examination of ASIC gene expression and biophysical function in each major bilaterian lineage of Xenacoelomorpha, Protostomia, and Deuterostomia suggests that the earliest bilaterian ASICs were probably expressed in the periphery, before being incorporated into the brain as it emerged independently in certain deuterostomes and xenacoelomorphs. The loss of certain peripheral cells from Ecdysozoa after they separated from other protostomes likely explains their loss of ASICs, and thus the absence of ASICs from model organisms Drosophila and Caenorhabditis elegans. Thus, our use of diverse bilaterians in the investigation of LGIC expression and function offers a unique hypothesis on the employment of LGICs in early bilaterian evolution.

Transcriptomic analysis reveals the genes involved in tetrodotoxin (TTX) accumulation, translocation, and detoxification in the pufferfish Takifugu rubripes

Tetrodotoxin (TTX) is a potent marine neurotoxin that exists in a variety of aquatic and terrestrial organisms. Pufferfish in different habitats show great variation in their TTX contents. Exploring the genes involved in TTX metabolism could contribute to our understanding of the molecular mechanisms underlying TTX accumulation, translocation, and detoxification in pufferfish. In this study, transcriptomic analysis was used to identify the functional genes related to TTX metabolism in the blood, liver, and muscle of the toxic and non-toxic tiger puffer (Takifugu rubripes). A total of 6101 differentially expressed genes (DEGs) were obtained after transcriptomic analysis; of these, 2401 were identified in the blood, 2262 in the liver, and 1438 in the muscle. After enrichment analysis, fourteen genes encoding glutathione S-transferases (GSTs), glutathione peroxidase (GPx), thioredoxins (TXNs), superoxide dismutase (SOD), ATP-binding cassettes (ABCs), apolipoproteins (APOs), inhibitors of apoptosis protein (IAP), and solute carrier (SLC), which are mainly antioxidant enzymes, membrane transporters, or anti-apoptotic factors, were revealed in the blood. Thirty-six genes encoding SLCs, ABCs, long-chain-fatty-acid-CoA ligases (ACSLs), interleukin 6 cytokine family signal transducer (IL6ST), endoplasmic reticulum (ER), and heat shock protein family A (Hsp70) were involved in transmembrane transporter activity and innate immune response. Notably, a large number of slc genes were found to play critical and diverse roles in TTX accumulation and translocation in the liver of T. rubripes. Nine genes from the slc, hsp70, complement C5 (c5), acsl, er, and serpin peptidase inhibitor (serpin) gene families were found to participate in the regulation of protein processing and anti-apoptosis. These results reflect the diverse functions of genes closely related to TTX accumulation, translocation, and detoxification in T. rubripes.

Transcriptome and Network Analyses Reveal the Gene Set Involved in PST Accumulation and Responses to Toxic Alexandrium minutum Exposure in the Gills of Chlamys farreri

Bivalve molluscs are filter-feeding organisms that can accumulate paralytic shellfish toxins (PST) through ingesting toxic marine dinoflagellates. While the effects of PST accumulation upon the physiology of bivalves have been documented, the underlying molecular mechanism remains poorly understood. In this study, transcriptomic analysis was performed in the gills of Zhikong scallop (Chlamys farreri) after 1, 3, 5, 10, and 15 day(s) exposure of PST-producing dinoflagellate Alexandrium minutum. Higher numbers of differentially expressed genes (DEGs) were detected at day 1 (1538) and day 15 (989) than that at day 3 (77), day 5 (82), and day 10 (80) after exposure, and most of the DEGs were only regulated at day 1 or day 15, highlighting different response mechanisms of scallop to PST-producing dinoflagellate at different stages of exposure. Functional enrichment results suggested that PST exposure induced the alterations of nervous system development processes and the activation of xenobiotic metabolism and substance transport processes at the acute and chronic stages of exposure, respectively, while the immune functions were inhibited by PST and might ultimately cause the activation of apoptosis. Furthermore, a weighted gene co-expression network was constructed, and ten responsive modules for toxic algae exposure were identified, among which the yellow module was found to be significantly correlated with PST content. Most of the hub genes in the yellow module were annotated as solute carriers (SLCs) with eight being OCTN1s, implying their dominant roles in regulating PST accumulation in scallop gills. Overall, our results reveal the gene set responding to and involved in PST accumulation in scallop gills, which will deepen our understanding of the molecular mechanism of bivalve resistance to PST.

Marine Neurotoxins' Effects on Environmental and Human Health: An OMICS Overview

Harmful algal blooms (HAB), and the consequent release of toxic metabolites, can be responsible for seafood poisoning outbreaks. Marine wildlife can accumulate these toxins throughout the food chain, which presents a threat to consumers’ health. Some of these toxins, such as saxitoxin (STX), domoic acid (DA), ciguatoxin (CTX), brevetoxin (BTX), tetrodotoxin (TTX), and β-N-methylamino-L-alanine (BMAA), cause severe neurological symptoms in humans. Considerable information is missing, however, notably the consequences of toxin exposures on changes in gene expression, protein profile, and metabolic pathways. This information could lead to understanding the consequence of marine neurotoxin exposure in aquatic organisms and humans. Nevertheless, recent contributions to the knowledge of neurotoxins arise from OMICS-based research, such as genomics, transcriptomics, proteomics, and metabolomics. This review presents a comprehensive overview of the most recent research and of the available solutions to explore OMICS datasets in order to identify new features in terms of ecotoxicology, food safety, and human health. In addition, future perspectives in OMICS studies are discussed.

Proteome Response of Meretrix Bivalves Hepatopancreas Exposed to Paralytic Shellfish Toxins Producing Dinoflagellate Gymnodinium catenatum

Paralytic shellfish toxins (PSTs) contamination of seafood has become a growing global problem. However, the molecular response of bivalves, some of the most popular seafoods, to PSP toxins has seldom been reported and the underlying molecular mechanisms of the interactions between Meretrix meretrix bivalves and PSTs-producing dinoflagellates are scarcely known. This study compared the protein expression profiles between PSP toxin-contaminated and non-PSP toxin contaminated M. meretrix, determined proteome responses and identified potential biomarkers based on feeding experiments. Results showed that the content of total PSP toxins in contaminated bivalves was 40.63 ± 4.08 μg saxitoxin (STX) equivalents per gram, with 95.3% in hepatopancreas, followed by gill (1.82%) and foot (1.79%). According to two-dimensional gel electrophoresis (2-DE), 15 differentially expressed proteins (at least 2-fold difference) between the hepatopancreas of bivalves with and without PSP toxins were detected. Eight of them were successfully identified by MALDI-TOF MS. These were catalase, protein ultraspiracle homolog, G2 and S phase-expression protein, paramyosin, Mn-superoxide dismutase, response regulator receiver domain-containing protein, sarcoplasmic calcium-binding protein and major facilitator superfamily transporters. The differences in the expression levels of the last three proteins involving in cell signaling, structure and membrane transport were 4.2, 5.3 and 4.9-fold, respectively. These proteins could be further developed as potential biomarkers. The other two up-regulated proteins, Mn-superoxide dismutase and catalase, were involved in cell defence mechanisms against oxidative stress, suggesting PSP toxin acts as xenobiotics and poses oxidative stress in bivalves. This study gives insights into the response of bivalves to PSP toxin-producing dinoflagellate at the proteomic level and the potential of using 2-DE to develop specific protein markers in bivalves.

Transcriptional Response in the Digestive Gland of the King Scallop (Pecten maximus) After the Injection of Domoic Acid

Some diatom species of the genus Pseudo-nitzschia produce the toxin domoic acid. The depuration rate of domoic acid in Pecten maximus is very low; for this reason, king scallops generally contain high levels of domoic acid in their tissues. A transcriptomic approach was used to identify the genes differentially expressed in the P. maximus digestive gland after the injection of domoic acid. The differential expression analysis found 535 differentially expressed genes (226 up-regulated and 309 down-regulated). Protein–protein interaction networks obtained with the up-regulated genes were enriched in gene ontology terms, such as vesicle-mediated transport, response to stress, signal transduction, immune system process, RNA metabolic process, and autophagy, while networks obtained with the down-regulated genes were enriched in gene ontology terms, such as response to stress, immune system process, ribosome biogenesis, signal transduction, and mRNA processing. Genes that code for cytochrome P450 enzymes, glutathione S-transferase theta-1, glutamine synthase, pyrroline-5-carboxylate reductase 2, and sodium- and chloride-dependent glycine transporter 1 were among the up-regulated genes. Therefore, a stress response at the level of gene expression, that could be caused by the domoic acid injection, was evidenced by the alteration of several biological, cellular, and molecular processes.

IAPs Gene Expansion in the Scallop Patinopecten yessoensis and Their Expression Profiles After Exposure to the Toxic Dinoflagellate

Inhibitors of apoptosis proteins (IAPs) are conserved regulators involved in cell cycle, cell migration, cell death, immunity and inflammation, should be due to the fact that they can assist with the ability to cope with different kinds of extrinsic or intrinsic stresses. Bivalve molluscs are well adapted to highly complex marine environments. As free-living filter feeders that may take toxic dinoflagellates as food, bivalves can accumulate and put up with significant levels of paralytic shellfish toxins (PSTs). PSTs absorption and accumulation could have a deleterious effect on bivalves, causing negative impact on their feeding and digestion capabilities. In the present study, we analyzed IAP genes (PyIAPs) in Yesso scallop (Patinopecten yessoensis), a major fishery and aquaculture species in China. Forty-seven PyIAPs from five sub-families were identified, and almost half of the PyIAP genes were localized in clusters on two chromosomes. Several sites under positive selection was revealed in the significantly expanded sub-families BIRC4 and BIRC5. After exposure to PST-producing dinoflagellates, Alexandrium catenella, fourteen PyIAPs showed significant responses in hepatopancreas and kidney, and more than eighty-five percent of them were from the expanded sub-families BIRC4 and BIRC5. The regulation pattern of PyIAPs was similar between the two tissues, with more than half exhibited expression suppression within three days after exposure. In contrast to hepatopancreas, more acute changes of PyIAPs expression could be detected in kidney, suggesting the possible involvement of these PyIAPs in tissue-specific PST tolerance. These findings also imply the adaptive expansion of bivalve IAP genes in response to algae derived biotoxins.

Molecular Identification and Toxin Analysis of Alexandrium spp. in the Beibu Gulf: First Report of Toxic A. tamiyavanichii in Chinese Coastal Waters

The frequency of harmful algal blooms (HABs) has increased in China in recent years. Information about harmful dinoflagellates and paralytic shellfish toxins (PSTs) is still limited in China, especially in the Beibu Gulf, where PSTs in shellfish have exceeded food safety guidelines on multiple occasions. To explore the nature of the threat from PSTs in the region, eight Alexandrium strains were isolated from waters of the Beibu Gulf and examined using phylogenetic analyses of large subunit (LSU) rDNA, small subunit (SSU) rDNA, and internal transcribed spacer (ITS) sequences. Their toxin composition profiles were also determined using liquid chromatography-tandem mass spectrometry (LC-MS/MS). All eight strains clustered in the phylogenetic tree with A. pseudogonyaulax, A. affine, and A. tamiyavanichii from other locations, forming three well-resolved groups. The intraspecific genetic distances of the three Alexandrium species were significantly smaller than interspecific genetic distances for Alexandrium species. Beibu Gulf isolates were therefore classified as A. pseudogonyaulax, A. affine, and A. tamiyavanichii. No PSTs were identified in A. pseudogonyaulax, but low levels of gonyautoxins (GTXs) 1 to 5, and saxitoxin (STX) were detected in A. tamiyavanichii (a total of 4.60 fmol/cell). The extremely low level of toxicity is inconsistent with PST detection above regulatory levels on multiple occasions within the Beibu Gulf, suggesting that higher toxicity strains may occur in those waters, but were unsampled. Other explanations including biotransformation of PSTs in shellfish and the presence of other PST-producing algae are also suggested. Understanding the toxicity and phylogeny of Alexandrium species provides foundational data for the protection of public health in the Beibu Gulf region and the mitigation of HAB events.

Genome-Wide Identification and Characterization of SODs in Zhikong Scallop Reveals Gene Expansion and Regulation Divergence after Toxic Dinoflagellate Exposure

As filter-feeding animals mainly ingesting microalgae, bivalves could accumulate paralytic shellfish toxins (PSTs) produced by harmful algae through diet. To protect themselves from the toxic effects of PSTs, especially the concomitant oxidative damage, the production of superoxide dismutase (SOD), which is the only eukaryotic metalloenzyme capable of detoxifying superoxide, may assist with toxin tolerance in bivalves. To better understand this process, in the present study, we performed the first systematic analysis of SOD genes in bivalve Chlamys farreri, an important aquaculture species in China. A total of six Cu/Zn-SODs (SOD1-6) and two Mn-SODs (SOD7, SOD8) were identified in C. farreri, with gene expansion being revealed in Cu/Zn-SODs. In scallops exposed to two different PSTs-producing dinoflagellates, Alexandrium minutum and A. catenella, expression regulation of SOD genes was analyzed in the top ranked toxin-rich organs, the hepatopancreas and the kidney. In hepatopancreas, which mainly accumulates the incoming PSTs, all of the six Cu/Zn-SODs showed significant alterations after A. minutum exposure, with SOD1, 2, 3, 5, and 6 being up-regulated, and SOD4 being down-regulated, while no significant change was detected in Mn-SODs. After A. catenella exposure, up-regulation was observed in SOD2, 4, 6, and 8, and SOD7 was down-regulated. In the kidney, where PSTs transformation occurs, SOD4, 5, 6, and 8 were up-regulated, and SOD7 was down-regulated in response to A. minutum feeding. After A. catenella exposure, all the Cu/Zn-SODs except SOD1 were up-regulated, and SOD7 was down-regulated in kidney. Overall, in scallops after ingesting different toxic algae, SOD up-regulation mainly occurred in the expanded Cu/Zn-SOD group, and SOD6 was the only member being up-regulated in both toxic organs, which also showed the highest fold change among all the SODs, implying the importance of SOD6 in protecting scallops from the stress of PSTs. Our results suggest the diverse function of scallop SODs in response to the PST-producing algae challenge, and the expansion of Cu/Zn-SODs might be implicated in the adaptive evolution of scallops or bivalves with respect to antioxidant defense against the ingested toxic algae.