Skip to main content

Insecticide resistance status of Aedes aegypti and Aedes albopictus mosquitoes in Papua New Guinea

Abstract

Background

Aedes aegypti and Ae. albopictus are important vectors of infectious diseases, especially those caused by arboviruses such as dengue, chikungunya and Zika. Aedes aegypti is very well adapted to urban environments, whereas Ae. albopictus inhabits more rural settings. Pyrethroid resistance is widespread in these vectors, but limited data exist from the Southwest Pacific Region, especially from Melanesia. While Aedes vector ecology is well documented in Australia, where incursion of Ae. albopictus and pyrethroid resistance have so far been prevented, almost nothing is known about Aedes populations in neighbouring Papua New Guinea (PNG). With pyrethroid resistance documented in parts of Indonesia but not in Australia, it is important to determine the distribution of susceptible and resistant Aedes populations in this region.

Methods

The present study was aimed at assessing Aedes populations for insecticide resistance in Madang and Port Moresby, located on the north and south coasts of PNG, respectively. Mosquitoes were collected using ovitraps and reared in an insectary. Standard WHO bioassays using insecticide-treated filter papers were conducted on a total of 253 Ae. aegypti and 768 Ae. albopictus adult mosquitoes. Subsets of samples from both species (55 Ae. aegypti and 48 Ae. albopictus) were screened for knockdown resistance mutations in the voltage-sensitive sodium channel (Vssc) gene, the target site of pyrethroid insecticides.

Results

High levels of resistance against pyrethroids were identified in Ae. aegypti from Madang and Port Moresby. Aedes albopictus exhibited susceptibility to pyrethroids, but moderate levels of resistance to DDT. Mutations associated with pyrethroid resistance were detected in all Ae. aegypti samples screened. Some genotypes found in the present study had been observed previously in Indonesia. No Vssc mutations associated with pyrethroid resistance were found in the Ae. albopictus samples.

Conclusions

To our knowledge, this is the first report of pyrethroid resistance in Ae. aegypti mosquitoes in PNG. Interestingly, usage of insecticides in PNG is low, apart from long-lasting insecticidal nets distributed for malaria control. Further investigations on how these resistant Ae. aegypti mosquito populations arose in PNG and how they are being sustained are warranted.

Background

Aedes aegypti and Ae. albopictus are vectors of important neglected infectious diseases, especially those caused by arboviruses such as dengue, chikungunya and Zika. They are known to be invasive and have expanded to several new regions relatively recently [1,2,3,4,5].

Insecticide resistance, especially against pyrethroids, is a major threat to vector-borne disease control worldwide. Aedes aegypti is well-adapted to urban habitats, and as a result, is usually more likely to be exposed to insecticides and develop resistance than is Ae. albopictus [6].

In the Asia-Pacific Region, reports have confirmed pyrethroid resistance in Aedes populations in many countries [7], and there is a geographical divide between regions where Aedes spp. are still susceptible (mainly Australia [8]) and those where Aedes spp. are resistant [8]. Aedes albopictus has not established populations in Australia so far [9].

While there are studies from parts of Indonesia indicating pyrethroid resistance in Ae. aegypti [10], Papua New Guinea (PNG) is one of the missing links in the region, where vector distribution and insecticide resistance status are almost completely unknown. While information about the vectors is very scarce [5], it is known that arboviruses represent a significant health burden in PNG; dengue viruses 1–4 circulate in the country and chikungunya outbreaks have been documented [11,12,13]. Due to the size of PNG, as well as its central location in the South Pacific Region, rapid population growth and increasing economic importance, understanding Aedes vector biology and resistance status in PNG is likely to be relevant to vector-control efforts in the South Pacific Region more generally. For instance, Aedes mosquitoes (and their resistance genes) may be transported from PNG to neighbouring countries, including Australia. Mosquito incursions to Australia are intercepted on a regular basis [14]. Aedes aegypti is present both in the Torres Strait and on the Australian mainland, but no phenotypic resistance or voltage-sensitive sodium channel (Vssc) resistance mutations have been detected in mosquitoes from these locations [8]. Although Ae. albopictus have not yet colonised the Australian mainland, they are present in the Torres Strait and PNG likely provides an important incursion route for this species [9, 15].

PNG exhibits a highly complex geography and a low level of urbanisation, with an estimated 13% of the population residing in urban areas [16]. However, this is rapidly changing as the expansion of urbanised areas, urban drift and anthropogenic environmental transformation are resulting in an increased proportion of humans at high risk from Aedes vector borne diseases. There is no widespread usage of insecticides in PNG in the public sector, apart from mass distribution of long-lasting insecticidal nets (LLIN) for malaria control since around 2009 [17], although large cities such as Port Moresby are not included in LLIN mass distributions. The LLIN distributed in PNG are treated with deltamethrin, which is a type 2 pyrethroid. In the private sector, insecticide usage is also limited, but not well documented. Large-scale commercial agriculture is not widespread and only a very small proportion of the available land (c.2.5%) is commercially farmed. Mining companies use fogging and indoor residual spraying (IRS) to protect their workforce from malaria. Hotels and businesses also engage in small-scale vector control activities, usually by hiring local pest control services. As such, there is potential for insecticide pressures in some areas in PNG to be underestimated.

The aim of this study was to generate initial information about Aedes vectors in PNG and their resistance levels to begin to address the knowledge gaps around these species. We collected Aedes mosquitoes from two population centres, Madang and Port Moresby, to determine species distribution and insecticide resistance status using bioassays to determine phenotype and DNA-based techniques to detect genetic markers of resistance.

Methods

Sampling sites and sample collection

Mosquitoes were collected as eggs and larvae in ovitraps placed in different locations around the provincial capital of Madang and the national capital of Port Moresby. In Madang, ovitraps (n = 15–20) were placed in the same locations every month for 12 months from February 2018 to January 2019. In Port Moresby traps were only placed once in December 2018 (n = 20). In both locations, traps were spread out across a large area in residential properties and hotel grounds, protected from rain. Trap contents were collected 5–7 days after trap placement.

Mosquito rearing and bioassays

Mosquitoes were reared to adult stage in a permanent (Madang) or temporary (Port Moresby) insectary. Bioassays were conducted as previously described [18] using WHO standard methodology, and limiting concentration filter papers obtained from School of Biological Sciences, Universiti Sains Malaysia, which is the regional centre manufacturing these commodities. Assays were always accompanied by at least one control of 20 mosquitoes (20 per cylinder). Where possible each test included 20 mosquitoes per assay (20 per cylinder). Ideally, we conducted 4 assays and 2 controls per insecticide in parallel; however, this was not always possible given the limited mosquito numbers. The average humidity in the assays was 76.1% (range: 67.8–82.3%) and the average temperature was 27.03 °C (range: 22.6–32.0 °C). In Madang, PNG, day length varies from approximately 12.35 h in winter to 13.10 h in summer (southern hemisphere).

Insecticides used in this study were deltamethrin (0.05%), DDT (4%) lambdacyhalothrin (0.05%), bendiocarb (0.1%) and malathion (5%). Aedes aegypti and Ae. albopictus were exposed to the filter papers as per standard protocol for 60 min, and then removed, and placed into holding cups with access to sugar solution for 24 h. Outcome variables were the proportion of mosquitoes knocked down after 60 min and the proportion of mosquitoes dead after the 24 h holding period. Since DDT has been shown to act more slowly on anophelines, we also assessed mortality 48 h after exposure [19].

In PNG, Ae. albopictus and Ae. aegypti occur mixed in containers and separation of live larvae or adults prior to conducting the bioassays is not possible. Therefore, we conducted the assays with mixed species and identified the adults morphologically after the assays.

Molecular analysis of kdr mutations

DNA was extracted from whole female mosquitoes using the High Pure PCR Template Preparation Kit (Roche, Millers Point, NSW, Australia) and suspended in 200 µl elution buffer. A 1:5 DNA dilution in PCR-grade H2O was prepared.

TaqMan® (Life Technologies Corporation, Mulgrave, VIC, Australia) assays were used to identify mutations in codons 989, 1016 and 1534 in the Vssc gene (numbered according to sequence of the most abundant splice variant of the house fly, Musca domestica) of Ae. aegypti. Primer and probe sequences used are as follows: Codon 989, forward primer (5′-TTC ATG ATC GTG TTC CGG GTA TT-3′), reverse primer (5′-ACG TCA CCC ACA AGC ATA CAA T-3′), probe (wildtype) (5′-CCC ACA TGG ATT CGA T-3′), probe (mutant) (5′-CCA CAT GGG TTC GAT-3′); codon 1016, forward primer (5′-CGT GCT AAC CGA CAA ATT GTT TCC-3′), reverse primer (5′-ATG AAC CGA AAT TGG ACA AAA GCA A-3′), probe (wildtype) (5′-AGA AAA GGT TAA GTA CCT GTG CG-3′), probe (mutant) (5′-AAG GTT AAG TCC CTG TGC G-3′); codon 1534, forward primer (5′-TCT ACA TGT ACC TCT ACT TTG TGT TCT TCA-3′), reverse primer (5′-GAT GAT GAC ACC GAT GAA CAG ATT C-3′), probe (wildtype) (5′-AAC GAC CCG AAG ATG A-3′), probe (mutant) (5′-ACG ACC CGC AGA TGA-3′). A 7 µl PCR reaction contained 40× TaqMan® assay (0.17 µl), 2× KAPA Fast PCR Probe Force qPCR Master Mix (KAPA Biosystems Inc.) (3.50 µl), ddH2O (1.33 µl) and genomic DNA (1:5 dilution) (2 µl).

Assays for each Vssc mutation were run in triplicate on a LightCycler® II 480 (Roche, Millers Point, NSW, Australia) instrument in 384-well plates with a pre-incubation of 3 min at 98 °C (ramp rate 4.8 °C/s) followed by 40 cycles of amplification at 95 °C for 10 s (2.5 °C/s ramp rate) and 60 °C for 20 s (2.5 °C/s ramp rate) (acquisition mode: single) with a final cooling step of 37 °C for 1 min (2.5 °C/s ramp rate). Endpoint genotyping was conducted for each mutation site (Roche LightCycler® 480 Software Version 1.5.1.62).

Mutations at Vssc codon 1534 in S6, Domain III of Ae. albopictus were screened using a forward primer of our own design (Alb171F: 5’-CCG ATT CGC GAG ACC AAC AT-3’) and the reverse primer of Kasai et al. (2011) (aegSCR8) [20]. Primer Alb171F was designed in the exon that contains Vssc codon 1534 in order to alleviate problems we encountered in sequencing across the intron. A 25 µl PCR mix included final concentrations of ThermoPol buffer Mg-free (1×) (New England Biolabs, Ipswich, MA, USA), dNTPs (200 µM each), MgCl2 (1.5 mM), 0.5 µM each of forward and reverse primers, 0.625 units of Immolase™ Taq polymerase (Bioline, London, UK), 2 µl genomic DNA (Chelex® extraction) diluted 1:10 and PCR-grade H2O to a final volume of 25 µl.

PCR conditions used to amplify the region were an initial denaturation at 95 °C for 10 min, 35 cycles of 95 °C for 30 s, annealing at 52 °C for 45 s and extension at 72 °C for 45 s, followed by a final extension step of 5 min at 72 °C and a hold at 10 °C.

PCR amplicons (220 bp) were sent to Macrogen Inc. in Seoul, Korea, for sequencing on a 3730xl DNA analyser. Sequences (up to 180 bp) were aligned and analysed using Geneious® 11.1.4 (Biomatters Ltd.) and mapped to a reference sequence (GenBank: KC152046.1).

Results

Overall species distribution

In Madang Ae. albopictus was observed to be the dominant species, whereas in Port Moresby Ae. aegypti was more abundant. Relative abundance of Ae. aegypti vs Ae. albopictus in December 2018, when mosquitoes were trapped in both locations simultaneously, was 19.6 vs 80.4% in Madang and 81.1 vs 18.9 % in Port Moresby, respectively, as shown in Fig. 1.

Fig. 1
figure 1

Location of the study sites and relative distribution of Aedes spp. in December 2018

Bioassays

Overall, 1021 female Ae. aegypti and Ae. albopictus mosquitoes were exposed to insecticides in the bioassays (n = 253 Ae. aegypti and n = 768 Ae. albopictus, excluding negative controls done alongside each assay). A total of 26 bioassays were conducted: with deltamethrin (0.05%) (n = 10); with DDT (4%) (n = 5); with lambda-cyhalothrin (0.05%) (n = 5); with bendiocarb (0.1%) (n = 3); and with malathion (5%) (n = 3). In Port Moresby, only two bioassays with deltamethrin were conducted. The results of the bioassays are presented in Tables 1 and 2 for Ae. aegypti and Ae albopictus, respectively.

Table 1 Bioassay results for Aedes aegypti, including samples collected in Madang and Port Moresby. Values are given as proportions and 95% CI [35]
Table 2 Bioassay results for Aedes albopictus, including samples collected in Madang and Port Moresby. Values are given as proportions and 95% CI [35]

These data clearly show high levels pyrethroid resistance in both Madang and Port Moresby Ae. aegypti populations. The Ae. aegypti population in Port Moresby seems to exhibit a higher level of resistance than that in Madang, with less than 10% of mosquitoes dying during the 24 h holding period after exposure to deltamethrin. Similarly, susceptibility to lambda-cyhalothrin and DDT was greatly reduced (although the number of Ae. aegypti exposed to DDT was low). Susceptibility against carbamates (bendiocarb) and organophosphates (malathion) was present in the Madang population (Port Moresby population not tested), but since the number of exposed mosquitoes was low, this result should be treated cautiously.

Susceptibility to DDT of Ae. albopictus also appeared to be reduced. The cut-off for confirmed resistance is < 90% mortality over the course of the 24 h holding period [21]. We observed a 24-h mortality of 0.79 (95% CI: 0.72–0.85). Aedes albopictus were found to be susceptible to pyrethroids, carbamates and organophosphates.

Genetic analyses

A total of 55 randomly selected Ae. aegypti mosquitoes from Madang were tested for sodium channel mutations using TaqMan® probe assays designed for these sites. The para sodium channel or voltage-sensitive sodium channel (Vssc) is the target site for both pyrethroid insecticides and DDT [22]. In addition, 48 Ae. albopictus mosquitoes were screened by sequencing for Vssc mutations at codon 1534, the main site of knockdown resistance mutations in this species [23].

For Ae. aegypti, three sodium channel mutations, V1016G, F1534C and S989P were found resulting in four composite genotypes. Composite genotype frequencies are summarised in Table 3.

Table 3 Frequencies of Vssc genotypes of 55 Aedes aegypti sampled from Madang, Papua New Guinea

The most common genotype (frequency = 0.65) consisted of homozygous mutants at codons 1016 and 989 and the wildtype homozygote at codon 1534. This is a common genotype also found in Bali and other locations throughout southeast Asia and the Pacific [24]. It confers resistance to both Type I and Type II pyrethroids.

The second most frequent genotype was the triple heterozygote at V1016G, F1534C and S989P (frequency = 0.16). This is also found commonly in Ae. aegypti in other countries in the region and confers a low level of resistance to Type I and II pyrethroids [25].

The remaining two genotypes found were (i) homozygous mutant at codon 1016, wildtype at codon 1534 and heterozygous mutant at codon 989 (frequency: 0.11); and (ii) heterozygous mutants at codons 1016 and 1534, and homozygous wildtype at codon 989 (frequency: 0.09). These genotypes are unusual, but have previously been noted in samples of Ae. aegypti from Yogyakarta, Indonesia [24]. Both these genotypes are expected to confer some level of pyrethroid resistance [26]. There were no susceptible wildtype individuals found in the sample and there is no indication that a wildtype haplotype exists in the population.

DNA sequences around codon 1534 were obtained for 36 Ae. albopictus specimens. None of the specimens showed a mutation within codon 1534. There were some synonymous mutations in other codons which comprised heterozygotes (IUPAC code Y for C/T) in nine individuals. Two individuals were homozygous synonymous mutants at codon 1528 and contained a single base mutation from C to T. The same mutation at this site was observed by Kasai et al. [20], but is not expected to affect susceptibility to pyrethroids. The sequences for the 36 Ae. albopictus are provided in Additional file 1.

Discussion

Insecticide resistance is a threat to vector borne disease control, and pyrethroids are still a widely used insecticide class, e.g. in LLIN for malaria control [27]. Pyrethroid resistance in the highly competent arbovirus vector Ae. aegypti is widely spread [28], and has also been confirmed in Ae. albopictus [28]. In the Asia-Pacific region, reports have confirmed pyrethroid resistance in Aedes populations in many countries including Indonesia, but not in Australia [8]. This implies that there is a geographical divide somewhere in the region, between areas where Aedes populations are still susceptible and those where Aedes populations are resistant [8, 9]. The present data from Port Moresby and Madang indicate that resistance in urban PNG populations may be widely spread, given that we found strong pyrethroid resistance in urban Ae. aegypti populations from the north and south coast of the country, and not a single wild type Ae. aegypti among the 55 samples tested. Australia has engaged in programs to prevent an incursion of Ae. albopictus via the Torres Strait for some time [9]. These efforts now seem even more important and should include Ae. aegypti, as imported Ae. aegypti originating in PNG may well carry resistant alleles. Other incursion routes from Asia also carry this risk [29].

In contrast to the Ae. aegypti data, we found no genetic mutations known to convey pyrethroid resistance in the Ae. albopictus mosquitoes tested. Moreover, Ae. albopictus mosquitoes were susceptible to deltamethrin and lambda-cyhalothrin in bioassays. We did observe reduced susceptibility of Ae. albopictus against DDT in Madang using bioassays. Similar results have been found in Ae. albopictus from India where DDT resistance is detected in bioassays, but is not associated with Vssc mutations [30]. DDT resistance in Ae. albopictus from PNG either may be persisting following previous use of the compound in the country during the malaria eradication era or may be due to incursion of resistant mosquitoes, an occurrence known from other parts of the world [31, 32].

Aedes aegypti are more prominent in urban than in rural areas, as these mosquitoes are well adapted to habitats found in urban settings [3]. This may account for the much higher proportion of Ae. aegypti found in the highly urbanised environment of the national capital Port Moresby compared with the proportion found in Madang (medium sized town, semi-urban). Aedes albopictus is known to use a wider range of habitats than Ae. aegypti and is better adapted to the natural environment than it is to urban areas.

The observed level of pyrethroid resistance in Ae. aegypti in the two PNG populations is surprising and unexpected, especially as LLIN usage in cities is lower than in rural areas. For example, no LLIN are distributed in Port Moresby (Tim Freeman, Rotary Against Malaria PNG, primary LLIN distribution agency in PNG, personal communication). No detailed analyses that quantify insecticide usage in PNG in the private sector exist, but to the best of our knowledge, insecticide usage is low. Indirect exposure to pesticides through the agricultural sector seems unlikely, particularly as subsistence farming usually does not involve the use of insecticides. Other sources of insecticide use in PNG are mining companies and other large commercial enterprises focused on resource extraction, as well as hotels and smaller businesses. Little is known about private household-based usage of insecticides, but this could provide further selection pressures. It is therefore possible that the observed pyrethroid resistance is focused on larger population centres or other commercial areas.

Another hypothesis arising from this work is that pyrethroid-resistant Ae. aegypti has spread to PNG from neighbouring regions. Human-mediated movement of pyrethroid-resistant Ae. aegypti is known to occur on a regular basis and mosquitoes with resistance alleles and identified, exotic origins are often intercepted at Australian and New Zealand air and sea ports [33]. It is possible that mosquitoes are transported on marine vessels used, e.g. for fishing and logging, as well as container transport [1] and so incursions to PNG could be possible. Pyrethroid resistance in Ae. aegypti is thought to be spreading in some parts of Indonesia [34], so potential sources of resistant mosquitoes in the region around PNG are likely to be on the increase. Further surveillance and sampling in PNG is required to map the level of resistance and population genetics of Ae. aegypti in different parts of the country and identify any potential pathways of incursion.

Limitations of this study include the absence of a susceptible Aedes laboratory strain to be used as a control, one of the challenges of working in a remote area with resource limitations. However, our observations were made on samples collected from multiple sites, tested in multiple assays. Bioassay results are also well supported by the genetic analyses presented. The concentrations on the test papers we used were those recommended for anopheline mosquitoes by WHO, whereas the concentrations (unofficially) recommended for Aedes mosquitoes are actually lower than those for anophelines [21]. As such our study is unlikely to overestimate the level of resistance but may represent an underestimation.

Conclusions

The present study identified pyrethroid resistance in Ae. aegypti in two locations in Papua New Guinea. To our knowledge, this is the first Aedes bioassay and kdr genotyping data from mosquito populations in Papua New Guinea. Due to presumably low insecticide usage in PNG, it is currently unclear how pyrethroid resistance in PNG arose and is maintained. Studies that investigate and quantify insecticides in the private sector in PNG are urgently needed. Population genetic studies that investigate the relationship between Aedes vectors in PNG and neighbouring countries would be useful to determine pathways of spread and incursion routes of resistant populations.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request. Genetic sequences generated during this study are included in Additional file 1.

Abbreviations

AE:

Aedes

DDT:

dichlorodiphenyltrichloroethane

IRS:

indoor residual spraying

KDR:

knockdown resistance

LLIN:

long-lasting insecticidal nets

PNG:

Papua New Guinea

VSSC:

voltage-sensitive sodium channel

WHO:

World Health Organisation

References

  1. Fonzi E, Higa Y, Bertuso AG, Futami K, Minakawa N. Human-mediated marine dispersal influences the population structure of Aedes aegypti in the Philippine Archipelago. PLoS Negl Trop Dis. 2015;9:e0003829.

    Article  Google Scholar 

  2. Hill MP, Axford JK, Hoffmann AA. Predicting the spread of Aedes albopictus in Australia under current and future climates: multiple approaches and datasets to incorporate potential evolutionary divergence. Aust Ecol. 2014;39:469–78.

    Article  Google Scholar 

  3. Kraemer MU, Sinka ME, Duda KA, Mylne AQ, Shearer FM, Barker CM, et al. The global distribution of the arbovirus vectors Aedes aegypti and Ae. albopictus. Elife. 2015;4:e08347.

    Article  Google Scholar 

  4. Lounibos LP, Juliano SA. Where vectors collide: the importance of mechanisms shaping the realized niche for modeling ranges of invasive Aedes mosquitoes. Biol Invasions. 2018;20:1913–29.

    Article  Google Scholar 

  5. Maynard AJ, Ambrose L, Cooper RD, Chow WK, Davis JB, Muzari MO, et al. Tiger on the prowl: invasion history and spatio-temporal genetic structure of the Asian tiger mosquito Aedes albopictus (Skuse, 1894) in the Indo-Pacific. PLoS Negl Trop Dis. 2017;11:e0005546.

    Article  Google Scholar 

  6. Ayorinde A, Oboh B, Oduola A, Otubanjo O. The insecticide susceptibility status of Aedes aegypti (Diptera: Culicidae) in farm and nonfarm sites of Lagos State, Nigeria. J Insect Sci. 2015;15:75.

    Article  Google Scholar 

  7. Amelia-Yap ZH, Chen CD, Sofian-Azirun M, Low VL. Pyrethroid resistance in the dengue vector Aedes aegypti in Southeast Asia: present situation and prospects for management. Parasit Vectors. 2018;11:332.

    Article  Google Scholar 

  8. Endersby-Harshman NM, Wuliandari JR, Harshman LG, Frohn V, Johnson BJ, Ritchie SA, et al. Pyrethroid susceptibility has been maintained in the dengue vector, Aedes aegypti (Diptera: Culicidae), in Queensland, Australia. J Med Entomol. 2017;54:1649–58.

    Article  CAS  Google Scholar 

  9. Muzari MO, Devine G, Davis J, Crunkhorn B, van den Hurk A, Whelan P, et al. Holding back the tiger: successful control program protects Australia from Aedes albopictus expansion. PLoS Negl Trop Dis. 2017;11:e0005286.

    Article  Google Scholar 

  10. 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:3688–92.

    Article  CAS  Google Scholar 

  11. Senn N, Luang-Suarkia D, Manong D, Siba PM, McBride WJ. Contribution of dengue fever to the burden of acute febrile illnesses in Papua New Guinea: an age-specific prospective study. Am J Trop Med Hyg. 2011;85:132–77.

    Article  Google Scholar 

  12. Horwood PF, Reimer LJ, Dagina R, Susapu M, Bande G, Katusele M, et al. Outbreak of chikungunya virus infection, Vanimo, Papua New Guinea. Emerg Infect Dis. 2013;19:1535–8.

    Article  Google Scholar 

  13. Moore PR, van den Hurk AF, Mackenzie JS, Pyke AT. Dengue viruses in Papua New Guinea: evidence of endemicity and phylogenetic variation, including the evolution of new genetic lineages. Emerg Microbes Infect. 2017;6(12):e114.

    PubMed  PubMed Central  Google Scholar 

  14. Webb C, Doggett SL, Russell R. A guide to mosquitoes of Australia. Clayton: CSIRO Publishing; 2016.

    Book  Google Scholar 

  15. van den Hurk AF, Nicholson J, Beebe NW, Davis J, Muzari OM, Russell RC, et al. Ten years of the tiger: Aedes albopictus presence in Australia since its discovery in the Torres Strait in 2005. One Health. 2016;2:19–24.

    Article  Google Scholar 

  16. Jones P. Managing urbanisation in Papua New Guinea: planning for planning’s sake? The Alfred Deakin Research Institute Working Papers; 2012. p. 33.

  17. Reimer LJ, Thomsen EK, Tisch DJ, Henry-Halldin CN, Zimmerman PA, Baea ME, et al. Insecticidal bed nets and filariasis transmission in Papua New Guinea. N Engl J Med. 2013;369:745–53.

    Article  CAS  Google Scholar 

  18. Koimbu G, Czeher C, Katusele M, Sakur M, Kilepak L, Tandrapah A, et al. Status of insecticide resistance in Papua New Guinea: an update from nation-wide monitoring of Anopheles mosquitoes. Am J Trop Med Hyg. 2017;98:162–5.

    Article  Google Scholar 

  19. Lukwa N, Sande S, Munosiyei P, Zimba M. Insecticide susceptibility tests conducted in Kamhororo, Masakadza and Chilonga villages in Zimbabwe during the 2011 malaria period. J Entomol Acarol Res. 2012;44:e19.

    Article  Google Scholar 

  20. Kasai S, Ng LC, Lam-Phua SG, Tang CS, Itokawa K, Komagata O, et al. First detection of a putative knockdown resistance gene in major mosquito vector, Aedes albopictus. Jpn J Infect Dis. 2011;64:217–21.

    CAS  PubMed  Google Scholar 

  21. WHO. Test procedures for insecticide resistance monitoring in malaria vector mosquitoes. Geneva: World Health Organization; 2013.

    Google Scholar 

  22. Soderlund DM, Bloomquist JR. Neurotoxic actions of pyrethroid insecticides. Annu Rev Entomol. 1989;34:77–96.

    Article  CAS  Google Scholar 

  23. Xu J, Bonizzoni M, Zhong D, Zhou G, Cai S, Li Y, et al. Multi-country survey revealed prevalent and novel F1534S mutation in voltage-gated sodium channel (VGSC) gene in Aedes albopictus. PLoS Negl Trop Dis. 2016;10:e0004696.

    Article  Google Scholar 

  24. Wuliandari JR, Lee SF, White VL, Tantowijoyo W, Hoffmann AA, Endersby-Harshman NM. Association between three mutations, F1565C, V1023G and S996P, in the voltage-sensitive sodium channel gene and knockdown resistance in Aedes aegypti from Yogyakarta, Indonesia. Insects. 2015;6:658–85.

    Article  Google Scholar 

  25. Du Y, Nomura Y, Satar G, Hu Z, Nauen R, He SY, et al. Molecular evidence for dual pyrethroid-receptor sites on a mosquito sodium channel. Proc Natl Acad Sci USA. 2013;110:11785–90.

    Article  CAS  Google Scholar 

  26. Plernsub S, Saingamsook J, Yanola J, Lumjuan N, Tippawangkosol P, Sukontason K, et al. Additive effect of knockdown resistance mutations, S989P, V1016G and F1534C, in a heterozygous genotype conferring pyrethroid resistance in Aedes aegypti in Thailand. Parasit Vectors. 2016;9:417.

    Article  Google Scholar 

  27. Hetzel MW, Pulford J, Ura Y, Jamea-Maiasa S, Tandrapah A, Tarongka N, et al. Insecticide-treated nets and malaria prevalence, Papua New Guinea, 2008–2014. Bull World Health Organ. 2017;95:695–705.

    Article  Google Scholar 

  28. Smith LB, Kasai S, Scott JG. Pyrethroid resistance in Aedes aegypti and Aedes albopictus: important mosquito vectors of human diseases. Pestic Biochem Physiol. 2016;133:1–12.

    Article  CAS  Google Scholar 

  29. Whelan P, Nguyen H, Hajkowicz K, Davis J, Smith D, Pyke A, et al. Evidence in Australia for a case of airport dengue. PLoS Negl Trop Dis. 2012;6:e1619.

    Article  Google Scholar 

  30. Kushwah RB, Mallick PK, Ravikumar H, Dev V, Kapoor N, Adak TP, et al. Status of DDT and pyrethroid resistance in Indian Aedes albopictus and absence of knockdown resistance (kdr) mutation. J Vector Borne Dis. 2015;52:95–8.

    CAS  PubMed  Google Scholar 

  31. Kamgang B, Marcombe S, Chandre F, Nchoutpouen E, Nwane P, Etang J, et al. Insecticide susceptibility of Aedes aegypti and Aedes albopictus in Central Africa. Parasit Vectors. 2011;4:79.

    Article  Google Scholar 

  32. Marcombe S, Farajollahi A, Healy SP, Clark GG, Fonseca DM. Insecticide resistance status of United States populations of Aedes albopictus and mechanisms involved. PLoS ONE. 2014;9:e101992.

    Article  Google Scholar 

  33. Schmidt TL, van Rooyen AR, Chung J, Endersby-Harshman NM, Griffin PC, Sly A, et al. Tracking genetic invasions: genome-wide single nucleotide polymorphisms reveal the source of pyrethroid-resistant Aedes aegypti (yellow fever mosquito) incursions at international ports. Evol Appl. 2019. https://doi.org/10.1111/eva.12787.

    Article  Google Scholar 

  34. Hamid PH, Ninditya VI, Prastowo J, Haryanto A, Taubert A, Hermosilla C. Current status of Aedes aegypti insecticide resistance development from Banjarmasin, Kalimantan, Indonesia. Biomed Res Int. 2018;2018:1735358.

    Article  CAS  Google Scholar 

  35. Agresti A, Coull BA. Approximate is better than “Exact” for interval estimation of binomial proportions. Am Stat. 1998;52:119–26.

    Google Scholar 

Download references

Acknowledgements

We would like to thank PNGIMR and Rotary Against Malaria employees for their assistance with placing mosquito traps in their premises in Madang and Port Moresby. We would like to thank Thomas Burkot, Tanya Russell, Scott Ritchie, Uriel Kitron and Tim Freeman for helpful discussions. We thank the Strive PNG partners for their support.

Funding

This study received no dedicated funding except that molecular assays were supported by an NHMRC Program grant and AH was supported by an NHMRC Fellowship. SK was supported by an NHMRC Career Development Fellowship GNT1141441, LT was supported by a JCU PhD Scholarship and LJR is supported by an NHMRC Career Development Fellowship. SD was supported, in part, by Strive PNG (Australian Department of Foreign Affairs and Trade).

Author information

Authors and Affiliations

Authors

Contributions

SD, RV and LT: field sample collection by ovitrapping and conduction phenotypic bioassays; data analysis. NEH and AH: molecular assays and analysis of molecular data, provided expertise with molecular techniques. LR provided expertise with field work and data analysis. SK and ML conceived the study; conducted data analysis; prepared the first manuscript draft; provided expertise with field work. SD, NEH, RV, LT, MS, LM, LJR, ML, AH and SK reviewed and contributed to final manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Stephan Karl.

Ethics declarations

Ethics approval and consent to participate

The present research was approved by the PNGIMR Institutional Review Board and the National Medical Research Advisory Committee of Papua New Guinea (Protocol number 18.10).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional file

Additional file 1.

DNA sequences for 36 individuals of Aedes albopictus for a small section of the voltage-sensitive sodium channel gene (Vssc) from S6, domain III.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Demok, S., Endersby-Harshman, N., Vinit, R. et al. Insecticide resistance status of Aedes aegypti and Aedes albopictus mosquitoes in Papua New Guinea. Parasites Vectors 12, 333 (2019). https://doi.org/10.1186/s13071-019-3585-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13071-019-3585-6

Keywords