Skip to main content

Voltage-sensitive sodium channel (Vssc) mutations associated with pyrethroid insecticide resistance in Aedes aegypti (L.) from two districts of Jeddah, Kingdom of Saudi Arabia: baseline information for a Wolbachia release program

Abstract

Background

Dengue suppression often relies on control of the mosquito vector, Aedes aegypti, through applications of insecticides of which the pyrethroid group has played a dominant role. Insecticide resistance is prevalent in Ae. aegypti around the world, and the resulting reduction of insecticide efficacy is likely to exacerbate the impact of dengue. Dengue has been a public health problem in Saudi Arabia, particularly in Jeddah, since its discovery there in the 1990s, and insecticide use for vector control is widespread throughout the city. An alternative approach to insecticide use, based on blocking dengue transmission in mosquitoes by the endosymbiont Wolbachia, is being trialed in Jeddah following the success of this approach in Australia and Malaysia. Knowledge of insecticide resistance status of mosquito populations in Jeddah is a prerequisite for establishing a Wolbachia-based dengue control program as releases of Wolbachia mosquitoes succeed when resistance status of the release population is similar to that of the wild population.

Methods

WHO resistance bioassays of mosquitoes with deltamethrin, permethrin and DDT were used in conjunction with TaqMan® SNP Genotyping Assays to characterize mutation profiles of Ae. aegypti.

Results

Screening of the voltage-sensitive sodium channel (Vssc), the pyrethroid target site, revealed mutations at codons 989, 1016 and 1534 in Ae. aegypti from two districts of Jeddah. The triple mutant homozygote (1016G/1534C/989P) was confirmed from Al Safa and Al Rawabi. Bioassays with pyrethroids (Type I and II) and DDT showed that mosquitoes were resistant to each of these compounds based on WHO definitions. An association between Vssc mutations and resistance was established for the Type II pyrethroid, deltamethrin, with one genotype (989P/1016G/1534F) conferring a survival advantage over two others (989S/1016V/1534C and the triple heterozygote). An indication of synergism of Type I pyrethroid activity with piperonyl butoxide suggests that detoxification by cytochrome P450s accounts for some of the pyrethroid resistance response in Ae. aegypti populations from Jeddah.

Conclusions

The results provide a baseline for monitoring and management of resistance as well as knowledge of Vssc genotype frequencies required in Wolbachia release populations to ensure homogeneity with the target field population. Vssc mutation haplotypes observed show some similarity with those from Ae. aegypti in southeast Asia and the Indo-Pacific, but the presence of the triple mutant haplotype in three genotypes indicates that the species in this region may have a unique population history.

Graphical Abstract

Background

Target-site resistance to pyrethroids in Aedes aegypti, also known as knockdown resistance (kdr), is an autosomal, incompletely recessive trait [1]. Voltage-sensitive sodium channel (Vssc) mutations at codons 1016 and 1534 occur in Ae. aegypti within the pyrethroid receptor sites in Domains II (S6) and III (S6) of the protein molecule [2]. A third mutation, S989P, which is often in perfect linkage with V1016G, is not known to reduce the sensitivity of the sodium channel [2], but confers some additive pyrethroid resistance in the homozygous state in combination with 1016G [3].

Pyrethroid resistance and Vssc mutations have been identified in Ae. aegypti from the Kingdom of Saudi Arabia and discussed in several studies [4,5,6,7]. Three mutation sites are being screened routinely in samples of Ae. aegypti from Jeddah and are labeled as S989P, V1016G and F1534C according to the codon numbering in the sequence of the most abundant splice variant of the house fly, Musca domestica, Vssc (GenBank accession nos. AAB47604 and AAB47605) [8]. However, direct links between these mutations and resistance phenotypes in this locality are still unclear.

Dengue has been a public health problem in Saudi Arabia and particularly in Jeddah since its discovery there in 1994 [9]. Dengue incidence is confined to the western and south western regions of Saudi Arabia, and cases have increased in Jeddah and Makkah from 2008 to 2012, possibly as a result of increased stagnant water due to construction, but also because of international travel associated with the Hajj [10]. Dengue cases peaked in 2013 with 6512 cases, and central Jeddah has been identified as a site of major dengue concentration [10]. Currently, dengue suppression relies on control of the mosquito vector, Ae. aegypti, through the applications of insecticides of which the pyrethroid group has played a dominant role [5]. Bifenthrin (Type I) and alpha-cypermethrin (Type II) are pyrethroids that have commonly been used in Jeddah.

Insecticide resistance is prevalent in Ae. aegypti around the world [11, 12], and the resulting reduction of insecticide efficacy is likely to be exacerbating the impact of dengue and of chikungunya, a second flavivirus that circulates in Saudi Arabia [10]. Insecticide resistance threatens the long-term utility of vector control in this region. Hence, an alternative approach based on blocking dengue transmission in mosquitoes by the endosymbiont Wolbachia is being trialed following the success of this approach in suppressing dengue in other countries [13, 14]. Preparation for a release of Wolbachia-infected mosquitoes for control of dengue transmission has commenced in the Al Safa and Al Rawabi districts of Jeddah in the Kingdom of Saudi Arabia. Aljamaa Municipality, in which Al Rawabi is located, had the highest numbers of both Ae. aegypti and Culex mosquitoes out of all Municipalities of Jeddah over a 2-month collection period in 2015–2016 [15]. Al Safa is located in the Almatar Municipality, where moderate to high numbers of Ae. aegypti were caught in the same study [15].

Knowledge of the insecticide resistance status of mosquito populations in Jeddah is one prerequisite for establishing a Wolbachia-based dengue control program. The presence and frequency of the mutants can be used as a measure of resistance in the local population of Ae. aegypti that can be compared with Wolbachia-infected mosquitoes destined for release. Mosquitoes are screened for mutations in the Vssc, the target site of pyrethroid insecticides and DDT [16], to characterize the population from the field for comparison with the future population for release. Similar releases of Wolbachia mosquitoes have been shown to succeed when insecticide resistance status of the release population is similar to that in the field [13, 17] and to fail or strike difficulties if the release population of mosquitoes does not carry equivalent resistance traits [18, 19]. Although a Wolbachia program is not based on insecticide use, it is a reality that released mosquitoes will be exposed to insecticides in the field and therefore must initially be able to survive in such an environment.

In this study, we aimed:

  1. 1.

    To characterize the Vssc mutations 989/1016/1534 and haplotypes in the Ae. aegypti populations in Jeddah in two areas (Al Safa and Al Rawabi) to compare with known haplotype patterns in other parts of the world.

  2. 2.

    To compare the distribution of these Vssc mutations in dead and surviving mosquitoes from bioassays with insecticides that target the Vssc (Type I and Type II pyrethroids as well as DDT, which shows a similar mode of action to Type I pyrethroids [20]).

  3. 3.

    To screen for the presence of a fourth mutation at codon 410 known to confer pyrethroid resistance to Ae. aegypti in a restricted number of locations [21,22,23] and, if found, to characterize this mutation in bioassay dead and surviving mosquitoes by sequencing a region of Vssc Domain I where this mutation may be located.

  4. 4.

    To screen for a mutation in codon 1763, known from Ae. aegypti in Taiwan [24, 25], which does not affect the sensitivity of the Vssc [2], but may play a similar role to S989P in enhancing resistance conferred by the other mutations [24].

  5. 5.

    To use a synergist bioassay to look for evidence of a metabolic component of pyrethroid resistance in Ae. aegypti from Jeddah.

Knowledge about insecticide resistance in Ae. aegypti from the potential Wolbachia release sites is needed to ensure that similar traits are found in the strains of mosquitoes to be released to facilitate their survival under continued pressure of insecticide applications. Comparison of Vssc resistance haplotypes with those found in Ae. aegypti in other parts of the world gives preliminary information about population history of the species in Saudi Arabia. This knowledge will benefit the planned Wolbachia release programs and the development of long-term insecticide resistance management strategies for Jeddah.

Methods

Sample collection

Two study districts were chosen for their status as dengue hotspots with large populations of Ae. aegypti, along with physical attributes such as highway boundaries to facilitate implementation of a Wolbachia program. Mosquito eggs were sampled in Al Safa-9 district (21.58343 N, 39.20830 E) located in central Jeddah and Al Rawabi district (21.47574 N, 39.26798 E) in the southern part of the city (Fig. 1) on three occasions (August 2018, October 2019 and February 2020) using oviposition bucket traps. The field sampling was performed over a 2-week period for each operation, and felt pieces containing eggs were collected twice and brought to the laboratory for rearing in separate containers, from which we tested an average of 10–12 mosquito samples (larval and adult stages, F0). These samples were used for screening the field populations for Vssc mutations.

Fig. 1
figure 1

Map of Jeddah, Kingdom of Saudi Arabia, showing two mosquito sampling districts, Al Safa-9 and Al Rawabi

Mosquito colonies were initiated and maintained in the Jeddah laboratory using larvae collected from various collection points in the Safa-9 district. Mosquitoes (F0) from this wild-type colony were used for WHO bioassays. Mosquito rearing was conducted at room temperature (26–27 °C) with around 70% relative humidity. Larvae were fed with fish food (Hikari Sinking Wafers® Item # 21503, Kyorin Food Industries, Ltd., Japan).

WHO bioassays

The standard WHO insecticide resistance bioassay involving insecticide-impregnated papers was used for adult mosquitoes [26, 27]. Twenty-five female adults aged < 5 days post-eclosion and non-blood fed were used as the test subjects. Five replicates and one control were used for each bioassay. Controls comprised the relevant solvent used for the insecticide papers (silicone oil for pyrethroids). The mosquitoes were exposed in WHO tubes to impregnated paper with diagnostic concentrations of deltamethrin (0.05%), permethrin (0.75%) or DDT (4%). Bioassays to test for synergism of permethrin 0.75% were conducted by exposing mosquitoes to piperonyl butoxide (PBO) (4%) for 1 h followed by exposure to permethrin (0.75%) for an additional hour. Insecticide impregnated papers were purchased from the Vector Control Research Unit, Universiti Sains Malaysia, Penang. Knockdown was scored every 10 min for a 1-h exposure period at room temperature (26–27 °C). After exposure, the mosquitoes were transferred into clean tubes and were supplied with cotton soaked in 10% sugar solution. Mortality was recorded after 24 h. Classification of resistance status was made based on pre-determined WHO guidelines in which mortality of < 90% characterizes a resistant population. Bioassays were repeated until at least 40 dead and 40 survivors could be collected for each insecticide to be used for analysis of Vssc mutations. Mosquitoes from the bioassays were stored in absolute ethanol at − 20 °C prior to Vssc screening.

DNA extraction

DNA was extracted from adult mosquitoes or late instar larvae using the DNeasy® Blood and Tissue kit (QIAGEN Sciences, Germantown, MD, USA), the Roche High Pure PCR template kit (Roche Molecular Systems, Inc., Pleasanton, CA, USA) according to the instructions of the manufacturer or Chelex® 100 resin (Bio-Rad Laboratories Inc., Hercules, CA, USA). Two final elutions of DNA from field mosquitoes (reared from field-collected eggs, F0) were made from the kit extractions with the first being used for construction of genomic libraries for a related study and the second being used for screening of Vssc mutations after being diluted 1:5 with water. The same dilution factor was used for DNA extracted using Chelex® 100 resin on samples which were not to be used for construction of genomic libraries (dead and surviving adult mosquitoes from insecticide bioassays, reared from field collected larvae, F0).

Screening of Vssc mutations

989/1016/1534

Custom TaqMan® SNP Genotyping Assays (Life Technologies, Carlsbad, CA, USA) were developed for each of the three target site mutations (codons 989, 1016, 1534) and were run on the Roche LightCycler® 480 [28]. Endpoint genotyping was conducted using the Roche LightCycler® 480 Software, version 1.5.1.62.

One hundred thirty mosquito samples were screened from each of the two districts of Jeddah; 117 dead and 121 surviving mosquitoes were screened from insecticide bioassays (Table 2).

410 and 1763

A subset of 32 samples was screened for a mutation at codon 410 in Domain I of the Vssc and 93 samples for codon 1763 in Domain IV by PCR and Sanger sequencing. Primers used to amplify the region around codon 410 (exon 10 in Domain I) of the Vssc were aegSCF10 (5ʹ-GTGTTACGATCAGCTGGACC-3ʹ and aegSCR10 (5ʹ-AAGCGCTTCTTCCTCGGC-3ʹ) from Tancredi et al. [29]. Primers to amplify codon 1763 in Vssc Domain IV were albSCF6 (5ʹ-TCGAGAAGTACTTCGTGTCG-3ʹ) and albSCR8 (5ʹ-AACAGCAGGATCATGCTCTG-3ʹ) [30]. Amplicons were sequenced with albSCF7 (5ʹ-AGGTATCCGAACGTTGCTGT-3ʹ) by Macrogen Inc., South Korea for Sanger sequencing.

A 25-µl PCR master mix was set up as follows: 2.5 µl Standard Taq (Mg-free) Reaction Buffer (10×) (B9015 New England BioLabs Inc., Ipswich, MA, USA), 2.0 µl dNTP mix (10 mM total) (Bioline (Aust) Pty, Ltd, Eveleigh, NSW, Australia), MgCl2 (50 mM) (Bioline (Aust) Pty Ltd, Eveleigh, NSW, Australia), 1.25 µl of each primer (10 µM) (aegSCF10 and aegSCR10), 0.125 µl IMMOLASE™ DNA Polymerase (5 u/µl) (Bioline (Aust) Pty Ltd, Eveleigh, NSW, Australia), 15.125 µl PCR-grade water and 2 µl DNA (1:10 dilution of Chelex®-extracted DNA). A PCR reaction of 95 °C for 10 min, 35 cycles of 95 °C 30 s, 55 °C 45 s, 72 °C 45 s and a final extension at 72 °C for 5 min was run on an Eppendorf Mastercycler® (Eppendorf AG, Hamburg, Germany). Amplicons were observed on a 1.5% agarose gel (Bioline (Aust) Pty, Ltd, Eveleigh, NSW, Australia) run for 40 min at 90 V stained with SYBR™ Safe DNA Gel Stain (Thermo Fisher Scientific Inc., Waltham, MA, USA) and viewed with a GelDoc XR Gel Documentation System (Bio-Rad Laboratories Inc., Hercules, CA, USA). Samples were sent to Macrogen Inc., South Korea, for Sanger sequencing. Amplicons for the codon 410 screen were sequenced with aegSCR10, and those for codon 1763 were sequenced with albSCF7 (5ʹ-AGGTATCCGAACGTTGCTGT-3ʹ). Sequences were analyzed with Geneious Prime Version 2020 (Biomatters Ltd.).

Samples screened were those collected in Al Rawabi and Al Safa in February 2020. Sequences were first aligned to KY747529 [21] for verification as Aedes aegypti Vssc and then to KY747530.1 for codon 410 (mutant L) [21] and MK495874 (wildtype D) and MK495875 (mutant Y) [25] for codon 1763 for identification of a possible mutation.

Statistical analyses

Vssc mutation data from Ae. aegypti from Al Safa and Al Rawabi collected over 3 years were analyzed for site differences in genotype frequencies using contingency tables with significance tested through the Chi-squared statistic. Monte Carlo tests were used to determine significance because expected values in some cells were particularly low, even when data were combined across collection dates. These analyses were performed in IBM SPSS Statistics (IBM Corp., Armonk, NY, USA; 2013).

Odds ratios (with 95% confidence intervals) [31] were calculated to indicate the odds of an individual surviving exposure to an insecticide if it carried one particular genotype compared with another. An odds ratio of 1 indicates that there is no relationship between survival and the genotype under investigation. If 95% confidence intervals of the odds ratio do not span the value “1” then this suggests that the genotype is associated with survival.

Results

Vssc mutations in temporal field samples

Data from multiple studies of Ae. aegypti [28, 32,33,34,35] suggest that certain haplotypes of the three mutation sites predominate in a population, and there is little evidence of crossing over to disrupt the phase patterns found. Four haplotypes and eight genotypes were found overall in the samples of Ae. aegypti tested from Jeddah (Fig. 2). An initial sample of ten individuals per district identified four Vssc genotypes (A, B, C, E) seen in samples from Asia and the Indo-Pacific (Table 1, Fig. 3). One individual from Safa had a unique genotype (a homozygous mutant at codon 1016 and 989, but a heterozygote at 1534—genotype I) (Table 1). The linkage patterns we have observed for these mutations in mosquitoes from Asia and the Indo-Pacific cannot produce individuals of genotype I found in Safa (Fig. 2). To obtain this genotype, one parent had to contribute a triple mutant haplotype (H5), and the other contributed a second haplotype (H1) that we see in Asia/Indo-Pacific samples (Fig. 2).

Fig. 2
figure 2

Vssc haplotypes and genotypes in Aedes aegypti from Jeddah from the current study and [5] compared with those found in southeast Asia [17] and the Indo-Pacific [28], including Australia [39]

Table 1 Voltage-sensitive sodium channel (Vssc) mutations in Aedes aegypti from Al Safa and Al Rawabi districts, Jeddah, 2018, 2019, 2020 (order of mutations is 1016/1534/989, upper genotype is the base configuration: T = wildtype, G or C = mutant; lower genotype is the amino acid configuration)
Fig. 3
figure 3

Comparison of frequency of Vssc genotypes in Aedes aegypti from Al Safa and Al Rawabi districts, Jeddah, Kingdom of Saudi Arabia, 2018–2020 (order of mutations 1016/1534/989, T = wildtype, G or C = mutant; AA = amino acid configuration)

The second screen consisted of 60 samples from each study site (Al Rawabi and Al Safa Districts, Jeddah), collected in October 2019. Three of the Vssc genotypes (A–C—Fig. 2) found in the mosquitoes collected from these sites in Saudi Arabia were those found in samples from the Indo-Pacific region [28]. However, there was again one individual from Al Safa with genotype I (Table 1).

The third sample of Ae. aegypti was taken in February 2020, and again 60 samples from each site were screened for Vssc mutations. Genotype I was again found in mosquitoes from Al Safa (Table 1). Genotype J, which also contains the triple mutant haplotype (Fig. 2), was found in both Al Safa and Al Rawabi, and triple mutant homozygotes (genotype K—Fig. 2) were found in both Al Safa (two individuals) and Al Rawabi (one individual) (Fig. 3). Two genotypes found only in the Al Rawabi sample (D and E—Fig. 2—one individual of each) indicate that the triple wildtype haplotype is present in that district, but was not found as a homozygote. The frequency of each of the common genotypes (A, B, C) is similar between Al Safa and Al Rawabi (Table 1), but there are differences in distribution of the rare genotypes (Fig. 3). When data for all sampling dates were combined, there was no significant difference in genotypes between mosquitoes sampled in Al Safa and Al Rawabi (χ2 = 7.27, df = 7, P = 0.401, 0.388–0.413, 99% confidence intervals).

WHO bioassays

  1. 1.

    Mortality

    None of the bioassays conducted showed mortality in the range 98–100%, which would indicate susceptibility according to WHO guidelines [36]. Permethrin 0.75% induced only 44.8% mortality, and DDT 4% bioassays showed mortality ranging from 3.2 to 10.4%, confirming a high level of resistance to these compounds in Ae. aegypti from Jeddah. Four bioassays with deltamethrin 0.05% showed an average mortality just under 90%, which indicates that the population of Ae. aegypti from Jeddah is resistant to this compound as well (Table 2).

  2. 2.

    Synergist bioassay

    Odds of mosquitoes being alive if exposed to permethrin 0.75% alone are just over twice those if exposed to permethrin 0.75% + PBO 4% (Table 3). PBO 4% alone caused some mortality, but the odds of being alive are not greater if exposed to permethrin compared with exposure to PBO alone (Table 3).

  3. 3.

    Comparison of Vssc genotypes in dead and survivors

    Five Vssc genotypes (A, B, C/L, J, K) were identified in the surviving pool of mosquitoes screened, and four genotypes (A, B, C/L, J) occurred in the dead mosquitoes exposed to 0.75% permethrin (Table 4). Two triple mutants were found in the survivors. There was no significant difference (α = 0.05) in the odds of being alive after exposure to permethrin 0.75% if carrying one genotype over any other (Additional file 1: Table S1).

Table 2 WHO bioassays of Aedes aegypti from Jeddah, KSA using permethrin 0.75%, deltamethrin 0.05%, DDT 4%—% mortality and numbers in Vssc screen
Table 3 WHO bioassays with permethrin 0.75% and synergist, piperonyl butoxide (PBO) and Aedes aegypti from Jeddah, KSA (OR odds ratio, α = 0.05, * = significance)
Table 4 Frequency (%) of Vssc mutations in dead and surviving Ae. aegypti mosquitoes from Jeddah from WHO bioassays with permethrin (0.75%), deltamethrin (0.05%) or DDT (4%) (order of mutations is 1016/1534/989, upper genotype is the base configuration: T = wildtype, G or C = mutant; lower genotype is the amino acid configuration; letters in parentheses refer to genotype code assigned in Fig. 2)

Survivors of deltamethrin 0.05% fell into four genotypes as did the dead individuals, though the fourth and least frequent genotype differed between them (I for alive and J for dead) (Table 4). Odds ratios were significant (α = 0.05) for two genotype comparisons indicating that the odds of surviving were higher for individuals of genotype A than genotype B and also higher for genotype A than the triple heterozygote (genotype C or L) (Table 5).

Table 5 Odds of Ae. aegypti surviving with one Vssc genotype compared with another 24 h after a 1 h exposure to deltamethrin 0.05%

Five genotypes (A, B, C/L, l, J) occurred in both survivors of DDT 4% and in the dead individuals (A, B, C/L, I, J) (Table 4). There was no significant difference (α = 0.05) in the odds of being alive after exposure to DDT 4% if carrying one genotype over any other (Additional file 1: Table S2).

Codons 410 and 1763

Thirty-two sequences of 93–95 bp in length (30 from mosquitoes from Al Rawabi and two from Al Safa) were obtained for codon 410. Each sequence was wildtype at codon 410, coding for valine (V) (Additional file 1: Table S3). Ninety-three sequences of 116 bp in length were obtained for codon 1763 (47 from mosquitoes from Al Rawabi and 46 from Al Safa). Each sequence was identical and wildtype, coding for aspartic acid (D) at codon 1763 (Additional file 1: Table S4).

Discussion

The consequences of insecticide resistance in insect vectors of disease have severe implications for human health. Multiple studies characterise mutations in the Vssc gene of the dengue vector, Ae. aegypti, around the world and attempt to relate mutation status to local resistance phenotype (for example, [3, 32, 37]) or, in other cases, to investigate physiological effects of the mutations, which leads to identification of those that play a causal role (for example, [12, 21]). Such studies are useful in aiding prediction about which types of pyrethroid insecticide might succeed or fail in controlling the mosquito if applied alongside knowledge of metabolic resistance status.

Our study takes a different approach in that after characterizing Vssc mutations in Ae. aegypti from two districts of Jeddah, we compare the mutations with those in Ae. aegypti from other world regions by analyzing as haplotypes and discuss their utility as both comparative markers in field and laboratory populations to ensure quality control for Wolbachia release programs and as indicators of whether resistance has developed by local selection or invasion by resistant individuals.

Vssc mutations at codons 989, 1016 and 1534 are present in Ae. aegypti from two districts of Jeddah representing potential future release sites for Wolbachia mosquitoes for control of dengue. The populations are not fixed for one genotype, but genotypes A, B, C are common, and no homozygote wildtype individuals were found in the sample. A wildtype haplotype exists, at least in Al Rawabi, so it is possible that homozygous wildtype individuals (genotype F) may occur in the population, albeit rarely. There is no indication that this haplotype still exists in Al Safa, but increased sampling might still find it. No instances of a mutation at codons 410 or 1763 were observed in the samples screened, indicating lack of independent selection or contact of the mosquito populations with those from Taiwan [24], Brazil [21], Mexico [23] or West and Central Africa [22].

The triple mutant homozygote (1016G/1534C/989P) (genotype K—Fig. 3) can now be confirmed from Al Safa and Al Rawabi districts, having been found previously only as a heterozygote in this region [5], suggesting that Ae. aegypti from Saudi Arabia may have a unique population history which could be further explored in a full genomic analysis. The implications of the triple mutant and genotypes I–L on pyrethroid resistance levels are not well understood, but their presence provides an opportunity to study effects of these genotypes in more detail in bioassays and fitness experiments. Although the triple mutant genotype is rare in the populations of Ae. aegypti sampled from Jeddah (2–3%), the frequency is higher than that reported from Myanmar (0.98%) [34] and is high enough to facilitate collection of adequate numbers to isolate individuals for crossing experiments which could allow further characterization of the H5 haplotype and I–L genotypes in bioassays. Two triple homozygous mutants were found in the live pool of our bioassay with permethrin 0.75%.

The Vssc genotypes found in Al Safa and Al Rawabi may all be constructed from haplotypes 1, 2, 3 and 5. Haplotype 4, found only in Taiwan of the countries sampled in the Indo-Pacific [28] and in Indonesia [3], is not found in the Saudi Arabian samples we have screened, which is consistent with results of Al Nazawi et al. [5]. It is possible that the triple heterozygous genotype in Saudi Arabia could be composed of either H1 + H2 (usual condition in the Indo-Pacific—genotype C) or H3 + H5 (triple wildtype plus triple mutant—genotype L) (Fig. 2). Both H3 and H5 are rare compared with H1 and H2, but there is a possibility that some of the heterozygotes have this rare configuration. It is not known whether there would be any difference in insecticide resistance in mosquitoes which have one or the other putative haplotype combination of the triple heterozygote.

Potential implications of the genotypes A–E and I–K on control of Ae. aegypti with pyrethroid insecticides have been surmised (Table 6) [5, 32, 38], but not all genotypes have been tested for effects on susceptibility. It is likely that efficacy of both Type I and Type II pyrethroids is compromised by many of these modifications to the Vssc, though genotypes have differential effects. As well as categorizing effects on resistance of each of these genotypes, a useful area of future research would be to pinpoint their fitness costs on Ae. aegypti both in wild and Wolbachia-infected populations.

Table 6 Pyrethroid resistance implications of voltage-sensitive sodium channel (Vssc) mutations in Aedes aegypti from the Kingdom of Saudi Arabia (Al Safa and Al Rawabi districts of Jeddah)

Although resistance to three insecticides which target the Vssc was detected in bioassays in the Ae. aegypti from Jeddah in this study, the only clear association between Vssc mutations and resistance was established for the Type II pyrethroid, deltamethrin, with genotype A conferring a survival advantage compared with genotypes B (1534C) or C (triple heterozygote). Genotype A is known to reduce the sensitivity of the Vssc for both Type I and Type II pyrethroids [2], but in our study there was no obvious effect for Type I (permethrin 0.75%) or DDT. It may be that this lack of effect is concentration related and a different concentration of permethrin or DDT would reveal a segregation of genotypes between dead and surviving mosquitoes as reported from the Jazan region of Saudi Arabia [4]. In a similar study in Ae. aegypti from Malaysia [17], a survival advantage over wildtype was conferred to mosquitoes by genotypes B and C when exposed to permethrin 0.25%.

Comparison of Vssc mutations in Ae. aegypti throughout the world and their association with phenotypic resistance levels reveals that their effects have some dependence on context which may arise from differences in genetic background. It is important to note that target-site resistance may not be the only mechanism of resistance to pyrethroids that has been selected in the Jeddah populations of Ae. aegypti. A role of metabolic resistance in Ae. aegypti from Al Safa and Al Rawabi is likely. Al Nazawi et al. [5] saw an increase in mortality of Ae. aegypti from Jeddah with deltamethrin, using the synergist (PBO), which indicates that detoxification by oxidases (cytochrome P450 mixed function oxidase system) accounts for some of the pyrethroid resistance response. We also found an indication of such synergism in bioassays with the Type I pyrethroid, permethrin 0.75%.

To ensure similarity in genetic background in mosquitoes from the field in Jeddah with those containing Wolbachia and reared in the laboratory for release, a backcrossing scheme is being employed. Male mosquitoes from the field release sites are crossed with females from the Wolbachia colony. At least three rounds of backcrossing will be completed before the colony is closed for mass-rearing. The backcrossing has the function of introducing resistance mechanisms into the release population. As well as Vssc mutations and metabolic resistance mechanism against pyrethroids, mechanisms for resistance to other chemical groups would also be incorporated, if present. Vssc resistance allele status will be monitored in the release population, and further backcrossing will be carried out as required to maintain their frequency. Screening of Vssc mutations provides a point of reference for when another round of backcrossing would be beneficial, as we have shown that these mutations can be lost over time in a closed population [17].

Conclusions

The continuing presence of the wildtype haplotype and minor difference in genotypes between mosquitoes in Al Safa and Al Rawabi suggest that selection for pyrethroid resistance may be patchy between districts in Jeddah. Changes in genotype frequencies over time also suggest that selection for pyrethroid resistance may be ongoing. A future area of our research, in which we look at fine-scale population genomics of Ae. aegypti in our study areas of Jeddah, will provide further information as to whether the resistance in the populations has arisen from local selection or by invasion of resistant mosquitoes from other regions. The presence of the triple mutant, which is only found in mosquito populations from a few regions throughout the world, suggests that local selection for resistance is likely in Jeddah. Our results provide preliminary insight into differences between Ae. aegypti populations in Jeddah and those from some other parts of the world. They add to the information available for ongoing monitoring of resistance particularly with implementation of resistance management programs and indicate Vssc genotypes required in Wolbachia mosquito release populations to ensure homogeneity with the target field population, a prerequisite for establishment and subsequent reductions in transmission of dengue.

Availability of data and materials

All data generated or analyzed during this study are included in this published article [and its additional information files].

References

  1. Chang C, Huang X-Y, Chang P-C, Wu H-H, Dai S-M. Inheritance and stability of sodium channel mutations associated with permethrin knockdown resistance in Aedes aegypti. Pestic Biochem Physiol. 2012;104(2):136–42.

    Article  CAS  Google Scholar 

  2. Du Y, Nomura Y, Satar G, Hu Z, Nauen R, He SY, et al. Molecular evidence for dual pyrethroid-receptor sites on a mosquito sodium channel. Proc Natl Acad Sci. 2013;110(29):11785–90.

    Article  CAS  Google Scholar 

  3. Wuliandari J, Lee S, White V, Tantowijoyo W, Hoffmann A, Endersby-Harshman N. Association between three mutations, F1565C, V1023G and S996P, in the voltage-sensitive sodium channel gene and knockdown resistance in Aedes aegypti from Yogyakarta, Indonesia. Insects. 2015;6(3):658.

    Article  Google Scholar 

  4. Dafalla O, Alsheikh A, Mohammed W, Shrwani K, Alsheikh F, Hobani Y, et al. Knockdown resistance mutations contributing to pyrethroid resistance in Aedes aegypti population, Saudi Arabia. East Mediterr Health J. 2019;25:905–13.

    Article  Google Scholar 

  5. Al Nazawi AM, Aqili J, Alzahrani M, McCall PJ, Weetman D. Combined target site (kdr) mutations play a primary role in highly pyrethroid resistant phenotypes of Aedes aegypti from Saudi Arabia. Parasit Vectors. 2017;10(1):161.

    Article  Google Scholar 

  6. Alsheikh AA, Mohammed WS, Noureldin EM, Daffalla OM, Shrwani YA, Hobani KJ, et al. Studies on Aedes aegypti resistance to some insecticides in the Jazan District, Saudi Arabia. J Egypt Soc Parasitol. 2016;46(1):209–16.

    PubMed  Google Scholar 

  7. Aziz AT, Dieng H, Hassan AA, Satho T, Miake F, Salmah MRC, et al. Insecticide susceptibility of the dengue vector Aedes aegypti (Diptera: Culicidae) in Makkah City, Saudi Arabia. Asian Pac J Trop Dis. 2011;1(2):94–9.

    Article  Google Scholar 

  8. Kasai S, Komagata O, Itokawa K, Shono T, Ng LC, Kobayashi M, et al. Mechanisms of pyrethroid resistance in the dengue mosquito vector, Aedes aegypti: target site insensitivity, penetration, and metabolism. PLoS Negl Trop Dis. 2014;8(6):e2948.

    Article  Google Scholar 

  9. Fakeeh M, Zaki AM. Virologic and serologic surveillance for dengue fever in Jeddah, Saudi Arabia, 1994–1999. Am J Trop Med Hyg. 2001;65(6):764–7.

    Article  CAS  Google Scholar 

  10. Al-Tawfiq JA, Memish ZA. Dengue hemorrhagic fever virus in Saudi Arabia: a review. Vector-borne Zoonotic Dis. 2018;18(2):75–81.

    Article  Google Scholar 

  11. Smith LB, Kasai S, Scott JG. Pyrethroid resistance in Aedes aegypti and Aedes albopictus: important mosquito vectors of human diseases. Pestic Biochem Physiol. 2016;133:1–12.

    Article  CAS  Google Scholar 

  12. Du Y, Nomura Y, Zhorov BS, Dong K. Sodium channel mutations and pyrethroid resistance in Aedes aegypti. Insects. 2016;7(4):1–11.

    Article  Google Scholar 

  13. Nazni WA, Hoffmann AA, NoorAfizah A, Cheong YL, Mancini MV, Golding N, et al. Establishment of Wolbachia strain wAlbB in Malaysian populations of Aedes aegypti for dengue control. Curr Biol. 2019;29(24):4241-8.e5.

    Article  CAS  Google Scholar 

  14. Ryan P, Turley A, Wilson G, Hurst T, Retzki K, Brown-Kenyon J, et al. Establishment of wMel Wolbachia in Aedes aegypti mosquitoes and reduction of local dengue transmission in Cairns and surrounding locations in northern Queensland, Australia [version 2; peer review: 2 approved]. Gates Open Res. 2020;3:1547.

    Article  Google Scholar 

  15. Khan MA, Elhossary S, Khan IA, Zahrani MHA, Zahrani FSA, Bashri FMA. The impact of climatic variables with GIS application on the abundance of medically important mosquitoes (Diptera: Culicidae) in Jeddah, Saudi Arabia. Int J Mosq Res. 2018;5(5):12–8.

    CAS  Google Scholar 

  16. Davies TGE, Field LM, Usherwood PNR, Williamson MS. DDT, pyrethrins, pyrethroids and insect sodium channels. IUBMB Life. 2007;59(3):151–62.

    Article  CAS  Google Scholar 

  17. Ahmad NA, Endersby-Harshman NM, Mohd Mazni NR, Mohd Zabari NZ, Amran SN, Ridhuan Ghazali MK, et al. Characterization of sodium channel mutations in the dengue vector mosquitoes Aedes aegypti and Aedes albopictus within the context of ongoing Wolbachia releases in Kuala Lumpur, Malaysia. Insects. 2020;11(8):529.

    Article  Google Scholar 

  18. Garcia GdA, Sylvestre G, Aguiar R, da Costa GB, Martins AJ, Lima JBP, et al. Matching the genetics of released and local Aedes aegypti populations is critical to assure Wolbachia invasion. PLoS Negl Trop Dis. 2019;13(1):e0007023.

    Article  Google Scholar 

  19. Garcia GA, Hoffmann AA, Maciel-de-Freitas R, Villela DAM. Aedes aegypti insecticide resistance underlies the success (and failure) of Wolbachia population replacement. Sci Rep. 2020;10(1):63.

    Article  CAS  Google Scholar 

  20. Du Y, Nomura Y, Zhorov BS, Dong K. Evidence for dual binding sites for 1,1,1-trichloro-2,2-bis(p-chlorophenyl)ethane (DDT) in insect sodium channels. J Biol Chem. 2016;291(9):4638–48.

    Article  CAS  Google Scholar 

  21. Haddi K, Tomé HVV, Du Y, Valbon WR, Nomura Y, Martins GF, et al. Detection of a new pyrethroid resistance mutation (V410L) in the sodium channel of Aedes aegypti: a potential challenge for mosquito control. Sci Rep. 2017;7:46549.

    Article  CAS  Google Scholar 

  22. Ayres CFJ, Seixas G, Borrego S, Marques C, Monteiro I, Marques CS, et al. The V410L knockdown resistance mutation occurs in island and continental populations of Aedes aegypti in West and Central Africa. PLoS Negl Trop Dis. 2020;14(5):e0008216.

    Article  CAS  Google Scholar 

  23. Villanueva-Segura OK, Ontiveros-Zapata KA, Lopez-Monroy B, Ponce-Garcia G, Gutierrez-Rodriguez SM, Davila-Barboza JA, et al. Distribution and frequency of the kdr mutation V410L in natural populations of Aedes aegypti (L.) (Diptera: Culicidae) from eastern and southern Mexico. J Med Entomol. 2020;57(1):218–23.

    Article  CAS  Google Scholar 

  24. Chang C, Shen W-K, Wang T-T, Lin Y-H, Hsu E-L, Dai S-M. A novel amino acid substitution in a voltage-gated sodium channel is associated with knockdown resistance to permethrin in Aedes aegypti. Insect Biochem Mol Biol. 2009;39(4):272–8.

    Article  CAS  Google Scholar 

  25. Chung H-H, Cheng IC, Chen Y-C, Lin C, Tomita T, Teng H-J. Voltage-gated sodium channel intron polymorphism and four mutations comprise six haplotypes in an Aedes aegypti population in Taiwan. PLoS Negl Trop Dis. 2019;13(3):e0007291.

    Article  CAS  Google Scholar 

  26. WHO. Instructions for determining the susceptibility or resistance of mosquito larvae to insecticides. World Health Organization. 1981. http://www.who.int/iris/handle/10665/69615.

  27. WHO. Monitoring and managing insecticide resistance in Aedes mosquito populations: interim guidance for entomologists. World Health Organization. 2016. http://www.who.int/iris/handle/10665/204588.

  28. Endersby-Harshman NM, Schmidt TL, Chung J, van Rooyen A, Weeks AR, Hoffmann AA. Heterogeneous genetic invasions of three insecticide resistance mutations in Indo-Pacific populations of Aedes aegypti (L.). Mol Ecol. 2020;29:1628–41.

    Article  CAS  Google Scholar 

  29. Tancredi A, Papandrea D, Marconcini M, Carballar-Lejarazu R, Casas-Martinez M, Lo E, et al. Tracing temporal and geographic distribution of resistance to pyrethroids in the arboviral vector Aedes albopictus. PLoS Negl Trop Dis. 2020;14(6):e0008350.

    Article  CAS  Google Scholar 

  30. Kasai S, Ng LC, Lam-Phua SG, Tang CS, Itokawa K, Komagata O, et al. First detection of a putative knockdown resistance gene in major mosquito vector, Aedes albopictus. Jpn J Infect Dis. 2011;64(3):217–21.

    CAS  PubMed  Google Scholar 

  31. Sokal RR, Rohlf FJ. Biometry: the principles and practice of statistics in biological research. 3rd ed. New York: W.H. Freeman; 1995.

    Google Scholar 

  32. Plernsub S, Saingamsook J, Yanola J, Lumjuan N, Tippawangkosol P, Sukontason K, et al. Additive effect of knockdown resistance mutations, S989P, V1016G and F1534C, in a heterozygous genotype conferring pyrethroid resistance in Aedes aegypti in Thailand. Parasit Vectors. 2016;9(1):417.

    Article  Google Scholar 

  33. Ishak IH, Jaal Z, Ranson H, Wondji CS. Contrasting patterns of insecticide resistance and knockdown resistance (kdr) in the dengue vectors Aedes aegypti and Aedes albopictus from Malaysia. Parasit Vectors. 2015;8(1):1–13.

    Article  Google Scholar 

  34. Kawada H, Oo SZM, Thaung S, Kawashima E, Maung YNM, Thu HM, et al. Co-occurrence of point mutations in the Voltage-Gated Sodium Channel of pyrethroid-resistant Aedes aegypti populations in Myanmar. PLoS Negl Trop Dis. 2014;8(7):1–8.

    Article  Google Scholar 

  35. Stenhouse SA, Plernsub S, Yanola J, Lumjuan N, Dantrakool A, Choochote W, et al. Detection of the V1016G mutation in the voltage-gated sodium channel gene of Aedes aegypti (Diptera: Culicidae) by allele-specific PCR assay, and its distribution and effect on deltamethrin resistance in Thailand. Parasit Vectors. 2013;6(1):253.

    Article  Google Scholar 

  36. WHO. Test procedures for insecticide resistance monitoring in malaria vector mosquitoes. Geneva: World Health Organisation; 2013.

    Google Scholar 

  37. Li C-X, Kaufman PE, Xue R-D, Zhao M-H, Wang G, Yan T, et al. Relationship between insecticide resistance and kdr mutations in the dengue vector Aedes aegypti in Southern China. Parasit Vectors. 2015;8(1):325.

    Article  Google Scholar 

  38. Hirata K, Komagata O, Itokawa K, Yamamoto A, Tomita T, Kasai S. A single crossing-over event in voltage-sensitive Na+ channel genes may cause critical failure of dengue mosquito control by insecticides. PLoS Negl Trop Dis. 2014;8(8):e3085.

    Article  Google Scholar 

  39. Endersby-Harshman NM, Wuliandari JR, Harshman LG, Frohn V, Johnson BJ, Ritchie SA, et al. Pyrethroid susceptibility has been maintained in the dengue vector, Aedes aegypti (Diptera: Culicidae), in Queensland, Australia. J Med Entomol. 2017;54(6):1649–58.

    Article  CAS  Google Scholar 

  40. Abbott W. A method of computing the effectiveness of an insecticide. J Econ Entomol. 1925;18:265–7.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

We are deeply grateful for the outstanding logistical support of the municipality of Jeddah in providing workforce, management, and local expertise during the field operations in Jeddah districts. We also thank Dr. Abdulghaffar Azhari for team management and all members involved in DNA extractions and bioassays. We thank Dan Pagendam for producing the map of the study sites in Jeddah. Photography in the graphical abstract is courtesy of Perran Stott-Ross.

Funding

The study was funded under the KACST-CSIRO co-investment collaborative research agreement [ETSC&KACST-CSIRO-2018-12-30-21] on “Management strategies of vector-borne disease in Saudi Arabia: feasibility of the Wolbachia-based approach as an alternative to chemical pesticides.”

Author information

Authors and Affiliations

Authors

Contributions

NMEH, AAH and SE designed the study and wrote the paper. AEA conducted all WHO bioassays. NMEH conducted DNA extraction, SNP genotyping and sample preparation for DNA sequencing. NMEH and AAH analyzed the data. BA, MAA, MBAF, AAM and MSA provided research services for the study and reviewed the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Nancy M. Endersby-Harshman.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

Frequency of Vssc mutations in dead and surviving Ae. aegypti mosquitoes from WHO insecticide paper bioassays with Type I pyrethroid, permethrin (0.75%) (OR = odds ratio with 95% confidence intervals) (order of mutations is 1016/1534/989, genotype given is the base configuration: T = wildtype, G or C = mutant). Table S2. Frequency of Vssc mutations in dead and surviving Ae. aegypti mosquitoes from WHO insecticide paper bioassays with DDT (4%) (OR = odds ratio with 95% confidence intervals) (order of mutations is 1016/1534/989, genotype given is the base configuration: T = wildtype, G or C = mutant). Table S3. DNA sequence around Vssc Domain I codon 410 in Aedes aegypti from Jeddah, Kingdom of Saudi Arabia. Table S4. DNA sequence around Vssc Domain IV codon 1763 in Aedes aegypti from Jeddah, Kingdom of Saudi Arabia (codon 1763 highlighted in yellow) aligned to GenBank sequences MK495874 (Ae. aegypti D1763—wildtype) and MK495875 (Ae. aegypti 1763Y—mutant).

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Endersby-Harshman, N.M., Ali, A., Alhumrani, B. et al. Voltage-sensitive sodium channel (Vssc) mutations associated with pyrethroid insecticide resistance in Aedes aegypti (L.) from two districts of Jeddah, Kingdom of Saudi Arabia: baseline information for a Wolbachia release program. Parasites Vectors 14, 361 (2021). https://doi.org/10.1186/s13071-021-04867-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13071-021-04867-3

Keywords