Minghui Zhao1,2
Xin Ran3Dan Xing1Yun Liao2Wei Liu2Yu Bai2Qiang Zhang2Kan Chen2Lan Liu2Mingyu Wu1Zu Ma1
Jian Gao4
Hengduan Zhang1*
Tongyan Zhao1*- 1State Key Laboratory of Pathogen and Biosecurity, Beijing Institute of Microbiology and Epidemiology, Beijing, China
- 2Jiangxi International Travel Healthcare Center, Nanchang, China
- 3Jiangxi Provincial Center for Disease Control and Prevention, Nanchang, China
- 4Jiangsu Provincial Center for Disease Control and Prevention, Nanjing, China
Background: Aedes aegypti and Aedes albopictus are important vectors of human arboviruses, transmitting arboviral diseases such as yellow fever, dengue, chikungunya and Zika. These two mosquitoes coexist on Hainan Island and the Leizhou Peninsula in China. Over the past 40 years, the distribution of Ae. albopictus has gradually expanded in these areas, while the distribution of Ae. aegypti has declined dramatically mainly due to the ecological changes and some other factors such as heavy use of insecticide indoor based on endophagic bloodfeeding of the species.
Methods: This study focused on the knockdown resistance (kdr) genes of both mosquitoes, investigated their mutations, and analyzed their haplotype and evolutionary diversity combined with population genetic features based on the ND4/ND5 genes to further elucidate the molecular mechanisms underlying the development of insecticide resistance in both mosquitoes.
Results: Three mutations, S989P, V1016G and F1534C, were found to be present in Ae. aegypti populations, and the three mutations occurred synergistically. Multiple mutation types (F1534C/S/L/W) of the F1534 locus are found in Ae. albopictus populations, with the three common mutations F1534C, F1534S and F1534L all having multiple independent origins. The F1534W (TTC/TGG) mutation is thought to have evolved from the F1534L (TTC/TTG) mutation. The F1534S (TTC/TCG) mutation has evolved from the F1534S (TTC/TCC) mutation. The most common form of mutation at the F1534 locus found in this study was S1534C, accounting for 20.97%, which may have evolved from the F1534C mutation. In addition, a new non-synonymous mutation M1524I and 28 synonymous mutations were identified in Ae. albopictus populations. Correlation analysis showed that the genetic diversity of Ae. aegypti and Ae. albopictus populations did not correlate with their kdr haplotype diversity (P>0.05), but strong gene flow between populations may have contributed to the evolution of the kdr gene.
Conclusion: The study of kdr gene evolution in the two mosquito species may help to identify the evolutionary trend of insecticide resistance at an early stage and provide a theoretical basis for improving the efficiency of biological vector control and subsequent research into new insecticides.
Introduction
Ae. aegypti and Ae. albopictus are vectors of four major arboviruses, dengue, chikungunya, yellow fever and Zika viruses, which have caused severe health burdens and economic losses worldwide. Many years ago, these vectors and viruses were endemic but confined to specific regions; however, in just over 200 years, they have spread to new tropical, subtropical and temperate regions, expanding their global distribution and likely contributing to coordinated outbreaks of arboviral disease epidemics (Kamgang et al., 2017; Leta et al., 2018; Lwande et al., 2020). In the absence of specific drugs or safe and effective vaccines against these viral infections, the control of arboviral diseases relies heavily on integrated mosquito control measures (Kushwah et al., 2020).
Chemical control of vector insects is one of the major initiatives of integrated control, using five main types of insecticides, including organochlorines, organophosphates, carbamates, pyrethroids and insect growth regulators (Li and Zhong, 2008). Pyrethroids have far outstripped the use of organophosphate and carbamate insecticides with their low toxicity, high efficacy and environmental friendliness. However, the long-term and widespread use of these insecticides has led to the intensification of mosquito resistance (Zhou et al., 2019; Li et al., 2021; Egid et al., 2022). As early as the 1980s, Ae. aegypti and Ae. albopictus in China developed high levels of resistance to insecticides (Wu and Chen, 1988). With frequent iterations and widespread use of insecticides, resistance to insecticides has since developed more rapidly in both species. In addition, Ae. aegypti and Ae. albopictus have different habitat preferences. Ae. aegypti is a domesticated mosquito species that prefers to biting and rest indoors, while Ae. albopictus is a semi-domesticated mosquito species that feeds mostly outdoors, so insecticides or electro-thermal repellents commonly used in homes will have a greater effect on Ae. aegypti than on Ae. albopictus. For example, a survey of insecticide resistance in Ae. aegypti and Ae. albopictus in Yunnan Province from 2021 to 2023 found that both mosquitoes had developed resistance to permethrin, deltamethrin and cyhalothrin, with Ae. aegypti populations in Jinghong have a resistance 825.61 times greater to cyhalothrin. This is also related to the widespread use of insecticide space spraying during the dengue epidemic (Yang et al., 2021; Zheng et al., 2022; Yang et al., 2023). High levels of resistance to pyrethroid insecticides of Ae. albopictus have also been reported in 13 dengue-endemic provinces in China (Zhao et al., 2022). High levels of resistance to pyrethroid insecticides all point to the need for rotation (Hayd et al., 2020).
Knockdown resistance (kdr) is a form of mosquito resistance to pyrethroid insecticides, mainly due to changes in one or more amino acids in the mosquito’s voltage-gated sodium channel (VGSC), which leads to a change in the conformation of the VGSC, reducing the ability of the pyrethroid molecule to bind as efficiently and further avoiding the toxic effects of the insecticide (Dong et al., 2014; Kamgang et al., 2020). Several kdr gene mutations in VGSC structural domains I-IV have been identified in Ae. aegypti and Ae. albopictus. These include the V410L mutation in domain I; the G923V, L982W, S989P, A1007G, I1011M/V and V1016G/I mutations in domain II; the T1520I, I1532T and F1534C/S/L/W/R mutations in domain III; and the D1763Y mutation in domain IV (Auteri et al., 2018; Gan et al., 2021; Akhir et al., 2022). Most of these mutations have been reported to be closely associated with the development of pyrethroid resistance (Kasai et al., 2011; Moyes et al., 2017; Kasai et al., 2019; Rahman et al., 2021b). In particular, the F1534C mutation form was reported to have first appeared in the Aedes mosquito to adapt to the selective pressure of pyrethroid insecticides, and also to provide a platform for the subsequent generation of other mutations, which in turn produced a broader spectrum of pyrethroid resistance (Chen et al., 2019). In addition, the synergistic occurrence of multiple mutations was found to significantly increase resistance in mosquitoes (Smith et al., 2016; Chen et al., 2020; Djiappi-Tchamen et al., 2021). For example, the V410L+V1016I+F1534C mutation in Ae. aegypti significantly increased resistance to pyrethroids (Saavedra-Rodriguez et al., 2018). The synergistic mutation V410L+F1534C was significantly more resistant to permethrin than the V410L mutation alone, and the synergistic mutation S989P+V1016G+F1534C was significantly more resistant to deltamethrin and permethrin than the S989P mutation alone (Hirata et al., 2014). S989P + V1016G and V1016G + F1534C synergistic mutations were also found in Ae. aegypti populations in China (Li et al., 2015; Lan et al., 2020), and have confirmed that the S989P + V1016G mutation is closely associated with pyrethroid resistance (Chung et al., 2019).
The increasing number of reports of resistance to these two mosquito insecticides worldwide has prompted a major effort by the scientific community to identify the key resistance mechanisms and respond to potential public health emergencies (Auteri et al., 2018). For example, the F1534C mutation in Ae. aegypti has been shown to have at least two independent origins, suggesting that selection pressure from DDT and pyrethroid insecticides is responsible for the origin of the F1534C mutation (Cosme et al., 2020). The V410L+F1534C allele mutation may have been caused by the addition of the V410L mutation to the F1534C allele or by a crossover event; the V410L+V1016I+F1534C allele mutation was formed by one or two mutations or recombination steps in the context of the F1534C mutation (Saavedra-Rodriguez et al., 2018; Fan et al., 2020). The S989P+V1016G+F1534C mutant form may be a further mutation type of S989P+V1016G in addition to the V1016G mutation, culminating in a combined form of three mutations (Scott, 2019). It has also been suggested that population migration and hybridization may play an important role in the transmission and maintenance of insecticide resistance alleles in Ae. aegypti and Ae. albopictus, influencing the generation and further evolution of mutations in the kdr gene (Chen, 2019; Rahman et al., 2021a). Studying mutations closely related to kdr, analyzing their mutation frequency and haplotyping and the factors influencing their evolution can provide a comprehensive understanding of the distribution and evolutionary pattern of kdr mutations, help to identify the development trend of insecticide resistance at an early stage and provide a theoretical basis for subsequent research on new insecticides (Cosme et al., 2020).
Hainan Island and the Leizhou Peninsula are located at the southernmost tip of mainland China and have a tropical (subtropical) climate with high temperatures and rainfall throughout the year. The natural conditions of high temperature and humidity are conducive to the breeding of vector mosquito including of Ae. aegypti and Ae. albopictus, and which is an important reason why tropical and subtropical regions are prone to insect-borne infectious diseases (Dong and Tan, 2005). There have been several outbreaks of dengue in these areas, and large-scale spraying during outbreaks is bound to cause the further development of mosquito resistance in the area (Wang et al., 2019; Li et al., 2021). In 1988, these two mosquito species were found not to have developed resistance to pyrethroid insecticides (Wu and Chen, 1988). However, in 1994, studies reported that, Ae. aegypti and Ae. albopictus on Hainan Island had developed resistance to both permethrin and deltamethrin (Li et al., 1994). In 2005, the resistance of these two Aedes species to deltamethrin and cyhalothrin was measured again, and the results showed that Ae. albopictus was less resistant to deltamethrin and cyhalothrin, while Ae. aegypti was susceptible to these two pyrethroids (Zeng et al., 2005). In 2010, a resistance survey of Ae. aegypti on Hainan Island and the Leizhou Peninsula conducted showed all the Ae. aegypti populations were resistant to deltamethrin, permethrin and cypermethrin (Meng et al., 2015). In 2015, studies reported high levels of resistance to both cypermethrin and cyhalothrin in Ae. aegypti in Hainan Island and the Leizhou Peninsula (Li et al., 2015). In recent years, resistance testing of Ae. albopictus populations in these areas has revealed high levels of resistance to pyrethroid insecticides, and mutations such as F1534S, F1534C and F1534 L associated with resistance development have been identified (Chen et al., 2021; Li et al., 2021). However, studies on Ae. aegypti resistance and kdr gene mutation types in these areas have not been reported in recent years, but pyrethroid resistance in Ae. aegypti is a global problem (Hirata et al., 2014). In particular, the evolutionary characteristics and factors influencing the kdr genes of these two mosquito species in these areas have not been studied. Therefore, in this study, it was investigated the mutations and evolutionary mechanisms of VGSC domain IIS6 and IIIS6 fragments, which are closely associated with the development of pyrethroid resistance in Ae. aegypti and Ae. albopictus, important vector mosquitoes in Hainan Island and Leizhou Peninsula, China. The study aims to provide theoretical guidance to improve the efficiency of vector control and to guide the scientific use of insecticides and the development of new insecticides in the region.
Materials and methods
Sample collection and DNA extraction
From June to October 2021, the group visited Hainan Island and the Leizhou Peninsula to collect samples of Ae. aegypti and Ae. albopictus. In order to collect more samples of both mosquito species, areas where the presence of Ae. aegypti had been reported were selected for extensive sampling. Five geographic strains of Ae. aegypti and twenty-two geographic populations of Ae. albopictus were collected. Of these, five geographical populations were sampled for Ae. albopictus sympatric with Ae. aegypti (Table 1). Aedes larvae were collected from various types of water containers such as buckets, tanks, foam boxes, plastic containers, etc. and reared to adults. Post-fledging mosquitoes were morphologically identified and individually stored in anhydrous ethanol for subsequent nucleic acid extraction. Nucleic acid extraction was performed using the DNeasy® Blood and Tissue Kit (Kadjer, Germany) according to the instructions and the extracted DNA was stored at -20°C.
Table 1 Sampling locations and number of Aedes mosquitoes.
Kdr gene amplification
Mutations closely related to pyrethroid resistance have been found in both II and III of the VGSC domain in various Ae. aegypti populations worldwide, with the majority of mutations in Ae. albopictus populations occurring in domain III (Auteri et al., 2018; Gan et al., 2021). Therefore, in this study, the domain IIS6 and IIIS6 fragments in the VGSC of Ae. aegypti and the domain IIIS6 fragment of Ae. albopictus were amplified to determine the mutations in the kdr gene of both mosquitoes. The primers for amplification were as follows: (1) The primers for the 585 bp IIS6 fragment of Ae. aegypti were named AII-F:5’GGTGGAACTTCACCGACTTC3’ and AII-R:5’ GGACGCAATCTGGCTTGTTA3’ (Hamid et al., 2017); (2) The 395bp IIIS6 fragments were named AIII-F:5’ GTGGGAAAGCAGCCGATTCGC3’ and AIII-R:5’ TGTTGAACCCGATGAACAAC3’ (Cosme et al., 2020); (3) The primers for the 780 bp IIIS6 fragment of Ae. albopictus were named BIII-F:5’ GAGAACTCGCCGATGAACTT3’ and BIII-R:5’ GACGACGAAATCGAACAGGT 3’, sequenced with the following primer 5’ TAGCTTTCAGCGGCTTCTTC 3’ (Kasai et al., 2011). The amplification reaction system 25 ul consisted of 12.5 ul of PCR reaction solution (TaKaRa, Japan), 8.5 ul of ddH2O (TaKaRa, Japan), 1 ul of forward primer and 1 ul of reverse primer (synthesized by Biotech, concentration 10uM), followed by 2 ul of DNA template. The PCR program consisted of one cycle at 94°C for 5 min pre-denaturation, followed by 35 cycles of amplification (94°C, 30 s; 63/60/59°C(AII/AIII/BIII), 60 s; 72°C, 60 s) and a final extension at 72°C for 5 min (BIORAD PCR instrument, USA). Amplification products were subjected to 2% agarose electrophoresis. Positive amplification samples were sent to BGI for sequencing. Samples containing intron insertion/deletion fragments of Ae. aegypti IIS6 were cloned, and positive clones selected for sequencing.
Mitochondrial ND4/ND5 gene amplification
The ND4 genes of Ae. aegypti and ND5 genes of Ae. albopictus were amplified for population genetic diversity and dispersal analysis, and the primers were as follows: (1) The 380 bp ND4 primers were named AE-F:5’ GTDYATTTATGATTRCCTAA3’ and AE-R:5’ CTTCGDCTTCCWADWCGTTC3’ (Gao, 2021); (2) the 400bp ND5 primers were named AL-F:5’ TCCTTAGAATAAAATCCCGC3’ and AL-R:5’ GTTTCTGCTTTAGTTCATTCTTC 3’ (Kamgang et al., 2011). The amplification reaction system 25 ul consisted of 12.5 ul of PCR reaction solution (TaKaRa, Japan), 8.5 ul of ddH2O (TaKaRa, Japan), 1 ul of forward primer and 1 ul of reverse primer (synthesized by Biotech), followed by 2 ul of DNA template. The PCR program consisted of one cycle at 94°C for a 5 min pre-denaturation, followed by 35 cycles of amplification (94°C, 30 s; 48/56°C (AE/AL), 60 s; 72°C, 60 s) and a final extension at 72°C for 5 min (BIORAD PCR instrument, USA). Amplification products were subjected to 2% agarose electrophoresis. Positive amplification samples were sent to BGI for sequencing.
Data analysis
To better compare the evolutionary characteristics and differences in the kdr genes of Ae. aegypti and Ae. albopictus in the sympatric areas, it was divided all samples into two large groups for analysis. First, five geographic populations of Ae. aegypti and five geographic populations of Ae. albopictus that sympatric with them were analyzed as a group; Second, the remaining 17 geographic populations of Ae. albopictus were analyzed as a group. All sequences were reviewed in BioEdit software and mutation sites were identified. Sequences were then collated separately in Mega v7 (Kumar et al., 2016). The haplotype diversity index (haplotype number-H, haplotype diversity-Hd, nucleotide diversity-π, average number of nucleotide differences-k) was calculated for both Aedes species by DnaSP v6.12.03 (Rozas et al., 2003). The TCS network maps were plotted using Popart v1.7 (Clement et al., 2002) and phylogenetic trees were plotted in Mega v7 software using a maximum likelihood model. In addition, molecular variance (AMOVA) and genetic differentiation Fst and Nm gene flow were calculated in ARLEQUIN v 3.1 (Excoffier and Lischer, 2010). Correlation analysis between kdr gene haplotype diversity and genetic diversity was assessed using SPSS Statistics.
Results
Sampling results and species identification
A comprehensive survey of potential breeding sites for Ae. aegypti on Hainan Island and the Leizhou Peninsula were conducted from June to October 2021. The findings were consistent with previous reports that the distribution of Ae. aegypti has been significantly reduced (Xie et al., 2011; Chen et al., 2018; Li et al., 2020). Only five populations of Ae. aegypti sympatric with Ae. albopictus were collected, one from the Leizhou Peninsula and four from Hainan Island. In contrast, Ae. albopictus was collected from all areas surveyed (Table 1). Morphological identification of adult mosquitoes was carried out 3 days after fledging. Some studies have reported no significant differences between male and female mosquitoes in kdr gene mutations, and genetic evolutionary analyses generally need to include samples from both sexes (Artigas et al., 2021; Yang et al., 2021). Therefore, half of the samples from each geographical population of female and male mosquitoes were selected for subsequent molecular experiments.
Mutations in kdr gene
A total of three non-synonymous mutations (S989P, V1016G and F1534C) were detected in the kdr gene in the five Ae. aegypti populations, with the F1534C mutation detected in all populations with 100% mutation frequency, in both heterozygous (F1534C) and homogeneous (C1534C) forms, with the highest mutation frequency in the homogeneous form (86.85%). The S989P and V1016G mutations were synergistic and were only present in the HW and WS populations, with both heterozygous (S989P+V1016G) and homogeneous (P989P+G1016G) forms. Both forms were present in the HW population and only one heterozygous form (S989P+V1016G) was present in the WS population (Figure 1A). These mutant forms also occur in association with the F1534C mutation, which was found in only 18 Aedes aegypti samples, including 13 in the HW population (S989P/P989P+V1016G/G1016G+F1534C) and 5 in the WS population (S989P+V1016G+F1534C). These 18 samples were sequenced after cloning due to the presence of both of the above mutations, with base insertions or deletions at intron positions. A clonogenic sequence was randomly selected for data analysis. In addition, as the F1534 locus was mutated in all Ae. aegypti populations and was the only SNP locus for the IIIS6 fragment. Therefore, subsequent Ae. aegypti kdr genotyping studies analyzed only the IIS6 fragment and focused on the S989P and V1016G mutant loci.
Figure 1 Kdr gene mutation frequency of (A) Sympatric Ae. aegypti populations, (B) Sympatric Ae. albopictus populations and (C) 17 Ae. albopictus populations.
Only one non-synonymous mutation (F1534C/S/L) was detected in the kdr gene of the five Ae. albopictus populations sympatric with Ae. aegypti, but the mutation frequencies were significantly different (Figure 1B). The F1534S mutation form was detected in all five Ae. albopictus populations, whereas the F1534L mutation form was only detected in the BS Ae. albopictus population, which had the most diverse mutation types at the F1534 locus (F1534S, F1534C, S1534S and F1534L). The YGH and WS Ae. albopictus populations had only one mutation from F1534S.
In addition to the non-synonymous mutation (F1534C/S/L/W) at the F1534 locus, a non-synonymous mutation M1524I (ATG/ATT) was also detected in 17 other Ae. albopictus populations. This mutation form was found in the FHJC Ae. albopictus population, which occurred in concert with the F1534S mutation. The most common type of mutation at the F1534 locus was S1534C, where there were two base mutations at the same locus, TTC/TTC and TCC/TGC. It was found in all 14 Ae. albopictus populations except ZJJC, SSXX and HQMT populations, with a total mutation frequency of 20.97% (99/472). Among them, the LPC Ae. albopictus population had the most types of kdr gene mutations, with almost all types of mutation forms. In addition, another F1534S mutation (TTC/TCG) was detected in several populations; another F1534L mutation (TTC/CTC) was detected in the NYSS Ae. albopictus population; and the F1534W mutation (TTC/TTG) was found in three populations (TSSS, LPC and CLZ). In addition, an additional 28 synonymous mutations were detected in the IIIS6 segment of the kdr gene in all Ae. albopictus populations. The P1516 locus had the highest mutation frequency at 14.12%, followed by the L1540 and N1541 loci at 10.30% and 10.47%, respectively, and the mutations at the latter two loci occurred synonymously in all populations except one sample; in addition, the F1528 locus also had a high mutation frequency of 9.97%.
Haplotype diversity of the kdr gene
The haplotype diversity of the kdr gene for the five Ae. aegypti, 5 sympatric Ae. albopictus and the 17 Ae. albopictus populations are displayed in Table 2, which showed that the haplotype diversity of the kdr gene was low in the Ae. aegypti population (Hd=0.000-0.6980, total value 0.2441), but very high in the sympatric Ae. albopictus population (Hd= 0.8529-0.9117, total value 0.9224). In addition, the remaining 17 Ae. albopictus populations had a lower haplotype diversity than the sympatric Ae. albopictus populations (Hd=0.6032-0.9035, total value 0.8726). The kdr gene haplotype TCS network showed that the most frequent haplotype in the Ae. aegypti populations was Hap_1, from which 12 other haplotypes evolved (Figure 2B) (Genebank accession number: OR195130). In contrast, the total number of haplotypes in the sympatric Ae. albopictus populations was higher (34 haplotypes) and the evolutionary relationships were more complex, with the more common haplotypes being Hap_5, Hap_3 and Hap_11 (Figure 2A) (Genebank accession number: OR195129). In addition, the total number of haplotypes in the 17 populations of Ae. albopictus was 56, with the most common haplotype being Hap_1, followed by Hap_2 and Hap_5 (Figure 2C) (Genebank accession number: OR166252).
Table 2 Haplotype diversity of Ae. aegypti and Ae. albopictus populations.
Figure 2 TCS network among haplotypes based on kdr gene for (A) Sympatric Ae. albopictus, (B) Sympatric Ae. aegypti and (C) 17 populations of Ae. albopictus. Each line segment represents a single mutation.
Phylogenetic tree analysis of these haplotypes and identification of the different kdr mutations in the tree (Figure 3) showed that the 13 kdr haplotypes of Ae. Aegypti populations could be divided into two branches, with Hap_1 being the S989S+V1016V wild type with the C1534C homogeneous mutation, while the other 12 haplotypes formed a large branch containing the mutation types S989P+V1016G+F1534C (Figure 3A). This also indicated that the S989P+V1016G mutation evolved from the wild-type mutation, which was consistent with the TCS network (Figure 2B). The 34 haplotypes of 5 sympatric Ae. albopictus populations and the 56 haplotypes of the 17 Ae. albopictus populations formed a similar evolutionary pattern, with the wild-type F1534F distributed in all genetic branches (Figures 3B, C). Figure 3B showed the phylogenetic patterns of five sympatric Ae. albopictus populations, from which it could be seen that F1534C, F1534S and F1534L mutations all evolved from the wild type. For example, the F1534S mutant haplotypes Hap_23, Hap_5, Hap_30, Hap_15, Hap_25 and Hap_18 all evolved from the wild haplotype Hap_4; the F1534C mutant haplotypes Hap_3 and Hap_34 all evolved from the wild haplotype Hap_11; the F1534L mutant haplotype Hap_26 evolved from wild haplotype Hap_6, and these results were consistent with the TCS network (Figure 2A). Figure 3C showed the phylogenetic pattern of the 17 Ae. albopictus populations. Similar to Figure 3B, the mutation haplotypes all evolved from the wild haplotype, and in addition, several noncommon mutation types could also be seen in the evolutionary direction from the phylogenetic tree and TCS network. For example, the F1534L mutation haplotype (TTC/CTC) Hap_8 evolved from the wild haplotype Hap_43; the F1534S mutation haplotype (TTC/TCG) Hap_24 evolved from the F1534S mutant haplotype (TTC/TCC) Hap_5; F1534W, the only mutation haplotype (TTC/TGG) Hap_3, evolved from the F1534L mutation haplotype (TTC/TTG) Hap_9. In addition, the most common mutation forms S1534C (TCC/TGC) found in this study were all clustered with the F1534C mutation according to the results of phylogenetic tree. For example, S1534C mutation haplotypes Hap_1 and Hap_42 clustered with F1534C mutant haplotypes Hap_31, Hap_32, Hap_33, Hap _36 and Hap_44 clustered together; S1534C mutation haplotype Hap_2 clustered with F1534C mutant haplotypes Hap_34 and Hap_30. In addition, the newly identified M1524I mutation belongs to the Hap_56 haplotype, which could be seen to have evolved from Hap_52 with one base reversal, consistent with the results of the TCS network diagram analysis (Figure 2C).
Figure 3 Maximum likelihood tree of the kdr haplotypes in (A) Sympatric Ae. aegypti, (B) Sympatric Ae. albopictus and (C) 17 Ae. albopictus populations. ML bootstrap percentages are noted for clades 90%. The different kdr mutations are indicated by different symbols.
Population genetic diversity based on ND4/ND5 genes
Table 2 shows the haplotype diversity for the ND4 gene of Ae. aegypti and the ND5 gene of Ae. albopictus. A relatively low haplotype diversity was found both in Ae. aegypti (Hd=0.37011-0.54924) and sympatric Ae. albopictus populations (Hd=0.00000-0.51282). Four haplotypes and five haplotypes were detected from 140 Ae. aegypti and 130 Ae. albopictus individuals, respectively (Figures 4A, B). The dominant haplotypes of Ae. aegypti were Hap_4 and Hap_3, with the other two haplotypes evolving from Hap_3. The dominant haplotype in Ae. albopictus was Hap_1, with the other four haplotypes evolving from Hap_1. There was only one haplotype in the HT Ae. albopictus population and no intra-population variation. The other 17 Ae. albopictus populations had a higher haplotype diversity (Hd=0.07692-0.72874) and 19 haplotypes were detected, with Hap_1 and Hap_4 being the dominant haplotypes and the other haplotypes all evolved from these two haplotypes (Figure 4C).
Figure 4 TCS network among haplotypes based on ND4 gene for (A) Sympatric Ae. aegypti, ND5 gene for (B) Sympatric Ae. albopictus and (C) 17 populations of Ae. albopictus. Each line segment represents a single mutation.
The AMOVA results indicated that the main variation across populations (Table S1) was consistent with other studies (Shi et al., 2017; Zhang et al., 2017; Wei et al., 2019; Lv et al., 2020). The results of the interpopulation paired Fst analysis are presented in Tables S2 and S3. Table S2 showed that the YGH Ae. aegypti population was more differentiated from the HT, BS and WS populations with less gene flow, while the HW and WS Ae. aegypti populations are less differentiated from the other populations with more gene flow. The WS population of sympatric Ae. albopictus showed significant genetic differentiation from the other populations with low gene flow. While other populations of Ae. albopictus were less differentiated from each other, with high gene flow. In addition, population differentiation varied among the 17 populations of Ae. albopictus, with most populations showing significant genetic differentiation (Fst > 0.05), accounting for 72.26% (101/136) (Guo et al., 2016). However, there were also populations with very low genetic differentiation (Fst < 0) and strong gene flow, accounting for 12.5% (17/136).
To verify whether population genetic diversity affected the evolution of the kdr gene, Pearson’s correlation analyses of the genetic diversity and kdr gene haplotype diversity was conducted, and found that there was no correlation between these two indices among the Ae. aegypti populations (r=0.653, P=0.232>0.05). The same results were also found in populations of sympatric Ae. albopictus (r=0.039, P=0.951>0.05) and the 17 Ae. albopictus (r=-0.003, P=0.990>0.05).
Discussion
The dramatic decline of Ae. aegypti over the past 40 years was closely linked to vector control efforts in Hainan Island and the Leizhou Peninsula, and the heavy use of insecticides had inevitably led to a reduction in mosquito populations and the development of insecticide resistance (Wang et al., 2015; Li et al., 2021). In the present study, five populations of Ae. aegypti and 22 populations of Ae. albopictus were collected from different habitats on Hainan Island and Leizhou Peninsula in 2021, focusing on Ae. aegypti and Ae. albopictus populations in these regions and analyses the mutations and genotyping characteristics of their resistance kdr genes, combined with population genetic and dispersal characteristics based on ND4/ND5 genes, to further reveal the differences in the evolutionary mechanisms and factors influencing the development of resistance in the two mosquito species. Detecting specific pyrethroid resistance mutations could help track and map the spread of mosquito resistance, as well as assess the response of mosquito populations to future insecticide interventions to improve the efficiency of mosquito control in Hainan Island and the Leizhou Peninsula (Ishak et al., 2015).
In the present study, only one mutation type at the F1534 locus was found in various populations of Ae. aegypti, but its mutation frequency reached 100%, with a frequency of 86.85% for the homozygous form (C1534C). On the other hand, the other two mutations (S989P and V1016G) accompanied the heterozygous form of the F1534 locus mutation (F1534C); no cases of S989P or V1016G alone were found, and the mutation frequency was low (13.14%), which was consistent with the report by Cosme et al. (Cosme et al., 2020). No S989P and V1016G mutations were found in the Ae. aegypti kdr gene on the Leizhou Peninsula in 2015, but the F1534C mutation frequency was already 100% (Li et al., 2015). In the present study, it was found that Ae. aegypti populations in the Leizhou Peninsula had evolved into S989P and V1016G mutation types, in addition to the 100% F1534C mutation frequency, which may be related to the frequent use of insecticides in these areas in recent years. In addition, the mutation frequency of F1534C in Ae. aegypti populations in Hainan Island reported in 2015 was only 87.9% (Li et al., 2015), indicating that after eight years of resistance selection, the F1534C mutation frequency of Ae. aegypti in Hainan Island has increased significantly. Some populations have also evolved the S989P and V1016G mutation types (HW). The S989P and V1016G homozygous mutation types were also found in the HW Ae. aegypti population. In addition, the HW population had the highest kdr haplotype diversity of the Ae. aegypti populations, suggesting that the HW Ae. aegypti population may have undergone greater insecticide selection. The genetic differentiation data also showed that only the HW Ae. aegypti population was less genetically distinct from all four other populations, with a strong gene flow (Table S2) (Guo et al., 2016). The strong gene flow may be responsible for the high diversity of kdr haplotypes, and also suggests that gene flow between populations may have contributed to the evolution of kdr genes (Chen, 2019; Rahman et al., 2021a).
Chung et al. also found a S989P+V1016G mutation pattern in Ae. aegypti populations in Taiwan China, and showed their hypothesized evolutionary pathway and closely association with pyrethroid resistance (Chung et al., 2019). Lan et al. found the V1016G+F1534C mutation pattern in Ae. aegypti populations in Yunnan province China (Lan et al., 2020). Chen et al. conducted a literature review and found 6 combined mutation patterns in Ae. aegypti populations from Lancang-Mekong River Basin, including S989P+V1016G, V1016G+F1534C and S989P+V1016G+ F1534C mutation patterns. S989P+ V1016G+F1534C are mainly found in countries such as Thailand, Myanmar, Cambodia and Laos (Chen and Jiang, 2023), no such mutation pattern has been found in China. It was first found in Ae. aegypti populations from China, but it was difficult to infer the evolutionary pattern of the S989P+ V1016G+F1534C mutation pattern in this population because there have been no previous studies of kdr gene mutations in Ae. aegypti in this region and all S989P+ V1016G+F1534C mutations occurred synergistically. The TCS network and phylogenetic tree of this study could not identify the mechanism by which the S989P and V1016G mutation types arose, and further studies are needed in the future.
In contrast to Ae. aegypti, the sympatric Ae. albopictus populations had a more complex kdr gene evolution, with more mutation types and evolutionary forms in the F1534 locus. In the five Ae. albopictus populations, only the F1534 locus had a non-synonymous mutation, but there were three mutation types, F1534C, F1534S and F1534L, with the latter two mutations first identified in Ae. albopictus populations in Hainan Island in 2016 (Chen et al., 2016). These mutations, as seen in the TCS network and phylogenetic tree, all evolved from the wild type and underwent multiple independent origins (Figures 2A, 3B). Compared to Ae. aegypti, all populations of Ae. albopictus have a higher kdr haplotype diversity index, with a lower genetic differentiation between populations, except for WS Ae. albopictus. In particular, the YGH, TW and HT populations show little genetic differentiation between populations, suggesting strong gene flow between these three populations. The WS population is located on the Leizhou Peninsula and is geographically distant from the other populations, resulting in significant genetic differentiation. To further test whether population genetic diversity affected the evolution of the kdr gene, the correlation analyses of two indices (population genetic diversity and the kdr gene haplotype diversity for Aedes aegypti and Aedes albopictus, respectively) was performed. It was found that neither was correlated. The suggested that the high genetic diversity of these two mosquito species did not affect the evolution of the kdr genes.
The genetic diversity of the 17 geographic Ae. albopictus populations was significantly higher compared to the five sympatric Ae. albopictus populations (Table 2), suggesting that the presence of Ae. aegypti may compete with Ae. albopictus for resources and ecological niches, affecting Ae. albopictus genetic diversity (Juliano and Lounibos, 2005; Lounibos and Juliano, 2018). However, the present study showed that there was no correlation between the level of genetic diversity and the genotype diversity of the kdr gene. In addition to the above three mutation types, F1534S (TTC/TCG), F1534L (TTC/CTC), F1534W (TTC/TGG) and S1534C (TCC/TGC) mutation forms were also found at the F1534 locus. Among them, the F1534S (TTC/TCG) mutation was first identified in the Ae. albopictus population in Haikou, Hainan Island, in 2018 (Zhao et al., 2019). This mutation type was found in several populations in Hainan Island and the Leizhou Peninsula in the present study (Figure 1). Evolutionary analysis revealed that this mutation evolved from the F1534S (TTC/TCC) mutation type. The F1534L (TTC/CTC) mutation was first identified in the Sanya Ae. albopictus population in Hainan in 2021 (Chen et al., 2021). It was found that the mutation occurred in only one population in the Leizhou Peninsula (NYSS) and none in Hainan populations in this study. An evolutionary analysis revealed that this mutation type evolved from the wild type, a similar evolutionary mechanism to the F1534L (TTC/TTG) mutation type. The F1534W mutation was first reported in the Shenzhen population of Ae. albopictus in Guangdong Province in 2021 (Chen et al., 2021), but has not yet been reported in Ae. albopictus populations in Hainan Island and the Leizhou Peninsula. The F1534W mutation type was found in three populations of TSSS, LPC and CLZ in this study (Figure 1). This mutant was found to have evolved from the F1534L (TTC/TTG) mutation based on phylogenetic trees and TCS network analysis, with a further base inversion. The heterozygous S1534C (TCC/TGC) mutation form was first identified in the Ae. albopictus population in Haikou in 2018 (Gao et al., 2018), and it was the most common mutation form in this study. Combined with mutational evolutionary analysis, it was clear that this type of mutation was most closely related to the F1534C mutation and was thought to have arisen from further base inversions based on F1534C. There was a fitness cost to the evolution of insecticide resistance in mosquitoes (Kliot and Ghanim, 2012; Deng et al., 2021); the impact of the complex evolutionary form of the F1534 locus in Ae. albopictus on its own fitness and vector competence is a very interesting area of research and deserves further study.
The LPC Ae. albopictus population had the most diverse mutation types at the F1534 locus. The population was located around county road X820 in Jiyang District, Sanya, which was the only convenient route connecting two major national highways G98 and G224 in Sanya, bridging the urban and rural areas with frequent trade and logistics. The rapid development of trade and logistics, high population mobility and frequent use of insecticides were all factors contributing to the development of resistance in mosquito populations (Sa et al., 2019). The WDGC and ASZX Ae. albopictus populations had the highest haplotype diversity in the kdr gene, suggesting that these two populations of Ae. albopictus have greater variation in the kdr gene. Both the WDGC and ASZX Ae. albopictus populations were in central urban areas; large-scale annual insecticide spraying in urban areas could significantly increase the selection pressure for resistance in Ae. albopictus populations (Xu, 2016; Su, 2019). In addition, there are also highly developed transport and logistics in urban areas. The Fst and Nm values also indicated that the WDGC Ae. albopictus population and the other two populations in the city (BSG and SSXX) have low genetic differentiation and strong gene flow, and the ASZX Ae. albopictus population has high genetic exchange (Nm>1) with the other three populations in the city (TSSS, NYSS and ZJJC) (Table S3). Some studies on the evolution of the kdr gene in Culex quinquefasciatus suggested that the distribution of kdr mutations may be the result of widespread insecticide use and frequent population movements (Zang et al., 2023), and the findings of the present study also supported these conclusions.
It was identified a new non-synonymous M1524I in the FHJC Ae. albopictus population. This mutation type has not been previously reported. It had been reported that any non-synonymous mutation on a protein targeted by a chemical insecticide could reduce the affinity of the insecticide for its target protein and lead to resistance (Zeng, 2020). However, as the mutation was only found in one sample, and resistance data for this population were lacking. Therefore, whether this mutation type is associated with insecticide resistance needs to be further investigated. The mutation type T1520I had been reported in Ae. aegypti populations at the nearby 1520 locus and occurred in concert with F1534C, but its association with pyrethroid insecticide resistance had not been reported (Kushwah et al., 2015).
In addition to the two non-synonymous mutations detected in Ae. albopictus populations, twenty-eight synonymous mutations were also detected, six of which were consistent with studies by Wei and Xu (Xu et al., 2016; Wei et al., 2021). Another study also identified synonymous mutations in VGSC domain III in Ae. albopictus populations (Chatterjee et al., 2018; Su, 2019). The presence of these synonymous mutations has increased the genetic diversity of the kdr gene in mosquitoes, but their impact on the development of resistance needs to be further investigated.
Conclusion
This study analyzed the mutations and haplotypes in the kdr genes of Ae. aegypti and Ae. albopictus in Hainan Island and the Leizhou Peninsula, and the influence of population genetic characteristics on them. Three mutations, S989P, V1016G and F1534C, were found in Ae. aegypti populations, and the three mutations were first found to occur synergistically in China. Multiple mutation types at the F1534 locus (F1534C/S/L/W) were found in Ae. albopictus populations, with the three common mutations F1534C, F1534S and F1534L all having multiple independent origins, and the F1534W mutation likely evolving from the F1534L mutation. The most common F1534 locus mutation found in this study was S1534C, which may have evolved further from the F1534C mutation. In addition, a new non-synonymous mutation M1524I and 28 synonymous mutations were identified in the Ae. albopictus population. A comparative analysis with genetic diversity showed that the genetic diversity of Ae. aegypti and Ae. albopictus populations did not correlate with their kdr gene haplotype diversity, but a strong gene flow between populations could contribute to the evolution of the kdr gene. The study of kdr gene mutations in the two mosquito species can help to identify the development trend of insecticide resistance at an early stage and provide a theoretical basis for improving the efficiency of biological vector control and subsequent research on new insecticides.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Materials, further inquiries can be directed to the corresponding authors.
Author contributions
MZ: Conceptualization, Formal Analysis, Investigation, Methodology, Software, Writing – original draft, Writing – review & editing. XR: Formal Analysis, Investigation, Writing – review & editing. DX: Investigation, Resources, Writing – review & editing. YL: Methodology, Writing – review & editing. WL: Resources, Writing – review & editing. YB: Methodology, Writing – review & editing. QZ: Methodology, Writing – review & editing. KC: Resources, Writing – review & editing. LL: Methodology, Writing – review & editing. MW: Resources, Writing – review & editing. ZM: Software, Writing – review & editing. JG: Software, Writing – review & editing. HZ: Conceptualization, Formal Analysis, Resources, Software, Writing – original draft, Writing – review & editing. TZ: Conceptualization, Formal Analysis, Resources, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was funded by grants from the Infective Diseases Prevention and Cure Project of China (no. 2017ZX10303404).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2023.1265873/full#supplementary-material
References
Akhir, M., Wajidi, M., Lavoue, S., Azzam, G., Jaafar, I. S., Awang, B. N., et al. (2022). Knockdown resistance (kdr) gene of Aedes aegypti in Malaysia with the discovery of a novel regional specific point mutation A1007G. Parasit. Vectors. 15 (1), 122. doi: 10.1186/s13071-022-05192-z
Artigas, P., Reguera-Gomez, M., Valero, M. A., Osca, D., Da, S. P. R., Rosa-Freitas, M. G., et al. (2021). Aedes albopictus diversity and relationships in south-western europe and Brazil by rdna/mtdna and phenotypic analyses: its-2, a useful marker for spread studies. Parasit. Vectors. 14 (1), 333. doi: 10.1186/s13071-021-04829-9
Auteri, M., La Russa, F., Blanda, V., Torina, A. (2018). Insecticide resistance associated with kdr mutations in aedes albopictus: an update on worldwide evidences. BioMed. Res. Int. 2018, 3098575. doi: 10.1155/2018/3098575
Chatterjee, M., Ballav, S., Maji, A. K., Basu, N., Sarkar, B. C., Saha, P. (2018). Polymorphisms in voltage-gated sodium channel gene and susceptibility of aedes albopictus to insecticides in three districts of Northern West Bengal, India. PloS Negl. Trop. Dis. 12 (1), e6192. doi: 10.1371/journal.pntd.0006192
Chen, H. (2019). Detection of genetic structure and knockdown resisitance mutations in aedes albopictus (diptera: culicidae) populations from China (Shanghai: Naval Medical University).
Chen, M., Du, Y., Nomura, Y., Zhorov, B. S., Dong, K. (2020). Chronology of sodium channel mutations associated with pyrethroid resistance in aedes aegypti. Arch. Insect Biochem. Physiol. 104 (2), e21686. doi: 10.1002/arch.21686
Chen, M., Du, Y., Wu, S., Nomura, Y., Zhu, G., Zhorov, B. S., et al. (2019). Molecular evidence of sequential evolution of ddt- and pyrethroid-resistant sodium channel in aedes aegypti. PloS Negl. Trop. Dis. 13 (6), e7432. doi: 10.1371/journal.pntd.0007432
Chen, L., Jiang, J. (2023). A review of knock down resistance in voltage-gated sodium channel of Aedes aegypti in Lancang-Mekong River Basin. J. Pathogen Biol. 18 (3), 354–357. doi: 10.13350/j.cjpb.230322
Chen, H., Li, K., Wang, X., Yang, X., Lin, Y., Cai, F., et al. (2016). First identification of kdr allele f1534s in vgsc gene and its association with resistance to pyrethroid insecticides in aedes albopictus populations from Haikou city, Hainan island, China. Infect. Dis. Poverty. 5, 31. doi: 10.1186/s40249-016-0125-x
Chen, Z., Xing, F., Zhang, L., Deng, H., Zhou, J., Huang, J., et al. (2018). Investigation of distribution of aedes aegypti and Ae. albopictus in Leizhou, Guangdong province. Chin. J. Vector Biol. Control. 29 (06), 46–49. doi: 10.11853/j.issn.1003.8280.2018.06.009
Chen, H., Zhou, Q., Dong, H., Yuan, H., Bai, J., Gao, J., et al. (2021). The pattern of kdr mutations correlated with the temperature in field populations of aedes albopictus in China. Parasit. Vectors. 14 (1), 406. doi: 10.1186/s13071-021-04906-z
Chung, H., Cheng, I., Chen, Y., Lin, C., Tomita, T., Teng, H. (2019). Voltage-gated sodium channel intron polymorphism and four mutations comprise six haplotypes in an aedes aegypti population in Taiwan. PloS Neglect. Trop. D. 13 (3), e7291. doi: 10.1371/journal.pntd.0007291
Clement, M. J., Snell, Q., Walker, P., Posada, D., Crandall, K. A. (2002). Tcs: estimating gene genealogies. IEEE Comput. Society. doi: 10.1109/IPDPS.2002.1016585
Cosme, L. V., Gloria-Soria, A., Caccone, A., Powell, J. R., Martins, A. J. (2020). Evolution of kdr haplotypes in worldwide populations of aedes aegypti: independent origins of the F1534C kdr mutation. PloS Negl. Trop. Dis. 14 (4), e8219. doi: 10.1371/journal.pntd.0008219
Deng, J., Guo, Y., Su, X., Liu, S., Yang, W., Wu, Y., et al. (2021). Impact of deltamethrin-resistance in aedes albopictus on its fitness cost and vector competence. PloS Neglect. Trop. D. 15 (4), e9391. doi: 10.1371/journal.pntd.0009391
Djiappi-Tchamen, B., Nana-Ndjangwo, M. S., Mavridis, K., Talipouo, A., Nchoutpouen, E., Makoudjou, I., et al. (2021). Analyses of insecticide resistance genes in aedes aegypti and aedes albopictus mosquito populations from Cameroon. Genes (Basel) 12 (6), 828. doi: 10.3390/genes12060828
Dong, K., Du, Y., Rinkevich, F., Nomura, Y., Xu, P., Wang, L., et al. (2014). Molecular biology of insect sodium channels and pyrethroid resistance. Insect Biochem. Mol. Biol. 50, 1–17. doi: 10.1016/j.ibmb.2014.03.012
Dong, B., Tan, Y. (2005). Prevention and control of important insect-borne infectious diseases in tropical zone. China Trop. Med. 05 (08), 1718–1721. doi: 10.3969/j.issn.1009-9727.2005.08.068
Egid, B. R., Coulibaly, M., Dadzie, S. K., Kamgang, B., McCall, P. J., Sedda, L., et al. (2022). Review of the ecology and behaviour of aedes aegypti and aedes albopictus in western africa and implications for vector control. Curr. Res. Parasitol. Vector Borne Dis. 2, 100074. doi: 10.1016/j.crpvbd.2021.100074
Excoffier, L., Lischer, H. E. (2010). Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under linux and windows. Mol. Ecol. Resour. 10 (3), 564–567. doi: 10.1111/j.1755-0998.2010.02847.x
Fan, Y., O’Grady, P., Yoshimizu, M., Ponlawat, A., Kaufman, P. E., Scott, J. G. (2020). Evidence for both sequential mutations and recombination in the evolution of kdr alleles in aedes aegypti. PloS Negl. Trop. Dis. 14 (4), e8154. doi: 10.1371/journal.pntd.0008154
Gan, S. J., Leong, Y. Q., Bin, B. M., Wong, S. T., Wong, S. F., Mak, J. W., et al. (2021). Dengue fever and insecticide resistance in aedes mosquitoes in Southeast Asia: a review. Parasit. Vectors. 14 (1), 315. doi: 10.1186/s13071-021-04785-4
Gao, J. (2021). Study on genetic differentiation and infection of zika virus among aedes aegypti and aedes albopictus populations (Beijing: Academy of Military Sciences).
Gao, J., Chen, H., Shi, H., Peng, H., Ma, Y. (2018). Correlation between adult pyrethroid resistance and knockdown resistance (kdr) mutations in aedes albopictus (diptera: culicidae) field populations in China. Infect. Dis. Poverty. 7 (1), 86. doi: 10.1186/s40249-018-0471-y
Guo, S., Ling, F., Wang, J., Wu, Y., Hou, J., Gong, Z. (2016). Genetic polymorphism analysis of cytochrome c oxidase submit i gene in aedes albopictus from Zhejiang province, China. Chin. J. Zoonoses 32 (2), 133–136. doi: 10.39697/j.issn.1002-2694.2016.02.007
Hamid, P. H., Prastowo, J., Widyasari, A., Taubert, A., Hermosilla, C. (2017). Knockdown resistance (kdr) of the voltage-gated sodium channel gene of aedes aegypti population in Denpasar, Bali, Indonesia. Parasit. Vectors. 10 (1), 283. doi: 10.1186/s13071-017-2215-4
Hayd, R., Carrara, L., de Melo, L. J., de Almeida, N., Lima, J., Martins, A. J. (2020). Evaluation of resistance to pyrethroid and organophosphate adulticides and kdr genotyping in aedes aegypti populations from roraima, the northernmost Brazilian state. Parasit. Vectors. 13 (1), 264. doi: 10.1186/s13071-020-04127-w
Hirata, K., Komagata, O., Itokawa, K., Yamamoto, A., Tomita, T., Kasai, S. (2014). 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. 8 (8), e3085. doi: 10.1371/journal.pntd.0003085
Ishak, I. H., Jaal, Z., Ranson, H., Wondji, C. S. (2015). Contrasting patterns of insecticide resistance and knockdown resistance (kdr) in the dengue vectors aedes aegypti and aedes albopictus from Malaysia. Parasit. Vectors. 8, 181. doi: 10.1186/s13071-015-0797-2
Juliano, S. A., Lounibos, L. P. (2005). Ecology of invasive mosquitoes: effects on resident species and on human health. Ecol. Lett. 8 (5), 558–574. doi: 10.1111/j.1461-0248.2005.00755
Kamgang, B., Brengues, C., Fontenille, D., Njiokou, F., Simard, F., Paupy, C. (2011). Genetic structure of the tiger mosquito, aedes albopictus, in Cameroon (central Africa). PloS One 6 (5), e20257. doi: 10.1371/journal.pone.0020257
Kamgang, B., Wilson-Bahun, T. A., Yougang, A. P., Lenga, A., Wondji, C. S. (2020). Contrasting resistance patterns to type i and ii pyrethroids in two major arbovirus vectors aedes aegypti and aedes albopictus in the republic of the Congo, Central Africa. Infect. Dis. Poverty. 9 (1), 23. doi: 10.1186/s40249-020-0637-2
Kamgang, B., Yougang, A. P., Tchoupo, M., Riveron, J. M., Wondji, C. (2017). Temporal distribution and insecticide resistance profile of two major arbovirus vectors aedes aegypti and aedes albopictus in yaounde, the capital city of Cameroon. Parasit. Vectors. 10 (1), 469. doi: 10.1186/s13071-017-2408-x
Kasai, S., Caputo, B., Tsunoda, T., Cuong, T. C., Maekawa, Y., Lam-Phua, S. G., et al. (2019). First detection of a vssc allele v1016g conferring a high level of insecticide resistance in aedes albopictus collected from europe (Italy) and asia (Vietnam), 2016: a new emerging threat to controlling arboviral diseases. Euro Surveill. 24 (5), 1700847. doi: 10.2807/1560-7917.ES.2019.24.5.1700847
Kasai, S., Ng, L. C., Lam-Phua, S. G., Tang, C. S., Itokawa, K., Komagata, O., et al. (2011). First detection of a putative knockdown resistance gene in major mosquito vector, aedes albopictus. Jpn. J. Infect. Dis. 64 (3), 217–221. doi: 10.7883/yoken.64.217
Kliot, A., Ghanim, M. (2012). Fitness costs associated with insecticide resistance. Pest Manage. Sci. 68 (11), 1431–1437. doi: 10.1002/ps.3395
Kumar, S., Stecher, G., Tamura, K. (2016). Mega7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 33 (7), 1870–1874. doi: 10.1093/molbev/msw054
Kushwah, R., Dykes, C., Kapoor, N., Adak, T., Singh, O. (2015). Pyrethroid-resistance and presence of two knockdown resistance (kdr) mutations, F1534C and a novel mutation T1520I, in Indian aedes aegypti. PloS Neglect. Trop. D. 9 (1), e3332. doi: 10.1371/journal.pntd.0003332
Kushwah, R., Kaur, T., Dykes, C. L., Ravi, K. H., Kapoor, N., Singh, O. P. (2020). A new knockdown resistance (kdr) mutation, f1534l, in the voltage-gated sodium channel of aedes aegypti, co-occurring with f1534c, s989p and v1016g. Parasit. Vectors. 13 (1), 327. doi: 10.1186/s13071-020-04201-3
Lan, X., Yang, M., Yang, R., Jiang, J. (2020). Analysis of knockdown resistance gene mutations in pyrethroid insecticide resistant populations of aedes aegypti in Yunnan province. Chin. J. Zoonoses 36 (12), 993–999. doi: 10.3969/j.issn.1002-2694.2020.00.163
Leta, S., Beyene, T. J., De Clercq, E. M., Amenu, K., Kraemer, M., Revie, C. W. (2018). Global risk mapping for major diseases transmitted by aedes aegypti and aedes albopictus. Int. J. Infect. Dis. 67, 25–35. doi: 10.1016/j.ijid.2017.11.026
Li, C. X., Kaufman, P. E., Xue, R. D., Zhao, M. H., Wang, G., Yan, T., et al. (2015). Relationship between insecticide resistance and kdr mutations in the dengue vector aedes aegypti in southern China. Parasit. Vectors. 8, 325. doi: 10.1186/s13071-015-0933-z
Li, M., Li, S., Yang, X., Liang, Q., Lan, X., Zeng, L. (1994). Report on the monitoring of vector mosquito resistance in Hainan province. Hainan Med. J. 04), 206–208.
Li, Z., Zhong, J. (2008). Advance in the research of insecticied resistance in dengue fever vectors. Chin. Bull. Entomol. 45 (06), 857–862. doi: 10.3969/j.issn.0452-8255.2008.06.004
Li, Y., Zhou, G., Zhong, D., Wang, X., Hemming-Schroeder, E., David, R. E., et al. (2021). Widespread multiple insecticide resistance in the major dengue vector aedes albopictus in Hainan province, China. Pest Manage. Sci. 77 (4), 1945–1953. doi: 10.1002/ps.6222
Li, Y., Zhou, G., Zhong, S., Wang, X., Zhong, D., Hemming-Schroeder, E., et al. (2020). Spatial heterogeneity and temporal dynamics of mosquito population density and community structure in Hainan island, China. Parasit. Vectors. 13 (1), 444. doi: 10.1186/s13071-020-04326-5
Lounibos, L. P., Juliano, S. A. (2018). Where vectors collide: the importance of mechanisms shaping the realized niche for modeling ranges of invasive aedes mosquitoes. Biol. Invasions 20 (8), 1913–1929. doi: 10.1007/s10530-018-1674-7
Lv, R., Zhu, C. Q., Wang, C., Ai, L., Lv, H., Zhang, B., et al. (2020). Genetic diversity and population structure of aedes aegypti after massive vector control for dengue fever prevention in Yunnan border areas. Sci. Rep.-UK 10 (1), 12731. doi: 10.1038/s41598-020-69668-7
Lwande, O. W., Obanda, V., Lindstrom, A., Ahlm, C., Evander, M., Naslund, J., et al. (2020). Globe-trotting aedes aegypti and aedes albopictus: risk factors for arbovirus pandemics. Vector Borne Zoonotic Dis. 20 (2), 71–81. doi: 10.1089/vbz.2019.2486
Meng, F., Wang, Y., Feng, L., Liu, Q. (2015). Review on dengue prevention and control and integrated mosquito management in China. Chin. J. Vector Biol. Control. 26 (1), 4–10. doi: 10.11853/j.issn.1003.4692.2015.01.002
Moyes, C. L., Vontas, J., Martins, A. J., Ng, L. C., Koou, S. Y., Dusfour, I., et al. (2017). Contemporary status of insecticide resistance in the major aedes vectors of arboviruses infecting humans. PloS Negl. Trop. Dis. 11 (7), e5625. doi: 10.1371/journal.pntd.0005625
Rahman, R. U., Cosme, L. V., Costa, M. M., Carrara, L., Lima, J. B. P., Martins, A. J. (2021a). Insecticide resistance and genetic structure of aedes aegypti populations from Rio de Janeiro state, Brazil. PloS Neglect. Trop. D. 15 (2), e8492. doi: 10.1371/journal.pntd.0008492
Rahman, R. U., Souza, B., Uddin, I., Carrara, L., Brito, L. P., Costa, M. M., et al. (2021b). Insecticide resistance and underlying targets-site and metabolic mechanisms in aedes aegypti and aedes albopictus from Lahore, Pakistan. Sci. Rep. 11 (1), 4555. doi: 10.1038/s41598-021-83465-w
Rozas, J., Sanchez-DelBarrio, J. C., Messeguer, X., Rozas, R. (2003). Dnasp, dna polymorphism analyses by the coalescent and other methods. Bioinformatics 19 (18), 2496–2497. doi: 10.1093/bioinformatics/btg359
Sa, E., Rodovalho, C. M., Sousa, N., Sa, I., Bellinato, D. F., Dias, L., et al. (2019). Evaluation of insecticide resistance in aedes aegypti populations connected by roads and rivers: the case of tocantins state in Brazil. Mem Inst Oswaldo Cruz 114, e180318. doi: 10.1590/0074-02760180318
Saavedra-Rodriguez, K., Maloof, F. V., Campbell, C. L., Garcia-Rejon, J., Lenhart, A., Penilla, P., et al. (2018). Parallel evolution of vgsc mutations at domains is6, iis6 and iiis6 in pyrethroid resistant aedes aegypti from Mexico. Sci. Rep. 8 (1), 6747. doi: 10.1038/s41598-018-25222-0
Scott, J. G. (2019). Life and death at the voltage-sensitive sodium channel: evolution in response to insecticide use. Annu. Rev. Entomol. 64, 243–257. doi: 10.1146/annurev-ento-011118-112420
Shi, Q., Zhang, H., Wang, G., Guo, X., Xing, D., Dong, Y., et al. (2017). The genetic diversity and population structure of domestic aedes aegypti (diptera: culicidae) in Yunnan province, southwestern China. Parasite Vector 10 (1), 292. doi: 10.1186/s13071-017-2213-6
Smith, L. B., Kasai, S., Scott, J. G. (2016). Pyrethroid resistance in aedes aegypti and aedes albopictus: important mosquito vectors of human diseases. Pestic Biochem. Physiol. 133, 1–12. doi: 10.1016/j.pestbp.2016.03.005
Su, X. (2019). Study on insecticide resistance of aedes albopictus (Guangzhou: Southern Medical University).
Wang, X., Chen, H., Yang, X., Lin, Y., Cai, F., Zhong, W., et al. (2015). Resistance to pyrethroid insecticides and analysis of knockdown resistance (kdr) gene mutations in aedes albopictus from Haikou city. Acadmic J. Second Military Med. University 36 (8), 832–838. doi: 10.3724/sp.j.1008.2015.00832
Wang, Y., Meng, F., Han, Z., Jiang, Z., Su, T., Xu, E., et al. (2019). Study on health risk assessment for adult mosquito resistance to public health pesticides. Chin. J. Vector Biol. Control. 30 (2), 117–122. doi: 10.11853/j.issn.1003.8280.2019.02.001
Wei, Y., Wang, J., Song, Z., He, Y., Zheng, Z., Fan, P., et al. (2019). Patterns of spatial genetic structures in aedes albopictus (diptera: culicidae) populations in China. Parasite Vector 12 (1), 552. doi: 10.1186/s13071-019-3801-4
Wei, Y., Zheng, X., He, S., Xin, X., Zhang, J., Hu, K., et al. (2021). Insecticide susceptibility status and knockdown resistance (kdr) mutation in aedes albopictus in China. Parasit. Vectors. 14 (1), 609. doi: 10.1186/s13071-021-05095-5
Wu, N., Chen, D. (1988). Insecticide resistance of aedes albopictus and aedes aegypti in China. J. Med. Pest Contrl 4 (2), 5–8.
Xie, H., Zhou, H., Yang, Y. (2011). Advances in the research on the primary dengue vector aedes aegypti in China. Chin. J. Vector Biol. Control. 22 (02), 194–197.
Xu, J. (2016). Mechanism study of pyrethroid resistance in aedes albopictus (Guangzhou: Southern Medical University).
Xu, J., Bonizzoni, M., Zhong, D., Zhou, G., Cai, S., Li, Y., et al. (2016). Multi-country survey revealed prevalent and novel f1534s mutation in voltage-gated sodium channel (vgsc) gene in aedes albopictus. PloS Negl. Trop. Dis. 10 (5), e4696. doi: 10.1371/journal.pntd.0004696
Yang, M., Lan, X., Jiang, J. (2021). Comparative study of the biological characteristics of lambda-cyhalothrin resistant aedes aegypti and the susceptible population in the Lancang-Meikong river region. J. pathogen Biol. 16 (10), 1220–1222. doi: 10.13350/j.cjpd.211022
Yang, Y., Lin, Z., Na, A., Zheng, Y. (2023). Surveillance of the important ecological habits and insecticedes resistance to dengue vector aedes mosquito in Gengma county, China-Myanmar border. J. pathogen Biol. 18 (04), 447–450. doi: 10.13350/j.cjpd.230415
Yang, L., Liu, D., Chen, J., Wei, X., Li, T., Qiu, D., et al. (2021). Detection and analysis of the knockdown resistance gene in the field populations of aedes albopictus in Zhongshan. Chin. J. Zoonoses 37 (02), 171–175. doi: 10.3969/j.issn.1002-2694.2021.00.005
Zang, C., Wang, X., Cheng, P., Liu, L., Guo, X., Wang, H., et al. (2023). Evaluation of the evolutionary genetics and population structure of culex pipiens pallens in shandong province, China based on knockdown resistance (kdr) mutations and the mtdna-coi gene. BMC Genomics 24 (1), 145. doi: 10.1186/s12864-023-09243-2
Zeng, D. (2020). Relationship between aedes albopictus density and meteorological factors and its insecitide resistance mechanism in Pingshan district, Shenzhen (Guangzhou: Guangdong Pharmaceutical University).
Zeng, L., Zhao, W., Wang, Z., Yang, X., Li, S. (2005). Determination of sensitivity of aedes albopictus and aedes egypy to pyrethroid indescticides in Hainan province. China Trop. Med. 05 (06), 1396–1399. doi: 10.3969/j.issn.1009-9727.2005.06.125
Zhang, R., Leng, P., Wang, X., Zhang, Z. (2017). Molecular analysis and genetic diversity of aedes albopictus (diptera, culicidae) from China. Mitochondrial DNA Part A. 1–6. doi: 10.1080/24701394.2017.1325481
Zhao, C., Zhou, X., Li, W., Lun, X., Liu, X., Wu, H., et al. (2022). Insecticide resistance surveillance and characteristic analysis of dengue vector aedes albopictus in 13 provinces of China in 2020. Chin. J. Vector Biol. Control. 33 (01), 30–37. doi: 10.11853/j.issn.1003.8280.2022.01.006
Zhao, C., Zhu, C., Kai, W., Liu, G., Liu, Q., Lin, L., et al. (2019). Genotypes of knockdown resistance gene and their distribution in aedes albopictus in Haikou, China, in 2018. Chin. J. Vector Biol. Control. 30 (01), 7–11. doi: 10.11853/j.issn.1003.8280.2019.01.002
Zheng, Y., Yang, C., Yang, M., Jiang, J. (2022). Ecological and insecticide resistance surveillance of dengue vector aedes in Yunnan border region of China. China J. Vector Biol. Control. 33 (01), 38–43. doi: 10.11853/j.issn.1003.8280.2022.01.007
Zhou, X., Yang, C., Liu, N., Li, M., Tong, Y., Zeng, X., et al. (2019). Knockdown resistance (kdr) mutations within seventeen field populations of aedes albopictus from Beijing China: first report of a novel v1016g mutation and evolutionary origins of kdr haplotypes. Parasit. Vectors. 12 (1), 180. doi: 10.1186/s13071-019-3423-x
Keywords: Ae. aegypti, Ae. albopictus, kdr, evolution, population genetics
Citation: Zhao M, Ran X, Xing D, Liao Y, Liu W, Bai Y, Zhang Q, Chen K, Liu L, Wu M, Ma Z, Gao J, Zhang H and Zhao T (2023) Evolution of knockdown resistance (kdr) mutations of Aedes aegypti and Aedes albopictus in Hainan Island and Leizhou Peninsula, China. Front. Cell. Infect. Microbiol. 13:1265873. doi: 10.3389/fcimb.2023.1265873
Received: 24 July 2023; Accepted: 05 September 2023;
Published: 21 September 2023.
Edited by:
Michail Kotsyfakis, Foundation for Research and Technology Hellas (FORTH), GreeceReviewed by:
Andreas Kudom, University of Cape Coast, GhanaAlden Estep, United States Department of Agriculture (USDA), United States
Copyright © 2023 Zhao, Ran, Xing, Liao, Liu, Bai, Zhang, Chen, Liu, Wu, Ma, Gao, Zhang and Zhao. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Hengduan Zhang, 453149452@qq.com; Tongyan Zhao, tongyanzhao@126.com