Next Article in Journal
The Extreme Ends of the Treatment Response Spectrum to Botulinum Toxin in Cervical Dystonia
Next Article in Special Issue
Transcriptional Response in the Digestive Gland of the King Scallop (Pecten maximus) After the Injection of Domoic Acid
Previous Article in Journal
The Anemonia sulcata Toxin BDS-I Protects Astrocytes Exposed to Aβ1–42 Oligomers by Restoring [Ca2+]i Transients and ER Ca2+ Signaling
Previous Article in Special Issue
Differentiating Two Closely Related Alexandrium Species Using Comparative Quantitative Proteomics
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Toxins
Volume 13
Issue 1
10.3390/toxins13010021
toxins-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Tissue-Biased and Species-Specific Regulation of Glutathione Peroxidase (GPx) Genes in Scallops Exposed to Toxic Dinoflagellates

by
Sein Moh Moh Hlaing
1,†,
Jiarun Lou
1,†,
Jie Cheng
1,2,
Xiaogang Xun
1,2,
Moli Li
1,
Wei Lu
1,*,
Xiaoli Hu
1,2,* and
Zhenmin Bao
1,2
1
Key Laboratory of Marine Genetics and Breeding, College of Marine Life Sciences, Ocean University of China, Ministry of Education, 5 Yushan Road, Qingdao 266003, China
2
Laboratory for Marine Fisheries Science and Food Production Processes, Pilot National Laboratory for Marine Science and Technology (Qingdao), 1 Wenhai Road, Qingdao 266237, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally.
Toxins 2021, 13(1), 21; https://doi.org/10.3390/toxins13010021
Submission received: 15 November 2020 / Revised: 20 December 2020 / Accepted: 28 December 2020 / Published: 31 December 2020
(This article belongs to the Special Issue Omic Technologies Applied to the Study of Marine Shellfish Toxins)
Browse Figures
Versions Notes

Abstract

:
Marine bivalves could accumulate paralytic shellfish toxins (PSTs) produced by toxic microalgae, which might induce oxidative stress. Glutathione peroxidases (GPxs) are key enzymes functioning in the antioxidant defense, whereas our understanding of their roles in PST challenge in bivalves is limited. Herein, through genome-wide screening, we identified nine (CfGPx) and eight (PyGPx) GPx genes in Zhikong scallop (Chlamys farreri) and Yesso scallop (Patinopecten yessoensis), respectively, and revealed the expansion of GPx3 sub-family in both species. RNA-Seq analysis revealed high expression of scallop GPx3s after D stage larva during early development, and in adult hepatopancreas. However, in scallops exposed to PST-producing dinoflagellates, no GPx was significantly induced in the hepatopancreas. In scallop kidneys where PSTs were transformed to higher toxic analogs, most CfGPxs were up-regulated, with CfGPx3s being acutely and chronically induced by Alexandrium minutum and A. catenella exposure, respectively, but only one PyGPx from GPx3 subfamily was up-regulated by A. catenella exposure. Our results suggest the function of scallop GPxs in protecting kidneys against the oxidative stresses by PST accumulation or transformation. The tissue-, species-, and toxin-dependent expression pattern of scallop GPxs also implied their functional diversity in response to toxin exposure.
Keywords:
Chlamys farreri; Patinopecten yessoensis; glutathione peroxidase (GPx); paralytic shellfish toxin (PST); antioxidant defense; dinoflagellates Alexandrium
Key Contribution: We systematically identified the glutathione peroxidase (GPxs) genes in two scallop species, revealing the expansion of GPx3 sub-family. We further found that scallop GPx3s were the dominantly expressed GPx genes, but their regulation after exposure to PST-producing Alexandrium was dependent on scallop organs, scallop species, and PSTs from the ingested Alexandrium.

1. Introduction

Reactive oxygen species (ROS) are produced in all living organisms as by-products of aerobic metabolism [1]. The excessive production of ROS can cause lipid peroxidation, protein oxidation, DNA damage, membrane disruption, cellular damage, immune dysfunction, as well as metabolic malfunction [2,3]. Low levels of ROS production are required to maintain the physiological functions, including proliferation, host defense, signal transduction, and gene expression that are involved in the antioxidant defense [4,5]. Aerobic organisms have developed various non-enzymatic and enzymatic antioxidant defense systems to maintain the appropriate level of ROS and protect against the oxidative stress which occurs imbalance between the production of ROS and the antioxidant capacity of the cells [6,7]. Among these enzymatic defense systems, the superoxide dismutase (SOD), catalase (CAT), glutathione S-transferase (GST), and glutathione peroxidase (GPx) are the important enzyme families involved in the detoxification of ROS [6,7].
GPxs are important members of the enzymatic antioxidant system that protects organisms from oxidative damage and scavenge peroxides generated in the cells [8]. They catalyze the reduction of H2O2 and organic hydro-peroxides to water or the corresponding alcohol using glutathione (GSH) as an electron donor [9]. In this way, they maintain a balance between the level of ROS and antioxidants [10]. GPxs are generally classified into two sub-groups: selenium-dependent glutathione peroxidase (SeGPx) and non-selenium glutathione peroxidase (non-SeGPx) based on the presence of selenocysteine (Sec, U) or cysteine (Cys, C) residues at their active sites [11]. The difference between these two enzyme types is that SeGPx catalyzes the reduction of both organic and inorganic peroxides, while non-SeGPx acts only on organic peroxides [12,13,14]. Until now, four major SeGPx isozymes have been reported in mammals [9]: (a) classical GPx (GPx1), (b) gastrointestinal GPx (GPx2), (c) plasma GPx (GPx3), and (d) phospholipid hydroperoxide GPx (PHGPx or GPx4). Moreover, there are also four non-SeGPx isozymes: epididymal GPx (GPx5), olfactory epithelium GPx (GPx6), and two endoplasmic reticulum phospholipid hydroperoxide GPxs (GPx7 and GPx8) [8,11,15,16,17].
GPxs have also been reported in mollusks, which live in diverse aquatic environments and face dynamic stresses. For example, the induction of GPxs was observed in response to various bacterial infections, like in Yesso scallop Patinopecten yessoensis [18], Zhikong scallop Chlamys farreri [19], Manila clam Venerupis philippinarum [20], disk abalone Haliotis discus discus [21], and clam Meretrix meretrix [22], all suggesting the significant role of GPxs in innate immune system of mollusks. Moreover, the up-regulation of GPxs under different chemical contaminant exposure have also been reported in Manila clam V. philippinarum [23], freshwater bivalve Unio tumidus [17], mussel Mytilus galloprovincialis [16], freshwater bivalve Anodonta woodiana [24], thick shell mussel Mytilus coruscus [25], Pacific oyster Crassostrea gigas [26], and abalone Haliotis discus hannai Ino [27], indicating the involvement of mollusk GPxs in the oxidative stress response caused by environmental pollutants.
Bivalves such as mussels, clams, oysters, as well as scallops, are all filter-feeders that can accumulate paralytic shellfish toxins (PSTs) from microalgae, especially the dinoflagellates of the genus Alexandrium [28]. PST exposure is known to induce oxidative stress in bivalve species through over-production of ROS [29,30,31]. The activity of antioxidant enzymes, such as SOD, CAT, GST, and GPx, could be enhanced and contributed to the removal of excessive ROS. For example, scallop SOD [32] and GST [33] genes were both induced after Alexandrium exposure; The CAT and GPx enzymes could be oxidative stress biomarkers [34] and were involved in detoxification in responding to PST accumulation in bivalve tissues [35], and the high enzymatic activity also mirrored the gene expression results [36]. Moreover, different organisms may present varied levels of tolerance to PSTs, and the activities of antioxidant enzymes could respond specifically with cell types [37], which indicate that the antioxidant system may be responsible for varied sensitivity to PST exposure. Like many other enzymes, the GPx activities are contributed by the expression of several GPx genes. Revealing the expression changes of all the GPx members during PSTs exposure could be helpful to our understanding of the antioxidant response variation mediated by GPx in bivalves.
Scallops represent bivalve species in which PSTs could be accumulated at high concentration and with long retaining time [38,39], but systematic analysis of scallop GPx genes has not been reported. In this study, we performed whole-genome identification of GPx genes in two scallops, Zhikong scallop, C. farreri and Yesso scallop, P. yessoensis. Their expression profiles in normal scallops and those exposed to PST-producing dinoflagellates were analyzed via RNA-Seq data analysis. We revealed the expansion of scallop GPx genes, and their tissue- and species-specific expression pattern during toxic algae challenge. These data expanded our understanding of GPx functions in the detoxification and antioxidant defense of shellfish.

2. Results and Discussion

2.1. Genome-Wide Identification of GPx Genes in C. farreri and P. yessoensis

A total of nine and eight GPx genes were identified in C. farreri (CfGPx) and P. yessoensis (PyGPx) genomes, respectively. General information regarding their genome position, intron number, encoding protein length, GPx domain position, isoelectric point (pI) value, and molecular weight was summarized in Table 1, and the coding nucleotide and deduced amino acid sequences of CfGPxs and PyGPxs were shown in Figure S1. In C. farreri genome, eight SeGPx genes were identified, including two GPx1s (CfGPx1-1 and CfGPx1-2), five GPx3s (CfGPx3-1, CfGPx3-2, CfGPx3-3a, CfGPx3-3b, and CfGPx3-4), and one GPx4 (CfGPx4), whereas only one non-SeGPx gene, CfGPx7, was found. Similarly, in P. yessoensis genome, one PyGPx1, five GPx3s (PyGPx3-1, PyGPx3-2, PyGPx3-3, PyGPx3-4, and PyGPx3-5), one PyGPx4, and one PyGPx7 were obtained. Moreover, the coding sequence of CfGPxs ranged from 366 to 843 bp in length and encoded proteins from 121 to 280 amino acids (aa), whereas the PyGPxs varied from 327 to 834 bp in length and encoded proteins between 108 and 277 aa. The CfGPxs had an average theoretical pI from 4.35 to 8.90 and an average molecular weight from 14.22 kDa to 24.68 kDa, while the theoretical pI of PyGPxs ranged from 4.43 to 9.41 and the average molecular masses varied from 12.73 kDa to 31.41 kDa.
Notably, the GPx3 sub-family has more duplicated copies in scallop genomes than in mammal, teleost, and amphibian, indicating the GPx3 expansion in scallop genome, but members of four GPx sub-families, including GPx2, GPx5, GPx6, and GPx8, were all absent from the scallop genomes (Table 2). Moreover, it was found that the amino acid sequences of GPx genes in C. farreri and P. yessoensis shared great sequence identity of more than 40% with vertebrates and other mollusks (Table S1), suggesting the conservation of the GPx family between vertebrates and invertebrates.

2.2. Conserved Structures of GPx Genes in C. farreri and P. yessoensis

The GPx sequence alignment among selected species revealed the presence of conserved amino acid residues, Cys (C), and the active site residues, Gln (Q) and Trp (W) (Figure 1), which were responsible for the fixation of selenium (Sec) [18,19,40]. In addition, the GPx family members were characterized by the presence of a conserved motif (GPx signature 1), (G[K/R]X[L/V][I/L]I[V/E/T]NVA[S/T/A][E/Q/L/Y][C/U]G[L/T]T) with a conserved Cys or Sec residue [41]. Two other conserved GPx domains: GPx signature motif 2, LAFPCNQF and active site motif WNF(S/T)KF, are critical sites for the catalytic activity of this enzyme, which are totally dependent on GSH for its regeneration [41]. It was clear from the sequence alignment that CfGPxs and PyGPxs shared the conserved features found in well-characterized GPxs from other vertebrates. For example, the typical GPx signature motif 2, L(G/A) (F/V)PC(N/D)QF and the active site motif, WNF(E/T/W)K(F/Y) were also highly conserved in scallop (Figure 1 and Figure S1), although the GPx signature motif 2 was missing in some scallop GPx3 sub-family members. Moreover, the GPx signature motif 1 (GQVSLVVNVASECGFT) was only predicted in the non-SeGPx, GPx7, of both scallops (Figure 1 and Figure S1). Qu et al. have also identified the signature sequence motif, LGFPCNQF, an extra active site motif, WNFEKF, and the GPx-1 active site motif, GKVILVENVASLUGTT, in the N-terminal region of thick shell mussel, M. coruscus [25].
In addition, two arginine (Arg, R) residues which direct the GSH donor substrate towards the catalytic center [42] were also well conserved in both scallop GPxs (Figure 1 and Figure S1). Two arginine residues (R, R) involved in binding GSH, and two amino acid residues glutamine and tryptophan (Q and W) responded for fixation of Sec in GPx were also observed in SeGPx of some mollusks, M. coruscus [25], P. yessoensis [18], C. farreri [19], and H. discus discus [21]. Moreover, the structural stability of GPx was reported to depend on three-loop structures where the first loop is Asn (N) to Tyr (Y), the second is Leu (L) to Gln (Q), and the third is Trp (W) to Phe (F) [43]. These three loop structures were also conserved in both scallop GPx sequences and were essential for maintaining the biochemical properties and enzymatic activities of GPxs [17,18]. In addition, three other conserved domains and four highly potential catalytic sites (Cys, Gln, Trp, and Asn) were also present in most amino acid sequences of both scallop GPxs (Figure 1 and Figure S1). These conserved amino acids were likely to stabilize GPx structure and function in evolutionary terms.

2.3. Phylogenetic Relationship of GPxs between Bivalves and Other Organisms

Phylogenetic analysis of GPx amino acid sequences from 19 selected species was conducted (Figure 2). As the result, the scallop GPx members can be classified into four major clades. One clade contained vertebrate GPx1 and GPx2 members, as well as mollusk GPx1s clustered with vertebrate GPx2. CfGPx1s and PyGPx1 formed a sub-cluster with GPx1s from other bivalve species: pearl oyster (P. fucata), Pacific oyster (C. gigas), and Mediterranean mussel (M. galloprovincialis). In another clade, the GPx3s of both scallops initially formed a marine bivalve GPx3 group and clustered with other mollusk GPx3s, such as Pacific abalone (H. discus hannai), zebra mussel (D. polymorpha), freshwater mussel (U. tumidus), Pacific oyster (C. gigas), and Manila clams (R. philippinarum). All these invertebrate GPx3s formed a sub-clade and clustered with vertebrate GPx3, GPx5, and GPx6 members, which indicated that vertebrate GPx3, 5, and 6 shared the same ancestor with mollusk GPx3 sub-family. Moreover, the GPx4 of both scallops positioned as an independent cluster with a direct linking to the vertebrate GPx4 clade. Besides, the scallop GPx7s were in a clade composed of GPx7s from Pacific oyster, C. gigas and then clustered with vertebrate GPx7s and GPx8s as the sisters of the well-supported GPx7/8 sub-group, which indicated that the mollusk GPx8 was probably lost in their common ancestor during evolution. The phylogenetic analysis provided evidence that CfGPxs and PyGPxs were derived from a common ancestor with other GPx family proteins as reported by De Zoysa et al. [21]. Zhang et al. also stated that the GPx phylogenetic relationships were in good agreement with traditional taxonomy, suggesting that this protein family might have a primarily similar functional role [20].

2.4. Spatio-Temporal Expression of Scallop GPxs During Development and in Adult Tissues

The spatio-temporal expression profiles of scallop GPxs were analyzed using transcriptome data [38,44]. In scallop embryos and larvae, CfGPx1-1, CfGPx3-1, CfGPx3-3b, and CfGPx4 in C. farreri, and PyGPx1, PyGPx3-1, PyGPx3-5, and PyGPx4 in P. yessoensis exhibited relatively higher expression level than other GPxs, suggesting the importance of these GPxs in scallop development (Figure 3a). During the entire developmental process, high expression of GPx genes (CfGPx3-1, CfGPx3-3b, PyGPx1, PyGPx3-1) was mainly present at the larval stages, especially after the D stage veliger, and in the following umbo larvae and juvenile stages (Figure 3a), implying the involvement of these GPxs in antioxidation or detoxification in scallop metamorphosis and post-larval development, which may be due to the elevation of oxygen consumption to meet high demands of energy reserve utilization during organ initiation and structural re-modelling [32]. Moreover, both CfGPx4 and PyGPx4 were expressed from the very beginning of fertilization and multicellular stages to the larvae and juvenile stages, indicating their maternal origin to play protective roles and to maintain a redox balance during development.
In adult individuals, the identified scallop GPx genes were expressed in all the examined tissues/organs, except that CfGPx1-2 transcript was not detected in any tissue (Figure 3b). Notably, higher expression of GPxs was mainly detected in the hepatopancreas of both scallops. These highly expressed GPxs included CfGPx3-1, CfGPx3-2, CfGPx3-3a, CfGPx3-4, and CfGPx4 in C. farreri, and PyGPx1, PyGPx3-1, PyGPx3-2, PyGPx3-3, PyGPx3-4, and PyGPx4 in P. yessoensis, most of which were from the GPx3 sub-family. The dominant GPx3s expression in scallop hepatopancreas was consistent with its role as the main metabolic and defense organs against oxidative stress caused by excessive ROS [27], which suggested that scallop GPx3s might play crucial roles in molluscan metabolism and antioxidative defense system. The high expression of SeGPx in hepatopancreas was also reported in abalone H. discus discus [21] and H. discus hannai Ino [27], freshwater bivalve U. tumidus [17,45] and D. polymorpha [40], and clam M. meretrix [22] and V. philippinarum [14,20]. In addition, GPx4s of both scallops, as well as PyGPx1, were transcribed in almost all the examined tissues. Similar expression patterns were also reported in other invertebrates. For example, the GPxs of M. galloprovincialis were mainly expressed in the gills, hemocytes, and digestive gland, moderately expressed in adductor muscle and gonad [16], and the GPx of L. vannamei was widely distributed in the hemocytes, hepatopancreas, gills, muscles, and intestine tissues [46].

2.5. Expression Regulation of Scallop GPxs in Response to PST-Producing Dinoflagellates

As filter-feeders, scallops have a high ability to accumulate PSTs [37]. Previous studies have suggested that the scallop hepatopancreas and kidneys were the two major organs involved in PST bioaccumulation and transformation, respectively [38,47,48]. Moreover, dinoflagellates of the genus Alexandrium are the major PST producers containing varied toxic analogs among species and strains. Therefore, to understand the scallop defense mechanism against toxic algae, RNA-Seq analysis was performed to determine the expression profiles of GPx genes in scallop kidneys and hepatopancreas in response to the ingestion of PST-producing algae, A. minutum (AM-1 strain) and A. catenella (ACDH strain) (Tables S2 and S3).
In C. farreri hepatopancreas, the major organ for PST accumulation from algae [38,48], most CfGPxs’ expression was not significantly altered after toxic algae exposure, with only CfGPx7 down-regulated at the 1st day after A. catenella exposure (log2FC: −2.06) (Figure 4a). This result implied that CfGPxs mediated antioxidant defense in hepatopancreas was not sensitive to PST challenge, and the highly expressed CfGPx3s might provide stable antioxidant capacity for redox homeostasis during toxic algae ingestion.
In C. farreri kidneys, where the ingested PSTs could be transformed to higher toxic analogs [38], a prominent up-regulation of CfGPx1-1, CfGPx3-1, CfGPx3-2, CfGPx3-3a, and CfGPx3-4 was detected after A. catenella exposure (log2FC: 2.18–5.61) (Figure 4a). Similarly, A. minutum exposure significantly induced the up-regulation of CfGPx1-1, CfGPx1-2, CfGPx3-1, CfGPx3-2, CfGPx3-3a, and CfGPx3-4 in this organ (log2FC: 2.19–5.88) (Figure 4a). For both toxic algae challenge, the majority of up-regulated CfGPxs in kidney were from the GPx1 and GPx3 sub-families, suggesting the importance of GPx1 and GPx3 members in PST detoxification in scallop kidneys, especially during the transformation of PSTs to more potent analogs [38]. The duplication of CfGPx3s and CfGPx1s might be related to the detoxification of C. farreri with respect to PST metabolism in the toxin transformation center kidney.
However, the regulation of CfGPxs in scallop kidneys exhibited different patterns between the two toxic algae exposure (Figure 4a). For example, the induction of CfGPx3s was mostly chronic (5th–10th days, log2FC: 2.18–5.61) with A. catenella challenge, whereas acute (1st–3rd days, log2FC: 2.19–5.45) with A. minutum exposure. Moreover, the induction of CfGPx1-2 and the down-regulation of CfGP3xs were only observed in the kidney of scallops exposed to A. minutum (Figure 4a and Table S3). Since the two micro-algae produce different PST analogs, with A. minutum mainly synthesizing gonyautoxins (GTXs, mainly GTX1-4) and A. catenella synthesizing N-sulfocarbamoyl toxins (C1 and C2) [38,48], it therefore suggested that the difference in time dependent pattern of GPx3s response between A. catenella and A. minutum challenge in C. farreri kidneys was probably due to the different toxin species accumulated (Figure S2). This diverse expression profiles of CfGPxs indicated that, with the challenges of different PST analogs, the CfGPxs functioning in PST detoxification and antioxidation during PST transformation were similar in scallop kidneys, but the gene stimulation or inhibition procedure/speed might vary.
To further compare the regulation of GPxs between scallop species, we examined PyGPxs transcripts in both hepatopancreas and kidneys of P. yessoensis ingesting A. catenella (Figure 4b). There were less regulated GPxs than those identified in C. farreri, with only PyGPx3-4 being up-regulated in P. yessoensis kidney (10th day, log2FC: 2.88) and three PyGPx3s down-regulated (3rd and 5th days, log2FC: −2.42–−2.95) (Figure 4b). In P. yessoensis hepatopancreas, only PyGPx3-5 was down-regulated at the 3rd day (log2FC: −2.60). The inhibition of PyGPxs might be resulted from the decline of GSH contents under cellular redox imbalance, to maintain the stable GSH/GSSG ratio which is crucial in preventing oxidative damage [49]. Moreover, the species-specific and tissue-biased expression of genes from the expanded GPx3 sub-family may also indicate the different sensitivity of scallops to varied PSTs (Figure S2) and the possible association of GPx3 expansion in the adaptation of scallops to the environment with different toxic algae.
The antioxidant metabolism is crucial for preventing the cellular oxidative damage in bivalves exposed to toxic dinoflagellates, such as Alexandrium sp. [30]. Stress response is usually modulated at the molecular level and shows significant changes in the transcripts expressed in the target tissues [50]. Many studies have described the impact of PST-producers on gene expression of diverse bivalve organs, majorly in C. gigas. For example, both repression and activation of genes in C. gigas gills by A. minutum induction was involved in key metabolic pathways such as oxidative stress, mitochondrial metabolism, endogenous clock, immunity, and detoxification processes, including SODs, GSTs, and GPxs [30,51]. The expression of defense-related genes modulated in the hemocytes of oyster after exposure to A. minutum could demonstrate the impact of direct contacting of some Alexandrium species on different cell types due to the production of extracellular compounds [52]. Expression of genes involved in antioxidant stress, detoxification, immune, and stress defense of C. gigas juveniles increased in a short exposure time by the toxic dinoflagellate, G. catenatum [53,54] and P. lima exposure [55]. Therefore, the high expression level of CfGP3xs maintained in scallop hepatopancreas, as well as the significantly regulated expression of CfGPxs in kidneys reasonably supported the fact that these two organs were the main centers involved in detoxification and antioxidant defense, but with different roles. These findings also indicated that the expression of GPx genes varied among scallop species and tissues, dinoflagellate strains, as well as the exposure time. The species- and organ-specific regulation of GPxs imply the diverse function of scallop GPxs in the defense against the harmful effects of PST accumulation and metabolism, which might contribute to the scallop adaptation to diverse marine environments.

3. Conclusions

In conclusion, a total of nine CfGPx and eight PyGPx genes were identified in scallops C. farreri and P. yessoensis, respectively, with the expansion of GPx3 sub-family in both species. GPx3s showed higher expression than other GPx genes, and hepatopancreas was the major organ for scallop GPx expression. After ingesting PST-producing algae, no GPx was induced in hepatopancreas of either scallop, but most GPxs were up-regulated in C. farreri kidneys, with acute and chronic regulation by A. minutum and A. catenella exposure, respectively. Meanwhile, only one GPx3 was induced in kidneys of P. yessoensis exposed to A. catenella. These results indicated that the response of scallop GPxs to PST-producing dinoflagellates was dependent on scallop species and tissues, as well as dinoflagellates strains. Our data revealed the involvement of GPxs in scallops against the toxic stress and highlighted their functional importance, which provided a better understanding of the detoxification systems in scallops and their regulation mechanism.

4. Materials and Methods

4.1. Screening GPx Genes from the C. farreri and P. yessoensis Genomes

The GPx genes were identified using the C. farreri and P. yessoensis transcriptome and genome data [38,44]. The coding sequences were translated using the ORF finder (https://www.ncbi.nlm.nih.gov/orffinder/), and the predicted ORFs were validated using the NCBI non-redundant (Nr) protein sequence database. The protein sequences of GPxs were further verified with BLASTP in NCBI (http://blast.ncbi.nlm.nih.gov/) and UniProt (http://www.uniprot.org/) with the setting of e-value to 1E-05. All possible GPx proteins from the BLASTP results were then tested using SMART (http://smart.embl-heidelberg.de/). Protein sequences longer than 100 amino acids were used for further analysis. The translated sequences were submitted to the SMART program for identification of the signal peptide and other conserved GPx structures. Compute pI/MW tool (https://www.expasy.org/) was used to predict the pI value and molecular weight (kDa).

4.2. Multiple Sequence Alignment and Phylogenetic Analysis of the GPx Gene Family

Multiple sequence alignment of GPx sequences was performed using Clustal W [56] and then was edited by Genedoc software (version 2.7.0) [57]. The amino acid similarity matrix was made using the Clustal Omega software (http://www.ebi.ac.uk/Tools/msa/clustalo/) with GPx protein sequences of both scallops and other selected species. The GPx amino acid sequences of selected species included Homo sapiens, Mus musculus, Bos taurus, Rattus norvegicus, Xenopus tropicalis, Anolis carolinensis, Gallus gallus, Danio rerio, and mollusk species such as Pinctada fucata, U. tumidus, H. discus hannai, Dreissena polymorpha, Ruditapes philippinarum, C. gigas, and M. galloprovincialis, all obtained from Ensembl (http://www.ensembl.org/) and NCBI (https://www.ncbi.nlm.nih.gov/) (Table S4). The phylogenetic tree was constructed using the GPx amino acid sequences of C. farreri and P. yessoensis, as well as other selected species using MEGA 7.0 with the maximum-likelihood (ML) method [58] and the tested best model of LG + G. The bootstrap analysis with 1000 replicates was applied to test the relative support for the branches produced by ML method.

4.3. Expression Analysis of GPxs During Scallop Development and in Adult Organs/Tissues

C. farreri and P. yessoensis were collected from Xunshan Group Co., Ltd. (Rongcheng, Shandong Province, China) and Zhangzidao Group Co., Ltd. (Dalian, Liaoning Province, China), respectively. The RNA-Seq data for scallop embryo/larva including zygote, 2–8 cells, blastula, gastrula, trochophore, D stage larvae, early umbo larvae, mid umbo larvae, late umbo larvae, creeping larvae, and juvenile scallops, and the adult organs/tissues including mantle, gill, eye, kidney, hepatopancreas, hemocyte, female gonad, male gonad, smooth muscle, and adductor muscle, were analyzed. The expression profiles of CfGPxs and PyGPxs during development and in adult organs/tissues were normalized and represented in the form of RPKM (Reads Per Kilobase of exon model per Million mapped reads) using the data from our previous studies [38,44]. To visualize the expression patterns of GPx genes on both scallops, the heatmap was generated via heatmap package under the R software environment with the log10(RPKM + 1) value.

4.4. Expression Analysis of GPx Genes in Scallops Exposed to Toxic Dinoflagellates

The scallops C. farreri were exposed to PST-producing dinoflagellates, A. minutum (AM-1 strain) and A. catenella (ACDH strain) [38,48]. Briefly, 2-year-old adult scallops were collected and acclimated in filtered and aerated seawater at 12–13 °C for three weeks for depuration by feeding the non-toxic algae, Isochrysis galbana (7.5 × 105 cells/mL) as the control group [59]. 18 scallops were separated into independent tanks with aeration, and 6 groups with 3 individuals for each were taken at different test time points (0, 1, 3, 5, 10, and 15 days) during the challenging experiments. The two Alexandrium dinoflagellate strains were cultured independently in F/2 medium using artificial seawater, and were harvested at the late exponential growth phase through centrifugation of 2500 g/10 min [54,60]. Then, each scallop was fed once a day with 3 L (2.5 × 103 cells/mL) A. minutum or A. catenella. For comparison, 2-year-old P. yessoensis were also exposed to A. catenella (ACDH strain) using the same procedure described above. Scallops were sampled at 0 day (control), 1 and 3 (acute response) days, 5, 10, and 15 (sub-chronic reaction) days with three individuals at each time point after AM-1 and ACDH exposure. The kidney and hepatopancreas of these scallops were dissected, washed with sterile seawater, and immediately frozen at −80 °C for subsequent RNA extraction. PST concentrations in the tissues of sampled scallops were determined by high-performance liquid chromatography with tandem mass spectrometry (HPLC-MS/MS) method [38].
Total RNA was isolated from scallop kidneys and hepatopancreas, and the RNA quantity and quality were checked by Nanovue Plus and gel electrophoresis before the following reverse transcription and sequencing analysis. RNA-Seq data obtained from these tissues were used to examine the expression profiles of GPx genes in response to different toxic dinoflagellates. Briefly, the RNA-Seq libraries of these tissues from C. farreri and P. yessoensis fed with AM-1 and ACDH were constructed independently using the NEB Next mRNA Library Prep Kit following the manufacturer’s instructions and were subjected to PE125 sequencing on the Illumina HiSeq 2000 platform. The RNA-Seq reads were mapped to the C. farreri and P. yessoensis genomes, respectively, using Tophat (ver 2.0.9), and the expression of all GPx genes were normalized and represented in the form of RPKM (Table S2). The fold change (FC) between the toxin-exposed and control groups was calculated and the significant differences were determined using the edgeR package with cutoff of |log2FC| > 2 and p value < 0.05, and the very significant difference with cutoff of |log2FC| > 2 and corrected FDR value < 0.05. The heat maps were drawn using pheatmap package with the log2FC values (Table S3).

Supplementary Materials

The following are available online at https://www.mdpi.com/2072-6651/13/1/21/s1, Table S1. The homology percentages of CfGPx and PyGPx amino acid sequences identity compared with other animals, Table S2. RPKM values of GPxs from RNA-Seq analysis after toxic A. catenella and A. miutum exposure in hepatopancreas and kidneys at different time points, Table S3. Different expression of GPxs in hepatopancreas and kidneys after A. catenella and A. minutum exposure analyzed using edgeR package, Table S4. Information of GPxs used in this study from vertebrates and mollusk species, Figure S1. Nucleotide and deduced amino acid sequences of (1) CfGPxs and (2) PyGPxs, Figure S2. The temporal abundance of PSTs in scallop kidneys and hepatopancreas after exposure to the toxic alga A. minutum and A. catenella.

Author Contributions

X.H., Z.B. and W.L. designed the study, S.M.M.H. and J.L. carried out the experiments. S.M.M.H., J.L., J.C., M.L. and X.X. analyzed the data. X.H., J.C., S.M.M.H. and W.L. contributed to the writing of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the National Key R&D Program of China (2019YFC1605704), the National Natural Science Foundation of China (31630081 and 31802292), Taishan Industry Leading Talent Project (LJNY201816), and National Infrastructure of Fishery Germplasm Resources (2020DKA30470).

Institutional Review Board Statement

This study was conducted in accordance with the Institutional Animal Care and Use Committee of Ocean University of China, and it does not contain any studies with human participants.

Data Availability Statement

The RNA-Seq data for exposure experiments in this study were deposited in NCBI Sequence Read Archive, with the accession numbers of SRX2445405-SRX2445440, and the RPKM values and differential expression results of GPxs were in supplementary Tables S2 and S3.

Acknowledgments

The authors wish to thank Rencheng Yu, Fanzhou Kong and Yang Liu (Institute of Oceanology, Chinese Academy of Sciences) for providing algae A. catenella, and their kind assistance in scallop challenge.

Conflicts of Interest

The authors declare that they have no conflict of interest.

References

  1. Lesser, M.P. Oxidative stress in marine environments: Biochemistry and Physiological Ecology. Annu. Rev. Physiol. 2006, 68, 253–278. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Nakamura, H.; Nakamura, K.; Yodoi, J. Redox regulation of cellular activation. Annu. Rev. Immunol. 1997, 15, 351–369. [Google Scholar] [CrossRef] [PubMed]
  3. Malhotra, J.D.; Kaufman, R.J. Endoplasmic Reticulum Stress and Oxidative Stress: A Vicious Cycle or a Double-Edged Sword? Antioxid. Redox Signal. 2007, 9, 2277–2294. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Dröge, W. Free Radicals in the Physiological Control of Cell Function. Physiol. Rev. 2002, 82, 47–95. [Google Scholar] [CrossRef] [PubMed]
  5. Nita, M.; Grzybowski, A. The Role of the Reactive Oxygen Species and Oxidative Stress in the Pathomechanism of the Age-Related Ocular Diseases and Other Pathologies of the Anterior and Posterior Eye Segments in Adults. Oxidative Med. Cell. Longev. 2016, 2016, 1–23. [Google Scholar] [CrossRef] [Green Version]
  6. Birben, E.; Sahiner, U.M.; Sackesen, C.; Erzurum, S.; Kalayci, O. Oxidative Stress and Antioxidant Defense. World Allergy Organ. J. 2012, 5, 9–19. [Google Scholar] [CrossRef] [Green Version]
  7. Migdal, C.; Serres, M. Reactive oxygen species and oxidative stress. Med. Sci. 2011, 27, 405–412. [Google Scholar]
  8. Herbette, S.; Roeckel-Drevet, P.; Drevet, J.R. Seleno-independent glutathione peroxidases: More than simple antioxidant scav-engers. FEBS J. 2007, 274, 2163–2180. [Google Scholar] [CrossRef]
  9. Margis, R.; Dunand, C.; Teixeira, F.K.; Margis-Pinheiro, M. Glutathione peroxidase family–an evolutionary overview. FEBS J. 2008, 275, 3959–3970. [Google Scholar] [CrossRef]
  10. Drevet, J.R. The antioxidant glutathione peroxidase family and spermatozoa: A complex story. Mol. Cell. Endocrinol. 2006, 250, 70–79. [Google Scholar] [CrossRef]
  11. Arthur, J.R. The glutathione peroxidases. Cell. Mol. Life Sci. 2001, 57, 1825–1835. [Google Scholar] [CrossRef] [PubMed]
  12. Almar, M.; Otero, L.; Santos, C.; Gallego, J.G. Liver glutathione content and glutathione-dependent enzymes of two species of freshwater fish as bioindicators of chemical pollution. J. Environ. Sci. Heal. Part. B 1998, 33, 769–783. [Google Scholar] [CrossRef] [PubMed]
  13. Ren, Q.; Sun, R.-R.; Zhao, X.-F.; Wang, J.-X. A selenium-dependent glutathione peroxidase (Se-GPx) and two glutathione S-transferases (GSTs) from Chinese shrimp (Fenneropenaeus chinensis). Comp. Biochem. Physiol. Part. C Toxicol. Pharmacol. 2009, 149, 613–623. [Google Scholar] [CrossRef] [PubMed]
  14. Zhang, L.; Liu, X.; Chen, L.; You, L.; Pei, D.; Cong, M.; Zhao, J.; Li, C.; Liu, N.; Yu, J.; et al. Transcriptional regulation of selenium-dependent glutathione peroxidase from Venerupis philippinarum in response to pathogen and contaminants challenge. Fish. Shellfish. Immunol. 2011, 31, 831–837. [Google Scholar] [CrossRef]
  15. Nguyen, V.D.; Saaranen, M.J.; Karala, A.-R.; Lappi, A.-K.; Wang, L.; Raykhel, I.B.; Alanen, H.I.; Salo, K.E.; Wang, C.-C.; Ruddock, L.W. Two Endoplasmic Reticulum PDI Peroxidases Increase the Efficiency of the Use of Peroxide during Disulfide Bond Formation. J. Mol. Biol. 2011, 406, 503–515. [Google Scholar] [CrossRef]
  16. Wang, Q.; Ning, X.; Zhang, Q.; Liu, F.; Wu, H.; Zhang, Y.; Zhao, J. Molecular characterization of two glutathione peroxidase genes in Mytilus galloprovincialis and their transcriptional responses to sub-chronic arsenate and cadmium exposure. Invertebr. Surviv. J. 2014, 11, 149–162. [Google Scholar]
  17. Doyen, P.; Vasseur, P.; Rodius, F. Identification, sequencing and expression of selenium-dependent glutathione peroxidase transcript in the freshwater bivalve Unio tumidus exposed to Aroclor 1254. Comp. Biochem. Physiol. Part. C Toxicol. Pharmacol. 2006, 144, 122–129. [Google Scholar] [CrossRef]
  18. Shan, Z.; Li, H.; Bao, X.; He, C.; Yu, H.; Liu, W.; Hou, L.; Wang, J.; Zhu, D.; Sui LZhu, B. A selenium-dependent glutathione peroxidase in the Japanese scallop, Mi-zuhopecten yessoensis: cDNA cloning, promoter sequence analysis and mRNA expression. Comp. Biochem. Physiol. Part B Biochem. Mol. Biol. 2011, 159, 1–9. [Google Scholar] [CrossRef]
  19. Mu, C.; Ni, D.; Zhao, J.; Wang, L.; Song, L.; Li, L.; Zhang, H.; Qiu, L.; Cong, M. cDNA cloning and mRNA expression of a selenium-dependent glutathione peroxidase from Zhikong scallop Chlamys farreri. Comp. Biochem. Physiol. Part. B Biochem. Mol. Biol. 2010, 157, 182–188. [Google Scholar] [CrossRef]
  20. Zhang, L.; Ning, X.; Chen, L.; Li, C.; Liu, F.; Wang, Q.; Wu, H.; Zhao, J. Molecular cloning and expression analysis of a selenium-independent glutathione peroxidase identified from Manila clam Venerupis philippinarum. Aquac. Res. 2011, 43, 1176–1183. [Google Scholar] [CrossRef]
  21. De Zoysa, M.; Pushpamali, W.A.; Oh, C.; Whang, I.; Kim, S.J.; Lee, J. Transcriptional up-regulation of disk abalone selenium dependent glutathione peroxidase by H2O2 oxidative stress and Vibrio alginolyticus bacterial infection. Fish. Shellfish. Immunol. 2008, 25, 446–457. [Google Scholar] [CrossRef] [PubMed]
  22. Wang, C.; Huan, P.; Yue, X.; Yan, M.; Liu, B. Molecular characterization of a glutathione peroxidase gene and its expression in the selected Vibrio-resistant population of the clam Meretrix meretrix. Fish. Shellfish. Immunol. 2011, 30, 1294–1302. [Google Scholar] [CrossRef] [PubMed]
  23. Cong, M.; Zhang, L.; Zhang, L.; Zhao, J.; Wu, H.; Chen, H.; Kong, J. Molecular characterization of a Se-containing glutathione peroxidases gene and its expressions to heavy metals compared with non-Se-containing glutathione peroxidases in Venerupis philippinarum. Agri Gene 2016, 1, 46–52. [Google Scholar] [CrossRef]
  24. Xia, X.; Hua, C.; Xue, S.; Shi, B.; Gui, G.; Zhang, N.; Wang, X.; Guo, L. Response of selenium-dependent glutathione peroxidase in the freshwater bivalve Anodonta woodiana exposed to 2,4-dichlorophenol,2,4,6-trichlorophenol and pentachlorophenol. Fish. Shellfish. Immunol. 2016, 55, 499–509. [Google Scholar] [CrossRef]
  25. Qu, C.; Liu, S.; Tang, Z.; Li, J.; Liao, Z.; Qi, P. Response of a novel selenium-dependent glutathione peroxidase from thick shell mussel Mytilus coruscus exposed to lipopolysaccharide, copper and benzo[α]pyrene. Fish. Shellfish. Immunol. 2019, 89, 595–602. [Google Scholar] [CrossRef]
  26. Jo, P.G.; Choi, Y.K.; Choi, C.Y. Cloning and mRNA expression of antioxidant enzymes in the Pacific oyster, Crassostrea gigas in response to cadmium exposure. Comp. Biochem. Physiol. Part. C Toxicol. Pharmacol. 2008, 147, 460–469. [Google Scholar] [CrossRef]
  27. Wu, C.; Mai, K.; Zhang, W.; Ai, Q.; Xu, W.; Wang, X.; Ma, H.; Liufu, Z. Molecular cloning, characterization and mRNA expression of seleni-um-dependent glutathione peroxidase from abalone Haliotis discus hannai Ino in response to dietary selenium, zinc and iron. Comp. Biochem. Physiol. Part. C Toxicol. Pharmacol. 2010, 152, 121–132. [Google Scholar] [CrossRef]
  28. Deeds, J.R.; Landsberg, J.H.; Etheridge, S.M.; Pitcher, G.C.; Longan, S.W. Non-traditional vectors for paralytic shellfish poisoning. Mar. Drugs 2008, 6, 308–348. [Google Scholar] [CrossRef]
  29. Estrada, N.; de Jesús Romero, M.; Campa-Córdova, A.; Luna, A.; Ascencio, F. Effects of the toxic dinoflagellate, Gymnodinium cat-enatum on hydrolytic and antioxidant enzymes, in tissues of the giant lions-paw scallop Nodipecten subnodosus. Comp. Biochem. Physiol. Part. C Toxicol. Pharmacol. 2007, 146, 502–510. [Google Scholar] [CrossRef]
  30. Fabioux, C.; Sulistiyani, Y.; Haberkorn, H.; Hégaret, H.; Amzil, Z.; Soudant, P. Exposure to toxic Alexandrium minutum activates the detoxifying and antioxidant systems in gills of the oyster Crassostrea gigas. Harmful Algae 2015, 48, 55–62. [Google Scholar] [CrossRef] [Green Version]
  31. Qiu, J.; Ma, F.; Fan, H.; Li, A. Effects of feeding Alexandrium tamarense, a paralytic shellfish toxin producer, on antioxidant enzymes in scallops (Patinopecten yessoensis) and mussels (Mytilus galloprovincialis). Aquaculture 2013, 76–81. [Google Scholar] [CrossRef]
  32. Lian, S.; Zhao, L.; Xun, X.; Lou, J.; Li, M.; Li, X.; Wang, S.; Zhang, L.; Hu, X.; Bao, Z. Genome-Wide Identification and Characterization of SODs in Zhikong Scallop Reveals Gene Expansion and Regulation Divergence after Toxic Dinoflagellate Exposure. Mar. Drugs 2019, 17, 700. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Lou, J.; Cheng, J.; Xun, X.; Li, X.; Li, M.; Zhang, X.; Li, T.; Bao, Z.; Hu, X. Glutathione S-transferase genes in scallops and their diverse expression patterns after exposure to PST-producingdinoflagellates. Mar. Life Sci. Technol. 2020, 2, 252–261. [Google Scholar] [CrossRef]
  34. Freitas, R.; Marques, F.; De Marchi, L.; Vale, C.; Botelho, M. Biochemical performance of mussels, cockles and razor shells con-taminated by paralytic shellfish toxins. Environ. Res. 2020, 188, 109846. [Google Scholar] [CrossRef] [PubMed]
  35. Mat, A.M.; Klopp, C.; Payton, L.; Jeziorski, C.; Chalopin, M.; Amzil, Z.; Tran, D.; Wikfors, G.H.; Hégaret, H.; Soudant PHuvet, A. Oyster transcriptome response to Alexandrium exposure is related to saxitoxin load and characterized by disrupted digestion, energy balance, and calcium and sodium signaling. Aquat. Toxicol. 2018, 199, 127–137. [Google Scholar] [CrossRef] [PubMed]
  36. Prego-Faraldo, M.V.; Vieira, L.R.; Eirin-Lopez, J.M.; Mendez, J.; Guilhermino, L. Transcriptional and biochemical analysis of anti-oxidant enzymes in the mussel Mytilus galloprovincialis during experimental exposures to the toxic dinoflagellate Prorocentrum lima. Mar. Environ. Res. 2017, 129, 304–315. [Google Scholar] [CrossRef]
  37. Zhu, M.Y.; Zou, Y.L.; Wu, R.J.; Hall, S. Accumulation and depuration of paralytic shellfish poisons (PSP) in Chinese scallop Chlamys farreri. Acta Oceanol. Sin. 2003, 25, 75–83. [Google Scholar]
  38. Li, Y.; Sun, X.; Hu, X.; Xun, X.; Zhang, J.; Guo, X.; Jiao, W.; Zhang, L.; Liu, W.; Wang, J.; et al. Scallop genome reveals molecular adaptations to semi-sessile life and neu-rotoxins. Nat. Commun. 2017, 8, 1–11. [Google Scholar] [CrossRef]
  39. Liu, Y.; Kong, F.-Z.; Xun, X.-G.; Dai, L.; Geng, H.-X.; Hu, X.-L.; Yu, R.-C.; Bao, Z.-M.; Zhou, M.-J. Biokinetics and biotransformation of paralytic shellfish toxins in different tissues of Yesso scallops, Patinopecten yessoensis. Chemosphere 2020, 261, 128063. [Google Scholar] [CrossRef]
  40. Doyen, P.; Bigot, A.; Vasseur, P.; Rodius, F. Molecular cloning and expression study of pi-class glutathione S-transferase (pi-GST) and selenium-dependent glutathione peroxidase (Se-GPx) transcripts in the freshwater bivalve Dreissena polymorpha. Comp. Biochem. Physiol. Part. C Toxicol. Pharmacol. 2008, 147, 69–77. [Google Scholar] [CrossRef]
  41. Churin, Y.; Schilling, S.; Börner, T. A gene family encoding glutathione peroxidase homologues in Hordeum vulgare (barley). FEBS Lett. 1999, 459, 33–38. [Google Scholar] [CrossRef] [Green Version]
  42. Aumann, K.D.; Bedorf, N.; Brigelius-Flohé, R.; Schomburg, D.; Flohé, L. Glutathione peroxidase revisited—Simulation of the catalytic cycle by computer-assisted molecular modelling. Biomed. Environ. Sci. 1997, 10, 136–155. [Google Scholar] [PubMed]
  43. Maiorino, F.M.; Brigelius-Flohé, R.; Aumann, K.; Roveri, A.; Schomburg, D.; Flohé, L. [5] Diversity of glutathione peroxidases. Methods Enzymol. 1995, 252, 38–53. [Google Scholar] [CrossRef]
  44. Wang, S.; Zhang, J.; Jiao, W.; Li, J.; Xun, X.; Sun, Y.; Guo, X.; Huan, L.; Dong, B.; Zhang, L.; et al. Scallop genome provides insights into evolution of bilaterian karyotype and development. Nat. Ecol. Evol. 2017, 1, 120. [Google Scholar] [CrossRef]
  45. Cossu, C.; Doyotte, A.; Jacquin, M.; Babut, M.; Exinger, A.; Vasseur, P. Glutathione reductase, selenium-dependent glutathione peroxidase, glutathione levels, and lipid peroxidation in freshwater bivalves, Unio tumidus, as biomarkers of aquatic con-tamination in field studies. Ecotoxicol. Environ. Saf. 1997, 38, 122–131. [Google Scholar] [CrossRef]
  46. Liu, C.; Tseng, M.-C.; Cheng, W. Identification and cloning of the antioxidant enzyme, glutathione peroxidase, of white shrimp, Litopenaeus vannamei, and its expression following Vibrio alginolyticus infection. Fish. Shellfish. Immunol. 2007, 23, 34–45. [Google Scholar] [CrossRef]
  47. Cheng, J.; Xun, X.; Kong, Y.; Wang, S.; Yang, Z.; Li, Y.; Kong, D.; Wang, S.; Zhang, L.; Hu, X.; et al. Hsp70 gene expansions in the scallop Patinopecten yessoensis and their expression regulation after exposure to the toxic dinoflagellate Alexandrium catenella. Fish. Shellfish. Immunol. 2016, 58, 266–273. [Google Scholar] [CrossRef]
  48. Hu, B.; Li, M.; Yu, X.; Xun, X.; Lu, W.; Li, X.; Li, Y.; Lou, J.; Wang, S.; Zhang, L.; et al. Diverse expression regulation of Hsp70 genes in scallops after exposure to toxic Alexandrium dinoflagellates. Chemosphere 2019, 234, 62–69. [Google Scholar] [CrossRef]
  49. Zhou, Z.; Tang, X.; Chen, H.; Wang, Y. Comparative studies of saxitoxin (STX) -induced cytotoxicity in Neuro-2a and RTG-2 cell lines: An explanation with respect to changes in ROS. Chemosphere 2018, 192, 66–74. [Google Scholar] [CrossRef]
  50. Liu, W.; He, C.; Li, W.; Zhou, Z.; Gao, X.; Fu, L. Discovery of host defence genes in the Japanese scallop Mizuhopecten yessoensis Jay by expressed sequence tag analysis of kidney tissue. Aquac. Res. 2010, 41, 1602–1613. [Google Scholar] [CrossRef]
  51. Mat, A.M.; Haberkorn, H.; Bourdineaud, J.-P.; Massabuau, J.-C.; Tran, D. Genetic and genotoxic impacts in the oyster Crassostrea gigas exposed to the harmful alga Alexandrium minutum. Aquat. Toxicol. 2013, 140, 458–465. [Google Scholar] [CrossRef] [PubMed]
  52. Mello, D.F.; Da Silva, P.M.; Barracco, M.A.; Soudant, P.; Hégaret, H. Effects of the dinoflagellate Alexandrium minutum and its toxin (saxitoxin) on the functional activity and gene expression of Crassostrea gigas hemocytes. Harmful Algae 2013, 26, 45–51. [Google Scholar] [CrossRef]
  53. García-Lagunas, N.; Romero-Geraldo, R.D.J.; Hernández-Saavedra, N.Y. Changes in gene expression and histological injuries as a result of exposure of Crassostrea gigas to the toxic dinoflagellate Gymnodinium catenatum. J. Molluscan Stud. 2015, 82, 193–200. [Google Scholar] [CrossRef] [Green Version]
  54. García-Lagunas, N.; Romero-Geraldo, R.; Hernandez-Saavedra, N.Y. Genomics Study of the Exposure Effect of Gymnodinium catenatum, a Paralyzing Toxin Producer, on Crassostrea gigas’ Defense System and Detoxification Genes. PLoS ONE 2013, 8, e72323. [Google Scholar] [CrossRef] [Green Version]
  55. Romero-Geraldo, R.D.J.; Hernández-Saavedra, N.Y. Stress Gene Expression in Crassostrea gigas (Thunberg, 1793) in response to experimental exposure to the toxic dinoflagellate Prorocentrum lima (Ehrenberg) Dodge, 1975. Aquac. Res. 2014, 45, 1512–1522. [Google Scholar] [CrossRef]
  56. Larkin, M.A.; Blackshields, G.; Brown, N.P.; Chenna, R.; Mcgettigan, P.A.; McWilliam, H.; Valentin, F.; Wallace, I.M.; Wilm, A.; Lopez, R.; et al. Clustal W and Clustal X version 2.0. Bioinformatics 2007, 23, 2947–2948. [Google Scholar] [CrossRef] [Green Version]
  57. Nicholas, K.B. GeneDoc: Analysis and visualization of genetic variation. Embnew. News 1997, 4, 14. [Google Scholar]
  58. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [Green Version]
  59. Escobedo-Lozano, A.Y.; Estrada, N.; Ascencio, F.; Contreras, G.; Alonso-Rodriguez, R. Accumulation, biotransformation, histo-pathology and paralysis in the Pacific calico scallop Argopecten ventricosus by the paralyzing toxins of the dinoflagellate Gymnodinium catenatum. Mar. Drugs 2012, 10, 1044–1065. [Google Scholar] [CrossRef]
  60. Navarro, J.; Muñoz, M.; Contreras, A. Temperature as a factor regulating growth and toxin content in the dinoflagellate Alex-andrium catenella. Harmful Algae 2006, 5, 762–769. [Google Scholar] [CrossRef]
Toxins 13 00021 g001 550
Figure 1. Multiple sequences alignments of the deduced amino acid sequence of C. farreri and P. yessoensis GPxs: (a) GPx1, (b) GPx3, (c) GPx4, and (d) GPx7 with other GPx orthologs deposited in GenBank. The GPx signature motif 1 (blue), GPx signature motif 2 (red) and active site motif (black) were marked by colored frames. The catalytically important residues were indicated by colored down-arrows, as active sites with red, arginine residues with green, and the loop structures with black. Gaps are indicated by dashes to improve the alignments.
Figure 1. Multiple sequences alignments of the deduced amino acid sequence of C. farreri and P. yessoensis GPxs: (a) GPx1, (b) GPx3, (c) GPx4, and (d) GPx7 with other GPx orthologs deposited in GenBank. The GPx signature motif 1 (blue), GPx signature motif 2 (red) and active site motif (black) were marked by colored frames. The catalytically important residues were indicated by colored down-arrows, as active sites with red, arginine residues with green, and the loop structures with black. Gaps are indicated by dashes to improve the alignments.
Toxins 13 00021 g001
Toxins 13 00021 g002 550
Figure 2. Phylogenetic analysis of GPx family genes of C. farreri and P. yessoensis with other selected organisms using the maximum-likelihood (ML) method. The numbers next to the branches indicated the bootstrap value of each internal branch in the tree nodes from 1000 replicates. The CfGPxs and PyGPxs were highlighted using red and blue dots, respectively. The colored branches represent GPx1 clade (yellow), mollusk GPx3 sub-clade (red), vertebrate GPx3 sub-clade (black), GPx4 clade (blue), and GPx7/8 clade (green).
Figure 2. Phylogenetic analysis of GPx family genes of C. farreri and P. yessoensis with other selected organisms using the maximum-likelihood (ML) method. The numbers next to the branches indicated the bootstrap value of each internal branch in the tree nodes from 1000 replicates. The CfGPxs and PyGPxs were highlighted using red and blue dots, respectively. The colored branches represent GPx1 clade (yellow), mollusk GPx3 sub-clade (red), vertebrate GPx3 sub-clade (black), GPx4 clade (blue), and GPx7/8 clade (green).
Toxins 13 00021 g002
Toxins 13 00021 g003 550
Figure 3. Heatmap of the GPx gene expression of C. farreri and P. yessoensis during embryonic developmental stages (a) and in different adult tissues/organs (b) with the log10(RPKM + 1) values.
Figure 3. Heatmap of the GPx gene expression of C. farreri and P. yessoensis during embryonic developmental stages (a) and in different adult tissues/organs (b) with the log10(RPKM + 1) values.
Toxins 13 00021 g003
Toxins 13 00021 g004 550
Figure 4. Temporal expression of scallop GPxs in hepatopancreas and kidneys of C. farreri after A. minutum and A. catenella exposure (a), and of P. yessoensis after A. catenella exposure (b). The heatmaps were drawn based on log2FC of the exposed group to the control group (0 day) (Table S3). Exposure time (1, 3, 5, 10, or 15 days) was displayed above each heatmap. Significant regulation (|log2FC| > 2 and p < 0.05) was indicated with “*”, and very significant regulation (|log2FC| > 2 and FDR < 0.05) was indicated with “**”.
Figure 4. Temporal expression of scallop GPxs in hepatopancreas and kidneys of C. farreri after A. minutum and A. catenella exposure (a), and of P. yessoensis after A. catenella exposure (b). The heatmaps were drawn based on log2FC of the exposed group to the control group (0 day) (Table S3). Exposure time (1, 3, 5, 10, or 15 days) was displayed above each heatmap. Significant regulation (|log2FC| > 2 and p < 0.05) was indicated with “*”, and very significant regulation (|log2FC| > 2 and FDR < 0.05) was indicated with “**”.
Toxins 13 00021 g004
Table
Table 1. General information of the GPx genes in C. farreri and P. yessoensis genomes.
Table 1. General information of the GPx genes in C. farreri and P. yessoensis genomes.
GeneGenomic PositionNo. IntronAmino Acid (aa) GPx Domain Position Isoelectric Point (pI) Molecular Weight (KDa)
CfGPx1-120345.14:285273–29019611451–717.6316.92
CfGPx1-2723835.1:769–175411211–477.6614.22
CfGPx3-116535.52:1096572–110280821301–635.8616.33
CfGPx3-241509.77:1599940–160607021451–735.6216.23
CfGPx3-3a41509.75:1596971–159876221731–246.7019.79
CfGPx3-3b41509.78:1608479–1610027222196–1438.2424.68
CfGPx3-441509.74:1587475–159166121921–246.2322.01
CfGPx457527.11:265310–27418441413–838.9015.00
CfGPx761639.7:153594–1579753280102–2104.3531.83
PyGPx17781.4:46294–5564611431–718.7316.67
PyGPx3-12921.10:304644–31570021411–635.0316.09
PyGPx3-27441.4:20641–2494621411–735.4015.89
PyGPx3-37441.1:18051–1912411081–949.4112.73
PyGPx3-49581.1:21483–2325411081–1008.7612.83
PyGPx3-510997.2:56916–6098822033–1088.3123.38
PyGPx42091.24:1003767–102230741739–979.2620.10
PyGPx77907.21:716289–720446327799–2074.4331.41
Table
Table 2. Copy number of GPx genes in vertebrates and mollusk genomes.
Table 2. Copy number of GPx genes in vertebrates and mollusk genomes.
GPxH. sapiensM. musculusA. carolinensisG. gallusX. tropicalisD. rerioL. giganteaC. gigasC. farreriP. yessoensis
GPx11111120121
GPx21100000000
GPx31101112355
GPx41101021111
GPx51100000000
GPx61100000000
GPx71111111111
GPx81111110000
Total8835474698
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hlaing, S.M.M.; Lou, J.; Cheng, J.; Xun, X.; Li, M.; Lu, W.; Hu, X.; Bao, Z. Tissue-Biased and Species-Specific Regulation of Glutathione Peroxidase (GPx) Genes in Scallops Exposed to Toxic Dinoflagellates. Toxins 2021, 13, 21. https://doi.org/10.3390/toxins13010021

AMA Style

Hlaing SMM, Lou J, Cheng J, Xun X, Li M, Lu W, Hu X, Bao Z. Tissue-Biased and Species-Specific Regulation of Glutathione Peroxidase (GPx) Genes in Scallops Exposed to Toxic Dinoflagellates. Toxins. 2021; 13(1):21. https://doi.org/10.3390/toxins13010021

Chicago/Turabian Style

Hlaing, Sein Moh Moh, Jiarun Lou, Jie Cheng, Xiaogang Xun, Moli Li, Wei Lu, Xiaoli Hu, and Zhenmin Bao. 2021. "Tissue-Biased and Species-Specific Regulation of Glutathione Peroxidase (GPx) Genes in Scallops Exposed to Toxic Dinoflagellates" Toxins 13, no. 1: 21. https://doi.org/10.3390/toxins13010021

APA Style

Hlaing, S. M. M., Lou, J., Cheng, J., Xun, X., Li, M., Lu, W., Hu, X., & Bao, Z. (2021). Tissue-Biased and Species-Specific Regulation of Glutathione Peroxidase (GPx) Genes in Scallops Exposed to Toxic Dinoflagellates. Toxins, 13(1), 21. https://doi.org/10.3390/toxins13010021

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop