Transcriptional responses in the hepatopancreas of Eriocheir sinensis exposed to deltamethrin

Abstract

Deltamethrin is an important pesticide widely used against ectoparasites. Deltamethrin contamination has resulted in a threat to the healthy breeding of the Chinese mitten crab, Eriocheir sinensis. In this study, we investigated transcriptional responses in the hepatopancreas of E. sinensis exposed to deltamethrin. We obtained 99,087,448, 89,086,478, and 100,117,958 raw sequence reads from control 1, control 2, and control 3 groups, and 92,094,972, 92,883,894, and 92,500,828 raw sequence reads from test 1, test 2, and test 3 groups, respectively. After filtering and quality checking of the raw sequence reads, our analysis yielded 79,228,354, 72,336,470, 81,859,826, 77,649,400, 77,194,276, and 75,697,016 clean reads with a mean length of 150 bp from the control and test groups. After deltamethrin treatment, a total of 160 and 167 genes were significantly upregulated and downregulated, respectively. Gene ontology terms “biological process,” “cellular component,” and “molecular function” were enriched with respect to cell killing, cellular process, other organism part, cell part, binding, and catalytic. Pathway analysis using the Kyoto Encyclopedia of Genes and Genomes showed that the metabolic pathways were significantly enriched. We found that the CYP450 enzyme system, carboxylesterase, glutathione-S-transferase, and material (including carbohydrate, lipid, protein, and other substances) metabolism played important roles in the metabolism of deltamethrin in the hepatopancreas of E. sinensis. This study revealed differentially expressed genes related to insecticide metabolism and detoxification in E. sinensis for the first time and will help in understanding the toxicity and molecular metabolic mechanisms of deltamethrin in E. sinensis.

About the Authors:

Zongying Yang

Contributed equally to this work with: Zongying Yang, Yiliu Zhang, Yingying Jiang

Roles Methodology, Writing – original draft, Writing – review & editing

Affiliations National Pathogen Collection Center for Aquatic Animals, Shanghai Ocean University, Shanghai, China, Nanchang Academy of Agricultural Sciences, Nanchang, China

ORCID http://orcid.org/0000-0001-5634-648X

Yiliu Zhang

Contributed equally to this work with: Zongying Yang, Yiliu Zhang, Yingying Jiang

Roles Methodology, Visualization

Affiliation: National Pathogen Collection Center for Aquatic Animals, Shanghai Ocean University, Shanghai, China

Yingying Jiang

Contributed equally to this work with: Zongying Yang, Yiliu Zhang, Yingying Jiang

Roles Methodology, Visualization

Affiliation: National Pathogen Collection Center for Aquatic Animals, Shanghai Ocean University, Shanghai, China

Fengjiao Zhu

Roles Methodology, Visualization

Affiliation: National Pathogen Collection Center for Aquatic Animals, Shanghai Ocean University, Shanghai, China

Liugen Zeng

Roles Formal analysis

Affiliation: Nanchang Academy of Agricultural Sciences, Nanchang, China

Yulan Wang

Roles Formal analysis

Affiliation: Nanchang Academy of Agricultural Sciences, Nanchang, China

Xiaoqing Lei

Roles Formal analysis

Affiliation: Nanchang Academy of Agricultural Sciences, Nanchang, China

Yi Yao

Roles Software

Affiliation: Nanchang Academy of Agricultural Sciences, Nanchang, China

Yujie Hou

Roles Software

Affiliation: Nanchang Academy of Agricultural Sciences, Nanchang, China

Liangqing Xu

Roles Software

Affiliation: Nanchang Academy of Agricultural Sciences, Nanchang, China

Chunxian Xiong

Roles Software

Affiliation: Nanchang Academy of Agricultural Sciences, Nanchang, China

Xianle Yang

Roles Validation

Affiliation: National Pathogen Collection Center for Aquatic Animals, Shanghai Ocean University, Shanghai, China

Kun Hu

Roles Validation

* E-mail: [email protected]

Affiliation: National Pathogen Collection Center for Aquatic Animals, Shanghai Ocean University, Shanghai, China

Introduction

Pyrethroids are synthetic pesticides, which were first derived from the extract of the flowers of Chrysanthemum cinerariaefolium, and they are widely used for agricultural and residential control of pest insects because of their selectivity and low toxicity against non-target organisms such as mammals and birds [1–4]. However, pyrethroids are highly toxic to fishes: LD50 for fishes is 10 to 1000 times less than the corresponding values for mammals and birds [5–8]; this means that pyrethroid contamination can make aquatic animals the targets of pyrethroid intoxication.

Pyrethroids have 2 subgroups because of their chemical structure: type I and II [9]. Deltamethrin,[(S)-cyano-(3-phenoxyphenyl)-methyl](1R,3R)-3-(2,2-dibromoethenyl)-2,2-dimethyl-cyclopropane-1-carboxylate, is a type II pyrethroid. Because of its short half-life in the environment and organisms, deltamethrin has become one of the most commonly used pesticides [10–12]. Many studies on the ecotoxicology of deltamethrin in fishes have been conducted [13–19], and the World Health Organization [20] determined that the lethal concentration of deltamethrin in fishes exposed for 96 h is 0.4–2 μg·L-1. Crustaceans are more sensitive to pyrethroids than fishes [21], and comparative tests have shown that deltamethrin is most toxic to crustaceans [22]. The 96-h LC50 value of deltamethrin for the pink shrimp (Penaeus duorarum) was 0.35 μg·L-1 [23], and Smith and Stratton [24] concluded that shrimps and lobsters are susceptible to all pyrethroids.

The Chinese mitten crab, Eriocheir sinensis, is an important crustacean species in China, and its culture in facilities was started in the early 1980s [25]. Because of an increasing demand for E. sinensis in the food market, commercial production is rapidly expanding, with an annual output worth approximately US$ 4 billion in Jiangsu Province, China [26]. The wide use of deltamethrin has posed a threat to crab breeding, as pyrethroids may be used for pond cleaning [27]. Biodegradation of pyrethroids in crustaceans has been reported [28]; however, there is currently no information on the effects of deltamethrin on E. sinensis, and the exact mode of action of deltamethrin in E. sinensis is unknown. Since alterations in gene expression after external stimuli are rapid, transcriptional responses to pharmaceutical drugs may help us to understand how an organism responds to a particular drug [29]. High-throughput RNA-sequencing (RNA-Seq) has become the conventional and highly effective technology for analyzing gene expression and identifying novel transcripts and differentially expressed genes [30]. RNA-Seq has been widely used to study various invertebrates, such as E. sinensis [31, 32], Litopenaeus vannamei [33], and Crassostrea gigas [34].

The hepatopancreas is not only a digestive gland but also an immune organ, and it plays an important role in the innate immune system and metabolism of xenobiotics. It is the primary site for synthesizing and excreting immune molecules, such as lectin or lectin-related proteins [35], antibacterial peptides [36] and beta-1,3-glucan-binding protein [37]. Transcriptome analysis of L. vannamei and Litopenaeus setiferus by expressed sequence tag analysis and novel gene discovery has revealed that the hepatopancreas plays a vital role in nonspecific immunity, and the cDNA library for the hepatopancreas is more representative than the library for hemocytes [35].

Li et al. [38] used RNA-Seq to construct a non-normalized cDNA library for the hepatopancreas in E. sinensis and identify immune-associated genes after bacterial infection. Huang et al. [39] assembled and annotated a comprehensive de novo transcriptome of E. sinensis hepatopancreas and characterized differentially expressed genes (DEGs) enriched at different molting stages. These studies were performed without a reference genome. However, Song et al. [40] performed, for the first time, genome sequencing, assembly, and annotation for E. sinensis. All the above-mentioned studies provide valuable information for studying important processes in E. sinensis or molecular mechanisms underlying the adaptive changes made by E. sinensis after stimulation with xenobiotics. The present study was conducted to analyze the transcriptome of the hepatopancreas in E. sinensis exposed to deltamethrin by using Illumina sequencing and bioinformatic analysis with the reference genome. The objective of this study was to annotate functional genes by using transcriptome analysis, analyze the short-term molecular effects of deltamethrin on E. sinensis, and evaluate the transcriptomic response of E. sinensis after exposure to deltamethrin. Our study helped in understanding the biological functions of the hepatopancreas of E. sinensis and provided a reference for further research on E. sinensis. Characterization of immune molecules and analysis of insecticide metabolism are crucial for healthy crab breeding and establishing guidelines for farmers to avoid inappropriate use of insecticides.

Materials and methods

Statement of ethics

This study was strictly performed in accordance with the guidelines on the care and use of animals for scientific purposes established by the Institutional Animal Care and Use Committee of Shanghai Ocean University, Shanghai, China.

Maintenance and treatment of E. sinensis

In October 2016, healthy mature male crabs (average weight, 114.04 ± 9.14 g) were collected from the crab breeding base (32.864256110°N, 119.866714447°E) of Tao Huadao Agricultural Development Co., Ltd. in Anfeng Town, Xinghua City, northern Jiangsu Province, China. The crabs were transported to the laboratory in polystyrene boxes filled with cultivation water, which was aerated during transportation with an aeration pump. Upon arrival, the crabs were cultured in experimental tanks (capacity, 90 L) with 20 L of natural water under controlled temperature (23–25°C), pH (7.25 ± 0.25), and dissolved oxygen (6.85 ± 0.15 mg·L-1) for 2 weeks. During the acclimation period, the crabs were fed with a commercial crab feed (37% crude protein) 2 times a day and were starved 1 day before the experiment. After 2 weeks of acclimation, 60 healthy crabs were separated into 6 groups (3 test and 3 control groups); the crabs in the test groups were given a 40 min bath in 3 ppb of deltamethrin (95% purity, Sigma; dissolved in acetone to prepare a stock solution of 25 mg·L-1), and the crabs in the control groups were exposed to the same concentration of acetone as that of deltamethrin in the test groups. Each experiment was performed in triplicate. After 40 min, the crabs were dissected on ice, and the hepatopancreas was collected in liquid nitrogen and stored at -80°C for RNA extraction.

RNA isolation, cDNA library construction, and sequencing

Total RNA was isolated from the hepatopancreas with TRIzol reagent (Invitrogen, USA), according to the manufacturer’s instructions. DNA contaminants were removed using RNase-free DNase I (TaKaRa Biotechnology, Dalian, China), and total RNA was eluted in 100 μL of RNase-free MilliQ H2O. Then, the RNA was stored at -80°C before the next process. The RNA quality was checked with the NanoDrop ND-1000 UV-vis Spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA), and RNA integrity was determined using the RNA 6000 Nano LabChip kit (Agilent Technologies) and Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA). We isolated poly(A) mRNA with oligo-dT beads and Oligotex mRNA kits (Qiagen). The mRNA was treated with fragmentation buffer, and the cleaved RNA fragments were used as templates to synthesize the first-strand cDNA with reverse transcriptase and random hexamer primers. The second-strand cDNA was synthesized with RNase H and DNA polymerase I. These double-stranded cDNA fragments were end-repaired with T4 DNA polymerase, Klenow fragment, and T4 polynucleotide kinase, followed by the addition of a single “A” base with Klenow 3′-5′ exopolymerase. Then, the fragments were ligated with an adapter or index adapter by using T4 quick DNA ligase. The adaptor-ligated fragments were selected on the basis of their size by using agarose gel electrophoresis. The desired cDNA fragment was excised from the gel, and PCR was performed to amplify the fragments. After validation using the Agilent 2100 Bioanalyzer and ABI StepOnePlus Real-Time PCR System (ABI, California, USA), the cDNA library was finally sequenced and constructed on a flow cell by using the high-throughput mode on the Illumina HiSeq 2500 unit (Illumina, San Diego, USA).

Illumina sequencing, data processing, and quality control

The low-quality reads were filtered out and 3′ adapter sequences were removed using Trim Galore. We cleaned the obtained reads by using FastQC software (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/), and content and quality of the remaining clean reads were evaluated. Then, we performed a comparative analysis with the reference genome of E. sinensis. For each sample belonging to the test and control groups, sequence alignment with the reference genome sequences was conducted using TopHat [41].

Identification of DEGs

To evaluate the expression levels of the transcripts in different groups, we used RSEM software with default parameter settings [42] to estimate the expression level (relative abundance) of a specific transcript with fragments per kilobase of transcript per million fragments mapped (FRKM) [43]. The expression level of each transcript was transformed using base log2 (FPKM+1). We used DESeq software to screen the DEGs and calculate the fold change of the transcript [44]. Two-fold changes were used to investigate the expression, and p-values < 0.05 were considered statistically significant.

GO functional annotation and enrichment analysis for DEGs

To analyze the potential functions of the DEGs, we first annotated the DEGs against the UniProt database (http://www.uniprot.org/). Then, we analyzed the functional annotation by using gene ontology terms (GO; http://www.geneontology.org) with Blast2GO (https://www.blast2go.com/) [45]. All the DEGs were mapped to GO terms in the GO database, and the number of genes for every term was calculated. The ultra-geometric test was used to detect the enriched GO terms of DEGs when compared with the transcriptome background. The formula used was as follows:where, N represents the number of genes with GO annotation added; n is the number of DEGs in N; M is the number of genes annotated to specific GO terms; and m is the number of DEGs in M. The calculated p-value was subjected to Bonferroni correction. A corrected p-value of 0.05 was determined as the threshold of statistical significance, and GO terms were considered significantly enriched in the DEGs when the corrected p-value was <0.05.

Pathway analysis of DEGs

Pathways of the DEGs were annotated using blastall (http://nebc.nox.ac.uk/bioinformatics/docs/blastall.html) against the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. We identified the enriched DEG pathways with the same formula as that used for the GO analysis. N represents the number of genes added with KEGG annotation; n represents the number of DEGs in N; M represents the number of genes annotated to specific pathways; and m is the number of DEGs in M.

Quantitative reverse-transcription PCR verification

We used quantitative reverse-transcription (qRT)-PCR to verify the expression levels of DEGs identified using RNA-Seq analysis. The primers were designed with Primer 5 software, and β-actin of E. sinensis was used as the internal control to normalize the expression levels. All the experiments were performed in triplicate. The reaction was performed using a 25 μL volume comprising 2 μL of cDNA, 12.5 μL of SYBR Premix Ex Taq (TaKaRa), 9.5 μL of RNase-free H2O, and 0.5 μL of the forward and reverse primers (10 μmol·L-1). The thermal cycling program was as follows: 95°C for 30 s, followed by 40 cycles of 95°C for 5 s, 60°C for 30 s, and 72°C for 30 s. The melting curve analysis was conducted at the end of qRT-PCR to confirm PCR specificity. The expression levels were analyzed with the 2—ΔΔ CT method.

Results

Illumina sequencing and quality assessment

RNA-Seq was conducted using Illumina sequencing to examine the effects of deltamethrin on E. sinensis transcriptome. The raw sequence reads in the test samples were 92,094,972, 92,883,894, and 92,500,828, and those in the control samples were 99,087,448, 89,086,478, and 100,117,958. After filtering and checking the quality, the trimmed reads in the test samples were 77,649,400, 77,194,276, and 75,697,016, and those in the control samples were 79,228,354, 72,336,470, and 81,859,826. The trim rates of the test samples were 84.31%, 83.11%, and 81.83%, and those of the control samples were 79.96%, 81.20%, and 81.76%; the mean length was 150 bp (Table 1). The results showed that the transcriptome sequencing of E. sinensis was successful, and the trimmed reads could be used for further analyses.

[Figure omitted. See PDF.]

Table 1. Summary of reads obtained from E. sinensis transcriptome sequencing.

https://doi.org/10.1371/journal.pone.0184581.t001

Comparative analysis with the reference genome

The trimmed reads from the hepatopancreas transcriptome of E. sinensis were compared with the reference genome sequence. The total mapped rates of the reads were 63.00%, 64.00%, and 63.00% in the 3 test groups and 63.00%, 64.00%, and 60.00% in the 3 control groups. Multiple mapped reads were about 3 million (3,030,201), 2 million (2,430,849), and 3 million (3,366,219) for the 3 test groups, and 3 million (3,791,979), 2 million (2,430,849), and 3 million (3,201,152) for the 3 control groups, accounting for 4.00%, 3.00%, 4.00%, 5.00%, 5.00%, and 4.00% of the total reads. There were 28,972,914, 28,138,541, and 28,200,925 unmapped reads in the test groups, and 29,285,518, 26,016,203, and 32,444,243 unmapped reads in the control groups (Table 2).

[Figure omitted. See PDF.]

Table 2. Statistical results for mapping the trimmed reads with the reference genome.

https://doi.org/10.1371/journal.pone.0184581.t002

Analysis of DEGs

To investigate the DEGs of E. sinensis exposed to deltamethrin, the Cuffdiff program was used to generate the gene expression profiles (Figs 1 and 2). This program identified 160 significantly upregulated DEGs and 167 significantly downregulated DEGs. The results indicated that deltamethrin affected gene expression in E. sinensis.

[Figure omitted. See PDF.]

Fig 1. Effects of deltamethrin treatment on the gene expression profile of E. sinensis.

Volcanic plot of the degree of differences in the expression profile of E. sinensis. X-axis, log2 (fold change); Y-axis, -log2 (p-value). Red, significantly upregulated genes; green, significantly downregulated genes. Each dot represents one gene.

https://doi.org/10.1371/journal.pone.0184581.g001

[Figure omitted. See PDF.]

Fig 2. Effects of deltamethrin treatment on the gene expression profile using pattern clustering.

Red, upregulated genes; green, downregulated genes. Each line represents one gene.

https://doi.org/10.1371/journal.pone.0184581.g002

GO annotation of DEGs

We classified 327 non-redundant unigenes into 65 GO terms under 3 domains: biological process, cellular component, and molecular function. The dominant subcategories of biological process were “cellular process” (164 genes), “single-organism process” (150 genes), “metabolic process” (131 genes), and “biological regulation” (99 genes). The dominant subcategories of cellular component were “cell part” (201 genes), “organelle” (121 genes), “organelle part” (95 genes), and “membrane” (82 genes). The dominant subcategories of molecular function were “binding” (160 genes), “catalytic” (123 genes), “transporter” (21 genes), and “molecular transducer” (15 genes) (Fig 3).

[Figure omitted. See PDF.]

Fig 3. Histogram of the enriched subcategories after gene ontology (GO) annotation of the differentially expressed genes (DEGs) in E. sinensis specimens that received deltamethrin treatment.

The GO terms (X-axis) were grouped into 3 main ontologies: biological process (BP), cellular component (CC), and molecular function (MF). The Y-axis indicates the number of DEGs.

https://doi.org/10.1371/journal.pone.0184581.g003

KEGG pathway analysis of DEGs

The DEGs were mapped to the KEGG database to further investigate the biological functions and important pathways on the basis of the whole transcriptome background; 327 unigenes were mapped to 262 pathways. Some genes were found in multiple pathways, and some genes were restricted to a single pathway. The first 20 most-enriched KEGG pathways of DEGs were noted; metabolic pathways (31 genes) were the most significantly enriched pathways, followed by glycine, serine, and threonine metabolism (8 genes); ribosome biogenesis in eukaryotes (6 genes); chemical carcinogenesis (5 genes); pyruvate metabolism (4 genes); drug metabolism-cytochrome P450 (3 genes); metabolism of xenobiotics by cytochrome P450 (3 genes); and glycerolipid metabolism (3 genes) (Fig 4).

[Figure omitted. See PDF.]

Fig 4. Histogram of the top 20 most-enriched KEGG pathways of DEGs in E. sinensis after deltamethrin treatment.

Y-axis, KEGG pathway categories; X-axis, statistical significance of the enrichment.

https://doi.org/10.1371/journal.pone.0184581.g004

DEGs involved in metabolic pathways

Metabolic pathways, biosynthesis of amino acids, and biosynthesis of secondary metabolites are all metabolism-related biological pathways. In our study, the lipid metabolism-related pathways included glycerophospholipid metabolism (1 DEG), glycerolipid metabolism (3 DEGs), steroid biosynthesis (2 DEGs), fatty acid biosynthesis (1 DEG), fatty acid elongation (2 DEGs), biosynthesis of unsaturated fatty acids (2 DEGs), synthesis and degradation of ketone bodies (1 DEG), fatty acid degradation (2 DEGs), fat digestion and absorption (1 DEG), and steroid hormone biosynthesis (2 DEGs). The DEGs involved in the lipid metabolism-related pathways were downregulated, except for VN_GLEAN_10002420|K11262 (fatty acid biosynthesis) and VN_GLEAN_10004144|K13368 (steroid hormone biosynthesis). The DEGs involved in glutathione metabolism (3 DEGs) were all downregulated. Besides the effect of deltamethrin on lipid and glutathione metabolism, we found that deltamethrin treatment affected carbohydrate metabolism, including glycolysis/gluconeogenesis (3 DEGs), pentose phosphate pathway (1 DEG), amino sugar and nucleotide sugar metabolism (3 DEGs), and fructose and mannose metabolism (2 DEGs). We should point out that the DEGs of the pentose phosphate pathway were upregulated, which may indicate that the metabolism of deltamethrin is an energy consumption process. We also found that deltamethrin had an effect on protein metabolism: the DEGs involved in ubiquitin-mediated proteolysis (2 DEGs) were downregulated, while the DEGs involved in protein digestion and absorption (2 DEGs) were upregulated.

The drug metabolic pathways were composed of metabolism of xenobiotics by cytochrome P450 (3 DEGs) and drug metabolism-other enzymes (2 DEGs). Interestingly, all the DEGs were downregulated. These findings indicated that deltamethrin treatment not only affected material (carbohydrate, lipid, protein, and other substances) metabolism but also detoxification mechanism of the hepatopancreas.

DEGs involved in signal transduction

We found that many DEGs were associated with signal transduction pathways. The DEGs involved in the calcium signaling pathway (1 DEG), axon guidance (3 DEGs), MAPK signaling pathway (2 DEGs), ABC transporters (1 DEG), glutamatergic synapse (2 DEGs), and GABAergic synapse (2 DEGs) were all upregulated.

Deltamethrin is neurotoxic, and it affects the nervous system by modifying the Na+ channels to remain open, resulting in prolonged Na+ transport along the membranes of nerve cells. The DEGs of axon guidance were all upregulated and may be associated with the toxic mechanism of deltamethrin.

The MAPK signaling pathway can be found in all eukaryotes, and mitogen-activated protein kinase signaling cascades include extracellular signal-regulated kinase, MAPK kinase, and MAPK kinase kinase. Most of the information on the MAPK signaling pathway is focused on vertebrates. In vertebrates, this pathway is significant for anti-stress, inflammation, cell development, reproduction, and differentiation. The MAPK signaling pathway transducts the extracellular signal to the cytoplasm and nucleus by substrate phosphorylation, and the physiological process is controlled by a conserved kinase cascade [46]. However, there is limited information on the MAPK signaling pathway in aquatic invertebrates. To date, Lin et al. [47] have found that the MAPK signaling pathway influenced Trichomonas vaginalis-induced proinflammatory cytokines in the shrimp Penaeus monodon; Feld et al. [48] found that the MAPK signaling pathway participated in neural plasticity in the crab Chasmagnathus; Zhu et al. [49] reported that the MAPK signaling pathway was involved in enrofloxacin metabolism in E. sinensis; and Li et al. [38] found that the MAPK signaling pathway played an important role in microbe-challenged E. sinensis. In our study, we found that the MAPK signaling pathway was associated with deltamethrin metabolism, and the expression levels of the genes involved in the pathway were upregulated. The results showed that the MAPK signaling pathway is involved in xenobiotic metabolism; however, the mechanism of action needs to be studied further.

Differential expression verification of DEGs

qRT-PCR was used to verify the gene expression profiles, and 15 genes were suggested to be related to the deltamethrin treatment after GO and KEGG analyses. The primer sequences for all the examined genes are listed in Table 3. The qRT-PCR verification results for the 15 genes were consistent with the results of RNA-Seq, except for 1 gene (VN_GLEAN_10005348; Fig 5). These results indicate that the RNA-Seq results are generally reliable, and further experiments are required to verify the results of this study.

[Figure omitted. See PDF.]

Fig 5. Comparison of 15 gene expression levels by using RNA-Seq and RT-PCR.

Negative values indicate that the gene expression in E. sinensis was downregulated after deltamethrin treatment; positive values indicate that the gene expression was upregulated.

https://doi.org/10.1371/journal.pone.0184581.g005

[Figure omitted. See PDF.]

Table 3. Oligonucleotide primers of qRT-PCR for validation of DEGs.

https://doi.org/10.1371/journal.pone.0184581.t003

Discussion

Deltamethrin (3 ppb) is often applied using a 30–40 min bath [10, 50]. We designed this study to evaluate the effects of deltamethrin on E. sinensis, which is a very economically important crab species in China. Transcriptome sequencing has become a powerful technique for studying the mechanisms underlying changes in the biological characteristics of an organism. When compared with the traditional methods of analyzing the hepatopancreas of E. sinensis, our study produced more sequencing reads [51, 52] and was useful in yielding genomic resources and molecular information on E. sinensis exposed to deltamethrin.

Detoxification mechanism of the hepatopancreas

Deltamethrin is a type II pyrethroid used against ectoparasites. Biodegradation of pyrethroids is principally catalyzed by P450 enzymes, carboxylesterases, and glutathione-S-transferase (GST) [10]. Biodegradation of deltamethrin occurs through hydrolysis of the central ester via carboxylesterases [9] and oxidation by cytochrome P450 enzymes [53]; however, the metabolites are different in different species [54]. Wheelock et al. [55] found that carboxylesterase activity was important for removing pyrethroid toxicity and could be used as an indicator for pyrethroids during toxicity identification.

Biodegradation of pyrethroids through oxidation by increasing the levels of reactive oxygen species (ROS) has been reported in crustaceans and fishes [15, 28]. Hence, elimination of ROS can alleviate oxidative stress and subsequently decrease toxicological damage [56]. Transcriptomic expression of the genes related to the antioxidant system, including non-enzymatic and enzymatic antioxidants, plays an important role in protecting organisms from damage. It is important to note that GST belongs to the group of transferases responsible for the biodegradation of pesticides, and the activity of this enzyme is associated with the metabolism of insecticides [57]. These transferases catalyze the conjugation of electrophilic compounds to the thiol group of reduced glutathione, producing more soluble products for excretion and, hence, directly decreasing the levels of this antioxidant [58]. We found that the expression levels of CYP450 and glutathione metabolism genes were downregulated, which means that CYP450 and glutamathione are relevant to deltamethrin metabolism. Our results were consistent with those of previous studies: Banka et al. [59] found that CYP450 enzymes were inhibited by deltamethrin in vivo and in vitro, and Abdel-Daim et al. [60] revealed that deltamethrin could induce lipid peroxidation and glutathione reduction. However, Erdogan et al. [61] found a significant increase in CYP450 in the muscle of rainbow trout exposed to deltamethrin, and Olsvik et al. [62] reported that the levels of non-enzymatic antioxidant glutathione were increased in all the tissues exposed to deltamethrin. The reasons for these discrepancies are unclear, but differences in exposure duration, deltamethrin concentration, and different species and tissues may explain the findings.

Candidate genes involved in energy and material metabolism

In our study, we found that the gene involved in the pentose phosphate pathway was upregulated, which indicated that deltamethrin disrupted carbohydrate metabolism. Abdel-Daim et al. [60] showed that serum glucose and cholesterol levels were significantly increased because of deltamethrin toxicity in the Nile tilapia, and El-Sayed & Saad [63] found that a reduction in hepatic glycogen levels accompanied by an increase in plasma glucose levels was a common reaction against toxic insult followed by metabolic stress in fishes. The hyperglycemic effect after pyrethroid treatment in the Nile tilapia suggests the effects of deltamethrin on the glycolytic and glycogenesis pathways [64]. Our study also confirmed that deltamethrin interfered with carbohydrate metabolism. Lipids are stored in droplets and consist of a triglyceride core surrounded by a layer of phospholipids and embedded proteins, and they represent a major component of the fat body and main source of metabolic fuel [65]. Phospholipids are a major component of biomembranes, and most pesticide preparations are characterized by high lipophilicity. The cell membrane is the principal site where pesticides act, and Cengiz et al. [66] found that deltamethrin increased the peroxidation of structurally important unsaturated fatty acids within the phospholipid structure of the membrane, resulting in cell membrane damage. In the enrichment analysis in our study, “Glycerophospholipid metabolism” and “Glycerolipid metabolism” were the important pathways related to lipid metabolism in E. sinensis, and the genes involved in the 2 pathways and biosynthesis of unsaturated fatty acids were all downregulated. We found that deltamethrin not only disrupted lipid storage and mobilization but also influenced protein metabolism. The genes involved in protein digestion and absorption were upregulated, and the genes involved in ubiquitin-mediated proteolysis were downregulated. Thus, deltamethrin disrupted carbohydrate, lipid, and protein metabolism, which would induce abnormalities in metabolism and nutrient absorption.

Effects on signal transduction in deltamethrin

The main mechanism of deltamethrin as a pesticide is believed to result from its binding to a specific receptor site on voltage-dependent Na+ channels, leading to prolonged Na+ transport along the membranes of nerve cells and neural hyperexcitation [5]. In this study, we found that the genes of the Na+- and Cl—dependent GABA transporter were upregulated. Besides, the genes involved in glutamatergic synapse, calcium signaling pathway, and cholinergic synapse were all upregulated.

Previous studies have revealed that many target sites other than the Na+ channel may be linked to the toxicity of deltamethrin. The voltage-dependent Cl- channel has been proposed as a target of pyrethroids [67]. Voltage-sensitive Cl- channels are found in many tissues, and their function is to control membrane excitability. Na+ and Cl- conductance has reciprocal effects on membrane excitability [68]. A previous study has shown that pyrethroids can act on GABA-gated Cl- channels [69], and this effect probably contributes to the seizures that accompany severe type II poisoning. Several other studies have suggested the role of the GABAA receptor–ionophore complex in type II pyrethroid toxicity [70, 71]. Glutamate is an important excitatory neurotransmitter; l-glutamate activates several subtypes of receptors, leading to an increase in intracellular calcium-ion concentration. Many responses to glutamate in the central nervous system can be directly attributed to these changes in membrane polarization and calcium ions. Radcliffe & Dani [72] reported that strong, brief stimulation of nicotinic acetylcholine receptors enhanced hippocampal glutamatergic synaptic transmission on 2 independent time scales and altered the relationship between consecutively evoked synaptic currents. In this study, we also found that deltamethrin could increase glutamatergic synaptic transmission. Acetylcholinesterase activity (AChE) has also been used as a biomarker for pyrethroids [73]; AChE plays a dominant role in cholinergic neurotransmission, hydrolyzing the neurotransmitter acetylcholine at the cholinergic synapses [74]. Tu et al. [28] reported inhibition of AChE in the black tiger shrimp exposed to deltamethrin; we also found that the genes involved in cholinergic synapse were upregulated. Above all, deltamethrin affected the nerve impulse transmissions.

Pyrethroids primarily act on the central nervous system. Researchers may be confused about their exact mechanism of action because pyrethroids interact with not only receptors but also a wide range of ion channels. Nevertheless, it is clear that pyrethroids are likely to interact reversibly with ion channels by their phosphorylation state, and the Na+ channels are a major target.

Conclusions

Because of the rapid metabolism and low toxicity of deltamethrin in non-target animals and humans, it has been widely used in many countries against pests [18, 56]. However, because of its properties, deltamethrin can cause enormous damage to aquatic animals. Alterations in the expression levels of CYP450, carboxylesterase, and GST could disturb the metabolization of deltamethrin, which may cause imbalances in detoxification processes and modifications in the synthesis and degradation of endogenous molecules with key biological activities.

In the present study, deltamethrin suppressed the gene expression levels of the antioxidant enzymes and glutathione. Deltamethrin is neurotoxic, and it also disrupted the metabolism of carbohydrates, lipids, and proteins in E. sinensis, affecting nutrient absorption and metabolization and making E. sinensis vulnerable to pathogen infections and environmental changes. In this study, we have elucidated the genes that were regulated when E. sinensis was exposed to deltamethrin; however, further studies are required to reveal the specific mechanism underlying deltamethrin detoxification, especially with respect to post-transcriptional processes.

Data archiving

The sequencing reads are available in the NCBI SRA database (SRP105235).

Supporting information

[Figure omitted. See PDF.]

S1 Table. Differentially expressed genes.

https://doi.org/10.1371/journal.pone.0184581.s001

(XLS)

S2 Table. GO analysis of the differentially expressed genes.

https://doi.org/10.1371/journal.pone.0184581.s002

(XLS)

S3 Table. KEGG pathways of all unigenes.

https://doi.org/10.1371/journal.pone.0184581.s003

(XLS)

S4 Table. Top 20 most-enriched KEGG pathways of the differentially expressed genes.

https://doi.org/10.1371/journal.pone.0184581.s004

(XLS)

S5 Table. Validation results of 15 differentially expressed genes by RT-PCR.

https://doi.org/10.1371/journal.pone.0184581.s005

(XLS)

Acknowledgments

We would like to thank the native English speaking scientists of Elixigen Company (Huntington Beach, California) for editing our manuscript.

Citation: Yang Z, Zhang Y, Jiang Y, Zhu F, Zeng L, Wang Y, et al. (2017) Transcriptional responses in the hepatopancreas of Eriocheir sinensis exposed to deltamethrin. PLoS ONE 12(9): e0184581. https://doi.org/10.1371/journal.pone.0184581

References

1. Velísek J, Jurcíková J, Dobsíková R, Svobodová Z, Piacková V, Máchová J, et al. Effects of deltamethrin on rainbow trout (Oncorhynchus mykiss). Environ Toxicol Pharmacol. 2007; 23(3): 297–301. pmid:21783771

2. Milam CD, Farris JL, Wilhide JD. Evaluating mosquito control pesticides for effect on target and nontarget organisms. Environ Contam Toxicol. 2000; 39(3): 324–328.

3. Parvez S, Raisuddin S. Copper modulates non-enzymatic antioxidants in the freshwater fish Channa punctata (Bloch) exposed to deltamethrin. Chamosphere. 2006; 62(8): 1324–1332.

4. Prasanthi K, Muralidhara , Rajini PS. Fenvalerate-induced oxidative damage in rat tissues and its attenuation by dietary sesame oil. Food Chem Toxicol. 2005; 43(2): 299–306. pmid:15696646

5. Bradbury SP, Coats JR. Comparative toxicology of the pyrethroid insecticides. Rev Environ Contam Toxicol. 1989; 108(1): 133–177.

6. Casida JE, Gammon DW, Glickman AH, Lawrence LJ. Mechanisms of selective action of pyrethroid insecticides. Annu Rev Pharmacol Toxicol. 1983; 23(23): 413–438.

7. Edwards R, Millburn P, Hutson DH. Comparative toxicity of cis-cypermethrin in rainbow trout, frog, mouse and quail. Toxicol Appl Pharmacol. 1986; 84(3): 512–522. pmid:3487848

8. Paul EA, Simonin HA. Toxicity of three mosquito insecticides to crayfish, Bull. Environ Contam Toxicol. 2006; 76(4): 614–621.

9. Sogorb MA, Vilanova E. Enzymes involved in the detoxification of organophosphorus, carbamate and pyrethroid insecticides through hydrolysis. Toxicology Letters. 2002; 128(3): 215–228.

10. Jacqueline CM, Cristian GE. Deltamethrin (AlphaMax™) reveals modulation of genes related to oxidative stress in the ectoparasite Caligus rogercresseyi: Implications on delousing drug effectiveness. Aquaculture. 2014; 433(37): 421–429.

11. Kan Y, Cengiz EI, Ugurlu P, Yanar M. The protective role of vitamin E on gill and liver tissue histopathology and micronucleus frequencies in peripheral erythrocytes of Oreochromis niloticus exposed to deltamethrin. Environ Toxicol Pharmacol. 2012; 34(2): 170–179. pmid:22534510

12. Villarini M, Moretti M, Pasquini R, Scassellati-Sforzolini G, Fatigoni C, Marcarelli M, et al. In vitro genotoxic effects of the insecticide deltamethrin in human peripheral blood leukocytes: DNA damage (‘comet’ assay) in relation to the induction of sister-chromatid exchanges and micronuclei. Toxicology. 1998; 130(2–3): 129–139. pmid:9865480

13. Guardiola FA, Gónzalez-Párraga P, Meseguer J, Cuesta A, Esteban MA. Modulatory effects of deltamethrin-exposure on the immune status, metabolism and oxidative stress in gilthead seabream (Sparus aurata L.). Fish Shellfish Immunol. 2014; 36(1): 120–129. pmid:24176818

14. Ceyhun SB, Senturk M, Ekinci D, Erdogan O, Ciltas A, Kocaman EM. Deltamethrin attenuates antioxidant defense system and induces the expression of heat shock protein 70 in rainbow trout. Comp Biochem Physiol C Toxicol Pharmacol. 2010; 152(2): 215–223. pmid:20417719

15. Dinu D, Marinescu D, Munteanu MC, Staicu AC, Costache M, Dinischiotu A. Modulatory effects of deltamethrin on antioxidant defense mechanisms and lipid peroxidation in Carassius auratus gibelio liver and intestine. Arch Environ Contam Toxicol. 2010; 58(3): 757–764. pmid:19789906

16. Sakin F, Ispir U, Yonar SM, Yonar ME, Taysi MR. Effect of short-termcypermethrin exposure on oxidant–antioxidant balance in the whole body of rainbow trout fry (Oncorhynchus mykiss). Fresenius Environ Bull. 2011; 20(10): 2806–2809.

17. Yonar ME, Sakin F. Ameliorative effect of lycopene on antioxidant status in Cyprinus carpio during pyrethroid deltamethrin exposure. Pestic Biochem Physiol. 2011; 99(3): 226–231.

18. Amin KA, Hashem KS. Deltamethrin-induced oxidative stress and biochemical changes in tissues and blood of catfish (Clarias gariepinus): antioxidant defense and role of alpha-tocopherol. BMC Vet Res. 2012; 8(1): 45.

19. Ensibi C, Perez-Lopez M, Rodriguez FS, Miguez-Santiyan MP, Yahya MND, Hernandez-Moreno D. Effects of deltamethrin on biometric parameters and liver biomarkers in common carp (Cyprinus carpio). Environ Toxicol Pharmacol. 2013; 36(2): 384–391. pmid:23736096

20. World Health Organization (WHO). Deltamethrin, Environmental Health Criteria 97. World Health Organization, Geneva; 1990. 133 p.

21. Clark JR, Goodman LR, Borthwick PW, Patrick JM, Cripe GM Jr, Moody PM, et al. Toxicity of pyrethroids to marine invertebrates and fish: a literature review and test results with sediment-sorbed chemicals. Environ Toxicol Chem. 1989; 8(5): 393–401.

22. Haya K. Toxicity of pyrethroid insecticides to fish. Environ Toxicol Chem. 1989; 8(5): 381–391.

23. L’ Hotellier M, Vincent P. Assessment of the impact of deltamethrin on aquatic species. In: British Crop Protection Conference-Pests and Disease; 1986. p. 1109–1116.

24. Smith TM, Stratton GW. Effects of synthetic pyrethroid insecticides on nontarget organisms. Res Rev. 1986; 97(97): 93–119.

25. Li XD, Lei Y, Gao XD, Ma CY, Dong SL. Calcium carbonate supersaturation and precipitation in Chinese mitten crab (Eriocheir japonica sinensis) larval ponds in china: mass mortality, crystal form analysis, and safety saturation index. Aquaculture. 2007; 272(1–4): 361–369.

26. Meng Q, Hou L, Zhao Y, Huang X, Huang Y, Xia S, et al. iTRAQ-based proteomic study of the effects of Spiroplasma eriocheiris on Chinese mitten crab Eriocheir sinensis hemocytes. Fish Shellfish Immunol. 2014; 40(1): 182–189. pmid:25017370

27. Manunza B, Deiana S, Pintore M. A molecular modeling study on the interaction between B-cyclodextrin and synthetic pyrethroids. Carbohydrate Research. 1997; 300(1): 89–93.

28. Tu HT, Silvestre F, Meulder BD, Thome JP, Phuong NT, Kestemont P. Combined effects of deltamethrin, temperature and salinity on oxidative stress biomarkers and acetylcholinesterase activity in the black tiger shrimp (Penaeus monodon). Chamosphere. 2012; 86(1): 83–91.

29. Olsvik PA, Lie KK, Hevroy EM. Do anesthetics and sampling strategies affect transcription analysis of fish tissues. BMC Mol Biol. 2007; 8(1): 1–9.

30. Garber M, Grabherr MG, Guttman M, Trapnell C. Computational methods for transcriptome annotation and quantification using RNA-seq. Nature Methods. 2011; 8(6): 469–477. pmid:21623353

31. He L, Wang Q, Jin X, Wang Y, Chen L. Transcriptome profiling of testis during sexual maturation stages in Eriocheir sinensis using Illumina sequencing. PLoS One. 2012; 7(3): e33735. pmid:22442720

32. Ou J, Meng Q, Li Y, Xiu Y, Du J. Identification and comparative analysis of the Eriocheir sinensis microRNA transcriptome response to Spiroplasma eriocheiris infection using a deep sequencing approach. Fish Shellfish Immunol. 2012; 32(2): 345–352. pmid:22166732

33. Djami-Tchatchou AT, Straker CJ, Allie F. 454 Sequencing for the Identification of Genes Differentially Expressed in Avocado Fruit (cv. Fuerte) Infected by Colletotrichum gloeosporioides. Journal of Phytopathology. 2012; 160(9): 449–460.

34. Gavery MR, Roberts SB. Characterizing short read sequencing for gene discovery and RNA-Seq analysis in Crassostrea gigas. Comp Biochem Physiol Part D Genomics Proteomics. 2012; 7(2): 94–99. pmid:22244882

35. Gross P, Bartlett T, Browdy C, Chapman R, Warr G. Immune gene discovery by expressed sequence tag analysis of hemocytes and hepatopancreas in the Pacific White Shrimp, Litopenaeus vannamei, and the Atlantic White Shrimp, L. setiferus. Developmental and Comparative Immunology. 2001; 25(7): 565–577. pmid:11472779

36. Ried C, Wahl C, Miethke T, Wellnhofer G, Landgraf C. High affinity endotoxin-binding and neutralizing peptides based on the crystal structure of recombinant Limulus anti-lipopolysaccharide factor. Journal of Biological Chemsitry. 1996; 271(45): 28120–28127.

37. Roux MM, Pain A, Klimpel KR, Dhar AK. The lipopolysaccharide and beta-1,3-glucan binding protein gene is upregulated in white spot virus-infected shrimp (Penaeus stylirostris). Journal of Virology. 2002; 76(14):7140–7149. pmid:12072514

38. Li X, Cui Z, Liu Y. Transcriptome analysis and discovery of genes involved in immune pathways from hepatopancreas of microbial challenged mitten crab Eriocheir sinensis. PLoS One. 2013; 8(7): e68233. pmid:23874555

39. Song L, Bian C, Luo Y. Draft genome of the Chinese mitten crab, Eriocheir sinensis. GigaScience. 2016; 5(1): 5.

40. Huang S, Wang J, Yue WC, Chen J, Sarah G, Lv WQ, et al. Transcriptomic variation of hepatopancreas reveals the energy metabolism and biological processes associated with molting in Chinese mitten crab, Eriocheir sinensis. Scientific Reports. 2015; 5: 4015.

41. Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009; 25(9): 1105–1111. pmid:19289445

42. He L, Jiang H, Cao D, Liu L, Hu S, Wang Q. Comparative transcriptome analysis of the accessory sex gland and testis from the Chinese mitten crab (Eriocheir sinensis). PLoS One. 2013; 8(1): e53915. pmid:23342039

43. Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc. 2012; 7(3): 562–578. pmid:22383036

44. Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010; 11(10): R106. pmid:20979621

45. Liu B, Jiang GF, Zhang YF, Li JL, Li XJ, Yue JS. Analysis of transcriptome differences between resistant and susceptible strains of the citrus red mite Panonychus citri (Acari: Tetranychidae). PLoS One. 2011; 6(12): e28516. pmid:22162774

46. Ma A, Wang Y, Zhang Z, Han K. MAPK signal transduction pathway and it’s research advance in aquatic invertebrates. Chinese Bulletin of Life Sciences. 2010; 22: 978–984.

47. Lin MC, Hui CF, Chen JY, Wu JL. The antimicrobial peptide, shrimp anti-lipopolysaccharide factor (SALF), inhibits proinflammatory cytokine expressions through the MAPK and NF-kB pathways in Trichomonas vaginalis adherent to HeLa cells. Peptides. 2012; 38(2): 197–207. pmid:23088922

48. Feld M, Dimant B, Delorenzi A, Coso O, Romano A. Phosphorylation of extra-nuclear ERK/MAPK is required for long-term memory consolidation in the crab Chasmagnathus. Behav Brain Res. 2005; 158(2): 251–261. pmid:15698891

49. Zhu FJ, Hu K, Yang ZY, Yang XL. Comparative transcriptome analysis of the hepatopancreas of Eriocheir sinensis following oral gavage with enrofloxacin. Canadian Journal of Fisheries and Aquatic Sciences. 2017; 74(4): 435–444.

50. Cárcamo JG, Aguilar MN, Carreño CF, Vera T, Arias-Darraz L, Figueroa JE, et al. Consecutive emamectin benzoate and deltamethrin treatments affect the expressions and activities of detoxification enzymes in the rainbow trout (Oncorhynchus mykiss). Comparative Biochemistry & Physiology Toxicology & Pharmacology Cbp. 2016; 191:129–137.

51. Jiang H, Cai YM, Chen LQ, Zhang XW, Hu SN, Wang Q. Functional annotation and analysis of expressed sequence tags from the hepatopancreas of mitten crab (Eriocheir sinensis). Mar Biotechnol (NY). 2009; 11(3): 317–326.

52. Jiang H, Yin Y, Zhang X, Hu S, Wang Q. Chasing relationships between nutrition and reproduction: A comparative transcriptome analysis of hepatopancreas and testis from Eriocheir sinensis. Comp Biochem Physiol Part D Genomics Proteomics. 2009; 4(3): 227–234. pmid:20403758

53. Ray DE, Fry JR. A reassessment of the neurotoxicity of pyrethroid insecticides. Pharmacol Ther. 2006; 111(1): 174–193. pmid:16324748

54. Soderlund DM, Clark JM, Sheets LP, Mullin LS, Piccirillo VJ, Sargent D, et al. Mechanisms of pyrethroid neurotoxicity: implications for cumulative risk assessment. Toxicology. 2002; 171(1): 3–59. pmid:11812616

55. Wheelock CE, Miller JL, Miller MJ, Philips BM, Huntley SA, Gee SJ, et al. Use of carboxylesterase activity to remove pyrethroid-associated toxicity to Ceriodaphnia dubia, and Hyalella azteca, in toxicity identification evaluations. Environmental Toxicology and Chemistry. 2006; 25(4): 973–984. pmid:16629136

56. Abdel-Daim MM, Abuzead SM, Halawa SM. Protective role of Spirulina platensis against acute deltamethrin-induced toxicity in rats. PLoS One. 2013; 8(9): e72991. pmid:24039839

57. Hemingway J, Ranson H. Insecticide resistance in insect vectors of human disease. Annual Review of Entomology. 2000; 45(1): 371–391.

58. Berdikova Bohne VJ, Hamre K, Arukwe A. Hepatic metabolism, phase I and II biotransformation enzymes in Atlantic salmon (Salmo Salar, L) during a 12 week feeding period with graded levels of the synthetic antioxidant, ethoxyquin. Food & Chemical Toxicology An International Journal Published for the British Industrial Biological Research Association. 2007; 45(5): 733–46.

59. Banka L, Deér KA, Nemcsók J, Abraham M. In vivo and in vitro effects of deltamethrin on cytochrome P450 monooxygenase activity in carp (Cyprinus carpio L.) liver. Journal of Environmental Science & Health part b Pesticides Food Contaminants & Agricultural Wastes. 1997; 32(5): 789–802.

60. Abdel-Daim MM, Abdelkhalek NK, Hassan AM. Antagonistic activity of dietary allicin against deltamethrin-induced oxidative damage in freshwater Nile tilapia; Oreochromis niloticus. Ecotoxicology & Environmental Safety. 2015; 111: 46–152.

61. Erdogan O, Ceyhun SB, Ekinci D, Aksakal E. Impact of deltamethrin exposure on mRNA expression levels of metallothionein A, B and cytochrome P450 1A in rainbow trout muscles. Gene. 2011; 484(1–4): 13–17.

62. Olsvik PA, Ørnsrud R, Lunestad BT, Steine N, Fredriksen BN. Transcriptional responses in Atlantic salmon (Salmo salar) exposed to deltamethrin, alone or in combination with azamethiphos. Comparative Biochemistry & Physiology Toxicology & Pharmacology Cbp. 2014; 162(1): 23.

63. El-Sayed YS, Saad TT. Subacute Intoxication of a Deltamethrin-Based Preparation (Butox ®; 5% EC) in Monosex Nile Tilapia, Oreochromis niloticus L. Basic Clin Pharmacol Toxicol. 2008; 102(3): 293–299. pmid:18053029

64. El-Sayed YS, Saad TT, El-Bahr SM. Acute intoxication of deltamethrin in monosex Nile tilapia, Oreochromis niloticus with special reference to the clinical, biochemical and haematological effects. Environmental Toxicology & Pharmacology. 2007; 24(3): 212–217.

65. Arrese EL, Soulages JL. Insect fat body: energy, metabolism, and regulation. Annu Rev Entomol. 2010; 55: 207–225. pmid:19725772

66. Cengiz EI, Kan Y, Kizmaz V, Bashan M, Yanar M. The protective role of vitamin E on the fatty acid composition of phospholipid structure in gill and liver tissues of Oreochromis niloticus exposed to deltamethrin. Ecotoxicol Environ Safety. 2012; 80(2): 381–385.

67. Forshaw PJ, Ray DE. A novel action of deltamethrin on membrane resistance in mammalian skeletal-muscle and non-myelinated nerve-fibers. Neuropharmacology. 1990; 29(1): 75–81.

68. Adrian RH, Marshall MW. Action potentials reconstructed in normal and myotonic muscle-fibers. J Physiol. 1976; 258(1): 125–143. pmid:940049

69. Bloomquist JR, Adams PM, Soderlund DM. Inhibition of gamma-aminobutyric acid-stimulated chloride flux in mouse brain vesicles by polychloroalkane and pyrethroid insecticides. Neurotoxicology. 1986; 7(3): 11–20. pmid:2434890

70. Lawrence LJ, Gee KW, Yamamura HI. Interactions of pyrethroid insecticides with chloride ionophore-associated binding sites. Neurotoxicology. 1985; 6(2): 87–98. pmid:2410833

71. Abalis IM, Eldefrawi ME, Eldefrawi AT. Effects of insecticides on GABA-induced chloride influx into rat-brain microsacs. J Toxicol Environ Health. 1986; 18(1): 13–23. pmid:3009836

72. Radcliffe KA, Dani JA. Nicotinic stimulation produces multiple forms of increased glutamatergic synaptic transmission. Journal of Neuroscience the Official Journal of the Society for Neuroscience. 1998; 18(18): 7075–7083.

73. Tu HT, Silvestre F, Scippo ML, Thome JP, Phuong NT, Kestemont P. Acetylcholinesterase activity as a biomarker of exposure to antibiotics and pesticides in the black tiger shrimp (Penaeus monodon). Ecotoxicology & Environmental Safety. 2009; 72(5): 1463.

74. Talesa V, Contenti S, Principato GB, Pascolini R, Giovannini E, Rosi G. Cholinesterases from Maia verricosa and Palinurus vulgaris: a comparative study. Comp Biochem Physiol. 1992; 101C: 499–503.