Cas9-mediated deletion of one carotenoid isomerooxygenase gene (EcNinaB-X1) from Exopalaemon carinicauda

Cas9-mediated deletion of one carotenoid isomerooxygenase gene (EcNinaB-X1) from Exopalaemon carinicauda

Journal Pre-proof CRISPR/Cas9-mediated deletion of one carotenoid isomerooxygenase gene (EcNinaB-X1) from Exopalaemon carinicauda Yuying Sun, Congcong...

2MB Sizes 0 Downloads 11 Views

Journal Pre-proof CRISPR/Cas9-mediated deletion of one carotenoid isomerooxygenase gene (EcNinaB-X1) from Exopalaemon carinicauda Yuying Sun, Congcong Yan, Mengfei Liu, Yujie Liu, Wenzheng Wang, Wei Cheng, Fusheng Yang, Jiquan Zhang PII:

S1050-4648(19)31166-0

DOI:

https://doi.org/10.1016/j.fsi.2019.12.037

Reference:

YFSIM 6679

To appear in:

Fish and Shellfish Immunology

Received Date: 23 October 2019 Revised Date:

9 December 2019

Accepted Date: 13 December 2019

Please cite this article as: Sun Y, Yan C, Liu M, Liu Y, Wang W, Cheng W, Yang F, Zhang J, CRISPR/ Cas9-mediated deletion of one carotenoid isomerooxygenase gene (EcNinaB-X1) from Exopalaemon carinicauda, Fish and Shellfish Immunology (2020), doi: https://doi.org/10.1016/j.fsi.2019.12.037. This is a PDF file of an article that has undergone enhancements after acceptance, such as the addition of a cover page and metadata, and formatting for readability, but it is not yet the definitive version of record. This version will undergo additional copyediting, typesetting and review before it is published in its final form, but we are providing this version to give early visibility of the article. Please note that, during the production process, errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain. © 2019 Published by Elsevier Ltd.

Graphical Abstract

1

CRISPR/Cas9-mediated deletion of one carotenoid isomerooxygenase

2

gene (EcNinaB-X1) from Exopalaemon carinicauda

3 Yuying Sun a, Congcong Yan a, Mengfei Liu a, Yujie Liu a, Wenzheng Wang a, Wei Cheng a,

4

Fusheng Yang b, Jiquan Zhang a,

5

b, *

6 a

7

Laboratory of Zoological Systematics and Application of Hebei Province, College of Life

8

Sciences, Hebei University, Baoding 071002, China

9

b

Xiaoshan Donghai Aquaculture Co., Ltd, Xiaoshan 310012, China

11

*

Corresponding author:

12

E-mail: [email protected]

10

13 14 15

Running title: CRISPR/Cas9-mediated deletion of EcNinaB-X1 gene

16

Abstract During the immune defense reaction of invertebrate, a plenty of reactive oxygen species

17

(ROS) could be induced to product. Though ROS can kill foreign invaders, the accumulation of

18

these reactive molecules in animals will cause serious cell damage. Carotenoids could function as

19

scavengers of oxygen radicals. In this research, cDNA and genomic DNA of one carotenoid

20

isomerooxygenase gene (named EcNinaB-X1) were cloned from Exopalaemon carinicauda.

21

EcNinaB-X1 gene was composed of 12 exons and 11 introns. EcNinaB-X1 knock-out (KO) prawns

22

were produced via CRISPR/Cas9 technology and the change of their phenotypes were analyzed.

23

Of the 400 injected one-cell stage embryos with cas9 mRNA and one sgRNA targeting the first

24

exon of EcNinaB-X1 gene, 26 EcNinaB-X1-KO prawns were generated and the mutant rate

25

reached 6.5% after embryo injection. The EcNinaB-X1-KO prawns had significant lower mortality

26

than those in wild-type group when the prawns were challenged with Vibrio parahaemolyticus or

27

Aeromonas hydrophila. In conclusion, we first demonstrate the function of the carotenoid

28

isomerooxygenase gene in immune defense of E. carinicauda by performing directed, heritable

29

gene mutagenesis.

30 31 32

Keywords: Exopalaemon carinicauda; carotenoid isomerooxygenase; CRISPR/Cas9

33

1. Introduction

34

In 2013, clustered regularly interspaced short palindromic repeats (CRISPR)/ CRISPR

35

associated protein (Cas9) technology was reported to be used in the genome editing of organisms

36

[1, 2]. Since then, CRISPR/Cas9 technology has been widely used in the genome editing of many

37

organisms including mammals [2-4], fly fruit [5, 6], zebrafish [7, 8], prawn [9-11] and so on. It

38

can produce targeted double-strand breaks (DSBs) in the genome and the resulting DSBs are

39

repaired by nonhomologous end-joining (NHEJ) or homology-directed repair (HDR) pathway,

40

thereby causing mutations [12].

41

A lot of crustaceans, such as crayfish, crabs, lobsters, prawns, and shrimp, are important

42

aquaculture species. At present, the genome draft of some crustaceans including Litopenaeus

43

vannamei [13], Exopalaemon carinicauda [14], Neocaridina denticulate [15], and Eriocheir

44

sinensis [16] has been reported. CRISPR/Cas9 technology was a useful tool to clarify the function

45

of genes by loss-of-function approaches in vivo. The ridgetail white prawn E. carinicauda, one of

46

the important commercial shrimp species, naturally distributed in the coasts of China [9]. In our

47

previous research, we used it as a model animal of crustacean in basic research and successfully

48

performed the genome editing using CRISPR/Cas9 technology [11]. We developed a highly

49

efficient microinjection method in E. carinicauda and deleted some interest genes using

50

CRISPR/Cas9 technology [9, 10, 17]. In addition, the low-coverage sequencing and de novo

51

assembly of E. carinicauda genome had also been finished, which made it possible to promote the

52

basic research of crustaceans [14].

53

Carotenoids such as beta-carotene, lycopene, lutein and β-cryptoxanthine are produced in

54

plants, certain bacteria, algae and fungi, where they function as accessory photosynthetic pigments

55

and as scavengers of oxygen radicals for photoprotection [18]. Recently, since carotenoid

56

accumulation played a key role in the formation of colorful animals, more researches have focused

57

on the change of body color in aquaculture animals [19]. All animals cannot produce these

58

naturally-occurring carotenoids, so they must obtain them from their diet. There are some

59

carotenoid oxygenases, including carotenoid isomerooxygenase (NinaB), β, β-carotene 15,

60

15'-monooxygenase (BCMO) and β, β-carotene 9', 10'-oxygenase (BCO2), which can cleave a

61

variety of carotenoids into a range of biologically important products [20, 21]. Then, the cleaved

62

products function as hormones, pigments, flavors, floral scents and defense compounds in animals.

63

In addition, during the immune defense reaction of invertebrate, a plenty of reactive oxygen

64

species (ROS) could be induced to product. Though ROS can kill foreign invaders, the

65

accumulation of these reactive molecules in animals will cause serious cell damage. Carotenoids

66

could function as scavengers of oxygen radicals [22-24].

67

At present, there is no report about the function of carotenoid oxygenase in decapods. In this

68

research, the full-length cDNA sequence of one carotenoid isomerooxygenase gene (named

69

EcNinaB-X1) was cloned from E. carinicauda. Then, EcNinaB-X1 knock-out (EcNinaB-X1-KO)

70

prawns were produced via CRISPR/Cas9 technology and the change of phenotypes was analyzed.

71

Furtherly, the EcNinaB-X1-KO prawns were challenged with different pathogenic bacteria and the

72

function of EcNinaB-X1 in immune defense was also clarified.

73 74

2. Materials and methods

75

2.1. Ethical Statement

76

All efforts were made to minimize animal suffering. This article does not contain any studies

77

with human participants.

78 79

2.2 Cultivation of the experimental animals The ridgetail white prawns, E. carinicauda, were bred in plastic tanks filled with aerated fresh

80 81

seawater at 26 ºC, and fed twice per day in our laboratory.

82

Referring to our previous research [11], the one-cell stage embryos were collected from the

83

randomly selected spawning prawns and transferred to a clean petri dish loaded up with

84

appropriate filter-sterilized seawater before microinjection. After being injected, the experimental

85

embryos were put in the clean petri dishes containing sterilized seawater and incubated on a

86

shaking bath at 26 ºC with the speed of 100 rpm. Fresh clean seawater was changed per day. The

87

mysis larval prawns were fed with Artemia salina larvae and the juvenile prawns were fed with

88

bait.

89 90

2.3. RNA isolation, cDNA synthesis and bioinformatic analysis

91

Fifteen healthy prawns with a body length of 4.0 ± 0.5 cm were collected for tissue

92

distribution analysis. Hepatopancreas, muscle, eyestalk, gills, nerve, intestines, epidermis and

93

heart were dissected out and immediately preserved in liquid nitrogen for RNA extraction.

94

Total RNA was extracted from the collected samples with Trizol® reagent (Thermo, USA).

95

Then, the extracted RNA was treated with RQI RNase-Free DNase (Promega, USA). Two

96

micrograms of total RNA and 0.2 µM random hexamer primers were used to synthesize cDNA by

97

M-MLV reverse transcriptase (Promega, USA).

98

Based on the transcriptomic and genomic data of E. carinicauda [14], the full-length

99

NinaB-X1

sequence

of

E.

carinicauda

(EcNinaB-X1)

was

confirmed

by

reverse

100

transcription-polymerase chain reaction (RT-PCR). The nucleotide sequence and deduced amino

101

acid

102

(http://www.ncbi.nlm.nih.gov/BLAST/). The characteristic structure of deduced EcNinaB-X1 was

103

predicted by SMART program (http://smart.embl-heidelberg.de/). The multiple sequence

104

alignments and phylogenetic analysis were performed using CLUSTAL W and MEGA 7.0 [25].

sequence

of

EcNinaB-X1

were

analyzed

by

BLAST

on-line

105 106

2.4. EcNinaB-X1 expression in different tissues analyzed by quantitative real-time PCR

107

Epidermis, heart, gill, eyestalk, hepatopancreas, intestine, stomach, nerve and muscle were

108

separated from 15 healthy prawns with a body length of 4.0 ± 0.5 cm for tissue distribution.

109

Quantitative real-time PCR (qRT-PCR) [26] was used to analyze EcNinaB-X1 distribution in

110

different tissues of E. carinicauda using Mastercycler ep realplex (Eppendorf). 18S rRNA was

111

used as the internal control. Primers were shown in Table 1. The expected size of EcNinaB-X1 and

112

18S rRNA was 131 bp and 147 bp, respectively. Samples were run in triplicate on the PCR System.

113

The data were analyzed using the comparative CT method and then subjected to one-way ANOVA

114

using SPSS 19.0. The p values less than 0.05 were considered statistically significant.

115 116 117 118

2.5. Designation and synthesis of gRNA specialized for EcNinaB-X1 According to our previous research [11], we designed one sgRNA target site for EcNinaB-X1 with the online tool ZiFiT (http://zifit.partners.org/ZiFiT/ChoiceMenu.aspx).

119

The gRNA of EcNinaB-X1 was synthesized using the Thermo Scientific TranscriptAid T7

120

High Yield Transcription Kit (Thermo, USA). Then, it was purified by phenol chloroform

121

extraction. The gRNA concentration was assessed by Nanodrop 2000 (Thermo Fisher Scientific,

122

USA) and the quality was assessed by electrophoresis on 1% agarose gel. Then it was preserved at

123

-80 ºC in portions for microinjection.

124 125

2.6. Preparation of Cas9 mRNA

126

According to our previous research [11], the pCMV-Cas9 plasmid (Sigma-Aldrich, USA)

127

was linearized by Xba I (Takara, Dalian) and purified by ethanol precipitation. The linearized

128

product was used as the template to synthesize Cas9 mRNA that have both 5’cap and 3’poly (A)

129

tail in vitro with mMACHINE® T7 Ultra Kit (Ambion, USA). Then, it was purified by phenol

130

chloroform extraction and preserved at -80 ºC in portions for microinjection.

131 132

2.7. Microinjection and Indels detection by Sanger sequencing

133

The microinjection liquid contains 200 ng/µL Cas9 mRNA, 100 ng/µL gRNA and 0.05% of

134

the inert dye phenol red as the indicator in the buffer (100 mM HEPES, 1.5 M NaCl). Before

135

microinjection, the liquids were filtered through 0.22 µm filtering membranes. A Warner

136

PLI-100A Pico-Injector microinjector (Warner Instruments, USA) and a dissecting microscope

137

MN-152 micromanipulator (Narishige, Japan) were used for microinjection with standardized

138

Femtotip II sterile microcapillaries (Eppendorf, Germany). The injection volume was

139

approximately 0.5 nL. The prawns in experimental group were injected with Cas9 mRNA and

140

EcNinaB-X1 gRNA.

141

The genomic DNA of mysis larvae prawns was extracted and the genomic region flanking the

142

target site was amplified by MightyAmp® Genotyping Kit (Takara, Dalian) according to the

143

manufacturer’s instruction. For Sanger sequencing detection, the amplified PCR products were

144

purified using Gel Extraction Kit (Tiangen, China) and then cloned into pMD19-T Simple Vector

145

(Takara, Dalian). The detection primers used to amplify the target fragment are detEcNinaB-F and

146

detEcNinaB-R (Table 1).

147 148

2.8. Evaluation of CRISPR/Cas9 generated EcNinaB-X1-KO prawns

149

According to our previous research [17], the heterozygous EcNinaB-X1-KO progenies from

150

the same family with one wild-type allele and one allele harboring 4 bp deletions at exon 1 were

151

crossed to produce wild-type, heterozygous and homozygous EcNinaB-X1-KO prawns [11]. The

152

generated wild-type prawns were crossed to produce the filial generations (wild-type prawns) and

153

the generated homozygous EcNinaB-X1-KO prawns were also crossed to produce the filial

154

generations (EcNinaB-X1-KO prawns). After being hatched, the prawns were fed with Artemia

155

salina larvae. When the mysis larvae grew into juvenile prawns, they were fed with clam meat.

156

Difference of phenotype between the wild-type and EcNinaB-X1-KO prawns: the prawns and

157

dissected hepatopancreas were photographed and the color difference was calculated according to

158

the method described by reference [27].

159

Vibrio parahaemolyticus and Aeromonas hydrophila used in this study were isolated and

160

identified by Dr. Yuying Sun. The bacteria were routinely cultured in Tryptic Soy Broth (TSB,

161

Difco) or TSA medium supplemented with additional 1% NaCl at 28 ºC, 180 rpm.

162

The wild-type and EcNinaB-X1-KO prawns cultured with the same size were challenged with

163

V. parahaemolyticus and A. hydrophila according to previous research [17]. EcNinaB-X1-KO

164

group and wild-type group were set up for each sampling point and 200 prawns were sampled

165

from each group. For the bacterial challenge experiment, the experimental group was injected

166

individually with 10 µL phosphate buffer saline (PBS) containing V. parahaemolyticus or A.

167

hydrophila (107 CFU mL-1). At the same time, the prawns injected with 10 µL sterile PBS were

168

maintained as the control. The residual prawns were calculated at 0, 12, 24, 48, 72, 96, and 120 h.

169 170

3. Results

171

3.1. Amplification and characterization of EcNinaB-X1

172

Based on the transcriptomic and genomic data of E. carinicauda published by Yuan et al. [14],

173

the full-length cDNA sequence of EcNinaB-X1 was obtained with 3185 bp (GenBank accession no.

174

MN583041). As shown in Fig.1, the nucleotide sequence of EcNinaB-X1 contained a 2070 bp

175

open reading frame (ORF) encoding EcNinaB-X1 of 689 amino acids with a predicted molecular

176

weight (MW) about 78286.67 Da and theoretical isoelectric point (pI) of 7.37. No putative signal

177

peptide was found. The domain architecture of deduced EcNinaB-X1 predicted by SMART

178

software showed that there was a RPE65 domain (residues 152-670) (Fig. 2). In addition, the

179

genomic DNA fragment of EcNinaB-X1 with the corresponding cDNA sequence was obtained,

180

which showed that it was composed of 12 exons and 11 introns (Fig.3). All intron-exon boundaries

181

are consistent with the consensus splicing junctions at both the 5’ splice donor site (GT) and the 3’

182

splice acceptor site (AG) of each intron.

183

A multiple sequence alignment showed that EcNinaB-X1 displayed high identities with

184

carotenoid isomerooxygenase of Armadillidium vulgare (AvNinaB-X1, 60%), and Nilaparvata

185

lugens (NlNinaB-X1, 51%) (Fig. 4).

186

Amino acid sequences of carotenoid isomerooxygenase (NinaB), β, β-carotene 15,

187

15'-monooxygenase (BCMO) and β, β-carotene 9', 10'-oxygenase (BCO2) from different species

188

were collected from the NCBI database and a phylogenetic tree was constructed using

189

Neighbor-joining method. The phylogenetic analysis showed that EcNinaB-X1 was divided into

190

the same branch of NinaB (Fig.5).

191 192

3.2 Tissue distribution of EcNinaB-X1 mRNA

193

Expression profiling of EcNinaB-X1 mRNA in different tissues of E. carinicauda was

194

detected by quantitative real-time PCR. The expression of EcNinaB-X1 was expressed

195

constitutively in all of the detected tissues and highly expressed in the hepatopancreas, gill, and

196

stomach (Fig. 6).

197 198

3.3. Knockout of EcNinaB-X1 using CRISPR/Cas9 technology

199

Ten days after injection to the one-cell embryo, ten embryos were selected randomly in each

200

group and extracted the genomic DNA. Then the target fragment was amplified with detection

201

primers. The sequencing results of the purified PCR products showed that multiple peaks occurred

202

initially after the PAM sites compared with blank and control groups (Fig. 7). It indicated that the

203

genome of prawns in the experimental group had been successfully edited. In order to identify the

204

result, cloning and Sanger sequencing were carried out. Results showed that seven types of

205

deletion mutation were generated in total, but no insertions were found (Fig. 8).

206

The mutation rate of the embryos injected with Cas9 mRNA and one sgRNA was analyzed

207

and compared with the results of our previous research (Table 2). Of the 400 injected one-cell

208

stage embryos with Cas9 mRNA and one sgRNA targeting the first exon of EcNinaB-X1 gene, 58

209

embryos could develop to postlarvae and the reproductive survival rate was 14.5%. Fifty days

210

later, when the 58 lived embryos developed to adult E. carinicauda, one leg of them was collected

211

respectively to detect indels and the results indicated that the number of mutant prawns was 26

212

and the mutant rate reached 6.5%. Furtherly, of the 26 mutant prawns, 5 kinds of mutant types

213

were identified in the targeted locus (Fig. 9).

214 215

3.4. Evaluation of CRISPR/Cas9 generated EcNinaB-X1-KO prawns

216

To determine the effect of the EcNinaB-X1-KO on the prawns, the phenotypes and dissected

217

hepatopancreas of wild-type and EcNinaB-X1-KO prawns were compared. The hepatopancreas

218

from EcNinaB-X1-KO prawns exhibited a red color compared with wild-type individuals (Fig.

219

10A, B). According to the reference [27], the industrial color criteria was used to evaluated the

220

color change in the hepatopancreas between wild-type and EcNinaB-X1-KO prawns, and the

221

EcNinaB-X1-KO prawns displayed a significant shift to red color compared to the wild-type

222

groups (Fig. 10C).

223

The effects of EcNinaB-X1-KO on the mortality of bacteria-challenged prawns by

224

challenging them with were also evaluated. The prawns in EcNinaB-X1-KO group had

225

significantly low mortality than those in wild-type group when the prawns were challenged with V.

226

parahaemolyticus or A. hydrophila (Fig. 11A, B). For the Vibrio-challenged prawns, the mortality

227

reached 50% in wild-type group and only 20% in EcNinaB-X1-KO group at 48 h post-challenge

228

and there was significant difference between wild-type and EcNinaB-X1-KO prawns (p < 0.01).

229

For the Aeromonas-challenged prawns, the mortality reached 65% in wild-type group and only 25%

230

in EcNinaB-X1-KO group at 48 h post-challenge and there was significant difference between

231

wild-type and EcNinaB-X1-KO prawns (p < 0.05). The negative control showed a cumulative

232

mortality of ~10% indicating that PBS-injection itself was non-toxic in prawns. Overall, these

233

results indicated that the deletion of EcNinaB-X1 in E. carinicauda decreased prawn’s mortality

234

following V. parahaemolyticus or A. hydrophila challenge.

235 236

4. Discussion

237

Carotenoids are naturally occurring red, orange and yellow pigments that are synthesized by

238

plants, alga and some microorganisms and fulfill many important physiological functions [28].

239

Carotenoids are fat-soluble compounds, and it is important to enhance the bioavailability of

240

β-carotene from food or supplements in animals [29]. β-carotene is involved in a number of

241

beneficial functions in animals such as taking part in crucial signaling functions of its metabolites

242

and working as an antioxidant [27]. In addition, the accumulation of carotenoids in animals could

243

change the fat color in some animals, which make it important in animal breeding [27, 30]. Li et al

244

[31] cultured a new variety of Yesso scallop (Patinopecten yessoensis) with orange adductor

245

muscle, the ‘Haida golden scallop’, which was caused by carotenoid accumulation.

246

In animals, carotenoids could be cleaved into a range of biologically important products and

247

carotenoid oxygenases play an important role in the reaction [32]. It is reported that there are some

248

carotenoid oxygenases, including carotenoid isomerooxygenase (NinaB), β, β-carotene 15,

249

15'-monooxygenase (BCMO) and β, β-carotene 9', 10'-oxygenase (BCO2) [20, 21]. In this

250

research, the full-length cDNA and genomic DNA of one carotenoid isomerooxygenase gene

251

(EcNinaB-X1) were obtained from E. carinicauda. There is one typical RPE65 domain at the

252

position of 152-670 in the deduced EcNinaB-X1 amino acid sequence. From the bioinformatic

253

analysis, we could find that the RPE65 domain played the oxidoreductase activity acting on single

254

donors with incorporation of molecular oxygen, or two atoms of oxygen [29]. Among three typical

255

carotenoid oxygenases, the phylogenetic analysis showed that EcNinaB-X1 was divided into the

256

same branch of NinaB, not belonging to β, β-carotene 15, 15'-monooxygenase (BCMO) and β,

257

β-carotene 9', 10'-oxygenase (BCO2).

258

Niu et al [27] reported that the biallelic modification of BCO2 from sheep resulted in yellow

259

fat, compared with the fat color in wild types (snow-flower white). Vage & Boman [33] reported

260

that a nonsense mutation in the BCO2 gene was tightly associated with accumulation of

261

carotenoids in adipose tissue of sheep (Ovis aries). Strychalski et al [34] reported that an

262

AAT-deletion mutation in the coding sequence of the BCO2 gene was related to the formation of

263

yellow-fat rabbit. Kyle-Little et al [35] reported that the mutant of BCMO1 gene affected macular

264

pigment optical density in young healthy caucasians. Amengual et al [36] reported that the

265

BCMO1 acted as a critical molecular player in deducing body adiposity of mice. However, there

266

was no report about the biological function of NinaB in regulating the physiology and

267

biochemistry of animals. In this research, we used CRISPR-Cas9 system to disrupt EcNinaB-X1

268

gene and obtained wild-type and homozygous mutants at EcNinaB-X1 loci in the offspring. In this

269

study, wild-type and homozygous EcNinaB-X1-KO prawns were selected as the experimental

270

animals to study the function of EcNinaB-X1 in immune defense. EcNinaB-X1-KO prawns

271

showed much lower mortality than those in wild-type group after V. parahaemolyticus or A.

272

hydrophila challenge, which indicated that EcNinaB-X1 might play a key role in immune defense

273

of prawns. In addition, the deletion of EcNinaB-X1 gene resulted in the color shift compared with

274

wild-type prawns. From above research, we concluded that NinaB-X1 gene could be used as a

275

candidate gene in marker-assistant selection in animal breeding.

276 277 278

Conflict of interest There is no conflict of interest.

279 280

Acknowledgments

281

This work was supported by the National Key R&D Program of China (No.

282

2018YFD0900205), National Natural Science Foundation of China (Grant Nos. 31872613,

283

41876196), The Natural Science Foundation of Hebei Province of China (Grant Nos.

284

C2019201236, D2019201239), Innovation & entrepreneurship training program for college

285

students of Hebei Province (Grant No. S201910075061), Innovation & entrepreneurship training

286

program for college students of Hebei University (Grant No. 2019177), and Hangzhou Qianjiang

287

Special Expert for Jiquan Zhang. We are grateful to Dr. Huan Gao for providing the prawns.

288 289

References

290 291 292 293 294 295 296 297 298 299 300 301 302 303

[1] L. Cong, F. Ran, D. Cox, S. Lin, R. Barretto, N. Habib, P. Hsu, X. Wu, W. Jiang, L. Marraffini, F. Zhang, Multiplex genome engineering using CRISPR/Cas systems, Science 339(6121) (2013) 819-23. [2] P. Mali, L. Yang, K. Esvelt, J. Aach, M. Guell, J. DiCarlo, J. Norville, G. Church, RNA-guided human genome engineering via Cas9, Science 339(6121) (2013) 823-6. [3] D. Li, Z. Qiu, Y. Shao, Y. Chen, Y. Guan, M. Liu, Y. Li, N. Gao, L. Wang, X. Lu, Y. Zhao, M. Liu, Heritable gene targeting in the mouse and rat using a CRISPR-Cas system, Nat. Biotechnol. 31(8) (2013) 681-3. [4] D. Harms, R. Quadros, D. Seruggia, M. Ohtsuka, G. Takahashi, L. Montoliu, C. Gurumurthy, Mouse genome editing using the CRISPR/Cas system, Curr Protoc Hum Genet 83 (2014) 15 7 1-27. [5] A. Bassett, C. Tibbit, C. Ponting, J. Liu, Highly efficient targeted mutagenesis of Drosophila with the CRISPR/Cas9 system, Cell Rep 4(1) (2013) 220-8. [6] S. Gratz, A. Cummings, J. Nguyen, D. Hamm, L. Donohue, M. Harrison, J. Wildonger, K. O'Connor-Giles, Genome engineering of Drosophila with the CRISPR RNA-guided Cas9 nuclease, Genetics 194(4) (2013) 1029-35.

304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347

[7] W. Hwang, Y. Fu, D. Reyon, M. Maeder, S. Tsai, J. Sander, R. Peterson, J. Yeh, J. Joung, Efficient genome editing in zebrafish using a CRISPR-Cas system, Nat. Biotechnol. 31(3) (2013) 227-9. [8] L. Jao, S. Wente, W. Chen, Efficient multiplex biallelic zebrafish genome editing using a CRISPR nuclease system, Proc Natl Acad Sci U S A 110(34) (2013) 13904-9. [9] J. Zhang, F. Song, Y. Sun, K. Yu, J. Xiang, CRISPR/Cas9-mediated deletion of EcMIH shortens metamorphosis time from mysis larva to postlarva of Exopalaemon carinicauda, Fish Shellfish Immunol 77 (2018) 244-51. [10] J. Zhang, Y. Liu, Y. Li, N. Su, Y. Zhou, J. Xiang, Y. Sun, Biological function of a gC1qR homolog (EcgC1qR) of Exopalaemon carinicauda in defending bacteria challenge, Fish Shellfish Immunol 82 (2018) 378-85. [11] T. Gui, J. Zhang, F. Song, Y. Sun, S. Xie, K. Yu, J. Xiang, CRISPR/Cas9-mediated genome editing and mutagenesis of EcChi4 in Exopalaemon carinicauda, G3-Genes Genom Genet 6(11) (2016) 3757-64. [12] P. Hsu, E. Lander, F. Zhang, Development and applications of CRISPR-Cas9 for genome engineering, Cell 157(6) (2014) 1262-78. [13] X. Zhang, J. Yuan, Y. Sun, S. Li, Y. Gao, Y. Yu, C. Liu, Q. Wang, X. Lv, X. Zhang, K. Ma, X. Wang, W. Lin, L. Wang, X. Zhu, C. Zhang, J. Zhang, S. Jin, K. Yu, J. Kong, P. Xu, J. Chen, H. Zhang, P. Sorgeloos, A. Sagi, A. Alcivar-Warren, Z. Liu, L. Wang, J. Ruan, K. Chu, B. Liu, F. Li, J. Xiang, Penaeid shrimp genome provides insights into benthic adaptation and frequent molting, Nature Communications 10 (2019) 356. [14] J. Yuan, Y. Gao, X. Zhang, J. Wei, C. Liu, F. Li, J. Xiang, Genome sequences of marine shrimp Exopalaemon carinicauda Holthuis provide insights into genome size evolution of Caridea, Mar. Drugs 15(7) (2017) 213. [15] N. Kenny, Y. Sin, X. Shen, Q. Zhe, W. Wang, T. Chan, S. Tobe, S. Shimeld, K. Chu, J. Hui, Genomic sequence and experimental tractability of a new decapod shrimp model, Neocaridina denticulata, Mar. Drugs 12(3) (2014) 1419-37. [16] L. Song, C. Bian, Y. Luo, L. Wang, X. You, J. Li, Y. Qiu, X. Ma, Z. Zhu, L. Ma, Z. Wang, Y. Lei, J. Qiang, H. Li, J. Yu, A. Wong, J. Xu, Q. Shi, P. Xu, Draft genome of the Chinese mitten crab, Eriocheir sinensis, Gigascience 5 (2016) 5. [17] Y. Sun, J. Zhang, J. Xiang, A CRISPR/Cas9-mediated mutation in chitinase changes immune response to bacteria in Exopalaemon carinicauda, Fish Shellfish Immunol 71 (2017) 43-9. [18] F. Karadas, E. Grammenidis, P. Surai, T. Acamovic, N. Sparks, Effects of carotenoids from lucerne, marigold and tomato on egg yolk pigmentation and carotenoid composition, Br Poult Sci 47(5) (2006) 561-6. [19] Y. Zhang, L. Zhang, J. Sun, J. Qiu, X. Hu, J. Hu, Z. Bao, Proteomic analysis identifies proteins related to carotenoid accumulation in Yesso scallop (Patinopecten yessoensis), Food Chem. 147 (2014) 111-6. [20] V. Oberhauser, O. Voolstra, A. Bangert, J. von Lintig, K. Vogt, NinaB combines carotenoid oxygenase and retinoid isomerase activity in a single polypeptide, Proc Natl Acad Sci U S A 105(48) (2008) 19000-5. [21] J. Amengual, M.A. Widjaja-Adhi, S. Rodriguez-Santiago, S. Hessel, M. Golczak, K. Palczewski, J. von Lintig, Two carotenoid oxygenases contribute to mammalian provitamin A metabolism, J Biol Chem 288(47) (2013) 34081-96. [22] A. Martinez, A. Melendez-Martinez, Lycopene, oxidative cleavage derivatives and antiradical

348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 381 382

activity, Computational and Theoretical Chemistry 1077 (2016) 92-8. [23] T. Papa, V. Pinho, E. do Nascimento, W. Santos, A. Burtoloso, L. Skibsted, D. Cardoso, Astaxanthin diferulate as a bifunctional antioxidant, Free Radical Res. 49(1) (2015) 102-11. [24] E. Rodrigues, L. Mariutti, A. Mercadante, Scavenging capacity of marine carotenoids against reactive oxygen and nitrogen species in a membrane-mimicking system, Mar. Drugs 10(8) (2012) 1784-98. [25] S. Kumar, G. Stecher, K. Tamura, MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets, Mol Biol Evol 33(7) (2016) 1870-74. [26] C. Yang, J. Zhang, F. Li, H. Ma, Q. Zhang, T. Jose Priya, X. Zhang, J. Xiang, A Toll receptor from Chinese shrimp Fenneropenaeus chinensis is responsive to Vibrio anguillarum infection, Fish Shellfish Immunol 24(5) (2008) 564-74. [27] Y. Niu, M. Jin, Y. Li, P. Li, J. Zhou, X. Wang, B. Petersen, X. Huang, Q. Kou, Y. Chen, Biallelic beta-carotene oxygenase 2 knockout results in yellow fat in sheep via CRISPR/Cas9, Anim Genet 48(2) (2016) 242-44. [28] J. Alcaino, M. Baeza, V. Cifuentes, Carotenoid Distribution in Nature, Subcell. Biochem. 79 (2016) 3-33. [29] J.D. Ribaya-Mercado, Influence of dietary fat on beta-carotene absorption and bioconversion into vitamin A, Nutr. Rev. 60(4) (2002) 104-10. [30] K.M. Sefc, A.C. Brown, E.D. Clotfelter, Carotenoid-based coloration in cichlid fishes, Comp. Biochem. Physiol. A Mol. Integr. Physiol. 173C (2014) 42-51. [31] X. Li, X. Ning, J. Dou, Q. Yu, S. Wang, L. Zhang, S. Wang, X. Hu, Z. Bao, An SCD gene from the Mollusca and its upregulation in carotenoid-enriched scallops, Gene 564(1) (2015) 101-8. [32] G. Giuliano, S. Al-Babili, J. von Lintig, Carotenoid oxygenases: cleave it or leave it, Trends Plant Sci. 8(4) (2003) 145-9. [33] D.I. Vage, I.A. Boman, A nonsense mutation in the beta-carotene oxygenase 2 (BCO2) gene is tightly associated with accumulation of carotenoids in adipose tissue in sheep (Ovis aries), BMC Genet. 11 (2010) 10. [34] J. Strychalski, P. Brym, U. Czarnik, A. Gugolek, A novel AAT-deletion mutation in the coding sequence of the BCO2 gene in yellow-fat rabbits, J Appl Genet 56(4) (2015) 535-37. [35] Z. Kyle-Little, A.J. Zele, C.P. Morris, B. Feigl, The effect of BCMO1 gene variants on macular pigment optical density in young healthy caucasians, Front Nutr 1 (2014) 22. [36] J. Amengual, E. Gouranton, Y.G. van Helden, S. Hessel, J. Ribot, E. Kramer, B. Kiec-Wilk, U. Razny, G. Lietz, A. Wyss, A. Dembinska-Kiec, A. Palou, J. Keijer, J.F. Landrier, M.L. Bonet, J. von Lintig, Beta-carotene reduces body adiposity of mice via BCMO1, PLoS ONE 6(6) (2011) e20644.

Table 1 Primers mentioned in the paper

383

Primers

384 385

Sequences (5’- 3’)

Sequence information

EcNinaF

ACAGCCCCTTGTCGTGTCCGT

Real-time PCR

EcNinaR

TCGTCGCTCCCGTGTCCCT

Real-time PCR

18S-F

TATACGCTAGTGGAGCTGGAA

Real-time PCR

18S-R

GGGGAGGTAGTGACGAAAAAT

Real-time PCR

DetEcNinaF

TTTCAGTAACTCTTACGATTC

Detection of mutation

DetEcNinaR

CCCGCAGTGGTAGGAATAGG

Detection of mutation

Note: F and R stand for forward primers and reverse ones, respectively.

Table 2 Mutation frequencies induced by microinjection of Cas9 mRNA and EcNinaB-X1 gRNA

386

Injected Group

Survival

Mutant

Survival

Mutant

RNA concentration

100 ng/µL

Embryos

Postlarvae

Postlarvae

Rate

Rate

247

35

18

14.17%

7.29%

[11]

250

23

12

9.20%

4.80%

[9]

400

58

26

14.50%

6.50%

In this paper

200 ng/µL

EcChi4 (gRNA-EcChi4)

(pCMV-Cas9)

100 ng/µL

200 ng/µL

EcMIH (gRNA-EcMIH)

(pCMV-Cas9)

100 ng/µL

200 ng/µL

EcNinaB-X1 (gRNA-EcNinaB-X1)

References

(pCMV-Cas9)

387

Legends of Figures

388

Fig.1 Nucleotide and deduced amino acid sequences of EcNinaB-X1 gene.

389

Nucleotides are numbered on the both sides of the sequence. The letters marked with double

390

underline represented the RPE65 domains.

391

Fig.2 The domain architecture of deduced EcNinaB-X1.

392

Fig. 3 Schematic representation of the structure of EcNinaB-X1 gene depicting exons, gRNA

393

target sequence (sequence in blue at the bottom) and PAM site (in red). The numbers indicate the

394

exact lengths of the exons.

395

Fig. 4 Alignment of the amino acid sequence of EcNinaB-X1 with other known NinaBs.

396

The identical residues are shown in solid boxes. Sequences start at the first methionine residue.

397

Armadillidium vulgare (AvNinaB-X1, GenBank accession no. RXG54508.1); Nilaparvata lugens

398

(NlNinaB-X1, XP_022185479.1); E. carinicauda (EcNinaB-X1, MN583041, in this research).

399

Fig.5 Phylogenetic tree of carotenoid isomerooxygenase (NinaB), β, β-carotene 15,

400

15'-monooxygenase (BCMO) and β, β-carotene 9', 10'-oxygenase (BCO2) from different species

401

based on the amino acid sequence comparisons.

402

Aedes albopictus (AaNinaB-X1, GenBank accession no. XP_019545564.1); Aethina tumida

403

(AtNinaB-X1, XP_019870879.1); Anoplophora glabripennis (AgNinaB-X1, XP_023310515.1);

404

Apis cerana (AcNinaB-X1, XP_016913671.1); A. dorsata (AdNinaB-X1, XP_006611100.1); A.

405

florea

406

Armadillidium vulgare (AvNinaB-X1, GenBank accession no. RXG54508.1); Bombus terrestris

407

(BtNinaB-X1, XP_003400564.1); Branchiostoma belcheri (BbBCO2-X1, XP_019636508.1);

408

Callorhinchus milii (CmBCMO-X1, XP_007887463.1); Ceratina calcarata (CcNinaB-X1,

(AfNinaB-X1,

XP_003695005.1);

A.

mellifera

(AmNinaB-X1,

XP_394000.4);

409

XP_017875448.1); Clupea harengus (ChBCO2-X1, XP_012685829.1); Dendroctonus ponderosae

410

(DpNinaB-X1, XP_019773593.1); Diachasma alloeum (DaNinaB-X1, XP_015113077.1);

411

Dinoponera quadriceps (DqNinaB-X1, XP_014471901.1); Dufourea novaeangliae (DnNinaB-X1,

412

XP_015439200.1); Eufriesea mexicana (EmNinaB-X1, XP_017759905.1); Fopius arisanus

413

(FaNinaB-X1, XP_011310428.1); Fundulus heteroclitus (FhBCMO-X1, XP_012709399.1);

414

Habropoda laboriosa (HlNinaB-X1, XP_017798623.1); Harpegnathos saltator (HsNinaB-X1,

415

XP_011149992.1); Ictalurus punctatus (IpBCO2-X1, XP_017316112.1); Latimeria chalumnae

416

(LcBCMO-X1,

417

XP_014345540.1;

418

XP_013396242.1); Megachile rotundata (MrNinaB-X1, XP_003702796.1); Nilaparvata lugens

419

(NlNinaB-X1, XP_022185479.1); Pseudomyrmex gracilis (PgNinaB-X1, XP_020278007.1);

420

Salmo salar (SsBCMO-X1, XP_014030705.1; SsBCMO-X2, XP_014030708.1; SsBCMO-X3,

421

XP_014030708.1); Sarcophilus harrisii (ShBCMO-X1, XP_003757903.1); Trichogramma

422

pretiosum (TpNinaB-X1, XP_014224940.1); E. carinicauda (EcNinaB-X1, MN583041, in this

423

research). Values on the line are bootstrap values showing percentage confidence of relatedness.

424

Fig.6 Detection of EcNinaB-X1 transcripts in different tissues of E. carinicauda. Tissues were

425

shown in the abscissa. The amount of EcNinaB-X1 mRNA was normalized to the 18S rRNA

426

transcript level. Data are shown as means ± SD (standard deviation) of three separate individuals

427

in the tissues.

428

Fig.7 Sanger sequencing of the PCR products from 10 prawns indicated the indel mutations

429

caused by CRISPR/Cas9 genome editing system. WT means the wild-type prawns (blank group).

430

M1 and M2 mean injected embryos. The PAM site was represented in pink rectangles.

XP_006004911.1; LcBCO2-X3,

LcBCO2-X1, XP_014345541.1);

XP_014345536.1; Lingula

anatina

LcBCO2-X2, (LaBCO2-X1,

431

Fig.8 Sequencing of CRISPR/Cas9 treated prawns. The wild-type sequence (ref) is shown at the

432

top. The gRNA site is underlined and the sequence in pink frame (TGG) represents the PAM

433

sequence. The numbers on the right represents the number of mutations recovered by sequencing.

434

Within the sequences, deletions are indicated by dashed lines.

435

Fig. 9 Sanger sequencing of the PCR products from prawns individually indicated the indel

436

mutations caused by CRISPR/Cas9 genome editing system. WT means the wild-type prawn (blank

437

group). MT-1, 2, 3, 4, and 5 mean five mutant-type prawns. The PAM site was represented in

438

black rectangles.

439

Figure 10 Phenotypes and dissection of EcNinaB-X1 knockouts in the hepatopancreas of E.

440

carinicauda. (a) Phenotypes of wild-type and EcNinaB-X1-KO prawns. (b) Hepatopancreas

441

dissected from prawns showed in figure (a). (c) Color differences in the hepatopancreas of

442

wild-type and EcNinaB-X1-KO prawns; ‘R’ indicates red color, ‘G’ indicates green color, ‘B’

443

indicates blue color. According to the reference [27], we defined the (R+G)/B value to indicate the

444

red color.

445

Fig. 11 The mortality of the prawns in wild-type group and EcNinaB-X1-KO group after the

446

prawns were challenged with Vibrio parahaemolyticus and equal volume of PBS (A), or

447

Aeromonas hydrophila and equal volume of PBS (B) at 0, 12, 24, 48, 72, 96, and 120 h.

448

449

450 Fig.1

451 452 453

.

454 455 456

Fig.2

457 458 459 460

Fig. 3

461 462 463

Fig. 4

464 465 466

Fig.5

467 468 469

Fig.6

470

471 472 473

Fig.7

474 475 476

Fig. 8

477

478 479 480

Fig. 9

481

482

(a)

483

(b)

(c)

484 485

Figure 10

486

.

487

488 489

Fig. 11

Highlights

EcNinaB-X1 gene was cloned from E. carinicauda. The expression profiles of EcNinaB-X1 were demonstrated. CRISPR-Cas9 system efficiently generated indels in EcNinaB-X1 loci. EcNinaB-X1-KO prawns had lower mortality than wild-type after bacterial challenge.