Dengue is the most important mosquito-borne viral disease in the world [1]. Dengue is also a significant public health problem in Indonesia since the first reported dengue outbreak in 1968. All four dengue virus (DENV; Flavivirus) serotypes are detected, suggesting hyperendemicity in all of 34 provinces in Indonesia [2, 3]. Dengue virus is transmitted by A. aegypti, a mosquito that has also been linked with the transmission of other major arboviral diseases including chikungunya (CHIKV) and Zika (ZIKV). The virologically confirmed chikungunya was initially reported during an outbreak in Jambi in 1982 [4]. After a near 20 years of absence, twenty-four distinct outbreaks of probable chikungunya reemerged in South Sumatera, Aceh and West Java in the early 2000s [4]. To date, there has been no recorded fatality related to CHIKV infection in Indonesia [5]. Zika virus (ZIKV, family Flaviviridae) has become one of the global public health threats because of its association with Guillain-Barré syndrome and microcephaly. Until recently, most evidence for Zika virus infection in Asia, including in Indonesia, has been serologic [6].
Indonesia has made dengue a notifiable disease. On average, 136,670 DENV infections and 1,112 deaths were reported annually from 2013 to 2016, an incidence of about 54 dengue fever or dengue hemorrhagic fever cases per 100,000 population and a case fatality rate of approximately 1% [7]. Dengue cases in Indonesia are expected to be under-reported due to poor disease surveillance and a low level of reporting [8]. The surveillance database likely covers dengue probable cases with supportive dengue serology or with epidemiologic linkage, and dengue confirmed cases with confirmatory laboratory criteria [8] [9]. Dengue infections deprived of diagnostic tests may not end up in the surveillance database [8]. Despite the level of uncertainty on total case number, dengue visibly results in a considerable cost to the health sector, and a heavy economic and social impact is likely.
Indonesia has made progress in many areas of dengue prevention and control. On August 2016, Indonesia became the sixth country to approve the first licensed dengue vaccine, Dengvaxia® [10, 11]. To date, the vaccine so far has been approved in 20 countries and launched in 11 including Indonesia [11]. However in December 2017, the National Agency of Drug and Food (Badan POM) suspended the use of dengue vaccine in Indonesia over safety concerns [12] and later in February 2018, the BPOM revised recommendation that the use of vaccine be limited for people 9 – 16 years of age who had a dengue infection prior to vaccination and that the vaccine should not be administered to people who have never had a dengue infection before [13]. There is no rapid, reliable test for previous dengue infection existing, so the Dengvaxia® vaccine cannot be widely used [14]. Nevertheless, effective vector control methods are still essential, targeting A. aegypti in its immature and adult stages. It has been argued that even if the current vaccine is highly targeted and low-cost, sustained mosquito control will remain cost-effective [15].
The mainstay of the current vector control program in Indonesia is environmental management, which in recent years, has emphasized community participation to reduce container breeding sites [2, 16]. The country is renowned for its 3M plus campaigns, aimed at covering and cleaning water containers and burying discarded water containers, complemented with biological approaches using natural predators or pathogens as alternatives for vector control [2, 16]. The efficacy of community-based approaches is measured by a larva free index (percentage of houses free from A. aegypti larvae and pupae infestation); unfortunately, data from the last five years show a larvae free index ranging from 24% - 80%, less than the target of 95% [7].
Chemical insecticides still play a central role in dengue vector control in Indonesia. Vector control during outbreaks depends primarily on chemical insecticides to effect a rapid reduction in the number of infected mosquitoes and to break the dengue transmission cycle [2, 16, 17, 18]. Since the 1970s, the organophosphates malathion and temephos have been widely used to control dengue, and starting in the 1980s, dengue vector control has been highly reliant on pyrethroids. Pyrethroids are also widely used in public health for prevention and control of other mosquito-borne diseases such as malaria and filariasis, and as agricultural insecticides [19]. Pyrethroids are also approved for household protection against dengue [20, 21]. The abundant and prolonged use of pyrethroids has led to the development of resistance in A. aegypti populations in many countries including Indonesia. Resistance to pyrethroids, based on bioassay data, has been reported in some areas of Indonesia [22, 23, 24, 25, 26].
Two main mechanisms for pyrethroid resistance have been identified in A. aegypti, metabolic resistance and target site resistance. Metabolic resistance occurs when increased levels or modified activities of one or more detoxifying enzymes result in a more rapid detoxification of the insecticide, preventing the insecticide from reaching its target in the nervous system. Metabolic resistance involves three groups of enzymes: esterases, multi-function oxidases P450 and glutathione s-transferases (GST) [27, 28, 29]. Limited studies about metabolic resistance based on biochemical assays are available for Indonesian dengue vectors [26, 30, 31, 32].
Pyrethroids act on the insect nervous system, targeting the voltage sensitive sodium channel (VSSC). They modify the gating kinetics of the channel by slowing both the activation and the inactivation, stimulating the nerve cells to produce repetitive discharge that lead to paralysis and death of insects, an effect known as knockdown. Target site resistance is caused by point mutations in the VSSC gene, resulting in amino acid substitutions that affect pyrethroid binding sites, being known as kdr (knockdown resistance) mutations[33, 34]. A total of 12 point mutations in the Vssc of A. aegypti have been identified to be associated with kdr resistance to pyrethroids. Only five of these mutations have been functionally confirmed to reduce the sensitivity of mosquito sodium channels to pyrethroids including S989P, I1011M, V1016G/I, F1534C, and recently V410L [35]The mutations L982W, G923V, I1011M, and V1016G were the first reported sodium channel mutations in permethrin/DDT-resistant populations of A. aegypti from various countries [36], located in the domain IIS5 (for L982W and G923 mutations) and IIS6 (for I1011M and V1016G mutations) of the VSSC [36]. Further studies have reported novel mutations, including I1011V and V1016I in Latin American populations [37]. In Asian countries, the F1534C mutation (in domain IIIS6) was detected in A. aegypti mosquitoes from Thailand [38]and Vietnam [39]. The S989P mutation, located in linker between domain IIS5–S6, was first reported in Ae aegypti populations in Thailand [40], the D1763Y mutation (in linker IVS5-S6) was reported in Taiwan population [41], and mutation T1520I (in domain IIS6) was detected in an Indian population [42]. In 2017, the mutation V410L was firstly identified in Brazilian strains of A. aegypti [43], and in 2018, a novel mutation V419L was found in populations from Colombia [44] (both mutations are located in domain IS6).
Some of these mutations, V1016G, F1534C, and S989P, are widely distributed and detected in pyrethroid-resistant populations in Southeast Asian countries including Thailand, Indonesia, Malaysia, Singapore, Vietnam, Cambodia, and Laos [45]. Co-occurrence of kdr mutations has been a common phenomenon observed and, for some combinations, has been shown to confer a higher level of resistance than singly occurring mutations [46]. In A. aegypti from Southeast Asia, at least three patterns of mutational associations have been identified: V1016G/F1534C, V1016G/S989P, and V1016G/F1534C/S989P. Co-occurrence of V1016G and F1534C was reported in populations from Thailand [47], Myanmar [46], Malaysia [48], and Indonesia [49, 50, 51, 52]. Meanwhile, co-occurrence of V1016G and S989P point mutations was detected in Thailand [47], Myanmar [46], Indonesia, [49, 50, 52], and Papua New Guinea [53]. In addition, co-occurrence of triple mutations V1016G/F1534C/S989P in heterozygous form has been identified commonly in A. aegypti from Thailand [54], Myanmar [46], and Indonesia [49] [50, 52]. However, co-occurrence of triple homozygous point mutations (homozygous mutation for each of V1016G, F1534C and S989P) is very scarce having only been reported from Myanmar at a frequency of 0.98% [46] and in Indonesia in one individual (at a frequency of 0.34%) [50]
Despite all of the strategies implemented in Indonesia, the existing methods of controlling dengue have limited success. The development of a novel strategy of vector control that can be incorporated into the existing vector control strategy is essential. One of the innovative approaches to preventing transmission of dengue virus involves introduction of strains of the bacterium Wolbachia into A. aegypti, which has both life-shortening effects on the mosquito and direct transmission-blocking effects on dengue virus. This new strategy together with more traditional approaches for vector control including insecticide application may provide promising results to reduce dengue transmission. With the support of communities and approval from regulators, Indonesia’s first field trial of Wolbachia-infected mosquitoes began in Yogyakarta in 2014 [55]. The first field trial in two areas in the outer city of Yogyakarta has yielded results that point to local invasion [56], and release of Wolbachia mosquitoes on a broader scale in the inner-city areas of Yogyakarta has since been initiated in August 2016 [57].
To assist the spread of Wolbachia, insecticide resistance in the strain of A. aegypti being released needs to match that of the background population to ensure that released individuals persist and reproduce, allowing a Wolbachia invasion to take place [58]In the inner city of Yogyakarta, resistance is expected because of the local heavy application of chemical insecticides. It is possible that insecticide usage is higher in these areas than in the outer rim of Yogyakarta, as inner city areas have more dengue cases [59] and are more densely populated [60]. In our previous study, we screened samples from Yogyakarta outer areas for kdr mutations and obtained a high frequency of kdr alleles in samples collected from the outer area of Yogyakarta [49]. The most common kdr co-occurrence was V1016G homozygous mutant / F1534 wildtype, combined either with S989P heterozygous, or S989P wild type, or S989P homozygous mutant respectively. This co-occurrence was more common in the surviving than in the dead mosquitoes [49] and a strong association between the V1016G mutation with pyrethroid resistance (type I and type II) was confirmed. In contrast, F1534C mutant homozygotes were rare and there was only a weak association between heterozygote individuals for the F1534C mutation and resistance to a type I pyrethroid. The S989P mutation, in addition to the V101G mutation was found to have an additive effect in resistance to Type II pyrethroids [49]. In this paper, we have screened Yogyakarta inner city areas for kdr mutations, based on a total of 1314 individuals collected from 27 localities (Figure 1, Table 1) and aim to compare frequency and occurrence of kdr mutations between the inner and outer city.
Table 1. Collection sites, localities and number of A. aegypti collected from inner city areas of Yogyakarta.
Sites
|
Localities
|
No. of samples
|
Geolocation coordinates
(longitude, latitude)
|
C1
|
Suryodinigratan
|
112
|
110.35971
|
-7.81972
|
C2
|
Mantrijeron
|
25
|
110.36539
|
-7.82453
|
C3
|
Brontokusuman
|
25
|
110.36996
|
-7.82397
|
C4
|
Sorosutan
|
124
|
110.37840
|
-7.82886
|
C5
|
Kadipaten,
Patehan
|
34
|
110.35872
110.36015
|
-7.80587
-7.81104
|
C6
|
Panembahan,
Keparakan,
Prawirodirjan
|
25
|
110.36486
110.37327
110.37287
|
-7.81063
-7.81464
-7.80581
|
C7
|
Notoparajan
Ngampilan
Ngupasan
|
140
|
110.35560
110.35577
110.36199
|
-7.80621
-7.79728
-7.80199
|
C8
|
Wirogunan,
Sorosutan
|
25
|
110.37810
110.37872
|
-7.81130
-7.81484
|
C9
|
Tahunan, Semaki
|
19
|
110.38273
|
-7.81056
|
C10
|
Bausasran,
Gunung Ketur,
Baciro,
Purwokinanti
|
25
|
110.37448
110.37883
110.37875
110.37382
|
-7.79253
-7.80142
-7.79241
-7.79717
|
C11
|
Pringgokkusuman,
Sosromenduran
|
21
|
110.36969
110.36458
|
-7.79702
-7.79287
|
C12
|
Suryatmajan,
Tegal Panggung
|
124
|
110.36969
110.36983
|
-7.79702
-7.79350
|
C13
|
Bumijo,
Sosromenduran,
Gowongan
|
25
|
110.35676
110.36460
110.36467
|
-7.78798
-7.78835
-7.78353
|
C14
|
Gowongan,
Kotabaru
|
25
|
110.36921
110.36935
|
-7.78485
-7.78920
|
C15
|
Cokrodiningaratan,
Terban
|
88
|
110.36914
110.37384
|
-7.77882
-7.77877
|
C16
|
Semaki,
Baciro
|
68
|
110.38617
110.38211
|
-7.79694
-7.79685
|
C17
|
Klitren,
Demangan
|
25
|
110.38212
110.38219
|
-7.78394
-7.78927
|
C18
|
Terban,
Klitren
|
25
|
110.37876
110.38331
|
-7.77835
-7.78029
|
R1
|
Bener,
Kricak,
Karangwaru
|
31
|
110.352325
110.360298
110.364609
|
-7.77969
-7.77916
-7.77532
|
R2
|
Tegalrejo,
Pakuncen
|
38
|
110.35648
110.35070
|
-7.78339
-7.79298
|
R3
|
Wirobrajan,
Patangpuluhan
|
49
|
110.35088
110.34646
|
-7.80139
-7.80694
|
R4
|
Gedongkiwo
|
32
|
110.35344
|
-7.82552
|
B1
|
Muja Muju
|
27
|
110.39291
|
-7.79768
|
B2
|
Muja Muju,
Warung Boto,
Pandeyan
|
87
|
110.39259
110.38804 110.38647
|
-7.80637
-7.81086
-7.81439
|
B3
|
Pandeyan,
Giwangan
|
24
|
110.38669 110.39137
|
-7.82032
-7.83290
|
Co1
|
Rejojwinangun
|
39
|
110.40043
|
-7.81836
|
Co2
|
Prenggan, Purbayan
|
36
|
110.39775
|
-7.82482
|
|
|
|
|
|
|
Note that localities are mapped in Figure 1. The clusters refer to areas that are being used in a controlled intervention design to investigate the impact of Wolbachia on dengue incidence by the World Mosquito Program (Indonesia). Areas marked by an R are not included as part of this design.