Transcriptome-Wide Identification of Differentially Expressed Genes in

Sep 27, 2017 - Jiangsu Key Laboratory for Bioresources of Saline Soils, Jiangsu Synthetic Innovation Center for Coastal Bio-agriculture, Jiangsu Provi...
2 downloads 12 Views 2MB Size
Subscriber access provided by PEPPERDINE UNIV

Article

Transcriptome-Wide Identification of Differentially Expressed Genes in Chinese Oak Silkworm Antheraea pernyi in Response to Lead Challenge Zhao-Zhe Xin, Qiu-Ning Liu, Yu Liu, Dai-Zhen Zhang, Zheng-Fei Wang, HuaBin Zhang, Bao-Ming Ge, Chun-Lin Zhou, Xin-Yue Chai, and Bo-Ping Tang J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.7b03391 • Publication Date (Web): 27 Sep 2017 Downloaded from http://pubs.acs.org on September 28, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Journal of Agricultural and Food Chemistry is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 32

Journal of Agricultural and Food Chemistry

1

Transcriptome-Wide Identification of Differentially Expressed

2

Genes in Chinese Oak Silkworm Antheraea pernyi in Response

3

to Lead Challenge

4

Zhao-Zhe Xin1, Qiu-Ning Liu1, *, Yu Liu, Dai-Zhen Zhang, Zheng-Fei Wang, Hua-Bin

5

Zhang, Bao-Ming Ge, Chun-Lin Zhou, Xin-Yue Chai, Bo-Ping Tang*

6

Jiangsu Key Laboratory for Bioresources of Saline Soils, Jiangsu Synthetic

7

Innovation Center for Coastal Bio-agriculture, Jiangsu Provincial Key Laboratory of

8

Coastal Wetland Bioresources and Environmental Protection, School of Ocean and

9

Biological Engineering, Yancheng Teachers University, Yancheng, 224007, PR China.

10 11 12 13 14

* Corresponding author: Bo-Ping Tang and Qiu-Ning Liu

15

Tel/fax: +86 515 88233991

16 17

E-mail: [email protected] (Bo-Ping Tang), [email protected], [email protected] (Qiu-Ning Liu)

18

Address: Jiangsu Key Laboratory for Bioresources of Saline Soils, Jiangsu Synthetic

19

Innovation Center for Coastal Bio-agriculture, Jiangsu Provincial Key Laboratory of

20

Coastal Wetland Bioresources and Environmental Protection, School of Ocean and

21

Biological Engineering, Yancheng Teachers University, Yancheng 224007, PR China.

22

1

These authors contribute to this work equally.

23 24 25 26 27 28 29 1

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

30

ABSTRACT

31

Antheraea pernyi is a commercially-cultivated silk moth and a source of insect food

32

with high-quality protein. Insects suffer oxidative stress on exposure to heavy metals,

33

and reactive oxygen species are cleared by antioxidant enzymes. To gain better

34

understanding of the antioxidant defense system of A. pernyi, we analyzed

35

transcriptomes of pupae after stimulation with lead and phosphate-buffered saline

36

(control). In total, 72,367 unigenes were identified. Gene ontology analysis revealed

37

these DEGs were in 20 biological process subcategories, 19 cellular component

38

subcategories and 18 molecular function subcategories. Clusters of orthologous

39

groups of proteins annotation placed a total of 528 DEGs into 25 categories. Kyoto

40

Encyclopedia of Genes and Genomes enrichment analysis identified antioxidant

41

defense pathways, including “Peroxisome” and “Glutathione metabolism”, which are

42

reported for the first time in A. pernyi. Our study enriches A. pernyi transcriptome

43

databases and provides insight into the heavy metals responses of antioxidant systems

44

of this insect fat bodies.

45

KEYWORDS: Antheraea pernyi, silkworm, heavy metals, lead, antioxidant defense

46

system, transcriptome

47 48 49 50 51

2

ACS Paragon Plus Environment

Page 2 of 32

Page 3 of 32

52

Journal of Agricultural and Food Chemistry

INTRODUCTION

53

Heavy metal pollution is one of the most serious ecological problem on account

54

of burning of mining and fossil fuels, the chaos of urban garbage, the overuse of

55

fertilizers.1 Heavy metals such as arsenic, cadmium, chromium and lead are very

56

common in the soil and water, which could result in growth inhibition and toxicity in

57

insects parts.2 The changes in lifestyles and environments make it easier for insects to

58

be exposed to heavy metals pollution. Exposure to heavy metals can promote the

59

generation of reactive oxygen species (ROS) in insects which is known to induce

60

oxidative stress.3 Therefore, insects have evolved antioxidant defense systems,

61

including antioxidant enzymes. Insects can clear the ROS through the antioxidant

62

enzymes to avoid oxidative damage. The general consensus is that insect larvae

63

contain catalase (CAT), superoxide dismutase (SOD), glutathione reductase (GR) and

64

glutathione S transferase (GSTs) activities.4-7 There have been a number of published

65

studies about antioxidant enzymes present in the species of lepidopteran insect

66

larvae.8-10 However, there are very few reports regarding the heavy metals responses

67

of antioxidant systems in the Antheraea pernyi fat bodies. In addition, transcriptome

68

data gives us a better understanding of insect immune mechanisms.11,12

69

A. pernyi is among the best-known species of wild silkworms. It is widely cultured

70

in China, India and Korea for silk production.13 A. pernyi has been used as an

71

experimental model insect to study breeding technique, development, metabolism

72

regulation, and immune response.14-23 It is also used as a bioreactor to produce

73

recombinant proteins.24 A. pernyi has four life-cycle stages: egg, larva, pupa and adult.

3

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

74

Recently, it is mainly used as a source of insect food (its larva, pupa and adult) in

75

China.25 The pupae of A. pernyi are regarded as a source of high-quality protein-rich

76

food containing all the essential amino acids.26 The larvae of A. pernyi are raised on

77

oak leaves in the field. Thus, the growth and development of larva are highly

78

susceptible to environmental conditions, including heavy metal pollution, foliage,

79

climate, and disease.

80

In this study, transcriptome sequencing libraries were constructed from the fat

81

bodies of A. pernyi from a lead-treated group and a phosphate-buffered saline (PBS)

82

control group. The differentially expressed genes (DEGs) were analyzed and

83

identified in different functional databases. The results enrich the information in A.

84

pernyi gene databases and provide insights into the potential antioxidant mechanisms

85

of A. pernyi.

86

MATERIALS AND METHODS

87

Experimental Insects and Challenge

88

A. pernyi were originally obtained from the Sericultural Research Institute of

89

Henan and reared on the leaves of oak until pupation. The pupae were collected and

90

kept at room temperature. Three individuals were used for lead challenge. Fat bodies

91

was used for experiment materials. Fat bodies were collected at 24 h through

92

dissection. In order to study the influencing mechanism of lead on A. pernyi, the

93

pupae of diapause were chosen to be injected with 10 μl Pb(NO3)2 for lead exposure.

94

The concentration of lead was 20 mg/l. The control pupae were injected with 10 μl

95

PBS solution at 24 h.

4

ACS Paragon Plus Environment

Page 4 of 32

Page 5 of 32

96

Journal of Agricultural and Food Chemistry

RNA Preparation and Library Construction

97

RNA was extracted using Trizol (Invitrogen, USA), based on the manufacturer’s

98

instructions. RNA degradation and contamination were monitored on 1% agarose gels.

99

RNA purity was checked using a NanoPhotometer® spectrophotometer (Implen,

100

Westlake Village, CA, USA). RNA concentration was measured using the Qubit®

101

RNA Assay Kit in a Qubit® 2.0 Fluorometer (Life Technologies, Carlsbad, CA, USA).

102

RNA integrity was assessed using the RNA Nano 6000 Assay Kit and the Agilent

103

Bioanalyzer 2100 system (Agilent Technologies, Santa Clara, CA, USA).

104

A total of 3 μg RNA per sample were used as input material for RNA sample

105

preparation. Sequencing libraries were generated using the NEBNext® Ultra™ RNA

106

Library Prep Kit for Illumina® (New England Biolabs (NEB), Ipswich, MA, USA)

107

following the manufacturer’s recommendations, and index codes were added to

108

attribute sequences to each sample. Briefly, mRNA was purified from total RNA using

109

poly-T oligo-attached magnetic beads. Fragmentation was carried out using divalent

110

cations at elevated temperature in NEBNext First Strand Synthesis Reaction Buffer

111

(5×). First strand cDNA was synthesized using random hexamer primers and

112

M-MuLV Reverse Transcriptase (RNase H). Second strand cDNA synthesis was

113

subsequently performed using DNA Polymerase I and RNase H. Remaining

114

overhangs were converted into blunt ends via exonuclease/polymerase activities. After

115

adenylation of 3ʹ - ends of DNA fragments, NEBNext Adaptors with hairpin loop

116

structure were ligated in preparation for hybridization. To select cDNA fragments of

117

preferentially 150–200 bp in length, the library fragments were purified with the

5

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

118

AMPure XP system (Beckman Coulter, Beverly, MA, USA). Then 3 μl

119

Uracil-Specific Excision Reagent Enzyme (NEB) was applied to size-selected,

120

adaptor-ligated cDNA at 37°C for 15 min, followed by incubation for 5 min at 95°C

121

before PCR. PCR was performed with Phusion High-Fidelity DNA polymerase,

122

universal PCR primers and Index (X) Primer. Finally, PCR products were purified and

123

library quality was assessed.

124

Transcriptome Assembly, Annotation and Function Enrichment

125

Raw data (raw reads) in FASTQ format were firstly processed through in-house

126

Perl scripts. In this step, clean data (clean reads) were obtained by removing reads

127

containing adapter, reads containing poly-N and low-quality reads from the raw data.

128

At the same time, Q20, Q30, G+C-content and sequence duplication level of the clean

129

data were calculated. All downstream analyses were based on clean data with high

130

quality. The left files (read1 files) from all libraries/samples were pooled into one big

131

left.fq file, and right files (read2 files) into one big right.fq file. Transcriptome

132

assembly was accomplished based on the left.fq and right.fq using Trinity software.27

133

Gene function was annotated based on the following databases: NCBI non-redundant

134

protein sequences (NR); Protein family (Pfam); Clusters of Orthologous Groups of

135

proteins (KOG/COG); Swiss-Prot (a manually annotated and reviewed protein

136

sequence database);28 Kyoto Encyclopedia of Genes and Genomes (KEGG);29 and

137

Gene Ontology (GO).30

138

Quantification of Gene Expression Levels and Differential Expression Analysis

139

Gene expression levels were estimated by RSEM31 for each sample. Differential

6

ACS Paragon Plus Environment

Page 6 of 32

Page 7 of 32

Journal of Agricultural and Food Chemistry

140

expression analysis of two conditions/groups was performed using the DESeq R

141

package. DESeq provides statistical routines for determining differential expression in

142

digital gene expression data using a model based on negative binomial distribution.

143

The resulting P-values were adjusted using the Benjamini and Hochberg approach for

144

controlling the false discovery rate (FDR). Genes with an adjusted P-value < 0.05

145

found by DESeq were assigned as differentially expressed.

146

GO and KEGG Pathway Enrichment Analysis

147

GO enrichment analysis of the DEGs was implemented by the topGO R package

148

based on the Kolmogorov-Smirnov test. KEGG32 is a database resource for

149

understanding high-level functions of biological systems, such as the cell, the

150

organism and the ecosystem, from molecular-level information, especially large-scale

151

molecular datasets generated by genome sequencing and other high-throughput

152

experimental technologies (http://www.genome.jp/kegg/). We used KOBAS33

153

software to test the statistical enrichment of DEGs in KEGG pathways.

154

RESULTS AND DISCUSSION

155

Assembly and Splicing

156

A total of 34,909,205 clean reads were obtained from lead-treated samples and

157

25,795,510 clean reads from the control (PBS-treated) samples, giving rise to total

158

clean base numbers of 10,351,801,068 and 7,662,720,366, respectively. Q30

159

was >92%, and the G+C content was approximately 44% in lead-treated samples and

160

43% in PBS-treated samples (Table 1). Using Trinity software, 106,751 transcripts

161

were generated with average length 953.20 bp and N50 length 1,779 bp; 72,367

7

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

162

unigenes were identified with mean length 738.57 bp and N50 length 1,336bp.

163

Among these unigenes, 45,461 (62.82%) were in the range 200–500 bp, 13,202

164

(18.24%) were 500–1000 bp, 7,870 (10.88%) were 1–2 kbp, and 5,833 (8.06%)

165

unigenes were longer than 2 kbp (Table 2). These results show that the data are

166

high-quality and the unigenes can be used for annotation analysis.

167

Functional Annotation and Classification

168

To obtain comprehensive gene function information, we used seven databases to

169

annotate the unigenes, including NR,34 Swiss-Prot, GO, COG,35 KOG,36 eggNOG37

170

and KEGG. In total, 26,426 unigenes were annotated as follows: 22,620 in NR

171

(85.60%); 10,001 in KEGG (37.85%); 9,002 in COG (34.06%); 12,123 in Swiss-Prot

172

(45.88%); 23,914 in eggNOG (90.49%); 17,204 in Pfam (65.10%); 12,602 in GO

173

(47.69%); and 14,961 in KOG (56.61%) (Table 3).

174

Identification and Analysis of DEGs

175

Based on the differential expression analysis, 1,609 significant DEGs ((|log2(fold

176

change)|) > 2, FDR < 0.01) were identified between the lead-challenged and PBS

177

control groups, including 807 upregulated and 802 downregulated unigenes (Fig. 1).

178

The number of DEGs identified in different functional databases were: 369 in COG,

179

692 in GO, 458 in KEGG, 757 in KOG, 900 in Pfam, 694 in Swiss-Prot, 1,126 in

180

eggNOG, and 1,168 in NR (Table 4).

181

Enrichment Analysis of DEGs

182

The functions of the predicted unigenes were classified in KOG, GO and KEGG.

183

All DEGs were classified into three GO categories by topGO,38 including "Biological

8

ACS Paragon Plus Environment

Page 8 of 32

Page 9 of 32

Journal of Agricultural and Food Chemistry

184

Process" (BP), "Cellular Component" (CC) and "Molecular Function" (MF),

185

comprised of 20, 19 and 18 subcategories respectively. Among the various categories

186

of BP, the top three by frequency were “metabolic process” (350), “cellular process”

187

(326) and “single-organism process” (284). Some BP terms were closely related to

188

immune response, for example “immune system process” (22). Within the CC

189

category, “cell” (179), “cell part” (179) and “organelle” (134) were the dominant

190

groups. In the MF classes, “catalytic activity” (371), “binding” (327) and “transporter

191

activity” (55) constituted the top three clusters by frequency (Fig. 2 and Table 5).

192

COG annotation showed that a total of 528 DEGs could be classified into 25

193

categories. Among them, the classification "General function prediction only (R)"

194

represented the largest group (95, 18.00%), followed by "Replication, recombination

195

and repair (L)" (45, 8.52%) and "Carbohydrate transport and metabolism (G)" (44,

196

8.33%). Six DEGs (1.14%) were related to the category "Defense mechanisms (V)".

197

In addition, the function of six DEGs was unknown (Fig. 3 and Table 6).

198

All KEGG annotation pathways are shown in Fig. 4. Four pathways were

199

connected with cellular processes, five were related to environmental information

200

processing, and six were associated with genetic information processing. Thirty-five

201

pathways were relevant to metabolism. The top 20 KEGG pathways are shown in

202

Table 7, which lists the number of DEGs and the total number of genes in each

203

pathway. Among them, “Phagosome”, “Endocytosis”, “Peroxisome” and “Glutathione

204

metabolism” are related to immune responses and antioxidant defence.39-41

205

Peroxisomes contains CAT, SOD and peroxidase, which are antioxidant enzymes. To

9

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

206

avoid oxidative damage, insects clear reactive oxygen species (ROS) using

207

antioxidant enzymes. Glutathione will arise in the course of glutathione metabolism,

208

which is a mechanism for removing toxins.5 Classifications associated with genetic

209

information42 processing were “Ribosome”, “Ubiquitin mediated proteolysis”,

210

“Protein processing in endoplasmic reticulum” and “RNA transport”. “FoxO signaling

211

pathway” is associated with signal transduction. Signal transduction is activated by

212

complex protein-protein interactions between ligands, receptors and kinases, which

213

play an important role in regulating a variety of physiological and metabolic

214

processes.43-46 “Biosynthesis of amino acids”, “Carbon metabolism”, “Fatty acid

215

metabolism”, “Arginine and proline metabolism”, “Glycolysis/Gluconeogenesis”,

216

“Oxidative phosphorylation”, “Glycerolipid metabolism”, “Pentose and glucuronate

217

interconversions”, “Purine metabolism”, “Glycine, serine and threonine metabolism”,

218

and “Starch and sucrose metabolism” were connected with metabolism, which

219

controls production, maintenance and destruction of biomolecules and energy

220

balance.47-50

221

Pathways Related to Antioxidant Defense

222

Figure 5 shows the peroxisome pathway; this pathway is associated with the

223

antioxidant system. In this pathway, peroxin-14 (PEX14), ATP-binding cassette,

224

subfamily D (ALD), member 3 (PMP70), phytanoyl-CoA Hydroxylase (PHYH),

225

acyl-CoA oxidase (ACOX), ATP-binding cassette, subfamily D (ALD), member 2

226

(ABCD) and xanthine dehydrogenase (XDH) are associated with the upregulated

227

genes. The alcohol-forming fatty acyl-CoA reductase (FAR) protein is associated with

10

ACS Paragon Plus Environment

Page 10 of 32

Page 11 of 32

Journal of Agricultural and Food Chemistry

228

the upregulated and downregulated genes. The main function of peroxisome is fatty

229

acid-oxidation, ether phospholipid biosynthesis and antioxidation system. In

230

antioxidation system, peroxisomes perform important functions, including hydrogen

231

peroxide metabolism and glutathione metabolism. CAT and SOD were generated in

232

the process of hydrogen peroxide metabolism. SODs convert O2− into H2O2. CATs

233

and peroxidases convert H2O2 into H2O. Thus, these enzymes work in concert. Two

234

major types of SOD, copper-zinc SOD (CuZnSOD) and manganese SOD (MnSOD),

235

are found in mammalian cells. CuZnSOD is found in the nucleus and cytosol,

236

MnSOD is localized in mitochondria. CAT is localized to the peroxisomes and

237

cytosol.51 peroxisomes are also connected with glutathione metabolism. Glutathione S

238

transferase Kappa 1 (GSTK1) was generated in the process of glutathione metabolism.

239

GSTK1 is a kind of protein enzyme, belonging to GST gene family. Figure 6 shows

240

that pyrimidodiazepine synthase (1.5.4.1) and ornithine decarboxylase (4.1.1.17) are

241

associated with the upregulated genes; The aminopeptidase N (3.4.11.2) and cytosol

242

aminopeptidase (LAP3) are associated with the downregulated genes; The glutathione

243

S-transferase (2.5.1.18) is associated with the upregulated and downregulated genes.

244

Among molecular antioxidants, reduced/oxidized glutathione (GSH/GSSG) plays a

245

crucial role in direct scavenging of radiation-induced ROS or as cofactor for GR and

246

peroxidase in the glutathione metabolism pathway.52 Glutathione peroxidase (GSH-Px)

247

can catalyze GSH to GSSG. In conjunction with NADPH, glutathione reductase can

248

catalyze GSSG to GSH (Fig. 6). GSH can adjust the synthesis of protein and

249

ribonucleotide and is related to the antioxidant capacity of the body.53 GST can be

11

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 12 of 32

250

considered a primary antioxidant enzyme as it can also remove hydroperoxides from

251

cells.51 Molecular and enzymatic antioxidant systems are crucial in determining the

252

survival of insects after heavy metal exposure. Some pathways of antioxidant defense

253

identified here are reported for the first time in A. pernyi. Our results suggest that A.

254

pernyi fat bodies have a strong antioxidant system that protects against heavy metal

255

pollution.

256

AUTHOR INFORMATION

257

Corresponding Author

258

B.P.T. and Q.N.L., E-mail: [email protected] (B.P.T.), [email protected],

259

[email protected] (Q.N.L.), Tel/fax: +86 515 88233991

260

Funding

261

This work was supported by the National Natural Science Foundation of China

262

(31640074 and 31672267), the Natural Science Foundation of Jiangsu Province

263

(BK20160444), the Natural Science Research General Program of Jiangsu Provincial

264

Higher Education Institutions (15KJB240002, 16KJA180008 and 12KJA180009), the

265

Special Guide Fund Project of Agricultural Science and Technology Innovation of

266

Yancheng city (YKN2014022), the Jiangsu Provincial Key Laboratory of Coastal

267

Wetland

268

JLCBE14006), and the Jiangsu Provincial Key Laboratory for Bioresources of Saline

269

Soils (JKLBS2014013 and JKLBS2015004).

270

Notes

271

The authors declare no competing interests.

272

REFERENCES

Bioresources

and

Environmental

Protection

12

ACS Paragon Plus Environment

(JLCBE13002

and

Page 13 of 32

273 274

Journal of Agricultural and Food Chemistry

(1) Boyd, R. S. Heavy metal pollutants and chemical ecology: exploring new frontiers. J. Chem. Ecol. 2010, 36, 46–58.

275

(2) Dubey, S.; Shri, M.; Misra, P.; Lakhwani, D.; Bag, S. K.; Asif, M. H.; Trivedi, P. K.; Tripathi,

276

R. D.; Chakrabarty, D. Heavy metals induce oxidative stress and genome-wide modulation in

277

transcriptome of rice root. Funct. Integr. Genomic. 2014, 14, 401–417.

278 279

(3) Wu, Q. X.; Xia, Q. Effect of heavy metal on antioxidant enzymes in insects. J. Env. Entomol. 2014, 36, 247–251.

280

(4) Ahmad, S.; Pritsos, C. A.; Bowen, S. M.; Heisler, C. R.; Blomquist, G. J.; Pardini, R. S.

281

Antioxidant enzymes of larvae of the cabbage looper moth, Trichoplusia ni: subcellular

282

distribution and activities of superoxide dismutase, catalase and glutathione reductase. Free Radic.

283

Res. Commun. 1988, 4, 403–408.

284 285

(5) Ahmad, S.; Pardini, R. S. Mechanisms for regulating oxygen toxicity in phytophagous insects. Free Radic. Biol. Med. 1990, 8, 401–413.

286

(6) Ahmad, S.; Pritsos, C. A.; Bowen, S. M.; Kirkland, K. E.; Blomquist, G. J.; Pardini, R. S.

287

Activities of enzymes that detoxify superoxide anion and related toxic oxyradicals in Trichoplusia

288

ni. Arch. Insect Biochem. Physiol. 1987, 6, 85–96.

289

(7) Aucoin, R. R.; Philogene, B. J. R.; Arnason, J. T. Antioxidant enzymes as biochemical

290

defenses against phototoxin-induced oxidative stress in three species of herbivorous lepidoptera.

291

Arch. Insect Biochem. Physiol. 1991, 16, 139–152.

292

(8) Zheng, X. Y.; Long, W. M.; Guo, Y. P; Ma, E. B. Effects of cadmium exposure on lipid

293

peroxidation and the antioxidant system in fourth-instar larvae of Propsilocerus akamusi (Diptera:

294

Chironomidae) under laboratory conditions. J. Econ. Entomol. 2011, 104, 827–832.

13

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

295

(9) Augustyniak, M.; Babczyńska, A. Oxidative stress in newly-hatched chorthippus brunneus

296

-the effects of zinc treatment during diapause, depending on the female's age and its origins. Comp.

297

Biochem. Physiol. C. Toxicol. Pharmacol. 2011, 154, 172–179.

298

(10) Cervera, A.; Maymó, A. C.; Martínez-Pardo, R.; Garcerá, M. D. Antioxidant enzymes in

299

Oncopeltus fasciatus (Heteroptera: Lygaeidae) exposed to cadmium. Environ. Entomol. 2003, 32,

300

705–710.

301

(11) Huang, Y. Y.; Shen, C.; Chen, J. X.; He, C. T.; Zhou, Q.; Tan, X.; Yuan, J. G.; Yang, Z. Y.

302

Comparative transcriptome analysis of two Ipomoea aquatica forsk. cultivars targeted to explore

303

possible mechanism of genotype-dependent accumulation of cadmium. J. Agr. Food Chem. 2016,

304

64, 5241–5250.

305

(12) Huang, R.; Huang, Y. J.; Sun, Z. C.; Huang, J. Q.; Wang, Z. J. Transcriptome analysis of

306

genes involved in lipid biosynthesis in the developing embryo of pecan (Carya illinoinensis). J.

307

Agr. Food Chem. 2017, 65, 4223–4236.

308

(13) Liu, Y. Q.; Li, Y. P.; Li, X. S.; Qin, L. The origin and dispersal of domesticated Chinese oak

309

silkworm Antheraea pernyi in China: a reconstruction based on ancient texts. J. Insect Sci. 2010,

310

10, 180.

311

(14) Liu, Q. N.; Zhu, B. J.; Wang, L.; Wei, G. Q.; Dai, L. S.; Lin, K. Z.; Sun, Y.; Qiu, J. F.; Fu, W.

312

W.; Liu, C. L. Identification of immune response-related genes in the Chinese oak silkworm,

313

Antheraea pernyi, by suppression subtractive hybridization. J. Invertebr. Pathol. 2013, 114,

314

313–323.

315

(15) Liu, Q. N.; Zhu, B. J.; Dai, L. S.; Fu, W. W.; Lin, K. Z.; Liu, C. L. Overexpression of small

316

heat shock protein 21 protects the Chinese oak silkworm Antheraea pernyi, against thermal stress.

14

ACS Paragon Plus Environment

Page 14 of 32

Page 15 of 32

317

Journal of Agricultural and Food Chemistry

J. Insect Physiol. 2013, 59, 848–854.

318

(16) You, Q. S.; Li, Q. Q.; Zheng, H. L.; Hu, Z. W.; Zhou, Y.; Wang, B. Discerning silk produced

319

by Bombyx mori from those produced by wild species using an enzyme-linked immunosorbent

320

assay combined with conventional methods. J. Agr. Food Chem. 2017, 65, 7805–7812.

321

(17) Kausar, S.; Abbas, M. N.; Qian, C.; Zhu, B. J.; Sun, Y.; Sun, Y. X.; Wang, L.; Wei, G. Q.;

322

Maqsood, I.; Liu, C. L. Serpin-14 negatively regulates prophenoloxidase activation and expression

323

of antimicrobial peptides in Chinese oak silkworm Antheraea pernyi. Dev. Comp. Immunol. 2017,

324

76, 45–55.

325

(18) Zhang, C. F.; Zhu, B. J.; Dai, L. S.; Liu, C. L.; Luo, X. G. Nucleopolyhedroviruses (NPV)

326

induce the expression of small heat shock protein 25.4 in Antheraea pernyi. Gene 2016, 591,

327

327–332.

328

(19) Qian, C.; Lin, K. Z.; Zhu, B. J.; Liu, C. L.; Wang, L.; Liu, Y.; Liang, D.; Wang, P.; Wei, G. Q.

329

Identification of a C-type lectin from Antheraea pernyi, which may associate with the immune

330

response to pathogens. J. Asia-Pac. Entomol. 2017, 20, 505–512.

331

(20) Qian, C.; Wang, F.; Zhu, B. J.; Wang, L.; Wei, G. Q.; Sun, Y.; Li, S.; Liu, C. L. Identification

332

of a hemolin protein from Actias selene mediates immune response to pathogens. Int.

333

Immunopharmacol. 2017, 42, 74–80.

334

(21) Fang, S. L.; Wang, L.; Fang, Q.; Chen, C.; Zhao, X. S.; Qian, C.; Wei, G. Q.; Zhu, B. J.; Liu,

335

C. L. Characterization and functional study of a cecropin-like peptide from the Chinese oak

336

silkworm, Antheraea pernyi. Arch. Insect Biochem. Physiol. 2017, 94, e21368.

337

(22) Liu, Q. N.; Lin, K. Z.; Yang, L. N.; Dai, L. S.; Wang, L.; Sun, Y.; Qian, C.; Wei, G. Q.; Liu, D.

338

R.; Zhu, B. J.; Liu, C. L. Molecular characterization of an apolipophorin-III gene from the Chinese

15

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

339

oak silkworm, Antheraea pernyi (Lepidoptera: Saturniidae). Arch. Insect Biochem. Physiol. 2015,

340

88, 155–167.

341

(23) Bian, D. D.; Zhao, X. M.; Chen, L.; Tian, J. W.; Liu, Q. N.; Zhou, C. L.; Tang, B. P.

342

Molecular cloning and expression analysis of the highly conserved eukaryotic translation initiation

343

factor 5A (eIF-5A) from Antheraea pernyi. Entomol. Res. 10.1111/1748-5967.12245.

344

(24) Huang, Y. J.; Kobayashi, J.; Yoshimura, T. Genome mapping and gene analysis of Antheraea

345

pernyi nucleopolyhedrovirus for improvement of baculovirus expression vector system. J. Biosci.

346

Bioeng. 2002, 93, 183–191.

347

(25) Chen, M.; Chen, M. M.; Yao, R.; Li, Y.; Wang, H.; Li, Y. P.; Liu, Y. Q. Molecular cloning

348

and characterization of two 12 kDa FK506-binding protein genes in the Chinese oak silkworm,

349

Antheraea pernyi. J. Agr. Food Chem. 2013, 61, 4599–4605.

350 351

(26) Zhou, J.; Han, D. X. Safety evaluation of protein of silkworm (Antheraea pernyi) pupae. Food Chem. Toxicol. 2006, 44, 1123–1130.

352

(27) Grabherr, M. G.; Haas, B. J.; Yassour, M.; Levin, J. Z.; Thompson, D. A.; Amit, I.;

353

Adiconis, X.; Fan, L.; Raychowdhury, R.; Zeng, Q. D.; Chen, Z. H.; Mauceli, E.; Hacohen, N.;

354

Gnirke, A.; Rhind, N.; Palma, F. D.; Birren, B. W.; Nusbaum, C.; Lindblad-Toh, K.; Friedman, N.;

355

Regev, A. Full length transcriptome assembly from RNA Seq data without a reference genome.

356

Nat. Biotechnol. 2011, 29, 644–652.

357 358

(28) Bairoch, A.; Apweiler, R. The SWISS-PROT protein sequence database and its supplement TrEMBL in 2000. Nucleic Acids Res. 2000, 28, 45–48.

359

(29) Kanehisa, M.; Araki, M.; Goto, S.; Hattori, M.; Hirakawa, M.; Itoh, M.; Katayama, T.;

360

Kawashima, S.; Okuda, S.; Tokimatsu, T.; Yamanishi, Y. KEGG for linking genomes to life and

16

ACS Paragon Plus Environment

Page 16 of 32

Page 17 of 32

361

Journal of Agricultural and Food Chemistry

the environment. Nucleic Acids Res. 2007, 36, D480–D484.

362

(30) Götz, S.; García-Gómez, J. M.; Terol, J.; Williams, T. D.; Nagaraj, S. H.; Nueda, M. J.;

363

Robles, M.; Talón, M.; Dopazo, J.; Conesa, A. High-throughput functional annotation and data

364

mining with the Blast2GO suite. Nucleic Acids Res. 2008, 36, 3420–3435.

365 366 367 368

(31) Li, B.; Dewey, C. N. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC bioinformatics 2011, 12, 323. (32) Kanehisa, M.; Goto, S.; Kawashima, S.; Okuno, Y.; Hattori, M. The KEGG resource for deciphering the genome. Nucleic Acids Res. 2004. 32, D277–D280.

369

(33) Mao, X.; Cai, T.; Olyarchuk, J. G.; Wei, L. Automated genome annotation and pathway

370

identification using the kegg orthology (KO) as a controlled vocabulary. Bioinformatics 2005, 21,

371

3787–3793.

372 373 374 375

(34) Deng, Y. Y.; Li, J. Q.; Wu, S. F.; Zhu, Y. P.; Chen, Y. W.; He, F. C. Integrated nr database in protein annotation system and its localization. Computer Engineering 2006, 32, 71–74. (35) Tatusov, R. L.; Galperin, M. Y.; Natale, D. A. The COG database: a tool for genome scale analysis of protein functions and evolution. Nucleic Acids Res. 2000, 28, 33–36.

376

(36) Koonin, E. V.; Fedorova, N. D.; Jackson, J. D.; Jacobs, A. R.; Krylov, D. M.; Makarova, K.

377

S.; Mazumder, R.; Mekhedov, S. L.; Nikolskaya, A. N.; Rao, B. S.; Rogozin, I. B.; Smirnov, S.;

378

Sorokin, A. V.; Sverdlov, A.V.; Vasudevan, S.; Wolf, Y. I.; Yin, J. J.; Natale, D. A. A

379

comprehensive evolutionary classification of proteins encoded in complete eukaryotic genomes.

380

Genome Biol. 2004, 5, 1–28.

381

(37) Huerta-Cepas, J.; Szklarczyk, D.; Forslund, K.; Cook, H.; Heller, D.; Walter, M. C.; Rattei,

382

T.; Mende, D. R.; Sunagawa, S.; Kuhn, M.; Jensen, L. J.; Mering, C. V.; Bork, P. eggNOG 4.5: a

17

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

383

hierarchical orthology framework with improved functional annotations for eukaryotic,

384

prokaryotic and viral sequences. Nucleic Acids Res. 2015, 44, D286–D293.

385

(38) Alexa, A.; Rahnenfuhrer, J. topGO: enrichment analysis for gene ontology. R package

386

version. 2010, 2.

387

(39) Blander, J. M.; Medzhitov, R. Regulation of phagosome maturation by signals from toll-like

388

receptors. Science 2004, 304, 1014–1018.

389

(40) Zhou, M. T.; Qin, Y.; Li, M.; Chen, C.; Chen, X.; Shu, H. B.; Guo, L. Quantitative

390

proteomics reveals the roles of peroxisome-associated proteins in antiviral innate immune

391

responses. Mol. Cell. Proteomics 2015, 14, 2535–2549.

392

(41) Nimmerjahn, F.; Ravetch, J. V. Fc [gamma] receptors as regulators of immune responses.

393 394 395

Nat. Rev. Immunol. 2008, 8, 34–47. (42) Patel, A. Why genetic information processing could have a quantum basis. J. Biosciences 2001, 26, 145–151.

396

(43) Asea, A.; Rehli, M.; Kabingu, E.; Boch, J. A.; Baré, O.; Auron, P. E.; Stevenson, M. A.;

397

Calderwood, S. K. Novel signal transduction pathway utilized by extracellular HSP70 role of

398

Toll-like receptor (TLR) 2 and TLR4. J. Biol. Chem. 2002, 277, 15028–15034.

399 400 401 402 403 404

(44) Kaplan, D. R.; Miller, F. D. Neurotrophin signal transduction in the nervous system. Curr. Opin. Neurobiol. 2000, 10, 381–391. (45) Smeekens, S. Sugar-induced signal transduction in plants. Annu. Rev. Plant Biol. 2000, 51, 49–81. (46) Simons, K.; Toomre, D. Lipid rafts and signal transduction. Nat. Rev. Mol. Cell Bio. 2001, 1, 31–39.

18

ACS Paragon Plus Environment

Page 18 of 32

Page 19 of 32

405 406

Journal of Agricultural and Food Chemistry

(47) Matarese, G.; La, C. A. The intricate interface between immune system and metabolism. Trends Immunol. 2004, 25,193–200.

407

(48) MacIver, N. J.; Jacobs, S. R.; Wieman, H. L.; Wofford, J. A.; Coloff, J. L.; Rathmell, J. C.

408

Glucose metabolism in lymphocytes is a regulated process with significant effects on immune cell

409

function and survival. J. Leukocyte Biol. 2008, 84, 949–957.

410

(49) Saliba, K. J.; Horner, H. A.; Kirk, K. Transport and metabolism of the essential vitamin

411

pantothenic acid in human erythrocytes infected with the malaria parasite Plasmodium falciparum.

412

J. Biol. Chem. 1998, 273, 10190–10195.

413 414 415 416

(50) Sivajothy, M. T.; Tsubaki, Y.; Hooper, R. E. Decreased immune response as a proximate cost of copulation and oviposition in a damselfly. Physiol. Entomol. 1998, 23, 274–277. (51) Wang, Y.; Oberley, L. W.; Murhammer, D. W. Antioxidant defense systems of two lipidopteran insect cell lines. Free Radical Biol. Med. 2001, 30, 1254–1262.

417

(52) Suman, S.; Khan, Z.; Zarin, M.; Chandna, S.; Seth, R. K. Radioresistant Sf9 insect cells

418

display efficient antioxidant defence against high dose γ-radiation. Int. J. Radiat. Biol. 2015, 91,

419

732–741.

420

(53) Pastore, A.; Federici, G.; Bertini, E.; Piemonte, F. Analysis of glutathione: implication in

421

redox and detoxification. Clin. Chim. Acta 2003, 333, 19–39.

422

Table 1. The quality of clean reads.

423

Sample

Clean Reads

Clean bases

GC (%)

Q30 (%)

PBS

25,795,510

7,662,720,366

43.12%

93.98%

Lead

34,909,205

10,351,801,068

43.95%

92.21%

Table 2. The distribution of splicing length.

424

19

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 20 of 32

Length Range

Transcript

Unigene

200-300

31,361(29.38%)

28,027(38.73%)

300-500

23,174(21.71%)

17,434(24.09%)

500-1000

21,887(20.50%)

13,202(18.24%)

1000-2000

16,671(15.62%)

7,870(10.88%)

2000+

13,656(12.79%)

5,833(8.06%)

Total Number

106,751

72,367

Total Length

101,755,331

53,447,744

N50 Length

1,779

1,336

Mean Length

953.20

738.57

425 426

Table 3. Summary statistics of the A. pernyi transcriptome annotation.

427 Annotated databases

Annotated_Number

Percentage (%)

9,002 12,602 10,001 14,961 17,204 12,123 23,914 22,620 26,426

34.06 47.69 37.85 56.61 65.10 45.88 90.49 85.60 100

COG_Annotation GO_Annotation KEGG_Annotation KOG_Annotation Pfam_Annotation Swissprot_Annotation eggNOG_Annotation NR_Annotation All_Annotated

300