Quantitative trait loci mapping of genome regions controlling permethrin resistance in the mosquito Aedes aegypti

Strong Positive Selection in Aedes aegypti and the Rapid Evolution of Insecticide Resistance

Aedes aegypti vectors the pathogens that cause dengue, yellow fever, Zika virus, and chikungunya and is a serious threat to public health in tropical regions. Decades of work has illuminated many aspects of Ae. aegypti's biology and global population structure and has identified insecticide resistance genes; however, the size and repetitive nature of the Ae. aegypti genome have limited our ability to detect positive selection in this mosquito. Combining new whole genome sequences from Colombia with publicly available data from Africa and the Americas, we identify multiple strong candidate selective sweeps in Ae. aegypti, many of which overlap genes linked to or implicated in insecticide resistance. We examine the voltage-gated sodium channel gene in three American cohorts and find evidence for successive selective sweeps in Colombia. The most recent sweep encompasses an intermediate-frequency haplotype containing four candidate insecticide resistance mutations that are in near-perfect linkage disequilibrium with one another in the Colombian sample. We hypothesize that this haplotype may continue to rapidly increase in frequency and perhaps spread geographically in the coming years. These results extend our knowledge of how insecticide resistance has evolved in this species and add to a growing body of evidence suggesting that Ae. aegypti has an extensive genomic capacity to rapidly adapt to insecticide-based vector control.

The role of voltage-gated sodium channel genotypes in pyrethroid resistance in Aedes aegypti in Taiwan

Background

Aedes aegypti is the major vector of dengue that threatens public health in tropical and subtropical regions. Pyrethroid-based control strategies effectively control this vector, but the repeated usage of the same insecticides leads to resistance and hampers control efforts. Therefore, efficient and prompt monitoring of insecticide resistance in local mosquito populations is critical for dengue control.

Methodology/Principal finding

We collected Ae. aegypti in southern Taiwan in March and October 2016. We analyzed the voltage-gated sodium channel (vgsc) genotypes of parentals (G0) and G1 adults after cypermethrin insecticide bioassay. Our results showed that four VGSC mutations (S989P, V1016G, F1534C, and D1763Y) associated with resistance were commonly detected in field-collected Ae. aegypti. The frequencies of these four mutations in the local mosquito population were significantly higher in October (0.29, 0.4, 0.27 and 0.11) than in March (0.09, 0.16, 0.18 and 0.03). Specific vgsc combined genotypes composed of the one to four such mutations (SGFY/SGFY, SVCD/SVCD, SGFY/PGFD, SVCD/SGFY, PGFD/PGFD, and SVCD/PGFD) shifted towards higher frequencies in October, implying their resistance role. In addition, the cypermethrin exposure bioassay data supported the field observations. Moreover, our study observed an association between the resistance level and the proportion of resistance genotypes in the population.

Conclusions/Significance

This is the first study to demonstrate the role of four-locus vgsc genotypes in resistance evaluation in a local Ae. aegypti population in Taiwan. This alternative method using resistance-associated genotypes as an indicator of practically insecticide resistance monitoring is a useful tool for providing precise and real-time information for decision makers.

Dengue outbreaks occur annually in Taiwan, and pyrethroid insecticides are commonly used to reduce mosquito density. Insecticide resistance of mosquitoes is commonly observed in the field and threatens vector control programs. Here, we analyzed the association between the combined vgsc genotype and resistance phenotype based on field surveillance data in March and October and a cypermethrin exposure bioassay. Resistance-attributable specific vgsc genotypes were proposed. Using the combined vgsc genotype rather than each vgsc allele is recommended for better resistance prediction to provide real-time information for control program managers.

Vector Specificity of Arbovirus Transmission

More than 25% of human infectious diseases are vector-borne diseases (VBDs). These diseases, caused by pathogens shared between animals and humans, are a growing threat to global health with more than 2.5 million annual deaths. Mosquitoes and ticks are the main vectors of arboviruses including flaviviruses, which greatly affect humans. However, all tick or mosquito species are not able to transmit all viruses, suggesting important molecular mechanisms regulating viral infection, dissemination, and transmission by vectors. Despite the large distribution of arthropods (mosquitoes and ticks) and arboviruses, only a few pairings of arthropods (family, genus, and population) and viruses (family, genus, and genotype) successfully transmit. Here, we review the factors that might limit pathogen transmission: internal (vector genetics, immune responses, microbiome including insect-specific viruses, and coinfections) and external, either biotic (adult and larvae nutrition) or abiotic (temperature, chemicals, and altitude). This review will demonstrate the dynamic nature and complexity of virus–vector interactions to help in designing appropriate practices in surveillance and prevention to reduce VBD threats.

Transcriptomic and proteomic analysis of pyrethroid resistance in the CKR strain of Aedes aegypti

Aedes aegypti is an important vector of human viral diseases. This mosquito is distributed globally and thrives in urban environments, making it a serious risk to human health. Pyrethroid insecticides have been the mainstay for control of adult A. aegypti for decades, but resistance has evolved, making control problematic in some areas. One major mechanism of pyrethroid resistance is detoxification by cytochrome P450 monooxygenases (CYPs), commonly associated with the overexpression of one or more CYPs. Unfortunately, the molecular basis underlying this mechanism remains unknown. We used a combination of RNA-seq and proteomic analysis to evaluate the molecular basis of pyrethroid resistance in the highly resistant CKR strain of A. aegypti. The CKR strain has the resistance mechanisms from the well-studied Singapore (SP) strain introgressed into the susceptible Rockefeller (ROCK) strain genome. The RNA-seq and proteomics data were complimentary; each offering insights that the other technique did not provide. However, transcriptomic results did not quantitatively mirror results of the proteomics. There were 10 CYPs which had increased expression of both transcripts and proteins. These CYPs appeared to be largely trans-regulated, except for some CYPs for which we could not rule out gene duplication. We identified 65 genes and lncRNAs as potentially being responsible for elevating the expression of CYPs in CKR. Resistance was associated with multiple loci on chromosome 1 and at least one locus on chromosome 3. We also identified five CYPs that were overexpressed only as proteins, suggesting that stabilization of CYP proteins could be a mechanism of resistance. Future studies to increase the resolution of the resistance loci, and to examine the candidate genes and lncRNAs identified here will greatly enhance our understanding of CYP-mediated resistance in A. aegypti.

A Two-Locus System with Strong Epistasis Underlies Rapid Parasite-Mediated Evolution of Host Resistance

Parasites are a major evolutionary force, driving adaptive responses in host populations. Although the link between phenotypic response to parasite-mediated natural selection and the underlying genetic architecture often remains obscure, this link is crucial for understanding the evolution of resistance and predicting associated allele frequency changes in the population. To close this gap, we monitored the response to selection during epidemics of a virulent bacterial pathogen, Pasteuria ramosa, in a natural host population of Daphnia magna. Across two epidemics, we observed a strong increase in the proportion of resistant phenotypes as the epidemics progressed. Field and laboratory experiments confirmed that this increase in resistance was caused by selection from the local parasite. Using a genome-wide association study, we built a genetic model in which two genomic regions with dominance and epistasis control resistance polymorphism in the host. We verified this model by selfing host genotypes with different resistance phenotypes and scoring their F1 for segregation of resistance and associated genetic markers. Such epistatic effects with strong fitness consequences in host–parasite coevolution are believed to be crucial in the Red Queen model for the evolution of genetic recombination.

Levels of cross-resistance to pyrethroids conferred by the Vssc knockdown resistance allele 410L+1016I+1534C in Aedes aegypti

Aedes aegypti is a primary vector of viral pathogens and is responsible for millions of human infections annually that represent critical public health and economic costs. Pyrethroids are one of the most commonly used classes of insecticides to control adult A. aegypti. The insecticidal activity of pyrethroids depends on their ability to bind and disrupt the voltage-sensitive sodium channel (VSSC). In mosquitoes, a common mechanism of resistance to pyrethroids is due to mutations in Vssc (hereafter referred as knockdown resistance, kdr). In this study, we found that a kdr (410L+V1016I+1534C) allele was the main mechanism of resistance in a pyrethroid-resistant strain of A. aegypti collected in Colombia. To characterize the level of resistance these mutations confer, we isolated a pyrethroid resistant strain (LMRKDR:RK, LKR) that was congenic to the susceptible Rockefeller (ROCK) strain. The full-length cDNA of Vssc was cloned from LKR and no additional resistance mutations were present. The levels of resistance to different pyrethroids varied from 3.9- to 56-fold. We compared the levels of resistance to pyrethroids, DCJW and DDT between LKR and what was previously reported in two other congenic strains that share the same pyrethroid-susceptible background (the ROCK strain), but carry different kdr alleles (F1534C or S989P + V1016G). The resistance conferred by kdr alleles can vary depending on the stereochemistry of the pyrethroid. The 410L+1016I+1534C kdr allele does not confer higher levels of resistance to six of ten pyrethroids, relative to the 1534C allele. The importance of these results to understand the evolution of insecticide resistance and mosquito control are discussed.

Aedes aegypti is a primary vector of viral pathogens and is responsible for millions of human infections and pyrethroids are one of the most commonly used classes of insecticides to control adult A. aegypti. The insecticidal activity of pyrethroids depends on their ability to bind and disrupt the voltage-sensitive sodium channel (VSSC). A common mechanism of resistance to pyrethroids is due to mutations in Vssc (hereafter referred as knockdown resistance mutations, kdr). In this study, we found that a kdr (410L+1016I+154C) allele was the main mechanism of resistance in a pyrethroid-resistant strain of A. aegypti collected in Colombia. We determined the levels of resistance to pyrethroids, DCJW and DDT conferred by this allele. The 410L+1016I+1534C kdr allele does not confer higher levels of resistance to six out of 10 pyrethroids relative to the 1534C allele. Resistance conferred by kdr alleles can vary between pyrethroids and can vary between the stereoisomers. Knowing the levels of resistance conferred by different alleles to different insecticides provides valuable information for pest control and for understanding the evolution of resistance.

Permethrin resistance in Aedes aegypti: Genomic variants that confer knockdown resistance, recovery, and death

Pyrethroids are one of the few classes of insecticides available to control Aedes aegypti, the major vector of dengue, chikungunya, and Zika viruses. Unfortunately, evolving mechanisms of pyrethroid resistance in mosquito populations threaten our ability to control disease outbreaks. Two common pyrethroid resistance mechanisms occur in Ae. aegypti: 1) knockdown resistance, which involves amino acid substitutions at the pyrethroid target site-the voltage-gated sodium channel (VGSC)-and 2) enhanced metabolism by detoxification enzymes. When a heterogeneous population of mosquitoes is exposed to pyrethroids, different responses occur. During exposure, a proportion of mosquitoes exhibit immediate knockdown, whereas others are not knocked-down and are designated knockdown resistant (kdr). When these individuals are removed from the source of insecticide, the knocked-down mosquitoes can either remain in this status and lead to dead or recover within a few hours. The proportion of these phenotypic responses is dependent on the pyrethroid concentration and the genetic background of the population tested. In this study, we sequenced and performed pairwise genome comparisons between kdr, recovered, and dead phenotypes in a pyrethroid-resistant colony from Tapachula, Mexico. We identified single-nucleotide polymorphisms (SNPs) associated with each phenotype and identified genes that are likely associated with the mechanisms of pyrethroid resistance, including detoxification, the cuticle, and insecticide target sites. We identified high association between kdr and mutations at VGSC and moderate association with additional insecticide target site, detoxification, and cuticle protein coding genes. Recovery was associated with cuticle proteins, the voltage-dependent calcium channel, and a different group of detoxification genes. We provide a list of detoxification genes under directional selection in this field-resistant population. Their functional roles in pyrethroid metabolism and their potential uses as genomic markers of resistance require validation.

Standing genetic variation in laboratory populations of insecticide-susceptible Phlebotomus papatasi and Lutzomyia longipalpis (Diptera: Psychodidae: Phlebotominae) for the evolution of resistance

Insecticides can exert strong selection on insect pest species, including those that vector diseases, and have led to rapid evolution of resistance. Despite such rapid evolution, relatively little is known about standing genetic variation for resistance in insecticide-susceptible populations of many species. To help fill this knowledge gap, we generated genotyping-by-sequencing data from insecticide-susceptible Phlebotomus papatasi and Lutzomyia longipalpis sand flies that survived or died from a sub-diagnostic exposure to either permethrin or malathion using a modified version of the Centers for Disease Control and Prevention bottle bioassay. Multi-locus genome-wide association mapping methods were used to quantify standing genetic variation for insecticide resistance in these populations and to identify specific alleles associated with insecticide survival. For each insecticide treatment, we estimated the proportion of the variation in survival explained by the genetic data (i.e., "chip" heritability) and the number and contribution of individual loci with measurable effects. For all treatments, survival to an insecticide exposure was heritable with a polygenic architecture. Both P. papatasi and L. longipalpis had alleles for survival that resided within many genes throughout their genomes. The implications for resistance conferred by many alleles, as well as inferences made about the utility of laboratory insecticide resistance association studies compared to field observations, are discussed.

Aedes aegypti miRNA-33 modulates permethrin induced toxicity by regulating VGSC transcripts

Aedes aegypti is a major vector of Zika, dengue, and other arboviruses. Permethrin adulticidal spraying, which targets the voltage-gated sodium channel (VGSC), is commonly done to reduce local mosquito populations and protect humans from exposure to arbovirus pathogens transmitted by this dangerous pest. Permethrin resistance, however, is a growing problem and understanding its underlying molecular basis may identify avenues to combat it. We identified a single G:C polymorphism in pre-miR-33 that was genetically associated with permethrin resistance; resulting isoforms had structural differences that may affect DICER-1/pre-miRNA processing rates. We then assessed the effects of overexpression of pre-miR-33 isoforms on permethrin toxicological phenotypes, VGSC transcript abundance and protein levels for two genetically related mosquito strains. One strain had its naturally high permethrin resistance levels maintained by periodic treatment, and the other was released from selection. VGSC protein levels were lower in the permethrin resistant strain than in the related permethrin-susceptible strain. Overexpression of the G-pre-miR-33 isoform reduced VGSC expression levels in both strains. To further elucidate changes in gene expression associated with permethrin resistance, exome-capture gDNA deep sequencing, genetic association mapping and subsequent gene set enrichment analysis revealed that transport genes, in particular, were selected in resistant versus susceptible mosquitoes. Collectively, these data indicate that miR-33 regulates VGSC expression as part of a nuanced system of neuronal regulation that contributes to a network of heritable features determining permethrin resistance.

Analysis of Anopheles messeae s.l. intron gene polymorphism associated with imidacloprid resistance

Due to their high solubility and stability, neonicotinoid insecticides are able to accumulate in water bodies, affecting aquatic organisms. The aims of this study were to evaluate resistance (LC₅₀ ) of Anopheles messeae s.l. (Anopheles messeae and An. daciae) to the neonicotinoid imidacloprid and to search for genetic markers associated with insecticide resistance. The LC₅₀ values of these species in the collections during 2017 and 2018 were indistinguishable and were in the range of 0.027-0.051 mg/l. In general, the LC₅₀ values of the mosquitoes were comparable with values of other mosquito species of the Anopheles and Culex genera. Gene polymorphisms of the variations in intron lengths and the presence of restriction sites in introns that were potentially associated with the metabolism of insecticides were studied. Polymorphisms of the studied genes in the pair of closely related species considered overlapped, but allele frequencies were different. Part of the genetic variants arose due to insertions of repetitive elements of the genome. Two variants of the cytochrome P450 gene Cyp6AG1 in An. daciae were associated with increased resistance to imidacloprid. Possible side effects of selection on insecticide resistance in blood-sucking mosquitoes are discussed.

Screening of insecticide resistance in Aedes aegypti populations collected from parishes in Eastern Jamaica

Owing to the increased reports in Aedes-borne diseases in the Caribbean and Latin America, the United States Agency for International Development assisted the Jamaican Ministry of Health and Wellness in conducting insecticide susceptibility tests on Aedes aegypti populations.

Sentinel sites were established in seven parishes of Jamaica (St. Catherine, Kingston and St. Andrew, St. Thomas, Portland, St. Mary and St. Ann) and Aedes aegypti eggs were collected, reared to adults per collected population and their susceptibility to varying pyrethroids and organophosphates were tested using the World Health Organization paper bioassays for these insecticides. The Centers for Disease Control and Prevention bottle bioassay was used to assess susceptibility to the carbamate, bendiocarb. The voltage gated sodium channel gene mutations V1016I and I1011V, normally associated with pyrethroid resistance, were also analysed.

The results showed that Aedes aegypti collected from all parishes exhibited resistance to pyrethroids at the following concentrations, permethrin 0.25–2.5%; deltamethrin 0.03–0.15%; lambda-cyhalothrin 0.03–0.3%; and etofenprox 0.5–2.5%. The insecticide deltamethrin at concentration 0.3% was the only pyrethroid tested that resulted in high mortality, 94.9 ± 0.34% knockdown within 1 hour of exposure and 98.95 ± 0.01% mortality (p <0.01) at 24 hours post exposure. The frequency of the voltage gated sodium channel gene mutation V1016I was high in the tested population, possibly accounting for the reduced sensitivity to pyrethroids. Organophosphate resistance was also observed in all populations tested. Mortality rates for 0.8% Malathion was 0.8 ± 0.70–60.68 ± 0.38% after 24 hour and 0.00–47.10 ± 3.02%, for pirimiphos-methyl 0.21%. Bendiocarb applied as 12.5 μg/ bottle resulted in mortality rates of 76.25 ± 4.30–100 ± 0.00% after 30 minutes of exposure.

The results showed that Ae. aegypti from the seven parishes analysed demonstrated resistance to the insecticides tested. Deltamethrin and bendiocarb at concentrations 0.3% and 12.5μg respectively, were considered most effective, causing high mortality in the local populations. Routine monitoring and evaluations of Ae. aegypti populations from the included parishes are recommended. Additionally, the study results represent the most comprehensive testing to date with local Aedes aegypti populations distributed across different parishes of Jamaica and should be useful to guide national and sub national strategies for vector control and surveillance.

The recent outbreaks of chikungunya, Zika and dengue viruses in the Caribbean and Latin America during 2013–2018, have motivated local health authorities in individual countries to seek international assistance from multiple agencies and adjust their vector management strategies in order to contain the transmission of Aedes-borne diseases. Countries within the region rely mainly on the use of commercially available insecticides to manage the mosquito populations and by extension, reduce the transmission of vector-borne diseases. However, this practice is stymied by the ability of these insects to develop resistance mechanisms to insecticides of frequent use. Many countries in the Caribbean and Latin America regions lack the facilities and supplies to conduct routine insecticide susceptibility testing. In an effort to ascertain knowledge on the presence of insecticide resistance in Aedes aegypti populations, and to determine which insecticide is more effective to reduce vector populations of Jamaica, we conducted an array of bioassays with three different types of insecticides (pyrethroid, organophosphate, carbamate) to use local evidence and inform public health authorities on vector control status and alternatives.

Impact of deltamethrin selection on kdr mutations and insecticide detoxifying enzymes in Aedes aegypti from Mexico

Background

Insecticide resistance is a serious problem for vector control programmes worldwide. Resistance is commonly attributed to mutations at the insecticide’s target site or increased activity of detoxification enzymes.

Methods

We determined the knockdown concentration (KC₅₀) and lethal concentration (LC₅₀) of deltamethrin in six natural populations of adult Aedes aegypti from southeastern Mexico. These populations were then selected over five generations using the LC₅₀ from the preceding generation that underwent selection, and the heritability of deltamethrin resistance was quantified. For each generation, we also determined the frequency of the kdr alleles L410, I1016 and C1534, and the levels of activity of three enzyme families (α- and β-esterases, mixed-function oxidases and glutathione S-transferases (GST)) associated with insecticide detoxification.

Results

There was an increase in KC₅₀ and LC₅₀ values in the subsequent generations of selection with deltamethrin (F_S5vs F_S0). According to the resistance ratios (RRs), we detected increases in LC₅₀ ranging from 1.5 to 5.6 times the values of the parental generation and in KC₅₀ ranging from 1.3–3.8 times the values of the parental generation. Triple homozygous mutant individuals (tri-locus, LL/II/CC) were present in the parental generations and increased in frequency after selection. The frequency of L410 increased from 1.18-fold to 2.63-fold after selection with deltamethrin (F_S5vs F_S0) in the populations analyzed; for I1016 an increase between 1.19-fold to 2.79-fold was observed, and C1534 was fixed in all populations after deltamethrin selection. Enzymatic activity varied significantly over the generations of selection. However, only α- esterase activity remained elevated in multiple populations after five generations of deltamethrin selection. We observed an increase in the mean activity levels of GSTs in two of the six populations analyzed.

Conclusions

The high levels of resistance and their association with high frequencies of kdr mutations (V410L, V1016I and F1534C) obtained through artificial selection, suggest an important role of these mutations in conferring resistance to deltamethrin. We highlight the need to implement strategies that involve the monitoring of kdr frequencies in insecticide resistance monitoring and management programmes.

Rodriguez-Ramirez A, Dzul F, Manrique-Saide P, Black WC. Loss of pyrethroid resistance in newly established laboratory colonies of Aedes aegypti

Background

Resistance to pyrethroid insecticides in Aedes aegypti has become widespread after almost two decades of the frequent use of these pesticides to reduce arbovirus transmission. Despite this resistance, pyrethroids continue to be used because they are relatively inexpensive and have low human toxicity. Resistance management has been proposed as a way to retain the use of pyrethroids in natural populations. A key component of resistance management is the assumption that negative fitness is associated with resistance alleles such that resistance alleles will decline in frequency when the insecticides are removed. At least three studies in Ae. aegypti have demonstrated a decrease in pyrethroid resistance once the insecticide has been removed.

Methods/Principal findings

The present study aims to evaluate variation in the loss of pyrethroid resistance among newly established laboratory populations of Ae. aegypti from Mexico. Eight field collections were maintained for up to eight generations, and we recorded changes in the frequencies of the mutations at the V1,016I locus and at the F1,534C locus in the voltage-gated sodium channel gene (VGSC). I1,016 and C1,534 confer resistance. We also examined resistance ratios (RR) with type 1 and 2 pyrethroids.

Conclusions/Significance

We demonstrate that, in general, the frequency of the Ae. aegypti pyrethroid-resistance alleles I1,016 and C1,534 decline when they are freed from pyrethroid pressure in the laboratory. However, the pattern of decline is strain dependent. In agreement with earlier studies, the RR was positively correlated with the frequencies of the resistance allele I1,016 and showed significant protection against permethrin, and deltamethrin, whereas F1,534C showed protection against permethrin but not against deltamethrin.

The author is interested in the evolution of genes that confer resistance to insecticides, especially when this evolution affects the binding of insecticides to their target sites. The Voltage-Gated Sodium Channel gene provides an excellent opportunity to understand how mutations at the target site(s) affect the evolution of resistance in many different pest insect species, including Aedes aegypti, which is the primary vector of the dengue, yellow fever, Zika and chikungunya viruses.

Using targeted next-generation sequencing to characterize genetic differences associated with insecticide resistance in Culex quinquefasciatus populations from the southern U.S

Resistance to insecticides can hamper the control of mosquitoes such as Culex quinquefasciatus, known to vector arboviruses such as West Nile virus and others. The strong selective pressure exerted on a mosquito population by the use of insecticides can result in heritable genetic changes associated with resistance. We sought to characterize genetic differences between insecticide resistant and susceptible Culex quinquefasciatus mosquitoes using targeted DNA sequencing. To that end, we developed a panel of 122 genes known or hypothesized to be involved in insecticide resistance, and used an Ion Torrent PGM sequencer to sequence 125 unrelated individuals from seven populations in the southern U.S. whose resistance phenotypes to permethrin and malathion were known from previous CDC bottle bioassay testing. Data analysis consisted of discovering SNPs (Single Nucleotide Polymorphism) and genes with evidence of copy number variants (CNVs) statistically associated with resistance. Ten of the seventeen genes found to be present in higher copy numbers were experimentally validated with real-time PCR. Of those, six, including the gene with the knock-down resistance (kdr) mutation, showed evidence of a ≥ 1.5 fold increase compared to control DNA. The SNP analysis revealed 228 unique SNPs that had significant p-values for both a Fisher’s Exact Test and the Cochran-Armitage Test for Trend. We calculated the population frequency for each of the 64 nonsynonymous SNPs in this group. Several genes not previously well characterized represent potential candidates for diagnostic assays when further validation is conducted.

Voltage-gated sodium channel intron polymorphism and four mutations comprise six haplotypes in an Aedes aegypti population in Taiwan

Background

Knockdown resistance (kdr) to dichlorodiphenyltrichloroethane (DDT) and pyrethroids is known to link amino acid substitutions in the voltage-gated sodium channel (VGSC) in Aedes aegypti. Dengue fever primarily transmitted by Ae. aegypti is an annual public health issue in Taiwan. Accordingly, pyrethroid insecticides have been heavily used for decades to control mosquito populations in the summer and autumn. In Taiwan, an Ae. aegypti population with two VGSC mutations, V1016G and D1763Y, was described previously.

Methodology/Principal finding

Aedes aegypti (G0) were collected in Tainan and Kaohsiung in southern Taiwan. The VGSC gene polymorphisms of the kdr mutations and the intron flanked by exons 20 and 21 were verified. The first generation offspring (G1) were used to measure the resistance level to cypermethrin, a pyrethroid insecticide currently used in Taiwan. In addition to V1016G and D1763Y, we describe two new mutations, S989P and F1534C, which have not been reported in Taiwan. Moreover, we also identify two types (groups A and B) of introns between exons 20 and 21. Intriguingly, the kdr mutations S989P, V1016G and D1763Y are strictly located on the haplotype harboring the group A intron, whereas F1534C links to the group B intron. When those data were taken together, we proposed the following six haplotypes for VGSC genes in Taiwan today: (i)S989-intron A-V1016-F1534-D1763, (ii)S989-intron A-V1016G-F1534-D1763, (iii)S989P-intron A-V1016G-F1534-D1763, (iv)S989-intron A-V1016G-F1534-D1763Y, (v)S989-intron B-V1016-F1534-D1763 and (vi)S989-intron B-V1016-F1534C-D1763. Triple heterozygous mutations of either S989P/V1016G/F1534C or V1016G/F1534C/D1763Y can be found in one single Ae. aegypti mosquito. The proportions of the VGSC mutations were relevant to cypermethrin resistance. Notably, the presence of S989P and V1016G in the population could be a helpful reference to predict the resistance level to cypermethrin. This is the first study to demonstrate the coexistence of four kdr mutations in a population of Ae. aegypti.

Conclusions/Significance

Four kdr mutations (S989P, V1016G, F1534C and D1763Y) and two intron forms (Group A and B) were commonly found in local Ae. aegypti populations in Taiwan.

VGSC mutations of Aedes aegypti threaten vector control programs and have been brought to attention in dengue endemic areas. Taiwan has suffered dengue outbreaks, which usually begin from an imported case. Pyrethroid insecticides were used to kill infectious females and adults in the surrounding environment of each suspected case. In Taiwan, V1016G and F1763Y mutations of VGSC have been described previously. Here, we further describe two additional amino acid substitutions (S989P, F1534C) and two forms of the intron between exon 20 and 21. We also propose six haplotypes for VGSC genes in Taiwan today. In conclusion, four kdr mutations (S989P, V1016G, F1534C and D1763Y) and two intron forms (Group A and B) are commonly found in local Ae. aegypti populations in Taiwan.

Long-Term Population Studies Uncover the Genome Structure and Genetic Basis of Xenobiotic and Host Plant Adaptation in the Herbivore Tetranychus urticae

Pesticide resistance arises rapidly in arthropod herbivores, as can host plant adaptation, and both are significant problems in agriculture. These traits have been challenging to study as both are often polygenic and many arthropods are genetically intractable. Here, we examined the genetic architecture of pesticide resistance and host plant adaptation in the two-spotted spider mite, Tetranychus urticae, a global agricultural pest. We show that the short generation time and high fecundity of T. urticae can be readily exploited in experimental evolution designs for high-resolution mapping of quantitative traits. As revealed by selection with spirodiclofen, an acetyl-CoA carboxylase inhibitor, in populations from a cross between a spirodiclofen-resistant and a spirodiclofen-susceptible strain, and which also differed in performance on tomato, we found that a limited number of loci could explain quantitative resistance to this compound. These were resolved to narrow genomic intervals, suggesting specific candidate genes, including acetyl-CoA carboxylase itself, clustered and copy variable cytochrome P450 genes, and NADPH cytochrome P450 reductase, which encodes a redox partner for cytochrome P450s. For performance on tomato, candidate genomic regions for response to selection were distinct from those responding to the synthetic compound and were consistent with a more polygenic architecture. In accomplishing this work, we exploited the continuous nature of allele frequency changes across experimental populations to resolve the existing fragmented T. urticae draft genome to pseudochromosomes. This improved assembly was indispensable for our analyses, as it will be for future research with this model herbivore that is exceptionally amenable to genetic studies.

Vgsc-interacting proteins are genetically associated with pyrethroid resistance in Aedes aegypti

Association mapping of factors that condition pyrethroid resistance in Aedes aegypti has consistently identified genes in multiple functional groups. Toward better understanding of the mechanisms involved, we examined high throughput sequencing data (HTS) from two Aedes aegypti aegypti collections from Merida, Yucatan, Mexico treated with either permethrin or deltamethrin. Exome capture enrichment for coding regions and the AaegL5 annotation were used to identify genes statistically associated with resistance. The frequencies of single nucleotide polymorphisms (SNPs) were compared between resistant and susceptible mosquito pools using a contingency χ² analysis. The -log₁₀(χ²p value) was calculated at each SNP site, with a weighted average determined from all sites in each gene. Genes with -log₁₀(χ²p value) ≥ 4.0 and present among all 3 treatment groups were subjected to gene set enrichment analysis (GSEA). We found that several functional groups were enriched compared to all coding genes. These categories were transport, signal transduction and metabolism, in order from highest to lowest statistical significance. Strikingly, 21 genes with demonstrated association to synaptic function were identified. In the high association group (n = 1,053 genes), several genes were identified that also genetically or physically interact with the voltage-gated sodium channel (VGSC). These genes were eg., CHARLATAN (CHL), a transcriptional regulator, several ankyrin-domain proteins, PUMILIO (PUM), a translational repressor, and NEDD4 (E3 ubiquitin-protein ligase). There were 13 genes that ranked among the top 10%: these included VGSC; CINGULIN, a predicted neuronal gap junction protein, and the aedine ortholog of NERVY (NVY), a transcriptional regulator. Silencing of CHL and NVY followed by standard permethrin bottle bioassays validated their association with permethrin resistance. Importantly, VGSC levels were also reduced about 50% in chl- or nvy-dsRNA treated mosquitoes. These results are consistent with the contribution of a variety of neuronal pathways to pyrethroid resistance in Ae. aegypti.

Life and Death at the Voltage-Sensitive Sodium Channel: Evolution in Response to Insecticide Use

The voltage-sensitive sodium channel (VSSC) is a critical component of the insect nervous system. Pyrethroids and DDT are insecticides that have been widely used, and they kill insects by perturbations of the VSSC. Decades of insecticide use selected for mutations in Vssc that give rise to resistance in almost all pest insects. However, the mutations responsible for the resistance are not always the same, and some unusual patterns have emerged. This review focuses on what pyrethroid/DDT selection has done, in terms of Vssc changes that have occurred, using four well-studied species as examples of the differences that have evolved. Information is provided about the mutations that occur, potential pathways by which alleles with multiple mutations arose, the relative fitness of the alleles, the levels of resistance conferred, and the geographic distribution of the mutations. The lessons learned and exciting new areas of research are discussed.

Restoration of pyrethroid susceptibility in a highly resistant Aedes aegypti population

Insecticide resistance has evolved in disease vectors worldwide, creating the urgent need to either develop new control methods or restore insecticide susceptibility to regain use of existing tools. Here we show that phenotypic susceptibility can be restored in a highly resistant field-derived strain of Aedes aegypti in only 10 generations through rearing them in the absence of insecticide.

Temporal and spatial trends in insecticide resistance in Anopheles arabiensis in Sudan: outcomes from an evaluation of implications of insecticide resistance for malaria vector control

Background

Long-lasting insecticidal nets (LLINs) (with pyrethroids) and indoor residual spraying (IRS) are the cornerstones of the Sudanese malaria control program. Insecticide resistance to the principal insecticides in LLINs and IRS is a major concern. This study was designed to monitor insecticide resistance in Anopheles arabiensis from 140 clusters in four malaria-endemic areas of Sudan from 2011 to 2014. All clusters received LLINs, while half (n = 70), distributed across the four regions, had additional IRS campaigns.

Methods

Anopheles gambiae (s.l.) mosquitoes were identified to species level using PCR techniques. Standard WHO insecticide susceptibility bioassays were carried out to detect resistance to deltamethrin (0.05%), DDT (4%) and bendiocarb (0.1%). TaqMan assays were performed on random samples of deltamethrin-resistant phenotyped and pyrethrum spray collected individuals to determine Vgsc-1014 knockdown resistance mutations.

Results

Anopheles arabiensis accounted for 99.9% of any anopheline species collected across all sites. Bioassay screening indicated that mosquitoes remained susceptible to bendiocarb but were resistance to deltamethrin and DDT in all areas. There were significant increases in deltamethrin resistance over the four years, with overall mean percent mortality to deltamethrin declining from 81.0% (95% CI: 77.6–84.3%) in 2011 to 47.7% (95% CI: 43.5–51.8%) in 2014. The rate of increase in phenotypic deltamethrin-resistance was significantly slower in the LLIN + IRS arm than in the LLIN-only arm (Odds ratio 1.34; 95% CI: 1.02–1.77). The frequency of Vgsc-1014F mutation varied spatiotemporally with highest frequencies in Galabat (range 0.375–0.616) and New Halfa (range 0.241–0.447). Deltamethrin phenotypic-resistance correlated with Vgsc-1014F frequency.

Conclusion

Combining LLIN and IRS, with different classes of insecticide, may delay pyrethroid resistance development, but the speed at which resistance develops may be area-specific. Continued monitoring is vital to ensure optimal management and control.

Electronic supplementary material

The online version of this article (10.1186/s13071-018-2732-9) contains supplementary material, which is available to authorized users.

Levels of insecticide resistance to deltamethrin, malathion, and temephos, and associated mechanisms in Aedes aegypti mosquitoes from the Guadeloupe and Saint Martin islands (French West Indies)

Background

In the Guadeloupe and Saint Martin islands, Aedes aegypti mosquitoes are the only recognized vectors of dengue, chikungunya, and Zika viruses. For around 40 years, malathion was used as a mosquito adulticide and temephos as a larvicide. Since the European Union banned the use of these two insecticide molecules in the first decade of the 21st century, deltamethrin and Bacillus thuringiensis var. israelensis are the remaining adulticide and larvicide, respectively, used in Guadeloupe. In order to improve the management of vector control activities in Guadeloupe and Saint Martin, we investigated Ae. aegypti resistance to and mechanisms associated with deltamethrin, malathion, and temephos.

Methods

Ae. aegypti mosquitoes were collected from six different localities of Guadeloupe and Saint Martin. Larvae were used for malathion and temephos bioassays, and adult mosquitoes for deltamethrin bioassays, following World Health Organization recommendations. Knockdown resistance (Kdr) genotyping for V1016I and F1534C mutations, and expression levels of eight enzymes involved in detoxification mechanisms were examined in comparison with the susceptible reference Bora Bora strain.

Results

Resistance ratios (RR₅₀) calculated for Ae. aegypti larvae showed high resistance levels to temephos (from 8.9 to 33.1-fold) and low resistance levels to malathion (from 1.7 to 4.4-fold). Adult females displayed moderate resistance levels to deltamethrin regarding the time necessary to affect 50% of individuals, varying from 8.0 to 28.1-fold. Molecular investigations on adult mosquitoes showed high resistant allele frequencies for V1016I and F1534C (from 85 to 96% and from 90 to 98%, respectively), as well as an overexpression of the glutathione S-transferase gene, GSTe2, the carboxylesterase CCEae3a, and the cytochrome genes 014614, CYP6BB2, CYP6M11, and CYP9J23.

Conclusions

Ae. aegypti populations from Guadeloupe and Saint Martin exhibit multiple resistance to organophosphates (temephos and malathion), and pyrethroids (deltamethrin). The mechanisms associated with these resistance patterns show strong frequencies of F1534C and V1016I Kdr mutations, and an over-expression of CCEae3a, GSTe2, and four cytochrome P450 genes (014614, CYP9J23, CYP6M11, CYP6BB2). These results will form the baseline for a deeper understanding of the insecticide resistance levels and associated mechanisms of Ae. aegypti populations and will be used to improve vector control strategies in Guadeloupe and Saint Martin.

Electronic supplementary material

The online version of this article (doi:10.1186/s40249-017-0254-x) contains supplementary material, which is available to authorized users.

Predicting Wolbachia invasion dynamics in Aedes aegypti populations using models of density-dependent demographic traits

BACKGROUND

Arbovirus transmission by the mosquito Aedes aegypti can be reduced by the introduction and establishment of the endosymbiotic bacteria Wolbachia in wild populations of the vector. Wolbachia spreads by increasing the fitness of its hosts relative to uninfected mosquitoes. However, mosquito fitness is also strongly affected by population size through density-dependent competition for limited food resources. We do not understand how this natural variation in fitness affects symbiont spread, which limits our ability to design successful control strategies.

RESULTS

We develop a mathematical model to predict A. aegypti-Wolbachia dynamics that incorporates larval density-dependent variation in important fitness components of infected and uninfected mosquitoes. Our model explains detailed features of the mosquito-Wolbachia dynamics observed in two independent experimental A. aegypti populations, allowing the combined effects on dynamics of multiple density-dependent fitness components to be characterized. We apply our model to investigate Wolbachia field release dynamics, and show how invasion outcomes can depend strongly on the severity of density-dependent competition at the release site. Specifically, the ratio of released relative to wild mosquitoes required to attain a target infection frequency (at the end of a release program) can vary by nearly an order of magnitude. The time taken for Wolbachia to become established following releases can differ by over 2 years. These effects depend on the relative fitness of field and insectary-reared mosquitoes.

CONCLUSIONS

Models of Wolbachia invasion incorporating density-dependent demographic variation in the host population explain observed dynamics in experimental A. aegypti populations. These models predict strong effects of density-dependence on Wolbachia dynamics in field populations, and can assist in the effective use of Wolbachia to control the transmission of arboviruses such as dengue, chikungunya and zika.

Pyrethroid resistance in Aedes aegypti and Aedes albopictus: Important mosquito vectors of human diseases

Aedes aegypti and A. albopictus mosquitoes are vectors of important human disease viruses, including dengue, yellow fever, chikungunya and Zika. Pyrethroid insecticides are widely used to control adult Aedes mosquitoes, especially during disease outbreaks. Herein, we review the status of pyrethroid resistance in A. aegypti and A. albopictus, mechanisms of resistance, fitness costs associated with resistance alleles and provide suggestions for future research. The widespread use of pyrethroids has given rise to many populations with varying levels of resistance worldwide, albeit with substantial geographical variation. In adult A. aegypti and A. albopictus, resistance levels are generally lower in Asia, Africa and the USA, and higher in Latin America, although there are exceptions. Susceptible populations still exist in several areas of the world, particularly in Asia and South America. Resistance to pyrethroids in larvae is also geographically widespread. The two major mechanisms of pyrethroid resistance are increased detoxification due to P450-monooxygenases, and mutations in the voltage sensitive sodium channel (Vssc) gene. Several P450s have been putatively associated with insecticide resistance, but the specific P450s involved are not fully elucidated. Pyrethroid resistance can be due to single mutations or combinations of mutations in Vssc. The presence of multiple Vssc mutations can lead to extremely high levels of resistance. Suggestions for future research needs are presented.

Reproductive Incompatibility Involving Senegalese Aedes aegypti (L) Is Associated with Chromosome Rearrangements

Aedes aegypti, the primary vector of dengue, yellow fever and Zika flaviviruses, consists of at least two subspecies. Aedes aegypti (Aaa) is light in color, has pale scales on the first abdominal tergite, oviposits in artificial containers, and preferentially feeds on humans. Aedes aegypti formosus (Aaf), has a dark cuticle, is restricted to sub-Saharan Africa, has no pale scales on the first abdominal tergite and frequently oviposits in natural containers. Scale patterns correlate with cuticle color in East Africa but not in Senegal, West Africa where black cuticle mosquitoes display a continuum of scaling patterns and breed domestically indoors. An earlier laboratory study did not indicate any pre- or postzygotic barriers to gene flow between Aaa and Aaf in East Africa. However, similar attempts to construct F₁ intercross families between Aaa laboratory strains and Senegal Ae. aegypti (SenAae) failed due to poor F₁ oviposition and low F₂ egg-to-adult survival. Insemination and assortative mating experiments failed to identify prezygotic mating barriers. Backcrosses were performed to test for postzygotic isolation patterns consistent with Haldane’s rule modified for species, like Aedes, that have an autosomal sex determining locus (SDL). Egg-pupal survival was predicted to be low in females mated to hybrid F₁ males but average when a male mates with a hybrid F₁ female. Survival was in fact significantly reduced when females mated to hybrid males but egg-pupal survival was significantly increased when males were mated to hybrid F₁ females. These observations are therefore inconclusive with regards to Haldane’s rule. Basic cytogenetic analyses and Fluorescent In Situ Hybridization (FISH) experiments were performed to compare SenAae strains with the IB12 strain of Aaa that was used for genome sequencing and physical mapping. Some SenAae strains had longer chromosomes than IB12 and significantly different centromeric indices on chromosomes 1 and 3. DAPI staining was used to identify AT-rich regions, chromomycin A3 following pretreatment with barium hydroxide stained for GC-rich regions and stained the ribosomal RNA locus and YOYO-1 was used to test for differential staining. Chromosome patterns in SenAae strains revealed by these three stains differed from those in IB12. For FISH, 40 BAC clones previously physically mapped on Aaa chromosomes were used to test for chromosome rearrangements in SenAae relative to IB12. Differences in the order of markers identified two chromosomal rearrangements between IB12 and SenAae strains. The first rearrangement involves two overlapping pericentric (containing the centromere) inversions in chromosome 3 or an insertion of a large fragment into the 3q arm. The second rearrangement is close to the centromere on the p arm of chromosome 2. Linkage analysis of the SDL and the white-eye locus identified a likely chromosomal rearrangement on chromosome 1. The reproductive incompatibility observed within SenAae and between SenAae and Aaa may be generally associated with chromosome rearrangements on all three chromosomes and specifically caused by pericentric inversions on chromosomes 2 and 3.

Aedes aegypti is one of the best studied mosquito species and it is the principal vector of dengue, Zika, and yellow fever flaviviruses and the Chikungunya alphavirus. Aedes aegypti occurs throughout all tropical and subtropical regions of the world, and previous population genetic studies have shown that the highest genetic diversity occurs in Africa. Aedes aegypti from Senegal, West Africa (SenAae) have a low oviposition rate; those that do oviposit have a low fecundity and poor egg-to-adult survival. Furthermore rearrangements were detected on all three chromosomes in SenAae. These observations are consistent with the presence of at least two cryptic subspecies of Ae. aegypti in Senegal arising from reproductive isolation due to chromosome rearrangements. Genetic control strategies are being considered for the suppression of Ae. aegypti populations worldwide. Barriers to gene flow in African Ae. aegypti populations could compromise these future control efforts.

Diagnostic doses and times for Phlebotomus papatasi and Lutzomyia longipalpis sand flies (Diptera: Psychodidae: Phlebotominae) using the CDC bottle bioassay to assess insecticide resistance

BACKGROUND

Insecticide resistance to synthetic chemical insecticides is a worldwide concern in phlebotomine sand flies (Diptera: Psychodidae), the vectors of Leishmania spp. parasites. The CDC bottle bioassay assesses resistance by testing populations against verified diagnostic doses and diagnostic times for an insecticide, but the assay has been used limitedly with sand flies. The objective of this study was to determine diagnostic doses and diagnostic times for laboratory Lutzomyia longipalpis (Lutz & Nieva) and Phlebotomus papatasi (Scopoli) to ten insecticides, including pyrethroids, organophosphates, carbamates, and DDT, that are used worldwide to control vectors.

METHODS

Bioassays were conducted in 1,000-ml glass bottles each containing 10-25 sand flies from laboratory colonies of L. longipalpis or P. papatasi. Four pyrethroids, three organophosphates, two carbamates and one organochlorine, were evaluated. A series of concentrations were tested for each insecticide, and four replicates were conducted for each concentration. Diagnostic doses were determined only during the exposure bioassay for the organophosphates and carbamates. For the pyrethroids and DDT, diagnostic doses were determined for both the exposure bioassay and after a 24-hour recovery period.

RESULTS

Both species are highly susceptible to the carbamates as their diagnostic doses are under 7.0 μg/ml. Both species are also highly susceptible to DDT during the exposure assay as their diagnostic doses are 7.5 μg/ml, yet their diagnostic doses for the 24-h recovery period are 650.0 μg/ml for Lu. longipalpis and 470.0 μg/ml for P. papatasi.

CONCLUSIONS

Diagnostic doses and diagnostic times can now be incorporated into vector management programs that use the CDC bottle bioassay to assess insecticide resistance in field populations of Lu. longipalpis and P. papatasi. These findings provide initial starting points for determining diagnostic doses and diagnostic times for other sand fly vector species and wild populations using the CDC bottle bioassay.

Quantitative trait locus mapping and functional genomics of an organophosphate resistance trait in the western corn rootworm, Diabrotica virgifera virgifera

The western corn rootworm, Diabrotica virgifera virgifera, is an insect pest of corn and population suppression with chemical insecticides is an important management tool. Traits conferring organophosphate insecticide resistance have increased in frequency amongst D. v. virgifera populations, resulting in the reduced efficacy in many corn-growing regions of the USA. We used comparative functional genomic and quantitative trait locus (QTL) mapping approaches to investigate the genetic basis of D. v. virgifera resistance to the organophosphate methyl-parathion. RNA from adult methyl-parathion resistant and susceptible adults was hybridized to 8331 microarray probes. The results predicted that 11 transcripts were significantly up-regulated in resistant phenotypes, with the most significant (fold increases ≥ 2.43) being an α-esterase-like transcript. Differential expression was validated only for the α-esterase (ST020027A20C03), with 11- to 13-fold greater expression in methyl-parathion resistant adults (P < 0.05). Progeny with a segregating methyl-parathion resistance trait were obtained from a reciprocal backcross design. QTL analyses of high-throughput single nucleotide polymorphism genotype data predicted involvement of a single genome interval. These data suggest that a specific carboyxesterase may function in field-evolved corn rootworm resistance to organophosphates, even though direct linkage between the QTL and this locus could not be established.

Rhodnius prolixus supergene families of enzymes potentially associated with insecticide resistance

Chagas disease or American trypanosomiasis, is a potentially life-threatening illness caused by the protozoan parasite, Trypanosoma cruzi. Once known as an endemic health problem of poor rural populations in Latin American countries, it has now spread worldwide. The parasite is transmitted by triatomine bugs, of which Rhodnius prolixus (Hemiptera, Reduviidae, Triatominae) is one of the vectors and a model organism. This species occurs mainly in Central and South American countries where the disease is endemic. Disease prevention focuses on vector control programs that, in general, rely intensely on insecticide use. However, the massive use of chemical insecticides can lead to resistance. One of the major mechanisms is known as metabolic resistance that is associated with an increase in the expression or activity of detoxification genes. Three of the enzyme families that are involved in this process - carboxylesterases (CCE), glutathione s-transferases (GST) and cytochrome P450s (CYP) - are analyzed in the R. prolixus genome. A similar set of detoxification genes to those of the Hemipteran Acyrthosiphon pisum but smaller than in most dipteran species was found in R. prolixus genome. All major CCE classes (43 genes found) are present but the pheromone/hormone processing class had fewer genes than usual. One main expansion was detected on the detoxification/dietary class. The phosphotriesterase family, recently associated with insecticide resistance, was also represented with one gene. One microsomal GST gene was found and the cytosolic GST gene count (14 genes) is extremely low when compared to the other hemipteran species with sequenced genomes. However, this is similar to Apis mellifera, a species known for its deficit in detoxification genes. In R. prolixus 88 CYP genes were found, with representatives in the four clans (CYP2, CYP3, CYP4 and mitochondrial) usually found in insects. R. prolixus seems to have smaller species-specific expansions of CYP genes than mosquitoes and beetles, among others. The number of R. prolixus CYP genes is similar to the hemipteran Ac. pisum, although with a bigger expansion in CYP3 and CYP4 clans, along with several gene fragments, mostly in CYP4 clan. Eleven founding members of new families were detected, consisting of ten genes in the CYP3 clan and 1 gene in the CYP4 clan. Members of these clans were proposed to have important detoxification roles in insects. The identification of CCE, GST and CYP genes is of utmost importance for directing detoxification studies on triatomines that can help insecticide management strategies in control programs.

"The Environment is Everything That Isn't Me": Molecular Mechanisms and Evolutionary Dynamics of Insect Clocks in Variable Surroundings

Circadian rhythms are oscillations in behavior, metabolism and physiology that have a period close to 24 h. These rhythms are controlled by an internal pacemaker that evolved under strong selective pressures imposed by environmental cyclical changes, mainly of light and temperature. The molecular nature of the circadian pacemaker was extensively studied in a number of organisms under controlled laboratory conditions. But although these studies were fundamental to our understanding of the circadian clock, most of the environmental conditions used resembled rather crudely the relatively constant situation at lower latitudes. At higher latitudes light-dark and temperature cycles vary considerably across different seasons, with summers having long and hot days and winters short and cold ones. Considering these differences and other external cues, such as moonlight, recent studies in more natural and semi-natural situations revealed unexpected features at both molecular and behavioral levels, highlighting the dramatic influence of multiple environmental variables in the molecular clockwork. This emphasizes the importance of studying the circadian clock in the wild, where seasonal environmental changes fine-tune the underlying circadian mechanism, affecting population dynamics and impacting the geographical variation in clock genes. Indeed, latitudinal clines in clock gene frequencies suggest that natural selection and demography shape the circadian clock over wide geographical ranges. In this review we will discuss the recent advances in understanding the molecular underpinnings of the circadian clock, how it resonates with the surrounding variables (both in the laboratory and in semi-natural conditions) and its impact on population dynamics and evolution. In addition, we will elaborate on how next-generation sequencing technologies will complement classical reductionist approaches by identifying causal variants in natural populations that will link genetic variation to circadian phenotypes, illuminating how the circadian clock functions in the real world.

Coevolution of the Ile1,016 and Cys1,534 Mutations in the Voltage Gated Sodium Channel Gene of Aedes aegypti in Mexico

Background

Worldwide the mosquito Aedes aegypti (L.) is the principal urban vector of dengue viruses. Currently 2.5 billion people are at risk for infection and reduction of Ae. aegypti populations is the most effective means to reduce the risk of transmission. Pyrethroids are used extensively for adult mosquito control, especially during dengue outbreaks. Pyrethroids promote activation and prolong the activation of the voltage gated sodium channel protein (VGSC) by interacting with two distinct pyrethroid receptor sites [1], formed by the interfaces of the transmembrane helix subunit 6 (S6) of domains II and III. Mutations of S6 in domains II and III synergize so that double mutants have higher pyrethroid resistance than mutants in either domain alone. Computer models predict an allosteric interaction between mutations in the two domains. In Ae. aegypti, a Ile1,016 mutation in the S6 of domain II was discovered in 2006 and found to be associated with pyrethroid resistance in field populations in Mexico. In 2010 a second mutation, Cys1,534 in the S6 of domain III was discovered and also found to be associated with pyrethroid resistance and correlated with the frequency of Ile1,016.

Methodology/Principal Findings

A linkage disequilibrium analysis was performed on Ile1,016 and Cys1,534 in Ae. aegypti collected in Mexico from 2000–2012 to test for statistical associations between S6 in domains II and III in natural populations. We estimated the frequency of the four dilocus haplotypes in 1,016 and 1,534: Val1,016/Phe1,534 (susceptible), Val1,016/Cys1,534, Ile1,016/Phe1,534, and Ile1,016/Cys1,534 (resistant). The susceptible Val1,016/Phe1,534 haplotype went from near fixation to extinction and the resistant Ile1,016/Cys1,534 haplotype increased in all collections from a frequency close to zero to frequencies ranging from 0.5–0.9. The Val1,016/Cys1,534 haplotype increased in all collections until 2008 after which it began to decline as Ile1,016/Cys1,534 increased. However, the Ile1,016/Phe1,534 haplotype was rarely detected; it reached a frequency of only 0.09 in one collection and subsequently declined.

Conclusion/Significance

Pyrethroid resistance in the vgsc gene requires the sequential evolution of two mutations. The Ile1,016/Phe1,534 haplotype appears to have low fitness suggesting that Ile1,016 was unlikely to have evolved independently. Instead the Cys1,534 mutation evolved first but conferred only a low level of resistance. Ile1,016 in S6 of domain II then arose from the Val1,016/Cys1,534 haplotype and was rapidly selected because double mutants confer higher pyrethroid resistance. This pattern suggests that knowledge of the frequencies of mutations in both S6 in domains II and III are important to predict the potential of a population to evolve kdr. Susceptible populations with high Val1,016/Cys1,534 frequencies are at high risk for kdr evolution, whereas susceptible populations without either mutation are less likely to evolve high levels of kdr, at least over a 10 year period.

Constant use of pyrethroid insecticides has driven mosquito populations to develop resistance. In Aedes aegypti, the primary mosquito vector of dengue, yellow Fever, and chikungunya viruses, pyrethroid resistance is primarily associated with mutations in the voltage-gated sodium channel protein. One mutation occurs in codon 1,016 and involves a replacement of valine with isoleucine (Ile1, 016), and a second located in subunit 6 of domain III in codon 1,534, replaces phenylalanine with cysteine (Cys1,534). In Mexico, we found that Cys1,534 was present in the same mosquito collections that were previously analyzed for Ile1,016. In this study, we performed a linkage disequilibrium analysis on both Ile1,016 and Cys1,534 in Mexican collections from 2000–2012. Our analysis suggests that pyrethroid resistance requires the sequential evolution of the two mutations and that Cys1,534 must occur first and appears to enable the Ile1,016 mutation to survive.

Identification of a major Quantitative Trait Locus determining resistance to the organophosphate temephos in the dengue vector mosquito Aedes aegypti

Organophosphate insecticides (OP) have extensively been used to control mosquitoes, such as the vector Aedes aegypti. Unfortunately, OP resistance has hampered control programs worldwide. We used Quantitative Trait Locus (QTL) mapping to evaluate temephos resistance in two F1 intercross populations derived from crosses between a resistant Ae. aegypti strain (RecR) and two susceptible strains (MoyoD and Red). A single major effect QTL was identified on chromosome 2 of both segregating populations, named rtt1 (resistance to temephos 1). Bioinformatics analyses identified a cluster of carboxylesterase genes (CCE) within the rtt1 interval. qRT-PCR demonstrated that different CCEs were up-regulated in F2 resistant individuals from both crosses. However, none exceeded the 2-fold expression. Primary mechanisms for temephos resistance may vary between Ae. aegypti populations, yet also appear to support previous findings suggesting that multiple linked esterase genes may contribute to temephos resistance in the RecR strain as well as other populations.

Identification of QTLs Conferring Resistance to Deltamethrin in Culex pipiens pallens

Culex pipiens pallens is the most abundant Culex mosquito species in northern China and is an important vector of bancroftian filariasis and West Nile virus. Deltamethrin is an insecticide that is widely used for mosquito control, however resistance to this and other insecticides has become a major challenge in the control of vector-borne diseases that appear to be inherited quantitatively. Furthermore, the genetic basis of insecticide resistance remains poorly understood. In this study, quantitative trait loci (QTL) mapping of resistance to deltamethrin was conducted in F₂ intercross segregation populations using bulked segregation analysis (BSA) and amplified fragment length polymorphism markers (AFLP) in Culex pipiens pallens. A genetic linkage map covering 381 cM was constructed and a total of seven QTL responsible for resistance to deltamethrin were detected by composite interval mapping (CIM), which explained 95% of the phenotypic variance. The major QTL in linkage group 2 accounted for 62% of the variance and is worthy of further study. 12 AFLP markers in the map were cloned and the genomic locations of these marker sequences were determined by applying the Basic Local Alignment Search Tool (BLAST) tool to the genome sequence of the closely related Culex quinquefasciatus. Our results suggest that resistance to deltamethrin is a quantitative trait under the control of a major QTL in Culex pipiens pallens. Cloning of related AFLP markers confirm the potential utility for anchoring the genetic map to the physical map. The results provide insight into the genetic architecture of the trait.

Assessing Insecticide Susceptibility of Laboratory Lutzomyia longipalpis and Phlebotomus papatasi Sand Flies (Diptera: Psychodidae: Phlebotominae)

Chemical insecticides are effective for controlling Lutzomyia and Phlebotomus sand fly (Diptera: Psychodidae) vectors of Leishmania parasites. However, repeated use of certain insecticides has led to tolerance and resistance. The objective of this study was to determine lethal concentrations (LCs) and lethal exposure times (LTs) to assess levels of susceptibility of laboratory Lutzomyia longipalpis (Lutz and Nieva) and Phlebotomus papatasi (Scopoli) to 10 insecticides using a modified version of the World Health Organization (WHO) exposure kit assay and Centers for Disease Control and Prevention (CDC) bottle bioassay. Sand flies were exposed to insecticides coated on the interior of 0.5-gallon and 1,000-ml glass bottles. Following exposure, the flies were allowed to recover for 24 h, after which mortality was recorded. From dose-response survival curves for L. longipalpis and P. papatasi generated with the QCal software, LCs causing 50, 90, and 95% mortality were determined for each insecticide. The LCs and LTs from this study will be useful as baseline reference points for future studies using the CDC bottle bioassays to assess insecticide susceptibility of sand fly populations in the field. There is a need for a larger repository of sand fly insecticide susceptibility data from the CDC bottle bioassays, including a range of LCs and LTs for more sand fly species with more insecticides. Such a repository would be a valuable tool for vector management.

The genetic architecture of a complex trait: Resistance to multiple toxins produced by Bacillus thuringiensis israelensis in the dengue and yellow fever vector, the mosquito Aedes aegypti

The bacterial insecticide Bacillus thuringiensis subsp. israelensis (Bti) is an increasingly popular alternative to chemical insecticides for controlling mosquito populations. Because Bti toxicity relies on the action of four main toxins, resistance to Bti is very likely a complex phenotype involving several genes simultaneously. Dissecting the underlying genetic basis thus requires associating a quantitative measure of resistance to genetic variation at many loci in a segregating population. Here, we undertake this task using the dengue and yellow fever vector, the mosquito Aedes aegypti, as a study model. We conducted QTL (Quantitative Trait Locus) and admixture mapping analyses on two controlled crosses and on an artificial admixed population, respectively, all obtained from resistant and susceptible lab strains. We detected 16 QTL regions, among which four QTLs were revealed by different analysis methods. These four robust QTLs explained altogether 29.2% and 62.2% of the total phenotypic variance in the two QTL crosses, respectively. They also all showed a dominant mode of action. In addition, we found six loci showing statistical association with Bti resistance in the admixed population. Five of the supercontigs highlighted in this study contained candidate genes as suggested by their function, or by prior evidence from expression and/or outlier analyses. These genomic regions are thus good starting points for fine mapping of resistance to Bti or functional analyses aiming at identifying the underlying genes and mutations. Moreover, for the purpose of this work, we built the first Ae. aegypti genetic map based on markers associated with genes expressed in larvae. This genetic map harbors 229 SNP markers mapped across the three chromosomes for a total length of 311.9cM. It brought to light several assembly discrepancies with the reference genome, suggesting a high level of genome plasticity in Ae. aegypti.

Selective sweep analysis in the genomes of the 91-R and 91-C Drosophila melanogaster strains reveals few of the 'usual suspects' in dichlorodiphenyltrichloroethane (DDT) resistance

Adaptation of insect phenotypes for survival after exposure to xenobiotics can result from selection at multiple loci with additive genetic effects. To the authors' knowledge, no selective sweep analysis has been performed to identify such loci in highly dichlorodiphenyltrichloroethane (DDT) resistant insects. Here we compared a highly DDT resistant phenotype in the Drosophila melanogaster (Drosophila) 91-R strain to the DDT susceptible 91-C strain, both of common origin. Whole genome re-sequencing data from pools of individuals was generated separately for 91-R and 91-C, and mapped to the reference Drosophila genome assembly (v. 5.72). Thirteen major and three minor effect chromosome intervals with reduced nucleotide diversity (π) were identified only in the 91-R population. Estimates of Tajima's D (D) showed corresponding evidence of directional selection in these same genome regions of 91-R, however, no similar reductions in π or D estimates were detected in 91-C. An overabundance of non-synonymous proteins coding to synonymous changes were identified in putative open reading frames associated with 91-R. Except for NinaC and Cyp4g1, none of the identified genes were the 'usual suspects' previously observed to be associated with DDT resistance. Additionally, up-regulated ATP-binding cassette transporters have been previously associated with DDT resistance; however, here we identified a structurally altered MDR49 candidate resistance gene. The remaining fourteen genes have not previously been shown to be associated with DDT resistance. These results suggest hitherto unknown mechanisms of DDT resistance, most of which have been overlooked in previous transcriptional studies, with some genes having orthologs in mammals.

Local evolution of pyrethroid resistance offsets gene flow among Aedes aegypti collections in Yucatan State, Mexico

The mosquito Aedes aegypti is the major vector of the four serotypes of dengue virus (DENV1-4). Previous studies have shown that Ae. aegypti in Mexico have a high effective migration rate and that gene flow occurs among populations that are up to 150 km apart. Since 2000, pyrethroids have been widely used for suppression of Ae. aegypti in cities in Mexico. In Yucatan State in particular, pyrethroids have been applied in and around dengue case households creating an opportunity for local selection and evolution of resistance. Herein, we test for evidence of local adaptation by comparing patterns of variation among 27 Ae. aegypti collections at 13 single nucleotide polymorphisms (SNPs): two in the voltage-gated sodium channel gene para known to confer knockdown resistance, three in detoxification genes previously associated with pyrethroid resistance, and eight in putatively neutral loci. The SNPs in para varied greatly in frequency among collections, whereas SNPs at the remaining 11 loci showed little variation supporting previous evidence for extensive local gene flow. Among Ae. aegypti in Yucatan State, Mexico, local adaptation to pyrethroids appears to offset the homogenizing effects of gene flow.

QTL mapping of genome regions controlling temephos resistance in larvae of the mosquito Aedes aegypti

INTRODUCTION

The mosquito Aedes aegypti is the principal vector of dengue and yellow fever flaviviruses. Temephos is an organophosphate insecticide used globally to suppress Ae. aegypti larval populations but resistance has evolved in many locations.

METHODOLOGY/PRINCIPAL FINDINGS

Quantitative Trait Loci (QTL) controlling temephos survival in Ae. aegypti larvae were mapped in a pair of F3 advanced intercross lines arising from temephos resistant parents from Solidaridad, México and temephos susceptible parents from Iquitos, Peru. Two sets of 200 F3 larvae were exposed to a discriminating dose of temephos and then dead larvae were collected and preserved for DNA isolation every two hours up to 16 hours. Larvae surviving longer than 16 hours were considered resistant. For QTL mapping, single nucleotide polymorphisms (SNPs) were identified at 23 single copy genes and 26 microsatellite loci of known physical positions in the Ae. aegypti genome. In both reciprocal crosses, Multiple Interval Mapping identified eleven QTL associated with time until death. In the Solidaridad×Iquitos (SLD×Iq) cross twelve were associated with survival but in the reciprocal IqxSLD cross, only six QTL were survival associated. Polymorphisms at acetylcholine esterase (AchE) loci 1 and 2 were not associated with either resistance phenotype suggesting that target site insensitivity is not an organophosphate resistance mechanism in this region of México.

CONCLUSIONS/SIGNIFICANCE

Temephos resistance is under the control of many metabolic genes of small effect and dispersed throughout the Ae. aegypti genome.

A single crossing-over event in voltage-sensitive Na+ channel genes may cause critical failure of dengue mosquito control by insecticides

The voltage-sensitive sodium (Na⁺) channel (Vssc) is the target site of pyrethroid insecticides. Pest insects develop resistance to this class of insecticide by acquisition of one or multiple amino acid substitution(s) in this channel. In Southeast Asia, two major Vssc types confer pyrethroid resistance in the dengue mosquito vector Aedes aegypti, namely, S989P+V1016G and F1534C. We expressed several types of Vssc in Xenopus oocytes and examined the effect of amino acid substitutions in Vssc on pyrethroid susceptibilities. S989P+V1016G and F1534C haplotypes reduced the channel sensitivity to permethrin by 100- and 25-fold, respectively, while S989P+V1016G+F1534C triple mutations reduced the channel sensitivity to permethrin by 1100-fold. S989P+V1016G and F1534C haplotypes reduced the channel sensitivity to deltamethrin by 10- and 1-fold (no reduction), respectively, but S989P+V1016G+F1534C triple mutations reduced the channel sensitivity to deltamethrin by 90-fold. These results imply that pyrethroid insecticides are highly likely to lose their effectiveness against A. aegypti if such a Vssc haplotype emerges as the result of a single crossing-over event; thus, this may cause failure to control this key mosquito vector. Here, we strongly emphasize the importance of monitoring the occurrence of triple mutations in Vssc in the field population of A. aegypti.

Pyrethroids are one of the major classes of insecticides that is widely used to control mosquito vectors. The target site of pyrethroids is found in the voltage-sensitive Na⁺ channel (Vssc) consisting of about 2100 amino acid residues. In this study we generated several types of Vssc with a single or multiple mutations, expressed in Xenopus oocytes, and examined their electrophysiological properties using two-electrode voltage clamp method. We confirmed that Aedes aegypti Vssc harboring a triple mutations exhibited extremely high levels of resistance to pyrethroid insecticides. This Vssc type can be generated by a single crossing-over event of two resistant Vssc genes that are widely distributed in Southeast Asia, one of the greatest dengue endemic areas. Our results highlight the importance of intensive monitoring for the triple mutations in Vssc in the A. aegypti mosquito.

Fine-mapping QTLs in advanced intercross lines and other outbred populations

Quantitative genetic studies in model organisms, particularly in mice, have been extremely successful in identifying chromosomal regions that are associated with a wide variety of behavioral and other traits. However, it is now widely understood that identification of the underlying genes will be far more challenging. In the last few years, a variety of populations have been utilized in an effort to more finely map these chromosomal regions with the goal of identifying specific genes. The common property of these newer populations is that linkage disequilibrium spans relatively short distances, which permits fine-scale mapping resolution. This review focuses on advanced intercross lines (AILs) which are the simplest such population. As originally proposed in 1995 by Darvasi and Soller, an AIL is the product of intercrossing two inbred strains beyond the F2 generation. Unlike recombinant inbred strains, AILs are maintained as outbred populations; brother-sister matings are specifically avoided. Each generation of intercrossing beyond the F2 further degrades linkage disequilibrium between adjacent makers, which allows for fine-scale mapping of quantitative trait loci (QTLs). Advances in genotyping technology and techniques for the statistical analysis of AILs have permitted rapid advances in the application of AILs. We review some of the analytical issues and available software, including QTLRel, EMMA, EMMAX, GEMMA, TASSEL, GRAMMAR, WOMBAT, Mendel, and others.

Mechanisms of pyrethroid resistance in the dengue mosquito vector, Aedes aegypti: target site insensitivity, penetration, and metabolism

Aedes aegypti is the major vector of yellow and dengue fevers. After 10 generations of adult selection, an A. aegypti strain (SP) developed 1650-fold resistance to permethrin, which is one of the most widely used pyrethroid insecticides for mosquito control. SP larvae also developed 8790-fold resistance following selection of the adults. Prior to the selections, the frequencies of V1016G and F1534C mutations in domains II and III, respectively, of voltage-sensitive sodium channel (Vssc, the target site of pyrethroid insecticide) were 0.44 and 0.56, respectively. In contrast, only G1016 alleles were present after two permethrin selections, indicating that G1016 can more contribute to the insensitivity of Vssc than C1534. In vivo metabolism studies showed that the SP strain excreted permethrin metabolites more rapidly than a susceptible SMK strain. Pretreatment with piperonyl butoxide caused strong inhibition of excretion of permethrin metabolites, suggesting that cytochrome P450 monooxygenases (P450s) play an important role in resistance development. In vitro metabolism studies also indicated an association of P450s with resistance. Microarray analysis showed that multiple P450 genes were over expressed during the larval and adult stages in the SP strain. Following quantitative real time PCR, we focused on two P450 isoforms, CYP9M6 and CYP6BB2. Transcription levels of these P450s were well correlated with the rate of permethrin excretion and they were certainly capable of detoxifying permethrin to 4'-HO-permethrin. Over expression of CYP9M6 was partially due to gene amplification. There was no significant difference in the rate of permethrin reduction from cuticle between SP and SMK strains.

Susceptibility to chlorpyrifos in pyrethroid-resistant populations of Aedes aegypti (Diptera: Culicidae) from Mexico

Resistance to the organophosphate insecticide chlorpyrifos was evaluated in females from six strains of Aedes aegypti (L.) that expressed high levels of cross-resistance to eight pyrethroid insecticides. Relative to LC50 and LC90 at 24 h of a susceptible New Orleans (NO) strain, three strains were highly resistant to chlorpyrifos (Coatzacoalcos, resistance ratio [RRLC90 = 11.97; Pozarica, RRLC90 = 12.98; and Cosoleacaque, RRLC50 = 13.94 and RRLC90 = 17.57), one strain was moderately resistant (Veracruz, RRLC90 = 5.92), and two strains were susceptible (Tantoyuca and Martinez de la Torre, RRLC50 and RRLC90 < 5) in bottle bioassays according to Centers for Disease Control and Prevention. Furthermore, high levels of alpha- or beta-esterase activity in the sample populations were correlated with resistance, suggesting that esterase activity may be a mechanism causing the development of organophosphate resistance in these populations. Overall, the populations in this study were less resistant to chlorpyrifos than to pyrethroids. Rotation of insecticides used in control activities is recommended to delay or minimize the occurrence of high levels of resistance to chlorpyrifos among local populations of Ae. aegypti. The diagnostic dose and diagnostic time for chlorpyrifos resistance monitoring was determined to be 85 microg per bottle and 30 min, respectively, using the susceptible NO strain.

Differential transcription profiles in Aedes aegypti detoxification genes after temephos selection

The mosquito Aedes aegypti is the main vector of Dengue and Yellow Fever flaviviruses. The organophosphate insecticide temephos is a larvicide that is used globally to control Ae. aegypti populations; many of which have in turn evolved resistance. Target site alteration in the acetylcholine esterase of this species has not being identified. Instead, we tracked changes in transcription of metabolic detoxification genes using the Ae. aegypti 'Detox Chip' microarray during five generations of temephos selection. We selected for temephos resistance in three replicates in each of six collections, five from Mexico, and one from Peru. The response to selection was tracked in terms of lethal concentrations. Uniform upregulation was seen in the epsilon class glutathione-S-transferase (eGST) genes in strains from Mexico prior to laboratory selection, while eGSTs in the Iquitos Peru strain became upregulated after five generations of temephos selection. While expression of many carboxyl/cholinesterase esterase (CCE) genes increased with selection, no single esterase was consistently upregulated and this same pattern was noted in the cytochrome P450 monooxygenase (CYP) genes and in other genes involved in reduction or oxidation of xenobiotics. Bioassays using glutathione-S-transferase (GST), CCE and CYP inhibitors suggest that various CCEs instead of GSTs are the main metabolic mechanism conferring resistance to temephos. We show that temephos-selected strains show no cross resistance to permethrin and that genes associated with temephos selection are largely independent of those selected with permethrin in a previous study.

Proof of concept for a novel insecticide bioassay based on sugar feeding by adult Aedes aegypti (Stegomyia aegypti)

Aedes aegypti L. (Stegomyia aegypti) (Diptera: Culicidae) is the principal vector of dengue and yellow fever viruses in tropical and subtropical regions of the world. Disease management is largely based on mosquito control achieved by insecticides applied to interior resting surfaces and through space sprays. Population monitoring to detect insecticide resistance is a significant component of integrated disease management programmes. We developed a bioassay method for assessing insecticide susceptibility based on the feeding activity of mosquitoes on plant sugars. Our prototype sugar-insecticide feeding bioassay system was composed of inexpensive, disposable components, contained minimal volumes of insecticide, and was compact and highly transportable. Individual mosquitoes were assayed in a plastic cup that contained a sucrose-permethrin solution. Trypan blue dye was added to create a visual marker in the mosquito's abdomen for ingested sucrose-permethrin solution. Blue faecal spots provided further evidence of solution ingestion. With the sugar-insecticide feeding bioassay, the permethrin susceptibility of Ae. aegypti females from two field-collected strains was characterized by probit analysis of dosage-response data. The field strains were also tested by forced contact of females with permethrin residues on filter paper. Dosage-response patterns were similar, indicating that the sugar-insecticide feeding bioassay had appropriately characterized the permethrin susceptibility of the two strains.

How many genetic options for evolving insecticide resistance in heliothine and spodopteran pests?

The widely accepted paradigm for the development of insecticide resistance in field populations of insects is of selection for one or a very few genes of major effect. Limited genetic mapping data for organophosphate and pyrethroid resistance in heliothine and spodopteran pests generally agrees with this paradigm. However, other biochemical and transcriptomic data suggest a more complex set of changes in multiple P450 and esterase gene/enzyme systems in resistant strains of these species. We discuss possible explanations for this paradox, including the likely embedding of these genes in regulatory cascades and emerging evidence for their arrangement in large clusters of closely related genes. We conclude that there could indeed be an unusually large number of genetic options for evolving resistance in these species.

Wide spread cross resistance to pyrethroids in Aedes aegypti (Diptera: Culicidae) from Veracruz state Mexico

Seven F1 strains of Aedes aegypti (L.) were evaluated by bottle bioassay for resistance to the pyrethroids d-phenothrin, permethrin, deltamethrin, lambda-cyalothrin, bifenthrin, cypermethrin, alpha-cypermethrin, and z-cypermethrin. The New Orleans strain was used as a susceptible control. Mortality rates after a 1 h exposure and after a 24 h recovery period were determined. The resistance ratio between the 50% knockdown values (RR(KC50)) of the F1 and New Orleans strains indicated high levels of knockdown resistance. The RR(KC50) with alpha-cypermethrin varied from 10 to 100 among strains indicating high levels of knockdown resistance. Most of the strains had moderate resistance to d-phenothrin. Significant but much lower levels of resistance were detected for lambda-cyalothrin, permethrin, and cypermethrin. For zeta-cypermethrin and bifenthrin, only one strain exhibited resistance with RR(KC50) values of 10- and 21-fold, respectively. None of the strains showed RR(KC50) >10 with deltamethrin, and moderate resistance was seen in three strains, while the rest were susceptible. Mosquitoes from all strains exhibited some recovery from all pyrethroids except d-phenothrin. Regression analysis was used to analyze the relationship between RR(LC50) and RR(KC50). Both were highly correlated (R2 = 0.84-0.97) so that the slope could be used to determine how much additional pyrethroid was needed to ensure lethality. Slopes ranged from 0.875 for d-phenothrin (RR(LC50) approximately equal to RR(KC50)) to 8.67 for lambda-cyalothrin (-8.5-fold more insecticide needed to kill). Both RR(LC50) and RR(KC50) values were highly correlated for all pyrethroids except bifenthrin indicating strong cross-resistance. Bifenthrin appears to be an alternative pyrethroid without strong cross-resistance that could be used as an alternative to the current widespread use of permethrin in Mexico.

Impact of environment on mosquito response to pyrethroid insecticides: facts, evidences and prospects

By transmitting major human diseases such as malaria, dengue fever and filariasis, mosquito species represent a serious threat worldwide in terms of public health, and pose a significant economic burden for the African continent and developing tropical regions. Most vector control programmes aiming at controlling life-threatening mosquitoes rely on the use of chemical insecticides, mainly belonging to the pyrethroid class. However, resistance of mosquito populations to pyrethroids is increasing at a dramatic rate, threatening the efficacy of control programmes throughout insecticide-treated areas, where mosquito-borne diseases are still prevalent. In the absence of new insecticides and efficient alternative vector control methods, resistance management strategies are therefore critical, but these require a deep understanding of adaptive mechanisms underlying resistance. Although insecticide resistance mechanisms are intensively studied in mosquitoes, such adaptation is often considered as the unique result of the selection pressure caused by insecticides used for vector control. Indeed, additional environmental parameters, such as insecticides/pesticides usage in agriculture, the presence of anthropogenic or natural xenobiotics, and biotic interactions between vectors and other organisms, may affect both the overall mosquito responses to pyrethroids and the selection of resistance mechanisms. In this context, the present work aims at updating current knowledge on pyrethroid resistance mechanisms in mosquitoes and compiling available data, often from different research fields, on the impact of the environment on mosquito response to pyrethroids. Key environmental factors, such as the presence of urban or agricultural pollutants and biotic interactions between mosquitoes and their microbiome are discussed, and research perspectives to fill in knowledge gaps are suggested.

Identifying permethrin resistance loci in malaria vectors by genetic mapping

Identification of the major loci responsible for insecticide resistance in malaria vectors would aid the development and implementation of effective resistance management strategies, which are urgently needed to tackle the growing threat posed by resistance to the limited insecticides available for malaria control. Genome-wide association studies in the major malaria vector, Anopheles gambiae, have been hindered by the high degree of within-population structuring and very low levels of linkage disequilibrium hence we revisited the use of quantitative trait loci (QTL) mapping to study resistance phenotypes in this vector species. Earlier work, identified two major QTL associated with pyrethroid resistance in A. gambiae s.s. from East Africa using genetic crossing of laboratory-colonized resistant and susceptible strains. In this study, we report the results from genetic mapping of pyrethroid resistance in three isofemale pedigrees established from wild-caught female A. gambiae s.s. mosquitoes from Benin. We identified two QTL on chromosomes 2L and 3R in these field populations, in similar genomic locations to the QTL identified in laboratory strains. The relative merits of two alternative study designs are discussed and suggestions made for future genetic mapping studies of insecticide resistance in mosquitoes.

Dengue research opportunities in the Americas

Dengue is a systemic arthropod-borne viral disease of major global public health importance. At least 2.5 billion people who live in areas of the world where dengue occurs are at risk of developing dengue fever (DF) and its severe complications, dengue hemorrhagic fever (DHF) and dengue shock syndrome (DSS). Repeated reemergences of dengue in sudden explosive epidemics often cause public alarm and seriously stress healthcare systems. The control of dengue is further challenged by the lack of effective therapies, vaccines, and point-of-care diagnostics. Despite years of study, even its pathogenic mechanisms are poorly understood. This article discusses recent advances in dengue research and identifies challenging gaps in research on dengue clinical evaluation, diagnostics, epidemiology, immunology, therapeutics, vaccinology/clinical trials research, vector biology, and vector ecology. Although dengue is a major global tropical pathogen, epidemiologic and disease control considerations in this article emphasize dengue in the Americas.

Comparative transcriptome analyses of deltamethrin-resistant and -susceptible Anopheles gambiae mosquitoes from Kenya by RNA-Seq

Malaria causes more than 300 million clinical cases and 665,000 deaths each year, and the majority of the mortality and morbidity occurs in sub-Saharan Africa. Due to the lack of effective vaccines and wide-spread resistance to antimalarial drugs, mosquito control is the primary method of malaria prevention and control. Currently, malaria vector control relies on the use of insecticides, primarily pyrethroids. The extensive use of insecticides has imposed strong selection pressures for resistance in the mosquito populations. Consequently, resistance to pyrethroids in Anopheles gambiae, the main malaria vector in sub-Saharan Africa, has become a major obstacle for malaria control. A key element of resistance management is the identification of resistance mechanisms and subsequent development of reliable resistance monitoring tools. Field-derived An. gambiae from Western Kenya were phenotyped as deltamethrin-resistant or -susceptible by the standard WHO tube test, and their expression profile compared by RNA-seq. Based on the current annotation of the An. gambiae genome, a total of 1,093 transcripts were detected as significantly differentially accumulated between deltamethrin-resistant and -susceptible mosquitoes. These transcripts are distributed over the entire genome, with a large number mapping in QTLs previously linked to pyrethorid resistance, and correspond to heat-shock proteins, metabolic and transport functions, signal transduction activities, cytoskeleton and others. The detected differences in transcript accumulation levels between resistant and susceptible mosquitoes reflect transcripts directly or indirectly correlated with pyrethroid resistance. RNA-seq data also were used to perform a de-novo Cufflinks assembly of the An. gambiae genome.

The effect of ivermectin in seven strains of Aedes aegypti (Diptera: Culicidae) including a genetically diverse laboratory strain and three permethrin resistant strains

Seven different strains of Aedes aegypti (L.), including a genetically diverse laboratory strain, three laboratory-selected permethrin-resistant strains, a standard reference strain, and two recently colonized strains were fed on human blood containing various concentrations of ivermectin. Ivermectin reduced adult survival, fecundity, and hatch rate of eggs laid by ivermectin-treated adults in all seven strains. The LC50 of ivermectin for adults and the concentration that prevented 50% of eggs from hatching was calculated for all strains. Considerable variation in adult survival after an ivermectin-bloodmeal occurred among strains, and all three permethrin-resistant strains were significantly less susceptible to ivermectin than the standard reference strain. The hatch rate after an ivermectin bloodmeal was less variable among strains, and only one of the permethrin-resistant strains differed significantly from the standard reference strain. Our studies suggest that ivermectin induces adult mortality and decreases the hatch rate of eggs through different mechanisms. A correlation analysis of log-transformed LC50 among strains suggests that permethrin and ivermectin cross-resistance may occur.

Insecticide resistance in the dengue vector Aedes aegypti from Martinique: distribution, mechanisms and relations with environmental factors

Dengue is an important mosquito borne viral disease in Martinique Island (French West Indies). The viruses responsible for dengue are transmitted by Aedes aegypti, an indoor day-biting mosquito. The most effective proven method for disease prevention has been by vector control by various chemical or biological means. Unfortunately insecticide resistance has already been observed on the Island and recently showed to significantly reduce the efficacy of vector control interventions. In this study, we investigated the distribution of resistance and the underlying mechanisms in nine Ae. aegypti populations. Statistical multifactorial approach was used to investigate the correlations between insecticide resistance levels, associated mechanisms and environmental factors characterizing the mosquito populations. Bioassays revealed high levels of resistance to temephos and deltamethrin and susceptibility to Bti in the 9 populations tested. Biochemical assays showed elevated detoxification enzyme activities of monooxygenases, carboxylesterases and glutathione S-tranferases in most of the populations. Molecular screening for common insecticide target-site mutations, revealed the presence of the "knock-down resistance" V1016I Kdr mutation at high frequency (>87%). Real time quantitative RT-PCR showed the potential involvement of several candidate detoxification genes in insecticide resistance. Principal Component Analysis (PCA) performed with variables characterizing Ae. aegypti from Martinique permitted to underline potential links existing between resistance distribution and other variables such as agriculture practices, vector control interventions and urbanization. Insecticide resistance is widespread but not homogeneously distributed across Martinique. The influence of environmental and operational factors on the evolution of the resistance and mechanisms are discussed.