Towards an understanding of the molecular basis of effective RNAi against a global insect pest, the whitefly Bemisia tabaci

Accepted Manuscript Towards an understanding of the molecular basis of effective RNAi against a global insect pest, the whitefly Bemisia tabaci Yuan Luo, Qingguo Chen, Junbo Luan, Seung Ho Chung, Joyce Van Eck, R. Turgeon, Angela E. Douglas PII:

S0965-1748(17)30106-6

DOI:

10.1016/j.ibmb.2017.07.005

Reference:

IB 2974

To appear in:

Insect Biochemistry and Molecular Biology

Received Date: 25 June 2017 Revised Date:

15 July 2017

Accepted Date: 19 July 2017

Please cite this article as: Luo, Y., Chen, Q., Luan, J., Chung, S.H., Van Eck, J., Turgeon, R., Douglas, A.E., Towards an understanding of the molecular basis of effective RNAi against a global insect pest, the whitefly Bemisia tabaci, Insect Biochemistry and Molecular Biology (2017), doi: 10.1016/ j.ibmb.2017.07.005. This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. 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.

ACCEPTED MANUSCRIPT

fate of dsRNA in insect gut

TE D

M AN U

SC

RI PT

In planta RNAi

insect performance

AC C

EP

dsRNA-insect gut nuclease dsRNA-essential insect gene insect gut nuclease

. .

. .

ACCEPTED MANUSCRIPT

1

Towards an understanding of the molecular basis of effective RNAi against a global insect

2

pest, the whitefly Bemisia tabaci

3 Yuan Luo1,5, Qingguo Chen2, Junbo Luan1, Seung Ho Chung1, Joyce Van Eck3, R. Turgeon2,

5

Angela E. Douglas1,4*

RI PT

4

6 7

1

8

2

9

NY14853, USA

SC

Department of Entomology, Cornell University, Ithaca, NY14853, USA

M AN U

Plant Biology Section, School of Integrative Plant Science, Cornell University, Ithaca,

3

Boyce Thompson Institute, Ithaca, NY14853, USA

11

4

Department of Molecular Biology & Genetics, Cornell University, Ithaca, NY14853, USA

12

5

Current address: Department of Molecular & Integrative Physiology, University of Michigan,

13

Ann Arbor, MI 48105

*Corresponding author. Email address: [email protected]

EP

15

AC C

14

TE D

10

ACCEPTED MANUSCRIPT

Abstract

17

In planta RNAi against essential insect genes offers a promising route to control insect crop

18

pests, but is constrained for many insect groups, notably phloem sap-feeding hemipterans, by

19

poor RNAi efficacy. This study conducted on the phloem-feeding whitefly Bemisia tabaci

20

reared on tomato plants investigated the causes of low RNAi efficacy and routes to ameliorate

21

the problem. Experiments using tomato transgenic lines containing ds-GFP (green fluorescent

22

protein) revealed that full-length dsRNA is phloem-mobile, ingested by the insects, and degraded

23

in the insect. We identified B. tabaci homologs of nuclease genes (dsRNases) in other insects

24

that degrade dsRNA, and demonstrated that degradation of ds-GFP in B. tabaci is suppressed by

25

administration of dsRNA against these genes. dsRNA against the nuclease genes was co-

26

administered with dsRNA against two insect genes, an aquaporin AQP1 and sucrase SUC1, that

27

are predicted to protect B. tabaci against osmotic collapse. When dsRNA constructs for AQP1,

28

SUC1, dsRNase1 and dsRNase2 were stacked, insect mortality was significantly elevated to 50%

29

over 6 days on artificial diets. This effect was accompanied by significant reduction in gene

30

expression of the target genes in surviving diet-fed insects. This study offers proof-of-principle

31

that the efficacy of RNAi against insect pests can be enhanced by using dsRNA to suppress the

32

activity of RNAi-suppressing nuclease genes, especially where multiple genes with related

33

physiological function but different molecular function are targeted.

35

SC

M AN U

TE D

EP

AC C

34

RI PT

16

Keywords: Bemisia tabaci, dsRNA degradation, RNA interference, RNAi efficacy

ACCEPTED MANUSCRIPT

36 37

1. Introduction RNA interference (RNAi) holds great promise as a novel strategy against insect pests of agricultural crops (Scott et al., 2013; Whyard et al., 2009; Xue et al., 2012). This is because, in

39

principle, the only information needed to target an essential insect gene with exquisite specificity

40

is the insect gene sequence; and the double-stranded RNA molecules (dsRNA) can be delivered

41

via transgenesis of the crop plant. Specifically, the plant is engineered to express dsRNA against

42

the insect gene of interest. On ingestion by the insect, the dsRNA is internalized into cells,

43

where it is cleaved by an insect dsRNA-specific enzyme, Dicer-2, into small interfering RNA

44

(siRNA, ca. 21 nt) that guides the Argonaute protein of the RNA induced silencing complex

45

(RISC) to degrade complementary mRNAs (Wilson and Doudna, 2013; Zamore et al., 2000). In

46

planta RNAi has yielded significant plant protection against the western corn rootworm

47

Diabrotica virgifera virgifera (Baum et al., 2007), cotton bollworm Helicoverpa armigera (Mao

48

et al., 2007) and Colorado potato beetle Leptinotarsa decemlineata (Zhang et al., 2015), and

49

insecticidal RNAi in transgenic crops are reported to be near commercial release.

SC

M AN U

TE D

50

RI PT

38

Despite these substantial advances, many RNAi studies on insects have yielded moderate or variable knock-down of gene expression, with limited effects on insect phenotype and

52

performance (Belles, 2010; Scott et al., 2013; Terenius et al., 2011). These problems apply

53

particularly to plant sap-feeding insects, such as whiteflies, aphids, psyllids and planthoppers,

54

including major pests and vectors of plant viruses, where RNAi against essential genes often

55

reduces growth or reproduction, but has small or no effect on survivorship (Coleman et al., 2015;

56

Pitino et al., 2011; Thakur et al., 2014; Tzin et al., 2015; Zha et al., 2011).

AC C

EP

51

57

A possible cause of the limited efficacy of in planta RNAi against many insects is that the

58

dsRNA molecules ingested by the insects may be delivered poorly to the intracellular RNAi

ACCEPTED MANUSCRIPT

machinery in the insect cells. The plants are engineered to produce long dsRNA molecules, but

60

the dsRNA can be processed by the plant RNAi machinery to generate small RNAs (sRNAs, 21-

61

24 nt) (Baulcombe, 2015; Melnyk et al., 2011b). It is widely believed that sRNAs are less

62

effective than long dsRNA molecules as triggers of insect RNAi (Huvenne and Smagghe, 2010).

63

However, this view is based largely on studies of RNAi in cultured insect cells (Saleh et al.,

64

2006; Ulvila et al., 2006), and appears to be contradicted by a study on the phloem-feeding aphid

65

Myzus persicae that obtained similar gene expression reduction with dsRNA and diced dsRNA

66

(Mutti et al., 2006). Further possible impediments to dsRNA delivery are, first, that the plant

67

dsRNA can be degraded by non-specific nucleases in the insect saliva, gut lumen or hemolymph

68

(Allen and Walker, 2012; Christiaens et al., 2014; Garbutt et al., 2013; Liu et al., 2012b; Luo et

69

al., 2013; Spit et al., 2017; Wang et al., 2016b; Wynant et al., 2014), and, second, that the

70

intracellular RNAi machinery in the insect may have low activity (Garbutt and Reynolds, 2012).

M AN U

SC

RI PT

59

The purpose of this study was to identify the factors in the plant and insect that limit the

72

efficacy of RNAi against phloem-feeding insects, and to use this information for improved

73

design of RNAi. Our experiments were conducted on the whitefly Bemisia tabaci, which is a

74

globally-important pest of many crops (Navas-Castillo et al., 2011). As in various other insects

75

(see above), multiple RNAi studies of B. tabaci have yielded incomplete knock-down of gene

76

expression (Ghanim et al., 2007; Luan et al., 2013; Raza et al., 2016; Thakur et al., 2014;

77

Upadhyay et al., 2011), constraining the application of RNAi as a control strategy for this pest.

78

We hypothesized that whitefly RNAi may be limited by processing of dsRNA to small RNAs

79

(sRNAs) in the plant and, additionally or alternatively, by nonspecific degradation of dsRNA in

80

the insect. We tested these hypotheses by following the fate of dsRNA constructed against a 370

81

nt fragment of the green fluorescent protein (GFP) gene of the jellyfish Aequorea victoria,

AC C

EP

TE D

71

ACCEPTED MANUSCRIPT

administered to the insects via artificial diets and plant transgenesis. These experiments led us to

83

identify non-specific degradation of dsRNA by B. tabaci, which we then reduced by RNAi

84

against two candidate B. tabaci nuclease genes. In our final experiments, we tested the efficacy

85

of stacking RNAi against the nuclease genes with RNAi against candidate essential genes of B.

86

tabaci. Our genes of choice were an aquaporin and a glucohydrolase of family GH-13,

87

previously identified as candidate osmoregulation genes that protect the insect from rapid

88

dehydration and death (Jing et al., 2016; Mathew et al., 2011). Following the nomenclature of

89

Jing et al. (2016), we call these genes BtAQP1 and BtSUC1, respectively. This study yields new

90

insights into the fate of dsRNA administered to phloem-feeding insects by in planta RNAi, and

91

the value of this technology as a novel insect pest control strategy.

M AN U

SC

RI PT

82

92 2. Material and methods

94

2.1 Plants and Insects

95

TE D

93

Tomato plants (Solanum lycopersicum cv. Florida Lanai) (seeds provided by C. L. McKenzie, USDA-ARS Fort Pierce, FL, USA) were grown in compost supplemented with Miracle-Gro®

97

Water Soluble All Purpose Plant Food in climate-controlled chambers at 25±2 °C with a

98

14L:10D light cycle at 400 µmol m-2 s-1 PAR. The Bemisia tabaci MEAM1 culture (mtCO1 GenBank accession no. KM507785) was derived

AC C

99

EP

96

100

from a collection from poinsettia (Euphorbia pulcherrima Willd. Ex Klotzsch) in Ithaca, NY,

101

USA in 1989. The insects were maintained on 5-6-week-old tomato plants at 25±2 °C with a

102

14L:10D light cycle at 400 µmol m-2 s-1 PAR. One-day-old adult males and female insects were

103

caged to test plant leaves using BugDorm-1 Insect Cages (Bio Quip, Rancho Dominguez, CA) or

104

custom clip-cages, and fed on a sterile artificial liquid diet (Douglas et al., 2001) containing 0.5

ACCEPTED MANUSCRIPT

105

M sucrose and 0.15 M amino acids in Parafilm sachets. Some experiments used isolated guts,

106

dissected with fine pins from adult female and male insects into phosphate-buffered saline.

107

109

2.2 RNA extractions

RI PT

108

For total RNA isolation, leaves or whiteflies were homogenized in TRIzol® Reagent (Cat# 15596-026, Thermo Fisher Scientific, Waltham, USA) with 1:1 (vol) Lysing Matrix D Bulk

111

beads (Cat# 116540434, MP Biomedicals, Santa Ana, USA) on a MP FastPrep-24™

112

homogenizer (Cat# 116004500, MP Biomedicals, Santa Ana, USA) with 5.5 M/S for 2 x 30

113

seconds. The lysed samples were centrifuged at 13,000 rpm for 3 min to remove the beads and

114

RNA was isolated from the supernatant following the manufacturer’s instructions. LMW (low

115

molecular weight) and HMW (high molecular weight) RNA isolation employed the same

116

protocol as above, except that the supernatant of the lysate homogenate was mixed with 70%

117

ethanol (1:1, by vol) and passed through a RNeasy mini column (Cat# 74104, Qiagen, Venlo,

118

Limburg, USA). RNAs longer than 200 nt were collected in the column and washed from the

119

column following the RNeasy mini kit protocol. RNAs shorter than 200 nt were precipitated

120

from the flowthrough with 100% isopropanol (1:1, by vol) and 3 M sodium acetate (pH 5.2,

121

10:1, by vol), followed by 75% ethanol wash and then dissolved in nuclease-free water. RNA

122

samples used for qRT-PCR were further processed to remove genomic DNA by incubation with

123

10 µl DNase 1, using the reagents and protocol in the RNase-free DNase set (Cat# 79254,

124

Qiagen, Valencia, USA). To recover dsRNA from diets, samples of the diet were combined with

125

isopropanol (1:1, by vol) and 3 M sodium acetate pH 5.2 (1:10 by vol). The concentration of

126

RNA was determined spectrophotometrically with a Nanodrop-2000 (Thermo Fisher Scientific,

AC C

EP

TE D

M AN U

SC

110

ACCEPTED MANUSCRIPT

127

Waltham, USA), and RNA integrity was verified by denaturing formaldehyde gel-

128

electrophoresis.

129

131

2.3 Synthesis of cDNA and dsRNA

RI PT

130

cDNA libraries for amplification of whitefly genes were prepared with oligo (dT) primers using the Reverse Transcription System (Cat# A3500, Promega, Madison, USA) following the

133

manufacturer’s instructions. Candidate osmoregulatory genes, comprising aquaporin AQP1

134

(NCBI Accession KX390870) (Jing et al., 2016; Mathew et al., 2011), sucrase SUC1 (also

135

known as α-glucohydrolase GH13-1 (Jing et al., 2016)) (NCBI Accession KX390872), and

136

candidate nuclease genes dsRNase1 (NCBI Accession KX390873) and dsRNase2 (NCBI

137

Accession KX390874) identified in this study, were amplified from the cDNA template in a

138

reaction mix containing 0.4 µM primers (Table S1), 2 U Invitrogen Platinum Taq DNA

139

polymerase (Cat# 10966018, Thermo Fisher Scientific, Waltham, USA), 1.5 mM MgCl2 and 200

140

ng cDNA template in 25 µl volume, using a Techne thermal cycler. The thermal profile

141

comprised 2 min at 94°C for initial denaturation, 30 cycles with 95°C for 30 sec, 55°C for 30

142

sec, 72°C for 1 min (depending on products length, 1 kb min-1) and final extension cycle of 72°C

143

for 5 min. Amplicon sequences were verified by Sanger sequencing, then introduced into

144

PGEM-T

145

competent cells. The plasmid was extracted and, following confirmation by sequencing, used as

146

template for dsRNA synthesis with primers listed in Table S1.

M AN U

TE D

EP

vector (Cat# A1360, Promega, Madison, USA) and transformed into DH5α™

AC C

147

SC

132

The dsRNA was synthesized using the AmpliScribe TM T7-Flash Transcription Kit (Cat#

148

ASF3257, Epicentre Biotechnologies, Madison, USA), according to the manufacturer’s

149

instructions. The templates were the pGFP2 plasmid for ds-GFP, and plasmids obtained for the

ACCEPTED MANUSCRIPT

150

whitefly genes identified above. The dsRNA product was quantified by Nanodrop, and run on a

151

gel with 1 kb plus molecular weight ladder (Cat# 10787026, Thermo Fisher Scientific, Waltham,

152

USA) to confirm the predicted size.

154

RI PT

153 2.4 Construction of dsRNA expression cassette

The dsRNA expression cassette was constructed in the pHANNIBAL vector. A 370 nt GFP

156

sequence was amplified and inserted in inverted orientation into pHANNIBAL using the PDK

157

intron as a spacer and different restriction sites, XhoI and EcoRI for sense GFP and HindIII and

158

XbaI for the inverted GFP sequence. Two phloem-companion cell-specific promoters

159

Galactinol Synthase from melon Cucumis melo, CmGAS (Ayre et al., 2003; Haritatos et al.,

160

2000) and sucrose-H+ symporter from Arabidopsis thaliana, AtSUC2 (Dunoyer et al., 2007;

161

Truernit and Sauer, 1995) were cloned and inserted into the binary vector pER8 using XhoI/SpeI

162

and XhoI restriction sites separately. The GFP-intron-rGFP cassette was assembled using the

163

Gibson assembly kit (New England Biolabs Cat# E2611S) following manufacturer’s

164

instructions, and inserted downstream of the promoter in pER8 vector using a XhoI restriction

165

site for AtSUC2 and SpeI site for CmGAS1.

168

M AN U

TE D

EP

167

2.5 Generation of tomato transgenics with dsRNA gene construct

AC C

166

SC

155

The binary vector pER8 was introduced via electroporation into Agrobacterium tumefaciens

169

strain LBA4404. Kanamycin selection at 500 µg ml-1 was used to select for transformants.

170

Transformation of the tomato plants was performed according to methods described by Van Eck

171

et al. (2006). Putative transgenic plants were transferred to soil and maintained in an incubator

172

at 25±2°C with a relative humidity of 60 – 70% and with a 14L:10D light cycle at 400 µmol m-2

ACCEPTED MANUSCRIPT

s-1 PAR. Total DNA was extracted from leaves of 5-6-week-old transgenic plants and verified

174

for transformation by PCR with sequence specific primers (Table S1). Eight transgenic lines

175

were confirmed for each promoter.

176 177 178

2.6 Administration of dsRNA to whiteflies

RI PT

173

dsRNA was administered to adult insects either via artificial diets (at 0.1-1 µg µl-1, varying with experiment with full details provided in the relevant figure legends) or via transgenic

180

tomato lines. For analysis of the fate of ds-GFP, 100 whiteflies were administered to each diet

181

cage and ca. 250 whiteflies were caged to each plant. For insect performance experiments on

182

diets, 10 replicate groups of 40 adult whiteflies (one day post-emergence) were applied to each

183

diet treatment and mortality was monitored daily over 6 days, with insects transferred to fresh

184

diet containing dsRNA every two days. At the end of each experiment, all live insects were

185

transferred to 500 µl TRIzol® Reagent (Cat# 15596-026, Thermo Fisher Scientific, Waltham,

186

USA) and stored at -80 °C prior to isolation of total RNA (as above).

187

189

2.7 Northern blots

EP

188

TE D

M AN U

SC

179

RNA extracted from whiteflies and plants was separated on denatured polyacrylamide-urea gels (SequaGel - UreaGel System, cat# EC-833, National Diagnostics, Atlanta, GA, USA)

191

containing 8% monomers. The gel was pre-run at 250 V in 0.5 X TBE buffer for 30 min, then

192

20 µg sample RNA was combined with an equal volume of Gel Loading Buffer II (Cat#

193

AM8547, Thermo Fisher Scientific, Waltham, USA) and heated at 95 °C for 4 minutes to

194

denature the RNA. Samples were loaded in urea-cleaned wells and run at 250 V in 0.5 X TBE

195

buffer until the loading dye migrated to the far end of the gel. Uniform sample loading was

AC C

190

ACCEPTED MANUSCRIPT

confirmed by staining of 5S rRNA with SYBR Gold, followed by transfer to Hybond-NX

197

membrane (Cat# RPN 203T, GE Healthcare, Wilkes-Barre, PA) with the Owl™ HEP Series

198

Semidry Electroblotting Systems at 0.4 A for 1 hour. The transferred RNA was cross-linked

199

using a UV crosslinker at 120 kJ for 30 sec for HMW RNA and by a published procedure (Pall

200

and Hamilton, 2008) for LMW RNA.

201

RI PT

196

The 370 nt GFP antisense probe used for northern blotting was generated using the

MAXIscript® In Vitro Transcription Kit (Cat# AM1308, Thermo Fisher Scientific, Waltham,

203

USA) with 3.125 µM alpha-32P UTP (10 mCi ml-1, 800 Ci mmol-1, Perkin Elmer, Waltham,

204

USA), 5 µM UTP and 100 µM ADP. Unincorporated 32P-label was removed using a RNeasy

205

mini column (Cat# 74104, Qiagen, Venlo, Limburg, USA). The 32P-labeled probe was brought

206

to 100 µl with nuclease-free water, mixed with 350 µl buffer RLT (RNeasy mini kit, Cat# 74104,

207

Qiagen, Venlo, Limburg, USA) and 250 µl 100% ethanol, collected onto a RNeasy mini column,

208

washed twice with buffer RPE (from RNeasy mini kit, Cat# 74104, Qiagen, Venlo, Limburg,

209

USA), and eluted with 40 µl nuclease-free water. For each membrane, half of the probe was

210

heated at 95 °C for 4 min and immediately added to the pre-hybridized membrane at 2 X 106 cpm

211

ml-1 final concentration. For HMW RNA detection, pre-hybridization and hybridization were

212

conducted at 68 °C using Ambion Ultrahyb buffer (CAT# AM8670, Thermo Fisher Scientific,

213

Waltham, USA) according to the manufacturer’s protocol. The blot was washed twice in 2X

214

SSC, 0.1% SDS buffer at 68°C, followed by two washes in 0.1X SSC, 0.1% SDS buffer at the

215

same temperature. For LMW RNA detection, the blot was pre-hybridized in hybridization buffer

216

(5 X SSC, 20 mM Na2HPO4 (pH 7.2), 7% SDS, 2 X Derhardt’s solution) at 50 °C for at least two

217

hours, then hybridized with a final concentration of 2.5 x 106 cpm ml-1 probe in the same buffer

218

at 50 °C overnight. The membrane was washed four times at 50 °C in non-stringent wash buffer

AC C

EP

TE D

M AN U

SC

202

ACCEPTED MANUSCRIPT

(3X SSC, 25 mM NaH2PO4 pH 7.5, 5% SDS) and once in stringent wash buffer (1X SSC, 0.1%

220

SDS), and then exposed for autoradiography. The signal was collected on phosphor screen

221

(Molecular Dynamics) and scanned using a Typhoon 9400 fluorescent imager.

222

ImageJ was used for ds-GFP band density analyses.

223 224

2.8 qRT-PCR

RI PT

219

To quantify the expression of target whitefly genes, qRT-PCR was performed with RNA

226

extracted from three biological replicates of whiteflies. cDNA was prepared using random

227

primers of High-Capacity cDNA Reverse Transcription Kit or SuperScript™ II Reverse

228

Transcriptase (Cat# 4368814 and 18064014, Thermo Fisher Scientific, Waltham, USA)

229

following the manufacturer’s instructions. For qRT-PCR, the 20 µl reaction mix comprised 10

230

µl Master Mix (Bio-Rad, Hercules, CA) or Power SYBR Green PCR Master Mix (Applied

231

Biosystems, Carlberg, CA, USA), precisely 1 µl cDNA template and 0.5-2 µl 10 µM primers

232

qRT-PCR primers (Table S1) designed with Primer Premier 5.0 software (Premier Biosoft

233

International, Palo Alto, CA). Amplifications were conducted in a C1000TM Thermal cycler

234

(Bio-Rad, Hercules, CA) with the following thermal profile: 95°C for 5 min, 40 amplification

235

cycles of 95°C for 15 sec, 55°C for 30 sec, and dissociation cycle of 95°C for 15 sec, 55°C for

236

15 sec then brought back to 95°C. Dissociation curves confirmed single peaks of the predicted

237

size without primer dimerization. All assays included three technical replicates with template-

238

free and non-RT as controls; and the relative expression was calculated using the 2- △△ Ct method

239

(Livak and Schmittgen, 2001), normalized to the whitefly 60S ribosomal protein L13a (RPL13)

240

gene. Mean Ct value of three technical replicates was calculated per sample.

241

AC C

EP

TE D

M AN U

SC

225

ACCEPTED MANUSCRIPT

242 243

2.9 Identification and phylogenetic analysis of candidate whitefly dsRNase genes To obtain an initial set of candidate nucleases in B. tabaci, the translated sequence of the nonspecific nuclease Bombyx mori, which has been demonstrated to degrade dsRNA and suppress

245

RNAi (Arimatsu et al., 2007; Liu et al., 2012b), was BLASTed (E value < 1.0 e-10) against the

246

translated RefSeq genes in the B. tabaci genome (Chen et al., 2016). The resultant B. tabaci

247

genes were analyzed by SignalP (Petersen et al., 2011) and NCBI conserved domain database for

248

signal peptide and conserved domain, respectively. For analysis of gene phylogenies, a

249

neighboring-joining tree was constructed of the conserved DNA/RNA non-specific nuclease

250

domain amino acid sequences (MEGA 6.06).

M AN U

251 252 253

SC

RI PT

244

2.10 Tomato plant grafting

For each graft, a scion with few mature leaves from 5-6-week-old wild-type tomato was cut and inserted to a transgenic dsRNA plant of the same age. The graft site was fastened with

255

Parafilm; and the plant was covered with a plastic bag and kept in darkness at 25°C for 48 h.

256

Lighting was gradually increased over the next 3-4 days, and the bag was removed at day-7.

257

Wildtype scions were grafted to two GAS:ds-GFP and two SUC2:ds-GFP plants. The plants

258

were harvested 3 weeks after grafting for RNA isolation and northern blotting analysis of ds-

259

GFP fragments.

261 262 263

EP

AC C

260

TE D

254

2.11 Statistical analysis

All data sets conformed to the expectations of normality by the Anderson Darling test and homogeneity of variance by the Levine and Bartlett tests. They were analyzed by one-way

ACCEPTED MANUSCRIPT

264

ANOVA with Fisher’s LSD post hoc test. Statistical analyses were conducted with JMP

265

software (SAS Institute, Miami, USA) and Minitab 17.

267

3. Results

268

3.1 Delivery of plant dsRNA to whiteflies

269

RI PT

266

The first experiments investigated the fate of dsRNA expressed under companion cell-specific promoters CmGAS from melon, Cucumis melo (expressed in companion cells of minor veins of

271

the leaf) and AtSUC2 from Arabidopsis thaliana (expressed in all companion cells) in transgenic

272

lines. The 370 bp ds-GFP construct (Fig. 1A) used in these experiments had no sequence

273

homology to tomato or whitefly genes, to ensure that processing of the dsRNA would not be

274

altered by sequence-specific effects of the dsRNA on gene expression in the plant or insect.

M AN U

275

SC

270

Northern blotting of ds-GFP fragments expressed under the CmGAS and AtSUC2 promoters in 5-6-week-old transgenic tomato lines yielded a band at ca. 21-25 nt in all lines and at 370 nt in

277

most lines (Fig 1B). These results are compatible with the interpretation that plant Dicer

278

enzyme(s) cleaved the full-length ds-GFP into sRNA(s). (Our methods could not discriminate

279

whether the sRNA product was one or multiple molecules in the range 21-25 nt.) Smearing of

280

the signal down the gel, indicative of non-specific degradation of ds-GFP, was very limited in

281

northern blots of plant RNA (Fig. S1A and B, whole gel picture of Fig 1B).

EP

AC C

282

TE D

276

The northern blots for the whiteflies that had fed on the ds-GFP plants bore either two or three

283

bands in the region of the predicted 370 nt ds-GFP band (Fig. S1C). The middle of these three

284

bands was the closest match to the 370 nt ds-GFP probe (Fig. S1C and Fig. 1C) and was detected

285

in the whiteflies feeding from ds-GFP plants but not the empty vector control plant. The other

286

two bands were evident in all the samples including the empty vector control, indicative of non-

ACCEPTED MANUSCRIPT

specific binding of whitefly RNA to the GFP probe. These results suggest that the whiteflies

288

feeding on the ds-GFP plants are ingesting ds-GFP. The blots for the whiteflies bore smeared

289

signal over an extended range of the gel, indicative of non-specific degradation, and no

290

detectable sRNA signal. (Fig. S1C&S1D, whole gel picture of Fig 1C).

RI PT

287

We postulated two alternative explanations for the apparent absence of GFP sRNA in the

292

whiteflies: that the sRNA, first, is not phloem-mobile and consequently not ingested by the

293

whiteflies; and, second, it is ingested but degraded in the whitefly. To test for the phloem

294

mobility of the full-length and sRNA, control scions were grafted onto transgenic plants

295

expressing ds-GFP under CmGAS or AtSUC2 promoters, with homografts onto wild-type (WT)

296

plants as negative control. Three weeks after grafting, entire scion apices were excised and

297

processed for GFP-RNA by northern blotting. The 370-nt band but not the sRNA band, was

298

detected (Fig. 1D and Fig. S1E). These data indicate that the processed sRNA derived from the

299

ds-GFP is not phloem-mobile, and therefore not available to the whiteflies feeding on phloem

300

sap from the transgenic plants.

303

M AN U

TE D

302

3.2 dsRNA degradation in the whiteflies

EP

301

SC

291

To investigate the fate of dsRNA molecules in the whiteflies further, adult insects were administered ds-GFP via artificial diet. For these experiments, we used the LMW hybridization

305

protocol to visualize the total RNA from diet and whitefly because this protocol enabled us to

306

detect LMW RNAs at low abundance as well as HMW RNA. The 370 nt ds-GFP was recovered

307

in northern blots of both the diet on which the whiteflies had fed and diet without whiteflies (Fig.

308

2A). Parallel analysis of the whiteflies yielded the 370 nt GFP band (as in the diet), a smear of

309

signal along the length of the gel, indicative of non-specific degradation, and a sRNA band at ca.

AC C

304

ACCEPTED MANUSCRIPT

21 nt (Fig. 2B). These data suggest that extra-oral degradation of dsRNA, as mediated by

311

salivary secretions in Lygus bugs (Allen and Walker, 2012), is not substantial in the whitefly, but

312

that an appreciable portion of ds-GFP ingested by the whiteflies is subjected to non-specific

313

degradation within the insect body. We hypothesized that dsRNA ingested by whiteflies is

314

subjected to non-specific degradation in the gut lumen, restricting the availability of dsRNA

315

molecules for uptake by cells of the gut epithelium and intracellular processing by the RNAi

316

machinery. To test this hypothesis, we applied phylogenetic methods to identify candidate

317

nuclease genes in the B. tabaci genome.

319 320

M AN U

318

SC

RI PT

310

3.3 Phylogenetic analysis of candidate dsRNase genes in the Bemisia tabaci genome Our strategy to identify candidate nuclease(s) in the B. tabaci genome was to identify orthologs of a Bombyx mori DNA/RNA non-specific nuclease gene (BmdsRNase) that has been validated

322

experimentally to cleave dsRNA and reduce the efficacy of RNAi (Liu et al., 2012b). The B.

323

mori dsRNase was BLASTed (E value < 1.0 e-10) against transcriptome databases of the whole

324

body, salivary gland and gut of B. tabaci, yielding three B. tabaci sequences (dsRNase1,

325

dsRNase2, dsRNase3 ) with a single DNA/RNA non-specific nuclease domain (NCBI conserved

326

domain database) and a predicted signal peptide (SignalP). A phylogenetic tree constructed by

327

the neighbor-joining method using the amino acid sequence of the conserved DNA/RNA non-

328

specific nuclease domain from multiple insect species aligned the B. tabaci dsRNase1 with aphid

329

nucleases with moderate bootstrap support, and dsRNase2 and dsRNase3 with nuclease-1 of

330

Tribolium castaneum, with excellent bootstrap support (Fig. 3A).

331 332

AC C

EP

TE D

321

We amplified part of the predicted full-length cDNA sequences of dsRNase1 and dsRNase2 from whole body and gut cDNA libraries of adult B. tabaci, but failed to amplify dsRNase3.

ACCEPTED MANUSCRIPT

Complementary searches of the B. tabaci whole body, gut and salivary gland transcriptome

334

databases (Su et al., 2012; Wang et al., 2010; Wang et al., 2012) yielded a few dsRNase3 reads

335

only in the salivary gland transcriptome, suggesting that this gene is weakly expressed in the

336

adult insects. Validating these transcriptome data, qRT-PCR analysis of the B. tabaci used in

337

this study confirmed that dsRNase1 and dsRNase2 are expressed, with three-fold enrichment of

338

dsRNase2 expression in the gut relative to the whole body (Fig. 3B).

RI PT

333

341

3.4 Effect of RNAi against dsRNase genes on ds-GFP ingested by the whiteflies We next asked whether inhibiting whitefly dsRNase genes could protect ds-GFP from

M AN U

340

SC

339

nonspecific degradation and thus improve RNAi. Enzymatically synthesized dsRNA against

343

dsRNase1 or dsRNase2 was fed to the whiteflies either synchronously or 3 days prior to adding

344

ds-GFP, with dsRNase-free treatment as the control (Fig 4A). qRT-PCR showed that

345

administration of dsRNA against dsRNase1 and dsRNase2 reduced their expression by 25-30%

346

(Fig. S2). Furthermore, northern blots from three independent experiments revealed that the

347

intensity of the full-length ds-GFP band and ca. 21 nt GFP was elevated in whiteflies

348

administered dsRNA against the dsRNase genes (Fig. S3), and this effect was significant for the

349

whiteflies pretreated with dsRNA against dsRNase2 and both dsRNases, but not dsRNase1(Fig.

350

4B). Because the RNAi-mediated reduction in expression of both dsRNase genes had the

351

greatest protective effect, our subsequent experiments used dsRNA against both genes.

353 354 355

EP

AC C

352

TE D

342

3.5 Effect of ds-dsRNase on efficacy of RNAi against whitefly osmoregulation genes The demonstration (Fig. 4) that RNAi against the whitefly dsRNase genes is protective for dsGFP provided the basis to investigate whether the efficacy of RNAi against whitefly genes of

ACCEPTED MANUSCRIPT

interest can be improved by combination with RNAi against the dsRNase genes. Our

357

experiments focused on whitefly genes AQP1 and SUC1, with the predicted function of

358

protecting the insect against osmotic dysfunction (Jing et al., 2016; Mathew et al., 2011).

359

Because these genes have been identified principally by bioinformatics methods, we first

360

conducted qRT-PCR experiments that confirmed the expression of these genes, including their

361

enriched expression in the gut (by 6-fold for SUC1 and 2.3-fold for AQP1), in the adult

362

whiteflies used for our experiments (Fig. S4).

SC

363

RI PT

356

We delivered the dsRNAs to adult whiteflies via artificial diet over a time-course of 6 days. The experimental treatments were dsRNA against the whitefly AQP1 and SUC1, either

365

separately or in combination, and with or without the dsRNA against dsRNase1 and dsRNase2.

366

Surviving insects on day-6 were used for gene expression analysis by qRT-PCR. The ds-

367

dsRNase treatments reduced expression of the cognate dsRNase genes by 30-35% (Fig. S5) but,

368

when administered as the sole dsRNA, did not affect expression of the target osmoregulation

369

genes (Fig. 5A & B). Expression of AQP1 and SUC1 was reduced in whiteflies feeding from

370

diets containing ds-AQP1 and dsSUC1, respectively, relative to the controls (diets with ds-GFP

371

or ds-dsRNase1&2, and dsRNA-free diets) and these effects were statistically significant apart

372

from ds-SUC1 without ds-dsRNase1&2. Inclusion of ds-RNase1&2 in the diet reduced mean

373

expression of AQP1 and SUC1 relative to diets without ds-RNase1&2 but the differences were

374

not statistically significant. However, this comparison may not be a reliable indicator of the

375

efficacy of suppression dsRNase1&2 expression. Specifically, AQP1 and SUC1 are expressed in

376

whitefly organs other than the gut, but the effect of ds-RNase1&2 on the RNAi-suppression of

377

their expression may be greater in the gut than other organs. This pattern could result from

378

limited systemic spread of RNAi (a likely constraint on the efficacy of RNAi in various insects

AC C

EP

TE D

M AN U

364

ACCEPTED MANUSCRIPT

379

(Scott et al., 2013)) with the consequence that our quantification of transcript abundance in the

380

whole insect bodies may have underestimated the impact of ds-RNase1&2 on expression of

381

AQP1 and SUC1 in the whitefly gut. The whiteflies in the two treatments containing both ds-AQP1 and ds-dsRNase1&2 displayed

383

elevated mortality, relative to the control diets (ds-GFP and dsRNA-free). In the absence of ds-

384

AQP1, ds-SUC1 had no discernible effect on whitefly mortality, but functioned synergistically

385

with ds-AQP1 in the presence of ds-dsRNase1&2 to yield mortality approaching 50% and

386

significantly greater than all other treatments.

388 389

M AN U

387

SC

RI PT

382

4. Discussion

The development of effective RNAi against plant sap-feeding insects cannot be achieved by trial-and-error. It requires a rational understanding of the fate of dsRNA molecules, in both the

391

plant tissues and in the feeding insects. This study demonstrates the complexity of these

392

processes, and identifies a strategy that improves the efficacy of RNAi against the globally-

393

important insect pest, the whitefly B.tabaci.

A critical factor determining the efficacy of in planta RNAi against phloem-feeding insects is

EP

394

TE D

390

the molecular identity of the phloem-mobile dsRNA molecules. It is inherently difficult to

396

identify phloem-mobile compounds with certainty. In particular, some widely-used methods for

397

direct sampling of the phloem contents of the plant sieve elements (incisions and EDTA

398

exudation) are now recognized to produce material that is contaminated by the contents of other

399

plant cells (Liu et al., 2012a; Zhang et al., 2012). However, reliable information about phloem-

400

mobile RNA molecules can be obtained from the RNA content of both wild-type tissue grafted

401

onto transgenic plants and phloem-feeding insects; and both of these approaches have been used

AC C

395

ACCEPTED MANUSCRIPT

previously (Dunoyer et al., 2010; Melnyk et al., 2011a; Yoo et al., 2004). When these two

403

methods were applied in this study to stable transgenic tomato containing ds-GFP under two

404

alternative phloem-specific promoters, they yielded the full-length 370 nt dsRNA, but not sRNA

405

in the 20-25 nt range. This result cannot be attributed to technical difficulties in detecting sRNA

406

because both full-length dsRNA and sRNA were detected in bulk transgenic leaf samples of the

407

same total RNA content. Detailed interpretation of this result requires consideration of the

408

structural relationships among plant cells. Specifically, the ds-GFP synthesis is restricted to

409

phloem companion cells, with predicted delivery via large diameter plasmodesmata to the sieve

410

elements and via more sparse and narrow-diameter plasmodesmata to non-vascular plant cells

411

that have little plasmodesmatal contact with the sieve elements (Kempers et al., 1993). Thus, the

412

recovery of dsRNA, but not sRNA, in the sieve elements can, in principle, be attributed to

413

multiple processes: dsRNA processing in the non-vascular tissue and not the companion cells; a

414

selective barrier reducing the delivery of sRNA (but not dsRNA) from the companion cells to the

415

sieve elements; degradation of sRNA in the sieve elements; and the efficient re-uptake of sRNA

416

from sieve elements back into companion cells. Further research is required to establish the

417

relative importance of these different processes. Interestingly, this study adds to the diversity of

418

dsRNA products reported in phloem sap, including short single-stranded RNA in cucurbits and a

419

legume (castor bean) (Yoo et al., 2004) and 21 nt or 24 nt siRNA in Arabidopsis (Dunoyer et al.,

420

2010; Melnyk et al., 2011a). Although the basis for these different results is uncertain because

421

the various studies use both different methods and plant systems, one possible interpretation is

422

that the phloem-mobile compounds vary with plant species. Generally, RNAi in insects is more

423

effective with long dsRNA molecules than sRNAs (either single- or double-stranded) (see

AC C

EP

TE D

M AN U

SC

RI PT

402

ACCEPTED MANUSCRIPT

Introduction), suggesting that the tomato transgenic lines containing dsRNA under phloem-

425

specific promoters, as used in this study, have the potential for effective RNAi against whiteflies.

426

Further insight into the delivery of in planta RNAi to phloem-feeding insects came from our

427

comparison of the ds-GFP products in whiteflies feeding from ds-GFP transgenic plants and ds-

428

GFP-supplemented artificial diet. Specifically, sRNA was detected in insects fed on ds-GFP via

429

the diet but not the plant. A parsimonious interpretation of this discrepancy is that non-specific

430

nuclease activity in the insect gut degraded the dsRNA in the phloem sap, quantitatively

431

preventing dsRNA delivery to gut cells, where siRNA is generated by the cytoplasmic RNAi

432

machinery; but that the high concentrations of dsRNA in the artificial diets saturated the gut

433

nuclease activity, such that a proportion of the ingested dsRNA was translocated to gut cells,

434

yielding detectable sRNA in the insect.

M AN U

SC

RI PT

424

It is becoming increasingly evident from this and other studies that non-specific nucleases in

436

the gut lumen or saliva are a key factor limiting the efficacy of orally-delivered dsRNA in insects

437

(Allen and Walker, 2012; Arimatsu et al., 2007; Liu et al., 2012b; Spit et al., 2017; Wang et al.,

438

2016b; Wynant et al., 2014). The cognate nuclease gene family is conserved across insects and

439

some crustaceans, and has also been reported in one viral genome (the shrimp white spot

440

syndrome virus in the Nimaviridae) but lacks homologs in other animals, raising the possibility

441

that it was acquired horizontally in the common ancestor of insects and crustaceans (Wynant et

442

al., 2014). Although these genes likely have an anti-viral function, as indicated by the increased

443

expression of the dsRNase gene in Bombyx mori infected with cytoplasmic polyhedrosis virus

444

(Wu et al., 2009), they are also strong candidates as suppressors of orally-delivered RNAi.

445

Consistent with this interpretation, hemolymph injections of dsRNA against the dsRNase genes

AC C

EP

TE D

435

ACCEPTED MANUSCRIPT

446

resulted in reduced transcript abundance of dsRNase genes and dsRNase enzymatic activity in

447

the midgut of in the locust Schistocerca gregaria (Wynant et al., 2014).

448

This study extends the application of RNAi in whiteflies to test whether RNAi-mediated suppression of dsRNase genes results in enhanced efficacy of RNAi against other insect genes,

450

specifically the two predicted osmoregulation genes, AQP1 and SUC1, identified previously by

451

Jing et al. (2016) and Mathew et al. (2011). As predicted, orally-delivered ds-dsRNase1&2 both

452

protected ds-GFP from non-specific degradation and increased the efficacy of RNAi against

453

osmoregulation genes, as quantified by gene expression and survivorship of insects administered

454

dsRNA via artificial diet. Interestingly, relatively modest suppression of expression of the

455

dsRNase genes resulted in substantial protection of orally-delivered dsGFP (Fig. 4). We

456

hypothesize that this effect may result from a very short half-life of the nuclease proteins coded

457

by the dsRNase genes (i.e. these proteins have a high turnover rate), resulting in a large reduction

458

in nuclease activity for relatively small changes in transcript abundance. More generally, our

459

results are paralleled by the recent demonstration that the efficacy of RNAi against the Colorado

460

potato beetle Leptinotarsa decemlineata is also enhanced by co-administration of dsRNA against

461

gut nuclease genes (Spit et al., 2017). These data suggest that the efficacy of RNAi in various

462

insects (Scott et al., 2013; Terenius et al., 2011; Wang et al., 2016b) may be significantly

463

improved by RNAi-mediated suppression of dsRNase genes, and potentially other RNAi-

464

suppressors.

SC

M AN U

TE D

EP

AC C

465

RI PT

449

It has been argued that the efficacy and durability of RNAi for insect pest control can be

466

enhanced by stacking dsRNA constructs against multiple targets in the pest species (Scott et al.,

467

2013; Tzin et al., 2015). This strategy is likely particularly valuable where the target genes have

468

unrelated molecular function but linked physiological function because the incomplete

ACCEPTED MANUSCRIPT

suppression of multiple contributing gene functions can combine synergistically to mediate

470

physiological collapse. This effect is the most likely explanation for the enhanced efficacy of

471

stacking dsRNA constructs against multiple osmoregulation genes in the phloem-feeding aphid

472

Myzus persicae and psyllid Bacteriocera cockerelli (Tzin et al., 2015). In our hands, however,

473

the performance of the whitefly B. tabaci is not significantly reduced by RNAi against

474

osmoregulation genes, either individually or stacked (this study; unpub results); and RNAi

475

efficacy is dependent on co-targeting of the dsRNase genes (this study).

SC

476

RI PT

469

Also pertinent to our understanding of the efficacy of RNAi in whiteflies are the results of Raza et al. (2016), who transformed tobacco plants with dsRNA against the B. tabaci AQP1 and

478

SUC2 genes, obtaining 80->98% expression knockdown of the target genes and 60-80% insect

479

mortality over 6 days. Further research is required to identify the factors contributing to the

480

differences between the two studies, potentially including variation in the expression of gut

481

nuclease genes among different genotypes of B. tabaci, and differences in the identity of phloem-

482

mobile dsRNA molecules between tobacco used by Raza et al. (2016) and tomato used in this

483

study. In addition, Raza et al. (2016) used a dsRNA construct against SUC2, which lacks a

484

signal peptide sequence (Jing et al., 2016). The SUC2 protein likely has an intracellular

485

function, as distinct from SUC1, used in this study, which has a signal peptide permitting export

486

to the gut lumen and osmoregulatory function (Jing et al., 2016).

TE D

EP

AC C

487

M AN U

477

In conclusion, this study has provided proof-of-principle that the efficacy of RNAi against the

488

whitefly B. tabaci can be enhanced by the dual strategies of, first, stacking dsRNA against

489

multiple genes with related physiological roles but distinct molecular functions and, second,

490

“using RNAi to suppress suppressors of RNAi”. These approaches create wider opportunities

491

for RNAi as a control strategy against B. tabaci, which is a globally important crop pest, and

ACCEPTED MANUSCRIPT

492

offer a template for comparable strategies against other insect pests with poor or variable RNAi

493

efficacy.

494 Acknowledgements

496

We thank Patricia Keen and Michelle Tjahjadi (Boyce Thompson Institute) for technical support

497

in preparing the transgenic plants. This work was funded by the Academic Venture Fund of the

498

Atkinson Center for a Sustainable Future at Cornell University and the Natural Resources

499

Institute, University of Greenwich from a grant provided by the Bill & Melinda Gates foundation

500

(Grant Agreement OPP1058938).

M AN U

501

SC

RI PT

495

Contributions of authors

503

YL: conducted the experiments using ds-GFP and osmoregulation genes on artificial diets, and

504

phylogenetic analysis of dsRNase genes. QC: designed and generated the constructs for plant

505

transgenesis. JL and SHC: conducted qPCR experiments. JVE: conducted the plant

506

transformations. YL, JL, SHC, JVE, RT and AED: devised the experimental plan, interpreted the

507

results and wrote the manuscript. All coauthors have commented on drafts of the manuscript,

508

and have read and agree to the submitted version.

EP

AC C

509

TE D

502

510

References

511

Allen, M.L., and Walker, W.B. 3rd, 2012. Saliva of Lygus lineolaris digests double stranded

512

ribonucleic acids. J. Insect Physiol. 58, 391-396.

ACCEPTED MANUSCRIPT

513

Arimatsu, Y., Kotani, E., Sugimura, Y., Furusawa, T., 2007. Molecular characterization of a

514

cDNA encoding extracellular dsRNase and its expression in the silkworm, Bombyx mori.

515

Insect Biochem. Mol. Biol. 37, 176-183. Ayre, B.G., Keller, F., Turgeon, R., 2003. Symplastic continuity between companion cells and

RI PT

516 517

the translocation stream: long-distance transport is controlled by retention and retrieval

518

mechanisms in the phloem. Plant Physiol. 131, 1518-1528.

521

SC

520

Baulcombe, D.C., 2015. VIGS, HIGS and FIGS: small RNA silencing in the interactions of viruses or filamentous organisms with their plant hosts. Curr. Opin. Plant Biol. 26, 141-146. Baum, J.A., Bogaert, T., Clinton, W., Heck, G.R., Feldmann, P., Ilagan, O., Johnson, S.,

M AN U

519

522

Plaetinck, G., Munyikwa, T., Pleau, M., et al., 2007. Control of coleopteran insect pests

523

through RNA interference. Nat. Biotechnol. 25, 1322-1326.

525 526

Belles, X., 2010. Beyond Drosophila: RNAi in vivo and functional genomics in insects. Annu. Rev. Entomol. 55, 111-128.

TE D

524

Chen, W., Hasegawa, D.K., Kaur, N., Kliot, A., Pinheiro, P.V., Luan, J., Stensmyr, M.C., Zheng, Y., Liu, W., Sun, H., et al., 2016. The draft genome of whitefly Bemisia tabaci MEAM1, a

528

global crop pest, provides novel insights into virus transmission, host adaptation, and

529

insecticide resistance. BMC Biol. 14, 110. Christiaens, O., Swevers, L., Smagghe, G., 2014. DsRNA degradation in the pea aphid

AC C

530

EP

527

531

(Acyrthosiphon pisum) associated with lack of response in RNAi feeding and injection assay.

532

Peptides 53, 307-314.

533

Coleman, A.D., Wouters, R.H., Mugford, S.T., Hogenhout, S.A., 2015. Persistence and

534

transgenerational effect of plant-mediated RNAi in aphids. J. Exp. Bot. 66, 541-548.

ACCEPTED MANUSCRIPT

535 536 537

Douglas, A.E., Minto, L.B., Wilkinson, T.L., 2001. Quantifying nutrient production by the microbial symbionts in an aphid. J. Exp. Biol. 204, 349-358. Dunoyer, P., Himber, C., Ruiz-Ferrer, V., Alioua, A., Voinnet, O., 2007. Intra- and intercellular RNA interference in Arabidopsis thaliana requires components of the microRNA and

539

heterochromatic silencing pathways. Nat. Genet. 39, 848-856.

RI PT

538

Dunoyer, P., Schott, G., Himber, C., Meyer, D., Takeda, A., Carrington, J.C., Voinnet, O., 2010.

541

Small RNA duplexes function as mobile silencing signals between plant cells. Science 328,

542

912-916.

SC

540

Garbutt, J.S., Belles, X., Richards, E.H., Reynolds, S.E., 2013. Persistence of double-stranded

544

RNA in insect hemolymph as a potential determiner of RNA interference success: evidence

545

from Manduca sexta and Blattella germanica. J. Insect Physiol. 59, 171-178.

546

M AN U

543

Garbutt, J.S., Reynolds, S.E., 2012. Induction of RNA interference genes by double-stranded RNA; implications for susceptibility to RNA interference. Insect Biochem. Mol. Biol. 42, 621-

548

628.

549

TE D

547

Ghanim, M., Kontsedalov, S., Czosnek, H., 2007. Tissue-specific gene silencing by RNA interference in the whitefly Bemisia tabaci (Gennadius). Insect Biochem. Mol. Biol. 37, 732-

551

738.

Haritatos, E., Ayre, B.G., Turgeon, R., 2000. Identification of phloem involved in assimilate

AC C

552

EP

550

553

loading in leaves by the activity of the galactinol synthase promoter. Plant Physiol. 123, 929-

554

937.

555 556

Huvenne, H., Smagghe, G., 2010. Mechanisms of dsRNA uptake in insects and potential of RNAi for pest control: a review. J. Insect Physiol. 56, 227-235.

ACCEPTED MANUSCRIPT

557

Jing, X., White, T.A., Luan, J., Jiao, C., Fei, Z., Douglas, A.E., 2016. Evolutionary conservation

558

of candidate osmoregulation genes in plant phloem sap-feeding insects. Insect Mol. Biol. 25,

559

251-258. Kempers, R., Prior, D.A.M., van Bel, A.J.E., Oparka, K.J., 1993. Plasmodesmata between sieve

RI PT

560 561

element and companion cell of extrafascicular stem phloem of Cucurbita maxima permit

562

passage of 3 kDa fluorescent probes. Plant J. 4, 567-575.

565

SC

564

Liu, D.D., Chao, W.M., Turgeon, R. 2012a. Transport of sucrose, not hexose, in the phloem. J. Exp. Bot. 63, 4315-4320.

Liu, J., Swevers, L., Iatrou, K., Huvenne, H., Smagghe, G., 2012b. Bombyx mori DNA/RNA

M AN U

563

566

non-specific nuclease: expression of isoforms in insect culture cells, subcellular localization

567

and functional assays. J. Insect Physiol. 58, 1166-1176.

569 570

Livak, K.J., Schmittgen, T.D., 2001. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods 25, 402-408.

TE D

568

Luan, J.B., Ghanim, M., Liu, S.S., Czosnek, H., 2013. Silencing the ecdysone synthesis and signaling pathway genes disrupts nymphal development in the whitefly. Insect Biochem. Mol.

572

Biol. 43, 740-746.

EP

571

Luo, Y., Wang, X., Wang, X., Yu, D., Chen, B., Kang, L., 2013. Differential responses of

574

migratory locusts to systemic RNA interference via double-stranded RNA injection and

575

feeding. Insect Mol. Biol. 22, 574-583.

576

AC C

573

Mao, Y.B., Cai, W.J., Wang, J.W., Hong, G.J., Tao, X.Y., Wang, L.J., Huang, Y.P., Chen, X.Y.,

577

2007. Silencing a cotton bollworm P450 monooxygenase gene by plant-mediated RNAi

578

impairs larval tolerance of gossypol. Nat. Biotechnol. 25, 1307-1313.

ACCEPTED MANUSCRIPT

579

Mathew, L.G., Campbell, E.M., Yool, A.J., Fabrick, J.A., 2011. Identification and

580

characterization of functional aquaporin water channel protein from alimentary tract of

581

whitefly, Bemisia tabaci. Insect Biochem. Mol. Biol. 41, 178-190. Melnyk, C.W., Molnar, A., Bassett, A., Baulcombe, D.C., 2011a. Mobile 24 nt small RNAs

RI PT

582 583

direct transcriptional gene silencing in the root meristems of Arabidopsis thaliana. Curr. Biol.

584

21, 1678-1683.

589 590 591 592 593 594 595 596 597

SC

M AN U

588

Mutti, N.S., Park, Y., Reese, J.C., Reeck, G.R., 2006. RNAi knockdown of a salivary transcript leading to lethality in the pea aphid, Acyrthosiphon pisum. J. Insect Sci. 6, 1-7. Navas-Castillo, J., Fiallo-Olive, E., Sanchez-Campos, S., 2011. Emerging virus diseases transmitted by whiteflies. Annu. Rev. Phytopathol. 49, 219-248.

Pall, G.S., Hamilton, A.J., 2008. Improved northern blot method for enhanced detection of small

TE D

587

RNA silencing signals. EMBO J. 30, 3553-3563.

RNA. Nat. Protoc. 3, 1077-1084.

Petersen, T.N., Brunak, S., von Heijne, G., and Nielsen, H., 2011. SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods 8, 785-786.

EP

586

Melnyk, C.W., Molnar, A., Baulcombe, D.C., 2011b. Intercellular and systemic movement of

Pitino, M., Coleman, A.D., Maffei, M.E., Ridout, C.J., Hogenhout, S.A., 2011. Silencing of aphid genes by dsRNA feeding from plants. PLoS One 6, e25709.

AC C

585

Raza, A., Malik, H.J., Shafiq, M., Amin, I., Scheffler, J.A., Scheffler, B.E., Mansoor, S., 2016.

598

RNA interference based approach to down regulate osmoregulators of whitefly (Bemisia

599

tabaci): potential technology for the control of whitefly. PLoS One 11, e0153883.

ACCEPTED MANUSCRIPT

600

Saleh, M.C., van Rij, R.P., Hekele, A., Gillis, A., Foley, E., O'Farrell, P.H., Andino, R., 2006.

601

The endocytic pathway mediates cell entry of dsRNA to induce RNAi silencing. Nat. Cell

602

Biol. 8, 793-802. Scott, J.G., Michel, K., Bartholomay, L.C., Siegfried, B.D., Hunter, W.B., Smagghe, G., Zhu,

RI PT

603 604

K.Y., Douglas, A.E., 2013. Towards the elements of successful insect RNAi. J. Insect Physiol.

605

59, 1212-1221.

Spit, J., Philips, A., Wynant, N., Santos, D., Plaetinck, G., Vanden Broeck, J., 2017. Knockdown

SC

606

of nuclease activity in the gut enhances RNAi efficiency in the Colorado potato beetle,

608

Leptinotarsa decemlineata, but not in the desert locust, Schistocerca gregaria. Insect Biochem.

609

Mol. Biol. 81, 103-116.

611 612

Su, Y.L., Li, J.M., Li, M., Luan, J.B., Ye, X.D., Wang, X.W., Liu, S.S., 2012. Transcriptomic analysis of the salivary glands of an invasive whitefly. PLoS One 7, e39303. Terenius, O., Papanicolaou, A., Garbutt, J.S., Eleftherianos, I., Huvenne, H., Kanginakudru, S.,

TE D

610

M AN U

607

Albrechtsen, M., An, C., Aymeric, J.L., Barthel, A., et al., 2011. RNA interference in

614

Lepidoptera: an overview of successful and unsuccessful studies and implications for

615

experimental design. J. Insect Physiol. 57, 231-245.

616

EP

613

Thakur, N., Upadhyay, S.K., Verma, P.C., Chandrashekar, K., Tuli, R., Singh, P.K., 2014. Enhanced whitefly resistance in transgenic tobacco plants expressing double stranded RNA of

618

v-ATPase A gene. PLoS One 9, e87235.

619

AC C

617

Truernit, E., Sauer, N., 1995. The promoter of the Arabidopsis thaliana SUC2 sucrose-H+

620

symporter gene directs expression of beta-glucuronidase to the phloem: evidence for phloem

621

loading and unloading by SUC2. Planta 196, 564-570.

ACCEPTED MANUSCRIPT

622

Tzin, V., Yang, X., Jing, X., Zhang, K., Jander, G., Douglas, A.E., 2015. RNA interference

623

against gut osmoregulatory genes in phloem-feeding insects. J. Insect Physiol. 79, 105-112.

624

Ulvila, J., Parikka, M., Kleino, A., Sormunen, R., Ezekowitz, R.A., Kocks, C., Ramet, M,. 2006. Double-stranded RNA is internalized by scavenger receptor-mediated endocytosis in

626

Drosophila S2 cells. J. Biol. Chem. 281, 14370-14375.

RI PT

625

Upadhyay, S.K., Chandrashekar, K., Thakur, N., Verma, P.C., Borgio, J.F., Singh, P.K., Tuli, R.,

628

2011. RNA interference for the control of whiteflies (Bemisia tabaci) by oral route. J. Biosci.

629

36, 153-161.

SC

627

Van Eck, J., Kirk, D.D., Walmsley, A.M., 2006. Tomato (Lycopersicon esculentum). In Methods

631

in Molecular Biology: Agrobacterium Protocols, K. Wang, ed. (Totowa, NJ: Humana Press

632

Inc.), pp. 459-473.

M AN U

630

Wang, F., Johnson, N.R., Coruh, C., Axtell, M.J., 2016a. Genome-wide analysis of single non-

634

templated nucleotides in plant endogenous siRNAs and miRNAs. Nucleic Acids Res. 44,

635

7395-7405.

TE D

633

Wang, K., Peng, Y., Pu, J., Fu, W., Wang, J., Han, Z. 2016b. Variation in RNAi efficacy among

637

insect species is attributable to dsRNA degradation in vivo. Insect Biochem. Mol. Biol. 77, 1-

638

9.

Wang, X.W., Luan, J.B., Li, J.M., Bao, Y.Y., Zhang, C.X., Liu, S.S., 2010. De novo

AC C

639

EP

636

640

characterization of a whitefly transcriptome and analysis of its gene expression during

641

development. BMC Genomics 11, 400.

642

Wang, X.W., Zhao, Q.Y., Luan, J.B., Wang, Y.J., Yan, G.H., Liu, S.S., 2012. Analysis of a

643

native whitefly transcriptome and its sequence divergence with two invasive whitefly species.

644

BMC Genomics 13, 529.

ACCEPTED MANUSCRIPT

646 647 648 649

Whyard, S., Singh, A.D., Wong, S., 2009. Ingested double-stranded RNAs can act as speciesspecific insecticides. Insect Biochem. Mol. Biol. 39, 824-832. Wilson, R.C., Doudna, J.A., 2013. Molecular mechanisms of RNA interference. Annu. Rev. Biophys. 42, 217-239.

RI PT

645

Wu, P., Li, M.-W., Wang, X., Zhao, P., Wang, X.-Y., Liu, T., Qin, G.-X., Guo, X.-J., 2009. Differentially expressed genes in the midgut of silkworm infected with cytoplasmic

651

polyhedrosis virus. Afr. J. Biotechnol. 8, 3711-3720.

652

SC

650

Wynant, N., Santos, D., Verdonck, R., Spit, J., Van Wielendaele, P., Vanden Broeck, J., 2014. Identification, functional characterization and phylogenetic analysis of double stranded RNA

654

degrading enzymes present in the gut of the desert locust, Schistocerca gregaria. Insect

655

Biochem. Mol. Biol. 46, 1-8.

656

M AN U

653

Xue, X.Y., Mao, Y.B., Tao, X.Y., Huang, Y.P., Chen, X.Y., 2012. New approaches to agricultural insect pest control based on RNA interference. Advances Insect Physiol. 42, 73-

658

113.

659

TE D

657

Yoo, B.C., Kragler, F., Varkonyi-Gasic, E., Haywood, V., Archer-Evans, S., Lee, Y.M., Lough, T.J., Lucas, W.J., 2004. A systemic small RNA signaling system in plants. Plant Cell 16,

661

1979-2000.

Zamore, P.D., Tuschl, T., Sharp, P.A., Bartel, D.P., 2000. RNAi: double-stranded RNA directs

AC C

662

EP

660

663

the ATP-dependent cleavage of mRNA at 21 to 23 nucleotide intervals. Cell 101, 25-33.

664

Zha, W., Peng, X., Chen, R., Du, B., Zhu, L., He, G., 2011. Knockdown of midgut genes by

665

dsRNA-transgenic plant-mediated RNA interference in the hemipteran insect Nilaparvata

666

lugens. PLoS One 6, e20504.

ACCEPTED MANUSCRIPT

667

Zhai, J., Zhao, Y., Simon, S.A., Huang, S., Petsch, K., Arikit, S., Pillay, M., Ji, L., Xie, M., Cao,

668

X., et al., 2013. Plant microRNAs display differential 3' truncation and tailing modifications

669

that are ARGONAUTE1 dependent and conserved across species. Plant Cell 25, 2417-2428.

671 672

Zhang, C., Yu, X., Ayre, B.G.,Turgeon, R., 2012. The origin and composition of cucurbit "phloem" exudate. Plant Physiol. 158, 1873-1882.

RI PT

670

Zhang, J., Khan, S.A., Hasse, C., Ruf, S., Heckel, D.G., Bock, R., 2015. Pest control. Full crop protection from an insect pest by expression of long double-stranded RNAs in plastids.

674

Science 347, 991-994.

AC C

EP

TE D

M AN U

SC

673

ACCEPTED MANUSCRIPT

Figure legends

676

Fig 1. Northern blot analyses of ds-GFP in transgenic tomato plants and whiteflies feeding on

677

the plants. (A) Map of transformation vectors for ds-GFP expression, designed to produce 370

678

bp hairpin-ds-GFP (hp-ds-GFP) under companion cell-specific promoter CmGAS or AtSUC2;

679

(B) 370 nt ds-GFP (top) and 21-25 nt small RNA (bottom) in the plants. (C) 370 nt ds-GFP

680

(horizontal arrowhead) in cohorts of ca. 250 whiteflies that had fed on transgenic plants for two

681

weeks. (D) RNA gel blot analysis of RNA extracted from the scion apex heterografted (GAS:4

682

and SUC2:11) plants. Numbers indicate transgenic tomato lines, with Ev comprising plants

683

transformed with the empty vector, and WT comprising non-transformed tomato plants. Twenty

684

µg total leaf or whitefly low molecular weight RNA (LMW) and high molecular weight RNA

685

(HMW) were separated on 8% denatured urea-acrylamide gel, with the HMW RNA

686

hybridization protocol for full-length ds-GFP detection and LMW RNA hybridization protocol

687

for small RNA detection. The rRNA loading control is SybrGold stained 5S rRNA.

SC

M AN U

TE D

688

RI PT

675

Fig 2. Northern blot analyses of ds-GFP administered to whiteflies via artificial diet. (A)

690

Northern blot of diet supplemented with 1 µg ds-GFP µl-1 and harvested immediately (0 h) and

691

after incubation for 72 h without whiteflies (-) or with whiteflies (+), with 100 ng RNA loaded

692

per well. (B) Northern blot of whiteflies that had fed on two diets containing ds-GFP (+) or ds-

693

GFP-free diets (-), with 20 µg whitefly total RNA loaded in each lane and Sybr-Gold stained

694

rRNA as loading control. The LMW RNA hybridization protocol was used for both blots.

AC C

695

EP

689

696

Fig 3. The dsRNase genes of Bemisia tabaci. (A) Neighbor-joining phylogenetic tree constructed

697

using the protein sequence of the conserved DNA/RNA non-specific nuclease domain of insect

ACCEPTED MANUSCRIPT

dsRNase genes. The numbers at the branches indicate the %bootstrap support, based on the

699

frequency of the clusters for 1000 bootstraps. The dsRNase sequences were: Bombyx mori 2

700

(NP_001091744.1), Papilio machaon (XP_014355571.1), Spodoptera littorallis (CAR92522.1),

701

Spodoptera frugiperda (CAR92521.1), Aedes aegypti (XP_001648469.1), Drosophila

702

melanogaster 1 (NM_140821.4), D. melanogaster 2 (NP_649078.1, CG3819), Tribolium

703

castaneum 1 (XP_973587.1), T. castaneum 2 (XP_970494.1), Schistocerca gregaria 1

704

(KJ135008), S.gregaria 2 (KJ135009), S.gregaria 3 (KJ135010), S.gregaria 4 (KJ135011),

705

Acyrthosiphon pisum (ACYPI008471), Myzus persicae

706

(MYZPE13164_0_v1.0_000125730.4_pep) and Bemisia tabaci 1(KX390872), B. tabaci 2

707

(KX390873) and B. tabaci 3 (Unigene11878_BT_Q_SG_ZJU). (B) qRT-PCR analysis of the

708

expression of BtdsRNase-1 (top) and BtdsRNase-2 (bottom) in dissected guts of B. tabaci,

709

relative to the whole body (Wb). Mean ± s.e. from 3 replicates are shown.

M AN U

SC

RI PT

698

TE D

710

Fig 4. Effect of RNAi against dsRNase on the abundance of ds-GFP administered orally to

712

whiteflies. (A) Experimental design. 200 ng µl-1 of each ds-dsRNase and 500 ng µl-1 of ds-GFP

713

were fed to the whiteflies. (B) Relative density of the 370 nt ds-GFP bands on northern blots

714

(mean + s.e, 3 replicate experiments), with ANOVA results and significantly different treatments

715

by post hoc test indicated by different letters. The northern blots are displayed in Fig. S3.

AC C

716

EP

711

717

Fig 5. RNAi against osmoregulation genes in adult B. tabaci on artificial diet. dsRNA (100 ng

718

µl-1 diet) were ds-GFP, ds-dsRNase1&2, and dsRNA against osmoregulation genes AQP1 and

719

SUC1 with (+) or without (-) ds-dsRNase1&2 administered over 6 days. A&B. qRT-PCR

720

analysis of AQP1 and SUC1 expression, respectively, relative to dsRNA-free diet, normalized to

ACCEPTED MANUSCRIPT

RPL13. Mean + s.e, 3 reps. The expression of SUC1 in insects fed on ds-AQP1, and of AQP1 in

722

insects fed on ds-SUC1, was quantified only in the presence of ds-dsRNase1&2. C. Number of

723

dead insects (mean + s.e. for 10 replicates of 40 insects) at day-6. Different letters indicate

724

significantly different treatments by Fisher’s LSD test.

725 726 Supplementary Material

SC

727 728

Table S1. Primers

M AN U

729

RI PT

721

730

Fig. S1. Northern blot analyses of ds-GFP in transgenic tomato plants and whiteflies feeding on

732

the plants. (A) Plant RNA with HMW RNA extraction and HMW RNA hybridization protocols.

733

(B) Plant RNA with LMW RNA extraction and LMW RNA hybridization protocols. (C)

734

Whitefly with HMW RNA extraction and HMW RNA hybridization protocols. (D) Whitefly

735

with LMW RNA extraction and LMW RNA hybridization protocols. (E) RNA extracted from

736

the scion apex heterografted (GAS:4 and SUC2:11) plants. RNA was extracted from tomato

737

leaves (A,B,E) and from ca. 250 adult whiteflies that had fed on transgenic plants for two weeks

738

(C,D), and the samples (20 µg RNA per lane) were separated on 8% denatured urea-acrylamide

739

gel. Numbers indicate transgenic tomato lines, with Ev, the negative control, comprising tomato

740

plants transformed with the empty vector. WT, wild-type tomato, MW, molecular weight. The

741

rRNA loading control is SybrGold stained 5S rRNA. The bands of molecular weight

742

intermediate between 370 nt GFP and candidate sRNAs (21-25 nt) were interpreted as non-

743

specific binding of plant and whitefly RNA to the GFP probe because they were detected in both

AC C

EP

TE D

731

ACCEPTED MANUSCRIPT

the experimental samples and Ev controls. The LMW RNA extraction protocol (B, D) reduced,

745

but did not exclude HMW RNA, particularly yielding multiple bands >> 21-25 bp in the plant

746

samples (B) that were detectable with the low stringency of the LMW hybridization protocol

747

required to visualize 21-25 bp bands.

RI PT

744

748

Fig. S2. Expression of dsRNase genes in whiteflies administered ds-dsRNase1, ds-dsRNase2 or

750

ds-dsRNase1&2. The experimental design is shown in Fig. 4A. (A) Relative expression level

751

of dsRNase1. (B) Relative expression level of dsRNase2. Gene expression was normalized to

752

RPL13, and mean ± s.e. from 3 replicates are shown. Different letters indicate significantly

753

different treatments by post hoc test following one-way ANOVA.

M AN U

SC

749

754

Fig. S3. Northern blots of ds-GFP fragment following treatment with ds-dsRNase1, ds-dsRNase2

756

and ds-dsRNase1&2, with SybrGold–stained 5S rRNA as loading control. The LMW RNA

757

hybridization protocol was used for all blots. The exposure time for each blot was optimized

758

independently, to ensure that none of the bands on each blot was saturated and so ensure the

759

reliability of the data in Fig. 4B. Biological variability and the criteria for optimizing exposure

760

time (above) contribute to the among-experiment variation. No: no treatment; Syn: synchronous

761

treatment with ds-GFP and ds-RNase; Pre: pretreatment with ds-dsRNase (see Fig. 4A for

762

details).

EP

AC C

763

TE D

755

764

Fig. S4. Transcription profiles of SUC1 (A) and AQP1 (B) in the whole body and dissected gut

765

of adult whiteflies. Gene expression was normalized to RPL13. Mean ± s.e. from 3 replicates are

766

shown. * indicate significantly different by Student’s t-test.

ACCEPTED MANUSCRIPT

767 Fig. S5. Expression of dsRNase genes in adult whiteflies reared for 6 days on artificial diets

769

supplemented with ds-dsRNase1&2 (+), with ds-dsRNase1&2-free diets (-) as control. (A)

770

Relative expression level of dsRNase1 (B) Relative expression level of dsRNase2. Gene

771

expression was normalized to RPL13, and mean ± s.e. from 3 replicates are shown. Different

772

letters indicate significantly different treatments by post hoc test following one-way ANOVA.

AC C

EP

TE D

M AN U

SC

RI PT

768

A

RI PT

ACCEPTED MANUSCRIPT

GFP

intron

GFP

terminator

dsGFP 0.37 kb

CmGAS1/AtSUC2

Ev

CmGAS1

3

4

5

AtSUC2 13

7

11 12

370 nt

M AN U

promoter

SC

B

rRNA

21-25 nt

rRNA 21-25 nt rRNA

Ev

AtSUC2 11

12

EP

370 nt

CmGAS1 3 4 5

AC C

C

TE D

rRNA

D 370 nt rRNA 21-25 nt rRNA

WT

CmGAS1

AtSUC2

4

11

4

11

ACCEPTED MANUSCRIPT

RI PT

M AN U

SC

370 nt

TE D

---

Whitefly Whitefly +/-+/dsGFP dsGFP + + ++ + + + 370 nt 370 nt

-- -

EP

370 nt

B

dsGFP +/- whitefly 0h 72h + + +

AC C

A A

21-25 nt 21-25 nt rRNA rRNA

AC C

EP

TE D

M AN U

SC

RI PT

ACCEPTED MANUSCRIPT

AC C

EP

TE D

M AN U

SC

RI PT

ACCEPTED MANUSCRIPT

AC C

EP

TE D

M AN U

SC

RI PT

ACCEPTED MANUSCRIPT

ACCEPTED MANUSCRIPT

Highlights

RI PT

SC M AN U

•

TE D

•

EP

• •

Tomato plants transformed with dsRNA deliver full-length dsRNA to the phloem-feeding whitefly Bemisia tabaci Orally-acquired dsRNA is degraded in the body of the whitefly B. tabaci The B. tabaci genome encodes three non-specific nucleases, two of which are expressed in the gut RNAi-mediated suppression of nuclease expression reduces degradation of other dsRNA molecules ingested by whiteflies Whitefly mortality on diets is enhanced by stacking oral RNAi against insect osmoregulation genes and nucleases

AC C

•