# RNAi-based reverse genetics in the chelicerate model Tetranychus urticae: A comparative analysis of five methods for gene silencing

Takeshi Suzuki, Maria Andreia Nunes, María Urizarna España, Hooman Hosseinzadeh Namin, Pengyu Jin, Nicolas Bensoussan, Vladimir Zhurov, Tawhid Rahman, Rebecca de Clercq, Pierre Hilson, et al.

<figure><img src="image_1.png" /></figure>

► To cite this version:

Takeshi Suzuki, Maria Andreia Nunes, María Urizarna España, Hooman Hosseinzadeh Namin, Pengyu Jin, et al. RNAi-based reverse genetics in the chelicerate model Tetranychus urticae: A comparative analysis of five methods for gene silencing. PLoS ONE, 2017, 12 (7), e02627337.

HAL Id: hal-02627337

https://hal.inrae.fr/hal-02627337

Submitted on 26 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# RNAi-based reverse genetics in the chelicerate model *Tetranychus urticae*: A comparative analysis of five methods for gene silencing

Takeshi Suzuki¹,*, Maria Andreia Nunes¹,*, María Urizarna España¹, Hooman Hosseinzadeh Namin¹, Pengyu Jin¹, Nicolas Bensoussan¹, Vladimir Zhurov¹, Tawhid Rahman¹,*, Rebecca De Clercq², Pierre Hilson²,³,⁴, Vojislava Grbic¹*, Miodrag Grbic¹,⁵*

1 Department of Biology, The University of Western Ontario, London, Ontario, Canada, 2 Department of Plant Systems Biology, VIB, Ghent, Belgium, 3 Department of Plant Biotechnology and Bioinformatics, Ghent University, Ghent, Belgium, 4 Institut Jean-Pierre Bourgin, UMR1318 INRA-AgroParisTech, Saclay Plant Science, Versailles, France, 5 University of La Rioja, Logrono, Spain

* These authors contributed equally to this work.

a Current address: Graduate School of Bio-Applications and Systems Engineering, Tokyo University of Agriculture and Technology, Koganei, Tokyo, Japan

b Current address: Centro de Citricultura Sylvio Moreira, Instituto Agronômico de Campinas, Cordeirópolis, São Paulo, Brazil

c Current address: Saskatoon Research Centre, Agriculture and Agri-Food Canada, Saskatoon, SK, Canada

* mgrbic@uwo.ca (MG); vgrbic@uwo.ca (VG)

<figure><img src="image_3.png" /></figure>

## Abstract

RNA interference (RNAi) can be used for the protection against agricultural pests through the silencing of genes required for pest fitness. To assess the potential of RNAi approaches in the two-spotted spider mite, *Tetranychus urticae*, we compared 5 methods for the delivery of double-stranded RNA (dsRNA). These methods include mite feeding on either (i) leaves floating on a dsRNA solution, (ii) dsRNA-expressing plants, (iii) artificial diet supplemented with dsRNA, or (iv) dsRNA-coated leaves, and (v) mite soaking in a dsRNA solution. In all cases, the gene targeted for method validation was the Vacuolar-type H⁺-ATPase (*TuVAT-Pase*), encoding a constitutively expressed ATP-driven proton pump located in the membrane. Down-regulation of *TuVATPase* increased mortality and/or reduced fecundity in all methods, but with variable efficiency. The most efficient methods for dsRNA delivery were direct soaking of mites in the dsRNA solution and mite feeding on dsRNA-coated leaves that mimics dsRNA application as a sprayable pesticide. Both resulted in a dark-body phenotype not observed in mites treated with a control dsRNA. Although with lower efficiency, dsRNA designed for *TuVATPase* silencing and expressed in transgenic *Arabidopsis* plants impacted the fitness of mites feeding on these plants. RNAi may thus be a valuable strategy to control spider mite populations, either as a sprayable pesticide or through transgenic crops. This comparative methodological study focusing on the induction of RNAi-based gene silencing in *T. urticae* paves the way for reverse genetics approaches in this model chelicerate system and prepares large-scale systematic RNAi screens as a first step towards the development of specific RNA-based pesticides. Such alternative molecules

<figure><img src="image_4.png" /></figure>

### OPEN ACCESS

**Citation:** Suzuki T, Nunes MA, España MU, Namin HH, Jin P, Bensoussan N, et al. (2017) RNAi-based reverse genetics in the chelicerate model *Tetranychus urticae*: A comparative analysis of five methods for gene silencing. PLoS ONE 12(7): e0180654. https://doi.org/10.1371/journal.pone.0180654

**Editor:** Mohammed Bendahmane, Ecole Normale Supérieure, FRANCE

**Received:** February 6, 2017

**Accepted:** June 19, 2017

**Published:** July 12, 2017

**Copyright:** © 2017 Suzuki et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper.

**Funding:** This work was supported by the Government of Canada through the Ontario Research Fund—Research Excellence Round 8 (RE08-067) (to MG and VG) and the University of Western Ontario through the Western Strategic Support program (to VG). TS has been supported by the Postdoctoral Fellowship for Research Abroad and KAKENHI (16K18661) funded by the
Japan Society for the Promotion of Science (JSPS) and MAN by the National Counsel of Technological and Scientific Development (CNPq/Brazil).

**Competing interests:** The authors declare that no competing interests exist.

may help control spider mites that cause significant damages to crops and ornamental plant species, as well as other chelicerates detrimental to agriculture and health.

## Introduction

The chelicerates represent the second largest group of terrestrial animals after insects [1]. This basal lineage of the phylum Arthropoda includes horseshoe crabs, scorpions, spiders, mites and ticks that display a plethora of different lifestyles and include economically important species for human health and agriculture. Mite species exhibit a large range of adaptations including herbivory, predation, parasitism, detritivory and symbiosis [2,3]. With almost 50,000 species described by the end of 20th century, and 0.5 to 1 million species estimated to exist, mites are one of the most diverse groups in the animal kingdom [4]. The two spotted spider mite (TSSM), *Tetranychus urticae* (Koch), is the first chelicerate whose complete genome was sequenced and annotated [5]. It is a compact genome of 90 Mbp (54% of which is protein coding sequence), with simple gene structure, complemented with several transcriptome databases. TSSM develops rapidly. It is easy to maintain in the laboratory and can be enriched in specific developmental stages. Because it is a major agricultural herbivorous pest, a strong research community supports the emergence of *T. urticae* as a versatile chelicerate model organism [6].

Most organisms whose genome has been sequenced in the last few years lack genetic tools available in established model systems (e.g. *Caenorhabditis elegans*, *Drosophila*, zebrafish and mouse). Thus, it is critical to develop high throughput platforms for reverse genetics, aimed at the functional characterization of unknown genes. Since its discovery, RNAi has been a useful tool to dissect gene function [7–9]. However, RNAi delivery is often a cumbersome process requiring the injection of dsRNA into the body of small animals to initiate gene silencing. In such cases, RNAi experiments remain limited to the silencing of individual genes, with the notable exception of a high throughput screen of over 5,000 genes in *Tribolium* [10]. Alternatively, streamlined dsRNA delivery methods resulted in relatively large assays. For example, a medium size screen (hundreds of genes) was performed in the coleopteran western corn rootworm *Diabrotica virgifera virgifera* based on an artificial diet laced with dsRNA [11,12]. Also, dsRNA delivery by soaking was implemented in a high throughput reverse genetic platform for *C. elegans* [13,14].

The genome of TSSM contains genes encoding for a complete RNAi processing machinery [5], including 2 Dicer homologs, pasha, Drosha and components of the RISC complex (including 7 orthologs of both *Argonaute* and *Piwi* genes). In contrast to other insects and similar to *C. elegans*, *T. urticae* possesses 5 copies of RNA-dependent RNA polymerase (RdRp) [5] required for the systemic spread of RNAi responses. It has already been shown that the material injection of either dsRNA or siRNAs in *T. urticae* induces RNAi in embryos and causes developmental aberrations consistent with the loss-of-function phenotype of the target gene [15]. RNAi is thus a valid reverse genetic approach in *T. urticae*. However, compared to other arthropods, for example *Tribolium* and *Oncopeltus* [16,17], it is difficult to inject dsRNA in TSSM females that are less than 0.5 mm in length. More recently, Kwon *et al.* [18] reported an alternative method for the oral administration of dsRNA via bean leaf discs floating on a dsRNA solution. Unfortunately, this protocol requires large quantities of dissolved dsRNA that must be replaced daily, making this method impractical for high throughput screens.

In the present study, we used the TSSM target gene *Vacuolar-type H$^+$-ATPase (VATPase)*, which encodes a protein that functions as a pump transferring protons across cellular
membranes using the energy released by the ATP hydrolysis [19]. *TuVATPase* is a good target for methodological comparative analysis for several reasons. It has been previously used as a target for RNAi silencing in multiple arthropod systems, e.g. the Western corn rootworm (*Diarrotica virgifera* *virgifera*), the pea aphid (*Acyrthosiphon pisum*), the red flour beetle (*Tribolium castaneum*), the tobacco hornworm (*Manduca sexta*), the whitefly (*Bemisia tabaci*) and the Colorado potato beetle (*Leptinotarsa decemlineata*) [12,18,20–22]. These studies established the essential role of *VATPase*, demonstrating that effective silencing of *VATPase* expression leads to significant and measurable reduction of fitness across arthropod species. In addition, these studies showed that orally-delivered dsRNAs against *VATPase* induce RNAi. Finally, *TuVATPase* has already been shown to yield modest but significant mite mortality when silenced in *T. urticae* through the leaf floating method [18].

We describe several methods for dsRNA delivery into spider mites towards the development of reverse genetics platforms: soaking, which may be applicable to high-throughput screens; leaf coating, that mimics the application of dsRNA as a sprayable pesticide; and artificial diet for the administration of quantifiable amounts of dsRNA to TSSM. In addition, we tested whether transgenic plants expressing dsRNA may relay silencing into the mites.

We show that all 4 tested dsRNA delivery methods had greater efficiency in triggering RNAi in comparison to the leaf floating benchmark. The leaf coating and soaking methods, compatible with the high throughput experimental setups, had the greatest efficiency in silencing *TuVATPase*. Furthermore, because dsRNA expressed in transgenic plants or applied on the leaf surface induced RNAi affecting mite fitness, they may become valuable tools for crop protection against spider mites.

## Materials and methods

### *T. urticae* strains and rearing

Stock populations of *T. urticae* London reference strain and inbred lines were maintained on California red kidney beans (*Phaseolus vulgaris* L.) (Stokes, Thorold, ON) in a climate controlled chamber at L:D 16:8, 26°C, and 50% relative humidity (RH). Inbred lines were derived from isofemale lines that underwent 6 consecutive generations of mother-son mating.

### Preparation of developmentally synchronized mites

A detailed protocol for the synchronization of mite development (larval and adult emergence) has been described in an accompanied article [23]. Briefly, adult female mites were allowed to lay eggs for 24 h in a climate controlled chamber at L:D 16:8, 26°C and 50% RH. After 24 h, the adult females were removed and 1-day-old eggs were submerged in water. Three days later, water was removed and the 4-day-old eggs were allowed to dry. Newly-hatched larvae were collected after 3 h, yielding a tightly synchronized cohort. For the preparation of adults, adult female mites were placed onto a fresh detached bean leaf and were allowed to lay eggs at L:D 16:8, 26°C and 50% RH, for 24 h. Adult females were then removed from the leaf and the 1-day-old eggs were allowed to develop until female mites became teleiochrysalis (7–8 days). Teleiochrysalis were collected and incubated at 18°C and 100% RH for 24 h, after which they were returned to 26°C and 50% RH environment. Adult female mites, emerged within the 3 h upon transfer to 50% RH environment, were collected for experimentation.

### dsRNA fragments

We used 2 dsRNAs targeting the *TuVATPase* transcript to ensure the reproducibility and the specificity of the RNAi phenotype. Fragments corresponding to the coding region of
<figure><img src="image_5.png" /><figcaption>**Fig 1.** Fragments used for the synthesis of dsRNAs. (A) Schematic of the *TuVATPase* locus. DNA sequences used for the generation of dsRNA-*TuVATPase* are located in the 3rd exon, fragment A (214bp), and in the 4th exon, fragment B (416bp). UTR and coding sequences are shown in yellow and blue, respectively. (B) Schematic of the part of the scaffold 12 of *T. urticae* genome depicting the location of the 382 bp fragment that was used to synthesize dsRNA-NC.</figcaption></figure>

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

*TuVATPase*, referred to as fragment A (an upstream fragment of 214 bp within the 3rd exon of the *TuVATPase* locus) and fragment B (a 416 bp fragment within the 4th exon) were chosen as templates for the preparation of dsRNAs (Fig 1). Furthermore, the fragment B corresponds to the dsRNA sequence previously used by Kwon *et al*. [18] enabling direct comparison of method efficiency between the previous report and this study. In addition, a 382 bp intergenic fragment spanning the region 1690614–1690995 of the genomic scaffold 12 [5] (Fig 1), was chosen as the template for the preparation of a negative control dsRNA, referred to as NC. A BLAST search against the *T. urticae* genome and transcriptome databases confirmed that the 382 bp intergenic sequence is unique and not transcribed [24,25].

dsRNA preparation

The nucleotide sequences of *TuVATPase* (tetur09g04140), and intragenic region (negative control (NC), genomic coordinates: scaffold 12, position 1690614–1690995) were obtained from the ORCAE database [24]. Total RNA was extracted from the frozen mite females with the RNeasy Mini Kit (Qiagen, Valencia, CA) and cDNA was synthesized from 3 µg of the extracted total RNA with the SuperScript II cDNA Synthesis Kit (Thermo Fisher Scientific, Waltham, MA) and stored at –20°C. Fragments (A and B) of *TuVATPase* and a control sequence (NC) targeting an intergenic region (Fig 1) were PCR amplified with primers shown in Table 1, using cDNA and genomic DNA as templates, respectively, with Phusion DNA polymerase (NEB, New England Biolabs, UK). Amplified DNA fragments were purified with the Gel/PCR DNA Fragments Extraction Kit (Geneaid Biotech, New Taipei, Taiwan). Purified fragments were cloned in pLITMUS 38i vector (NEB, New England Biolabs, UK) containing T7 polymerase promotor sequences flanking multiple cloning sites to streamline dsRNA production. Inserts selected for the synthesis of dsRNA were re-sequenced to confirm their identity. RNA fragments were synthesized using 1 µg of DNA template with the TranscriptAid T7 High Yield Transcription Kit (Thermo Fisher Scientific) in 1.5 mL centrifuge tubes, denaturated at 95°C for 5 min, followed by slow cool-down to room temperature to facilitate
<table><caption>**Table 1. Primers used in this study.**</caption><thead><tr><th colspan="3">Cloning primers</th></tr><tr><th>Primer names</th><th>Oligonucleotide sequence (5' to 3')</th><th>Fragment names (bp)</th></tr></thead><tbody><tr><td>Tetur-VATP-F-A</td><td>CAGTTCTCCGAACCGGTAAA</td><td>A (214)</td></tr><tr><td>Tetur-VATP-R-A</td><td>CCACCGGTAATGTGACTACCA</td><td></td></tr><tr><td>Tetur-VATP-F-B</td><td>CCGTGATATGGGTTACCATG</td><td>B (416)</td></tr><tr><td>Tetur-VATP-R-B</td><td>GAAGAGGTACGAAATCTGGG</td><td></td></tr><tr><td>Tetur-sc12-F</td><td>GCCCTCTCCTGGTTGTAAACTT</td><td>Control, NC (382)</td></tr><tr><td>Tetur-sc12-R</td><td>CGACCCCATCAGGCTATTGA</td><td></td></tr><tr><th colspan="3">Spider mite RT-qPCR assay primers</th></tr><tr><th>Primer names</th><th>Oligonucleotide sequence (5' to 3')</th><th>Primer efficiency</th></tr><tr><td>RP49 (tetur18g03590) F</td><td>CTTCAAGCGGCATCAGAGC</td><td>97.6%</td></tr><tr><td>RP49 (tetur18g03590) R</td><td>CGCATCTGACCCTTGAACTTC</td><td></td></tr><tr><td>VATPase qPCR F</td><td>GGGTACCATCACATTCCTCG</td><td>104.1%</td></tr><tr><td>VATPase qPCR R</td><td>AATCGGTCTGGTTTGACGAAC</td><td></td></tr><tr><th colspan="3">Transgenic plant RT-qPCR assay primers</th></tr><tr><th>Primer names</th><th>Oligonucleotide sequence (5' to 3')</th><th>Primer efficiency</th></tr><tr><td>PEX4 (AT5G25760) F</td><td>GCTCTTATCAAAGGACCTTCGG</td><td>99.2%</td></tr><tr><td>PEX4 (AT5G25760) R</td><td>CGAACTTGAGGAGGTTGCAAAG</td><td></td></tr><tr><td>VATPase qPCR F</td><td>GGGTACCATCACATTCCTCG</td><td>104.1%</td></tr><tr><td>VATPase qPCR R</td><td>AATCGGTCTGGTTTGACGAAC</td><td></td></tr></tbody></table>

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

formation of dsRNA. dsRNA was purified by phenol-chloroform extraction followed by ethanol precipitation [26]. dsRNA was dissolved in nuclease free water and quantified using Nano-Drop (Thermo Fisher Scientific, Waltham, MA).

## Delivery of dsRNA

**Leaf floating method.** Delivery of dsRNA into mites was performed as described by Kwon *et al.* [18] with minor modifications. Solutions of dsRNA for the *TuVATPase* A, B and control sequences were prepared at the concentration of 160 ng/µL. Four hundred microliters of dsRNA solution was added in each well of a 24-well polystyrene plate (Corning, NY) and bean leaf discs (10-mm diameter) were floated on its surface. Newly-hatched larvae or newly-molted adult females were inoculated onto the floated discs (5 larvae/disc or 1 adult/disc). After inoculation, the plate was covered with a lid and maintained at L:D 16:8 and 26°C. Survival and development of the larvae and survival and fecundity of adults were recorded over 8 and 5 days, respectively. For RT-qPCR analysis, the adult females were collected into a 1.5-mL tube (at least 30 adults/tube) after 3 days of treatment. The collected samples were frozen in liquid nitrogen and stored at -80°C until RNA extraction. The biological assays on larvae and adults were conducted on 5 independent experimental runs, and the RT-qPCR analysis on 3 independent experimental runs. No leaf withering was observed after 8 days of floating on the dsRNA solution.

**Leaf coating method.** Six microliters of dsRNA solution (160 ng/µL), transcribed from the *TuVATPase* A, B or control sequence, were evenly distributed onto the upper surface of bean leaf discs (10-mm diameter). The surface of dsRNA-coated leaf discs was allowed to dry and were subsequently placed on water-soaked cotton in a polystyrene cup (94 mm in diameter, 57 mm in depth) incubated at 26°C, L:D 16:8 and 50% RH. Newly-hatched larvae (5 larvae/disc) or newly-molted adult females (1 adult/disc) were inoculated onto the dsRNA-coated discs. Survival and development of the larvae, and the survival and fecundity of the
adults were observed for 10 days. For RT-qPCR analysis, the adult females were collected into a 1.5 mL tube at day 5 after transfer to dsRNA-coated leaf discs. The collected mites were frozen in liquid nitrogen and stored at -80°C until RNA extraction. The biological assays and the material collection for RT-qPCR analysis were conducted in 3 independent experimental runs.

**Transgenic plants expressing a hairpin dsRNA.** With the aim to induce RNAi through mite feeding on transgenic plants expressing a hairpin dsRNA, the fragment A (Fig 1), of the *TuVATPase* transgene was cloned as inverted copies into the pAGRIKOLA binary vector carrying the *bar* gene coding for BASTA resistance, and transformed into *Arabidopsis thaliana* plants [27]. Wild-type Kondara plants were transformed via the *Agrobacterium* floral dip method [28]. Seeds collected from the infiltrated plants were germinated and grown under a short-day photoperiod (L:D 10:14) with a photosynthetic photon flux density of 100–150 µmol m⁻² s⁻¹ at 24°C and 55% RH and the transformed seedlings were selected by daily spraying the BASTA herbicide (10 mg/L). Six independent T1 seedlings were self-pollinated. The collected T2 progenies were analyzed for the BASTA resistance to discriminate lines that segregate BASTA resistance as single or multiple genetic loci. Representative lines that segregated BASTA resistance as single (line 2–3) and multiple (line 2–1) genetic traits were used in further experiments. The progeny seedlings of the T2 2–3 line were subjected to BASTA treatment to select homozygous seed stocks with 100% resistant phenotype. The progeny of the T2 2–1 line was selected for BASTA resistance prior to experiments, to ensure the presence of a transgene in each experimental plant.

Detached youngest adult fully-expanded leaves were collected from 3-week-old transgenic and control (untransformed) *Arabidopsis* plants, and were infested with either 10 newly-hatched larvae (for determination of the developmental rate and mortality; total of 70 larvae/ experimental run) or with one mated newly-emerged adult female mite (for determination of fecundity over 10-day period; total of 10 females/experimental run). These experiments were conducted in 3 independent experimental repeats. Detached leaves were replaced every other day. We used a London mite population adapted to *Arabidopsis* for these experiments, in order to enhance mite feeding.

For RT-qPCR analysis of *TuVATPase* transcript level, detached *Arabidopsis* leaves collected from either transgenic or control plants were infested with about 250 newly-hatched larvae. About 150 adult spider mites (collected in 3 independent batches of approximately 50) were recovered from a detached leaf setup with an aspirator after 2 weeks. Mites were flash frozen in liquid nitrogen and stored in -80°C until RNA extraction.

**Artificial diet.** The artificial diet was formulated as previously described [23,29,30]. Parafilm® M hemispheres were assembled with a custom built vacuum device as described in Suzuki *et al.* [23] and Jonckheere *et al.* [30], and 40 µL of diet was added in each hemisphere before sealing with Scotch tape. Individual hemispheres were placed in wells of the 24-well culture plate prefilled with 1 mL of water to create a water barrier preventing mite escapes. Diet was supplemented with 160 ng/µL of the dsRNA fragments (*TuVATPase* fragment A, B or control). Newly-hatched larvae or newly-emerged female adult mites were placed on individual diet-filled Parafilm® M hemisphere, and incubated at L:D 16:8, 26°C and 50% RH. Larval and adult survivorship was observed over 5 days in 3 independent experimental runs.

For RT-qPCR analysis, 100 newly-emerged female adult mites were fed on 4 Parafilm® M hemispheres filled with the artificial diet for 5 days. Mites were collected and placed into a 1.5-mL tube for subsequent RNA extraction. The collection was performed in 3 independent experimental runs.

**Soaking.** Twenty five newly-hatched larvae or 50 newly-emerged adult females were soaked in 25 µL and 50 µL of dsRNA solution (160 ng/µL; 0.1% v/v Tween 20) respectively, and incubated at 20°C for 2 h (larvae) or 24 h (adults) according to Suzuki *et al.* [23]. After
soaking, mites were transferred onto bean leaf discs (10-mm diameter; 5 larvae/disc or 1 adult/disc) placed on water-soaked cotton in the cup with a polyethylene lid with 4 venting holes each covered with a gas-permeable filter (0.45 micron pore size; Milliseal; EMD Millipore, Billerica, MA), and incubated at L:D 16:8, 26°C and 50% RH. Survival of larvae and adults, larval development and fecundity of the adult females were recorded over 10 days. The biological assays were conducted in 3 independent experimental runs.

For RT-qPCR analysis, the adult females were collected into a 1.5-mL tube (at least 30 adults/tube) at 5 days after 24-h soaking. The collected samples were frozen in liquid nitrogen and stored at -80°C until RNA extraction. The collection and the RT-qPCR analysis were conducted in 3 independent experimental runs.

## Analysis of RNAi efficiency in inbred lines

Newly-emerged adult females from eleven inbred lines generated from the reference London population were soaked in a 1.5-mL tube (60–80 adults/tube) with 50 µL of dsRNA solution (160 ng/µL; 0.1% v/v Tween 20). Adult females soaked in the dsRNA solution were incubated at 20°C in a water bath for 24 h. After soaking, mites were washed in 100 µL of double distilled water and transferred onto a bean leaf that was placed on top of water-soaked cotton in a cup with vented lid, and incubated at 26°C, L:D 16:8 and 50% RH. After 5 days, mites with dark and normal body color were evaluated for mortality; surviving mites from each phenotypic group were separated and placed on new fresh bean leaves. After 3 additional days, the fecundity from each group was evaluated. Each experimental run was performed with 4 replicates and the reference London population was used as a positive control.

## RT-qPCR analysis

Plant total RNA was extracted from approximately 100 mg of leaf tissue harvested from 3-week-old transgenic and control *Arabidopsis* plants, with the RNeasy Plant Mini Kit, including a DNase treatment (Qiagen). To assess the stability of transgene expression, RNA samples were collected from 3 leaves at 1, 4 and 7 days after they were detached from the original plant.

Mite total RNA was extracted with the RNeasy Mini Kit, including a DNase treatment (Qiagen). Two micrograms of total RNA was reverse transcribed with the Maxima First Strand cDNA Synthesis Kit for RT-qPCR (Thermo Fisher Scientific).

qPCR reactions were performed in 3 technical replicates for each sample with the Maxima SYBR Green/ROX qPCR Master Mix (Thermo Fisher Scientific). The RT-qPCR was performed on an Agilent Mx3005P qPCR instrument (Agilent Technologies, Santa Clara, CA). The reference plant gene was *PEX4* (AT5G25760), coding for a ubiquitin conjugating enzyme [31], while *RP49* (tetur18g03590), a ribosomal protein reference gene [32] was the reference gene for mite samples. Primer sequences and amplification efficiencies (E) are listed in Table 1. Cycle Threshold (Ct) values from 3 technical replicates were averaged to calculate the Ct value of each independent experimental run [33]. For plotting, expression value for each target gene (T) was normalized to the reference gene (R) and normalized relative quantity (NRQ) was calculated as follows: $NRQ = (1+E_R)^{CtR}/(1+E_T)^{CtT}$. NRQ values were then normalized to a mean of those in the control. Differences in the mean of normalized NRQ values between the control and treatments were analyzed with the Dunnett's test (function glht, R package multcomp) in the R 3.2.5 software (R Core Team 2016).

## Data analysis of survival and life history parameters

Survival curves were calculated with the Kaplan-Meier method (function survfit, R package survival) with comparisons performed based on the log-rank test (function survdiff, R package
survival). Results for the developmental times and fecundity are displayed as box-plots where central lines (second quartile, Q2) indicate the median of data, the distance between the box bottom (first quartile, Q1) and top (third quartile, Q3) indicate interquartile ranges (IQRs), and the whisker bottom and top indicate the minimum and maximum of data (except outliers that are outside the range between the lower (Q1-1.5×IQR) and upper limits (Q3+1.5×IQR) that are plotted as a white circle). Significant differences in the median developmental time and the median number of eggs laid between the control and other treatments were analyzed with the Wilcoxon-Mann-Whitney test (function wilcox.exact, R package exactRankTests). A significant difference in the proportion of color phenotype (normal or dark-body) was analyzed with the Fisher's exact test (function fisher.test). For multiple comparisons based on the paired tests, the level of significance (α) was adjusted with the Bonferroni correction (α/k, where k is the number of pairs in the multiple comparison). Analysis was performed with the R 3.2.5 software (R Core Team 2016).

## Data analysis of *TuVATPase* RNAi effect on inbred lines

RNAi was induced by soaking newly-emerged adult females as described in the “Soaking” section. Number of alive and dead individuals, and phenotype (normal or dark-body) was assessed 5 days post soaking. Fecundity of surviving females was assessed over 3 days (days 5-8 after soaking). To build the heatmap, dead/alive and normal/dark-body counts were converted to proportion of total mites recovered after soaking procedure and hierarchical clustering analysis was performed based on Euclidean distance and with the average clustering method [34]. Fecundity data were scaled to the 0-1 range and the heatmap was organized according to survival and phenotype clusters. The Wilcoxon-Mann-Whitney test was applied to analyze the differences in mortality and normalized fecundity between phenotype and treatment classes. P-values were adjusted for multiple testing with Bonferroni correction and *k* = 9. Analysis was performed with the R 3.2.5 software (R Core Team 2016).

## Results

### Induction of RNAi in the leaf floating setup

Delivery of dsRNA to mites through floating leaf discs is a current benchmark for RNAi in *T. urticae* [18]. This method and dsRNA directed against the *TuVATPase-B* sequence thus served as our reference baseline in this comparative analysis of alternative dsRNA delivery methods. The effects of dsRNAs were assessed by following the expression of the target gene together with indicators of mite performance: larval and adult survivorship, larval and adult developmental timing, and fecundity (Fig 2). As over half of the adults escaped from the floated discs and drowned in the surrounding solution at the end of the treatment (>6 days), the survival curves and fecundity measurements were limited to the first 5 days, only including mites that continuously fed on the floating discs during this period.

Although some of the measured traits tended to decrease, no significant difference was observed in the larval or adult performance parameters when *T. urticae* fed on leaf discs floating on the solution containing dsRNA-*TuVATPase*-A or dsRNA-*TuVATPase*-B, relative to the dsRNA-NC control, (Fig 2A-2D). Furthermore, no significant effects of dsRNA-*TuVATPase*-A or dsRNA-*TuVATPase*-B were observed on the relative level of *TuVATPase* expression (Fig 2E). Thus, the effects of dsRNAs homologous to the *TuVATPase* transcripts and supplied to the mites through floating leaf discs were non-significant and in accordance with the previously published data [18].
A

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

B

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

C

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /></figure>

D

<figure><img src="image_13.png" /></figure>

E

<figure><img src="image_14.png" /></figure>

**Fig 2.** The effect of dsRNA-*TuVATPase* delivered by the leaf floating method. Larval (A) and adult female (B) survivorship when feeding on bean leaf discs floating on the solution of dsRNA-*TuVATPase*-A or
dsRNA-*TuVATPase*-B (dashed lines), or dsRNA-NC (solid line). Data were collected from 5 independent experimental runs. Survival curves were plotted using the Kaplan-Meier method and compared using the log-rank test with Bonferroni correction (no asterisks, *P* > 0.05/2). (C) Developmental times of the larval (left) and adult female (middle) and male (right) mites and (D) mite fecundity. Data were collected from 5 independent experimental runs and were compared using the Wilcoxon-Mann-Whitney test with Bonferroni correction (no asterisks, *P* > 0.05/2). (E) *TuVATPase* gene expression relative to the expression of *RP49* reference gene. For the analysis of the gene expression, data were represented as mean±SE and analysed using the Dunnett's test (no asterisks, corrected *P* > 0.05).

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

## Induction of RNAi through the artificial diet

Feeding of spider mite larvae and adults on the artificial diet supplemented with 160 ng/µL of dsRNA-*TuVATPase*-B for 5 days resulted in significantly higher mortality of mites at both developmental stages compared to the dsRNA-NC control (Fig 3A and 3B). On the other hand, the application of dsRNA-*TuVATPase*-A resulted in the significant increase of larval but not of adult mortality (Fig 3A). Ultimately, the dsRNA-*TuVATPase*-A or -B reduced larval survivorship down to 30% after a 5-day feeding. The effects of dsRNA-*TuVATPase* on mite mortality correlated with the significant reduction of *TuVATPase* transcript level in these mites, down to 75% relative to mites fed on the artificial diet supplemented with the dsRNA-NC control (Fig 3C).

## Plant transgene induced RNAi

*Arabidopsis* lines carrying a single (2–3) or multiple (2–1) transgenic loci were selected based on the segregation of the *bar* gene that confers resistance to the BASTA herbicide. This marker gene was co-integrated in *Arabidopsis* chromosomes with a transgene driving the constitutive expression of the hairpin dsRNA corresponding to the *TuVATPase*-A fragment. As expected, the dsRNA was transcribed in lines 2–1 and 2–3, but not in untransformed *Arabidopsis* control plants (Fig 4A). Since mite performance assays were done on detached leaves, we initially verified that the expression of the hairpin dsRNA-*TuVATPase*-A transgene was not affected by the leaf detachment from the plant in the course of the experiment (Fig 4B). Survivorship and transitions through developmental stages of the newly-hatched larvae were similar over a 10-day period irrespective of the leaf genotype on which they fed (Fig 4C and 4D). However, the fecundity was significantly lower in females fed on lines 2–1 and 2–3, respectively, compared to those that fed on Kondara control plants (Fig 4E), even though the *TuVATPase* transcripts levels were not significantly lower in spider mites fed on transgenic lines compared to the control (Fig 4F).

## Induction of RNAi in the leaf coating setup

In leaf coating method, the solution of dsRNA is applied on the surface of the leaf disc, mimicking spray delivery common in agricultural practice. Although the foliar application of dsRNA can be effective against chewing herbivores [35], it was uncertain whether dsRNA coated on the leaf surface would be effective against mites that feed on plant mesophyll cells located within the leaf [36]. However, as shown in Fig 5A, larvae that fed on leaves coated with either dsRNA-*TuVATPase*-A or dsRNA-*TuVATPase*-B showed significantly lower survivorship relative to mites that fed on leaves coated with dsRNA-NC. In contrast to larvae fed on dsRNA-*TuVATPase* leaf discs, that had no visible phenotype, a portion of the adult mites exposed to the same dsRNAs applied to the leaf surface developed a dark-body color (Fig 5B). This phenotype was seen in 34% and 32% of mites treated with dsRNA-*TuVATPase*-A and dsRNA-*TuVATPase*-B, respectively, but was not observed in control population treated with
<figure><img src="image_15.png" /><figcaption>**Fig 3.** The effect of dsRNA-*TuVATPase* delivered through an artificial diet on mite survivorship and the endogenous *TuVATPase* gene expression. Larval (A) and adult female (B) survivorship during a 5-day feeding on artificial diet supplemented with dsRNA-*TuVATPase*-A (dashed line), dsRNA-*TuVATPase*-B (dashed line) or dsRNA-NC (solid line). Survival curves were plotted using the Kaplan-Meier method and compared using the log-rank test with Bonferroni correction (no asterisk, *P* > 0.05/2; *, *P* < 0.05/2; **, *P* < 0.01/2). (C) *TuVATPase* gene expression relative to the expression of *RP49* reference gene. For the analysis of the gene expression, data were represented as mean±SE and analysed using the Dunnett's test (**, corrected *P* < 0.01).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0180654.g003
A

<figure><img src="image_16.png" /></figure>

B

<figure><img src="image_17.png" /></figure>

C

<figure><img src="image_18.png" /></figure>

<figure><img src="image_19.png" /></figure>

D

<figure><img src="image_20.png" /></figure>

E

<figure><img src="image_21.png" /></figure>

F

<figure><img src="image_22.png" /></figure>

**Fig 4.** Mite performance and the endogenous *TuVATPase* gene expression upon mite feeding on transgenic plants expressing dsRNA-*TuVATPase* transgene. (A) Levels of dsRNA-*TuVATPase*
transgene expression in 2 independent transgenic lines of *Arabidopsis* (accession Kondara), line 2–1 (carrying multiple transgene insertions) and line 2–3 (carrying a single transgene locus), relative to the expression of *PEX4* reference gene. Samples for the analysis of the gene expression were collected from 3 independent experimental runs. (B) Stability of dsRNA-*TuVATPase* transcripts in detached leaves of transgenic *Arabidopsis* plants (line 2–3) over 7 days. Levels of transgene expression were assessed by RT-qPCR and presented relative to the expression of *PEX4* reference gene. Samples for the analysis of the gene expression were collected from 3 independent experimental runs and data were analysed using the Dunnett's test (no asterisk, *P* > 0.05). (C) Larval survivorship over a 10-day period upon mite development on detached leaves of control (Kondara; solid line) and transgenic *Arabidopsis* lines expressing dsRNA-*TuVATPase* transgene (dashed line). Data were collected from 3 independent experimental runs. Survival curves were plotted using the Kaplan–Meier method and compared using the log-rank test with Bonferroni correction (no asterisk, *P* > 0.05/2). (D) Developmental time of adulthood in mites feeding on detached leaves of control (Kondara) and transgenic *Arabidopsis* lines expressing dsRNA-*TuVATPase* transgene. Data were collected from 3 independent experimental runs and data were analysed using the Wilcoxon–Mann–Whitney test with Bonferroni correction (no asterisk, *P* > 0.05/2). (E) Fecundity of mites feeding on detached leaves of control (Kondara) and transgenic *Arabidopsis* lines expressing dsRNA-*TuVATPase* transgene over a 10-day period. Data were collected from 3 independent experimental runs and data were analysed using the Wilcoxon–Mann–Whitney test with Bonferroni correction (**, *P* < 0.05/2; **, *P* < 0.001/2). (F) The expression of the endogenous *TuVATPase* gene relative to the expression of *RP49* reference gene. Data were represented as mean±SE and analysed using the Dunnett's test (no asterisk, corrected *P* > 0.05).

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

dsRNA-NC. To test the association between the body color phenotype and RNAi effects, we followed mite performance separately on mites with normal (labeled as *normal* in Fig 5) and dark color (labeled as *dark* in Fig 5). Survivorship of adult mites treated with dsRNA-*TuVATPase*-A (regardless of the body-color phenotype) and dsRNA-NC was similar (Fig 5C). Likewise, survivorship of mites with normal body color treated with dsRNA-*TuVATPase*-B did not significantly differ from the control. However, the survivorship of dsRNA-*TuVATPase*-B treated mites displaying a dark-body phenotype was dramatically reduced (Fig 5C). Further, the effect of dsRNA-*TuVATPase* on mite fecundity was differential in mites with normal and dark-body phenotype. dsRNA-*TuVATPase* either did not affect (dsRNA-*TuVATPase*-A) or had a modest effect (dsRNA-*TuVATPase*-B) on fecundity of mites with normal body color. But, dark-colored mites, in treatments with either dsRNA-*TuVATPase*-A or -B fragments, had dramatically reduced fecundity (Fig 5E). The relative expression level of *TuVATPase* was reduced by dsRNA-*TuVATPase* relative to the control (Fig 5F). Even though mites with dark body color had more dramatic RNAi phenotypes relative to their counterparts with normal body color, expression levels of the endogenous *TuVATPase* correlated with the effect of RNAi only when mites were treated with dsRNA-*TuVATPase*-B. However, when mites were treated with dsRNA-*TuVATPase*-A, the reduction of the relative expression level in dark mites was comparable to that in mites with normal color.

## Induction of RNAi in soaked mites

In the methods described above, dsRNA intake into mites occurred through feeding. In parallel, we developed a direct dsRNA delivery method in which mites are soaked in the solution of dsRNA [23]. We showed that tracer dyes introduced into mites by soaking could be detected in the mite gut, and possibly in a wider domain, as a result of direct penetration [23].

Larval survivorship of mites soaked in a dsRNA-*TuVATPase*-A or dsRNA-*TuVATPase*-B solution could not be distinguished from the control (Fig 6A). However, over 37% of the adults soaked in dsRNA-*TuVATPase*-A and 47% of adults soaked in dsRNA-*TuVATPase*-B developed the dark-body phenotype previously observed in mites that fed on bean leaf discs coated with a dsRNA-*TuVATPase* solution (Figs 5B and 6B). Adult mites with normal body color phenotype (*normal* in Fig 6) soaked in dsRNA-*TuVATPase*-A or dsRNA-*TuVATPase*-B had the same survivorship as mites soaked in the solution containing the control dsRNA-NC (Fig 6C). Although
<figure><img src="image_23.png" /><figcaption>**Fig 5.** The effect of dsRNA-*TuVATPase* delivered by the leaf coating method. (A) Larval survivorship when feeding on bean leaf discs coated with the solution of dsRNA-*TuVATPase*-A, dsRNA-*TuVATPase*-B (dashed lines) or dsRNA-NC (solid line). Data were collected from 3 independent experimental runs. Survival curves were plotted using the Kaplan-Meier method and compared using the log-rank test with Bonferroni correction (***, *P* < 0.001/2). (B) Normal and dark-body female mites, and frequency of phenotypes observed after feeding on leaf discs coated with the solution of dsRNA-*TuVATPase*-A, dsRNA-*TuVATPase*-B or dsRNA-NC. The bar graph represents dark-body mite frequencies collected from 3 independent experimental runs. (C) Survivorship of adult female mites upon feeding on bean leaf discs coated with the solution of dsRNA-*TuVATPase*-A, dsRNA-*TuVATPase*-B (dashed lines) or dsRNA-NC (solid line) for normal and dark-body mites respectively. Data were collected from 3 independent experimental runs. Survival curves were plotted using the Kaplan-Meier method and compared using the log-rank test with Bonferroni correction (no asterisk, *P* > 0.05/4; **, *P* < 0.01/4). (D) Developmental times of the larval and adult female and male mites. Data were collected from 3 independent experimental runs and were compared using the Wilcoxon-Mann-Whitney test with Bonferroni correction (no asterisk, *P* > 0.05/2; *, *P* < 0.05/2; **, *P* < 0.01/2). (E) Normal and dark-body mite fecundity. Data were collected from 3 independent experimental runs and were compared using the Wilcoxon-Mann-Whitney test with Bonferroni correction (no asterisk, *P* > 0.05/4; **, *P* < 0.01/4). (F) *TuVATPase* gene expression relative to the expression of *RP49* reference gene in normal and dark-body female mites. Data were represented as mean±SE and analysed using the Dunnett's test (*, corrected *P* < 0.05; ***, corrected *P* < 0.001).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0180654.g005

mites with the dark-body phenotype (*dark* in Fig 6) soaked with the dsRNA-*TuVATPase*-A fragment had a slightly lower survival than the control, the difference was not significant (Fig 6C). However, mites with the dark-body phenotype soaked with the dsRNA-*TuVATPase*-B fragment had a significantly lower survival than the control (Fig 6C). Developmental times of
<figure><img src="image_24.png" /><figcaption>**Fig 6.** The effect of dsRNA-*TuVATPase* delivered by the soaking method. (A) Larval survivorship after soaking treatment in solution of dsRNA-*TuVATPase*-A, dsRNA-*TuVATPase*-B (dashed lines) or dsRNA-NC (solid line). Data were collected from 3 independent experimental runs. Survival curves were plotted using the Kaplan–Meier method and compared using the log-rank test with Bonferroni correction (no asterisk, *P* > 0.05/2). (B) Normal and dark-body female mites, and frequency of phenotypes observed after application of dsRNA-*TuVATPase*-A, dsRNA-*TuVATPase*-B or dsRNA-NC through soaking. The bar graph represents dark-body mite frequencies collected from 3 independent experimental runs. (C) Adult survivorship after soaking treatment in solution of dsRNA-*TuVATPase*-A, dsRNA-*TuVATPase*-B or dsRNA-NC separately for normal and dark-body mites. Survival curves were plotted using the Kaplan–Meier method and compared using the log-rank test with Bonferroni correction (no asterisk, *P* > 0.05/4; ***, *P* < 0.001/4). (D) Developmental times of the larval and adult female and male mites. Data were collected from 3 independent experimental runs and were compared using the Wilcoxon–Mann–Whitney test with Bonferroni correction (no asterisk, *P* > 0.05/2). (E) Normal and dark-body mite fecundity after soaking in solution of dsRNA-*TuVATPase*-A, dsRNA-*TuVATPase*-B or dsRNA-NC. Data were collected from 3 independent experimental runs and were compared using the Wilcoxon–Mann–Whitney test with Bonferroni correction (**, *P* < 0.01/4; ***, *P* < 0.001/4). (F) *TuVATPase* gene expression relative to the expression of *RP49* reference gene in normal and dark-body female mites. Data were represented as mean±SE and analysed using the Dunnett's test relative to dsRNA-NC treatment (*, corrected *P* < 0.05; **, *P* < 0.01; ***, *P* < 0.001).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0180654.g006

the larvae, females or males were unaffected by treatment (Fig 6D). The fecundity was significantly affected by mite soaking in the solution containing the dsRNA-*TuVATPase* (Fig 6E) and was more prominently reduced in mites with dark than with normal body color. In addition, the endogenous *TuVATPase* transcripts were significantly reduced in the dsRNA-*TuVATPase* treated mites (Fig 6F). Like with mites that received dsRNA through leaf coating, the severity of
RNAi phenotypes correlated with the residual levels of the endogenous *TuVATPase* only upon treatment with one of the dsRNA fragments, the dsRNA-*TuVATPase*-A (Fig 6F).

## Differential RNAi response in *T. urticae* inbred lines

Application of dsRNA-*TuVATPase* by leaf coating and soaking resulted in the phenotypic change of mite body color that correlated with the reduced mite fitness (survivorship, fecundity). Given that changes in body color occurred only in mites treated with dsRNA-*TuVATPase* and not when mites were soaked in the control dsRNA-NC fragment, the observed phenotype is not a general response to dsRNA, but a specific response to the RNA interference with *TuVATPase* expression.

These assays involved mite larvae and adults that were tightly synchronized in their development (within a 3-hour interval) originating from the London mite population whose genetic variability is narrow. Nevertheless, the strong RNAi effects (characterized by dark-body color, and significant increase of mortality and reduced fecundity) were observed only in a portion of the treated mites. When dsRNA-*TuVATPase* was coated on the bean leaf, around 30% of adults showed the strong RNAi effects, while mite soaking in the solution of dsRNA-*TuVATPase* increased the proportion of the affected adults to 50%. Such partial penetrance of the RNAi-induced phenotypes may be explained in 2 ways: a) it is an intrinsic property of the delivery methods; or b) genetic factors control the RNAi response and are variable in the London mite population.

To test the latter hypothesis, we analyzed the RNAi response of thirteen independent inbred (over 6 generations) lines derived from the London population. For each inbred line, newly-emerged adult female mites were soaked in a solution containing either the dsRNA-*TuVATPase*-B or dsRNA-NC fragment. Mites with normal or dark-body phenotypes were then counted and the fitness (mortality and fecundity) of both classes was measured separately. Genotype/dsRNA combinations were clustered according to survivorship and body color (Fig 7).

The highest level clusters clearly partitioned treatments with dsRNA-*TuVATPase*-B (top) and the control dsRNA-NC (bottom) (Fig 7A). Within these 2 major groups, mite genotypes (London population and inbred lines) were intermixed (Fig 7A) and their performance values were similar (Fig 7B and 7C), indicating that inbred lines reconstitute genetic variability existing in the London population.

The RNAi response triggered by dsRNA-*TuVATPase*-B was characterized by dark mites that invariably had higher mortality (Fig 7B) and lower fecundity (Fig 7C) relative to control or normal mites within the same population, further strengthening the link between the silencing of the *TuVATPase* gene with the color body phenotype and performance parameters. However, mites with a normal body color (normal in Fig 7) treated with dsRNA-*TuVATPase*-B also displayed an effect of the RNAi, albeit at reduced level. Their mortality is indistinguishable from the control population (Fig 7B), however, their fecundity is significantly reduced (Fig 7C). Thus, mite populations treated with the dsRNA-*TuVATPase*-B partition in 2 phenotypic classes (dark and normal body color) that differ in the quantitative severity of the RNAi effects.

None of the dsRNA-*TuVATPase*-B treated populations displayed 100% of dark mites, indicating that the high RNAi-efficiency is not completely penetrant and that the variability is intrinsic to the treatment. Two inbred lines (L2 and L14) were peculiar outliers that clustered with the control dsRNA-NC treatments, yielding no or very few dark mites. However, despite a normal color phenotype, these mite lines had similar mortality and fecundity as other inbred lines and London population, indicating that they respond to RNAi but seemingly lack the ability to display a dark-body phenotype.
<figure><img src="image_25.png" /><figcaption>**Fig 7. RNAi response in soaked *T. urticae* inbred lines.** (A) Experimental combinations of 14 mite groups (13 inbred lines and the reference population) with 2 dsRNAs (dsRNA-*TuVATPase*-B (label VA) and dsRNA-NC control (label NC)) were classified using hierarchical average clustering and Euclidean distance based on normalized survivorship (dead/alive counts) and body color (proportion of dark mites). The mite/dsRNA combinations were tested in 8 experimental batches (indicated by colors at the right side of the panel) and analyzed together. Each run included a sample of the London reference population (label REF) and 1 or 2 inbred line samples (label L). Normalized fecundity over 3 days and total proportion of dark mites are shown in the same order as dead/alive counts but were not taken into account for clustering. In most dsRNA-NC control treatments dark-body mites were not present and their fecundity could not be measured (not applicable, NA). Fecundity for VA-REF5 samples was not recorded (NR) due to a technical issue. Data were collected from 4 independent experimental runs in 8 batches. All measurements are expressed as ratios relative to the total number of tested mites and represented accordingly to the heatmap color code (bottom). (B) Normal and dark-body mite mortality (proportion of dead mites) in inbred lines (label L) and the reference population (label REF) treated with 2 dsRNAs (dsRNA-*TuVATPase*-B (label VA) and dsRNA-NC control (label NC)). (C) Normal and dark-body mite fecundity in inbred lines and the reference population treated with dsRNA-*TuVATPase*-B (label VA) and dsRNA-NC control (label NC). Data from 8 experimental batches was analyzed together and separated at the level of mite line type</figcaption></figure>
(REF and L), treatment (NC and VA), and phenotype (normal and dark-body) classes. As in most dsRNA-NC control treatments dark-body mites were not present their mortality and fecundity could not be measured (not applicable, NA). All P-values shown are corrected for multiple comparisons.

https://doi.org/10.1371/journal.pone.0180654.g007

Discussion

Even though *T. urticae* reverse genetics protocols describing the maternal injection and oral delivery of dsRNA have been published [15,18,37], they are not suitable for high throughput screens. In this study, we tested 4 methods based on artificial diet, transgenic plants, leaf coating or soaking to deliver dsRNA into spider mites. The efficiency of these methods was compared to the floating leaf disc assay, with the same dsRNA-*TuVATPase* fragment and concentration previously reported by Kwon *et al.* [18].

Several mite fitness parameters (larval and adult survivorship, larval and adult developmental timing and fecundity) and the endogenous targeted transcript level were measured to determine RNAi efficacy. As our aim was to assess relative performance of delivery methods and RNAi efficiency, we utilized conservative methods of statistical analysis. In our hands, the 4 alternative methods yielded higher silencing and/or stronger phenotypes than the floating leaf disc assay (Figs 2–6 and 8). The superior efficiency of these methods is reflected in: a) the lower amount of dsRNA applied (80 µL per larva of dsRNA solution in the floating leaf disc assay vs. 1 to 4 µL in the others), and b) the severity of the RNAi effects on mite fitness and/or endogenous target silencing. dsRNA-*TuVATPase* applied in the floating leaf disc assay had no significant effects on any of the recorded fitness parameters (Fig 2). Induction of RNAi through the expression of the dsRNA-generating transgene also had limited effect on mites: of all parameters tested, only fecundity was significantly reduced by mite feeding on transgenic relative to the control plants (Fig 4). Lower larval feeding intensity on *Arabidopsis* and potentially lower concentration of dsRNA present in transgenic plant tissues can explain the reduced effect of transgene-expressed dsRNA. Overall, the fecundity was the most sensitive parameter indicative of RNAi that was significantly reduced upon dsRNA delivery through transgenic plants, leaf coating or soaking (Figs 4, 5 and 6; note that this parameter could not have been tested for delivery of dsRNA via artificial diet, as female mites are not fecund under this rearing method).

Our comparative analysis identified soaking and leaf coating as 2 most effective ways to induce RNAi in *T. urticae* and application of dsRNA-*TuVATPase* with either method resulted

<figure><img src="image_26.png" /><figcaption>**Fig 8.** dsRNA delivery methods for RNAi-based gene silencing in *T. urticae*. Graphical summary of required dsRNA amount, severity of phenotype, and overall efficiency for the 5 dsRNA delivery methods compared in this study.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0180654.g008
in reduction of mite survivorship, fecundity and target gene expression (Figs 5 and 6). Similar dsRNA-*VATPase* silencing phenotypes have been described across multiple experimental systems, e.g. whiteflies, western corn rootworms, western flower thrips, Colorado potato beetles, fruit flies, flour beetles, pea aphids, and tobacco hornworms [12,20–22,38,39]. They confirm the essential role of the VATPase pump in proton transport across cellular membranes of eukaryotic organisms [19]. However, in addition to these phenotypes, and unique to this study, the silencing of *TuVATPase* resulted in a dark-body phenotype (Figs 5 and 6). The effect of dsRNA-*TuVATPase* on body pigmentation may reflect distinctive features of mite physiology. It may also be a more general phenomenon that is only now observed in mites because they have semitransparent integument. In our study, the dark-body phenotype tightly correlated with dsRNA-*VATPase* treatment and high RNAi effectiveness (Figs 5–7). However, it could also be observed in aging mites in untreated mite populations (note that adult mites used in our study were newly molted adult females). Thus, the dark-body phenotype is not specific to the disruption of the *TuVATPase* gene function, but may reflect mite body color change that associates with stress that can be induced by multiple independent ways.

Only a portion of mites treated with the dsRNA-*VATPase*-B fragment developed a dark body phenotype (Figs 5–7). The incomplete phenotype penetrance is likely due to the differential ability to: a) deliver dsRNA into a portion of the treated mites, or b) to silence the target gene. We previously used tracer dyes to track the delivery of small molecules into the mite body and have observed differential dye accumulation in the treated mites [15,23]. In addition, we measured that siRNA designed to silence *Tu-Dll* induced a reduction of gene expression and caused the loss-of-function phenotype in ~60% of mite embryos that received the siRNA [15]. The partial RNAi efficiency was again observed in this study despite the homogeneity of the mite genotype (inbred lines) and the tight synchronization of mite development in the experimental population (Figs 5–7). The instability of dsRNA due to the presence of nucleases, variable expression levels, or mutations in genes encoding the RNAi machinery have been shown to affect RNAi in several experimental systems [40–42]. In addition, the sequence variability of gene regions targeted for interference and differential physiological backgrounds have been postulated as factors that affect RNAi efficiency [43]. Further studies will be required to understand the source of partial and variable RNAi efficiency in spider mites.

Our analysis further highlighted that genetic component(s) can modulate RNAi-induced responses in mites. Two out of 13 inbred lines soaked in the solution containing dsRNA-*VATPase*-B did not develop the dark-body phenotype. However, the fecundity of these mites was significantly reduced relative to the dsRNA-NC treated controls (Fig 7), indicating effective RNAi despite the inability of these mites to develop a dark-body phenotype. The dark-body phenotype in dsRNA-*VATPase*-B responsive mite populations may be used as an easily assessable trait to compare factors such as RNA chemical modifications, transfection reagents and carriers when optimizing dsRNA penetration into mites.

The analysis of the *T. urticae* genome sequence identified many genes that either belong to expanded known gene families or that are “orphan” genes, with no obvious orthologs, often restricted to certain taxons [5]. Efficient reverse genetics platforms are urgently needed to study the biological and evolutionary role of these uncharacterized sequences. In this perspective, the 4 delivery methods presented here could be adapted for the application of dsRNA as a reverse-genetics tool for spider mite that will be an important asset for both the fundamental and applied sciences.

*T. urticae* is a pest with a staggering host range and one of the pest arthropod species that is most resistant to chemical pesticides [44,45]. It is thus a prime model to study the evolution of host range, plant-herbivore interactions and mechanisms of xenobiotic resistance. RNAi itself, applied in fields or greenhouses, could prove to be a valuable and sustainable biotechnological
approach for pest control [18,35,37]. For this approach, future investigations on design of dsRNA consisting of pest-lethal sequences that are highly specific to the pest species and use of a carrier for increasing stability of dsRNA will be needed.

**Acknowledgments**

The authors are grateful to Deok Ho Kwon and Si Hyeock Lee from Seoul National University for sharing data prior to publication, and to the Globalink Research Interns Rosa Ochoa Martínez and María Cristina Meza Miranda for their help with the experiments.

**Author Contributions**

**Conceptualization:** PH VG MG.

**Data curation:** TS MAN MUE VZ VG MG.

**Formal analysis:** TS MAN MUE HN PJ NB VZ VG MG.

**Funding acquisition:** PH VG MG.

**Investigation:** TS MAN MUE HN PJ NB TR RDC.

**Methodology:** TS MAN MUE VZ VG MG.

**Project administration:** VG MG.

**Resources:** PH VG MG.

**Software:** TS MAN MUE VZ.

**Supervision:** PH VG MG.

**Validation:** TS MAN MUE HN PJ NB.

**Visualization:** TS MAN MUE HN VZ.

**Writing – original draft:** TS MAN MUE VZ PH VG MG.

**Writing – review & editing:** TS VZ PH VG MG.

**References**

1. Brusca RC, Brusca GJ. Invertebrates. 2003: Sinauer Associates.

2. Dunlop JA. Geological history and phylogeny of Chelicerata. Arthropod Struct Dev. 2010; 39: 124–42.
https://doi.org/10.1016/j.asd.2010.01.003 PMID: 20093195

3. Dunlop JA, Alberti G. The affinities of mites and ticks: a review. J Zool Syst Evol Res. 2008; 46: 1–18.
https://doi.org/10.1111/j.1439-0469.2007.00429.x

4. Halliday R, O'Connor B, Baker A. Global diversity of mites. 2000. In: Raven PH, Williams T, editors. Nature and Human Society. The Quest for a Sustainable World. Proceedings of the 1997 Forum on Bio-diversity. Washington DC: National Academy Press. pp. 192–203.

5. Grbic M, Van Leeuwen T, Clark RM, Rombauts S, Rouze P, Grbic V, et al. The genome of *Tetranychus urticae* reveals herbivorous pest adaptations. Nature. 2011; 479: 487–92. https://doi.org/10.1038/nature10640 PMID: 22113690

6. Grbic M, Khila A, Lee KZ, Bjelica A, Grbic V, Whistlecraft J, et al. Mity model: *Tetranychus urticae*, a candidate for chelicerate model organism. Bioessays. 2007; 29: 489–96. https://doi.org/10.1002/bies.20564 PMID: 17450600

7. Zhai Z, Sooksa-nguan T, Vatamaniuk OK. Establishing RNA interference as a reverse-genetic approach for gene functional analysis in protoplasts. Plant Physiol. 2009; 149: 642–52. https://doi.org/10.1104/pp.108.130260 PMID: 19005083
8. Cullen LM, Arndt GM. Genome-wide screening for gene function using RNAi in mammalian cells. *Immuno Cell Biol.* 2005; 83: 217–23. https://doi.org/10.1111/j.1440-1711.2005.01332.x PMID: 15877598

9. Kuttenkeuler D, Boutros M. Genome-wide RNAi as a route to gene function in *Drosophila*. *Brief Funct Genomic Proteomic*. 2004; 3: 168–76. https://doi.org/10.1093/bfgp/3.2.168 PMID: 15355598

10. Schmitt-Engel C, Schultheis D, Schwirz J, Ströhlein N, Troelenberg N, Majumdar U, et al. The iBeetle large-scale RNAi screen reveals gene functions for insect development and physiology. *Nat Commun.* 2015; 6: 7822. https://doi.org/10.1038/ncomms8822 PMID: 26215380

11. Van BE, Kubler L, Raemaekers R, Bogaert T, Plaetinck G. Methods for controlling pests using RNAi. 2011. *Google Patents*.

12. Baum JA, Bogaert T, Clinton W, Heck GR, Feldmann P, Ilagan O, et al. Control of coleopteran insect pests through RNA interference. *Nat Biotechnol.* 2007; 25: 1322–6. https://doi.org/10.1038/nbt1359 PMID: 17982443

13. Maeda I, Kohara Y, Yamamoto M, Sugimoto A. Large-scale analysis of gene function in *Caenorhabditis elegans* by high-throughput RNAi. *Curr Biol.* 2001; 11: 171–6. https://doi.org/10.1016/S0960-9822(01) 00052-5 PMID: 11231151

14. Tabara H. RNAi in *C. elegans*: soaking in the genome sequence. *Science*. 1998; 282: 430–1. https:// doi.org/10.1126/science.282.5388.430 PMID: 9841401

15. Khila A, Grbic M. Gene silencing in the spider mite *Tetranychus urticae*: dsRNA and siRNA parental silencing of the *Distal-less* gene. *Dev Genes Evol.* 2007; 217: 241–51. https://doi.org/10.1007/s00427- 007-0132-9 PMID: 17262226

16. Liu PZ, Kaufman TC. even-skipped is not a pair-rule gene but has segmental and gap-like functions in *Oncopeltus fasciatus*, an intermediate germband insect. *Development*. 2005; 132: 2081–92. https:// doi.org/10.1242/dev.01807 PMID: 15788450

17. Bucher G, Scholten J, Klingler M. Parental RNAi in *Tribolium* (Coleoptera). *Curr Biol.* 2002; 12: R85–R6. https://doi.org/10.1016/s0960-9822(02)00666-8 PMID: 11839285

18. Kwon DH, Park JH, Lee SH. Screening of lethal genes for feeding RNAi by leaf disc-mediated systematic delivery of dsRNA in *Tetranychus urticae*. *Pestic Biochem Physiol.* 2013; 105: 69–75. https://doi. org/10.1016/j.pestbp.2012.12.001 PMID: 24238293

19. Finbow ME, Harrison MA. The vacuolar H+-ATPase: a universal proton pump of eukaryotes. *Biochem J.* 1997; 324 (Pt 3): 697–712. https://doi.org/10.1042/bj3240697

20. Zhu F, Xu J, Palli R, Ferguson J, Palli SR. Ingested RNA interference for managing the populations of the Colorado potato beetle, *Leptinotarsa decemlineata*. *Pest Manage Sci.* 2011; 67: 175–82. https:// doi.org/10.1002/ps.2048 PMID: 21061270

21. Whyard S, Singh AD, Wong S. Ingested double-stranded RNAs can act as species-specific insecticides. *Insect Biochem Mol Biol.* 2009; 39: 824–32. https://doi.org/10.1016/j.ibmb.2009.09.007 PMID: 19815067

22. Upadhyay SK, Chandrashekar K, Thakur N, Verma PC, Borgio JF, Singh PK, et al. RNA interference for the control of whiteflies (*Bemisia tabaci*) by oral route. *J Biosci* (Bangalore). 2011; 36: 153–61. https://doi.org/10.1007/s12038-011-9009-1

23. Suzuki T, Urizarna M, Nunes A, Zhurov V, Osakabe M, Van Leeuwen T, et al. Protocols for the delivery of small molecules to two-spotted spider mite, *Tetranychus urticae*. *PLoS ONE*. 2017. https://doi.org/ 10.1371/journal.pone.0180658

24. Sterck L, Billiau K, Abeel T, Rouze P, Van de Peer Y. ORCAE: online resource for community annotation of eukaryotes. *Nat Methods*. 2012; 9: 1041. https://doi.org/10.1038/nmeth.2242 PMID: 23132114

25. Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. *Nucleic Acids Res.* 1997; 25: 3389–402. PMID: 9254694

26. Kirby KS. A new method for the isolation of ribonucleic acids from mammalian tissues. *Biochem J.* 1956; 64: 405–8. PMID: 13373784

27. Hilson P, Allemeersch J, Altmann T, Aubourg S, Avon A, Beynon J, et al. Versatile gene-specific sequence tags for *Arabidopsis* functional genomics: transcript profiling and reverse genetics applications. *Genome Res.* 2004; 14: 2176–89. https://doi.org/10.1101/gr.2544504 PMID: 15489341

28. Clough SJ, Bent AF. Floral dip: a simplified method for *Agrobacterium*-mediated transformation of *Arabidopsis thaliana*. *Plant J.* 1998; 16: 735–43. PMID: 10069079

29. Van Der Geest LPS, Bosse TC, Veerman A. Development of a meridic diet for the two-spotted spider mite *Tetranychus urticae*. *Entomol Exp Appl.* 1983; 33: 297–302. https://doi.org/10.1111/j.1570-7458. 1983.tb03272.x
30. Jonckheere W, Dermauw W, Zhurov V, Wybouw N, Van den Bulcke J, Villarroel CA, et al. The salivary protein repertoire of the polyphagous spider mite *Tetranychus urticae*: a quest for effectors. Mol Cell Proteomics. 2016. https://doi.org/10.1074/mcp.M116.058081 PMID: 27703040

31. Czechowski T, Stitt M, Altmann T, Udvardi MK, Scheible WR. Genome-wide identification and testing of superior reference genes for transcript normalization in *Arabidopsis*. Plant Physiol. 2005; 139: 5–17. https://doi.org/10.1104/pp.105.063743 PMID: 16166256

32. Demaeght P, Dermauw W, Tsakireli D, Khajehali J, Nauen R, Tirry L, et al. Molecular analysis of resistance to acaricidal spirocyclic tetronic acids in *Tetranychus urticae*: CYP392E10 metabolizes spirodiclofen, but not its corresponding enol. Insect Biochem Mol Biol. 2013; 43: 544–54. https://doi.org/10.1016/j.ibmb.2013.03.007 PMID: 23523619

33. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2^−ΔΔCT^ method. Methods. 2001; 25: 402–8. https://doi.org/10.1006/meth.2001.1262 PMID: 11846609

34. Everitt BS, Landau S, Leese M, Stahl D. Cluster analysis. 2011: Wiley. 352 p.

35. San Miguel K, Scott JG. The next generation of insecticides: dsRNA is stable as a foliar-applied insecticide. Pest Manage Sci. 2016; 72: 801–9. https://doi.org/10.1002/ps.4056 PMID: 26097110

36. Bensoussan N, Santamaria ME, Zhurov V, Diaz I, Grbic M, Grbic V. Plant-herbivore interaction: dissection of the cellular pattern of *Tetranychus urticae* feeding on the host plant. Front Plant Sci. 2016; 7. https://doi.org/10.3389/fpls.2016.01105 PMID: 27512397

37. Kwon DH, Park JH, Ashok PA, Lee U, Lee SH. Screening of target genes for RNAi in *Tetranychus urticae* and RNAi toxicity enhancement by chimeric genes. Pestic Biochem Physiol. 2015. https://doi.org/10.1016/j.pestbp.2015.11.005 PMID: 27155477

38. Thakur N, Upadhyay SK, Verma PC, Chandrashekar K, Tuli R, Singh PK. Enhanced whitefly resistance in transgenic tobacco plants expressing double stranded RNA of v-ATPase A gene. PLoS ONE. 2014; 9: e87235. https://doi.org/10.1371/journal.pone.0087235 PMID: 24595215

39. Badillo-Vargas IE, Rotenberg D, Schneweis BA, Whitfield AE. RNA interference tools for the western flower thrips, *Frankliniella occidentalis*. J Insect Physiol. 2015; 76: 36–46. https://doi.org/10.1016/j.jinsphys.2015.03.009 PMID: 25796097

40. Kennedy S, Wang D, Ruvkun G. A conserved siRNA-degrading RNase negatively regulates RNA interference in *C. elegans*. Nature. 2004; 427: 645–9. https://doi.org/10.1038/nature02302 PMID: 14961122

41. Garbutt JS, Belles X, Richards EH, Reynolds SE. Persistence of double-stranded RNA in insect hemolymph as a potential determiner of RNA interference success: evidence from *Manduca sexta* and *Blattella germanica*. J Insect Physiol. 2013; 59: 171–8. https://doi.org/10.1016/j.jinsphys.2012.05.013 PMID: 22664137

42. Wang K, Peng Y, Pu J, Fu W, Wang J, Han Z. Variation in RNAi efficacy among insect species is attributable to dsRNA degradation in vivo. Insect Biochem Mol Biol. 2016; 77: 1–9. https://doi.org/10.1016/j.ibmb.2016.07.007 PMID: 27449967

43. Chu CC, Sun W, Spencer JL, Pittendrigh BR, Seufferheld MJ. Differential effects of RNAi treatments on field populations of the western corn rootworm. Pestic Biochem Physiol. 2014; 110: 1–6. https://doi.org/10.1016/j.pestbp.2014.02.003 PMID: 24759044

44. Migeon A, Dorkeld F. Spider mites web: a comprehensive database for the Tetranychidae. 2006–2017. url: www.montpellier.inra.fr/CBGP/spmweb

45. Whalon M, Mota-Sanchez D, Hollingworth R. Arthropod pesticide resistance database. 2000–2017. url: www.pesticideresistance.org.