A selenium-dependent glutathione peroxidase in the Japanese scallop, Mizuhopecten yessoensis: cDNA cloning, promoter sequence analysis and mRNA expression

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.

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

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.

Response of a novel selenium-dependent glutathione peroxidase from thick shell mussel Mytilus coruscus exposed to lipopolysaccharide, copper and benzo[α]pyrene

Glutathione peroxidase (GPx) plays an important antioxidant role in cellular defense against environmental stress. In the present study, a novel selenium-dependent glutathione peroxidase termed McSeGPx firstly identified in thick shell mussel Mytilus coruscus. McSeGPx consists of 197 amino acid residues, characterized with one selenocysteine residue encoded by an opal stop codon TGA, one selenocysteine insertion sequence (SECIS) in the 3' untranslated region (UTR), two active site motifs and one signature sequence motif. McSeGPx transcripts were constitutively expressed in all examined tissues, and were significantly induced in gills and digestive glands with the stimulations of lipopolysaccharide (LPS), copper (Cu) and benzo[α]pyrene (B[α]P). Additionally, rough increases in McSeGPx activity were detected in both tissues under the challenge of LPS, Cu and B[α]P. Collectively, these results suggested that McSeGPx affiliate to selenocysteine dependent GPx (SeGPx) family and might play an important role in mediating the environmental stressors and antioxidant response in M. coruscus.

Molecular characterization and functional analysis of two phospholipid hydroperoxide isoforms from Larimichthys crocea under Vibrio parahaemolyticus challenge

Glutathione peroxidases family is a key role in the antioxidant system in oxybiotic organisms for cell redox homeostasis. One of their members, phospholipid hydroperoxide glutathione peroxidase (GPx4) have unique monomeric structure and can directly react with complex lipid and membrane-bound peroxides under the presence of glutathione(GSH). In this paper, two complete GPx4 cDNAs (designated as LcGPx4a and LcGPx4b) from Larimichthys crocea are identified by rapid amplification of cDNA ends. The cDNA of LcGPx4a was consisted of a 5'-untranslated region (UTR) of 258 bp, a 3'-UTR of 330 bp, and an open reading frame (ORF) of 561 bp encoding 186 amino acid (aa) polypeptides. And the full-length sequence of LcGPx4b was 1164 bp with a 5'-UTR of 34 bp, a 3'-UTR of 551 bp and an ORF of 576 bp encoding a polypeptide of 191 aa residues with a predicted signal peptide of 15 aa. The characteristic selenocysteine insertion (SECIS) sequence was detected in the 3'UTR of the two sequences with 78 bp in length. The conserved active site of selenocysteine (Sec) encoded by TGA was also identified and formed a tetrad functional structure with glutamine, tryptophan, and asparagine in LcGPx4a and LcGPx4b. Two signature site motifs ("LRILAFPSNQFGNQEPG" and "LRILGFPCNQFGGQEPG") were both conserved in the deduced amino acid of LcGPx4a and LcGPx4b. The genomic structure analysis revealed that the two sequences both had 7 exons and 6 introns, and the Sec opal codon and SECIS element were located at the third and seventh exons, respectively. LcGPx4a and LcGPx4b both have a wide distribution in 9 tissues with various relative expression levels and a highest expression pattern in the liver. Under Vibrio parahaemolyticus challenge, their relative expression levels were altered in the liver, spleen, kidney, and head kidney but with different magnitudes and response time. LcGPx4a and LcGPx4b showed a significantly up-regulated trend in the spleen during experimental period. Above results suggested that LcGPx4a and LcGPx4b were two conserved immune molecules and might play a role in the immune response of fish with a tissue-depemdent manners.

Identification and characterization of two selenium-dependent glutathione peroxidase 1 isoforms from Larimichthys crocea

Glutathione peroxidases, a vital family of antioxidant enzymes in oxybiotic organisms, are involved in anti-pathogen immune response. In this study, two complete selenium-dependent glutathione peroxidase 1 cDNAs (designated as LcGPx1a and LcGPx1b) were obtained from the large yellow croaker Larimichthys crocea by rapid amplification of cDNA ends. The full-length sequence of LcGPx1a was 917 bp with a 5'-untranslated region (UTR) of 52 bp, a 3'-UTR of 289 bp, and an open reading frame of 576 bp encoding 191 amino acid (aa) polypeptides. The cDNA of LcGPx1b was composed of 884 bp with a 5'-UTR of 59 bp, a 3'-UTR of 258 bp, and an open reading frame of 567 bp encoding 188 aa polypeptides. The conserved selenocysteine insertion sequence was detected in the 3'-UTR of both isoforms, which can classify types I and II. Protein sequence analysis revealed that both isoforms included a selenocysteine encoded by an opal codon (TGA) and formed the functioning tetrad site with glutamine, tryptophan, and asparagine. Three conservative motifs, including one active site motif ("GKVVLIENVASLUGTT") and two signature site motifs ("LVILGVPCNQFGHQENC" and "V(A/S)WNFEKFLI"), were conserved both in sequence and location. Multiple alignments revealed that they exhibited a high level of identities with GPx1 from other organisms, especially in the abovementioned conserved amino acid sequence motifs. Tissue expression analysis indicated that LcGPx1a and LcGPx1b had a wide distribution in nine tissues with various abundances. The transcript level of LcGPx1a was not significantly different among the nine tissues, whereas that of LcGPx1b was higher in the kidney and head kidney than in the other tissues. After Vibrio parahaemolyticus stimulation, the expression levels of LcGPx1a and LcGPx1b were unanimously altered in the liver, spleen, kidney, and head kidney but with different magnitudes and response time. LcGPx1a and LcGPx1b showed distinct expression trends in the liver, where LcGPx1b was induced and LcGPx1a was depressed in response to pathogen infection. These results indicate that LcGPx1a and LcGPx1b display functional diversities and play crucial roles in mediating the immune response of fish.

Characterization of the TRAF3IP1 gene in Yesso scallop (Patinopecten yessoensis) and its expression in response to bacterial challenge

Tumor necrosis factor receptor-associated factor 3 (TRAF3) is an important adaptor that transmits upstream activation signals to induce innate immune responses. TRAF3 interacting protein 1 (TRAF3IP1) interacts specifically with TRAF3, but its function in innate immunity remains unclear, especially in marine invertebrates. In this study, to better understand the functions of TRAFs in innate immune responses, we identified and characterized the first bivalve TRAF3IP1 gene, PyTRAF3IP1, from Yesso scallop (Patinopecten yessoensis), one of the most important mollusk species for aquaculture. The PyTRAF3IP1 cDNA is 2,367 bp, with an open reading frame of 1,629 bp encoding 542 amino acids. Phylogenetic and protein structural analysis confirmed the gene's identity and revealed that PyTRAF3IP1 was more similar to vertebrate TRAF3IP1s than to those of invertebrates. PyTRAF3IP1 was expressed in all the adult tissues and developmental stages sampled, implying that it plays versatile roles in many biological processes. Furthermore, PyTRAF3IP1 expression was dramatically induced in the acute phase (3-6 h) after infection with both Gram-positive (Micrococcus luteus) and Gram-negative (Vibrio anguillarum) bacteria, even stronger induction being observed after V. anguillarum challenge. This is the first report of the characterization and immune response involvement of TRAF3IP1 in marine invertebrates, and suggests that TRAF3IP1 contributes to innate immunity in bivalves.

Effects of Organic Selenium Supplementation on Growth, Accumulation, Haematology and Histopathology of Juvenile Barramundi (Lates calcarifer) Fed High Soybean Meal Diets

Soybean meal (SBM) has been commonly utilised as a substitute for fishmeal (FM) in the diets of several fish species. However, little is known regarding their effects on trace element availability and thus their importance to fish. The present study employed two feeding trials to evaluate the implications of dietary selenium (Se) on the growth, accumulation, antioxidant, and histopathological responses of juvenile barramundi (Lates calcarifer). In the first trial, each of three basal diets containing 0, 15 and 43 % SBM as replacements for 0, 25 and 75 % of FM protein on an isoproteic and isocalorific basis were either supplemented or not supplemented with 2 mg kg-1 organic Se (OS). In the second trial, the potential effect of OS supplementation in a high SBM diet was investigated in a feeding trial with five experimental diets: 75 % SBM protein as replacement of FM was supplemented with 2, 3, 4, 5 or 7 mg OS kg-1. Growth was independently influenced by the SBM level and the OS supplementation level but not by their interaction. Glutathione peroxidase (GPx) activity, haematocrit, Se accumulation and muscle tissue integrity were significantly enhanced in fish fed on OS-supplemented diets. Furthermore, when high SBM was included in diets, elevated Se tended to lower the barramundi's performance. These findings suggest that dietary supplementation of OS at 2-3 g kg-1 diet is necessary when high plant protein ingredients are incorporated in the diet, in order to maintain better growth and to afford protection against oxidative stress.

Two variants of selenium-dependent glutathione peroxidase from the disk abalone Haliotis discus discus: Molecular characterization and immune responses to bacterial and viral stresses

Glutathione peroxidase (GPx) is an essential member of the antioxidant systems of living organisms and may be involved in immune defense against pathogenic invasion. In the current study, two selenium-dependent glutathione peroxidases (AbSeGPxs) that shared 54.3% identity were identified from the disk abalone Haliotis discus discus. The open reading frames (ORFs) of AbSeGPx-a and AbSeGPx-b coded for 222 and 220 amino acids, respectively, with a characteristic selenocysteine residue encoded by an opal stop codon (TGA). The conserved selenocysteine insertion sequence (SECIS) element was predicted in the 3' untranslated region (UTR) of both isoforms, and they were found to form two stem-loop structures. Amino acid comparison and phylogenetic studies revealed that the AbSeGPxs were closely related to those in other mollusk species and were evolutionarily distinct from those of other taxonomic groups. The SYBR Green qPCR was employed in investigating the transcripts of AbSeGPxs. The expression of AbSeGPxs mRNA was examined in different embryonic developmental stages and differential expression patterns for AbSeGPx-a and AbSeGPx-b were noted. Meanwhile, the highest expression of AbSeGPxs was detected in the hepatopancreas of healthy adult animals. Next, transcriptional levels were profiled in hemocytes of adults to determine the immune responses of AbSeGPxs to microbial infections. The results revealed the significant up-regulation of AbSeGPx-a in a time-dependent manner after bacterial (Listeria monocytogenes and Vibrio parahaemolyticus) and viral (viral hemorrhagic septicemia virus) infections. Consequently, these findings indicate that AbSeGPx-a and AbSeGPx-b might be involved in the embryonic development of disk abalone and the regulation of immune defense system of adult animals.

Expression profiles of selenium dependent glutathione peroxidase and glutathione S-transferase from Exopalaemon carinicauda in response to Vibrio anguillarum and WSSV challenge

A selenium-dependent glutathione peroxidase cDNA was obtained from the ridgetail white prawn Exopalaemon carinicauda (EcGPx) by rapid amplification of cDNA ends (RACE) methods. The full-length cDNA of EcGPx was 946 bp, which contains a 5'-untranslated region (UTR) of 105 bp, 3'-UTR of 280 bp with a poly (A) tail, an open reading frame (ORF) of 561 bp, encoding a 186 amino-acid polypeptide with the predicted molecular weight of 21.35 kDa and estimated isoelectric point of 7.65. It involves a putative selenocysteine (U39) residue which is encoded by an opal codon, (220)TGA(222), and forms an active site with residues Q73 and W141. Sequence analysis revealed that a GPx signature motif 2 ((63)LAFPCNQF(70)), an extra active site motif ((151)WNFEKF(156)), two putative N-glycosylation site ((75)NNT(77) and (107)NGS(109)), and two arginine residues (R89 and R167) were observed in the EcGPx sequence. Comparison of amino acid sequences showed that white shrimp GPx is more closely related to GPx1 and GPx2 subgroups. Quantitative real-time RT-PCR analysis indicated that two glutathione antioxidant enzymes of E. carinicauda, glutathione peroxidase (designated EcGPx) and glutathione S-transferase (designated EcGST) were widely expressed in all the tested tissues, but showed different expression patterns. After Vibrio anguillarum and WSSV challenge, EcGPx and EcGST transcripts both in hemocytes and hepatopancreas increased in the first 6 h and 3 h, respectively. The results suggested that EcGPx and EcGST might be associated with the immune defenses to V. anguillarum and WSSV in E. carinicauda.

De novo characterization of Japanese scallop Mizuhopecten yessoensis transcriptome and analysis of its gene expression following cadmium exposure

Background

Japanese scallop has been cultured on a large-scale in China for many years. However, serious marine pollution in recent years has resulted in considerable loss to this industry. Moreover, due to the lack of genomic resources, limited research has been carried out on this species. To facilitate the understanding at molecular level immune and stress response mechanism, an extensive transcriptomic profiling and digital gene expression (DGE) database of Japanese scallop upon cadmium exposure was carried out using the Illumina sequencing platform.

Results

RNA-seq produced about 112 million sequencing reads from the tissues of adult Japanese scallops. These reads were assembled into 194,839 non-redundant sequences with open reading frame (ORF), of which 14,240 putative amino acid sequences were assigned biological function annotation and were annotated with gene ontology and eukaryotic orthologous group terms. In addition, we identified 720 genes involved in response to stimulus and 302 genes involved in immune-response pathways. Furthermore, we investigated the transcriptomic changes in the gill and digestive gland of Japanese scallops following cadmium exposure using a tag-based DGE system. A total of 7,556 and 3,002 differentially expressed genes were detected, respectively, and functionally annotated with KEGG pathway annotations.

Conclusion

This study provides a comprehensive transcripts sequence resource for the Japanese scallop and presents a survey of gene expression in response to heavy metal exposure in a non-model marine invertebrate via the Illumina sequencing platform. These results may contribute to the in-depth elucidation of the molecular mechanisms involved in bivalve responses to marine pollutants.

Molecular cloning and characterization of cat, gpx1 and Cu/Zn-sod genes in pengze crucian carp (Carassius auratus var. Pengze) and antioxidant enzyme modulation induced by hexavalent chromium in juveniles

Hexavalent chromium (Cr(6+)) is a common pollutant transient metal with high toxicity in the environment. The toxicological effects partly result from oxidative damage due to the production of excessive reactive oxygen species (ROS) in the reductive process of Cr(6+). To explore the influence of ROS induced directly by Cr(6+) on the oxidative stress generation and antioxidant system, the full length cDNAs of antioxidant-related genes cat, gpx1 and Cu/Zn-sod were successfully acquired from pengze crucian carp first and analyzed. Furthermore, the mRNA expression of the antioxidant genes encompassing catalase (cat), copper/zinc superoxide dismutase (Cu/Zn-sod) and glutathione peroxidase (gpx1), antioxidant enzyme activities of CAT, SOD, and GPx and total protein content were further studied in the gill, intestine and liver of pengze crucian carp (Carassius auratus var. Pengze) juveniles upon acute exposure to Cr(6+) at concentrations of 0.1, 1.0, 10 and 100 mg/L for 4 days. Differential significant changes of the antioxidant enzymes and gene expression were observed in different tissues. The findings contribute to better understanding the antioxidant mechanisms induced by Cr(6+) and selecting the organic-specific sensitive biomarkers to monitor the safety of the aquatic ecosystem.

[Cloning and expression analysis of lipopolysaccharide-induced TNF-a(LITAF) of Japanese scallop (Mizuhopecten yessoensis)]

The lipopolysaccharide-induced TNF-alpha factor (LITAF) is an inflammatory cytokine, which plays an important role in innate immunity system. Based on the expressed sequence tag (EST) of Japanese scallop (Mizuhopecten yessoensis), the cDNA of LITAF gene was amplified using rapid amplification of cDNA ends (RACE) approach. Results showed that the full-length cDNA of LITAF is 1 551 bp consisting of a 5' untranslated region (UTR) of 76 bp, a 3' UTR of 1 001 bp, and an open reading frame (ORF) of 474 bp encoding a polypeptide of 157 amino acids, and there is a conserved LITAF domain in amino acid sequences. The estimated molecular mass is 16.99 kDa and the theoretical isoelectric point is 6.24. The total length of LITAF is 3 698 bp, which includes three exons and two introns. Real-time quantitative PCR was carried out to measure LITAF mRNA expression in adult tissues and monitor mRNA expression patterns during embryonic development after bacteria (Vibrio anguillarum) challenged. The expression level of LITAF mRNA was detected in all the adult tissues with the highest in the kidneys. The trochophore owns the highest expression level of LITAF in embryonic development. LITAF expression showed significant difference(P<0.01)between the control and bacteria challenged specimens at 36 h. These results suggest that the LITAF should be a member of the LITAF family that perhaps involved in the innate immune response of Japanese scallop.

Selenium-dependent glutathione peroxidase gene expression during gonad development and its response to LPS and H₂O₂ challenge in Scylla paramamosain

A selenium-dependent glutathione peroxidase cDNA was obtained from green mud crab Scylla paramamosain (SpGPx) by homology PCR technique and rapid amplification of cDNA ends (RACE) methods. The 1135 bp full-length cDNA contains a 9 bp 5'-untranslated region (UTR), an open reading frame (ORF) of 564 bp encoded a deduced protein of 187 amino acids (aa), and a 562 bp 3'-UTR with a 100 bp conserved eukaryotic selenocysteine insertion sequence (SECIS). It involves a putative selenocysteine (Sec⁴⁰, or U⁴⁰) residue which is encoded by an opal codon, ¹²⁷TGA¹²⁹, and forms an active site with residues Q⁷⁴ and W¹⁴². Sequence characterization revealed that SpGPx contain a characteristic GPx signature motif 2 (⁶⁴LAFPCNQF⁷¹), an active site motif (¹⁵²WNFEKF¹⁵⁷), a potential N-glycosylation site (⁷⁶NTT⁷⁸), and two residues (R⁹⁰ and R¹⁶⁸) which contribute to the electrostatic architecture by directing the glutathione donor substrate. Multiple sequence alignment and phylogenetic analysis showed that SpGPx share a high level of identities and closer relationship with other selected invertebrate GPxs and vertebrate GPx1 and GPx2. Molecular modelling analysis results also supported these observations. Real time quantitative PCR analysis revealed that SpGPx was constitutively expressed in 10 selected tissues, and its expression level in gill and testis was higher than that in the other tissues (p < 0.05). The SpGPx expression increased and then declined during ovarian and testicular development implying thatnscrpits yowed that SpGPx might play an important role in gonad development by protecting them from oxidative stress. The expression of SpGPx mRNA was induced by lipopolysaccharide (LPS) and hydrogen peroxide (H₂O₂) in hepatopancreas and haemocytes. The results suggested that SpGPx was implicated in the immune response induced by LPS and H₂O₂.