Expression and functional analysis of two odorant binding proteins

Mar 14, 2019 - Expression and functional analysis of two odorant binding proteins from Bradysia odoriphaga (Diptera: Sciaridae). Bowen Tang , Shulei T...
0 downloads 0 Views 2MB Size
Subscriber access provided by La Trobe University Library

Agricultural and Environmental Chemistry

Expression and functional analysis of two odorant binding proteins from Bradysia odoriphaga (Diptera: Sciaridae) Bowen Tang, Shulei Tai, Wu Dai, and Chunni Zhang J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.9b00568 • Publication Date (Web): 14 Mar 2019 Downloaded from http://pubs.acs.org on March 17, 2019

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

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 48

Journal of Agricultural and Food Chemistry

1

Expression and functional analysis of two odorant binding proteins from

2

Bradysia odoriphaga (Diptera: Sciaridae)

3

Bowen Tang, Shulei Tai, Wu Dai*, Chunni Zhang*

4

Key Laboratory of Plant Protection Resources and Pest Integrated Management of the

5

Ministry of Education, College of Plant Protection, Northwest A&F University,

6

Yangling, Shaanxi, China

7 8

*Corresponding

9

Wu Dai, Tel. 86-13389220985; email: [email protected];

10

Chunni

author.

Zhang,

Tel.

86-18729551958;

email:

1

ACS Paragon Plus Environment

[email protected]

Journal of Agricultural and Food Chemistry

11

ABSTRACT

12

Two OBP genes, BodoOBP1 and BodoOBP2 were cloned from Bradysia

13

odoriphaga, a major agricultural pest of Chinese chives. The amino acid sequence

14

alignment of both BodoOBPs showed high similarity. Fluorescence competitive

15

binding assays revealed that both BodoOBPs have a moderate binding affinity to

16

dipropyl trisulfide. Tissue expression profiles indicated that both BodoOBPs are

17

antennae-specific and more abundant in the male antennae than in the female

18

antennae. Developmental expression profile analysis indicated that expression levels

19

of both BodoOBPs were higher in the male adult stage than in the other

20

developmental stages. Both BodoOBPs also showed differential expression in pre-

21

and post-mating adults. RNAi assays indicated that ability of dsOBPs-males to detect

22

females was significantly reduced compared to controls. Attraction of plant volatile

23

dipropyl trisulfide to dsOBPs-treated adults was also significantly lower than in the

24

control. Our findings indicate that both BodoOBPs are involved in host-seeking

25

behavior and in detecting sex pheromones.

26 27

KEYWORDS: Bradysia odoriphaga, odorant binding protein, binding affinity,

28

expression profile, RNA interference

2

ACS Paragon Plus Environment

Page 2 of 48

Page 3 of 48

Journal of Agricultural and Food Chemistry

29

INTRODUCTION

30

The olfactory system is essential for multiple behaviors of insects, including

31

locating sexual partners and oviposition sites, host-seeking, and avoiding predators1-4.

32

The insect olfactory system is mediated by eight main types of proteins: soluble

33

binding proteins [odorant binding proteins (OBPs), chemosensory proteins (CSPs)

34

and Niemann-Pick C2 protein (NPC2)]5, 6, membrane-bound receptors [odorant

35

receptors (ORs), gustatory receptors (GRs), ionotropic receptors (IRs) and sensory

36

neuron membrane proteins (SNMPs)] and odorant-degrading enzymes (ODEs)4.

37

OBPs are small soluble proteins that are recognized as playing a central role in the

38

initial phases of olfactory perception. OBPs are highly concentrated in the sensillum

39

lymph of chemosensilla and are capable of binding and transporting hydrophobic

40

odorant molecules via the chemosensilla lymph to corresponding chemosensory

41

receptors1, 4, 7, 8.

42

To date, a large number of OBP sequences have been identified from different

43

insect orders including Hemiptera, Lepidoptera, Diptera, Hymenoptera, Coleoptera,

44

Orthoptera, Isoptera, Anoplura, Blattaria and Neuroptera1, 9-15. Based on the number

45

of cysteine residues they contain, OBPs are now classified as “Classic”, “Dimer”,

46

“Minus - C”, “Plus - C”, and “Atypical” OBPs15. Most insect OBPs are characterized

47

by six conserved cysteines (Classic OBPs) that form three disulphide bridges and play

48

an important role in stabilizing the secondary and tertiary structure16, 17. The number

49

of OBPs is highly variable even between closely related species. For example, there

50

are 33 OBPs in Lygus lineolaris (Hemiptera: Miridae)18, 28 OBPs in Nysius ericae 3

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 4 of 48

51

(Hemiptera: Lygaeidae)19, 49 OBPs in Episyrphus balteatus (Diptera: Syrphidae), 44

52

OBPs in Eupeodes corollae (Diptera: Syrphidae)20. It has been demonstrated that

53

many of the identified OBP genes are highly expressed in the antennae11,

54

OBPs were also reportedly expressed in other chemosensory organs, including

55

mouthpart palps, maxillae, labium, legs, and wings, suggesting that some OBPs might

56

also be involved in recognition of taste compounds18,

57

studies have demonstrated that several OBPs are specifically distributed in

58

non-chemosensory organs, such as the male reproductive organs35, 36. For example, in

59

Bactrocera dorsalis, OBP19c was abundantly expressed in the thorax and abdomen37.

60

All of these studies revealed that OBPs play significant roles in insect olfactory

61

systems, and may have different functions in non-sensory organs of the insect body,

62

such as pheromone delivery, solubilisation of nutrients, development and insecticide

63

resistance8.

26-34.

18, 21-25.

However, some recent

64

The Chinese chive maggot, Bradysia odoriphaga Yang et Zhang (Diptera:

65

Sciaridae) is a major agricultural pest that causes severe economic damage in

66

liliaceous vegetable crops, edible mushrooms and humus38,

67

(Allium tuberosum) suffers the most severe damage in Northern China, with yield loss

68

of 30%–80% and economic losses of over 30%38, 40. B. odoriphaga larvae tend to

69

aggregate in the roots, bulbs, and immature stems of Chinese chives, resulting in

70

moisture loss and even death41. Since the larvae primarily damage the underground

71

portions of plants, infestation is difficult to prevent or control. Management practices

72

against B. odoriphaga still largely rely on the application of chemical insecticides. 4

ACS Paragon Plus Environment

39.

The Chinese chive

Page 5 of 48

Journal of Agricultural and Food Chemistry

73

However, insecticide use causes environmental pollution, may affect human health

74

and may also lead to pesticide resistance. Because B. ororiphaga presumably locate

75

their host plants by detecting volatile compounds released by the plants and also use

76

pheromones in courtship42-44, identification and functional characterization of OBP

77

and other olfactory genes are of potential importance for developing olfaction-based

78

and environmentally friendly techniques to control the pest. At present, 49 OBP genes

79

have been identified from B. odoriphaga transcriptomic data14, but their functions

80

have not been addressed.

81

In this study, in order to explore the possible olfactory roles of the two B.

82

odoriphaga OBPs, quantitative real-time PCR was used to assess the expression

83

profiles of the two OBPs genes in different tissues of both sexes, at different

84

developmental stages, and before and after mating. In addition, the binding abilities of

85

the two OBPs to host plant volatiles were measured using fluorescence competitive

86

binding assays and the roles of two BodoOBPs in the detection of odorants were

87

investigated by RNAi. Our systematic studies provide a better understanding of the

88

molecular basis of B. odoriphaga olfaction.

89

MATERIALS AND METHODS

90

Insect rearing and sample collection. A laboratory colony of B. odoriphaga was

91

established using specimens collected from a Chinese chive field in Shangcun Town,

92

Zhouzhi County, Xi’an City, China in 2014. The insects were reared on fresh

93

rhizomes of Chinese chive (1 cm in length) in petri dishes at 25 ± 1 °C, 70±5 %

94

relative humidity with a photoperiod of 16 : 8 h (L : D). For gene expression pattern 5

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

95

study, 1st to 4th instar larvae, pupae and adults were collected into RNase-free tubes

96

and immediately immersed in liquid nitrogen. The antennae, head, thorax, abdomen,

97

legs and wings were dissected under a microscope while soaking in Ringer’s Solution

98

(6.5g/L NaCl, 0.2g/L NaHCO3, 0.14g/L KCl, 0.01g/L NaH2PO4, 0.12g/L CaCl2) and

99

immediately transferred into RNase-free tubes and immersed into liquid nitrogen. The

100

individual virgin adults were reared in a separate plastic box after emergence. After

101

24h, 5 adults (male and female separately) were collected into RNase-free tubes;

102

meanwhile other mature virgin males and females were placed into one plastic box.

103

The mating behavioral assays were conducted at 25°C between 8:00 AM and 11:00

104

AM. 5 post-mating adults (male and female separately) were collected into

105

RNase-free tubes at 0h and 1h after the first successful mating. All samples were

106

collected in three biological replicates and stored at -80 °C until further experiments.

107

Total RNA isolation and cDNA synthesis. Total RNA for each sample was isolated

108

by using the RNAiso plus Reagent (TaKaRa, Japan). The quality and concentration of

109

total RNA were determined by agarose gel electrophoresis and spectrophotometer

110

analysis (Infinite M200 PRO, Tecan Group Ltd., Mannedorf, Switzerland). For gene

111

cloning, total RNA of adults was digested with DNase I to remove genomic DNA,

112

and then was used to synthesize first-strand cDNA using a RevertAid First Strand

113

cDNA Synthesis Kit (Thermo Scientific, USA). For quantitative real-time PCR

114

analysis, 1μg of each RNA sample was employed as template to synthesize

115

first-strand cDNA by PrimeScript™ RT reagent Kit with gDNA Eraser, Perfect Real

6

ACS Paragon Plus Environment

Page 6 of 48

Page 7 of 48

Journal of Agricultural and Food Chemistry

116

Time (TaKaRa, Japan). All operations strictly followed the manufacturer’s handbook.

117

The cDNA template was stored at -20 °C until use.

118

Cloning of the full-length OBPs cDNA. The partial sequence of two candidate OBPs

119

was identified based on the previously obtained transcriptome sequence. To obtain the

120

full‐length cDNA, 5’ and 3’ RACE was performed. Specific primers for RACE were

121

designed based on identified PCR fragment sequences (Table S1). The first strand

122

cDNA for 5’ and 3’ RACE was synthesized from 1μg of total RNA using a

123

SMARTer™ RACE cDNA Amplification Kit (TaKaRa, Japan). 5’RACE was carried

124

out by nested PCR using OBP1-5-GSP1/ OBP2-5-GSP1 and the outer Primer from

125

the primary PCR reactions and the second PCR reactions were carried out using

126

OBP1-5-GSP2/ OBP2-5-GSP2 combined with the inner Primer. 3’ RACE was

127

conducted by nested PCR using OBP1-3-GSP1/ OBP2-3-GSP1 and Long Primer for

128

the primary PCR and the second PCR was carried out using

129

OBP1-3-GSP2/OBP2-3-GSP2 combined with the Short Primer. For the second PCR

130

amplification, a 50-fold dilution of the primary reaction mixture was used as template.

131

All thermal cycling followed a standard program: 94 °C for 3 min; 35 cycles of 30 s at

132

94 °C, 30 s at 54-58 °C and 60 s at 72 °C; 72 °C for 10 min. All amplified products

133

were cloned in pMD18-T Vector (TaKaRa, Japan) and sequenced completely. Based

134

on the sequences obtained from the 3′- and 5′-RACE, the putative full-length OBP

135

genes were amplified with specific primers OBP-F and OBP-R. The amplified

136

products were gel-purified, and cloned into pMD18-T vectors (TaKaRa, Japan) and

137

sequenced completely. 7

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

138

Sequence analysis. The ExPASy Proteomics Server

139

(http://cn.expasy.org/tools/pi_tool.html) was used to compute the isoelectric points

140

and molecular weights of deduced protein sequences. The signal peptide was

141

predicted by SignalP 3 Serve (http://www.cbs.dtu.dk/service/SignalP/). Alignments of

142

nucleotide and amino acid sequences were made using DNAMAN software. 3D

143

structure analysis was performed with the program SWISS MODEL

144

(https://swissmodel.expasy.org/).

145

Expression and purification of recombinant BodoOBP1 and BodoOBP2 proteins.

146

The coding sequences of mature BodoOBPs were amplified by PCR using specific

147

primers carrying endonuclease restriction sites (Table S1). The purified PCR products

148

of two OBPs were subcloned into pMD18-T vector (TaKaRa, Japan ), then the

149

plasmids containing the correct sequence were digested by specific restriction

150

enzymes. The digested fragments were purified and ligated into the expression vectors

151

pET28a (+) (Novagen, Madison, WI) previously digested with the same enzymes

152

(TaKaRa, Japan). These constructed vectors were confirmed by sequencing. The

153

recombinant plasmids were transformed into Escherichia coli BL21(DE3) competent

154

cells. The expression of recombinant proteins was induced at 28 °C for 8 h with a

155

final concentration of 0.5 mM isopropyl β-d-1-thiogalactopyranoside (IPTG).

156

Bacteria cells were collected by centrifugation at 6,000 g for 10 min, suspended in

157

balance buffer (300 mM NaCl, 50 mM NaH2PO4, 10 mM imidazole, 10 mM Tris

158

base, pH8.0) and sonicated in ice (15 s, 20 passes), then centrifuged again at 12,000 g

159

for 20 min. Sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) 8

ACS Paragon Plus Environment

Page 8 of 48

Page 9 of 48

160

Journal of Agricultural and Food Chemistry

showed that the OBPs were expressed mainly in inclusion bodies.

161

The inclusion bodies were washed by 50 mM Tris buffer (pH 6.8) with 0.2%

162

Triton X-100 and dissolved in 5mL of 8M urea. The proteins were refolded according

163

to the protocol of Prestwich45. Protein purification was accomplished by Ni-NTA

164

Resin (TransGen Biotech, Beijing, China). Purified recombinant proteins were

165

analyzed by SDS-PAGE and dialyzed overnight on ice with 10 mM

166

phosphate-buffered saline (PBS) (pH 7.4).

167

Fluorescence competitive binding assay. 10 typical chemical components of volatile

168

oils from the Chinese chive46 were purchased from Aladdin (China Shanghai), HX-R

169

(China, Chengdu), Sigma-Aldrich (USA) and TCI (Japan). The purity of all ligands

170

and the fluorescent probe N-phenyl-1-naphthylamine(1-NPN) is higher than 95% and

171

all of these were dissolved independently in GC grade methanol (Aladdin, China) to

172

create 1mM stock solutions. Fluorescence competitive binding assays were performed

173

at 25 °C on a Hitachi F-4500 spectrofluorimeter with a 1 cm light path fluorimeter

174

quartz cuvette and 10-nm slits for both excitation and emission. The fluorescent probe

175

was excited at 337 nm and emission spectra were recorded between 350 and 600 nm.

176

The affinity of the fluorescent probe to each OBP was measured by titrating a 2 μM

177

solution of the protein in 50 mM Tris-HCl (pH 6.8) with aliquots of 1 mM 1-NPN to

178

final concentrations of 2-20 μM. In competitive binding assays, the ligands from 0 to

179

50μM were added into each protein (2 μM) and 1-NPN (2 μM).

180 181

The fluorescence intensities at the maximum fluorescence emission between 350 and 600 nm were plotted against the free ligand concentration to determine the 9

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

182

binding constants. All data were analysed based on the assumption that the protein

183

was 100% active and had a stoichiometry of 1 : 1 (protein : ligand) at saturation. The

184

binding curves with 1-NPN were linearized using a Scatchard plot. The dissociation

185

constants of the competitors were calculated from the corresponding IC50 values

186

based on the following equation: Ki=[IC50]/(1+[1-NPN]/K1-NPN), where [1-NPN] is

187

the free concentration of 1-NPN and K1-NPN is the dissociation constant of the

188

complex protein/1-NPN47. If the IC50 value of the candidate competitive ligand

189

exceeded 30 μM, the further calculation of Ki was not considered.

190

Quantitative real-time PCR. Quantitative real-time PCR (qRT-PCR) was

191

performed using a Cycler IQ Real-Time PCR detection system (Bio-Rad, Hercules,

192

CA, USA). 1μg of each RNA sample was employed as template to synthesize

193

first-strand cDNA using a PrimeScript™ RT Reagent Kit with gDNA Eraser (Perfect

194

Real Time) (TaKaRa, Japan). Each reaction included 1μL of template, 1 μL (0.6 μL

195

for BodoOBP2) of each primer (10 mM), 10μL 2× UltraSYBR Mixture (KWBIO,

196

Beijing, China) and added RNase-free water to a total volume of 20μL. Amplification

197

conditions were: 95 °C for 10 min; followed by 40 cycles of 95 °C for 15 s, and 60 °C

198

for 60 s. A melting curve analysis was conducted by a protocol that varied

199

temperature from 55 to 95 °C with increments of 0.5 °C every 30 s with continuous

200

fluorescence acquisition to ensure consistency and specificity of the amplified product.

201

The amplification efficiency (E) of all primer pairs was validated by constructing a

202

standard curve using a 10-fold serial dilution of cDNA template using the equation: E

203

= (10[-1/slope]-1) ×100. No-template controls were included. Each experiment contained 10

ACS Paragon Plus Environment

Page 10 of 48

Page 11 of 48

Journal of Agricultural and Food Chemistry

204

three biological replicates (three technical replications for each biological replication).

205

Relative expression levels were determined using the 2−ΔΔCT method48. The 18S,

206

RPS15 and RPS7 was used as a reference to normalize the expression of the target

207

genes across different developmental stage. The 18S, RPS7 and Tubulin was used as a

208

reference to normalize the expression of the target genes across various tissues.

209

Results are shown as means ± SE. Statistical analyses of the gene expression data

210

were performed using ANOVA in SPSS 20.0 (IBM Inc., Chicago, Illinois, USA).

211

dsRNA synthesis and RNA interference assay. Primers of dsRNA fragments were

212

designed according to the two OBP sequences, and the selected region for dsRNA

213

synthesis was different from the region for qRT-PCR. These primers contained a T7

214

polymerase binding sequence required for dsRNA synthesis. PCR conditions were

215

94 °C for 3 min, followed by 35 cycles of 94 °C for 30 s, 60 °C for 30 s and 72 °C for

216

60 s, with a final extension at 72 °C for 10 min. The amplified PCR products were

217

separated on an agarose gel and purified by a Biospin Gel Extraction Kit (Bioer

218

Technology, Hangzhou, China). dsRNAs were synthesized using the T7

219

RiboMAXTM Express RNAi System (Promega) following the manufacturer’s

220

instructions. dsRNA was purified and precipitated using 3M sodium acetate (pH 5.2)

221

and isopropanol, then washed using 0.5ml of cold 70% ethanol. The purified dsRNA

222

was eluted into diethyl pyrocarbonate (DEPC) -treated nuclease-free water and stored

223

at -80 °C until use. As controls, double strand Green fluorescent protein (dsGFP,

224

GenBank accession number: KF410615.1) was synthesized using the protocol

225

described above. The dsRNA concentration was determined by a Nanodrop™ 1000 11

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

226

(Thermoscientific, Lithuania). The primers used for dsRNA synthesis are shown in

227

Table S1.

228

30 nL dsRNA (6 μg/μL ) was injected into the abdomen of 1-day-old pupae

229

using a Nanoject III (Drummond Scientific Company, USA) and each insect was

230

injected only once. Pupae were injected with an equivalent amount of dsGFP as a

231

control. The injected pupae were reared under the conditions described above and

232

adults emerged on the third day of the pupal stage. 5 adults (male and female

233

separately) were collected into RNase-free tubes and immersed into in liquid nitrogen

234

at 72 h after dsRNA injection. RNA was extracted with RNAiso plus Reagent

235

(TaKaRa, Japan) to measure the transcript levels of BodoOBP1 and BodoOBP2 by

236

qRT-PCR as described above.

237

Y-tube olfactometer assay. In order to study the role of OBPs in B. odoriphaga

238

olfactory recognition behavior, a Y-tube olfactometer with diameter 2.5 cm (each arm

239

with length of 15 cm and stem with length of 15 cm) was used in this study. The air

240

was first filtered through active carbon and then humidified with doubly distilled

241

water before entering the Y-tube. The air that travelled through both arms of the

242

olfactometer was pumped in at a flow rate of 300 mL min-1. Each individual virgin

243

adult was placed at the middle position of the olfactometer stem, and its choice was

244

recorded when it passed the middle position of an arm and stayed for 30 s. Virgin

245

adults that made no choice within 5 min were considered nonresponders. In these

246

assays, 10 μL of the dipropyl trisulfide (10 μg/μL in methanol (HPLC)) applied to

247

6mm diameter filter paper or virgin females (15 individuals, 24 hours after emergence) 12

ACS Paragon Plus Environment

Page 12 of 48

Page 13 of 48

Journal of Agricultural and Food Chemistry

248

were placed in the incoming airstream of one of the arms of the Y-tube olfactometer.

249

Each group included at least 30 individuals with 3 replicates. The number of adults

250

that travelled upwind towards the odour sources in each arm was recorded. The time

251

spent by individuals moving toward the odorant source was measured using a digital

252

chronometer. The positions of the odour sources were exchanged after every 15 tests.

253

The Y-tube olfactometer was cleaned with 75% ethanol and dried after every

254

exchange. An attraction index was calculated from each replicate using the formula:

255

(number of individuals in the treatment arm – number of individuals in the control

256

arm)/total number of tested individuals49.

257

RESULTS

258

cDNA sequence analysis of BodoOBP1 and BodoOBP2

259

The full length cDNA sequences of BodoOBP1 and BodoOBP2 were obtained

260

by PCR and submitted to GenBank (accession numbers: KX099613 and KX099614).

261

BodoOBP1 and BodoOBP2 are 622 bp and 690 bp in length, respectively, and both

262

contain open reading frames (ORF) of 453 bp which encode a protein of 150 amino

263

acids. The calculated molecular masses of BodoOBP1 and BodoOBP2 proteins are

264

17.3 kDa and 17.5 kDa, and the predicted isoelectric points are 5.67 and 6.03,

265

respectively. Both proteins all have predicted signal peptides of 17 amino acids at the

266

N-terminus.

267

BodoOBP1 and BodoOBP2 proteins have the typical six cysteine signature

268

(C1-X25−30-C2-X3-C3-X36−42-C4-X8−14-C5-X8-C6) and belong to the Classic OBP

269

subfamily (Fig. 1). The amino acid sequence alignment revealed that both BodoOBPs 13

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

270

share very high similarities (identity >42%) with OBPs from other Dipteran insects

271

(Fig. 1). The amino acid sequence of BodoOBP1 showed 86.7% identity with

272

BodoOBP2. Both BodoOBPs have approximately 58% and 55% identities with

273

AgamOBP1 (GenBank accession number: AAL84179.1), respectively (Fig. 2).

274

AgamOBP1 was used as template to predict the 3D structure of BodoOBP1 and

275

BodoOBP2 proteins according to SWISS MODEL (Fig. 3). 3D models of BodoOBP1

276

and BodoOBP2 proteins showed that both have the typical predicted structure of

277

OBPs, with six α-helices (labeled α1 to α6) surrounding a central hydrophobic cavity.

278

Four helices (α1, α4, α5 and α6) converge to form a binding pocket. The narrow end

279

of the binding pocket is open with the opposite end being capped by the α3 helix. The

280

binding pocket is delimited by the hydrophobic sides of these helices, and therefore

281

well suited to being a binding site for hydrophobic ligands. The 3D structures of

282

BodoOBP1 and BodoOBP2 are similar, suggesting that they may have similar

283

functions.

284

Bacterial expression and purification of BodoOBP1 and BodoOBP2

285

The genes encoding BodoOBP1 and BodoOBP2 were amplified from cDNA and

286

ligated into the expression vector pET-28a (+), as described in the Materials and

287

Methods section. The two OBPs were mainly present in inclusion bodies and highly

288

expressed, yielding approximately 30 mg/L culture. The denaturation, renaturation

289

and purification were performed as described in the Materials and Methods. The

290

SDS-PAGE analysis showed that the molecular weights of the recombinant proteins

291

(contain His-Tag) were approximately 19 kDa (Fig. 4). 14

ACS Paragon Plus Environment

Page 14 of 48

Page 15 of 48

292 293

Journal of Agricultural and Food Chemistry

Ligand-binding assays To investigate the ligand binding properties of BodoOBP1 and BodoOBP2 to

294

different host volatiles, the fluorescence competitive binding assay was further

295

performed using 1-NPN as the fluorescence probe. The binding curves with 1-NPN

296

and Scatchard plots indicated that the bindings between the proteins and 1-NPN were

297

saturable with a single population of binding sites, without an apparent allosteric

298

effect (Fig. 5). The results showed that both recombinant BodoOBPs could well

299

interact with 1-NPN with dissociation constants of 4.03±0.46 μM and 6.91±0.61 μM.

300

The above results indicated that using 1-NPN as the fluorescence probe is suitable for

301

a fluorescence competitive binding assay.

302

Ten chemicals were selected based on previous reports (Table S2) and used in

303

the competitive binding assay. Competitive binding curves of these ligands to

304

BodoOBP1 and BodoOBP2 are shown in Fig. 6. When the concentration of dipropyl

305

trisulfide reached 15.46 μM and 12.06 μM, respectively, the fluorescence intensity of

306

the BodoOBP1-2 /1-NPN complex rapidly decreased to approximately 50%. These

307

results showed that the both BodoOBPs have moderate binding affinity to dipropyl

308

trisulfide, with Ki values of 12.39±0.83 μM and 10.53±0.36 μM, respectively (Table

309

S2). For other chemicals, the two BodoOBPs showed weak or no binding affinity

310

(Ki > 30 μM).

311

Tissue Expression Patterns of BodoOBPs

312 313

The quantitative real-time PCR (qRT-PCR) was used to measure transcript levels of BodoOBP1 and BodoOBP2 in different tissues. The results revealed that both 15

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

314

genes were significantly more highly expressed in the antennae of both sexes than in

315

other tissues (Fig. 7). Furthermore, expression levels of both genes were much higher

316

in male than in female antennae, approximately 3-fold higher for BodoOBP1 and

317

9-fold higher for BodoOBP2. The expressions of the two OBP genes were also

318

detected in the head of both sexes.

319

Expression patterns of BodoOBPs in different life stages

320

The expression patterns of BodoOBP1 and BodoOBP2 in different life stages

321

were measured by RT-qPCR. The transcript levels of both BodoOBPs were higher in

322

the male adult stage than in the other developmental stages (Fig. 8). Compared to

323

females, expression of both BodoOBPs was significantly higher in the adult and pupal

324

stage of males.

325

The results showed that expression of BodoOBP1 and BodoOBP2 in males was

326

significantly higher in the 24 h after eclosion than in the female (Fig. 9). Expression

327

of both BodoOBPs was down-regulated in 0h post-mating males but up-regulated in

328

1h post-mating males compared to virgin males. Compared to virgin females, the

329

expression of both BodoOBPs was up-regulated in 1h post-mating females. No

330

significant differences were detected in the expression of OBP genes between mature

331

virgin females and 0h post-mating females.

332

Effect of dsRNA treatment on B. odoriphaga behavior

333

To test the effect of reducing transcript levels of BodoOBP1 and BodoOBP2

334

genes on B. odoriphaga behavior, we injected pupae with dsOBPs, and measured

16

ACS Paragon Plus Environment

Page 16 of 48

Page 17 of 48

Journal of Agricultural and Food Chemistry

335

transcript levels of both BodoOBPs at 72 h post-injection. After treatment with

336

dsOBPs, BodoOBP expression was reduced by 85% or more in adults (Fig. 10).

337

To detect the role of BodoOBP1 and BodoOBP1 in response to sex pheromones,

338

we used virgin female adults as odour sources. The Y‐tube assay indicated that ability

339

of dsOBPs-males to detect females was significantly reduced when compared to

340

control insects (Fig. 11A). The attraction indices of males treated with dsOBP1 to

341

females were significantly lower than for those injected with dsOBP2. Virgin males

342

responded to virgin females with an attraction index of 53%, whereas the attraction

343

index was only 4% and 31%, respectively, for dsOBP1-treated and dsOBP2-treated

344

virgin males.

345

The attraction of plant volatile dipropyl trisulfide to dsGFP-treated male or

346

female adults was 32% and 31%, respectively, which was not significantly different

347

(Fig. 11B). However, the attraction of dipropyl trisulfide to dsOBP1-treated males or

348

dsOBP1-females was significantly lower than dsGFP-treated insects, with attraction

349

index of 9% and 7%, respectively, which was similar to results for attraction of

350

dipropyl trisulfide to dsOBP2-treated male or female adults (Fig. 11B).

351

The time spent by dsOBPs-treated virgin male adults moving toward females or

352

dipropyl trisulfide was significantly longer than for dsGFP-treated groups (Fig.12A,

353

B). The response time of dsOBPs-treated female individuals to dipropyl trisulfide

354

increased, but no significant difference was observed (Fig. 12C).

355

DISCUSSION

356

The olfactory recognition system plays a vital role in interactions of insects with 17

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

357

their environment, and is important for their survival and reproduction4. Our study

358

demonstrates that OBPs play a central role in the initial phases of olfactory perception

359

in B. odoriphaga. OBPs are responsible for transporting hydrophobic ligands to

360

specific ORs4. This suggests that OBPs could be useful for controlling insect

361

populations by blocking the first step of the olfactory recognition process. Although

362

the chemosensory genes of B. odoriphaga have been detected previously in antennal

363

and whole-body transcriptomes14, there was no previous information on expression

364

profiling of OBPs in B. odoriphaga across different developmental stages, which is

365

necessary to determine the role of these proteins on reproduction. Identification and

366

functional analysis of odorant-binding proteins is needed to provide new and more

367

environmentally friendly pest control strategies utilizing olfactory interference.

368

In this study, two OBP genes, BodoOBP1 and BodoOBP2 were cloned and

369

characterized from B. odoriphaga, and recombinant proteins were expressed in E.

370

coli. The amino acid sequence alignment revealed that both BodoOBPs are similar

371

(identity >42%) to OBPs from other dipteran insects. Both BodoOBPs have high

372

homology with Drosophila melanogaster OBPs OS-E and OS-F, which clustered

373

together with DmelOBP83a/83b14, an OBP group associated with the detection of

374

volatile pheromones in D. melanogaster27, 50. Binding assays showed that dipropyl

375

trisulfide, which has been isolated and identified from the B. odoriphaga host plant

376

Allium tuberosum, had the highest binding affinity of the tested ligands to BodoOBP1

377

and BodoOBP2. For other chemicals, these two BodoOBPs showed weak binding or

378

no binding affinity. Thus, we infer that both BodoOBPs are involved in host-seeking 18

ACS Paragon Plus Environment

Page 18 of 48

Page 19 of 48

379

Journal of Agricultural and Food Chemistry

behavior.

380

Previous work has reported that the transcripts of BodoOBP1 and BodoOBP2

381

were expressed specifically at very high levels in the female and male antennae14, as

382

also observed in our study, which implies that the two proteins are involved in

383

olfactory chemoreception. Interestingly, BodoOBP1 and BodoOBP2 were expressed

384

much more highly in male antennae than female antennae (Fig. 7), suggesting that

385

they could also play a role in male chemosensory behavior, such as sex pheromone

386

detection. This has been reported for the mosquito Aedes aegypti, in which OBP10

387

enriched in antennae may play a role in adult male chemosensory behavior51. The

388

results of qRT-PCR indicated that expression of BodoOBP1-2 differs during different

389

developmental stages of B. odoriphaga, with higher expression in male adults than in

390

other developmental stages (Fig. 8). Previous study of D. melanogaster showed that

391

Obp19a-d was expressed only in adults52. A similar result was observed in

392

Chrysomya megacephala with qPCR results showing that higher expression of two

393

OBP genes in adults than in larvae53. The male adult-biased expression observed for

394

BodoOBP1-2 suggests putative roles for these proteins in mate location.

395

For B. odoriphaga, we found a difference in expression levels between virgin

396

and post-mating adults for BodoOBP1. Both BodoOBPs were down-regulated in 0h

397

post-mating males but up-regulated in 1h post-mating males compared to virgin males

398

(Fig. 9). Similar results were observed in Anastrepha fraterculus post-mating males,

399

in which OBP56d and OBP57c were up-regulated in the 6h post-mating54. Higher

400

expression of OBPs in post-mating males would suggest that these genes could be 19

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

401

related to mate seeking including responding to the female pheromones55. Due to the

402

habit of multiple mating in B. odoriphaga male adults, we speculate that the higher

403

expression of both BodoOBPs in the 1h post-mating males could signify preparation

404

for re-mating. The expression of both BodoOBPs was up-regulated in 1h post-mating

405

females compared to virgin females (Fig. 9). Although in A. fraterculus females, most

406

OBP genes were down-regulated post-mating, in A. obliqua females the majority of

407

OBPs were up-regulated post-mating54, the same pattern observed in D. melanogaster

408

females52. Previous research indicated that behavioral changes in adults due to mating

409

were mainly related to responses of female moths to plant volatiles and of male moths

410

responding to sex pheromones56-59. Mature virgin females of the Mediterranean fruit

411

fly, Ceratitis capitata, have a strong preference for a male-produced pheromone,

412

whereas mated females are more attracted to host plant odors60, 61.

413

To test the effect of reducing transcript levels of BodoOBP1 and BodoOBP2

414

genes on B. odoriphaga behavior, we performed RNA interference of both

415

BodoOBPs. Our Y-tube olfactometer assay indicated that ability of dsOBP1-males to

416

detect females was significantly reduced when compared to control insects (Fig. 11).

417

Similar results were obtained for dsOBP2-males. Our findings are consistent with the

418

previous finding that in Rhodnius prolixus, RproOBP27 is significantly expressed in

419

male antennae33. After treatment with dsOBP27, ability of male adults to detect

420

females was reduced, and adult males spent significantly less time close to females as

421

compared to controls, indicating that RproOBP27 may be involved in the perception

422

of semiochemicals related to mate finding33. Additionally, our results indicated that 20

ACS Paragon Plus Environment

Page 20 of 48

Page 21 of 48

Journal of Agricultural and Food Chemistry

423

attraction of plant volatile dipropyl trisulfide to dsOBPs-treated adults was

424

significantly lower than dsGFP-treated insects (Fig. 11). Therefore, BodoOBP1 and

425

BodoOBP2 may be required for host plant recognition. In Holotrichia oblita,

426

dsOBP-treated female beetles showed significantly lower attraction to (E)-2-hexenol

427

and phenethyl alcohol than those treated with GFP-dsRNA, indicating that

428

HoblOBP13 and HoblOBP9 are essential for H. oblita perception of these volatile

429

compounds49. Taken together, the results presented here strongly suggest that both

430

BodoOBPs are likely involved in the perception of sex pheromone(s) and

431

host-seeking.

21

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

432

ACKNOWLEDGMENT

433

We thank Dr. Chris Dietrich (University of Illinois, IL, USA) and Dr. John Richard

434

Schrock (Emporia State University, Emporia, KS, USA) for their comments on an

435

earlier draft of this paper.

436

AUTHOR CONTRIBUTIONS

437

W.D. and C.N.Z. designed the experiment. B.W.T., S.L.T. and C.N.Z., performed the

438

experiments. B.W.T., W.D. and C.N.Z. prepared the manuscript. All authors

439

contributed to data analysis.

440

FUNDING

441

The research was supported by the National Natural Science Foundation of China

442

(31672037) and the Special Fund for Agro-scientific Research in the Public Interest

443

from the Ministry of Agriculture of China (201303027).

444

NOTES

445

The authors declare no competing financial interest.

446

SUPPORTING INFORMATION

447

The Supporting Information is available free of charge on the ACS Publications

448

website.

449

Table S1. Primers used for cloning of genes, qRT-PCR and RNAi. Table S2. Binding

450

affinities of selected host plant compounds to the recombinant proteins.

22

ACS Paragon Plus Environment

Page 22 of 48

Page 23 of 48

Journal of Agricultural and Food Chemistry

451

REFERENCES

452

(1) Pelosi, P.; Zhou, J. J.; Ban, L. P.; Calvello, M., Soluble proteins in insect chemical

453

communication. Cellular & Molecular Life Sciences 2006, 63, 1658-1676.

454

(2) Field, L. M.; Pickett, J. A.; Wadhams, L. J., Molecular studies in insect olfaction.

455

Insect Molecular Biology 2000, 9, 545-551.

456

(3) Hansson, B.; Stensmyr, M., Evolution of Insect Olfaction. Neuron 2011, 72,

457

698-711.

458

(4) Leal, W. S., Odorant reception in insects: roles of receptors, binding proteins, and

459

degrading enzymes. Annual Review of Entomology 2013, 58, 373-391.

460

(5) Pelosi, P.; Iovinella, I.; Felicioli, A.; Dani, F. R, Soluble proteins of chemical

461

communication: an overview across arthropods. Frontiers in physiology 2014, 5, 320.

462

(6) Zhu, J.; Guo, M. B.; Ban, L. P.; Song, L. M.; Liu, Y.; Pelosi, P.; Wang, G. R.,

463

Niemann-Pick C2 Proteins: A New Function for an Old Family. Frontiers in

464

physiology 2018, 9, 52.

465

(7) Pelosi, P.; Calvello, M.; Ban, L., Diversity of odorant-binding proteins and

466

chemosensory proteins in insects. Chemical Senses 2005, 30 Suppl 1, i291-2.

467

(8) Pelosi, P.; Iovinella, I.; Zhu, J.; Wang, G.; Dani, F. R., Beyond chemoreception:

468

diverse tasks of soluble olfactory proteins in insects. Biological Reviews 2018, 93,

469

184-200.

470

(9) Fan, J.; Francis, F.; Liu, Y.; Chen, J. L.; Cheng, D. F., An overview of

471

odorant-binding protein functions in insect peripheral olfactory reception. Genetics &

472

Molecular Research 2011, 10, 3056-69. 23

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

473

(10) Li, Z. Q.; Zhang, S.; Ma, Y.; Luo, J. Y.; Wang, C. Y.; Lv, L. M.; Dong, S. L.;

474

Cui, J. J., First Transcriptome and Digital Gene Expression Analysis in Neuroptera

475

with an Emphasis on Chemoreception Genes in Chrysopa pallens (Rambur). PloS one

476

2013, 8, e67151.

477

(11) Niu, D. J.; Liu, Y.; Dong, X. T.; Dong, S. L., Transcriptome based identification

478

and tissue expression profiles of chemosensory genes in Blattella germanica

479

(Blattaria: Blattidae). Comparative biochemistry and physiology. Part D, Genomics &

480

proteomics 2016, 18, 30-43.

481

(12) Pelosi, P.; Maida, R., Odorant-binding proteins in insects. Comparative

482

Biochemistry & Physiology Part B Biochemistry & Molecular Biology 1995, 111,

483

503-514.

484

(13) Xu, Y. L.; He, P.; Zhang L; Fang, S. Q.; Dong, S. L.; Zhang, Y. J.; Li, F.,

485

Large-scale identification of odorant-binding proteins and chemosensory proteins

486

from expressed sequence tags in insects. BMC genomics 2009, 10, 632.

487

(14) Zhao, Y. H.; Ding, J. F.; Zhang, Z. Q.; Liu, F.; Zhou, C. G.; Mu, W., Sex- and

488

Tissue-Specific Expression Profiles of Odorant Binding Protein and Chemosensory

489

Protein Genes in Bradysia odoriphaga (Diptera: Sciaridae). Frontiers in Physiology

490

2018, 9, 107.

491

(15) Zhou, J. J., Odorant-binding proteins in insects. Vitamins & Hormones 2010, 83,

492

241.

493

(16) Leal, W. S.; Nikonova, L.; Peng, G., Disulfide structure of the pheromone

494

binding protein from the silkworm moth, Bombyx mori. FEBS letters 1999, 464, 24

ACS Paragon Plus Environment

Page 24 of 48

Page 25 of 48

Journal of Agricultural and Food Chemistry

495

85-90.

496

(17) Scaloni, A.; Monti, M.; Angeli, S.; Pelosi, P., Structural analysis and

497

disulfide-bridge pairing of two odorant-binding proteins from Bombyx mori.

498

Biochemical and biophysical research communications 1999, 266, 386-391.

499

(18) Hull, J. J.; Perera, O. P.; Snodgrass, G. L., Cloning and expression profiling of

500

odorant-binding proteins in the tarnished plant bug, Lygus lineolaris. Insect Molecular

501

Biology 2014, 23, 78-97.

502

(19) Zhang, Y. N.; Zhu, X. Y.; Zhang, Q.; Yin, C. Y.; Dong, Z. P.; Zuo, L. H.; Deng,

503

D. G.; Sun, L.; Li, X. M., De novo assembly and characterization of antennal

504

transcriptome reveal chemosensory system in Nysius ericae. Journal of Asia-Pacific

505

Entomology 2016, 19, 1077-1087.

506

(20) Wang, B.; Liu, Y.; Wang, G. R., Chemosensory genes in the antennal

507

transcriptome of two syrphid species, Episyrphus balteatus and Eupeodes corollae

508

(Diptera: Syrphidae). BMC genomics 2017, 18, 586.

509

(21) Zhou, J. J.; Vieira, F. G.; He, X. L.; Smadja, C.; Liu, R.; Rozas, J.; Field, L. M.,

510

Genome annotation and comparative analyses of the odorant ‐ binding proteins and

511

chemosensory proteins in the pea aphid Acyrthosiphon pisum. Insect Molecular

512

Biology 2010, 19, 113-122.

513

(22) Gu, S. H.; Wang, S. P.; Zhang, X. Y.; Wu, K. M.; Guo, Y. Y.; Zhou, J. J.; Zhang,

514

Y. J., Identification and tissue distribution of odorant binding protein genes in the

515

lucerne plant bug Adelphocoris lineolatus (Goeze). Insect Biochemistry & Molecular

516

Biology 2011, 41, 254-263. 25

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

517

(23) Gu, S. H.; Wu, K. M.; Guo, Y. Y.; Field, L. M.; Pickett, J. A.; Zhang, Y. J.;

518

Zhou, J. J., Identification and expression profiling of odorant binding proteins and

519

chemosensory proteins between two wingless morphs and a winged morph of the

520

cotton aphid Aphis gossypii glover. PloS one 2013, 8, e73524.

521

(24) He, P.; Zhang, J.; Liu, N. Y.; Zhang, Y. N.; Yang, K.; Dong, S. L., Distinct

522

expression profiles and different functions of odorant binding proteins in Nilaparvata

523

lugens Stal. PloS one 2011, 6, e28921.

524

(25) Li, Z. Q.; Zhang, S.; Luo, J. Y.; Wang, S. B.; Dong, S. L.; Cui, J. J.,

525

Odorant-binding proteins display high affinities for behavioral attractants and

526

repellents in the natural predator Chrysopa pallens. Comparative biochemistry and

527

physiology. Part A, Molecular & integrative physiology 2015, 185, 51-57.

528

(26) Galindo, K.; Smith, D. P., A large family of divergent Drosophila

529

odorant-binding proteins expressed in gustatory and olfactory sensilla. Genetics 2001,

530

159, 1059-1072.

531

(27) Shanbhag, S. R.; Park, S. K.; Pikielny, C. W.; Steinbrecht, R. A., Gustatory

532

organs of Drosophila melanogaster: fine structure and expression of the putative

533

odorant-binding protein PBPRP2. Cell and tissue research 2001, 304, 423-37.

534

(28) Ishida, Y.; Ishibashi, J.; Leal, W. S., Fatty acid solubilizer from the oral disk of

535

the blowfly. PloS one 2013, 8, e51779.

536

(29) Jeong, Y. T.; Shim, J.; Oh, S. R.; Yoon, H. I.; Kim, C. H.; Moon, S. J.; Montell,

537

C., An odorant-binding protein required for suppression of sweet taste by bitter

538

chemicals. Neuron 2013, 79, 725-737. 26

ACS Paragon Plus Environment

Page 26 of 48

Page 27 of 48

Journal of Agricultural and Food Chemistry

539

(30) Sparks, J. T.; Bohbot, J. D.; Dickens, J. C., The genetics of chemoreception in

540

the labella and tarsi of Aedes aegypti. Insect biochemistry and molecular biology

541

2014, 48, 8-16.

542

(31) Sun, L.; Wang, Q.; Wang, Q.; Dong, K.; Xiao, Y.; Zhang, Y. J., Identification

543

and Characterization of Odorant Binding Proteins in the Forelegs of Adelphocoris

544

lineolatus (Goeze). Frontier in Physiology 2017, 8, 735.

545

(32) He, P.; Li, Z. Q.; Zhang, Y. F.; Chen, L.; Wang, J.; Xu, L.; Zhang, Y. N.; He, M.,

546

Identification of odorant-binding and chemosensory protein genes and the ligand

547

affinity of two of the encoded proteins suggest a complex olfactory perception system

548

in Periplaneta americana. Insect Molecular Biology 2017, 26, 687-701.

549

(33) Oliveira, D. S.; Brito, N. F.; Franco, T. A.; Moreira, M. F.; Leal, W. S.; Melo, A.

550

C. A., Functional characterization of odorant binding protein 27 (RproOBP27) from

551

Rhodnius prolixus antennae. Frontiers in Physiology 2018, 9, 1175.

552

(34) Pregitzer, P.; Zielonka, M.; Eichhorn, A. S.; Jiang, X.; Krieger, J.; Breer, H.,

553

Expression of odorant ‐ binding proteins in mouthpart palps of the desert locust

554

Schistocerca

555

https://doi.org/10.1111/imb.12548

556

(35) Ban, L. P.; Napolitano, E.; Serra, A.; Zhou, X. H.; Iovinella, I.; Pelosi, P.,

557

Identification of pheromone ‐ like compounds in male reproductive organs of the

558

oriental locust Locusta migratoria. Biochemical and Biophysical Research

559

Communications 2013, 437, 620–624.

560

(36) Sun, Y. L.; Huang, L. Q.; Pelosi, P.; Wang, C. Z., Expression in antennae and

gregaria.

Insect

Molecular

27

ACS Paragon Plus Environment

Biology

2018,

Journal of Agricultural and Food Chemistry

561

reproductive organs suggests a dual role of an odorant-binding protein in two sibling

562

Helicoverpa species. PLoS One 2012, 7, e30040.

563

(37) Zhang, J.; Luo, D. Y.; Wu, P.; Li, H. Z.; Zhang, H. Y.; Zheng, W. W.,

564

Identification and expression profiles of novel odorant binding proteins and functional

565

analysis of OBP99a in Bactrocera dorsalis. Archives of Insect Biochemistry and

566

physiology 2018, https://doi.org/10.1002/arch.21452

567

(38) Li, W. X.; Yang, Y. T.; Xie, W.; Wu, Q. J.; Xu, B. Y.; Wang, S. L.; Zhu, X.;

568

Wang, S. J.; Zhang, Y. J., Effects of Temperature on the Age-Stage, Two-Sex Life

569

Table of Bradysia odoriphaga (Diptera: Sciaridae). Journal of Economic Entomology

570

2015, 108, 126-134.

571

(39) Yang, Y. T.; Li, W. X.; Xie, W.; Wu, Q. J.; Xu, B. Y.; Wang, S. L.; Li, C. R.;

572

Zhang, Y. J., Development of Bradysia odoriphaga (Diptera: Sciaridae) as affected

573

by humidity: an age–stage, two-sex, life-table study. Applied Entomology & Zoology

574

2015, 50, 3-10.

575

(40) Zhang, H.; Mallik, A.; Zeng, R. S., Control of Panama disease of banana by

576

rotating and intercropping with Chinese chive (Allium tuberosum Rottler): role of

577

plant volatiles. Journal of Chemical Ecology 2013, 39, 243-252.

578

(41) Ma, J.; Moens, M.; Han, R.; Clercq, P. D., Efficacy of entomopathogenic

579

nematodes (Rhabditida: Steinernematidae and Heterorhabditidae) against the chive

580

gnat, Bradysia odoriphaga. Journal of Pest Science 2013, 86, 551-561.

581

(42) Li, H. J.; He, X. K.; Zeng, A. J.; Liu, Y. J.; Jiang S. R., Bradysia odoriphaga

582

copulatory behavior and evidence of a female sex pheromone. Journal of Agricultural 28

ACS Paragon Plus Environment

Page 28 of 48

Page 29 of 48

Journal of Agricultural and Food Chemistry

583

and Urban Entomology, 2007, 24, 27-34.

584

(43) Li, H. J., He, X. K.; Zeng, A. J.; Liu, Y. J.; Jiang, S. R., Evidence of a female sex

585

pheromone in Bradysia odoriphaga (Diptera: Sciaridae). The Canadian Entomologist,

586

2008, 140, 324 -326.

587

(44) Kafil Uddin, M. D.; Yin, X. W.; Zhang, L., Courtship and mating behavior of the

588

Chinese chive fly, Bradysia odoriphaga (Diptera: Sciaridae) and evidence of female

589

sex pheromone. Pakistan Journal of Zoology 2016, 48 ,1543-1548.

590

(45) Prestwich, G. D., Bacterial expression and photoaffinity labeling of a pheromone

591

binding protein. Protein Science 1993, 2, 420-428.

592

(46) Wang, H. M.; Feng, J., Identification on chemical components of volatile oil of

593

chive by GC/MS. Journal of Tianjin Medical University 2002, 8, 191-192.

594

(47) Campanacci, V.; Krieger, J.; Bette, S.; Sturgis, J. N.; Lartigue, A.; Cambillau, C.;

595

Breer, H.; Tegoni, M., Revisiting the specificity of Mamestra brassicae and Antheraea

596

polyphemus pheromone-binding proteins with a fluorescence binding assay. Journal

597

of Biological Chemistry 2001, 276, 20078.

598

(48) Livak, K. J.; Schmittgen, T. D., Analysis of relative gene expression data using

599

real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 2001, 25,

600

402-408.

601

(49) Yin, J.; Wang, C. Q.; Fang, C. Q.; Zhang, S.; Cao, Y. Z.; Li, K. B.; Leal, W. S.,

602

Functional characterization of odorant-binding proteins from the scarab beetle

603

Holotrichia oblita based on semiochemical-induced expression alteration and gene

604

silencing. Insect Biochemistry and Molecular Biology, 2019, 104, 11-19. 29

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

605

(50) Siciliano, P.; He, X. L.; Woodcock, C.; Pickett, J. A.; Field, L. M.; Birkett, M.

606

A.; Kalinova, B.; Gomulski, L. M.; Scolari, F.; Gasperi, G., Identification of

607

pheromone components and their binding affinity to the odorant binding protein

608

CcapOBP83a-2 of the Mediterranean fruit fly, Ceratitis capitata. Insect Biochemistry

609

& Molecular Biology 2014, 48, 51-62.

610

(51) Bohbot, J.; Vogt, R. G., Antennal expressed genes of the yellow fever mosquito

611

(Aedes aegypti L.); characterization of odorant-binding protein 10 and takeout. Insect

612

Biochemical and Molocular Biology 2005, 35, 961–979.

613

(52) Zhou, S.; Stone, E. A.; Mackay, T. F.; Anholt, R. R., Plasticity of the

614

chemoreceptor repertoire in Drosophila melanogaster. PLoS Genetics 2009, 5,

615

e1000681-235.

616

(53) Wang, X. Y.; Xiong, M.; Lei, C. L.; Zhu, F., The developmental transcriptome of

617

the synanthropic fly Chrysomya megacephala and insights into olfactory proteins.

618

BMC Genomics, 2015, 16, 20.

619

(54) Campanini, E. B.; Congrains, C.; Torres, F. R.; Brito, R. A. D., Odorant-binding

620

proteins expression patterns in recently diverged species of Anastrepha fruit flies.

621

Scientific Reports 2017, 7, 2194.

622

(55) Nakamura, A. M.; Chahad-Ehlers, S.; Lima, A. L.; Taniguti, C. H.; Jr, S. I.;

623

Torres, F. R.; de Brito, R. A., Reference genes for accessing differential expression

624

among developmental stages and analysis of differential expression of OBP genes in

625

Anastrepha obliqua. Scientific Reports 2015, 6, 17480.

626

(56) Kromann, S. H.; Saveer, A. M.; Binyameen, M.; Bengtsson, M.; Birgersson, G.; 30

ACS Paragon Plus Environment

Page 30 of 48

Page 31 of 48

Journal of Agricultural and Food Chemistry

627

Hansson, B. S.; Schlyter, F.; Witzgall, P.; Ignell, R.; Becher, P. G., Concurrent

628

modulation of neuronal and behavioural olfactory responses to sex and host plant cues

629

in a male moth. Proceedings of the Royal Society B: Biological Sciences 2015, 282,

630

20141884.

631

(57) Barrozo, R. B.; Gadenne, C.; Anton, S., Switching attraction to inhibition:

632

mating-induced reversed role of sex pheromone in an insect. Communicative &

633

Integrative Biology 2010, 213, 2933-2939.

634

(58) Barrozo, R. B.; Jarriault, D.; Simeone, X.; Gaertner, C.; Gadenne, C.; Anton, S.,

635

Mating-induced transient inhibition of responses to sex pheromone in a male moth is

636

not mediated by octopamine or serotonin. Journal of Experimental Biology 2010, 213,

637

1100.

638

(59) Saveer, A. M.; Kromann, S. H.; Birgersson, G.; Bengtsson, M.; Lindblom, T.;

639

Balkenius, A.; Hansson, B. S.; Witzgall, P.; Becher, P. G.; Ignell, R., Floral to green:

640

mating switches moth olfactory coding and preference. Proceedings of the Royal

641

Society B: Biological Sciences 2012, 279, 2314-2322.

642

(60) Jang, E. B., Effects of mating and accessory gland injections on

643

olfactory-mediated behavior in the female mediterranean fruit fly, Ceratitis capitata.

644

Journal of Insect Physiology 1995, 41, 705-710.

645

(61) Kamala Jayanthi, P. D.; Kempraj, V.; Aurade, R. M.; Venkataramanappa, R. K.;

646

Nandagopal, B.; Verghese, A.; Bruce, T. J., Specific volatile compounds from mango

647

elicit oviposition in gravid Bactrocera dorsalis females. Journal of Chemical Ecology

648

2014, 40, 259-266. 31

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

649

(61) Yin, J.; Wang, C. Q.; Fang, C. Q.; Zhang, S.; Cao, Y. Z.; Li, K. B.; Leal, W. S.,

650

Functional characterization of odorant-binding proteins from the scarab beetle

651

Holotrichia oblita based on semiochemical-induced expression alteration and gene

652

silencing. Insect Biochemistry and Molecular Biology, 2019, 104, 11-19.

32

ACS Paragon Plus Environment

Page 32 of 48

Page 33 of 48

Journal of Agricultural and Food Chemistry

653

FIGURE CAPTIONS

654

Figure 1. Alignment of amino acid sequences of OBPs from dipteran insects. The six

655

conserved cysteines are marked with a black star.

656

Figure 2. Alignment of amino acid sequences of BodoOBP1, BodoOBP2 and

657

AgamOPB1(AAL84179.1). Conserved amino acids in all OBPs are marked with

658

black triangles. Predicted signal peptides are indicated by bold underline.

659

Figure 3. 3D structure of BodoOBP1 and BodoOBP2. A: BodoOBP1. B: BodoOBP2.

660

3D structure analysis was performed with the program SWISS MODEL using the

661

structure of AgamOBP1 as a template.

662

Figure 4. Expression and purification of recombinant BodoOBP1 and BodoOBP2. All

663

proteins were expressed in E. coli BL21(DE3) using the vector pET28a(+) and

664

purified by Ni-NTA Resin. Lane 1 and Lane 15: Molecular weight markers. Lane 2:

665

supernatant of pET28a(+) (no insert). Lane 3: inclusion body of pET28a(+) (no

666

insert). Lane 4: supernatant of BodoOBP1 without induction. Lane 5: inclusion body

667

of BodoOBP1 without induction. Lane 6: supernatant of BodoOBP1after induction

668

with 0.4mM IPTG. Lane 7: inclusion body of BodoOBP1 after induction with 0.4mM

669

IPTG. Lane 8: purified BodoPBP1. Lane 9: supernatant of BodoOBP2 without

670

induction. Lane 10: inclusion body of BodoOBP2 without induction. Lane 11:

671

supernatant of BodoOBP2 after induction with 0.4mM IPTG. Lane 12: inclusion body

672

of BodoOBP2 after induction with 0.4mM IPTG. Lane 13: purified BodoPBP2.

673

Figure 5. Binding curves and Scatchard plots of 1-NPN to BodoOBP1 and 33

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

674

BodoOBP2.

675

Figure 6. Competitive binding curves of some ligands to BodoOBP1 and BodoOBP2.

676

A and B: BodoOBP1. C and D: BodoOBP2.

677

Figure 7. Expression patterns of BodoOBP1 and BodoOBP2 in different tissues. The

678

data and error bars represent the means and standard errors of three biological

679

replicates. Different letters above the bars indicate significant differences (P< 0.05,

680

Student–Newman–Keuls (SNK) in one way ANOVA). “*” indicates significant

681

differences, * p < 0.05; ** p < 0.01; *** p < 0.001 (Student's t test).

682

Figure 8. Expression patterns of BodoOBP1 and BodoOBP2 (relative to the 1st) in

683

different developmental stages. M/F-Pupae: Male or female pupae. M/F-Adults: Male

684

or female adults. The data and error bars represent the means and standard errors of

685

three biological replicates. Different letters above the bars indicate significant

686

differences (P< 0.05, Student–Newman–Keuls (SNK) in one way ANOVA). “*”

687

indicates significant differences, * p < 0.05; ** p < 0.01; *** p < 0.001 (Student's t

688

test).

689

Figure 9. Expression levels of BodoOBP1 and BodoOBP2 (relative to 24h after

690

eclosion in female) in pre- and post-mating adults. The data and error bars represent

691

the means and standard errors of three biological replicates. Different letters above the

692

bars indicate significant differences (P< 0.05, Student–Newman–Keuls (SNK) in one

693

way ANOVA).

694

Figure 10. Effect of RNAi treatment on the transcript levels of BodoOBP1 and 34

ACS Paragon Plus Environment

Page 34 of 48

Page 35 of 48

Journal of Agricultural and Food Chemistry

695

BodoOBP2. The data and error bars represent the means and standard errors of three

696

biological replicates.

697

Figure 11. The attraction indexes of B. odoriphaga (A) virgin female adults to male

698

adults and (B) dipropyl trisulfide to adults. The data and error bars represent the

699

means and standard errors of three biological replicates. “*” means significant

700

difference between dsOBP1 or dsOBP2 treated group and dsGFP treated group, * p