- Open Access
Larval application of sodium channel homologous dsRNA restores pyrethroid insecticide susceptibility in a resistant adult mosquito population
© The Author(s). 2016
- Received: 7 January 2016
- Accepted: 7 June 2016
- Published: 14 July 2016
Mosquitoes host and pass on to humans a variety of disease-causing pathogens such as infectious viruses and other parasitic microorganisms. The emergence and spread of insecticide resistance is threatening the effectiveness of current control measures for common mosquito vector borne diseases, such as malaria, dengue and Zika. Therefore, the emerging resistance to the widely used pyrethroid insecticides is an alarming problem for public health. Herein we demonstrated the use of RNA interference (RNAi) to increase susceptibility of adult mosquitoes to a widely used pyrethroid insecticide.
Experiments were performed on a field-collected pyrethroid resistant strain of Ae. aegypti (Rio de Janeiro; RJ). Larvae from the resistant Ae. aegypti population were soaked with double-stranded RNAs (dsRNAs) that correspond either to voltage-gate sodium channel (VGSC), P-glycoprotein, or P450 detoxification genes and reared to adulthood. Adult mortality rates in the presence of various Deltamethrin pyrethroid concentrations were used to assess mosquito insecticide susceptibility.
We characterized the RJ Ae. aegypti strain with regard to its level of resistance to a pyrethroid insecticide and found that it was approximately 6 times more resistant to Deltamethrin compared to the laboratory Rockefeller strain. The RJ strain displayed a higher frequency of Val1016Ile and Phe1534Cys substitutions of the VGSC gene. The resistant strain also displayed a higher basal expression level of VGSC compared to the Rockefeller strain. When dsRNA-treated mosquitoes were subjected to a standard pyrethroid contact bioassay, only dsRNA targeting VGSC increased the adult mortality of the pyrethroid resistant strain. The dsRNA treatment proved effective in increasing adult mosquito susceptibility over a range of pyrethroid concentrations and these results were associated with dsRNA-specific small interfering RNAs in treated adults, and the corresponding specific down regulation of VGSC gene expression level. Finally, we demonstrated that the efficiency of our approach was further improved by ‘tiling’ along the VGSC gene in order to identify the most potent dsRNA sequences.
These results demonstrate that dsRNA applied to mosquito larvae retains its biological activity into adulthood. Thus, the RNAi system reported here could be a useful approach to control the widespread insecticide resistance in mosquitoes and other insect vectors of human diseases.
- Insecticide resistance
- Voltage-gated sodium channel
- Aedes aegypti
Mosquitoes are responsible for a huge detrimental impact on global public health. They are by far the deadliest animals to humans on the planet . Approximately 3.3 to 3.6 billion people are at risk of being infected by dengue or malaria, with around 21,000 dengue related deaths and 600,000 malaria related deaths annually . Vector control still constitutes the critical element in the current global strategies for the control of vector-borne diseases such as malaria, Yellow fever, filariasis, Japanese encephalitis, West Nile fever, Rift Valley fever, Saint Louis encephalitis, chikungunya, dengue and Zika . The Zika virus belongs to the family Flaviviridae and is transmitted to humans by Aedes aegypti mosquitoes . Zika virus infection in pregnant women has been associated with various embryonic defects and abortion. Recently, two independent studies demonstrated that Zika causes placental damage as well as fetal intrauterine growth restriction, microcephaly and death in mice [5, 6]. The expanding Aedes aegypti spread and the emerging Zika virus disease are of worldwide concern.
For controlling the Zika and dengue vector Aedes (Stegomyia) aegypti (Linnaeus, 1762), several different approaches have been recently developed, such as the use of transgenic mosquitoes and Wolbachia-infected mosquitoes . However, the traditional methods, i.e. environmental sanitation measures and chemical control are still the most widely used. For dengue control, efficient programs against Ae. aegypti were implemented in the late 1940s and early 1950s, resulting in the elimination of this species in almost all of the Americas, including Brazil . Indoor residual spraying (IRS) and pyrethroid impregnated bed-nets (ITN) both targeting the Anopheles vector, led to unprecedented success in malaria prevention in Africa, preventing hundreds of millions of malaria cases over the last decade . Although this chemical control has proven effective, it is fickle. The intensive use of insecticides and the constant exposure to them is resulting in strong selection pressure and development of resistance . Indeed, recent monitoring of malaria vectors showed that resistance levels are increasing, 53 countries reported resistance to at least one insecticide, and 41 countries reported resistance to two or more insecticide classes [2, 11].
There are four main groups of neurotoxic insecticides registered for use for public health purposes, but pyrethroids are considered the most effective and safest insecticide class to combat the adult vectors [12, 13]. All pyrethroids act on the same target protein, the voltage-gated sodium channel (VGSC). This adulticide keeps VGSC neuron membranes open, stimulating the nerve cells to produce repetitive discharges which eventually cause paralysis and mosquito death . To date, two main mechanisms are involved in pyrethroid resistance: enhanced metabolic detoxification and target site insensitivity [15, 16]. Indeed, single amino acid substitutions resulting from mutations in the VGSC gene have been associated with target site insensitivity to this class of insecticides. In natural populations of Ae. aegypti the most common mutations associated with resistance are the V1016G, F1534C [10, 17–19] and in Anopheles sp. is the L1014F [20, 21].
Among the new approaches implemented for pest control, one of the most promising is RNA interference (RNAi). RNAi is an evolutionarily conserved regulatory mechanism found in almost all eukaryotic organisms, which downregulates gene expression in a sequence-specific manner. RNAi is induced by the presence of double-stranded RNA (dsRNA) that is then processed by Dicer into small interfering RNAs (siRNAs). Later, the siRNAs recruit and guide the RNA-induced silencing complex (RISC) to the complementary RNA sequence, which is then cleaved and degraded [22–25]. The recruitment of this mechanism has led to the development of new methods for insect pest control, manipulation of insect disease vectors and disease management of beneficial insects [26–33]. Importantly, large-scale field implementation of this technology has proven to be successful in honey bees .
In mosquitoes, dsRNA delivered to larvae, resulted in increased toxicity of the larvicide temephos . Additionally, dsRNA fed directly to adults achieved a gene knockdown effect at 48 h post-feeding . However, these findings were focused on delivering dsRNA and assessing gene silencing and related-phenotype in the same mosquito developmental stage. Recent observations suggest that RNAi triggered in mosquito larvae could be maintained until adulthood [37, 38]. Nevertheless, these studies did not establish molecular evidence for a RNAi mechanism that persists across developmental stages in mosquitoes.
Herein we describe a novel RNAi-based approach for vector control. Through treating Ae. aegypti mosquito larvae with dsRNA targeting the VGSC gene, we were able to induce and characterize long-lasting gene silencing, which persisted into the mosquito adult. Moreover, the treatment successfully restored pyrethroid susceptibility in a field-collected resistant strain of Ae. aegypti. This study shows that dsRNA delivered during the mosquito larval stage can be potentially useful in integrated vector control programs, by helping to break down the resistance of adult mosquitoes to the still widely used pyrethroid insecticides.
Ae. aegypti eggs of the Rockefeller strain (Rock), Rio de Janeiro strain (RJ, Nova Iguaçu, 23°S, 43°W) were obtained from the Laboratório de Fisiologia e Controle de Artrópodes Vetores of the Instituto Oswaldo Cruz (Rio de Janeiro). Adult mosquitoes were maintained as described by Rutledge et al. , with modifications. Briefly, mosquitoes were fed on a 10 % sugar solution and were reared at 26–28 °C, 70–80 % relative humidity and a 12:12 photoperiod. Adult females were fed with artificial feeders, using natural pork gut membrane and defibrinated sheep blood (Laborclin®).
Total cDNA from Ae. aegypti adults was employed as the template for synthesis of dsRNA sequences using specific primers containing T7 promoter sequences. The dsRNAs were synthetized by in vitro transcription using Megascript RNAi kit (Ambion) according to the manufacturer’s instructions. Integrity and purity of dsRNA were verified by non-denaturing 1 % agarose gel electrophoresis, and concentration was determined spectrophotometrically. We produced one dsRNA sequence targeting the Ae. aegypti CYP9J26 gene [NCBI: XM_001649047.2, VectorBase: AAEL014609-RA], one targeting the CYP9J32 gene [NCBI: XM_001653404.2, VectorBase: AAEL008846-RA], one targeting P-glycoprotein gene [NCBI: XM_001654442.1, VectorBase: AAEL010379-RA], and sequences targeting the VGSC gene [NCBI: KC107440.1, VectorBase: AAEL006019-RD]: the VGSC422 dsRNA and 8 different tiling sequences (Tile 1, Tile 2, Tile 3, Tile 4, Tile 5, Tile 6, Tile 7 and Tile 8). The primer sequences, as well as the position of each dsRNA in the VGSC gene are indicated in Additional file 1. Artificial dsRNA (Random-dsRNA) sequence, which has no perfect homology greater than 19 bp with the Ae. aegypti genome, was generated by “Random DNA Sequence Generator” (http://www.faculty.ucr.edu/~mmaduro/random.htm). The Random-dsRNA’s full sequence (357 bp excluding the T7 promoters) is shared in Additional file 2. Additional dsRNA (Random, VGSC422 and VGSCT6) was purchased from AgroRNA.
DsRNA soaking assay
For each treatment, 100 third instar larvae were placed inside a 50 ml plastic cup with 50 μl of autoclaved food solution (6 g/100 ml – 40 % cat food, 40 % rabbit food, 10 % yeast, 10 % fish flakes) and the desired amount of indicated dsRNA ranging from 0.05 to 0.5 μg/μl in a total of 3 ml autoclaved water. The untreated group was kept in 3 ml autoclaved water as the control group. One hundred third instar larvae were soaked for 24 h and after that transferred to a new tray with 1 l of aged water. Larvae were allowed to develop until they became adults. All procedures and larval development were performed at 26–28 °C. Adult females were used in the adulticide bioassays as described below.
Deltamethrin (PESTANAL®; Sigma-Aldrich) diluted in acetone was used in the adulticide bioassays, performed as previously described . Bottles with acetone only served as a control. The resistance status of the RJ population (without dsRNA treatment of larvae) was measured and compared to the Rock strain using qualitative assays (diagnostic dose: 2 μg/ml; diagnostic time: 30 min). The knockdown effect was evaluated after 24 h into the pyrethroid bioassay. Populations with a mortality rate above 98 % were considered susceptible; whereas populations with mortality rates below 80 % were considered resistant. Populations with mortality rates between 80 and 98 % were considered populations with modifications in susceptibility .
For the quantitative bioassays, female adults (with dsRNA treatment of larval stage) were exposed for 30 min at different deltamethrin concentrations, ranging from 0 to 2 μg/ml, or at the indicated concentration. Three to four bottles (250 ml) were impregnated with each insecticide concentration. For each bottle, 10 to 15 unfed females (5–7 day-old maximum) were used. Females (live and dead) were immediately collected following the bioassays. After flash freezing, samples were stored at -80 °C for analysis of gene silencing.
Mating assay and F1 generation effect
Third instar larvae were treated with VGSC dsRNA or left untreated and grown until the pupal stage where the pupae were separated by sex within each treatment. The sex separation was performed at the pupal stage in order to ensure that all the females remain virgin. After their emergence, the adults were crossed as follows: X1 males untreated (water) crossed with females untreated (water); X2 males treated with dsRNA VGSC crossed with females treated with dsRNA VGSC; X3 males treated with dsRNA VGSC crossed with females untreated (water); X4 males untreated (water) crossed with females treated with dsRNA VGSC. After mating, the females within each crossing were blood fed and allowed to lay their eggs. These eggs (F1 generation) were grown until the adult stage and then were used for the pyrethroid bioassay.
Individual mosquitoes from either RJ or Rock strain were genotyped at position 1016 and 1534 of the genomic DNA using allele-specific PCR (AS-PCR). Following DNA extraction, PCR reactions were performed using GoTaq® Green Mastermix (Promega Corporation), 100 ng genomic DNA, 0.24 μM of the reverse primer and 0.12 μM of each one of the forward primer with a final volume of 12.5 μl. The amplification reactions were performed according to reference . The sequences of primers used were described by Linss et al.  and are shown in Additional file 3.
The amplified alleles were analyzed on a 10 % polyacrylamide gel. For the mutation Val1016Ile, a 78 bp band refers to the mutant allele (1016Ile) and a band with 98 bp represents the wild type allele (1016Val). For the mutation Phe1534Cys, the 110 bp band is in linkage with the mutant allele (1534Cys) and the 90 bp band is linked to the wild type allele (1534Phe).
Reverse transcription and relative mRNA quantification
VGSC expression levels were measured using Real time PCR following the 2-ΔΔCT method . Total RNA was extracted using TRIzol reagent (Invitrogen Life Technologies, Carlsbad, CA, USA), according to the manufacturer’s instructions. cDNA synthesis was performed in a final volume of 20 μl, using ImProm-II Reverse Transcriptase (Promega Corporation, Madison, WI, USA). The reaction mixture contained 2 μg of total RNA, 20 pmol of oligo dT primer (Invitrogen Life Technologies, Carlsbad, CA, USA), 40 U of RNAsin, 1 mM of dNTP mix, and 1 U of reverse transcriptase. PCR amplification and analysis were achieved using an ABI Prism 7500 sequence detector (Applied Biosystems, Foster City, CA, USA). All reactions were performed with SYBR Green PCR Master Mix (Applied Biosystems) using a 10 μl final volume in each reaction, which contained 1 μl of template cDNA, 2.5 pmol of each primer and 5 μl of SYBR Green. The cycles were processed according to the manufacturer’s instructions. Each sample was tested in duplicate and all quantifications were normalized using α-tubulin as the endogenous control. The primers used for all PCR amplifications are described in Additional file 1.
Alternatively, mRNA quantification was determined using the QuantiGene Plex 2.0 branched DNA assay. Briefly, total RNA from each sample (400 ng) was analyzed by using a Bio-Plex 200 system array reader with Luminex 100 X-MAP technology, and data was acquired and analyzed by using Bio-Plex Data Manager 5.0 software (Bio-Rad Laboratories, Hercules, CA). Individual bead-based probe sets specific to the VGSC gene and the endogenous control gene S7 were designed by Affymetrix and are available in the company’s website. Assays were performed as recommended by the manufacturer. Data are shown as relative light units specific to mRNA expression, with normalization to levels of the endogenous control gene.
Small RNA sequencing
Total RNA was extracted from a pool of 5 mosquitoes treated with Random or VGSC422 dsRNAs, using the Direct-zol™ RNA MiniPrep (Zymo Research), following the manufacturer’s instructions. RNA quality was verified using Agilent Total RNA Nano chip on Bioanalyzer 2100 instrument, and the quantification of total RNA for small RNA library preparation was performed using Qubit RNA BR Assay fluorimetric method (Invitrogen, USA).
One microgram (μg) RNA was used to generate a small RNA sequencing library using the TruSeq Small RNA Sample Prep Kit 7A (Illumina, USA), according to the manufacturer’s instructions. Briefly, T4 RNA ligase and T4 RNA ligase 2 Deletion mutant (Epicentre part #LR2D1132K, USA) were used to ligate RA5 and RA3 oligonucleotides to the 5′ and 3′ ends of the RNA, respectively. Adapter-ligated RNA was reverse-transcribed using a RTP primer and the resulting cDNA was amplified in an 11-cycle PCR that used RP1 and indexed RP1 primers. PCR products of 140–160 bp were isolated following electrophoresis through a 4 % agarose gel (Promega). Quality of the generated small RNA sequencing library was confirmed using Agilent High Sensitivity D1000 Screen Tape on TapeStation 2200 instrument. Sequencing was performed in a Hiseq 2500 (Illumina, USA) instrument, in a 100 bp single-end read run using HiSeq SBS clustering and sequencing kit v3-HS. Base-calling and generation of raw, de-multiplexed sequencing data in FASTQ format, was performed using BaseSpace (Illumina, USA) platform.
Processing and analyses of RNA sequencing data
Illumina Truseq small RNA 3′ adapters were trimmed using cutadapt  version 1.2.1. Untrimmed reads were discarded from the analysis. RNA-seq data were aligned to a broad institute assembly of Ae. aegypti genome  using STAR, an RNA-seq specific aligner, version 2.4.2a ; parameters were set for small RNA’s alignment: _outFilterMismatchNoverLmax = 0.05, outFilterMatchNmin = 16,_outFilterScoreMinOverLread = 0, outFilterMatchNminOverLread = 0, alignIntronMax = 1. Reads that were uniquely mapped to the VGSC gene were subjected to further downstream analysis.
Counting and normalization
The number of reads associated with each position along the silencing region and along the entire sodium channel mRNA were counted using coverageBed (bedtools version 2.17.0). Raw read counts were normalized to one million mapped reads per sample (i.e. RPM).
Statistical analysis was performed with GraphPad prism 6.0 (GraphPad Software Inc,San Diego, CA). The corresponding tests performed are indicated in the legend of each figure.
Characterization of pyrethroid resistance based on biological assays and molecular analysis of mutations in the VGSC
Based on WHO guidelines, the biological assay and the molecular analysis classifies the RJ field population as resistant to pyrethroid, whereas the Rock population is susceptible.
Soaking larvae with specific VGSC dsRNA results in gene silencing and increased pyrethroid susceptibility
Adult mosquitoes treated with VGSC422-dsRNA during the larval stage carry a siRNA population that corresponds only to the dsRNA target sequence
Targeting different sequences (‘tiling’) along the sodium channel gene generated a refined dsRNA with a greater pyrethroid efficacy
The emergence and geographical expansion of pyrethroid resistant mosquito populations is a great concern to vector-borne disease control measures. In the present study we report on a gene silencing approach that restores pyrethroid susceptibility of resistant Ae. aegypti mosquitoes. Specific dsRNA application to larvae downregulated the corresponding sodium channel gene in both larval and adult stages. This sustainable larva-to-adult silencing was sufficient to significantly reduce pyrethroid insensitivity of adults, demonstrating that reversal of insecticide resistance is a promising approach for pest control.
Resistance to a pyrethroid insecticide could evolve through various mechanisms including reduced uptake (behavioral resistance), enhanced metabolic detoxification or via gene modification that affects the insecticide’s affinity to- or activity on its target site [12, 52]. Indeed, genotyping a wild resistant strain (RJ) revealed two mutations in the VGSC gene that have been associated with increased resistance to pyrethroids (Fig. 1b) [10, 18]. Although our RNAi screen also included detoxification components, P450 and PgP, a significant phenotypic effect was observed only with dsRNA that targets the VGSC. This does not exclude the involvement of metabolic detoxification in our wild pyrethroid resistant strain. However, it further supports that the presence of mutated VGSC is involved in pyrethroid resistance in our field population.
It has been shown that enhanced detoxification activity is associated with increased expression of related metabolic components including the P450 enzymes [48, 54]. Interestingly, we detected enhanced VGSC expression in the resistant RJ strain compared to the sensitive Rock strain. Although such target expression variation might be irrelevant to pyrethroid susceptibility, it is tempting to propose a potential effect. At low pyrethroid concentrations, overexpression of the wild type VGSC might shift the toxin: target ratio, leaving enough pyrethroid-free channels to stabilize the sodium concentration within axonic cells. Alternatively, overexpression of VGSC could potentially compensate fitness costs derived from newly evolved structural modifications; hence, preserving housekeeping properties while acquiring new traits. Lastly, a cooperative interaction between highly expressed toxin chelating alleles and wild type functional alleles would be important when novel target gene modifications evolve. Nevertheless, in all scenarios, knockdown of VGSC is predicted to compromise pyrethroid insensitivity, as demonstrated by the RNAi experiments in this study.
The aim of this study was to provide a feasible RNAi solution that can be applied on wild pyrethroid resistant mosquito populations in the near future. To achieve this, high dsRNA efficacy at economic quantities is required. It is recognized that the sodium channel transcript variability governs its functional diversity including the emergence of insecticide resistance [55, 56]. Therefore, to maximize the RNAi effect we tiled a number of overlapping dsRNA constructs that together target about half of the VGSC transcript variants annotated in this work (Fig. 5a). Despite having only 45 % of the gene covered, following the I-score prediction of siRNA regions we were able to cover 15 out of 22 siRNA regions pointed out using the I-Score Designer software (Additional file 9). This strategy provided a refined dsRNA trigger (tile-6) that increased mortality with a three-fold decrease in dsRNA amounts compared to the primary VGSC dsRNA construct. Since tile-6 is not homologous to a unique VGSC sequence and not predicted to be enriched with high numbers of active siRNAs (Additional file 9), its relatively high potency is not well understood. We postulate that local secondary structure of target transcripts plays a role in dsRNA-derived siRNAs’ efficiency. It is acknowledged that in field applications dsRNA has to be packaged to preserve its stability, quantity and delivery. Further addressing such requirements, similar increased susceptibility to the adulticide effect occurred when dsRNA was complexed by a polyamine reagent (see Additional file 10), indicating that naked dsRNA is not crucial for efficient dsRNA delivery to Ae. aegypti larvae.
RNAi has been applied to modulate mosquito gene expression in vivo mostly in the adult stage. A few studies reported on environmental dsRNA uptake (soaking or ingestion) and systemic RNAi effect in the larval stage [35, 57]. Lately, Whyard et al. showed that treating larvae with dsRNA corresponding to testis-specific genes resulted in male sterility phenotype in adults . However, it has not been established molecularly whether RNAi that is induced in the larval stage can persist until adulthood. Our findings demonstrate an active prolonged cross-developmental stage, gene silencing that is effective for at least 35 days post-dsRNA introduction. The detection of VGSC-specific siRNAs (Fig. 4) coupled with VGSC silencing in adults (treated with dsRNA-VGSC as larvae) (Fig. 2e and 3b) is evidence that a sufficient amount of dsRNA was taken up by the larva and that active RNAi was sustained into adulthood. Interestingly, dsRNA-mediated VGSC knockdown had no apparent deleterious effect on treated larvae and adults in the absence of the pyrethroid insecticide. This implies that the VGSC gene is not a suitable RNAi target to induce mortality in Ae. aegypti possibly due to insufficient knockdown levels.
In conclusion, this work has demonstrated a novel strategy for the control of pyrethroid insecticide-resistant, disease vectoring mosquitoes. We applied dsRNA to Ae. aegypti larvae and showed long-lasting gene silencing that affects adults. The RNAi system was directed to impair resistance to pyrethroids by knocking down the well-established pyrethroid target, the voltage-gated sodium channel. Solving the pyrethroid resistance problem is of great importance as pyrethroids are the only class of insecticide approved to treat bed-nets . Moreover, the approach reported in this study has the potential to reduce the environmental use of chemicals with broad-spectrum insecticidal activity such as pyrethroids. It is proposed that the dsRNA-triggered long term RNAi could be applied into various complementary strategies to control mosquitoes at different developmental stages. Finally, it should be noted that based on the procedure reported herein, we are working on additional avenues to achieve disease vector control.
dsRNA, double-stranded RNA; IRS, indoor residual spraying; PgP, P-glycoprotein gene; RISC, RNA-induced silencing complex; RJ, Rio de Janeiro; RNAi, RNA interference; Rock, Rockefeller strain; RPM, reads per million; siRNAs, small interfering RNAs; VGSC, voltage-gate sodium channel; WT, wild type; μg, microgram
We are thankful to Helixxa - Genomics Service Provider (Campinas, Brazil) for the excellent work with miRNA-sequencing and to Dr. Vered Kunik for the all bioinformatic and RNA-sequencing analyses.
The study was supported by Forrest Innovations’ internal resources. Forrest Innovations’ funders had no role in the design of the study and collection, analysis, and interpretation of data as well as in writing the manuscript.
Availability of data and materials
The small RNA sequencing data reported in this publication have been deposited in NCBI's Gene Expression Omnibus, and are accessible through GEO Series accession number GSE84030 (http://0-www.ncbi.nlm.nih.gov.brum.beds.ac.uk/geo/query/acc.cgi?acc=GSE84030).
Conceived the study: ESB, NP, EM. Designed the experiments: ACDB, RFC, MLF, LCP, AW, JBPL, NP, ESB, JH, EM. Performed the experiments: ACDB, RFC, MLF, LCP. Analyzed the data: ACDB, RFC, MLF, LCP, AW, JBPL, NP, ESB, JH, EM. Wrote the paper: ACDB, RFC, MLF, LCP, NP, JH, EM. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Gates B. The Deadliest Animal in the World. Gates Fundation. 2014. http://www.gatesnotes.com/Health/Most-Lethal-Animal-Mosquito-Week. Accessed 15 Dec 2015.
- Organization WH. World Malaria Report. 2014. http://www.who.int/malaria/publications/world_malaria_report_2014/report/en/. Accessed 15 Dec 2015.Google Scholar
- van der Berg H, Zaim M, Yadav RS, Soares A, Ameneshewa B, Mnzava A, et al. Global trends in the use of insecticides to control vector-borne diseases. Environ Health Perspec. 2012;120:577–82.View ArticleGoogle Scholar
- Dick GW, Kitchen SF, Haddow AJ. Zika virus (I). Isolations and serological specificity. Trans R Soc Trop Med Hyg. 1952;46:509–20.View ArticlePubMedGoogle Scholar
- Miner JJ, Cao B, Govero J, Smith AM, Fernandez E, Cabrera OH, et al. Zika virus infection during pregnancy in mice causes placental damage and fetal demise. Cell. 2016;165:1081-1091.Google Scholar
- Cugola FR, Fernandes IR, Russo FB, Freitas BC, Dias JLM, Guimarães KP, Benazzato C, Almeida N, Pignatari GC, Romero S, Polonio CM, Cunha I, Freitas CL, Brandão WN, Rossato C, Andrade DG, Faria DP, Garcez AT, Buchpigel CA, Braconi CT, Mendes E, Sall AA, Zanotto PMA, Peron JPS, Muotri AR, Beltrão-Braga PCB. The Brazilian Zika virus strain causes birth defects in experimental models. Nature. 2016. doi:10.1038/nature18296.PubMedGoogle Scholar
- Maciel-de-Freitas R, Aguiar R, Bruno RV, Guimaraes MC, Lourenco-de-Oliveira R, Sorgine MH, et al. Why do we need alternative tools to control mosquito-borne diseases in Latin America? Mem Inst Oswaldo Cruz. 2012;107:828–9.View ArticlePubMedGoogle Scholar
- Soper FL. The 1964 status of Aedes aegypti erradication and yellow fever in the Americas. Am J Trop Med Hyg. 1965;14:887–91.PubMedGoogle Scholar
- Bhatt S, Weiss DJ, Cameron E, Bisanzio D, Mappin B, Dalrymple U, et al. The effect of malaria control on Plasmodium falciparum in Africa between 2000 and 2015. Nature. 2015;526:207–11.View ArticlePubMedPubMed CentralGoogle Scholar
- Linss JG, Brito LP, Garcia GA, Araki AS, Bruno RV, Lima JB, et al. Distribution and dissemination of the Val1016Ile and Phe1534Cys Kdr mutations in Aedes aegypti Brazilian natural populations. Parasit Vectors. 2014;7:25.View ArticlePubMedPubMed CentralGoogle Scholar
- Butler D. Mosquitoes score in chemical war. Nature. 2011;475:19.View ArticlePubMedGoogle Scholar
- Silva AP, Santos JM, Martins AJ. Mutations in the voltage-gated sodium channel gene of Anophelines and their association with resistance to pyrethroids - a review. Parasit Vectors. 2014;7:450.View ArticlePubMedPubMed CentralGoogle Scholar
- Brogdon WG, McAllister JC. Simplification of adult mosquito bioassays through use of time-mortality determinations in glass bottles. J Am Mosq Control. 1998;14:159–64.Google Scholar
- Davies TGE, Field LM, Usherwood PNR, Williamson MS. A comparative study of voltage-gated sodium channels in the Insecta: implications for pyrethroid resistance in Anopheline and other Neopteran species. Insect Mol Biol. 2007;16:361–75.View ArticlePubMedGoogle Scholar
- Kawada H, Higa Y, Komagata O, Kasai S, Tomita T, Thi Yen N, et al. Widespread distribution of a newly found point mutation in voltage-gated sodium channel in pyrethroid-resistant Aedes aegypti populations in Vietnam. PLoS Negl Trop Dis. 2009;3:e527.View ArticlePubMedPubMed CentralGoogle Scholar
- Martins AJ, Ribeiro CD, Bellinato DF, Peixoto AA, Valle D, Lima JB. Effect of insecticide resistance on development, longevity and reproduction of field or laboratory selected Aedes aegypti populations. PloS One. 2012;7:e31889.View ArticlePubMedPubMed CentralGoogle Scholar
- Lima EP, Paiva MH, de Araujo AP, da Silva EV, da Silva UM, de Oliveira LN, et al. Insecticide resistance in Aedes aegypti populations from Ceara. Brazil Parasit Vectors. 2011;4:5.View ArticlePubMedGoogle Scholar
- Martins AJ, Lima JB, Peixoto AA, Valle D. Frequency of Val1016Ile mutation in the voltage-gated sodium channel gene of Aedes aegypti Brazilian populations. Trop Med Int Health. 2009;14:1351–5.View ArticlePubMedGoogle Scholar
- Yanola J, Somboon P, Walton C, Nachaiwieng W, Somwang P, Prapanthadara LA. High-throughput assays for detection of the F1534C mutation in the voltage-gated sodium channel gene in permethrin-resistant Aedes aegypti and the distribution of this mutation throughout Thailand. Trop Med Int Health. 2011;16:501–9.View ArticlePubMedGoogle Scholar
- Martinez-Torres D, Chandre F, Williamson MS, Darriet F, Berge JB, Devonshire AL, et al. Molecular characterization of pyrethroid knockdown resistance (kdr) in the major malaria vector Anopheles gambiae s.s. Insect Mol Biol. 1998;7:179–84.View ArticlePubMedGoogle Scholar
- Kawada H, Futami K, Komagata O, Kasai S, Tomita T, Sonye G, et al. Distribution of a knockdown resistance mutation (L1014S) in Anopheles gambiae s.s. and Anopheles arabiensis in western and southern Kenya. PloS One. 2011;6:e24323.View ArticlePubMedPubMed CentralGoogle Scholar
- Almeida R, Allshire RC. RNA silencing and genome regulation. Trends Cell Biol. 2005;15:251–8.View ArticlePubMedGoogle Scholar
- Kim D, Rossi J. RNAi mechanisms and applications. BioTechniques. 2008;44:613–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Ding SW. RNA-based antiviral immunity. Nat Rev Immunol. 2010;10:632–44.View ArticlePubMedGoogle Scholar
- Wilson RC, Doudna JA. Molecular mechanisms of RNA interference. Annu Rev Biophys. 2013;42:217–39.View ArticlePubMedGoogle Scholar
- Maori E, Paldi N, Shafir S, Kalev H, Tsur E, Glick E, et al. IAPV, a bee-affecting virus associated with Colony Collapse Disorder can be silenced by dsRNA ingestion. Insect Mol Biol. 2009;18:55–60.View ArticlePubMedGoogle Scholar
- Garbian Y, Maori E, Kalev H, Shafir S, Sela I. Bidirectional transfer of RNAi between honey bee and Varroa destructor: Varroa gene silencing reduces Varroa population. PLoS Pathog. 2012;8:e1003035.View ArticlePubMedPubMed CentralGoogle Scholar
- Pridgeon JW, Zhao L, Becnel JJ, Strickman DA, Clark GG, Linthicum KJ. Topically applied AaeIAP1 double-stranded RNA kills female adults of Aedes aegypti. J Med Entomol. 2008;45:414–20.View ArticlePubMedGoogle Scholar
- Katoch R, Thakur N. Advances in RNA interference technology and its impact on nutritional improvement, disease and insect control in plant. Appl Biochem Biotechnol. 2013;169:1579–605.View ArticlePubMedGoogle Scholar
- Scott JG, Michel K, Bartholomay LC, Siegfried BD, Hunter WB, Smagghe G, et al. Towards the elements of successful insect RNAi. J Insect Physiol. 2013;59:1212–21.View ArticlePubMedGoogle Scholar
- Paldi N, Glick E, Oliva M, Zilberberg Y, Aubin L, Pettis J, et al. Effective gene silencing in a microsporidian parasite associated with honeybee (Apis mellifera) colony declines. Appl Environ Microbiol. 2010;76:5960–4.View ArticlePubMedPubMed CentralGoogle Scholar
- Bing ZS, Jiang P, Wang ZQ, Long SR, Liu RD, Zhang X, Yang W, Ren HJ, Cui J. DsRNA-mediated silencing of Nudix hydrolase in Trichinella spiralis inhibits the larval invasion and survival in mice. Exp Parasitol. 2016;162:35–42.View ArticleGoogle Scholar
- Pan ZH, Gao K, Hou CX, Wu P, Qin GX, Geng T, Guo XJ. dsRNA interference on expression of a RNA-dependent RNA polymerase gene of Bombyx mori cytoplasmic polyhedrosis virus. Gene. 2015;565:56–61.View ArticlePubMedGoogle Scholar
- Hunter W, Ellis J, Vanengelsdorp D, Hayes J, Westervelt D, Glick E, et al. Large-scale field application of RNAi technology reducing Israeli acute paralysis virus disease in honey bees (Apis mellifera, Hymenoptera: Apidae). PLoS Pathog. 2010;6:e1001160. http://www.who.int/malaria/publications/atoz/who_cds_cpc_mal_98_12/en/
- Figueira-Mansur J, Ferreira-Pereira A, Mansur JF, Franco TA, Alvarenga ES, Sorgine MH, et al. Silencing of P-glycoprotein increases mortality in temephos-treated Aedes aegypti larvae. Insect Mol Biol. 2013;22:648–58.View ArticlePubMedGoogle Scholar
- Coy MR, Sanscrainte ND, Chalaire KC, Inberg A, Maayan I, Glick E, et al. Gene silencing in adult Aedes aegypti mosquitoes through oral delivery of double-stranded RNA. J Appl Entomol. 2012;136:741–8.View ArticleGoogle Scholar
- Whyard S, Erdelyan CNG, Partridge AL, Singh AD, Beebe NW, Capina R. Silencing the buzz: a new approach to population suppression of mosquitoes by feeding larvae double-stranded RNAs. Parasite Vector. 2015;8:96.View ArticleGoogle Scholar
- Mysore K, Sun L, Tomchaney M, Sullivan G, Adams H, et al. siRNA-mediated silencing of doublesex during female development of the dengue vector mosquito Aedes aegypti. PLoS Negl Trop Dis. 2015;9:e0004213.View ArticlePubMedPubMed CentralGoogle Scholar
- Rutledge LC, Ward RA, Gould DJ. Studies on the feeding response of mosquitoes to nutritive solutions in a new membrane feeder. Mosq News. 1964;24:407–19.Google Scholar
- Brogdon W, Chan A. Guidelines for evaluating insecticide resistance in vectors using the CDC bottle bioassay. Methods in Anopheles Research. Second ed. 2010.Google Scholar
- Organization WH. Test procedures for insecticide resistance monitoring in malaria vectors, bio efficacy and persistence of insecticides on treated surfaces. WHO/CDS/CPC/MAL/9812. 1998. http://www.who.int/malaria/publications/atoz/who_cds_cpc_mal_98_12/en/.
- Saavedra-Rodriguez K, Urdaneta-Marquez L, Rajatileka S, Moulton M, Flores AE, Fernandez-Salas I, et al. A mutation in the voltage-gated sodium channel gene associated with pyrethroid resistance in Latin American Aedes aegypti. Insect Mol Biol. 2007;16:785–98.View ArticlePubMedGoogle Scholar
- Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001;25:402–8.View ArticlePubMedGoogle Scholar
- Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011;17:10–2.View ArticleGoogle Scholar
- National Center for Biotechnology Information, Accession: PRJNA12434, ID: 12434. https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/bioproject/12434 .
- Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29:15–21.View ArticlePubMedGoogle Scholar
- Montella IR, Martins AJ, Viana-Medeiros PF, Lima JB, Braga IA, Valle D. Insecticide resistance mechanisms of Brazilian Aedes aegypti populations from 2001 to 2004. J Trop Med Hyg. 2007;77:467–77.Google Scholar
- Stevenson BJ, Pignatelli P, Nikou D, Paine MJ. Pinpointing P450s associated with pyrethroid metabolism in the dengue vector Aedes aegypti: developing new tools to combat insecticide resistance. PLoS Negl Trop Dis. 2012;6:e1595.View ArticlePubMedPubMed CentralGoogle Scholar
- Davies TE, O'Reilly AO, Field LM, Wallace B, Williamson MS. Knockdown resistance to DDT and pyrethroids: from target-site mutations to molecular modelling. Pest Manag Sci. 2008;64:1126–30.View ArticlePubMedGoogle Scholar
- Price DRG, Gatehouse JA. RNAi-mediated crop protection against insects. Trends Biotechnol. 2008;26:393–400.View ArticlePubMedGoogle Scholar
- Hoa NT, Keene KM, Olson KE, Zheng L. Characterization of RNA interference in an Anopheles gambiae cell line. Insect Biochem Molec. 2003;33:949–57.View ArticleGoogle Scholar
- Ranson H, N'Guessan R, Lines J, Moiroux N, Nkuni Z, Corbel V. Pyrethroid resistance in African Anopheline mosquitoes: what are the implications for malaria control? Trends Parasitol. 2011;27:91–8.View ArticlePubMedGoogle Scholar
- Stenhouse SA, Plernsub S, Yanola J, Lumjuan N, Dantrakool A, Choochote W, et al. Detection of the V1016G mutation in the voltage-gated sodium channel gene of Aedes aegypti (Diptera: Culicidae) by allele-specific PCR assay, and its distribution and effect on deltamethrin resistance in Thailand. Parasit Vectors. 2013;6:253.View ArticlePubMedPubMed CentralGoogle Scholar
- David JP, Ismail HM, Chandor-Proust A, Paine MJ. Role of cytochrome P450s in insecticide resistance: impact on the control of mosquito-borne diseases and use of insecticides on Earth. Philos Trans R Soc Lond B Biol Sci. 2013;368:20120429.View ArticlePubMedPubMed CentralGoogle Scholar
- Dong K. Insect sodium channels and insecticide resistance. Invertebr Neurosci. 2007;7:17–30.View ArticleGoogle Scholar
- Dong K, Du Y, Rinkevich F, Nomura Y, Xu P, Wang L, et al. Molecular biology of insect sodium channels and pyrethroid resistance. Insect Biochem Molec Biol. 2014;50:1–17.View ArticleGoogle Scholar
- Zhang X, Zhang J, Zhu KY. Chitosan/double-stranded RNA nanoparticle-mediated RNA interference to silence chitin synthase genes through larval feeding in the African malaria mosquito (Anopheles gambiae). Cell. 2016;165:683–693.View ArticleGoogle Scholar
- Sarkies P, Miska EA. Small RNAs break out: the molecular cell biology of mobile small RNAs. Nat Rev Mol Cell Biol. 2014;15:525–35.View ArticlePubMedGoogle Scholar