Pyrethroid resistance status and co-occurrence of V1016G, F1534C and S989P mutations in the Aedes aegypti population from two dengue outbreak counties along the China-Myanmar border

Background Over the past two decades, dengue fever (DF) has emerged as a significant arboviral disease in Yunnan province, China, particularly in the China-Myanmar border area. Aedes aegypti, an invasive mosquito species, plays a crucial role in transmitting the dengue virus to the local population. Insecticide-based vector control has been the primary tool employed to combat DF, but the current susceptibility status of Ae. aegypti to commonly used insecticides is unknown. Assessment of Ae. aegypti resistance to pyrethroid insecticides and an understanding of the underlying mechanisms of this resistance in the China-Myanmar border region is of significant strategic importance for effectively controlling the DF epidemic in the area. Methods Aedes aegypti larvae collected from Ruili and Gengma counties in Yunnan Province were reared to adults in the laboratory and tested for susceptibility to three pyrethroid insecticides (3.20% permethrin, 0.08% lambda-cyhalothrin and 0.20% deltamethrin) by the standard WHO susceptibility bioassay. Genotyping of mutations in the knockdown gene (kdr), namely S989P, V1016G and F1534C, that are responsible for resistance to pyrethroid insecticides was performed using allele-specific PCR methods. A possible association between the observed resistant phenotype and mutations in the voltage-gated sodium channel gene (VGSC) was also studied. Results Aedes aegypti mosquitoes collected from the two counties and reared in the laboratory were resistant to all of the pyrethroids tested, with the exception of Ae. aegypti from Gengma County, which showed sensitivity to 0.20% deltamethrin. The mortality rate of Ae. aegypti from Ruili county exposed to 3.20% permethrin did not differ significantly from that of Ae. aegypti from Gengma County (χ2 = 0.311, P = 0.577). By contrast, the mortality rate of Ae. aegypti from Ruili County exposed to 0.08% lambda-cyhalothrin and 0.20% deltamethrin, respectively, was significantly different from that of Ae. aegypti from Gengma. There was no significant difference in the observed KDT50 of Ae. aegypti from the two counties to various insecticides. Four mutation types and 12 genotypes were detected at three kdr mutation sites. Based on results from all tested Ae. aegypti, the V1016G mutation was the most prevalent kdr mutation (100% prevalence), followed by the S989P mutation (81.6%) and the F1534C mutation (78.9%). The constituent ratio of VGSC gene mutation types was significantly different in Ae. aegypti mosquitoes from Ruili and those Gengma. The triple mutant S989P + V1016G + F1534C was observed in 274 Ae. aegypti mosquitoes (60.8%), with the most common genotype being SP + GG + FC (31.4%). The prevalence of the F1534C mutation was significantly higher in resistant Ae. aegypti from Ruili (odds ratio [OR] 7.43; 95% confidence interval [CI] 1.71–32.29; P = 0.01) and Gengma (OR 9.29; 95% CI 3.38–25.50; P = 0.00) counties than in susceptible Ae. aegypti when exposed to 3.20% permethrin and 0.08% lambda-cyhalothrin, respectively. No significant association was observed in the triple mutation genotypes with the Ae. aegypti population exposed to 3.20% permethrin and 0.20% deltamethrin resistance (P > 0.05), except for Ae. aegypti from Gengma County when exposed to 0.08% lambda-cyhalothrin (OR 2.86; 95% CI 1.20–6.81; P = 0.02). Conclusions Aedes aegypti from Ruili and Gengma counties have developed resistance to various pyrethroid insecticides. The occurrence of multiple mutant sites in VGSC strongly correlated with the high levels of resistance to pyrethroids in the Ae. aegypti populations, highlighting the need for alternative strategies to manage the spread of resistance. A region-specific control strategy for dengue vectors needs to be implemented in the future based on the status of insecticide resistance and kdr mutations. Graphical Abstract Supplementary Information The online version contains supplementary material available at 10.1186/s13071-024-06124-9.


Background
Dengue fever (DF) is a mosquito-borne disease transmitted by Aedes aegypti or Aedes albopictus mosquitoes that carry the dengue virus.It is prevalent throughout tropical and subtropical countries and regions, with particularly severe outbreaks recorded in Southeast Asia and the Western Pacific regions [1].China's Yunnan Province shares borders with Vietnam, Laos and Myanmar, with a border length of 2185.7 km with Myanmar [2].DF is a major public health concern in Myanmar, where a total of 58,461 cases were reported between 2017 and 2019 [3,4].This high incidence rate of DF in Myanmar has resulted in multiple local outbreaks in other regions, including in Ruili and Gengma counties, Yunnan Province [5].Approximately 5500 cases of DF have been reported in these counties, indicating that DF is a significant vectorborne disease with a substantial burden on public health [6][7][8][9].The dominant mosquito species inhabiting Ruili and Gengma counties is Ae.aegypti [8,[10][11][12].As a preventative measure aimed at combating the DF epidemic, pyrethroid insecticides have been applied consistently in an effort to effectively reduce the population density of Aedes mosquitoes and prevent the spread of the disease.However, this extensive and prolonged use of these insecticides has resulted in the emergence of resistance among local mosquito populations, posing a significant obstacle to the control of DF in the region.
Insecticide resistance in mosquitoes can be attributed to several mechanisms, including knockdown resistance (kdr), metabolic resistance, epidermal resistance and behavioral resistance [13].Resistance to pyrethroid insecticides in Ae. aegypti is primarily attributed to mutations in kdr genes (kdr) and to metabolic detoxification mechanisms [14].Kdr mutations result in alterations in the voltage-gated sodium channel proteins (VGSC), the target site of pyrethroids and organochlorines.These mutations reduce the binding affinity of pyrethroids to the VGSC, thereby reducing their effectiveness [15].The most commonly observed mutation associated with pyrethroid resistance, V1016G, is found in the IIS6 region of the Ae.aegypti VGSC gene [16], with mutations also found at other sites in IIS6 region, such as G923V, L982W, S989P, V1016I, A1007G and I1011V/M.Another mutation, F1534C, in the IIIS6 region of the VGSC gene, has also been identified and is associated with pyrethroid resistance [16][17][18].
The phenomenon of Ae. aegypti mosquitoes developing resistance to pyrethroid insecticides has been noted across multiple countries and regions worldwide.Progress in researching this topic has been made in both China and the Greater Mekong Subregion (GMS), with pyrethroid resistance gaining significant attention from researchers in China in recent years.Several studies have investigated the prevalence and mechanisms of resistance in different regions of China.One study conducted in southern China found high levels of pyrethroid resistance in Ae. aegypti populations and identified the S989P, V1016G and F1534C mutations [19].Another study conducted in Yunnan Province, southwest China, also reported high levels of pyrethroid resistance in Ae. aegypti populations, and identified the V1016G mutation as the primary resistance mechanism in this region.In addition to identifying single-point mutations, researchers have also detected new co-occurrence of mutations related to pyrethroid resistance in Ae.Aegypti [20,21].Several studies have shown that the presence of both V1016G and F1534C mutations in the same individual mutation sites.Based on results from all tested Ae. aegypti, the V1016G mutation was the most prevalent kdr mutation (100% prevalence), followed by the S989P mutation (81.6%) and the F1534C mutation (78.9%).The constituent ratio of VGSC gene mutation types was significantly different in Ae. aegypti mosquitoes from Ruili and those Gengma.The triple mutant S989P + V1016G + F1534C was observed in 274 Ae. aegypti mosquitoes (60.8%), with the most common genotype being SP + GG + FC (31.4%).The prevalence of the F1534C mutation was significantly higher in resistant Ae. aegypti from Ruili (odds ratio [OR] 7.43; 95% confidence interval [CI] 1.71-32.29;P = 0.01) and Gengma (OR 9.29; 95% CI 3.38-25.50;P = 0.00) counties than in susceptible Ae. aegypti when exposed to 3.20% permethrin and 0.08% lambdacyhalothrin, respectively.No significant association was observed in the triple mutation genotypes with the Ae.aegypti population exposed to 3.20% permethrin and 0.20% deltamethrin resistance (P > 0.05), except for Ae.aegypti from Gengma County when exposed to 0.08% lambda-cyhalothrin (OR 2.86; 95% CI 1.20-6.81;P = 0.02).
Conclusions Aedes aegypti from Ruili and Gengma counties have developed resistance to various pyrethroid insecticides.The occurrence of multiple mutant sites in VGSC strongly correlated with the high levels of resistance to pyrethroids in the Ae.aegypti populations, highlighting the need for alternative strategies to manage the spread of resistance.A region-specific control strategy for dengue vectors needs to be implemented in the future based on the status of insecticide resistance and kdr mutations.
or population in Southeast Asian countries results in a higher level of resistance to pyrethroids [22][23][24].
Additionally, a study found evidence of co-occurrence of triple mutations (S989P, V1016G and F1534C), which resulted in even higher resistance to pyrethroids, as observed in Ae. aegypti populations from Myanmar [25].These mutations in the VGSC gene of Ae. aegypti are believed to play a significant role in the development of pyrethroid resistance.However, it is important to note that resistance mechanisms can be complex and multifactorial, involving additional genetic variations and metabolic detoxification pathways.
In light of the development of insecticide resistance, data are asymmetric both within China and in the GSM countries.In-depth studies on insecticide resistance in Ae. aegypti are critical to controlling the spread of diseases transmitted by this mosquito species.Therefore, we have investigated the resistance status of two Ae.aegypti populations from Ruili and Gengma counties to three pyrethroid insecticides, namely permethrin, lambda-cyhalothrin, and deltamethrin.Molecular surveillance based on allele-specific PCR (AS-PCR) was utilized to confirm the presence of the commonly reported V1016G, F1534C and S989P mutations and to determine the presence of particular genetic variations or mutation combinations that contribute to pyrethroid resistance.Finally, the association between mutant and phenotype was statistically analyzed to explore the mechanisms of the resistance.The findings of the study will provide scientific evidence for the control of Ae. aegypti in Yunnan Province.

Mosquito collection and rearing
The study was conducted in two different urban areas of Yunnan province, namely Ruili County of Dehong prefecture (24.0108°N, 97.8556°E) and Gengma County of Lincang City (23.5540°N, 99.0819°E).From September to December 2022, larvae and pupae of Ae. aegypti were collected from various artificial outdoor breeding sites, such as abandoned car tires, discarded containers and cans (Fig. 1).The collected larvae were collected in plastic containers, brought back to the laboratory and reared until they emerged into adults (F 0 ).The rearing conditions consisted of a constant temperature of 26 ± 1 °C, relative humidity of 65 ± 5% and a dark:light cycles of 12:12 h.After emergence, the adult mosquitoes were identified to species level based on morphological characters.

WHO susceptibility test
The adult Ae. aegypti mosquito bioassay was conducted following the WHO susceptibility test protocol [26].A total of 20-30 non-blood-fed female Ae.aegypti mosquitoes aged 3-5 days were used in the bioassays.These adults were exposed to impregnated paper with 3.20% permethrin, 0.08% lambda-cyhalothrin and 0.20% deltamethrin, respectively (all insecticides were provided by the Institute of Infectious Disease Prevention and Control, Chinese Center for Disease Control and Prevention), for 1 h.The diagnostic doses reported in this study refer to the values established by the China Centers for Disease Control and Prevention [27,28].Following exposure, the mosquitoes were transferred to separate holding tubes, and the number of deceased mosquitoes was tallied after a 24-h period.Probit regression analysis was used to determine the time (in minutes) taken for 50% of the test mosquitoes to be knocked down (KDT 50 ) and 95% confidence interval (95% CI) of Ae. aegypti in Ruili County and Gengma County in response to different insecticides.After 24 h of recovery, the number of dead and live mosquitoes were recorded, and the mortality rate associated with each tested insecticide was calculated.Mosquitoes were considered to be deceased if they were unable to fly or had turned upside down regardless of slight leg tremors.Each insecticide was tested in three replicates, and control groups were set up.If the mortality rate in the control group was > 20%, the experiment was considered to be unsuccessful and to necessitate a restart.If the mortality rate in the control group was < 20%, the results had to be corrected using Abbott's formula.S indicates susceptible, based on a mortality rate of ≥ 98%; M indicates suspected resistant, based on a mortality rate 90% to 98%; R indicates resistant, based on a mortality rate < 90% [29].

Allele-specific PCR genotyping of kdr mutation
Following the bioassay, mosquito DNA was extracted according to the manufacturer's instructions (TianGen Biotech, Beijing; Art.No. 22011210T332).Two fragments of the VGSC gene (DII-S6 and DIII-S6) were amplified using allele-specific PCR (AS-PCR) [13,19] sets (Additional file 1: Table S1) as previously described.The thermal cycling conditions used for the experiment consisted of an initial denaturation at 94 °C for 3 min; followed by 35 cycles of denaturation at 94 °C for 30 s, annealing at 51 °C or 60 °C for 30 s (51 °C for S989P, 60 °C for V1016G and F1534C) and extension at 72 °C for 60 s; with a final extension at 72 °C for 5 min.A negative control PCR tube was run, which contained all of the components of the AS-PCR except for DNA.After we conducted a thorough analysis and made precise modifications to the various factors involved in the experiment, the AS-PCR protocol displayed exceptional efficiency in precisely distinguishing between homozygous and heterozygous mosquitoes with regards to three crucial mutations, namely S989P, V1016G, and F1534C, with the application of 1.2% agarose gel electrophoresis.

Data analysis
Microsoft Excel version 365 (Microsoft Corp., Redmond, WA, USA) was used to organize the data and calculate the mortality rates and kdr mutation rates of Ae. aegypti in response to the three types of pyrethroids.SPSS version 25.0 software (SPSS IBM, Armonk, NY, USA) was used for the statistical analysis that compared the mortality rates and kdr gene mutation rates at different sites, as well as the differences in resistance phenotype and sensitive phenotype kdr gene mutation rates of Ae. aegypti in Ruili and Gengma counties in response to different insecticides.Pearson's Chi-square test (χ 2 ) was applied to analyze differences in mortality rates and kdr gene mutation rates at different sites.Additionally, the kdr gene mutation rates for both resistant and susceptible phenotypes of Ae. aegypti to the three insecticides were assessed in Ruili and Gengma counties.P values < 0.05 were considered to indicate statistical significant difference.Odds ratios (ORs) and 95% CIs were estimated using Fisher's exact test to assess the associations between kdr mutations and resistant phenotypes; at an OR > 1.0 and P < 0.05, the kdr mutation type was considered to be significantly associated with insecticide resistance.

Species identification and sampling size
From September to December 2022, > 1000 larvae were collected from Ruili and Gengma; following emergence in the laboratory, 485 female mosquitoes were morphologically identified as Ae.aegypti and subsequently tested for insecticide susceptibility to 3.20% permethrin, 0.08% lambda-cyhalothrin and 0.20% deltamethrin, respectively.

Aedes aegypti insecticide bioassays
The mortality rates of Ae. aegypti from Ruili County upon exposure to 3.20% permethrin, 0.08% lambdacyhalothrin and 0.20% deltamethrin were 18.1%, 11.43% and 24.0%, respectively.The mortality rates of Ae. aegypti from Gengma County upon exposure to these three insecticides at the same concentrations were 21.9%, 36.1% and 100.0%, respectively.Based on the mortality rates, Ae. aegypti mosquitoes from the two counties were considered to be resistant to all of the tested pyrethroids, with the exception of Ae. aegypti from Gengma County, which showed sensitivity to 0.20% deltamethrin (Table 1).The mortality rate of Ae. aegypti from Ruili County exposed to 3.20% permethrin was not significantly different from that of Ae. aegypti from Gengma County (χ 2 = 0.311, df = 1, P = 0.577).In contrast, the mortality Table 1 The resistance status of Aedes aegypti from Ruili and Gengma counties to the three tested insecticides The KDT 50 and 95% CI of the two Ae.aegypti populations in Ruili and Gengma counties in response to different insecticides are shown in Fig. 2 and Table 1.There was no significant difference in the knockdown rates of Ae. aegypti from Ruili and those from Gengma showed following exposure to 3.20% permethrin, 0.08% lambda-cyhalothrin and 0.20% deltamethrin, respectively (Table 1).Exposure to 0.20% deltamethrin had a rapid knockdown effect on Ae. aegypti from both Gengma County and Ruili County.In contrast, the KDT 50 of Ae. aegypti from Gengma County was much higher than the KDT 50 of Ae. aegypti from Ruili County following exposure to permethrin (3.20%) and lambda-cyhalothrin (0.08%), respectively (Table 1).
Four kdr mutation types were observed in Ae. aegypti mosquitoes from Ruili and Gengma counties, resulting in the identification of 12 genotypes.Among these genotypes, the combination S989P + V1016G + F1534C was particularly prevalent, accounting for 60.8% of the observed genotypes.The dominant mutation type was SP + GG + FC (31.4%).Additionally, two genotypes of S989P + V1016G were identified, accounting for 20.8% of the observed genotypes, with PP + GG being the primary type (97.9%).Also, three genotypes of V1016G + F1534C were detected, accounting for 18.19% of observed genotypes, with GG + CC being the primary type (85.4%).Only one single mutation genotype of V1016G was identified as GG.The triple mutation of S989P + V1016G + F1534C was found to be dominant in both Ruili and Gengma counties, accounting for 64.6% and 57.4%, respectively.Althoug S989P + V1016G and V1016G + F1534C combined mutations were found at both locations, the V1016G + F1534C combination was higher in Ruili County, whereas the S989P + V1016G combination was higher in Gengma County (Fig. 3a,  b).These findings indicate that the mutation types of Ae. aegypti varied at the different locations, and that a comprehensive understanding of these mutations could Fig. 2 Knockdown times of Aedes aegypti from Ruili County (a) and Gengma County (b) following exposure to 3.20% permethrin, 0.08% lambda-cyhalothrin and 0.20% deltamethrin, respectively be helpful in devising effective strategies for controlling vector-borne diseases.

Association between kdr gene mutations and resistant phenotypes
To determine the impact of kdr mutations on resistance to the three pyrethroid insecticides, we first analyzed them separately for their associations with 3.20% permethrin, 0.08% lambda-cyhalothrin and 0.20% deltamethrin.Regarding exposure to 3.20% permethrin, we found that the frequency of the F1534C mutation in mosquitoes from Ruili County was significantly higher in resistant mosquitoes than in susceptible ones (Fisherʼs exact test, P = 0.01, OR = 7.43, 95% CI 1.71-32.29).Regarding exposure to 0.08% lambda-cyhalothrin, we found that the frequency of the F1534C mutation in mosquitoes from Gengma County was significantly higher in resistant mosquitoes than in susceptible ones (Chi-square test, χ 2 = 9.29, df = 1, P < 0.001).Overall, the F1534C mutation was found to provide increased protection against insecticides, whereas the S989P and V1016G mutations were found not to have a significant relationship with the insecticies (Table 2).Generalized linear model analysis was used to reveal the relationship between the frequencies of mutation types and resistance to each insecticide.The results showed that triple mutations conferred Ae. aegypti from Gengma County with significantly greater resistance (Chi-square test, χ 2 = 2.86, df = 1, P = 0.02) when exposed to 0.08% lambda-cyhalothrin; however, the triple mutation did not show a significant association with 3.20% permethrin and 0.20% deltamethrin resistance (P > 0.05) (Table 3).

Discussion
Aedes aegypti, commonly known as the yellow fever mosquito, has been identified as a significant transmission vector of DF since its first recorded presence in Ruili County in Yunnan Province in 2002 [30].Over the years, the species has caused numerous local DF outbreaks in border areas such as Ruili and Dehong prefectures [31].Consequently, the use of insecticides for controlling the Ae.aegypti vector is crucial for preventing DF in the region.However, the overuse and misuse of insecticides have led to the development of resistance in these mosquito populations, rendering these control measures increasingly ineffective [32,33].To tackle this issue effectively, it is crucial to have a clear understanding of the resistance status of Ae. aegypti populations in the region to the various insecticides being used.This study aimed to investigate the resistance of Ae. aegypti to the three main insecticides used in the two border counties of Ruili and Gengma in Yunnan Province.The   results showed that the Ae.aegypti populations from both Ruili and Gengma exhibited widespread resistance to all three insecticides, with the exception of Ae. aegypti from Gengma County, which showed sensitivity to 0.20% deltamethrin.The development of resistance in mosquitoes can be attributed to the extensive use of these insecticides due to outbreaks of imported and local DF and the rapid increase in population density in both regions.
To address this challenge, new strategies for vector control are urgently needed, including the use of alternative insecticides and implementation of novel approaches.The resistance of Ae. aegypti to different insecticides is an important indicator for evaluating the effectiveness of insecticides.In our study, we found that the resistance of Ae. aegypti from two different locations, Ruili and Gengma counties, to 3.20% permethrin was similar, with a mortality rate of 18.1% and 21.9%, respectively.We also found that the mortality rate of Ae. aegypti from Ruili County exposed to 0.08% lambda-cyhalothrin was significantly lower than that of Ae. aegypti from Gengma County, indicating that Ae. aegypti from Ruili County has a higher level of resistance to this insecticide than those from Gengma County.This was also observed for 0.20% deltamethrin.Previous studies have also shown that Ae. aegypti from Ruili County have developed resistance to 0.03% lambda-cyhalothrin, with a mortality rate of 20.1% [32].These results suggest that Ae. aegypti in different locations have developed varying levels of resistance to different insecticides.It should be noted that the concentration of insecticides used in various earlier studies is inconsistent, which limits the horizontal comparison of resistance levels and can only provide a reference level for resistance.Furthermore, there is a deficiency in the systematic monitoring of Ae. aegypti resistance at established monitoring points, which hinders the provision of comprehensive information on the evolution of resistance levels over time.Therefore, we recommend that local public health authorities systematically monitor insecticide resistance in different Ae.aegypti populations to provide a basis for practical decision-making regarding control measures for vector-borne diseases.
Several kdr mutations have been identified in Ae. aegypti, including V410L, G923V, L982W, S989P, A1007G, I1011M/V, V1016G/I, T1520I, F1534C/L and D1763Y [16,34].Previous research has demonstrated that specific individual mutation sites, such as V1016G and F1534C, can directly confer resistance to pyrethroids [35] and that combinations of double and triple mutations can further augment the level of resistance [16,35].Our analysis of 451 individual mosquitoes revealed mutation rates of 81.6% for S989P, 78.9% for F1534C and 100.0%for V1016G.In 2015, Li et al. [19] reported an 81.9% occurrence of the S989P mutation in Ae. aegypti in southern China, and in 2017, Shi Q-M et al. [36] found a 25.0% occurrence of the F1534C mutation in mosquitoes in Yunnan Province, China.These findings indicate that the S989P and F1534C mutations are prevalent mutations in Ae. aegypti populations in Yunnan and other regions of China.Notably, in 2017, Shi C-N et al. [28] observed the V1016G mutation in all 100 samples from a population in Yunnan Province, and in 2020, Lan et al. [37] reported a 100.0%mutation rate in four counties in Yunnan Province, with a mutation rate of 99.3% in a fifth county.These results suggest that the V1016G mutation has become a fixed trend in these populations.In addition, the distribution and origin of kdr mutations demonstrate unique global geographical characteristics.For example, the kdr mutation combination S989P + V1016G was found exclusively in Asia, while the F1534C mutation was present in samples from almost all continents [38,39].Understanding the distribution and origin of kdr mutations is crucial for developing effective strategies to combat insecticide resistance and control vector-borne diseases [40].
Our study found four mutation types and 12 genotypes in the tested Ae. aegypti specimens, with triple, double and single mutations accounting for 60.8%, 36.8% and 0.2% of all mutations, respectively.The most common genotypes were SP + GG + FC and PP + GG.According to Moyes et al. [41], the presence of two mutations (S989P + V1016G/V1016G + F1534C) could enhance the resistance of Ae. aegypti to pyrethroids.The presence of a triple mutation (S989P + V1016G + F1534C) has been linked to an even higher insecticide resistance, with the authors of one study reporting a resistance level that was 28-fold [42] and 90-fold that of the sensitive strain [24].Additionally, high frequencies of 1520I + 1534C and T1520 + 1534C in Ae. aegypti from Pakistan have been found to confer resistance to pyrethroids [43].These results suggest that multiple mutations have a strong synergistic effect and that combined mutations are more effective than a single mutation in conferring resistance to pyrethroids.Monitoring changes in the category and intensity of resistance is crucial for responding to DF epidemics and for their prevention and control, as well as for providing important references and evidence for the control of other vector-borne diseases.
Although this study provides valuable insight into the resistance of Ae. aegypti, certain limitations must be taken into consideration.Firstly, the VGSC gene comprises multiple loci, and only the three most common ones were studied in this research.The contribution of the remaining loci to resistance needs to be explored further.Additionally, other mechanisms related to resistance, such as enzymatic metabolism, have not yet been systematically studied.Secondly, the study only measured resistance to those insecticides commonly used for dengue vectors.In contrast, the potential effects of interactions and resistance levels of other insecticide classes, particularly agricultural insecticides, were not considered.Moreover, due to the rapid expansion of this invasive mosquito species, it is essential to conduct vector resistance surveillance (including bioassay and genetic assay) in other nearby regions, including both endemic and non-endemic countries, to provide comprehensive and systematic evaluation indicators for the prevention and control of DF in Yunnan Province.

Conclusions
In conclusion, this study provides substantial preliminary data on the resistance of Ae. aegypti to insecticides commonly used for dengue vector control.The co-occurrence of the kdr mutations resulted in a higher resistance level and conferred a complex resistance pattern.The development of resistance to synthetic pyrethroid insecticides in the Ae.aegypti populations from the two endemic areas studied here highlight the need for continuous monitoring of insecticide resistance and the rational selection of insecticides.Further research into the mechanism of insecticide resistance in vector mosquitoes is essential for effective decision-making and the development of novel strategies for vector control.

Fig. 1
Fig.1 Aedes aegypti mosquito collection sites from two different urban areas of Yunnan Province, China.Map was created with the software ArcMap 10.1 in ArcGIS 10.1 (ESRI, Redlands, CA, USA)
CI Confidence interval KDT 50 Time (in minutes) taken for 50% of the test mosquitoes to be knocked down, *Statistical significant difference at P < 0.05 a S indicates susceptible, based on a mortality rate of ≥ 98%; R indicates resistant, based on a mortality rate < 90% b -, Not determined c /, No results

Table 2
Correlation between the single kdr mutations S989P, V1016G and F1534C in Ae. aegypti from Ruili and Gengma counties with their resistance phenotypes to the three insecticides

Table 3
Correlation between the kdr mutation types in Ae. aegypti from Ruili and Gengma counties with their resistance phenotypes to the three insecticides