RESEARCH ARTICLE


Initial Finding of CYP9M10 Gene Haplotype in Culex Quinquefasciatus-Permethrin Resistant Isolate from Natural Population of Central Java Province, Indonesia



Irfanul Chakim*, Rahayu Astuti, Sayono Sayono
Department of Epidemiology, Faculty of Public Health, Universitas Muhammadiyah Semarang, 50273 Semarang, Indonesia


Creative Commons License
© 2018 Chakim et al.

open-access license: This is an open access article distributed under the terms of the Creative Commons Attribution 4.0 International Public License (CC-BY 4.0), a copy of which is available at: https://creativecommons.org/licenses/by/4.0/legalcode. This license permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

* Address correspondence to this author at the Department of Epidemiology, Faculty of Public Health, Universitas Muhammadiyah Semarang, Jl. Kedung Mundu Raya 18, 50273 Semarang, Indonesia; Tel: +6282235657394; E-mail: irfanul.c@student.chula.ac.th


Abstract

Introduction:

Culex is a mosquito genus which widely distributed in the tropical, subtropical and temperate climates, both in urban and rural areas all over the world. Culex mosquitoes, especially Culex quinquefasciatus, are the main vector of filariasis. Thus, public health prevention is mainly directed to control these species and the chemical control is the most widely used approach. However, unfavourable effects of such control to the C. quinquefasciatus population have been known in the form of resistance and emergence of this resistance to various insecticides has been reported in many countries.

Methods:

The metabolic resistance in mosquitoes occurs through the expression of P450 gene family. One of the P450 families is CYP9M10 gene. Partial analysis of the CYP9M10 gene sequences C. quinquefasciatus mosquito reveals the differences between susceptible and resistant alleles. This study was started from the bioassay test of C. quinquefasciatus mosquitoes from five filariasis-endemic areas in Central Java and was continued on the molecular analysis of CYP9M10 gene.

Results:

The results of resistance analysis using bioassay test showed that the C. quinquefasciatus mosquito which has been isolated in Central Java showed high levels of resistance against permethrin 0.75% when mortality rates ranged from 4.8% to 21.6%.

Conclusion:

On the other hand, the molecular analysis revealed three types of haplotypes. The third haplotype had the highest frequency and it exhibited sequence pattern change in the CYP9M10 gene altering the susceptible strain into a resistant one.

Keywords: Culex quinquefasciatus, CYP9M10 gene, Permethrin resistant, Natural population, Central Java, Filariasis.



1. INTRODUCTION

Culex is a mosquito genus which is widespread in tropical, subtropical and temperate climates, both in urban and rural environments around the world [1]. Culex mosquitoes, especially C. pipiens and C. quinquefasciatus, are the main vectors of filariasis in many regions of the world including the Middle East and Eastern Mediterranean countries [2-4]. Filariasis affects more than 25 million people with the disease in the genital area and more than 15 million people with lymphoedema. The distribution of patients is 65% lives in the region of Southeast Asia, 30% is in African region, and the rest is in the other tropical regions [5]. Japanese Encephalitis (JE) is also a disease which is transmitted by Culex mosquitoes and it is spread across South Asia, Southeast Asia, and Pacific [6]. The population at risk is estimated to reach 3 billion people and those live in the endemic countries of JE virus [7]. Culex is also responsible for the transmission of some viral diseases such as Rift Valley Fever [8], the West Nile virus [9-11], Saint Louis Encephalitis (SLE) and Eastern Equine Encephalitis (EEE) [12].

Filariasis is a disease which causes public health problems in Indonesia. According to the cumulative data up to the year 2016, there are 13,009 cases of filariasis clinical cases which spread over 278 of 514 Districts/Municipalities in Indonesia [13]. In Central Java, the number of filariasis cumulative cases by 2012 is 565 patients [14]. In 2013 and 2014, new cases of filariasis were found in 7 districts/municipalities, namely Semarang (11 cases), Demak (3 cases), Grobogan (1 case), the district of Jepara (2 cases), Pekalongan (2 cases), Sukoharjo (1 case), and Wonogiri (1 case) [15].

The distribution of Culex mosquitoes in various regions of the world, including Indonesia and its relationship with the incidence of filariasis has been reported [16]. The survey of vector density showed that Culex is the most dominant mosquitoes, both in Malaysia and Indonesia [17-19]. The density of C. quinquefasciatus mosquitoes reaches 5.91 mosquitoes/person/hour in the house and 4.75 mosquitoes/person/hour in the outdoors [20]. Culex mosquitoes have the highest infectivity against microfilariae compared with other genera [21]. Furthermore, the existence of C. quinquefasciatus mosquito habitat was associated with the incidence of filariasis [22].

Vector control is an important element of the strategy used to control vector-borne disease, and the chemical control is the most widely used approach [23]. The effectiveness of interventions using insecticide has been shown to reduce the incidence of malaria, and it is used in a variety of vector-borne disease control efforts in the epidemiology [24]. Something similar may also occur on the epidemiology of filariasis which is transmitted by Culex sp. The increase of insecticides utilization for vector control in the last decade has been associated with the increase of resistance due to the continuous pressure of the resistance genes. The data about vector resistance to insecticides is still limited and difficult to be found because many countries have not done any adequate routine susceptibility testing [5].

The resistance of Culex mosquito to various insecticides has been reported in many countries. The data from India indicates that C. quinquefasciatus is highly resistant to DDT and Malathion [25]. This mosquito is also resistant to Pyrethroids, Permethrin and Deltamethrin in Zambia [26]. Larvae of C. quinquefasciatus was resistant to Malathion and Permethrin in Kuala Lumpur [27], and it is also resistant to Pyrethroid in Zanzibar [28]. Our previously published paper found that this species also possess a high level of resistance against permethrin insecticide [29].

The mechanism of resistance to insecticides is very complex, including the behavioral change to avoid insecticides, physiological changes, reducing penetration and absorption, target site changes and metabolic resistance. Generally, the resistance in mosquitoes is associated with modification of the target site and metabolic resistance. The changes in the target site resistance include mutations that are resistant to chemical insecticides [30]. On the other hand, the metabolic resistance includes the subtler changes in the expression of a complex array of enzymes and detoxification pathways, the mechanism is still not well understood [31, 32]. Metabolic resistance occurs because of the increase of biological degradation to the insecticide, which is generally the mechanism occurs through overproduction of detoxification enzymes such as P450s, Glutathione S-Transferase (GSTs) and carboxyl/cholinesterase (CCE) [32].

P450s enzyme is expressed by CYP450 gene family. This gene has a role as the resistance marker to most of Pyrethroid Insecticides. Pyrethroid is the most widely known insecticide to be used in vector control. The increased activity level of CYP450 has been observed very often in the Pyrethroid resistance of malaria vector in Africa [33-35]. CYP450 is well-known as the main role in the resistance to Pyrethroid in Culex spp. [36, 37].

In the C. quinquefasciatus mosquitoes, four genes from CYP450 family have emerged as the candidates for resistance to permethrin i.e. CYP6AA7, CYP9J40, CYP9J34 and CYP9M10. The level of overexpression is closely related to the level of resistance in the mosquito laboratory strain of Alabama, USA [37-40]. The overexpression of CYP9M10 has been reported in the Culex-resistant mosquito laboratory strain in Saudi Arabia [41], which also metabolize permethrin and deltamethrin [39]. CYP9M10 is overexpressed in larvae of C. quinquefasciatus permethrin resistant in the laboratory strain [42]. CYP9M10 gene-resistant strain is expressed about 260-fold to pyrethroids compared with the susceptible strain [43].

Besides the data about the level of overexpression of detoxification enzymes, Single Nucleotide Polymorphisms (SNP) which encodes these enzymes have also been observed. Partial analysis of CYP9M10 sequences of C. quinquefasciatus mosquito revealed eight polymorphic locations distinguishing between susceptible alleles and resistant alleles [42]. Yet, there is no information related to the analysis of CYP9M10 genes in the C. quinquefasciatus mosquito population in the field. Thus, this study analyzed the site of CYP9M10 gene in the population of C. quinquefasciatus mosquitoes in the field that are resistant to permethrin.

2. MATERIALS AND METHODS

2.1. Mosquito Specimens and Bioassay Test

The specimens in this study are C. quinquefasciatus mosquitoes identified morphologically and genetically which are originally from 5 filariasis endemic Districts/Municipalities in Central Java in the year 2015 based on our previous report [29, 44]. Those areas are Semarang, Demak, Pekalongan, Jepara, and Grobogan. The adult mosquitoes in each district/municipality have been contacted using the insecticide of permethrin 0.75% for 1 hour with 24 hours holding period. The mosquitoes were distinguished between dead (susceptible) and alive (resistant), then continued for molecular detection of CYP9M10 gene in Eijkman Institute for Molecular Biology, Jakarta.

2.2. Polymerase Chain Reaction (PCR) of CYP9M10 Gene

Genomic DNA of C. quinquefasciatus was extracted using Chelex-100. Then, the DNA was amplified using Forward primers (9M-F): 5'- GAGGCGGAT CCAGTGTTAG -3' and Reverse primer (9M-R): 5'-TCAGTAG CTTCTTTAG GGATTATG-3' [42]. The detail of PCR reaction mixture can be found in our previously published papers [44]. The temperature for PCR conditions was 95ºC (5') pre-denaturation, continued for 40 cycles of 95ºC (30”), 50ºC (30”) and 72ºC (30”) and elongation period of 72ºC (5'). The PCR results were visualized using the electrophoresis method with 2% of agarose gel (2-gram agarose in 100 ml 1X TAE) in 10% ethidium bromide. Subsequently, the PCR products confirmed with gel electrophoresis were sent to bio service unit to be sequenced.

2.3. Sequencing Data Analysis

Improvement of sequencing results in the form of a DNA chromatogram was performed using Finch TV program. Blast identification method was performed at the site ncbi.glm.nih. DNA sequence alignment used Bio edit 7.1.9

Geographic location was determined using eTrex® 10 global positioning system (Garmin, USA) and was analyzed using Arc GIS 10.2

Statistical analysis of ӽ2 and graphic representation of calculation was performed by GraphPad Prism 7.0 program. Haplotype frequency calculations and haplotype visualization used were R (haplo.stats) and TCS 2.1 program, respectively.

3. RESULTS

3.1. The Result of Bioassay

Fig. (1) showed a graphic representation of bioassay result indicating Grobogan was the highest number of died mosquitoes with the lower KDT50 and KDT95. However, all the sites were categorized to be highly resistant to permethrin. A mapping of the spread of the resistant can be seen in Fig. (2). Apparently, due to the most remote area, Grobogan possessed a relatively low level of resistant.

3.2. The Results of the Molecular Analysis of CYP9M10 Gene

Based on the previous report, there are four associated codons that discriminate between susceptible and resistant isolates. These codons are at position of 96, 134, 288 and 390. Our sequences have 99% similarities to reference sequence of SLAB (sensitive) and ISOP450 (resistant). No mutation detected at codon 96, only sensitive and heterozygote at codon 134 and 288 (Fig. 3). However, codon 390 exhibited all types of sensitive homozygote, heterozygote and resistant homozygote. Of each codon, there are 25 samples of mosquitoes, either dead (susceptible) or alive (resistant).

Fig. (1). A graphic representation of the number of dead and tested mosquitoes and KDT50&KDT95 of bioassay test.

Fig. (2). A map exhibits geographic location of the study sites and the phenotypic status of Culex quinquefasciatus towards permethrin 0.75.

Fig. (3). A graphic representation of each indicated alleles of CYP9M10 gene in Culex quenquefasciatus. A ӽ2 test result indicated by *** is considered to be statistically significant < 0.001.

CYP9M10 gene at codon 390 showed polymorphism. The polymorphisms do not change the amino acid (silent mutation). There are three alleles at codon 390 i.e. CAA (wild-type), CAG (mutant allele) and CAA/G (heterozygote allele). Genotype frequencies of CAA, CAG / A, CAG and the allele frequency 390 CAA and 390 CAG respectively are 14%, 23%, 63%, 27% and 73% (Table 1). The tendency of changes occurred in the resistant allele (CAG). The allele frequencies in the resistant group are 56% for the CAG (mutant allele) and 44% for the CAA (wild-type). While on the susceptible group, the percentage increased to 85% for the CAG (mutant allele) and 15% for the CAA (wild-type). Additionally, there are a strong association of ӽ2 test result which showed P-value < 0.001 with a respective odd ratio as 4.706 (Table 1).

Table 1. Allele frequency of codon 96, 134, 288, and 390 of CYP9M10 gene of Culex quinquefasciatus.
CYP9M10
Phenotype
Codon 96 P-value of ӽ2 Codon 134 P-value of ӽ2 Codon 288 P-value of ӽ2 Codon 390 P-value of ӽ2 OR
CAA/CAG
AGT AGA/T AGA ACG ATG/ACG ATG GTC GTC/T GTT CAA CAG/A CAG
Resistant 1,00 0,00 0,00 - 1,00 0,00 0,00 - 1,00 0,00 0,00 - 0,16 0,25 0,59 0.0004 4.706
Sensitive 1,00 0,00 0,00 - 0,90 0,10 0,00 - 0,90 0,10 0,00 - 0,1 0,2 0,7

Haplotype analysis showed that there were three types of haplotypes which were detected based on specific SNP. The calculations of haplotype analysis used R 3.2.1 program with haplo.stats packages. The three of haplotypes are ATTG (resistant strain), TCCA (susceptible strain) and TCCG (strain mixture). The frequencies of each haplotype are ATTG (8%), TCCA (36%), and TCCG (56%) (Table 2). It showed that the strained mixture and susceptible are dominant.

Table 2. Haplotype frequency analysis based on specific SNP.
Haplotype Code Loc-1 Loc-2 Loc-2 Loc-3 Haplotype Frequency
ISOP450/ISOJPAL A T T G 0.08
SLAB T C C A 0.36
C T C C G 0.56

The relationships among each haplotype were analyzed using statistical haplotype kinship networks (Fig. 4). This haplotype kinship network statistics has a function for visualizing the pattern of linkages between each haplotype. Haplotype C which is the most commonly found haplotypes in this study indicates the closeness to the haplotype B / SLAB. While the haplotype A / ISOP450 indicates a distance of 2 small circles, which means that there are genetic distances related to the base substitution in specific SNP.

Fig. (4). Kinship network statistics for haplotype of SLAB, ISO P450 / ISO PAL and C haplotype. The line illustrates the close relationship between haplotypes the probability of > 95%, the each line shows the stages of mutation. The small circle indicates the distance among haplotypes. The relative size of the box or oval shows the haplotypes frequency.

4. DISCUSSION

The results of resistance analysis using a bioassay test denote that C. quinquefasciatus mosquito originally from Central Java has a high level of resistance to permethrin 0.75%. These findings are different from the previous studies in the other countries which found that KDT50 ranged between 40-50 minutes [45] and 33.1 minutes [46]. High levels of resistance to permethrin among C. quinquefasciatus from Central Java, in contrast to some of the isolates from various countries exhibiting a vulnerable response to permethrin 0.75% [45, 47]. It is relatively higher than the isolates from Pathumthani, Nonthaburi and Bangkok of Thailand with a mortality rate of each permethrin with 67.4 ± 2.4%, 72.4 ± 0.4% and 80.6 ± 1.7% [48]. Although in some areas also found high levels of resistance in C. quinquefasciatus resistant to permethrin 0.75% i.e. the isolates of Benin with a range of 4-24% [49] and isolates of Baan Suan, Thailand, with an average mortality of 10.1% [50].

The high number of resistance of C. quinquefasciatus mosquito in Central Java can be the important information related to insecticide-based vector control program. In line with this study, some findings have described the distribution of vector resistance in Indonesia on Aedes and the degree of resistance is increasing [51-53]. In Central Java, a high resistance has also been reported in Aedes with mortality ranged from 1, 6-15, 2% and the KDT50 ranged from 120.2 to 1337.7 minutes [54].

The distribution of insecticide resistance to pyrethroid is various, it depends on the region. The four of the regions are the city of Semarang, Pekalongan district, Demak district and Jepara district, showed a very high level of resistance. This contrasts with the Grobogan district which showed a lower number than the other districts/cities. Geographically, Grobogan district is the most remote region compared to the other regions. However, it still shows a relatively high level of resistance. The phenomenon of the high level of resistance may be caused by an exposure of the local insecticides in dengue vector control program for a long time. Moreover, the breeding places of C. quinquefasciatus mosquitoes may also have been contaminated with insecticides [55].

Insecticides in Central Java have been very widely used with a high frequency. Moreover, the household insecticides and the repellant containing active compounds of pyrethroid have been increasingly used [54]. Interestingly, C. quinquefasciatus species may relatively develop the resistance process to the insecticide faster than other mosquitoes [56]. The close interrelatedness between the resistance and the intensity of insecticide exposure makes clear how the relationship between them is. So, it is necessary to evaluate the use of insecticides in a vector control program. Besides, it is needed for the society to understand the effects of the use of household insecticides.

CYP9M10 is one of gene families encoding P450s enzymes, and it is responsible toward the mechanisms of metabolic resistance to insecticides. The scientific publications about the mechanisms of metabolic resistance related to P450s enzymes are presented in the form of overexpression in the relevant genes, one of them is CYP9M10 [36, 38-43, 57]. The expression of the CYP450 gene is shown by the pure C. quinquefasciatus strain (ISOP450) of mating results between individuals exposed to permethrin over eleven generations [42]. Sequence analysis comparing the laboratory strain (SLAB) and ISOP450 have found five specific codons in exon and three in introns. We replicate a partial analysis using the field strain at four districts and one municipality in Central Java.

The sequence analysis results showed three out of four the analyzed exon sites showed that the isolates from Central Java are still susceptible to molecularly. It is because there are similarities to laboratory sequences (SLAB). However, this result is not in line with the high level of resistance in the bioassay test. It may occur due to a broad resistance mechanism which includes the reduction of penetration and absorption, the modification of target site (i.e. kdr gene mutation) and the complex mechanism of metabolic resistance [58]. Especially in the complex metabolic mechanism, it is not only influenced by a single gene, but also other genes included in the CYP450 gene family. It has a role in the expression of metabolic detoxification [59]. Notwithstanding, at codon 390, all three types of a sensitive homozygote, heterozygote and resistant homozygote were detected with statistically significant and odd ratio of 4.706 which strengthening the fact of haplotype analysis.

CONCLUSION

The study revealed three haplotype patterns between the reference sequence and the isolates from Central Java. Two of them are susceptible and resistant types. But, this study revealed that the third haplotype with the highest frequency which is a mix between susceptible and resistant alleles. This third haplotype may inform the existence of a pattern of sequence changes in the CYP9M10 gene from susceptible strain becomes resistant showing the association between the genotype and phenotype [60]. The data is supported by the variation on the composition of introns.

ETHICS APPROVAL AND CONSENT TO PARTICIPATE

Not applicable.

HUMAN AND ANIMAL RIGHTS

No animals/humans were used for studies that are the basis of this research.

CONSENT FOR PUBLICATION

Not applicable.

CONFLICT OF INTEREST

The authors declare no conflict of interest, financial or otherwise.

ACKNOWLEDGEMENTS

The authors wish to thank people who have consented to take a part in this study, to Eijkman Institute for Molecular Biology in Jakarta; Health Office of Semarang City, Grobogan District, Jepara District, Demak District and Pekalongan District, and those have helped.

REFERENCES

[1] Smith JL, Fonseca DM. Rapid assays for identification of members of the Culex (Culex) pipiens complex, their hybrids, and other sibling species (Diptera: Culicidae). Am J Trop Med Hyg 2004; 70(4): 339-45.
[2] Omar MS. A survey of bancroftian filariasis among South-East Asian expatriate workers in Saudi Arabia. Trop Med Int Health 1996; 1(2): 155-60.
[3] Omar MS, Sheikha AK, Al-Amari OM, Abdalla SE, Musa RA. Field evaluation of two diagnostic antigen tests for Wuchereria bancrofti infection among Indian expatriates in Saudi Arabia. Southeast Asian J Trop Med Public Health 2000; 31(2): 415-8.
[4] Abdel-Hameed AA, Dura WT, Alkhalife IS. An inguinal mass with local vascular lesions induced by a lymphatic filaria. Saudi Med J 2004; 25(8): 1106-8.
[5] World Health Organization. A Global Brief on Vector-Borne Disease 2014.
[6] Solomon T. Control of Japanese encephalitis-within our grasp? N Engl J Med 2006; 355(9): 869-71.
[7] United Nations. World urbanization prospects 2005. Available from: http://www.un.org/esa/population/publications/WUP2005/wup.htm
[8] Balkhy HH, Memish ZA. Rift Valley fever: An uninvited zoonosis in the Arabian peninsula. Int J Antimicrob Agents 2003; 21(2): 153-7.
[9] Paramasivan R, Mishra AC, Mourya DT. West Nile virus: The Indian scenario. Indian J Med Res 2003; 118: 101-8.
[10] Turell MJ, Dohm DJ, Sardelis MR, Oguinn ML, Andreadis TG, Blow JA. An update on the potential of north American mosquitoes (Diptera: Culicidae) to transmit West Nile Virus. J Med Entomol 2005; 42(1): 57-62.
[11] Balenghien T, Vazeille M, Grandadam M, et al. Vector competence of some French Culex and Aedes mosquitoes for West Nile virus. Vector Borne Zoonotic Dis 2008; 8(5): 589-95.
[12] From the Centers for Disease Control and Prevention. Arboviral disease - United States 1994. JAMA 1995; 274(14): 1110-2.
[13] Ministry of Health of the Republic of Indonesia. Profil Kesehatan Indonesia Tahun 2016 2017. Available from: http://www.depkes.go.id/ resources/download/ pusdatin/profil-kesehatan-indonesia/Profil-Kesehatan-Indonesia-2016.pdf
[14] Ministry of Health of the Republic of Indonesia. 2014. Profil Kesehatan Provinsi Jawa Tengah Tahun 2013 [Health Profile of Central Java Province, Year 2013]. Ministry of Health, Jakarta, Republic of Indonesia (in Indonesian). Available from: http://www.depkes.go.id/ resources/download/ profil/PROFIL_KES_PROVINSI_2013/13_Prov_Jateng_2013.pdf
[15] Ministry of Health of the Republic of Indonesia. 2016. Profil Kesehatan Indonesia Tahun 2015 [Indonesia Health Profile, Year 2015]. Ministry of Health, Jakarta, Republic of Indonesia (in Indonesian). Available from: http://www.depkes.go.id/ resources/download/pusdatin/ profil-kesehatan-indonesia/Profil-Kesehatan-Indonesia-2016.pdf
[16] Barr AR. Occurrence and distribution of the Culex pipiens complex. Bull World Health Organ 1967; 37(2): 293-6. Available from: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2554335/pdf/bullwho00241-0123.pdf
[17] Low VL, Chen CD, Lee HL, Lim PE, Leong CS, Sofian-Azirun M. Nationwide distribution of Culex mosquitoes and associated habitat characteristics at residential areas in Malaysia. J Am Mosq Control Assoc 2012; 28(3): 160-9.
[18] Islamiyah M, Leksono SA, Gama PZ. Distribusi dan komposisi nyamuk di wilayah Mojokerto. Jurnal Biotropika 2013; 1: 80-5. [in Indonesian with summary in English].
[19] Juhanudin N, Leksono SA. Distribusi spasial nyamuk diurnal secara ekologi di Kabupaten Lamongan. Jurnal Biotropika 2013; 1(3): 124-8.
[20] Juwita R. 2011.Kepadatan nyamuk Culex quinquefasciatus tersangka vektor filariasis dan gambaran kondisi lingkungan di desa Jeruksari Kecamatan Tirto Kabupaten Pekalongan Available from: http://eprints.undip.ac.id/34975/1/4169.pdf
[21] Michael E, Snow LC, Bockarie MJ. Ecological meta-analysis of density-dependent processes in the transmission of lymphatic filariasis: Survival of infected vectors. J Med Entomol 2009; 46(4): 873-80.
[22] Windiastuti IKA, Suhartono S, Jazuli N. Hubungan kondisi lingkungan rumah, sosial ekonomi, dan perilaku masyarakat dengan kejadian filariasis di Kecamatan Pekalongan Selatan Kota Pekalongan. J Kesehat Lingkung Indones 2013; 12: 51-7.
[23] World Health Organization. 2011. Global Insecticide Use for Vector-borne Disease Control: a 10-year assessment (2000-2009). WHO, Geneva. Ref. number: WHO/HTM/NTD/VEM/WHOPES/2011.6
[24] World Health Organization. 2007.Insecticide-treated mosquito nets: a WHO position statement Available from: http://www.bvsde.paho.org/bvsacd/cd66/ITNspospaperfinal.pdf
[25] Kumar K, Sharma AK, Kumar S, Patel S, Sarkar M, Chauhan LS. Multiple insecticide resistance/susceptibility status of Culex quinquefasciatus, principal vector of bancroftian filariasis from filaria endemic areas of northern India. Asian Pac J Trop Med 2011; 4(6): 426-9.
[26] Norris LC, Norris DE. Insecticide resistance in Culex quinquefasciatus mosquitoes after the introduction of insecticide-treated bed nets in Macha, Zambia. J Vector Ecol 2011; 36(2): 411-20.
[27] Hamdan H, Sofian-Azirun M, Nazni W, Lee HL. Insecticide resistance development in Culex quinquefasciatus (Say), Aedes aegypti (L.) and Aedes albopictus (Skuse) larvae against malathion, permethrin and temephos. Trop Biomed 2005; 22(1): 45-52.
[28] Jones CM, Machin C, Mohammed K, et al. Insecticide resistance in Culex quinquefasciatus from Zanzibar: Implications for vector control programmes. Parasit Vectors 2012; 5: 78.
[29] Chakim I, Sayono S, Astuti R. High levels of resistance in a Culex quinquefasciatus population to the insecticide permethrin in filariasis endemic areas in central java. Makara J Sci 2017; 21: 149-54.
[30] Donnelly MJ, Corbel V, Weetman D, Wilding CS, Williamson MS, Black WC IV. Does kdr genotype predict insecticide-resistance phenotype in mosquitoes? Trends Parasitol 2009; 25(5): 213-9.
[31] Feyereisen R. Insect cytochrome P450. Insect Mol Biol 2005.
[32] Hemingway J, Ranson H. Insecticide resistance in insect vectors of human disease. Annu Rev Entomol 2000; 45: 371-91.
[33] Vulule JM, Beach RF, Atieli FK, et al. Elevated oxidase and esterase levels associated with permethrin tolerance in Anopheles gambiae from Kenyan villages using permethrin-impregnated nets. Med Vet Entomol 1999; 13(3): 239-44.
[34] Matambo TS, Paine MJ, Coetzee M, Koekemoer LL. Sequence characterization of cytochrome P450 CYP6P9 in pyrethroid resistant and susceptible Anopheles funestus (Diptera: Culicidae). Genet Mol Res 2010; 9(1): 554-64.
[35] Djouaka RF, Bakare AA, Coulibaly ON, et al. Expression of the cytochrome P450s, CYP6P3 and CYP6M2 are significantly elevated in multiple pyrethroid resistant populations of Anopheles gambiae s.s. from Southern Benin and Nigeria. BMC Genomics 2008; 9: 538.
[36] Kasai S, Weerashinghe IS, Shono T. P450 monooxygenases are an important mechanism of permethrin resistance in Culex quinquefasciatus Say larvae. Arch Inch Biochein Physiol 1998; 37(1): 47-56.
[37] Kasai S, Shono T, Yamakawa M. Molecular cloning and nucleotide sequence of a cytochrome P450 cDNA from a pyrethroid-resistant mosquito, Culex quinquefasciatus say. Insect Mol Biol 1998; 7(2): 185-90.
[38] Nauen R. Insecticide resistance in disease vectors of public health importance. Pest Manag Sci 2007; 63(7): 628-33.
[39] Wilding CS, Smith I, Lynd A, et al. A cis-regulatory sequence driving metabolic insecticide resistance in mosquitoes: Functional characterisation and signatures of selection. Insect Biochem Mol Biol 2012; 42(9): 699-707.
[40] Liu N, Li T, Reid WR, Yang T, Zhang L. Multiple Cytochrome P450 genes: their constitutive overexpression and permethrin induction in insecticide resistant mosquitoes, Culex quinquefasciatus. PLoS One 2011; 6(8): e23403.
[41] Amin AM, Hemingway J. Preliminary investigation of the mechanisms of DDT and pyrethroid resistance in Culex quinquefasciatus Say (Diptera, Culicidae) from Saudi Arabia. Bull Entomol Res 1989; 79: 361-6.
[42] Hardstone MC, Komagata O, Kasai S, Tomita T, Scott JG. Use of isogenic strains indicates CYP9M10 is linked to permethrin resistance in Culex pipiens quinquefasciatus. Insect Mol Biol 2010; 19(6): 717-26.
[43] Itokawa K, Komagata O, Kasai S, Okamura Y, Masada M, Tomita T. Genomic structures of Cyp9m10 in pyrethroid resistant and susceptible strains of Culex quinquefasciatus. Insect Biochem Mol Biol 2010; 40(9): 631-40.
[44] Chakim I, Tyasningrum WS, Wardani HC, Sayono S. Molecular identification, phylogeography, and genetic diversity of Culex quinquefasciatus in Central Java province, Indonesia. Ann Trop Med Public Health 2017; 10: 341-7.
[45] Wan-Norafikah O, Nazni WA, Lee HL, Zainol-Ariffin P, Sofian-Azirun M. Development of permethrin resistance in Culex quinquefasciatus Say in Kuala Lumpur, Malaysia. Saudi J Biol Sci 2013; 20(3): 241-50.
[46] Komalamisra N, Trongtokit Y, Palakul K, et al. Insecticide susceptibility of mosquitoes invading tsunami-affected areas of Thailand. Southeast Asian J Trop Med Public Health 2006; 37(Suppl. 3): 118-22.
[47] Vatandoost H, Ezeddinloo L, Mahvi AH, Abai MR, Kia EB, Mobedi I. Enhanced tolerance of house mosquito to different insecticides due to agricultural and household pesticides in sewage system of Tehran, Iran. Iran J Environ Health Sci Eng 2004; 1: 42-6.
[48] Thanispong K, Sathantriphop S, Chareonviriyaphap T. Insecticide resistance of Aedes aegypti and Culex quinquefasciatus in Thailand. J Pestic Sci 2008; 33: 351-6.
[49] Yadouléton A, Badirou K, Agbanrin R, et al. Insecticide resistance status in Culex quinquefasciatus in Benin. Parasit Vectors 2015; 8: 17.
[50] Sathantriphop S, Paeporn P, Supaphathom K. Detection of insecticides resistance status in Culex quinquefasciatus and Aedes aegypti to four major groups of insecticides. Trop Biomed 2006; 23(1): 97-101.
[51] Ahmad I, Astari S, Tan M. Resistance of Aedes aegypti (Diptera: Culicidae) in 2006 to pyrethroid insecticides in Indonesia and its association with oxidase and esterase levels. Pak J Biol Sci 2007; 10(20): 3688-92.
[52] Ahmad I, Astari S, Rahayu R, Hariani N. Status Kerentanan Aedes aegypti (Diptera: Culicidae) pada Tahun 2006-2007 terhadap Malathion di Bandung, Jakarta, Surabaya, Palembang dan Palu. Biosfera 2009; 26: 82-9.
[53] Ministry of Health of the Republic of Indonesia. 2011. Profil Kesehatan Indonesia Tahun 2010 [Indonesian Health Profile Year 2010]. Jakarta, Indonesia (in Indonesian).
[54] Sayono S, Hidayati APN, Fahri S, et al. Distribution of voltage-gated sodium channel (Nav) alleles among the aedes aegypti populations in central java province and its association with resistance to pyrethroid insecticides. PLoS One 2016; 11(3): e0150577.
[55] Canyon DV, Hii JLK. Insecticide susceptibility status of Aedes aegypti (Diptera: Culicidae) from Townsville. Aust Entomol 1999; 38: 40-3.
[56] Hamon J, Mouchet J. La résistance aux insecticides chez Culex pipiens fatigans Wiedemann. Bull World Health Organ 1967; 37(2): 277-86. [in French with summary in English].
[57] David JP, Ismail HM, Chandor-Proust A, Paine MJI. Role of cytochrome P450s in insecticide resistance: Impact on the control of mosquito-borne diseases and use of insecticides on Earth. Philos Trans R Soc Lond B Biol Sci 2013; 368(1612): 20120429.
[58] Yang T, Liu N. Genome analysis of cytochrome P450s and their expression profiles in insecticide resistant mosquitoes, Culex quinquefasciatus. PLoS One 2011; 6(12): e29418.
[59] Gong Y, Li T, Zhang L, Gao X, Liu N. Permethrin induction of multiple cytochrome P450 genes in insecticide resistant mosquitoes, Culex quinquefasciatus. Int J Biol Sci 2013; 9(9): 863-71.
[60] Crawford DC, Nickerson DA. Definition and clinical importance of haplotypes. Annu Rev Med 2005; 56: 303-20.