Frequency of kdr mutations in the voltage-sensitive sodium channel (VSSC) gene in Aedes aegypti from Yogyakarta and implications for Wolbachia-infected mosquito trials

DOI: https://doi.org/10.21203/rs.3.rs-17263/v3

Abstract

Background : In the inner city of Yogyakarta, Indonesia, insecticide resistance is expected in the main dengue vector, Aedes aegypti , because of the intensive local application of pyrethroid insecticides. However, detailed information about the nature of resistance in this species is required to assist the release of Wolbachia mosquitoes in a dengue control program so that we can ensure that insecticide resistance in the strain of Ae. aegypti being released matches that of the background population.

Methods: High-resolution melt genotyping was used to screen for kdr mutations associated with pyrethroid resistance in the voltage-sensitive sodium channel (V SSC ) gene in Ae. aegypti of some areas in the inner city of Yogyakarta.

Results: The results show that the V1016G mutation predominated, with individuals homozygous for the 1016G allele at a frequency of 82.1% and the mutant allele G at a frequency of 92%. Two patterns of co-occurrence of mutations were detected in this study, homozygous individuals V1016G/S989P; and heterozygous individuals V1016G/F1534C/S989P. We found the simultaneous occurrence of kdr mutations V1016G and F1534C at all collection sites, but not within individual mosquitoes. Homozygous mutants at locus 1016 were homozygous wildtype at locus 1534 and vice versa, and heterozygous V1016G were also heterozygous for F1534C. The most common tri-locus genotype co-occurrences were homozygous mutant 1016GG and homozygous wild-type FF1534, combined with homozygous mutant 989PP (GG/FF/PP) at a frequency of 38.28%.

Conclusions: Given the relatively small differences in frequency of resistance alleles across the city area, locality variations in resistance should have minor implications for the success of Wolbachia mosquito trials being undertaken in the Yogyakarta area.

Background

Dengue is the most important mosquito-borne viral disease in the world [1]. Dengue is also a significant public health problem in Indonesia since the first reported dengue outbreak in 1968. All four dengue virus (DENV; Flavivirus) serotypes are detected, suggesting hyperendemicity in all of 34 provinces in Indonesia [2, 3]. Dengue virus is transmitted by A. aegypti, a mosquito that has also been linked with the transmission of other major arboviral diseases including chikungunya (CHIKV) and Zika (ZIKV). The virologically confirmed chikungunya was initially reported during an outbreak in Jambi in 1982 [4].  After a near 20 years of absence, twenty-four distinct outbreaks of probable chikungunya reemerged in South Sumatera, Aceh and West Java in the early 2000s [4]. To date, there has been no recorded fatality related to CHIKV infection in Indonesia [5]. Zika virus (ZIKV, family Flaviviridae) has become one of the global public health threats because of its association with Guillain-Barré syndrome and microcephaly. Until recently, most evidence for Zika virus infection in Asia, including in Indonesia, has been serologic [6]. 

 

Indonesia has made dengue a notifiable disease. On average, 136,670 DENV infections and 1,112 deaths were reported annually from 2013 to 2016, an incidence of about 54 dengue fever or dengue hemorrhagic fever cases per 100,000 population and a case fatality rate of approximately 1% [7]. Dengue cases in Indonesia are expected to be under-reported due to poor disease surveillance and a low level of reporting [8].  The surveillance database likely covers dengue probable cases with supportive dengue serology or with epidemiologic linkage, and dengue confirmed cases with confirmatory laboratory criteria [8] [9]. Dengue infections deprived of diagnostic tests may not end up in the surveillance database [8]. Despite the level of uncertainty on total case number, dengue visibly results in a considerable cost to the health sector, and a heavy economic and social impact is likely.

 

Indonesia has made progress in many areas of dengue prevention and control. On August 2016, Indonesia became the sixth country to approve the first licensed dengue vaccine, Dengvaxia® [10, 11]. To date, the vaccine so far has been approved in 20 countries and launched in 11 including Indonesia [11]. However in December 2017, the National Agency of Drug and Food (Badan POM) suspended the use of dengue vaccine in Indonesia over safety concerns [12] and later in February 2018, the BPOM revised recommendation that the use of vaccine be limited for people 9 – 16 years of age who had a dengue infection prior to vaccination and that the vaccine should not be administered to people who have never had a dengue infection before [13]. There is no rapid, reliable test for previous dengue infection existing, so the Dengvaxia® vaccine cannot be widely used [14]. Nevertheless, effective vector control methods are still essential, targeting A. aegypti in its immature and adult stages. It has been argued that even if the current vaccine is highly targeted and low-cost, sustained mosquito control will remain cost-effective [15].

 

The mainstay of the current vector control program in Indonesia is environmental management, which in recent years, has emphasized community participation to reduce container breeding sites [2, 16]. The country is renowned for its 3M plus campaigns, aimed at covering and cleaning water containers and burying discarded water containers, complemented with biological approaches using natural predators or pathogens as alternatives for vector control [2, 16]. The efficacy of community-based approaches is measured by a larva free index (percentage of houses free from A. aegypti larvae and pupae infestation); unfortunately, data from the last five years show a larvae free index ranging from 24% - 80%, less than the target of 95% [7].

 

Chemical insecticides still play a central role in dengue vector control in Indonesia. Vector control during outbreaks depends primarily on chemical insecticides to effect a rapid reduction in the number of infected mosquitoes and to break the dengue transmission cycle [2, 16, 17, 18]. Since the 1970s, the organophosphates malathion and temephos have been widely used to control dengue, and starting in the 1980s, dengue vector control has been highly reliant on pyrethroids. Pyrethroids are also widely used in public health for prevention and control of other mosquito-borne diseases such as malaria and filariasis, and as agricultural insecticides [19]. Pyrethroids are also approved for household protection against dengue [20, 21]. The abundant and prolonged use of pyrethroids has led to the development of resistance in A. aegypti populations in many countries including Indonesia. Resistance to pyrethroids, based on bioassay data, has been reported in some areas of Indonesia [22, 23, 24, 25, 26].

 

Two main mechanisms for pyrethroid resistance have been identified in A. aegypti, metabolic resistance and target site resistance. Metabolic resistance occurs when increased levels or modified activities of one or more detoxifying enzymes result in a more rapid detoxification of the insecticide, preventing the insecticide from reaching its target in the nervous system. Metabolic resistance involves three groups of enzymes: esterases, multi-function oxidases P450 and glutathione s-transferases (GST) [27, 28, 29]. Limited studies about metabolic resistance based on biochemical assays are available for Indonesian dengue vectors [26, 30, 31, 32].

 

Pyrethroids act on the insect nervous system, targeting the voltage sensitive sodium channel (VSSC). They modify the gating kinetics of the channel by slowing both the activation and the inactivation, stimulating the nerve cells to produce repetitive discharge that lead to paralysis and death of insects, an effect known as knockdown. Target site resistance is caused by point mutations in the VSSC gene, resulting in amino acid substitutions that affect pyrethroid binding sites, being known as kdr (knockdown resistance) mutations[33, 34]. A total of 12 point mutations in the Vssc of A. aegypti have been identified to be associated with kdr resistance to pyrethroids. Only five of these mutations have been functionally confirmed to reduce the sensitivity of mosquito sodium channels to pyrethroids including S989P, I1011M, V1016G/I, F1534C, and recently V410L [35]The mutations L982W, G923V, I1011M, and V1016G were the first reported sodium channel mutations in permethrin/DDT-resistant populations of A. aegypti from various countries [36], located in the domain IIS5 (for L982W and G923 mutations) and IIS6 (for I1011M and V1016G mutations) of the VSSC [36]. Further studies have reported novel mutations, including I1011V and V1016I in Latin American populations [37]. In Asian countries, the F1534C mutation (in domain IIIS6) was detected in A. aegypti mosquitoes from Thailand [38]and Vietnam [39]. The S989P mutation, located in linker between domain IIS5–S6, was first reported in Ae aegypti populations in Thailand [40], the D1763Y mutation (in linker IVS5-S6) was reported in Taiwan population [41], and mutation T1520I (in domain IIS6) was detected in an Indian population [42]. In 2017, the mutation V410L was firstly identified in Brazilian strains of A. aegypti [43], and in 2018, a novel mutation V419L was found in populations from Colombia [44] (both mutations are located in domain IS6).

 

Some of these mutations, V1016G, F1534C, and S989P, are widely distributed and detected in pyrethroid-resistant populations in Southeast Asian countries including Thailand, Indonesia, Malaysia, Singapore, Vietnam, Cambodia, and Laos [45]. Co-occurrence of kdr mutations has been a common phenomenon observed and, for some combinations, has been shown to confer a higher level of resistance than singly occurring mutations [46]. In A. aegypti from Southeast Asia, at least three patterns of mutational associations have been identified: V1016G/F1534C, V1016G/S989P, and V1016G/F1534C/S989P. Co-occurrence of V1016G and F1534C was reported in populations from Thailand [47], Myanmar [46], Malaysia [48], and Indonesia [49, 50, 51, 52]. Meanwhile, co-occurrence of V1016G and S989P point mutations was detected in Thailand [47], Myanmar [46], Indonesia, [49, 50, 52], and Papua New Guinea [53]. In addition, co-occurrence of triple mutations V1016G/F1534C/S989P in heterozygous form has been identified commonly in A. aegypti from Thailand [54], Myanmar [46], and Indonesia [49] [50, 52]. However, co-occurrence of triple homozygous point mutations (homozygous mutation for each of V1016G, F1534C and S989P) is very scarce having only been reported from Myanmar at a frequency of 0.98% [46] and in Indonesia in one individual (at a frequency of 0.34%) [50]

 

Despite all of the strategies implemented in Indonesia, the existing methods of controlling dengue have limited success. The development of a novel strategy of vector control that can be incorporated into the existing vector control strategy is essential. One of the innovative approaches to preventing transmission of dengue virus involves introduction of strains of the bacterium Wolbachia into A. aegypti, which has both life-shortening effects on the mosquito and direct transmission-blocking effects on dengue virus. This new strategy together with more traditional approaches for vector control including insecticide application may provide promising results to reduce dengue transmission. With the support of communities and approval from regulators, Indonesia’s first field trial of Wolbachia-infected mosquitoes began in Yogyakarta in 2014 [55]. The first field trial in two areas in the outer city of Yogyakarta has yielded results that point to local invasion [56], and release of Wolbachia mosquitoes on a broader scale in the inner-city areas of Yogyakarta has since been initiated in August 2016 [57].

 

To assist the spread of Wolbachia, insecticide resistance in the strain of A. aegypti being released needs to match that of the background population to ensure that released individuals persist and reproduce, allowing a Wolbachia invasion to take place [58]In the inner city of Yogyakarta, resistance is expected because of the local heavy application of chemical insecticides. It is possible that insecticide usage is higher in these areas than in the outer rim of Yogyakarta, as inner city areas have more dengue cases [59] and are more densely populated [60]. In our previous study, we screened samples from Yogyakarta outer areas for kdr mutations and obtained a high frequency of kdr alleles in samples collected from the outer area of Yogyakarta [49]. The most common kdr co-occurrence was V1016G homozygous mutant / F1534 wildtype, combined either with S989P heterozygous, or S989P wild type, or S989P homozygous mutant respectively. This co-occurrence was more common in the surviving than in the dead mosquitoes [49] and  a strong association between the V1016G mutation with pyrethroid resistance (type I and type II) was confirmed. In contrast, F1534C mutant homozygotes were rare and there was only a weak association between heterozygote individuals for the F1534C mutation and resistance to a type I pyrethroid. The S989P mutation, in addition to the V101G mutation was found to have an additive effect in resistance to Type II pyrethroids [49]. In this paper, we have screened Yogyakarta inner city areas for kdr mutations, based on a total of 1314 individuals collected from 27 localities (Figure 1, Table 1) and aim to compare frequency and occurrence of kdr mutations between the inner and outer city.

Methods

Mosquito samples

A. aegypti samples (male and female adults) were collected in June – July 2015 using BioGent-Sentinel (BG-S) traps (Biogents AG, Regensburg, Germany), each locality was covered by several BG-traps. Samples were collected every week and identified to species, and A. aegypti were preserved in absolute ethanol.

 

HRM (High-resolution melt) genotyping of kdr mutations in A. aegypti samples from Yogyakarta

A total of 1293 samples from 27 collection sites were screened for kdr mutations at 1016, 1534, and 989 following the HRM assay protocol of our previous study [49]. The PCR primer pairs, the region amplified in Vssc, and the product size are shown in Table 2.

We used the Roche “High Pure PCR Template Preparation” kit (Mannheim, Germany; Cat. No. 11 796 828 001) to extract mosquito DNA in a final volume of 200 μL of elution buffer and prepared a 10-fold dilution of the template DNA. PCR reactions (10 μL) contained 2 μL of 1 in 10 diluted template DNA, 12.5 μL 0.4 μL each of the primers at 10 μM, 1 μL of the ThermoPol reaction buffer (NEB Inc., Ipswich, MA, USA; Cat. No. B9004S), 0.064 μL of dNTPs at 25 mM (Bioline, Alexandria, NSW, Australia; Cat. No. BIO-39029), 0.4 μL of MgCl2 (50 mM) (Bioline, Alexandria, NSW, Australia; Cat. No. BIO-21047), 0.25 μL of the LightCycler® 480 High Resolution Melting Master (Roche, Mannheim, Germany; Cat. No. 04909631001), 0.01 μL of IMMOLASETM DNA polymerase (10 u/μL) (Bioline, Alexandria, NSW, Australia; Cat. No. BIO-21047) and 5.476 μL of ddH2O (Honeywell, Burdick and Jackson; Muskegon, MI, USA; Cat. No. 365-4).

Amplification was performed on a Roche LightCycler® 480 system (384-well format) using temperature cycling conditions of: 95°C for 10 min, 20 cycles of 95°C for 5 s, 65°C (reduce 0.5°C each cycle) for 15 s, 72°C for 15 s, followed by an additional 20 cycles of 95°C for 5 s, 55°C for 15 s, and 72°C for 15 s. Fluorescence information was captured at the end of each 72°C step. Following amplification, PCR products were subjected to HRM analysis. The HRM step involved heating the PCR products to 95°C for 1 min, cooling to 40°C for 20 s and then increasing the temperature to 65°C, and then melting at 0.05°C/s with continuous acquisition of fluorescence as the temperature increased from 65 to 95°C. Melt curves were generated in the Gene Scanning module of the Roche LightCycler® 480 software package. Melting data were normalized, temperature shifted, and displayed as derivative curves compared to the known genotyped samples. We used the same PCR protocol and thermocycling conditions, and HRM step for all kdr mutations except the parameter settings for melt curve normalization.

 

Statistical Analysis

Statistical comparisons of locations for genotypic frequencies for the V1016G, F1534C and S989P loci were undertaken using a 𝜒2 contingency test (IBM SPSS version 22) with degrees of freedom dictated by the number of localities being compared. We also tested for spatial signature in the data using the average GPS points for the BG traps in a locality and then computing Moran’s I using the software Spatial Analysis in Macroecology [61].

Hardy–Weinberg equilibrium for observed genotyped frequencies for each kdr mutation was calculated using GenAlEx 6.5 (available from http://biology.anu.edu.au/GenAlEx/) [62, 63]. Linkage disequilibrium between the two kdr mutations, V1016G and F1534C, was also tested. Data were pooled across adjacent sites to increase sample size for detecting linkage disequilibrium.

We identified three patterns of mutational associations: nearly all homozygous mutant individuals 1016GG occurred in conjunction with homozygous wild-type 1534FF; homozygous mutant individuals 1534CC occurred with the wild-type 1016VV; and heterozygous individuals for V1016G were commonly heterozygous for F1534C. We used the Pearson correlation (MS Excel Ver. 16.23) to measure the strength of relationship of these mutational associations.

Results

HRM genotyping of kdr mutations in A. aegypti samples from Yogyakarta

A total of 1314 of A. aegypti samples were collected from 27 different sites in the inner city of Yogyakarta (Figure 1, Table 1). We genotyped 1293 out of 1314 samples collected for kdr mutations at V1016G, F1534C and S989P sites using the high-resolution melt (HRM) protocol [49] and compared these according to their geographic location (Figure 1). V1016G was the most widespread kdr point mutation and was detected in all collection sites. 82.06% individuals had the homozygous mutant genotype 1016GG, 16.86% were heterozygous 1016VG and the remaining 1.08 % were homozygous wild-type 1016VV (Table 3). The allelic frequencies of V1016 and 1016G were 0.105 and 0.895 respectively. All collection sites displayed a high frequency of the mutant allele 1016G, ranging from 0.760 to 1.000. The observed frequency of the mutant allele 1016G was statistically different between collection sites as indicated by a chi-squared contingency test (χ2(26) = 44.145, p = 0.015). For the 1016 allele frequencies, Moran’s I was not significant for any distance class, indicating locality differences unrelated to distance.

The second most common kdr mutation detected in our samples was S989P. The frequencies of genotypes detected overall are as follows: homozygous wild-type 989SS = 14.54%, heterozygous 989SP = 47.18%, and homozygous mutant 989PP = 38.28%. Allelic frequencies overall were 0.599 for S989, and 0.401 for 989P. The mutant allelic frequency 989P ranged from 0.360 to 0.700 across localities, and overall there was a statistically significant difference between localities (chi-squared contingency test χ2(26) = 43.473, p = 0.017). For the 989 allele frequencies, Moran’s I did not show a significant departure from randomness at the smallest distance class (0.85 km) but there was a positive association at 1.7 km and negative association at the next class up (2.2 km) suggesting a complex pattern not related directly to distance.

For the 1534 site, the homozygous wild-type genotype 1534FF was the most common at 83.06%, while the homozygous mutant genotype 1534CC was scarce at 1.00%, and heterozygous genotype 1534F/C was at 15.94%. A contingency test indicates no significant difference between the localities in allele frequencies (χ2(26) = 32.151, p = 0.188). Moran’s I was not significant for any distance classes for the 1534 alleles.

All genotypes at the three loci in all field collections were in Hardy – Weinberg equilibrium except in Co1 (Rejowinangun) (Table 3). When combined across sites, genotype frequencies were also in Hardy-Weinberg equilibrium (Table 4).

 

Mutation combinations

We observed that there were nine out of 27 possible genotype combinations (three genotypes at three loci) in our 1293 individual mosquitoes (Figure 2). Figure 2 shows the frequency of each of the nine tri-locus combinations. The most common tri-locus genotypes were the homozygous mutant 1016GG and homozygous wild-type FF1534, combined with either the homozygous mutant 989PP (GG/FF/PP) (38.28%), the heterozygous S989P (GG/FF/SP) (35.43%), or the homozygous wild-type SS989 (GG/FF/SS) (8.35%). These genotypes totalled 495, 458 and 108 mosquitoes respectively. Triple mutations were found only in the heterozygous genotype state (VG/FC/SP) at a frequency of 11.14% (144 individuals).

The V1016G homozygote mutant genotype was found alone (GG/FF/SS) in 108 (8.35%) individuals, and there were 8 (1.00%) individuals exhibiting the homozygous mutant genotype for F1534C (VVCCSS). S989P was only found in conjunction either with V1016G (as a homozygote or heterozygote) or F1534C (only as a heterozygote). We only obtained one wild type individual for the three mutations (VV/FF/SS), likewise there was no individual that expressed triple mutations (homozygous mutants for each site V1016G, F1534C, and S989P) among the 1293 mosquito individuals analysed.

Interestingly, we found that nearly all homozygous mutant 1016GG individuals occurred in conjunction with homozygous wild-type 1534FF, and homozygous mutant individuals 1534CC occurred with the wild-type 1016VV. Heterozygous individuals for V1016G were commonly heterozygous for F1534C, although there was a small number of combined homozygous/heterozygous genotypes (Table 3, Figure 2). A Pearson's correlation was run to assess these relationship patterns and there was a significant positive correlation between co-occurrence of 1016G homozygous mutants and F1534 wild type at all collection sites (r = 0.999, N = 27, p < 0.0005), with 1016G homozygous mutants explaining 99.8% of the variation in V1016G mutations (Figure 3). A strong positive correlation was also found between the co-occurrence of individuals with heterozygous V1016G and heterozygous F1534C (r = 0.988, N = 27, p <0.0005) (Figure 4). The most common genotype for V1016C and F1534C was the homozygous mutant genotype 1016GG combined with the homozygous wild-type 1534FF with a frequency of 82.1%, followed by the double heterozygous 1016VG with 1534FC which had a frequency of 9.74%, and finally the homozygous wild-type 1016VV combined with the homozygous mutant 1534CC which had a frequency of 1.01%. However, there were 12 individuals with heterozygous 1016VG that were combined with homozygous wild-type 1534FF. Table 4 shows the linkage disequilibrium coefficients (Rij), χ2 and associated probabilities obtained between pairwise loci. To calculate linkage disequilibrium, samples from adjacent collection sites were pooled. In all pooled collection sites, there was very strong linkage disequilibrium between V1016G and F1534C, with coefficients (Rij) ranged from 0.863 to 1.

We also compared V1016G, F1534C, and S989P from those sites where data were available from our previous study in 2012 [49] with this study in 2015. Most allele frequencies were similar between the two sampling periods and did not differ significantly by contingency tests, but there was a significant difference at Cokrodiningratan for S989 (X2(1) = 9.88, P = 0.002) where frequency of the mutant 989P increased (Table 5).

Discussion

The diagnostic identification of the kdr mutations in A. aegypti natural populations is nowadays an important tool to predict resistance to pyrethroids in the field. In this study, we genotyped kdr mutations (F1534C, V1016G, and S989P) that have been confirmed as being associated with pyrethroid resistance in the Vssc of A. aegypti [49, 50, 64, 65]. The homozygous V1016G mutation predominated, with individuals homozygous for the 1016G allele being at a frequency of 82.1% and the mutant allele (G) being at a frequency of 92%. The high frequency of the mutant allele is similar to previous findings for A. aegypti samples from the outer region of Yogyakarta [49], and some other provinces in Indonesia, including Denpasar (70% in a resistant group and 50% in susceptible groups [52], Jakarta (65% in a resistant group and 35% in a susceptible group) [51], and Central Java [50]. In Asia, a high frequency of the mutant allele V1016G is also common, including in Myanmar where it is at 84.4% [46], and in southern China where it is at 95.5% [66]. The V1016G mutation was detected at a lower frequency in Thailand [47], Vietnam [46], and Malaysia [48].

 

For locus 1534, the homozygous wild-type genotype was frequent at all collection sites. A high frequency of wild-type was also found in A. aegypti mosquitoes from Jakarta (70.70%), Central Java (91.16%) [50], and Denpasar (63.91%) [52]. The 1534C allele frequency in this study (8%) was lower compared to the neighbouring countries of southeast Asia, including Malaysia (40-80%) [48]Vietnam (22%) [46], Myanmar (21%) [46], and Thailand (51%) [67].

 

In the present study, there were three patterns of co-occurrence of point mutations, namely, V1016G/F1534C, V1016G/S989P, and V1016G/F1534C/S989P. We found the simultaneous occurrence of kdr mutations V1016G and F1534C in all collection sites. Homozygous mutants at locus 1016 had homozygous wildtype alleles at locus 1534 and vice versa, and heterozygous V1016G were also heterozygous for F1534C, apart from a small number of combined homozygous/heterozygous genotypes. The main co-occurrence was homozygous mutants at 1016 and homozygous wildtype alleles at 1534 (1016G/F1534 haplotype) at a frequency of 82.1%. Pearson’s correlation analysis strongly supports statistical associations between mutations in all collection sites. Correspondingly, a similar linkage association is also found in the A. aegypti population from our previous study where these genotypes occurred at 69.2% [49], and from Myanmar where the equivalent frequency was 59.3% [46]. In contrast, the more common association in A. aegypti from Thailand was homozygous F1534C mutation with V1016G wild-type at frequency 43.5% [47]. Double mutants of 1016G and 1534C were never found. This is consistent with other studies, which have noted that in natural populations of A. aegypti, the double homozygote for V1016G and F1534C mutations are rare [46, 47, 48, 49] or absent [47], suggesting a possible fitness cost related to this haplotype or absence of recombination to date to bring homozygotes for V1016G and F1534C mutations onto the same haplotype [54]. However, a small number of mosquitoes homozygous for 1016G and 1534C occur in Malaysia [48] and Myanmar [46]. Diverse haplotypes associated with V1016G and F1534C co-occurrences across the countries suggest that selection has acted differently across sites perhaps related to pyrethroid usage, or else that different genetic backgrounds have had an impact on evolutionary trajectories [68]. Given the linkage disequilibrium between V1016G and F1534C, there is a non-random association of alleles at the two kdr loci which may indicate a selection process acting simultaneously at these two loci.

 

The co-occurrence of the S989P mutation with the V1016G mutation has been noted in previous work from some other Asian countries [46, 66, 68, 69, 70]. However, it also appears that the S989P mutation can be found in the absence of the V1016G mutation. A small number of homozygous and heterozygous mutants at S989P were found to coexist with the F1534C mutation instead of the V1016G mutation in A. aegypti samples from Yogyakarta outer city areas [49]. and Central Java [50]. Additionally, the S989P mutation has been reported as occurring alone in samples from Central Java [50]. The most common tri-locus genotype co-occurrences were homozygous mutant 1016GG and homozygous wild-type FF1534, combined with homozygous mutant 989PP (GG/FF/PP); followed by heterozygous S989P (GG/FF/SP). These results are somewhat different from our previous findings, which detected the most common combination being the homozygous mutant 1016G and homozygous wild-type F1534, combined with heterozygous S989P which was recorded at 39.75% in the resistant group, 26.83% in the susceptible group, and 31% in the field collection [49]. These differences may reflect local selection processes for resistance that differ between the city and outer areas of Yogyakarta.

 

The spatial variation in resistance allele frequencies found in this study may have minor implications for Wolbachia mosquito trials being undertaken in the Yogyakarta area [56]. At present, some areas of the city have been invaded by Wolbachia whereas other areas are acting as uninvaded controls. The expectation is that dengue incidence will be lower in invaded areas because there is no local transmission or at least reduced transmission where Wolbachia is common, based on the fact that the wMel strain released in A. aegypti mosquitoes partially blocks dengue transmission [71]. However, if there are differences in pesticide resistance across areas, then this may confound interpretations because increased resistance may be associated with a higher local density of mosquitoes which could influence transmission as well. This problem will be alleviated to some extent by backcrossing the Wolbachia mosquito release strain to the local wild mosquitoes, so that the local insecticide resistance genetic profile is introgressed to the release strain [56]. Given the relatively small differences in frequency of resistance alleles, we do not anticipate large confounding effects of locality variation in resistance, but it is nevertheless worthwhile controlling for such an effect when comparing impacts of Wolbachia.

 

Finally, the data suggest that there is no spatial structure to variation in resistance mutations across the city area. This seems to point to local selection events driving resistance frequencies rather than migration at this scale. Still, any attempts to manage pesticide resistance such as through rotating active ingredients will need to be mindful of selection pressures that can locally increase resistance and aim to alter behaviour across the entire inner-city area.

Conclusion

This study identified kdr mutations in Aedes aegypti from some areas in the inner city of Yogyakarta managed for Wolbachia releases. We detected the highest frequency of V1016G mutation, and the lowest frequency of F1534C mutation. The low frequency of the F1565C mutations in the samples may indicate that Type I pyrethroids have not been used as extensively as Type II. We also observed the pattern of the simultaneous occurrence of kdr mutations V1016G and F1534C: the V1016/C1534 double homozygote (VV/CC), the G1016/C1534 double heterozygote (VG/FC) and the G1016/F1534 homozygote (GG/FF), and a small amount of combined heterozygous/ homozygous (VG/FF) genotypes. Given the linkage disequilibrium between V1016G and F1534C, there is a non-random association of alleles at the two kdr loci which may indicate a selection process acting simultaneously at these two loci. The most common tri-locus genotype co-occurrences were homozygous mutant 1016GG and homozygous wild-type FF1534, combined with homozygous mutant 989PP (GG/FF/PP). This is a bit different from the co-occurrence from Yogyakarta outer city samples which was homozygous mutant 1016G and homozygous wild-type F1534, combined with heterozygous S989P. These differences are perhaps due to their having different local selection processes for resistance between the city and outer areas of Yogyakarta. The relatively low spatial variation in resistance allele frequency across

Declarations

Ethics approval and consent to participate

No ethics approval was required for this study.

 

Consent for publication

Not applicable.

 

Availability of data and material

Data are presented in the main manuscript.

 

Competing interests

The authors declare that they have no competing interests.

 

Funding

The funding body had no role in the design of the study and collection, analysis, and interpretation of data or in writing the manuscript.

 

Authors' contributions

JRW, NMEH and AAH designed the study. WT collected the samples. JRW collected the data. JRW and AAH analysed the data. JRW, NMEH and AAH wrote the manuscript. All authors read and approved the final manuscript.

 

Acknowledgements

The study was funded by the National Health and Medical Research Council with a Programme Grant/Award Numbers: 1037003, 1132412 and a Fellowship to Ary A. Hoffmann, NHMRC Fellowship Grant no. 1118640; FNIH Vector-based Control of Transmission: Discovery Research (VCTR) program Eliminate Dengue; JJT Foundation and the Tahija Foundation.

References

  1. WHO: Global strategy for dengue prevention and control 2012–2020. France: World Health Organization; 2014.
  2. Indonesia KKR. Demam berdarah dengue di Indonesia tahun 1968 - 2009. Buletin Jendela Epidemiologi. 2010;2 Agustus:1 - 14. https://www.kemkes.go.id/download.php?file=download/pusdatin/buletin/buletin-dbd.pdf.
  3. Indonesia KKR: Situasi DBD di Indonesia. Edited by Indonesia KKR. Jakarta: Pusdatin Kementrian Kesehatan Republik Indonesia; 2016.
  4. Laras K, Sukri NC, Larasati RP, Bangs MJ, Kosim R, Djauzi, et al. Tracking the re-emergence of epidemic chikungunya virus in Indonesia. Trans R Soc Trop Med Hyg. 2005;99 2:128-41; doi: 10.1016/j.trstmh.2004.03.013. https://doi.org/10.1016/j.trstmh.2004.03.013.
  5. Harapan H, Michie A, Mudatsir M, Nusa R, Yohan B, Wagner AL, et al. Chikungunya virus infection in Indonesia: a systematic review and evolutionary analysis. BMC Infect Dis. 2019;19 1:243; doi: 10.1186/s12879-019-3857-y. https://doi.org/10.1186/s12879-019-3857-y.
  6. Yudhaputri FA, Trimarsanto H, Perkasa A, Yohan B, Haryanto S, Wiyatno A, et al. Genomic characterization of Zika virus isolated from Indonesia. Virology. 2017;510:248-51. https://doi.org/10.1016/j.virol.2017.07.026.
  7. Indonesia KKR. Profil kesehatan Indonesia tahun 2016. Jakarta: Kementerian Kesehatan Republik Indonesia; 2017.
  8. Karyanti MR, Uiterwaal CSPM, Kusriastuti R, Hadinegoro SR, Rovers MM, Heesterbeek H, et al. The changing incidence of dengue haemorrhagic fever in Indonesia: a 45-year registry-based analysis. BMC Infec Dis. 2014;14 1:1-7; doi: 10.1186/1471-2334-14-412. http://dx.doi.org/10.1186/1471-2334-14-412.
  9. WHO. Comprehensive guidelines for prevention and control of dengue and dengue haemorrhagic fever. Revised and expanded edition. India: World Health Organization Regional Office for South-East Asia; 2011.
  10. RI BP: Badan POM menyetujui izin edar vaksin dengue di Indonesia. http://pionas.pom.go.id/artikel/badan-pom-menyetujui-izin-edar-vaksin-dengue-di-indonesia (2017). Accessed May 20 2018.
  11. Pasteur S: First dengue vaccine approved in more than 10 countries. Lyon, France: Sanofi Pasteur; 2016.
  12. RI BP: Penjelasan Badan POM RI terkait isu keamanan vaksin dengue (demam berdarah). https://www.pom.go.id/new/view/more/klarifikasi/72/PENJELASAN-BADAN-POM-RI--TERKAIT-ISU-KEAMANAN-VAKSIN-DENGUE--DEMAM-BERDARAH-.html (2017). Accessed May 20 2018.
  13. RI BP: Penjelasan Badan POM RI terkait isu keamanan vaksin dengue (demam berdarah). http://www.pom.go.id/new/view/more/klarifikasi/79/PENJELASAN-BADAN-POM-RI-----TERKAIT--ISU-KEAMANAN-VAKSIN-DENGUE--DEMAM-BERDARAH-.html (2018). Accessed May 20 2018.
  14. "The dengue vaccine dilemma". Lancet Infect Dis. 2018;18 2:123. Editorial; doi: 10.1016/S1473-3099(18)30023-9. http://dx.doi.org/10.1016/S1473-3099(18)30023-9.
  15. Fitzpatrick C, Haines A, Bangert M, Farlow A, Hemingway J, Velayudhan R. An economic evaluation of vector control in the age of a dengue vaccine. PLoS Negl Trop Dis. 2017;11 8:e0005785; doi: 10.1371/journal.pntd.0005785. https://doi.org/10.1371/journal.pntd.0005785.
  16. RI KK. Petunjuk teknis implementasi PSN 3M-Plus dengan gerakan 1 rumah 1 jumantik. Jakarta: Kementerian Kesehatan RI; 2016.
  17. RI KK: Informasi umum demam berdarah dengue. 2011.
  18. RI KK. Modul pengendalian demam berdarah dengue. Jakarta: Direktorat Jenderal Pengendalian Penyakit Dan Penyehatan Lingkungan; 2011.
  19. Indonesia MoH: Review of national vector control policy in Indonesia. http://www.actmalaria.net/files/vector_control/i_resistance/Indonesia.pdf (2010). Accessed 27 March 2012 2012.
  20. Pratamawati DA, Irawan AS, Widiarti W. Relationship between knowledge of vector with household insecticide usage behavior in dengue hemorrhagic fever endemic areas in Bali Province. Vektora. 2012;4 2:99-116. http://ejournal.litbang.depkes.go.id/index.php/vk/article/view/3503.
  21. Setyawati DR: Pengetahuan, sikap, dan tindakan penggunaan pestisida rumah tangga di Jakarta dan Surabaya. http://repository.ipb.ac.id/handle/123456789/27154 (2010). Accessed 25 January 2016.
  22. Ikawati B, Sunaryo S, Widiastuti D. Peta status kerentanan Aedes aegypti (Linn.) terhadap insektisida cypermethrin dan malathion di Jawa Tengah. Aspirator. 2015;7 1 Juni:23-8. http://dx.doi.org/10.22435/aspirator.v7i1.3722.23-28.
  23. Widiarti W, Heriyanto B, Boewono DT, Widyastuti U, Mujiono M, Lasmiati L, et al. Peta resistensi vektor demam berdarah dengue Aedes aegypti terhadap insektisida kelompok organofosfat, karbamat dan pyrethroid di Propinsi Jawa Tengah dan Daerah Istimewa Yogyakarta. Buletin Penelitian Kesehatan. 2011;39 4:176-89. http://ejournal.litbang.kemkes.go.id/index.php/BPK/article/view/54.
  24. Ahmad I, Astari S, Rahardjo B, Tan M, Munif A: Resistance of Aedes aegypti from three provinces in Indonesia, to pyrethroid and organophosphate insecticides. In: International Conference on Mathematics and Natural Sciences (ICMNS)2006; Bandung, Indonesia.
  25. Mantolu Y, Kustiati K, Ambarningrum TB, Yusmalinar S, Ahmad I. Status dan perkembangan resistensi Aedes aegypti (Linnaeus) (Diptera: Culicidae) strain Bandung, Bogor, Makassar, Palu, dan VCRU terhadap insektisida Permetrin dengan seleksi lima generasi. Jurnal Entomologi Indonesia. 2016;13 1:1- 8; doi: 10.5994/jei.13.1.1. http://dx.doi.org/10.5994/jei.13.1.1.
  26. Putra RE, Ahmad I, Prasetyo DB, Susanti S, Rahayu R, Hariani N. Detection of insecticide resistance in the larvae of some Aedes aegypti (Diptera: Culicidae) strains from Java, Indonesia to temephos, malathion and permethrin. Int J Mosq Res. 2016;3 3:23-8. http://www.dipterajournal.com/vol3issue3/pdf/3-2-4.1.pdf.
  27. Hemingway J. The molecular basis of two contrasting metabolic mechanisms of insecticide resistance. Insect Biochem Mol Biol. 2000;30 11:1009-15. https://doi.org/10.1016/S0965-1748(00)00079-5.
  28. Hemingway J, Ranson H. Insecticide resistance in insect vectors of human disease. Ann Rev Entomol. 2000;45 1:371-91; doi: doi:10.1146/annurev.ento.45.1.371. https://doi.org/10.1146/annurev.ento.45.1.371.
  29. Hemingway J, Hawkes NJ, McCarroll L, Ranson H. The molecular basis of insecticide resistance in mosquitoes. Insect Biochem Mol Biol. 2004;34 7:653-65; doi: 10.1016/j.ibmb.2004.03.018. https://doi.org/10.1016/j.ibmb.2004.03.018.
  30. Ahmad I, Astari S, Tan M. Resistance of Aedes aegypti (Diptera: Culicidae) in 2006 to pyrethroid insecticides in Indonesia and its association with oxidase and esterase levels. Pakistan Journal of Biological Sciences. 2007;10 20:3688-92. http://docsdrive.com/pdfs/ansinet/pjbs/2007/3688-3692.pdf.
  31. Astari S, Ahmad I. Insecticide resistance and effect of piperonyl butoxide as a synergist in three strains of Aedes aegypti (Linn.)(Diptera: Culicidae) against insecticides permethrin, cypermethrin, and d-allethrin. Buletin Penelitian Kesehatan. 2005;33 2 73-9. http://ejournal.litbang.kemkes.go.id/index.php/BPK/article/view/210.
  32. Widiastuti D, Sunaryo S, Pramestuti N, Martini M. Aktivitas enzim monooksigenase pada populasi nyamuk Aedes aegypti di Kecamatan Tembalang, Kota Semarang. Aspirator-Journal of Vector-borne Disease Studies. 2015;7 1:1-6. http://ejournal.litbang.depkes.go.id/index.php/aspirator/article/view/3696.
  33. Soderlund DM. Pyrethroids, knockdown resistance and sodium channels. Pest management science. 2008;64 6:610-6; doi: 10.1002/ps.1574. https://doi.org/10.1002/ps.1574.
  34. Dong K, Du Y, Rinkevich F, Nomura Y, Xu P, Wang L, et al. Molecular biology of insect sodium channels and pyrethroid resistance. Insect Biochem Mol Biol. 2014;50:1-17. https://doi.org/10.1016/j.ibmb.2014.03.012.
  35. Saavedra-Rodriguez K, Maloof FV, Campbell CL, Garcia-Rejon J, Lenhart A, Penilla P, et al. Parallel evolution of VGSC mutations at domains IS6, IIS6 and IIIS6 in pyrethroid resistant Aedes aegypti from Mexico. Sci Rep. 2018;8 1:6747; doi: 10.1038/s41598-018-25222-0. https://doi.org/10.1038/s41598-018-25222-0.
  36. Brengues C, Hawkes NJ, Chandre F, McCarroll L, Duchon S, Guillet P, et al. Pyrethroid and DDT cross-resistance in Aedes aegypti is correlated with novel mutations in the voltage-gated sodium channel gene. Medical and veterinary entomology. 2003;17 1:87-94. https://doi.org/10.1046/j.1365-2915.2003.00412.x.
  37. Saavedra-Rodriguez K, Urdaneta-Marquez L, Rajatileka S, Moulton M, Flores AE, Fernandez-Salas I, et al. A mutation in the voltage-gated sodium channel gene associated with pyrethroid resistance in Latin American Aedes aegypti. Insect Mol Biol. 2007;16 6:785-98; doi: 10.1111/j.1365-2583.2007.00774.x. http://dx.doi.org/10.1111/j.1365-2583.2007.00774.x.
  38. Yanola J, Somboon P, Walton C, Nachaiwieng W, Prapanthadara L. A novel F1552/C1552 point mutation in the Aedes aegypti voltage-gated sodium channel gene associated with permethrin resistance. Pesticide Biochem Physiol. 2010;96 3:127-31; doi: 10.1016/j.pestbp.2009.10.005. https://doi.org/10.1016/j.pestbp.2009.10.005.
  39. Kawada H, Higa Y, Komagata O, Kasai S, Tomita T, Thi Yen N, et al. Widespread distribution of a newly found point mutation in voltage-gated sodium channel in pyrethroid-resistant Aedes aegypti populations in Vietnam. PLoS Negl Trop Dis. 2009;3 10:e527; doi: 10.1371/journal.pntd.0000527. https://doi.org/10.1371/journal.pntd.0000527.
  40. Srisawat R, Komalamisra N, Eshita Y, Zheng M, Ono K, Itoh TQ, et al. Point mutations in domain II of the voltage-gated sodium channel gene in deltamethrin-resistant Aedes aegypti (Diptera: Culicidae). Appl Entomol Zool. 2010;45 2:275-82; doi: 10.1303/aez.2010.275. https://doi.org/10.1303/aez.2010.275.
  41. Chang C, Shen WK, Wang TT, Lin YH, Hsu EL, Dai SM. 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; doi: 10.1016/j.ibmb.2009.01.001. https://doi.org/10.1016/j.ibmb.2009.01.001.
  42. Kushwah RBS, Dykes CL, Kapoor N, Adak T, Singh OP. Pyrethroid-resistance and presence of two knockdown resistance (kdr) mutations, F1534C and a novel mutation T1520I, in Indian Aedes aegypti. PLoS Negl Trop Dis. 2015;9 1:e3332; doi: 10.1371/journal.pntd.0003332. https://dx.doi.org/10.1371%2Fjournal.pntd.0003332.
  43. 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; doi: 10.1038/srep46549. http://dx.doi.org/10.1038/srep46549.
  44. Granada Y, Mejia-Jaramillo AM, Strode C, Triana-Chavez O. A point mutation V419L in the sodium channel gene from natural populations of Aedes aegypti is involved in resistance to lambda-cyhalothrin in Colombia. Insects. 2018;9 1; doi: 10.3390/insects9010023. https://doi.org/10.3390/insects9010023.
  45. Amelia-Yap ZH, Chen CD, Sofian-Azirun M, Low VL. Pyrethroid resistance in the dengue vector Aedes aegypti in Southeast Asia: present situation and prospects for management. Parasites Vectors. 2018;11 1:332; doi: 10.1186/s13071-018-2899-0. https://doi.org/10.1186/s13071-018-2899-0.
  46. 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:e3032; doi: 10.1371/journal.pntd.0003032. https://dx.doi.org/10.1371%2Fjournal.pntd.0003032.
  47. 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. Parasites Vectors. 2013;6 1:253; doi: 10.1186/1756-3305-6-253. https://doi.org/10.1186/1756-3305-6-253.
  48. 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. Parasites Vectors. 2015;8 1:181; doi: 10.1186/s13071-015-0797-2. https://doi.org/10.1186/s13071-015-0797-2.
  49. Wuliandari JR, Lee SF, White VL, Tantowijoyo W, Hoffmann AA, Endersby-Harshman NM. 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-85. https://doi.org/10.3390/insects6030658.
  50. Sayono S, Hidayati APN, Fahri S, Sumanto DD, E., Hadisaputro S, Asih PBS, et al. Distribution of voltage-gated sodium channel (Nav) alleles among the Aedes aegypti populations in Central Java Province and its association with resistance to pyrethroid insecticides. PLOS ONE. 2016;11 3:e0150577; doi: 10.1371/journal.pone.0150577. https://doi.org/10.1371/journal.pone.0150577.
  51. Hamid PH, Prastowo J, Ghiffari A, Taubert A, Hermosilla C. Aedes aegypti resistance development to commonly used insecticides in Jakarta, Indonesia. PLOS ONE. 2017;12 12:e0189680; doi: 10.1371/journal.pone.0189680. https://doi.org/10.1371/journal.pone.0189680.
  52. Hamid PH, Prastowo J, Widyasari A, Taubert A, Hermosilla C. Knockdown resistance (kdr) of the voltage-gated sodium channel gene of Aedes aegypti population in Denpasar, Bali, Indonesia. Parasites Vectors. 2017;10 1:283; doi: 10.1186/s13071-017-2215-4. https://dx.doi.org/10.1186%2Fs13071-017-2215-4.
  53. Demok S, Endersby-Harshman N, Vinit R, Timinao L, Robinson LJ, Susapu M, et al. Insecticide resistance status of Aedes aegypti and Aedes albopictus mosquitoes in Papua New Guinea. Parasites Vectors. 2019;12 1:333; doi: 10.1186/s13071-019-3585-6. https://doi.org/10.1186/s13071-019-3585-6.
  54. 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. Parasites Vectors. 2016;9:417; doi: 10.1186/s13071-016-1713-0. https://dx.doi.org/10.1186%2Fs13071-016-1713-0.
  55. Yogyakarta E: Communities show support for mosquito releases. http://www.eliminatedengue.com/id/progress/view/news/242/year/2014/month/1/pg/1 (2014).
  56. Tantowijoyo W, Andari B, Arguni E, Budiwati N, Nurhayati I, Fitriana I, et al. Stable establishment of wMel Wolbachia in Aedes aegypti populations in Yogyakarta, Indonesia. PLoS Negl Trop Dis. 2020;14 4:e0008157; doi: 10.1371/journal.pntd.0008157. https://doi.org/10.1371/journal.pntd.0008157.
  57. Yogyakarta EDP: Nyamuk Aedes aegypti ber-Wolbachia siap diasuh warga Kota Yogya untuk atasi DBD. Yogyakarta: Eliminate Dengue Project Yogyakarta; 2016: 2.
  58. 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; doi: 10.1371/journal.pntd.0007023. https://doi.org/10.1371/journal.pntd.0007023.
  59. Yogyakarta DKDI: Profil kesehatan Daerah Istimewa Yogyakarta tahun 2016 (data tahun 2015). Edited by Indonesia DKR. Yogyakarta: Dinas Kesehatan D.I. Yogyakarta; 2016.
  60. Yogyakarta BPDI: The population density by regency/city in D.I. Yogyakarta. https://yogyakarta.bps.go.id/dynamictable/2017/08/02/31/kepadatan-penduduk-menurut-kabupaten-kota-di-d-i-yogyakarta.html (2016).
  61. Rangel TF, Diniz-Filho JAF, Bini LM. SAM: a comprehensive application for Spatial Analysis in Macroecology. Ecography. 2010;33 1:46-50; doi: doi:10.1111/j.1600-0587.2009.06299.x. https://doi.org/10.1111/j.1600-0587.2009.06299.x.
  62. Peakall R, Smouse PE. GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research—an update. Bioinformatics. 2012;28 19:2537-9; doi: 10.1093/bioinformatics/bts460. https://doi.org/10.1093/bioinformatics/bts460.
  63. Peakall R, Smouse PE. GENALEX 6: Genetic Analysis in Excel. Population genetic software for teaching and research. vol. 6; 2006.
  64. 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. Proceedings of the National Academy of Sciences of the United States of America. 2013;110 29:11785-90; doi: 10.1073/pnas.1305118110. https://doi.org/10.1073/pnas.1305118110.
  65. 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; doi: 10.1371/journal.pntd.0003085. http://dx.doi.org/10.1371%2Fjournal.pntd.0003085.
  66. Li CX, Kaufman PE, Xue RD, Zhao MH, Wang G, Yan T, et al. Relationship between insecticide resistance and kdr mutations in the dengue vector Aedes aegypti in Southern China. Parasites Vectors. 2015;8:325; doi: 10.1186/s13071-015-0933-z. https://doi.org/10.1186/s13071-015-0933-z.
  67. Saingamsook J, Saeung A, Yanola J, Lumjuan N, Walton C, Somboon P. A multiplex PCR for detection of knockdown resistance mutations, V1016G and F1534C, in pyrethroid-resistant Aedes aegypti. Parasites Vectors. 2017;10 1:465; doi: 10.1186/s13071-017-2416-x. https://doi.org/10.1186/s13071-017-2416-x.
  68. Plernsub S, Saingamsook J, Yanola J, Lumjuan N, Tippawangkosol P, Walton C, et al. Temporal frequency of knockdown resistance mutations, F1534C and V1016G, in Aedes aegypti in Chiang Mai city, Thailand and the impact of the mutations on the efficiency of thermal fogging spray with pyrethroids. Acta Trop. 2016;162:125-32; doi: 10.1016/j.actatropica.2016.06.019. https://doi.org/10.1016/j.actatropica.2016.06.019.
  69. Srisawat R, Komalamisra N, Apiwathnasorn C, Paeporn P, Roytrakul S, Rongsriyam Y, et al. Field-collected permethrin-resistant Aedes aegypti from central Thailand contain point mutations in the domain IIS6 of the sodium channel gene (kdr). Southeast Asian J Trop Med Public Health. 2012;43 6:1380-6. https://www.tm.mahidol.ac.th/seameo/journal-43-6-2012.html.
  70. 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; doi: 10.1371/journal.pntd.0002948. http://dx.doi.org/10.1371%2Fjournal.pntd.0002948.
  71. Walker T, Johnson PH, Moreira LA, Iturbe-Ormaetxe I, Frentiu FD, McMeniman CJ, et al. The wMel Wolbachia strain blocks dengue and invades caged Aedes aegypti populations. Nature. 2011;476 7361:450-3; doi: 10.1038/nature10355. https://doi.org/10.1038/nature10355.

Tables

Table 1. Collection sites, localities and number of A. aegypti collected from inner city areas of Yogyakarta.

Sites

Localities

No. of samples

Geolocation coordinates

(longitude, latitude)

C1

Suryodinigratan

112

110.35971

-7.81972

C2

Mantrijeron

25

110.36539

-7.82453

C3

Brontokusuman

25

110.36996

-7.82397

C4

Sorosutan

124

110.37840

-7.82886

C5

Kadipaten,

Patehan

34

110.35872

110.36015

-7.80587

-7.81104

C6

Panembahan,

Keparakan,

Prawirodirjan

25

110.36486

110.37327

110.37287

-7.81063

-7.81464

-7.80581

C7

Notoparajan

Ngampilan

Ngupasan

140

110.35560

110.35577

110.36199

-7.80621

-7.79728

-7.80199

C8

Wirogunan,

Sorosutan

25

110.37810

110.37872

-7.81130

-7.81484

C9

Tahunan, Semaki

19

110.38273

-7.81056

C10

Bausasran,

Gunung Ketur,

Baciro,

Purwokinanti

25

110.37448

110.37883

110.37875

110.37382

-7.79253

-7.80142

-7.79241

-7.79717

C11

Pringgokkusuman,

Sosromenduran

21

110.36969

110.36458

-7.79702

-7.79287

C12

Suryatmajan,

Tegal Panggung

124

110.36969

110.36983

-7.79702

-7.79350

 

C13

Bumijo,

Sosromenduran,

Gowongan

25

110.35676

110.36460

110.36467

-7.78798

-7.78835

-7.78353

C14

Gowongan,

Kotabaru

25

110.36921

110.36935

-7.78485

-7.78920

C15

Cokrodiningaratan,

Terban

88

110.36914

110.37384

-7.77882

-7.77877

C16

Semaki,

Baciro

68

110.38617

110.38211

-7.79694

-7.79685

C17

Klitren,

Demangan

25

110.38212

110.38219

-7.78394

-7.78927

C18

Terban,

Klitren

25

110.37876

110.38331

-7.77835

-7.78029

R1

Bener,

Kricak,

Karangwaru

31

110.352325

110.360298

110.364609

-7.77969

-7.77916

-7.77532

R2

Tegalrejo,

Pakuncen

38

110.35648

110.35070

-7.78339

-7.79298

R3

Wirobrajan,

Patangpuluhan

49

110.35088

110.34646

-7.80139

-7.80694

R4

Gedongkiwo

32

110.35344

-7.82552

B1

Muja Muju

27

110.39291

-7.79768

B2

Muja Muju,

Warung Boto,

Pandeyan

87

110.39259

110.38804 110.38647

-7.80637

-7.81086

-7.81439

B3

Pandeyan,

Giwangan

24

110.38669 110.39137

-7.82032

-7.83290

Co1

Rejojwinangun

39

110.40043

-7.81836

Co2

Prenggan, Purbayan

36

110.39775

-7.82482

           

Note that localities are mapped in Figure 1. The clusters refer to areas that are being used in a controlled intervention design to investigate the impact of Wolbachia on dengue incidence by the World Mosquito Program (Indonesia). Areas marked by an R are not included as part of this design.

 

Table 2. Primers to amplify Vssc mutations of A. aegypti from the city of Yogyakarta, Indonesia.

Mutation

Region amplified

Sequence

Product size (bp)

V1016G

IIS6

first codon of exon 21

5’-GACAAATTGTTTCCCACCCGCACAG-3’ and 5’-AAGCAAGGCTAAGAAAAGGTTAAG-3’

52

 

 

 

 

F1534C

IIIS6

24th codon of exon 31

5’-TACCTCTACTTTGTGTTCTTCATCATC-3’and 5’-GATTCAGCGTGAAGAACGACCCG-3’

52

 

 

 

 

S989P

IIS5 S6

first codon of P region

5’-CGGGTATTATGCGGCGAGTGGATC-3’ and

5’-CCCACAAGCATACAAT CCCACATGG-3’

52

 

Table 3. V1016G, F1534C and S989P genotypes of Aedes aegypti samples from Yogyakarta and chi-squared tests for Hardy-Weinberg equilibrium.

Collection codes

Collection sites

N

kdr mutations

SS

RS

RR

Frequency of R allele (95%CI)

X2

P

C1

Suryodiningratan

110

V1016G

1

21

88

0.896 (0.848, 0.929)

0.042

0.837

 

 

 

F1534C

89

20

1

0.100 (0.067, 0.147)

0.011

0.916

 

 

 

S989P

23

47

40

0.577 (0.511, 0.641)

1.706

0.191

 

C2

Mantrijeron

25

V1016G

0

4

21

0.920 (0.812, 0.969)

0.189

0.664

 

 

 

F1534C

21

4

0

0.080 (0.032, 0.188)

0.189

0.664

 

 

 

S989P

3

17

5

0.540 (0.414, 0.670)

3.400

0.065

 

 

 

 

 

 

 

 

 

 

C3

Brontokusuman

25

V1016G

0

5

20

0.900 (0.786, 0.957)

0.309

0.664

 

 

 

F1534C

20

5

0

0.100 (0.044, 0.214)

0.309

0.664

 

 

 

S989P

6

13

6

0.500 (0.333, 0.634)

0.040

0.065

 

 

 

 

 

 

 

 

 

 

C4

Sorosutan

122

V1016G

0

9

113

0.963 (0.931, 0.981)

0.198

0.656

 

 

 

F1534C

114

8

0

0.033 (0.017, 0.063)

0.198

0.656

 

 

 

S989P

15

57

50

0.643 (0.582, 0.701)

2.258

0.133

 

 

 

 

 

 

 

 

 

 

C5

Kadipaten,

30

V1016G

0

6

24

0.900 (0.781, 0.953)

0.370

0.543

 

Patehan,

 

F1534C

24

6

0

0.100 (0.047, 0.216)

0.370

0.543

 

Panembahan

 

S989P

1

16

13

0.700 (0.575, 0.801)

2.184

0.139

 

 

 

 

 

 

 

 

 

 

C6

Panembahan,

24

V1016G

0

4

20

0.917 (0.805, 0.967)

0.023

0.878

 

Keparakan,

 

F1534C

20

4

0

0.083 (0.033, 0.196)

0.004

0.709

 

Prawirodirjan

 

S989P

2

15

7

0.604 (0.463, 0.730)

0.002

0.272

 

 

 

 

 

 

 

 

 

 

C7

Notoprajan,

136

V1016G

1

23

112

0.908 (0.868, 0.937)

0.140

0.709

 

Ngupasan,

 

F1534C

113

22

1

0.088 (0.060, 0.128)

0.140

0.709

 

Ngampilan

 

S989P

15

60

61

0.669 (0.611, 0.722)

0.851

0.356

 

 

 

 

 

 

 

 

 

 

C8

Wirogunan,

25

V1016G

2

8

15

0.760 (0.626, 0.857)

1.205

0.272

 

Sorosutan

 

F1534C

20

3

2

0.140 (0.070, 0.262)

1.205

0.272

 

 

 

S989P

6

16

3

0.440 (0.312, 0.577)

1.340

0.247

 

 

 

 

 

 

 

 

 

 

C9

Tahunan,

19

V1016G

0

3

16

0.921 (0.792, 0.921)

0.416

0.519

 

Semaki.

 

F1534C

16

3

0

0.079 (0.027, 0.208)

0.416

0.519

 

 

 

S989P

2

11

6

0.605 (0.447, 0.744)

0.012

0.914

 

 

 

 

 

 

 

 

 

 

C10

Bausasran,

25

V1016G

1

5

19

0.860 (0.738, 0.931)

n.a

 

 

Gn Ketur, Baciro,

 

F1534C

20

4

1

0.120 (0.056, 0.238)

n.a

 

 

Purwokinanti

 

S989P

11

10

4

0.360 (0.241, 0.499)

1.995

0.158

 

 

 

 

 

 

 

 

 

 

C11

Pringgokussuman,

21

V1016G

0

0

21

1.000 (0.916, 1.000)

1.574

0.210

 

Sosromenduran

 

F1534C

21

0

0

0.000 (0.000, 0.084)

1.574

0.210

 

 

 

S989P

5

7

9

0.595 (0.445, 0.730)

0.884

0.347

 

 

 

 

 

 

 

 

 

 

C12

Suryatmajan,

123

V1016G

0

25

98

0.898 (0.854, 0.930)

0.107

0.744

 

Tegal Panggung

 

F1534C

98

25

0

0.102 (0.070, 0.146)

0.107

0.744

 

 

 

S989P

13

61

49

0.646 (0.585, 0.703)

0.160

0.689

 

 

 

 

 

 

 

 

 

 

C13

Bumijo,

24

V1016G

0

3

21

0.938 (0.832, 0.979)

1.361

0.243

 

Sosromenduran,

 

F1534C

21

3

0

0.063 (0.022, 0.168)

1.361

0.243

 

Gowongan

 

S989P

7

11

6

0.479 (0.345, 0.617)

1.434

0.231

 

 

 

 

 

 

 

 

 

 

C14

Gowongan,

24

V1016G

1

4

19

0.875 (0.753, 0.941)

0.718

0.397

 

Kotabaru

 

F1534C

19

4

1

0.125 (0.059, 0.247)

1.469

0.225

 

 

 

S989P

2

14

8

0.625 (0.484, 0.748)

0.435

0.509

 

 

 

 

 

 

 

 

 

 

C15

Cokrodiningratan,

88

V1016G

1

12

75

0.921 (0.871, 0.952)

0.377

0.539

 

Terban

 

F1534C

75

12

1

0.080 (0.048, 0.129)

6.292

0.012

 

 

 

S989P

10

40

38

0.659 (0.586, 0.725)

2.231

0.135

 

 

 

 

 

 

 

 

 

 

C16

Semaki, Baciro

68

V1016G

1

12

55

0.897 (0.835, 0.938)

0.135

0.714

 

 

 

F1534C

55

12

1

0.103 (0.062, 0.165)

0.135

0.714

 

 

 

S989P

8

37

23

0.610 (0.526, 0.692)

1.408

0.235

 

 

 

 

 

 

 

 

 

 

C17

Klitren,

24

V1016G

0

5

19

0.896 (0.778, 0.955)

0.324

0.569

 

Demangan

 

F1534C

19

5

0

0.104 (0.045, 0.222)

0.324

0.569

 

 

 

S989P

6

10

8

0.542 (0.403, 0.674)

0.621

0.431

 

 

 

 

 

 

 

 

 

 

C18

Terban, Klitren

25

V1016G

0

9

16

0.820 (0.692, 0.921)

0.179

0.672

 

 

 

F1534C

16

9

0

0.180 (0.098, 0.308)

0.140

0.708

 

 

 

S989P

3

15

7

0.580 (0.442, 0.706)

0.041

0.840

 

 

 

 

 

 

 

 

 

 

R1

Karangwaru,

30

V1016G

0

5

25

0.920 (0.819, 0.964)

0.248

0.619

 

Kricak, Bener

 

F1534C

26

4

0

0.067 (0.026, 0.159)

0.153

0.696

 

 

 

S989P

4

13

13

0.605 (0.447, 0.744)

0.068

0.794

 

 

 

 

 

 

 

 

 

 

R2

Tegalrejo,

38

V1016G

2

7

29

0.855(0.759, 0.917)

2.489

0.115

 

Pakuncen

 

F1534C

29

7

2

0.145 (0.083, 0.241)

2.489

0.115

 

 

 

S989P

5

15

18

0.671 (0.660, 0.766)

0.426

0.514

 

 

 

 

 

 

 

 

 

 

R3

Wirobrajan,

46

V1016G

0

9

37

0.902 (0.825, 0.948)

0.541

0.462

 

Patangpuluhan

 

F1534C

37

9

0

0.098(0.052, 0.176)

0.541

0.462

 

 

 

S989P

6

19

21

0.663 (0.562, 0.751)

0.263

0.608

 

 

 

 

 

 

 

 

 

 

R4

Gedongkiwo

32

V1016G

0

4

28

0.938 (0.850, 0.975)

0.142

0.706

 

 

 

F1534C

28

4

0

0.063(0.025, 0.150)

0.142

0.706

 

 

 

S989P

7

9

16

0.641(0.518, 0.747)

4.847

0.028

 

 

 

 

 

 

 

 

 

 

B1

Muja Muju

27

V1016G

0

7

20

0.870 (0.756, 0.936)

0.600

0.439

 

 

 

F1534C

20

7

0

0.130 (0.064, 0.244)

0.600

0.439

 

 

 

S989P

4

12

11

0.630 (0.496, 0.746)

0.600

0.807

 

 

 

 

 

 

 

 

 

 

B2

Muja Muju,

85

V1016G

0

17

68

0.900 (0.846, 0.937)

1.049

0.306

 

Warungboto,

 

F1534C

68

17

0

0.100 (0.063, 0.154)

1.049

0.306

 

Pandeyan

 

S989P

12

41

32

0.618 (0.543, 0.687)

0.038

0.845

 

 

 

 

 

 

 

 

 

 

B3

Pandeyan,

24

V1016G

0

1

23

0.979 (0.891, 0.996)

0.011

0.917

 

Giwangan

 

F1534C

23

1

0

0.021 (0.004, 0.109)

0.011

0.917

 

 

 

S989P

0

15

9

0.688 (0.547, 0.801)

4.959

0.026

 

 

 

 

 

 

 

 

 

 

Co1

Rejowinangun

39

V1016G

4

4

31

0.846 (0.750, 0.910)

14.325

0.000

 

 

 

F1534C

33

3

3

0.115 (0.062, 0.205)

15.146

0.000

 

 

 

S989P

7

13

19

0.652 (0.543, 0.756)

2.710

0.100

 

 

 

 

 

 

 

 

 

 

Co2

Prenggan,

34

V1016G

0

6

28

0.912 (0.821, 0.959)

0.318

0.573

 

Purbayan

 

F1534C

28

6

0

0.088 (0.041, 0.179)

0.318

0.573

 

 

 

S989P

7

14

13

0.588 (0.470, 0.697)

0.765

0.382

 

 

 

 

 

 

 

 

 

 

Total

 

1293

V1016G

14

218

1061

0.915 (0.893, 0.916)

0.552

0.458

 

 

 

F1534C

1073

207

13

0.080 (0.026, 0.102)

0.721

0.396

 

 

 

S989P

190

608

495

0.612 (0.599, 0.637)

0.022

0.882

 

Table 4. Linkage disequilibrium coefficients (Rij) and χ2 (1 df) for pairwise tests of V1016G and F1534C in Aedes aegypti in pooled sites (listed in Table 1).

Pooled sites

Rij

χ2

P

C1+C2

0.979

129.331

0.001

C3+C4

0.960

135.481

0.001

C5+C6+C7

0.967

177.719

0.001

C8+C9+C10

0.863

51.336

0.001

C11+C12

1.000

144.000

0.001

C13+C14+C15

1.000

136.000

0.001

C16+C17+C18

1.000

117.000

0.001

R1+R2

0.969

63.886

0.001

R3+R4

1.000

78.000

0.001

B1+B2+B3

1.000

136.000

0.001

Co1+Co2

0.880

56.561

0.001

 

Table 5. Comparison of V1016G/F1534C/S989P in A. aegypti collected from three sites in 2012/2015 in Yogyakarta, Indonesia

Site

Year

Locus

SS

SR

RR

Frequency of R allele

Cokrodiningratan

2012

V1016

0

5

35

0.938

 

2015

V1016

1

12

75

0.921

 

2012

F1534

37

3

0

0.038

 

2015

F1534

75

12

1

0.080

 

2012

S989

12

20

8

0.450

 

2015

S989

10

40

38

0.659

 

 

 

 

 

 

 

Purwokinanti

2012

V1016

1

10

27

0.842

 

2015

V1016

1

5

19

0.860

 

2012

F1534

27

11

0

0.145

 

2015

F1534

20

4

1

0.120

 

2012

S989

18

18

4

0.488

 

2015

S989

11

10

4

0.360

 

Mantrijeron

2012

V1016

1

5

32

0.908

 

2015

V1016

0

4

21

0.920

 

2012

F1534

32

6

0

0.079

 

2015

F1534

21

4

0

0.080

 

2012

S989

14

17

9

0.438

 

2015

S989

3

17

5

0.540