Abstract
Mosquitoes are globally distributed and important vectors for the transmission of many human diseases. Mosquito control is a difficult task and the cost of preventing mosquito-borne diseases is much lower than that for curing the associated diseases. Thus, chemical control remains the most effective tool for mosquito. Due to the long-term intensive use of insecticides to control mosquito vectors, resistance to most chemical insecticides has been reported. This study aimed to investigate the relationship between insecticide resistance and target site mutation of L1014 kdr and G119 ace alleles in 5 species/species group of mosquitoes (Aedes vexans, Ae. albopictus, Anopheles spp., Culex pipiens complex, and Cx. tritaeniorhynchus) obtained from 6 collection sites. For Anopheles spp., the proportion of mosquitoes with mutated alleles in L1014 was 88.4%, homozygous resistant genotypes were observed in 46.7%, and heterozygous resistant genotypes were observed in 41.8%. For the Cx. pipiens complex and Cx. tritaeniorhynchus species, homozygous resistant genotypes were found in 25.9% and 9.8%, respectively. However, target site mutation of L1014 in the Ae. vexans nipponii and Ae. albopictus species was not observed. Anopheles spp., Cx. pipiens complex, and Cx. tritaeniorhynchus mosquitoes were resistant to deltamethrin and chlorpyriphos, whereas Ae. vexans nipponii and Ae. albopictus were clearly susceptible. We also found a correlation between the resistance phenotype and the presence of the L1014 kdr and G119 ace mutations only in the Anopheles spp. population. In this study, we suggest that insecticide resistance poses a growing threat and resistance management must be integrated into all mosquito control programs.
-
Key words: kdr, ace, insecticide resistance, mutation
INTRODUCTION
Global warming affects the survival and density of poikilothermic animals, including insects. The temperature in Korea is rising rapidly due to climate change [
1]. Climate change is a major factor in mosquito ecology and provides conditions that are more favorable for the prevalence of mosquito-borne diseases along with new urban development or population growth [
2]. In addition, countermeasures for infectious diseases must be prepared to prevent foreign mosquito-borne diseases from entering Korea [
2–
4].
In this sense, vector control of mosquitoes has become a pivotal portion of the global strategy to manage mosquito-associated diseases, especially malaria, Japanese encephalitis, West Nile encephalitis, and dengue fever. In this effort, insecticides, such as chlorpyriphos (
O,
O-diethyl
O-(3,5,6-trichloropyridin-2-yl) phosphorothioate) and pyrethroids, are the most important tool. Chlorpyriphos is a broad-spectrum organophosphate insecticide widely used as a pesticide to control a variety of insects [
4]. Pyrethroids, which are analogous to naturally occurring pyrethrins, are a class of insecticides commonly used for indoor spraying against mosquitoes due to their insecticidal potency and environmental stability as well as safety [
5,
6]. Although massive usage of insecticides was highly effective in controlling insect-borne diseases in the past decades, the widespread development of resistance in insects, especially mosquitoes, is now causing serious issues in many regions [
7–
9].
After prolonged exposure to an insecticide over several generations, mosquito populations may evolve resistance against the insecticide. Since mosquitoes can produce many generations per year, high levels of resistance can evolve very quickly. It has been shown that multiple resistance mechanisms are specific to individual mosquito species [
10,
11]. Overall, insecticide resistance is not only conferred via multiple mechanisms as well as reduced cuticle penetrance, but also mediated through the interaction of regulatory genes and resistance genes. However, it is not evident which and how various genes are involved.
The knockdown resistance
kdr gene is the primary mechanism responsible for resistance to dichlorodiphenyltrichloroethane (DDT) and pyrethroids [
12–
14]. The enzyme synaptic acetylcholinesterase (AChE1) is the target of organophosphorus and carbamate, which are competitive inhibitors of acetylcholine [
15,
16]. Genomic studies commonly focus on single nucleotide polymorphisms because they can be readily identified with short region sequencing. However, it is evident that structural variants, including inversions, duplications, and deletions may be involved in adaptation and speciation within mosquitoes [
17,
18]. In this study, we examined sequence variations in specific genes, including knockdown resistance (
kdr) (L1014 in domain II segment 6) and acetylcholine esterase (
ace) (G119 region) collected from mosquitoes at various sampling locations comprising
Ae. vexans nipponii,
Ae. albopictus,
Anopheles spp.,
Cx. pipiens complex, and
Cx. tritaeniorhynchus.
MATERIALS AND METHODS
Mosquito collection
In general, mosquitoes have different insecticide resistances depending on their habitat. Thus, mosquito collection sites for this study were considered depending on the type of mosquito habitat. Adult mosquitoes were collected using unbaited black-light (BL) traps (SC-2000, Shinyeong Co., Namyangju, Korea), BG-Sentinel
TM (BGS) traps (Biogents, Regenburg, Germany), or hand-held battery-powered aspirators at each sample collection site in summer 2017 and 2018. The mosquitoes were identified and classified according to morphology-based specialized taxonomic keys [
19–
22] and compared with standard specimens in the collection of the School of Medicine, Inha University, and the Korea Centers for Diseases Control and Prevention (KCDC). Collections were made over a minimum of 2 days at each site. The collected mosquitoes were pooled by species/species group, collection site, and time of survey, and quickly transferred to microcentrifuge tubes at 4°C where the species/species group were identified. The mosquitoes were assayed to examine target site sequence variations in specific genes, including knockdown resistance (
kdr, L1014 in domain II) and acetylcholine esterase (
ace, G119). For insecticide resistance bottle assays, living mosquitoes were transferred to the School of Medicine, Inha University, as quickly as possible and assayed according to the standard method described below in the section “insecticide resistance bottle assays”.
Table 1 lists the geographic locations, collection methods, and numbers of the collected species/species group of mosquitoes. A total of 7,720 adult mosquitoes comprising the 5 species/species group of interest (
Ae. vexans nipponii (n=1,708),
Ae. albopictus (n=955),
Anopheles spp. (n=2,744),
Cx. pipiens complex (n=1,685), and
Cx. tritaeniorhynchus (n=628), which were identified and classified according to morphology-based specialized taxonomic keys [
19,
20], were collected in summer 2017 and 2018 and the target site mutations were evaluated. L1014F
kdr and G119S
ace mutations were analyzed in 1,010 and 1,258 mosquitoes, respectively.
Adult mosquitoes were collected, were kept in cooled cages to increase survival and provided with 10% sucrose solution ad libitum until the time of testing. The bottle bioassays followed the KCDC and World Health Organization (WHO) guidelines [
23]. In brief, mosquitoes were placed in a 250 ml glass Wheaton bottle coated with 1.25 mg/ml of deltamethrin (Pyrethroid) and chlorpyrifos (Organophosphate) diluted in acetone. The negative control bottle was treated with acetone only. Bottles were swirled so that the glass bottom and inside cap were coated and rotated while rocking so that the side were evenly coated with insecticide. An aspirator was used to gently add 25 mosquitoes into each bottle. Afterward, knockdown mosquitoes were recorded every 30 min for 6 hr. The exposure time for diagnosis resistance of each insecticide was set to 30 min. Collected mosquitoes were transferred into Wheaton bottle coated with each insecticide and incubate for 30 min. After the incubation periods, all living mosquitoes were placed in a pint-sized paper cage with constant access to a 10% sucrose solution. The mortality was assessed after 24 hr. Surviving individuals were frozen at −70°C for molecular analysis.
Mosquito genomic DNA was extracted individually using the QIAamp DNA mini kit (Qiagen GmbH, Hilden, Germany) according to the manufacturer’s protocol. Briefly, the entire body of a single mosquito was homogenized in a microcentrifuge tube containing 80 μl of PBS using a pestle for 30 sec. The homogenized samples were incubated with 20 μl of proteinase K and 100 μl of ATL buffer at 56°C for 2 hr until completely lysed with occasional vortexing during incubation. After centrifuging quickly, 200 μl of AL buffer was added to the mixture followed by incubation at 70°C for 10 min and 200 μl of absolute ethanol was added with vigorous vortexing for thorough mixing. After brief centrifugation, the mixtures were applied to the QIAamp Mini Spin column. The mixture-applied column was washed with 1 ml of AW1 buffer followed by centrifugation at 6,000×g for 1 min and then with 1 ml of AW2 buffer followed by centrifuging at maximum speed for 3 min. The column-bound DNA was eluted with 150 μl of distilled water by centrifugation at 6,000×g for 1 min. The concentration of the extracted DNA was determined by measuring ultraviolet absorbance at 260 nm using a spectrophotometer (NanoDropTM ND-1000, Thermo Fisher Scientific, Lafayette, Colorado, USA), after which the samples were stored at −70°C before use.
Detection of L1014 and G119 mutations
The primer set sequences and reaction conditions for the amplification of
kdr and
ace genes for each targeted mosquito are listed in (
Supplementary Table S1). The specificity of the primers was confirmed using a BLAST search in the GenBank database from NCBI. PCR parameters were as follows: a denaturation step at 94°C for 5 min and then 30 cycles of a denaturation step at 94°C for 1 min, an annealing step at 52°C for 1 min, and extension at 72°C for 1 min, with a final extension at 72°C for 10 min. The PCR products were analyzed by electrophoresis in 1.5% agarose gel. The amplified DNA fragments were cloned into TA cloning vector (Promega Corporation, Madison, Wisconsin, USA) and sequenced by Macrogen Inc. (Seoul, Korea). The DNA sequences were analyzed using DNASTAR and the NCBI-BLAST service.
RESULTS
Genotyping of target site mutations: kdr and ace genes
In this study, we investigated the status of resistance in 5 common mosquitoes (
Ae. vexans nipponii,
Ae. albopictus,
Anopheles spp.,
Cx. pipiens complex, and
Cx. tritaeniorhynchus) collected from different locations in Korea (
Table 1). L1014F
kdr and G119S
ace mutations were analyzed in 1,010 and 1,258 mosquitoes, respectively. The frequency of target site mutations is presented in
Table 2.
kdr target site mutations in 1,010 mosquito samples were analyzed by sequencing of the Domain II region. For the L1014F
kdr target site mutations, homozygous resistant genotypes (L1014F/L1014F and L1014C/L1014C) were observed in 17.5% (177/1,010) and in 0.7% (7/1,010) of 1,010 tested samples, respectively. Homozygous resistant genotypes (L1014C/L1014F and L1014/L1014F) were found in 5.5% (56/1,010) and 5.2% (53/1,010), respectively. A homozygous susceptible genotype (L1014/L1014) was found in 195 (19.3%) mosquitoes from a total of 1,010 sequenced samples. For
Anopheles spp., the rate of mosquitoes with mutated alleles in the L1014
kdr gene was 88.4% (199/225 sequenced samples of
Anopheles spp.), homozygous resistant genotypes were observed in 46.7% (105/225), and heterozygous resistant genotypes (L1014/L1014F and L1014C/L1014F) were observed in 41.8% (94/225). In practice, since it is difficult to classify species group of the genus
Anopheles by morphological differences alone and it is hard to perform the PCR experiments for each
Anopheles subspecies mosquito, we categorized mosquitoes in this genus as simply just
Anopheles spp. It has been reported that there are various
Anopheles subspecies in Korea including
Anopheles lesteri,
An. sinensis,
An. pullus,
An. kleini, and
An. belenrae [
24].
For
Cx. pipiens complex and
Cx. tritaeniorhynchus, homozygous resistant genotypes were found in 25.9% (68/263) and in 9.8% (11/112) of samples, respectively. Heterozygous resistant genotypes were found in 13.4% (15/112) of
Cx. tritaeniorhynchus samples (
Table 2). In addition, various point mutations in the Domain II region were amplified and analyzed in this study, including M976V, L982P, C983R, W986R, N1013D, and L1025M. However, as expected, target site mutation for L1014 in the
kdr Domain II region was not observed in the
Ae. vexans nipponii and
Ae. albopictus species.
For the G119
ace target site mutation, the individuals with the homozygous susceptible genotype (G119) comprised 21% (50/238 sequenced samples in
Anopheles spp.) and those with the homozygous resistant genotype (G119S) comprised 79% (188/238) (
Table 2). However, target site mutations of G119 in the
ace gene in the
Ae. vexans nipponii,
Ae. albopictus,
Cx. pipiens complex, and
Cx. tritaeniorhynchus species were not observed in this study. In addition, various point mutations in the
ace gene were amplified and analyzed, including H103R, S109C, and G118S, especially in the
Culex genus.
Insecticide resistance bottle bioassays showed that the mosquitoes collected from the study areas exhibited 0–61.0% resistance to deltamethrin and chlorpyriphos.
Ae. vexans nipponii and
Ae. albopictus from the study areas showed no resistance to the insecticides tested (pyrethroids (deltamethrin) and organophosphates (chlorpyriphos)), which were used as target insecticides for L1014
kdr and G119
ace genes, respectively (
Table 3). However, resistance was observed in a mosquito species group,
Anopheles spp., with mortality rates of 86.1% (537/624 tested mosquitoes) for deltamethrin and 66.6% (384/577) for chlorpyriphos. Deltamethrin resistance was also observed in the mosquito species
Cx. pipiens complex and
Cx. tritaeniorhynchus, with mortality rates of 78.3% (267/341) and 53.8% (50/93), respectively, but resistance to chlorpyriphos was not observed at all. After insecticide bioassays, the insecticide-resistant mosquitoes exposed to deltamethrin or chlorpyriphos were tested using PCR to identify the target site mutations for L1014
kdr and G119
ace genes (
Table 3). The mosquitoes used for the bioassays with deltamethrin (pyrethroids) (204 specimens) and chlorpyriphos (organophosphate) (316 specimens) were evaluated. The species/species group with the highest frequency of insecticide-resistance was
Anopheles spp., with 78% (68/87 resistant mosquitoes) resistance to deltamethrin, followed by the same species group with 69% (133/193) resistance to chlorpyriphos, and
Cx. pipiens complex and
Cx. tritaeniorhynchus with 14% (10/74, and 6/43, respectively) resistance to deltamethrin. In this sense, the diagnostic test using
kdr and
ace as resistance markers in 5 species/species groups of mosquitoes revealed differences between deltamethrin and chlorpyriphos, with a high specificity for both (78% and 69%, respectively) especially in
Anopheles spp.
Cx. pipiens complex and
Cx. tritaeniorhynchus populations resistant to chlorpyriphos were not observed. These results indicate that the presence of the L1014
kdr and G119
ace gene mutations was significantly associated with resistance to the 2 tested insecticides in
Anopheles spp. However, with the exception of
Anopheles spp. resistance to the 2 insecticides tested, no significant correlation was observed between target site mutations in the L1014
kdr and G119
ace genes and the resistance phenotype in the tested mosquito species/species group.
DISCUSSION
Mosquitoes are important vectors of diseases. The absence of effective vaccines forces the use of insecticides as the major tool to control vector populations [
25,
26]. Mosquitoes were collected from 6 locations nationwide to ensure local diversity of the mosquito species/species group, and the habitat characteristics for each mosquito species/species group were taken into consideration. In general,
Ae. albopictus is considered a forest mosquito, and therefore has had minimal exposure to insecticides sprayed countless times in rice paddies and fields, resulting in low resistance against insecticides. However,
Anopheles spp. mosquitoes, which are found in rice paddies and are therefore exposed to agricultural pesticides sprayed each year, are thought to be more resistant than
Ae. albopictus. Insecticide resistance bottle bioassays showed that the mosquitoes collected from the study areas exhibited resistance to deltamethrin and chlorpyriphos. Especially,
Ae. vexans nipponii and
Ae. albopictus species showed 100% susceptibility to the 2 tested insecticides. However,
Cx. tritaeniorhynchus mosquitoes showed the highest resistance level of 46.2% to deltamethrin and 61.0% to chlorpyriphos.
Anopheles spp. exhibited a resistance level of 13.9% to deltamethrin and 33.4% to chlorpyriphos while the
Cx. pipiens complex exhibited 21.7% resistance to deltamethrin and 21.0% to chlorpyriphos (
Table 3).
Grasping the mechanisms underlying these insecticide resistances is critical to understanding the optimal method to apply these chemicals and conserve their efficacy as tools to control vector populations. In addition, we also examined the well-known and widely tested insecticide target site mutations in L1014 of domain II segment 6 in the
kdr and G119
ace genes. The L1014F substitution mutation confers target site insensitivity and has been reported in
Cx. pipiens mosquitoes [
27]. Other mutations including L1014S and L1014C were also reported in various regions [
28–
31]. Both L1014F
kdr and L1014H/C/S/W
kdr-type resistances have been identified in numerous species [
32]. However, it has been reported that many mosquitoes that present a resistant phenotype also carry the L1014F
kdr allele and the
kdr site mutation may not be a unique mechanism that confers resistance to pyrethroids and DDT [
33]. G119S
ace gene mutation is responsible for the reduction of AChE1 activity in cholinergic synapses and is one of the most common mutations detected in
Cx. pipiens mosquitoes [
34]. An obvious correlation between the degrees of resistance and the frequency of insecticide usage has been well-established in the
Cx. pipiens population [
35]. Our findings indicated that the G119S
ace allele was present only in the tested
Anopheles spp. population, which presented a resistant phenotype to chlorpyriphos. Chlorpyriphos, as well as temephos and pyrimiphos-methyl, are widely used in the control of mosquito larvae [
36].
There are 2 ways to control mosquito-borne diseases; the first is by administering medicine to infected individuals and the second is by controlling vector mosquitoes. Chemical control for mosquito elimination is still the most efficient primary tool. However, due to the large scale and intensive application of insecticides to control mosquito vectors over many years, resistance to most chemical insecticides has been reported. After the US military was stationed in the Southern part of the Korean Peninsula in 1940’s, DDT solution was sprayed in the region to prevent the outbreak of diseases carried by flies and mosquitoes in US military stations [
37]. After this period, extensive use of agricultural pesticides also contributed greatly to the decrease in disease incidence through the control of the mosquito vector of the disease. Due to the nature of their habitat environment,
Ae. vexans nipponii and
Ae. albopictus are thought to have a lower probability of exposure to insecticides than other species. The habitats of these species are usually forests, parks, or quarantine regions, and they generally live out of reach of people and are therefore more susceptible to insecticides. In contrast, other species of mosquitoes examined in this study are continuously exposed to various insecticides. This continuous and intensive exposure to insecticides is responsible for the high levels of resistance. The resistance levels are not only due to intensive prior use but also due to the acquisition of cross-resistance [
36]. Sinegre et al. [
35] reported that in
Cx. pipiens treated with chlorpyriphos, prolonged treatment led to the appearance of resistance to other organophosphates. In addition, it was confirmed that prolonged exposure to organophosphates always leads to the appearance of cross-resistance to other organophosphates and sometimes to certain carbamate-family insecticides [
38]. In addition, it is suggested that the high rates of resistance could be explained by the appearance of cross-resistance, especially when the area was treated only with temephos for larvae and pyrethroid for adults [
36,
39].
Our study demonstrated some correlations between the resistance phenotype and the presence of the L1014
kdr and G119
ace mutations in the
Anopheles spp. population. However, there was no clear correlation between the insecticide resistance phenotype and target site mutation in the other 4 mosquito species tested. Target site mutation coupled with enzyme detoxification has been reported in
An. gambiae and
Cx. quinquefasciatus from Benin [
40]. Moreover, the involvement of metabolic resistance such as carboxylesterase, cytochrome P450 monooxygenase, and glutathione S-transferase requires further investigation to elucidate the mechanisms of resistance to the insecticides currently used in mosquito-vector control.
Insecticide-resistance bottle assays and molecular identification of target-site mutation of L1014 kdr and G119 ace alleles clearly indicated that Anopheles spp., Cx. pipiens complex, and Cx. tritaeniorhynchus mosquitoes were resistant to 2 tested insecticides, deltamethrin and chlorpyriphos, but Ae. vexans nipponii and Ae. albopictus were clearly susceptible. The frequencies of resistant kdr and ace alleles containing L1014F or L1014C and G119S substitution mutations are high for Anopheles spp. populations collected at the monitoring sites. A correlation between the resistance phenotype and the presence of the L1014 kdr and G119 ace mutations was found only in the Anopheles spp. population. Thus, it is clear that insecticide resistance poses a growing threat and resistance management must be integrated into all mosquito control programs.
Notes
-
CONFLICT OF INTEREST
The authors declare that they have no conflicts of interest.
Supplementary Information
ACKNOWLEDGMENT
This work was supported by the Inha University Research Fund and the Korea Centers for Disease Control and Prevention (2017NER44030) in the Republic of Korea.
Table 1Collection sites, methods, and number of mosquitoes collected in summer 2017 and 2018
Table 1
|
Species/species group |
Aedes vexans nipponii
|
Aedes albopictus
|
Anopheles spp. |
Culex pipiens complex |
Culex tritaeniorhynchus
|
|
|
|
|
|
|
Collection site |
Hogok-ri, Hwaseong-si, Gyeonggi-do (Cowshed) |
Arang-dong, Jeju Special Self-Governing Province |
Incheon Grand Park, Incheon Metropolitan City |
Hogok-ri, Hwaseong-si, Gyeonggi-do (Cowshed) |
Hogok-ri, Hwaseong-si, Gyeonggi-do (Cowshed) |
Seochang-dong, Gwangju Metropolitan City |
Gorang-dong, Jeju Special Self-Governing Province |
|
Collection method |
Black-light traps |
Aspirator |
BG-SentinelTM traps |
Aspirator |
Black-light traps |
Aspirator |
Black-light traps |
Aspirator |
Black-light traps |
Aspirator |
|
|
No. of mosquitoes collected |
1,242 |
465 |
743 |
212 |
1,543 |
1,201 |
1,102 |
584 |
417 |
211 |
|
|
Total No. of mosquitoes collected |
1,707 |
955 |
2,744 |
1,686 |
628 |
Table 2Genetic variations in kdr and ace genes, which were found in this study
Table 2
|
Gene locus |
kdr (L1014) |
ace (G119) |
|
|
|
|
Species |
Aedes vexans nipponii
|
Aedes albopictus
|
Anopheles spp. |
Culex pipiens complex |
Culex tritaenior-hynchus
|
Aedes vexans nipponii
|
Aedes albopictus
|
Anopheles spp. |
Culex pipiens complex |
Culex tritaeniorhynchus
|
|
No. of total collected mosquitoes |
1,707 |
955 |
2,744 |
1,686 |
628 |
1,707 |
955 |
2,744 |
1,686 |
628 |
|
|
No. of sequenced samples |
138 |
272 |
225 |
263 |
112 |
134 |
272 |
238 |
441 |
173 |
|
|
No. of mosquitoes carrying mutated alleles (%) |
0 (0) |
0(0) |
199 (88.4) |
68 (25.9) |
26 (23.2) |
0 (0) |
0 (0) |
188 (80.0) |
0 (0) |
0 (0) |
|
|
Other point mutations in same allele or variation of target site (%) |
M976V
L982P
C983R
W986R
W991E
V996E
P1003S
F1005S
A1007T
N1013D |
F959S D960G
R960G R965Q
T969A F971S
M972T S974L
R980Q V999Q
V1010L/MC1011Y
L1025M D1039G
I104555T N1050D
I1063S |
L1014C (3.1)
L1014C/L1014F (24.9)
L1004F (43.6)
L1014/L1014F (16.9)
L1014 (11.5) |
L1014 (74.1)
L1014F (25.9) |
L1014F (9.8)
L1014/L1014F (13.4) |
H103R
S109C
G118S |
H103Y
M1123L
I115F
G117R
G123D
T126A |
G119 (21.0)
G119S (79.0) |
|
|
Table 3Mutation frequency of kdr and ace genes in the mosquito (resistant/susceptible) phenotype grouped using insecticide resistance bioassay
Table 3
|
Insecticide |
Deltamethrin (Pyrethroid, for kdr L1014)a
|
|
|
|
Species |
Aedes vexans nipponii
|
Aedes albopictus
|
Anopheles spp. |
Culex pipiens complex |
Culex tritaenior-hynchus
|
|
Susceptible |
194 |
94 |
537 |
267 |
50 |
|
|
Resistant |
0 |
0 |
87 |
74 |
43 |
|
|
Subtotal |
194 |
94 |
624 |
341 |
93 |
|
|
No. of genetic variation§
|
0 |
0 |
68/87 |
10/74 |
6/43 |
|
|
Frequency of mutated allele |
0 |
0 |
0.78 |
0.14 |
0.14 |
|
|
Chlorpyriphos (Organophosphate, for ace G119)b
|
|
|
Susceptible |
271 |
118 |
384 |
192 |
46 |
|
|
Resistant |
0 |
0 |
193 |
51 |
72 |
|
|
Subtotal |
271 |
118 |
577 |
243 |
118 |
|
|
No. of genetic variationc
|
0 |
0 |
133/193 |
0/51 |
0/72 |
|
|
Frequency of mutated allele |
0 |
0 |
0.69 |
0 |
0 |
References
- 1. Korea Meteorological Administration (KMA). 2019. Abnormal Weather Report [Internet]; [cited 2020 Feb 20]. Available from: http://www.climate.go.kr/home/bbs/view.php?code=93&bname=abnormal&vcode=6385&cpage=1&vNum=Notice&skind=&sword=&category1=&category2
- 2. Yeom JS. Current status and outlook of mosquito-borne diseases in Korea. J Korean Med Assoc 2017;60:468-474. (Korean).
- 3. Bahk YY, Lee HW, Na BK, Kim J, Jin K, Hong YS, Kim TS. Epidemiological characteristics of re-emerging vivax malaria in the Republic of Korea (1993–2017). Korean J Parasitol 2018;56:531-543.
- 4. Korea Centers for Disease Control and Prevention. Epidemiology and Management of Vaccine Preventable Diseases. 5th eds. Cheongju, Korea. Korea Centers for Disease Control and Prevention; 2017.
- 5. Kharkongor M, Hooroo RNK, Dey S. Effects of the insecticide chlorpyrifos, on hatching, mortality and morphology of Duttaphrynus melanostictus embryos. Chemosphere 2018;210:917-921.
- 6. Lu Z, Gan J, Cui X, Delgado-Moreno L, Lin K. Understanding the bioavailability of pyrethroids in the aquatic environment using chemical approaches. Environ Int 2019;129:194-207.
- 7. Hemingway J, Field L, Vontas J. An overview of insecticide resistance. Science 2002;298:96-97.
- 8. Ketterman AJ, Saisawang C, Wongsantichon J. Insect glutathione transferases. Drug Metab Rev 2011;43:253-265.
- 9. Zaim M, Guillet P. Alternative insecticides: an urgent need. Trends Parasitol 2002;18:161-163.
- 10. Liu N. Insecticide resistance in mosquitoes: impact, mechanisms, and research directions. Ann Rev Entomol 2015;60:537-559.
- 11. Main BJ, Everitt A, Cornel AJ, Hormozdiari F, Lanzaro GC. Genetic variation associated with increased insecticide resistance in the malaria mosquito, Anopheles coluzzii
. Parasit Vectors 2018;11:225.
- 12. Dong K. A single amino acid change in the para sodium channel protein is associated with knockdown-resistance (kdr) to pyrethroids insecticides in German cockroach. Insect Biochem Mol Biol 1997;27:93-100.
- 13. Busvine JR. Mechanism of resistance to insecticide in houseflies. Nature 1951;168:193-195.
- 14. Soderlund DM, Knipple DC. The molecular biology of knockdown resistance to pyrethroids insecticides. Insec Biochem Mol Biol 2003;33:563-577.
- 15. Weill M, Lutfalla G, Mogensen K, Chandre F, Berthomieu A, Berticat C, Pasteur N, Philips A, Fort P, Raymond M. Comparative genomics: Insecticide resistance in mosquito vectors. Nature 2003;423:136-137.
- 16. Weill M, Malcolm C, Chandre F, Mogensen K, Berthomieu A, Marquine M, Raymond M. The unique mutation in ace-1 giving high insecticide resistance is easily detectable in mosquito vectors. Insect Mol Biol 2004;13:1-7.
- 17. Lanzaro GC, Touré YT, Carnahan J, Zheng L, Dolo G, Traoré S, Petrarca V, Vernick KD, Taylor CE. Complexities in the genetic structure of Anopheles gambiae populations in west Africa as revealed by microsatellite DNA analysis. Proc Natl Acad Sci USA 1998;95:14260-14265.
- 18. Schmidt JM, Good RT, Appleton B, Sherrard J, Raymant GC, Bogwitz MR, Martin J, Daborn PJ, Goddard ME, Batterham P, Robin C. Copy number variation and transposable elements feature in recent, ongoing adaptation at the Cyp6g1 locus. PLoS Genet 2010;6:e1000998.
- 19. Tanaka K, Mizusawa K, Saugsad ES. A Revision of the Adult and Larval Mosquitoes of Japan (including the Ryukyu Archipelago and the Ogasawara Islands) and Korea (Diptera: Culicidae). Ann Arbor, USA. American Entomological Institute; 1979.
- 20. Lee KW. A Revision of the Illustrated Taxonomic Keys to Genera and Species of Female Mosquitoes of Korea. Seoul, Korea. Department of the Army, 5th Medical Detachment, 168th Medical Battalion, 18th Medical Command; 1998, p 40.
- 21. Jegal S, Jun H, Kim-Jeon MD, Park SH, Ahn SK, Lee J, Gong YW, Joo K, Kwon MJ, Roh JY, Lee WG, Lee W, Bahk YY, Kim TS. Three-year surveillance of culicine mosquiotes (Diptera: Culicidae) for flavivirus infections in Incheon metropolitan city and Hwaseong-si of Gyeonggi-do province, Republic of Korea. Acta Trop 2020;202:105258.
- 22. Bahk YY, Jun H, Park SH, Jung H, Jegal S, Kim-Jeon MD, Roh JY, Lee WG, Ahn SK, Lee J, Joo K, Gong YW, Kwon MJ, Kim TS. Flaviviral disease mosquito vector surveillance in Incheon Metropolitan City and the Hwaseong area, Gyeonggi-do, Republic of Korea. Entomol Res 2020;50:3-13.
- 23. World Health Organization. Test Procedures for Insecticide Resistance Monitoring in Malaria Vector Mosquitoes [Internet]; [cited 2020 Mar 4] Available from: http://www.who.int/malaria/publications/atoz/9789241511575/en/
- 24. Foley DH, Klein TA, Lee IY, Kim MS, Wilkerson RC, Harrison G, Rueeda LM, Kim HC. Mosquito species composition and Plasmodium vivax infection rates on Baengyeong-do (Island), Republic of Korea. Korean J Parasitol 2011;49:313-316.
- 25. Zaim M, Aitio A, Nakashima N. Safety of pyrethroid-treated mosquito nets. Med Vet Entomol 2000;14:1-5.
- 26. World Wildlife Fund. Hazards and Exposures Associated with DDT and Synthetic Pyrethroids used for Vector Control. Toronto, Canada. World Wildlife Fund; 1999.
- 27. Martinez-Torres D, Chevillon C, Brun-Barale A, Bergé JB, Pasteur N, Pauron D. Voltage-dependent Na+ channels in pyrethroid-resistant Culex pipiens L. mosquitoes. Pest Manag Sci 1999;55:1012-1020.
- 28. Xu Q, Tain L, Zhang L, Liu N. Sodium channel genes and their differential genotypes at the F-to-F kdr locus in the mosquito Culex quinquefasciatus
. Biochem Biophys Res Comm 2011;40:645-649.
- 29. Wang ZM, Li CX, Xing D, Yu YH, Liu N, Xue RD, Dong YD, Zhao TY. Detection and widespread distribution of sodium channel alleles characteristic of insecticide resistance in Culex pipiens complex mosquitoes in China. Med Vet Entomol 2012;26:228-232.
- 30. Diabate A, Brengues C, Baldet T, Dabiré KR, Hougard JM, Akogbeto M, Kengne P, Simard F, Guillet P, Hemingway J, Chandre F. The spread of the Leu-Phe kdr mutation through Anopheles gambiae complex in Burkina Faso: genetic introgression and de novo phenomena. Tropical Med Int Health 2004;9:1267-1273.
- 31. Scott JG, Yoshimizu M, Kasal S. Pyrethroid resistance in Culex pipiens mosquitoes. Pest Biochem Physiol 2015;120:68-76.
- 32. Brooke BD.
kdr: can a single mutation produce an entire insecticide resistance phenotype? Trans R Soc Trop Med Hyg 2008;102:524-525.
- 33. Rinkevich FD, Du Y, Dong K. Diversity and convergence of sodium channel mutations involved in resistance to pyrethroids. Pestic Biochem Physiol 2013;106:93-100.
- 34. Weill M, Lutfalla G, Mogensen K, Chandre F, Berthomieu A, Berticat C, Pasteur N, Philips A, Fort P, Raymond M. Comparative genomics: insecticide resistance in mosquito vectors. Nature 2003;423:136-137.
- 35. Sinegre G, Jullian J, Gaven B. Acquisition progressive de la résistance a chlorpyiphos chez les larves de Culex pipiens (L) dans le midi la France. Parasitologia 1977;19:1-2.
- 36. Tmimi FZ, Faraj C, Bkhache M, Mounaji K, Failloux AB, Sarih M. Insecticide resistance and target site mutations (G119S ace-1 and L1014F kdr) of Culex pipiens in Morocco. Parasit Vectors 2018;11:51.
- 37. Chavasse DC, Yop HH. Chemical Methods for the Control of Vectors and Pests of Public Health Importance. 2:Geneva, Switzerland. World Health Organization; 1997, p 129.
- 38. Yeo IS. U.S. military administration’s malaria control activities (1945–1948). Korean J Med Hist 2015;24:35-65.
- 39. Yadouléton A, Badirou K, Agbanrin R, Jöst H, Attolou R, Srinivasan R, Padonou G, Akogbéto M. Insecticide resistance status in Culex quinquefasciatus in Benin. Parasit Vectors 2015;8:17.
- 40. Antonio-Nkondjio C, Poupardin R, Tene BF, Kopya E, Costantini C, Awono-Ambene P, Wondji CS. Investigation of bendiocarb resistance in Anopheles gambiae populations from the city of Yaoundé, Cameroon. Malar J 2016;15:424.
, Byoung-Kuk Na4,5,*, Young Yil Bahk8,*, Tong-Soo Kim1,*
