Transcriptome analysis of hepatopancreas of Eriocheir sinensis with hepatopancreatic necrosis disease (HPND)

Abstract

Hepatopancreatic necrosis disease (HPND) is a newly emerging disease in the Chinese mitten crab, Eriocheir sinensis, which has resulted in large economic losses. However, the underlying cause of this disease remains unclear. To better understand the pathogenesis and pathogenic mechanism of HPND, we compared the transcriptome differences of the hepatopancreas of E. sinensis with and without HPND. The analysis yielded > 30 million reads for each sample of three test (with HPND) and three control groups (without HPND). We observed 978 downregulated genes and 644 upregulated genes. Among the gene ontology categories “biological process,” “cellular component,” and “molecular function”, the subcategories cellular process, single-organism process, biological regulation, metabolic process, cell part, organelle, organelle part, binding, and catalytic were enriched. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that “metabolism of xenobiotics by cytochrome P450,” “drug metabolism—cytochrome P450,” “chemical carcinogenesis,” and “material metabolism” were the “five” most significantly enriched pathways in the hepatopancreas of E. sinensis with HPND. The results revealed that material metabolic abnormalities and drug effects from the external environment might be associated with HPND in the Chinese mitten crab. Considering the wide use of pyrethroids for pond cleaning in Xinghua city, we speculated that pyrethroids might cause HPND in the Chinese mitten crab. Our study provided useful information about the cause and pathogenetic mechanisms of HPND and could help to prevent this disease in production practice.

About the Authors:

Zongying Yang

Contributed equally to this work with: Zongying Yang, Kun Hu

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

Affiliation: Nanchang Academy of Agricultural Sciences, Nanchang, China

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

Kun Hu

Contributed equally to this work with: Zongying Yang, Kun Hu

Roles Methodology, Software, Visualization

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

Yujie Hou

Roles Software, Visualization

Affiliation: Nanchang Academy of Agricultural Sciences, Nanchang, China

Yulan Wang

Roles Software, Visualization

Affiliation: Nanchang Academy of Agricultural Sciences, Nanchang, China

Yi Yao

Roles Methodology

Affiliation: Nanchang Academy of Agricultural Sciences, Nanchang, China

Xiaoqing Lei

Roles Validation

Affiliation: Nanchang Academy of Agricultural Sciences, Nanchang, China

Baohua Yan

Roles Validation

Affiliation: Nanchang Academy of Agricultural Sciences, Nanchang, China

Qinglong Jiang

Roles Validation

Affiliation: Nanchang Academy of Agricultural Sciences, Nanchang, China

Chunxian Xiong

Roles Formal analysis

Affiliation: Nanchang Academy of Agricultural Sciences, Nanchang, China

Liangqing Xu

Roles Formal analysis

Affiliation: Nanchang Academy of Agricultural Sciences, Nanchang, China

Liugen Zeng

Roles Formal analysis

* E-mail: [email protected]

Affiliation: Nanchang Academy of Agricultural Sciences, Nanchang, China

Introduction

The Chinese mitten crab, Eriocheir sinensis, is one of the most important economic crustacean species in China, and techniques to rear its larvae and its culture facilities were developed in the early 1980s. The Chinese mitten crab is regarded as a delicious and nutritious crustacean by many Chinese consumers because of its nutritional ingredients and delicate flavor [1]. With the increasing demand for the Chinese mitten crab in food markets, the commercial production of E. sinensis is rapidly expanding. As a province hosting large scale E. sinensis farming, the annual output of E. sinensis in Jiangsu Province, China is worth approximately US$ 4 billion [2]. With the development of crab breeding industry, many diseases, such as black gill syndrome, shell disease, edema disease [3], and tremor disease (TD) [4] have emerged and have resulted in serious economic losses. Fortunately, researchers have revealed the etiologies of these diseases and the corresponding prevention and control measures have been developed. However, the disease termed hepatopancreatic necrosis disease (HPND) in E. sinensis is prevalent in Xinghua city, Jiangsu province, and has spread slowly to other main areas of crab rearing nationwide since 2015 [5]. The clinical symptoms are a color change of the hepatopancreas from orange-yellow to light yellow or gray white, muscle atrophy, and edema. Although the diseased crabs do not die immediately, they are of no economic value.

Despite intensive research into HPND, its etiology has not been elucidated clearly. Ding et al. first claimed that the microsporidian fungus Hepatospora eriocheir was the pathogen of HPND in the Chinese mitten crab [5]. They proposed a pathogenic mechanism by which H. eriocheir starved its host by destroying the hepatopancreas by appropriating host resources, resulting in associated changes in crab metabolism and immunity [6]. However, Pan et al. observed hepatopancreatic cells, spermatogonium, gill tissues, and muscle cells of diseased crabs using transmission electron microscopy, and microorganisms such as bacteria, fungi, microsporidia, and viruses were not detected. In addition, when the healthy crabs were fed or injected with diseased tissues, the symptoms of HPND did not appear. Furthermore, HPND symptoms appeared when crabs were cultured in water with a pH of 9.5 to 10, but not when they were cultured in water with different concentrations of avermectin. Therefore, they concluded that HPND in the Chinese mitten crab was not caused by a virus or microsporidian, but caused by water at a high pH or other environmental factors [7]. Shen et al. [8] constructed the meta-transcriptomic libraries of the hepatopancreas from crabs with and without HPND, and no significant differences in viral and microsporidial communities were detected. However, the prevalence of bacteria belonging to the Tenericutes and Actinobacteria phyla increased, whereas the prevalence of bacteria belonging to the Bacteroidetes phylum decreased in crabs with HPND. In addition, the expression of carboxylesterase family genes was significantly upregulated; therefore, they claimed that HPND was not likely to be caused by a microsporidian and speculated that pyrethroid insecticides injured the hepatopancreas, thus leading to abnormalities in metabolism, nutrient absorption, and microbial dysbiosis; however, this required further study. Gao et al. [9] used liquid chromatography-mass spectrometry to analyze the metabolite profile of the hepatopancreas of crabs with HPND, and found that Hepatospora eriocheir was unlikely to be the cause of HPND in the Chinese mitten crab. They revealed that fatty acid metabolic abnormalities and high concentrations of propamocarb might be associated with HPND. Yang et al. [10] sought the etiology of HPND by parasites examination, pathogenic bacteria isolation, challenge experiment, and electron microscopy observation; however, no pathogenic microorganisms were found. The authors suggested that the disease is caused by non-living agents. Generally, there is no unanimous conclusion about the etiology of HPND in the Chinese mitten crab, and all the above speculations require further confirmation.

Gene expression variations in response to external stimuli are sensitive and rapid, and such transcriptional responses might help to explain how an organism responds to a particular situation [11]. Recently, as a highly effective technology to investigate gene expression and to analyze differentially expressed genes and novel transcripts [12], high-throughput RNA-sequencing has been applied widely to study a large variety of invertebrates, such as Crassostrea gigas [13], Litopenaeus vannamei [14], and E. sinensis [15, 16].

In crustaceans, the hepatopancreas is a digestive gland and an immune organ, which plays an important role in the metabolism of xenobiotics and the innate immune system [17–19]. Transcriptome analyses of Litopenaeus setiferus and L. vannamei indicated that the hepatopancreas has a vital role in non-specific immunity, and a cDNA library prepared from the hepatopancreas is more representative than that of hemocytes [17]. For the Chinese mitten crabs with HPND, the changes in the color and morphological characteristics of the hepatopancreas are the most important clinical symptoms. Therefore, the present study aimed to analyze the transcriptome of the hepatopancreas of E. sinensis with and without HPND using Illumina sequencing and bioinformatic analysis. The objective of the study was to annotate functional genes identified using transcriptome analysis, and to evaluate the transcriptomic response of the Chinese mitten crab affected with HPND. Our study will help to understand the etiology and pathogenic mechanism of HPND in E. sinensis and will provide a reference for further research.

Materials and methods

Statement of ethics

This study was carried out in strict accordance with the recommendations and 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. The study design and protocol were approved by the Ethics Committee of Shanghai Ocean University (Approval Number: SHOU-DW-2016-021; SHOU-DW-2019-053). The crabs were anesthetized and dissected on ice; all efforts were made to minimize their suffering.

Maintenance and treatment of E. sinensis

Healthy mature male crabs (average weight, 115.26 ± 6.84 g), with yellow hepatopancreas, and diseased mature male crabs (average weight, 113.28 ± 7.26 g), with clinical symptoms of HPND, 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, in December 2016. These crabs were transported to the laboratory in polystyrene boxes filled with cultivation water, which was aerated during transportation using an aeration pump. Upon arrival, the crabs were anesthetized and dissected on ice to collect the hepatopancreas, which was stored in liquid nitrogen at −80°C for RNA extraction. Healthy crabs and diseased crabs were separated into three groups (six groups in total). Each group consisted of a mixed collection of ten crabs, and the test and control group consisted of three replicates, respectively. A total of 60 crabs were used.

RNA isolation, cDNA library construction, and sequencing

According to the manufacturer’s instructions, total RNA was isolated from hepatopancreas using the TRIzol reagent (Invitrogen, Waltham, MA, USA). DNA contaminants were removed using RNase-free DNase I (Takara Biotechnology, Dalian, China), and the total RNA was eluted in 100 μL of RNase-free MilliQ H2O. The eluted RNA was stored in liquid nitrogen before the next step. A NanoDrop ND-1000 UV–vis Spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA) was used to check the RNA quality. RNA integrity was determined using an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA) and an RNA 6000 Nano LabChip kit (Agilent Technologies). Poly (A) mRNA was isolated using Oligotex mRNA kits (Qiagen) and oligo-dT beads. The mRNA was treated using fragmentation buffer, and the first-strand cDNA was synthesized using the cleaved RNA fragments as templates. RNase H and DNA polymerase I were used to synthesize the second-strand cDNA. A single “A” base was added to these double-stranded cDNA fragments using Klenow 3′–5′ exopolymerase, and the fragments were end-repaired using T4 polynucleotide kinase, T4 DNA polymerase, and the Klenow fragment. Then, the double-stranded fragments were ligated to an index adapter or adapter using T4 quick DNA ligase. Then, the adaptor-ligated fragments were screened out according to their size via agarose gel electrophoresis. The target cDNA fragments were excised from the gel and a PCR reaction was performed to amplify them. After validation using the ABI StepOnePlus Real-Time PCR System (ABI, Foster City, CA, USA) and the Agilent 2100 Bioanalyzer, the cDNA libraries of the corresponding samples were sequenced on a flow cell with the high-throughput mode on an Illumina HiSeq 2500 unit (Illumina, San Diego, CA, USA) [20].

Data processing and quality control

Trim Galore was used to filter out the low-quality reads and remove 3′ adapter sequences. The FastQC software (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/) was used to clean the obtained reads and to evaluate the quality and content of the remaining clean reads. Then, we conducted a comparative analysis with the reference genome of the Chinese mitten crab. For every sample belonging to the control and test groups, sequence alignment with the reference genome sequences was performed using TopHat (version 2.0.12), with default parameters [21].

Identification of the differentially expressed genes (DEGs)

The RSEM software with default parameter settings [22] was used to estimate the expression levels of specific transcripts using fragments per kilobase of transcript per million fragments mapped (FPKM) [23]. The expression levels of transcripts in the control and test groups were transformed using base log2 (FPKM+1) and evaluated. The DESeq software (version 1.14.0) was used to screen out the DEGs and the fold-change of each transcript was calculated [24]. A two-fold change in expression identified a DEG, and an adjusted p-value < 0.05 was considered statistically significant [20].

GO functional annotation and enrichment analysis for DEGs

To investigate the potential functions of these DEGs, the DEGs were firstly annotated using the UniProt database (http://www.uniprot.org/). The Blast2GO suite (https://www.blast2go.com/) was used to obtain the GO annotations against gene ontology terms (GO; http://www.geneontology.org) [25]. All DEGs were mapped to all GO terms in the GO database, and the number of DEGs mapped to each term was calculated. When compared with the transcriptome background, we used the hypergeometric test to identify the enriched GO terms of the DEGs. The formula used was as follows:where, N is the number of genes with GO annotation added; n represents the number of DEGs included in N; M represents the number of genes that were annotated to the specific GO term; and m represents the number of DEGs included in M. The p-value was assigned to each gene and adjusted using the Benjamini and Hochberg’s approach to control the false discovery rate. The corrected p-value of 0.05 was used as the threshold to determine statistical significance. When the corrected p-value was < 0.05, the GO term was considered significantly enriched in the DEGs [20].

Pathway analysis of DEGs

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

Quantitative real-time reverse-transcription PCR verification

Quantitative real-time reverse-transcription-PCR (qRT-PCR) was used to verify the expression levels of the DEGs identified by RNA-Seq analysis. The genes selected for validation by qRT-PCR were involved in enriched GO terms and enriched KEGG pathways. We designed the primers using Primer 5 software. The gene encoding β-actin of the Chinese mitten crab was selected as internal control to normalize the expression levels of the DEGs; all the experiments were conducted in triplicate. The PCR reaction was carried out in a 25 μL volume, comprising 0.5 μL of primers, 9.5 μL of RNase-free H2O, 12.5 μL of SYBR Premix Ex Taq, and 2 μL of cDNA. 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. Melting curve analysis was performed at the end of qRT-PCR to verify the PCR specificity. The expression levels of DEGs were analyzed using the 2—ΔΔ CT method [20].

Results

Illumina sequencing and quality assessment

RNA-Seq was performed using Illumina sequencing to investigate the transcriptional changes in E. sinensis with HPND. In the test groups, 48,219,524; 44,950,578; and 46,421,300 raw reads were obtained; and in the control groups 31,060,184; 31,016,354; and 32,478,004 raw reads were obtained. After the raw sequence reads were filtered and checked, we obtained 46,485,186; 43,436,726; and 44,832,770 trimmed reads in the test groups, and 30,506,938; 30,453,508; and 31,904,838 in the control groups. The trim rates of test groups were 96.40, 96.63, and 96.58%, and in the control groups were 98.22, 98.19, and 98.24%. The average length of the reads was 150 bp (Table 1). These results revealed that the transcriptome sequencing of the Chinese mitten crab was successful, and the trimmed reads of test and control groups could be used for the further analyses.

[Figure omitted. See PDF.]

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

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

Comparison with the reference genome

We compared the trimmed reads with the reference genome sequence of the Chinese mitten crab [26]. The total mapped rates of the reads were 55.92, 56.24, and 56.20% in the three test groups and 58.82, 60.14, and 56.83% in three control groups. The numbers of multiple mapped reads for the three test groups and three control groups were 2,664,557; 2,522,319; and 2,708,791; and 1,353,719; 1,374,223; and 1,162,749, respectively, which accounted for 5.73, 5.81, and 6.04%; and 4.44, 4.51, and 3.64% of the total reads, respectively. The numbers of unmapped reads in the three test groups were 20,491,399; 19,008,186; and 19,634,582; and in the three control groups were 12, 562,134; 12,138,743; and 13,772,587 (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.0228623.t002

Analysis of DEGs

To explore the DEGs of the Chinese mitten crab with HPND, we drew a volcano plot by comparing the differentially expressed genes (DEGs) of each group. The x-axis was the change in gene expression multiple in the different groups, and the y-axis was the statistically significant changes in gene expression (Fig 1). We identified 978 significantly downregulated DEGs and 644 significantly upregulated DEGs in the HPND group compared with those in the healthy group. The results showed that HPND changed the gene expression of the Chinese mitten crab.

[Figure omitted. See PDF.]

Fig 1. Volcano plot of the degree of differences in the expression profile of E. sinensis samples between the test and control groups.

The x-axis is log2 (fold-change), and the y-axis is -log2 (p-value). Red represents significantly upregulated genes, green represents significantly downregulated genes, and every dot represents one gene.

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

GO annotation of DEGs

The 1622 DEGs were classified into 65 GO terms belonging to the three domains: biological process, cellular component, and molecular function (Fig 2 and S1 Fig).

[Figure omitted. See PDF.]

Fig 2. Histogram of enriched subcategories after gene ontology (GO) analysis of the DEGs in E. sinensis samples with HPND.

The x-axis represents the percent of genes. The y-axis represents GO terms involved in the two main ontologies (biological process and molecular function).

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

Lipids are stored in droplets in crabs and represent the major component of the fat body and the main source of metabolic fuel [27]. When et al. [28] reported that the energy of mature crabs is mainly derived from fat, followed by carbohydrates, with proteins providing the smallest proportion of energy. Research showed that prawns mainly used carbohydrates as a source of energy in response to stress [29]. When the existing carbohydrates are insufficient to provide the energy needed by the body, prawns changed their metabolic pathways to gain energy by breaking down fat and proteins [30, 31]. Under long-term stress, the catabolism of the body is strengthened and anabolism is weakened, resulting in the exhaustion of certain nutrients, which can not only disturb the normal physiological state of the body, but also damage the tissues and organs of the body [32, 33]. In the present study, we found that the expression levels of DEGs involved in “lipid metabolic process” and “carbohydrate catabolic process” were mostly downregulated in the GO annotation. These changes to the two energy metabolism processes in crabs might indicate that HPND in Eriocheir sinensis was related to stress of some kind.

KEGG pathway analysis of DEGs

We mapped the DEGs to the KEGG database to further explore their related biological functions and important pathways. We found that some genes participated in multiple pathways, whereas other genes were included in a single pathway. The first nine most-enriched KEGG pathways of the DEGs were metabolism of xenobiotics by cytochrome P450; glycine, serine and threonine metabolism; pentose and glucuronate interconversions; drug metabolism—cytochrome P450; retinol metabolism; tryptophan metabolism; ascorbate and aldarate metabolism; chemical carcinogenesis; and histidine metabolism (Fig 3). Furthermore, we found that the downregulated DEGs were involved in the following pathways: Metabolism of xenobiotics by cytochrome P450; drug metabolism—cytochrome P450; and chemical carcinogenesis.

[Figure omitted. See PDF.]

Fig 3. Histogram of the most enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways of DEGs in E. sinensis samples with HPND.

The x-axis represents the statistical significance of the enrichment. The y-axis represents the KEGG pathway categories.

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

Verification of the differential expression of DEGs

We used qRT-PCR to verify the gene expression profiles of 10 DEGs selected randomly after the GO and KEGG analyses. All the primer sequences of the 10 examined genes are listed in Table 3. The qRT-PCR results revealed no differences compared with the results of RNA-Seq for the 10 genes (Fig 4). These results further confirmed that the results of RNA-Seq were reliable.

[Figure omitted. See PDF.]

Fig 4. Comparison of 10 gene expression levels between the RNA-sequencing (RNA-Seq) data and the quantitative real-time reverse transcription PCR (qRT-PCR) data.

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

[Figure omitted. See PDF.]

Table 3. Oligonucleotide primers for qRT-PCR to validate DEGs expression.

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

Positive values indicate that the expression of the DEG is upregulated. Negative values indicate that the expression of the DEG is downregulated.

Discussion

The GO annotation showed that compared with those in healthy crab, the expression levels of most DEGs involved in the lipid metabolic process and carbohydrate catabolic process in crabs with HPND were downregulated. Furthermore, the KEGG analysis showed that DEGs participating in the biosynthesis of unsaturated fatty acids were downregulated, and the DEGs involved in glycerolipid metabolism were also downregulated, except for two genes (VN_GLEAN_10001404 and VN_GLEAN_10006205). In the enrichment analysis of KEGG pathways, “glycerolipid metabolism” and “biosynthesis of unsaturated fatty acids” were the important pathways related to lipid metabolism in E. sinensis, and the expression levels of genes involved in these two pathways were mostly downregulated. The results indicated that abnormal lipid storage and mobilization occurred in crabs with HPND. In addition to lipid metabolism, DEGs that were enriched in pathways including glycine, serine and threonine metabolism and pentose and glucoronate interconversions were mainly downregulated. These results were consistent with the results of previous studies, in which nutrients metabolic abnormalities, especially fatty acid and carbohydrate metabolism, were proposed to be associated with HPND in the Chinese mitten crab [8, 9].

In the present study, the enriched drug metabolic pathways associated with HPND in the Chinese mitten crab were metabolism of xenobiotics by cytochrome P450, drug metabolism-other enzymes, drug metabolism-cytochrome P450, chemical carcinogenesis, and glutathione metabolism. Interestingly, all the DEGs associated with these pathways were downregulated. Furthermore, we found that the expression levels of carboxylesterase family genes were upregulated, and CYP450 family genes were downregulated. Previous studies have showed that the biodegradation of pyrethroids was mainly catalyzed by carboxylesterases, P450 enzymes, and glutathione-S-transferase (GST) [34]. Carboxylesterase activity could be used as a reliable indicator of pyrethroids in the process of toxicity identification, and plays an important part in removing the toxicity of pyrethroids [35], CYP450 enzymes were inhibited by pyrethroids in vivo and in vitro [36], and pyrethroids affect material metabolism [20]. Considering that previous research failed to find a relationship between HPND in the Chinese mitten crab and viral, bacterial, or other pathogenic organisms, we speculated that pyrethroids, which are widely used for pond cleaning [37] in Xinghua city, might cause HPND in the Chinese mitten crab.

Conclusions

Our study conducted a comparative transcriptome analysis between diseased crabs with HPND and healthy crabs. GO functional annotation analysis and KEGG pathway analysis showed that abnormal metabolism of nutrients occurred in the diseased crabs. Alterations in the expression levels of CYP450 and carboxylesterase family genes, which are involved in xenobiotics or drug associated pathways, suggested that pyrethroids or other drugs were related to HPND in crabs.

In conclusion, our study revealed that under the stress effects of some as-yet-unidentified toxic substances, abnormalities in metabolism and absorption of nutrients appeared in crabs, which subsequently led to HPND. This hypothesis is consistent with that of Shen et al. [8]. However, whether the toxic substances were pyrethroids alone requires further study. Our research provide new information about the etiology and pathogenesis of HPND in the Chinese mitten crab, and further supported the hypothesis that environment factors play important roles in the formation of HPND in the Chinese mitten crab [38].

Supporting information

[Figure omitted. See PDF.]

S1 Fig. Histogram of enriched subcategories after gene ontology (GO) analysis of the DEGs in E. sinensis samples with HPND.

The x-axis represents the GO terms involved in the three main ontologies (biological process, cellular component, and molecular function). The y-axis represents the percent of genes.

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

(TIF)

S1 Table. Differentially expressed genes.

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

(XLS)

S2 Table. GO analysis of the differentially expressed genes.

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

(XLS)

S3 Table. KEGG pathways of all unigenes.

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

(XLS)

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

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

(XLS)

S5 Table. Validation results of 10 differentially expressed genes using RT-PCR.

https://doi.org/10.1371/journal.pone.0228623.s006

(XLS)

Acknowledgments

We thank the anonymous reviewers and editors for their valuable comments and suggestions. We would like to thank the native English speaking scientists of Elixigen Company (Huntington Beach, California) for editing our manuscript.

Citation: Yang Z, Hu K, Hou Y, Wang Y, Yao Y, Lei X, et al. (2020) Transcriptome analysis of hepatopancreas of Eriocheir sinensis with hepatopancreatic necrosis disease (HPND). PLoS ONE 15(2): e0228623. https://doi.org/10.1371/journal.pone.0228623

References

1. Sui LY, Wille M, Cheng XY, Wu XG, Sorgeloos P. Larviculture techniques of Chinese mitten crab Eriocheir sinensis. Aquaculture. 2011; 315: 16–19.

2. 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

3. Wang W. Bacterial diseases of crabs: A review. Journal of Invertebrate Pathology. 2011; 106(1): 1–26.

4. Wang W, Chen J. Ultrastructural study on a novel microsporidian, Endoreticulatus eriocheir sp. nov. (Microsporidia, Encephalitozoonidae), parasite of Chinese mitten crab, Eriocheir sinensis (Crustacea, Decapoda). Journal of Invertebrate Patholog. 2007; 94: 77–83.

5. Ding Z, Meng Q, Liu H, Yuan S, Zhang F, Sun M, et al. First case of hepatopancreatic necrosis disease in pond-reared Chinese mitten crab, Eriocheir sinensis, associated with microsporidian. Journal of Fish Diseases. 2016; 39: 1043–1051. pmid:26850704

6. Ding ZF, Pan J, Huang H, Jiang GC, Chen JQ, Zhu XS, et al. An integrated metabolic consequence of, Hepatospora eriocheir, infection in the Chinese mitten crab, Eriocheir sinensis. Fish Shellfish Immunol. 2018; 72: 443–451. pmid:29146449

7. Pan ZH, Song XH, Hu XL, Xue RY, Cao GL, Zar MS, et al. Pathological Changes and Risk Factors of Hepatopancreas Necrosis Disease of Mitten Crab, Eriocheir Sinensis. Fisheries and Aquaculture Journal. 2017; 8: 3.

8. Shen HY, Zang YN, Song K, Ma YC, Dai TH, Serwadda A. A Meta-Transcriptomics Survey Reveals Changes in the Microbiota of the Chinese Mitten Crab Eriocheir sinensis Infected with Hepatopancreatic Necrosis Disease. Frontiers in Microbiology. 2017; 8: 732. pmid:28491058

9. Gao TH, Xu Y, Wang K, Deng YF, Yang YJ, Lu QP, et al. Comparative LC-MS based non-targeted metabolite profiling of the Chinese mitten crab, Eriocheir sinensis, suffering from hepatopancreatic necrosis disease (HPND). Aquaculture. 2018; 491: 338–345.

10. Yang ZY, Zhang YL, Hu K, Liu LS, Cai HG, Zhang FX, et al. Etiological and histopathological study on hepatopancreatic necrosis syndrome in Eriocheir sinesis. Acta Hydrobiologica Sinica. 2018; 42(1): 17–22.

11. 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.

12. 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

13. 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

14. 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.

15. 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

16. 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

17. 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

18. 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 Chemistry. 1996; 271(45): 28120–28127. pmid:8910426

19. 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

20. Yang ZY, Zhang YL, Jiang YY, Zhu FJ, Zeng LG, Wang YL, et al. Transcriptional responses in the hepatopancreas of Eriocheir sinensis exposed to deltamethrin. PLoS One. 2017; 12(9): e0184581. pmid:28910412

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

22. 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

23. 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

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

25. 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

26. Song LS, Bian C, Luo YJ, Wang LL, You XX, Li J, et al. Draft genome of the Chinese mitten crab, Eriocheir sinensis. GigaScience. 2016; 5(1): 5.

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

28. When XB, Chen LQ, Ai CX. Studies on Standard Metabolism of the Parent Crab Eriocheir Sinensis. Journal of East China Normal University (Natural Science). 2002; 3: 105–109.

29. Pascual C, Ariadna Sánchez, Adolfo Sánchez, Vargas-albores F, Lemoullac G, Rosas G. Haemolymph metabolic variables and immune response in Litopenaeus setiferus adult males: the effect of an extreme temperature. Aquaculture, 2003; 218(1–4): 0–650.

30. Racottais , Hernandezherrera R. Metabolic responses of the white shrimp, Penaeus vannamei, to ambient ammonia. Comparative Biochemistry and Physiology, 2000; 125(4): 437. pmid:10840218

31. Beis I, Newsholme EA. The contents of adenine nucleotides, phosphagens and some glycolytic intermediates in resting muscles from vertebrates and invertebrates. Biochemical Journal, 1975; 152(1): 23–32. pmid:1212224

32. Mackett D B, Tam W H, Fryer J N. Histological changes in insulin-immunoreactive pancreatic β-cells, and suppression of insulin secretion and somatotrope activity in brook trout (Salvelinus fontinalis) maintained on reduced food intake or exposed to acidic environment. Fish Physiology & Biochemistry, 1992; 10(3): 229–243.

33. Weerd J H V, Komen J. The effects of chronic stress on growth in fish: a critical appraisal. Comparative Biochemistry and Physiology, 1998; 120(1): 107–112.

34. 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.

35. 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

36. 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.

37. 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.

38. Chen XW, Wang J, Yue WC, Liu JS, Wang CH. Hepatopancreas transcriptome analysis of Chinese mitten crab (Eriocheir sinensis) with white hepatopancreas syndrome. Fish Shellfish Immunol. 2017; 70: 302–307. pmid:28860074