Research | Open | Published:
Behaviour and molecular identification of Anopheles malaria vectors in Jayapura district, Papua province, Indonesia
Malaria Journalvolume 15, Article number: 192 (2016)
Members of the Anopheles punctulatus group dominate Papua, Indonesia and Papua New Guinea (PNG), with a geographic range that extends south through Vanuatu. An. farauti and An. punctulatus are the presumed major vectors in this region. Although this group of species has been extensively studied in PNG and the southern archipelagoes within their range, their distribution, ecology and vector behaviours have not been well characterized in eastern Indonesia.
Mosquitoes were collected in five villages in Jayapura province, Papua, Indonesia using human-landing collections, animal-baited tents and backpack aspirators. Mosquitoes were morphologically typed and then molecularly distinguished based on ribosomal ITS2 sequences and tested for Plasmodium falciparum and P. vivax infection using circumsporozoite ELISA and PCR.
The presence and vector status of An. farauti 4 in Papua, Indonesia is confirmed here for the first time. The data indicate that this species is entering houses at a rate that increases its potential to come into contact with humans and act as a major malaria vector. An. farauti 4 was also abundant outdoors and biting humans during early evening hours. Other species collected in this area include An. farauti 1, An. hinesorum, An. koliensis, An. punctulatus, and An. tessellatus. Proboscis morphology was highly variable within each species, lending support to the notion that this characteristic is not a reliable indicator to distinguish species within the An. punctulatus group.
The vector composition in Papua, Indonesia is consistent with certain northern areas of PNG, but the behaviours of anophelines sampled in this region, such as early and indoor human biting of An. farauti 4, may enable them to act as major vectors of malaria. Presumed major vectors An. farauti and An. punctulatus were not abundant among these samples. Morphological identification of anophelines in this sample was often inaccurate, highlighting the importance of using molecular analysis in conjunction with morphological investigations to update keys and training tools.
The province of Papua, Indonesia has a high burden of malaria, with the highest reported malaria incidence among Indonesia’s 34 provinces. Both Plasmodium falciparum and P. vivax are prevalent in the province of Papua and throughout eastern Indonesia [1–3]. P. vivax, while often not a direct cause of mortality, largely contributes to the disease burden and morbidity in this region [4, 5]. The primary vectors of P. vivax are not well known, as mosquito surveys do not always include screening for P. vivax infection along with P. falciparum. While there have been many studies of the ecology of the An. punctulatus group in neighbouring Papua New Guinea (PNG), sampling of malaria vectors and characterization of their ecology and behaviour in Indonesian Papua is lacking. In particular, comprehensive morphological and molecular analysis of anophelines in this province is uncommon.
There are currently 13 recognized members of the An. punctulatus group and eight recognized members of the An. farauti complex [6–8]. Only An. farauti 1 (hereafter referred to as An. farauti) and An. punctulatus are considered to be the primary malaria vectors. Some members of the An. punctulatus group can be distinguished through different morphological features, such as spotting on the palpi [9, 10], but the members of the An. farauti complex are morphologically indistinguishable and can only be typed to species using molecular tools . There is significant morphological variation within even molecularly distinct species across their geographic range. There is also phenotypic variation within this complex, such as increased saline tolerance in An. farauti and An. farauti 7 and large body size in An. oreios (formerly An. farauti 6) [12–14]. An. farauti is considered anthropophilic throughout its range and can exploit slightly saline larval habitats along the coast, unlike many other members of the An. punctulatus group . Since most of these species have been shown to be capable of harbouring both P. falciparum and P. vivax parasites, the primary determinate of their status as a major vector of malaria is by their tendency to bite humans indoors or near human dwellings . An. farauti 4 (like most members of the An. punctulatus group) is known to be primarily outdoor-biting and outdoor-resting . An. farauti 4 and An. hinesorum (formerly An. farauti 2) have been found to have P. falciparum positivity rates comparable to An. farauti and An. punctulatus in Papua New Guinea [17, 18]. Species within the An. punctulatus group can also act as vectors of several different types of human filariasis .
Current malaria control efforts in Indonesia are primarily dependent upon the widespread use of long-lasting insecticide-treated bed nets (LLINs), which do not target outdoor-biting mosquito species [6, 20, 21]. Nonetheless, Indonesia Ministry of Health statistics show a marked decline in malaria incidence consequent to LLIN distribution and improved diagnosis and malaria case management over the past several years. A high prevalence of malaria was reported in this district in 2008 , and in 2014 was reported to have among the highest transmission rates in Indonesia (annual parasite incidence of 172 reported cases/1000 population in 2014) (Nyoman, Kabupatan, Jayapura, pers comm). Previous mosquito surveys in Papua, Indonesia, have not distinguished between cryptic members of the An. farauti complex  and most vector studies in this region pre-date the understanding of cryptic species complexes.
It is crucial to understand the bionomic traits of each vector species throughout their geographic range to be able to implement effective malaria control and elimination efforts, particularly since many of the species in the An. punctulatus group occur sympatrically . This study adds to what is known about the An. punctulatus group within Indonesia through a brief survey and characterization of Anopheles species in five villages in Jayapura Province, Papua, Indonesia. The study illustrates the limits of reliance upon morphological identifications only for vector surveys in this particularly complex region. There are implications of the results for malaria control strategies in Indonesian Papua.
Papua, Indonesia is the easternmost province of Indonesia, bordering PNG. The eastern half of the island of Papua comprises PNG, while the western half comprises the Indonesian provinces of Papua and West Papua. The Malaria Transmission Consortium (MTC) collection sites in Papua province are in a series of villages along the coast near the provincial capital city Jayapura, very near to the border of PNG. Activities in this region include fishing, subsistence farming and harvesting spices in the jungle. Houses in this area are made of bamboo or thin pieces of wood with thatched roofs, many of which are up on stilts. Few houses have metal roofs. The open construction of these houses allows for mosquito entry from the bottom, sides and top. The health system in Papua is generally weaker than in other parts of Indonesia, which impacts surveillance, monitoring and treatment of malaria. Malaria transmission in this area is stable and high [3, 25].
Mosquitoes were collected in five villages in Jayapura province, Papua, Indonesia during May of 2011: Demta village, Demta sub-district; Bunyom village, Nimbokrang sub-district; Kehiran village, Sentani Kota sub-district; Nolokla village, Sentani Timur (Harapan) sub-district, and Ongan Jaya village, Yapsi (Taja) sub-district (Fig. 1). The villages were spread across the district of Jayapura with Demta, Bunyom and Ongan Jaya villages being in forested areas close to rivers, while Kehiran and Nolokla villages were in agricultural areas.
Several methods were utilized to collect anophelines. Backpack aspirators  were utilized hourly (from 18:00 to 06:00 hours) for 10 min to capture mosquitoes on bushes and foliage within and surrounding housing areas. Hourly human-landing catches (HLC) were performed both indoors and outdoors (from 18:00 to 06:00 hours). There were two indoor and two outdoor HLC collectors per night per site. Net tents were constructed using untreated netting with two wide openings situated around pig or cow baits. Anopheles mosquitoes were collected from the inside of the tent for ten min every hour (from 18:00 to 06:00 hours) using a mouth aspirator. The baited tents were situated away from human gathering areas and the human landing collectors. Collectors were advised not to wear any insect repellant during collections. Collections were performed for one night in each village except Kehiran, which had three nights of collections. All mosquitoes were collected, stored in holding cups labelled by hour until they were processed for morphological identification using keys for Indonesian anophelines . High-quality adult specimens from the field collections were individually processed for detailed morphological descriptions of the palpi and proboscis, a key diagnostic feature of some members of the An. punctulatus group.
Molecular processing of samples
Genomic DNA was isolated from individual specimens using a CTAB DNA extraction. Species were molecularly identified using sequences of the ribosomal DNA internal transcribed spacer region two (rDNA ITS2). This region of rDNA was isolated using PCR with ITS2A and ITS2B primers . The amplified fragments were purified using an enzyme clean-up: 2U of Exonuclease 1 (USB Corporation, Cleveland, OH), 1U of Shrimp Alkaline Phosphatase (USB), and 1.8 μl of ddH20 were added to 8 μl of PCR product. This mixture was incubated at 37 °C for 15 min, followed by 15 min at 80 °C to inactivate the enzymes. PCR products were sequenced directly using Sanger sequencing on ABI 3730 xl DNA Analyzer platform (Applied Biosystems).
The ITS2 sequences were blasted against the NCBI database with BLASTn for confirmation of molecular species identification. Sequences were visually checked for sequence quality and potential contamination. Low quality or contaminated sequences were excluded from the analysis. Sequences with greater than 99 % sequence identity to voucher reference sequences were confirmed as that molecular species. Voucher sequences are available in the NCBI database for the species sampled in this study. The primary PCR diagnostics for these species is also based on the sequence of the rDNA ITS2 region.
The infection status of the mosquitoes was determined using the standard CDC sandwich CS-ELISA test for the detection of P. falciparum, P. vivax-210, and P. vivax-247 circumsporozoite (CS) proteins . A sub-set of specimens were analysed for Plasmodium infection using a multiplex PCR for P. falciparum and P. vivax .
Statistical analyses to evaluate indoor and outdoor biting of An. farauti 4 and An. koliensis were performed using SPSS.
A total of 1071 (of 1968) individual mosquito specimens were successfully molecularly identified from Jayapura district, Papua province, Indonesia. Five morphological species were identified to six molecular species with an accuracy rate of 51 % (Table 1). These molecular species include An. farauti, An. farauti 4, An. hinesorum, An. koliensis, An. punctulatus, and An. tessellatus and were determined with sequence identity, greater than 99 % with reference sequences from voucher specimens.
There was a low level of accuracy when the specimens were identified morphologically in the field, with only 51 % accuracy using morphology alone. The few An. farauti and An. punctulatus specimens collected were morphologically identified accurately, while three specimens of the even more distinctive An. tessellatus was identified to three different species. An. koliensis, abundant in these collections and morphologically distinguishable from the other species, was identified as An. koliensis 13 % of the time. The most abundant species in the collection, An. farauti 4, was morphologically identified to An. farauti s.l. with 48 % accuracy. An. koliensis may be a species complex based on rDNA ITS2 sequence variation of samples found in other areas [18, 30], but this level of sequence variation was not observed within any of the molecularly identified specimens in this collection. These species have very similar and overlapping morphological characteristics, which can be further complicated by the quality of field-collected adult specimens. This level of accuracy in morphologically identifying field specimens when compared to molecular species identification is consistent with field collections in other areas of Southeast Asia and even countries within high transmission areas of Africa  (MTC data). This low level of morphological identification accuracy underlies the importance of incorporating molecular tools to help distinguish vector species.
The HLC collections indicate indoor and outdoor exposure of residents to host-seeking anophelines. Here, the majority of the species were collected both outdoors and early in the night, around dusk, at a time when people are extremely active, cooking and socializing outdoors. About 752 (74 %) of the total 1017 specimens collected were captured before midnight, with 394 (39 %) captured before 22.00. Only 299 (29 %) were caught between midnight and 06.00 (Fig. 2). This early evening biting behaviour is consistent with observations in PNG . The total Anopheles species collected were more evenly distributed through the night. The An. farauti 4 in this sample set were mostly collected in outdoor HLCs with 88 % of the total An. farauti 4 collected captured in outdoor HLCs and 9 % in indoor HLCs. This species may be much more prone to biting humans than was previously suspected (Fig. 3). The second-most abundant species, An. koliensis, was captured primary in indoor HLCs, 45 and 35 % in outdoor HLCs. Backpack aspiration captured all six species sampled, but with relatively low numbers, and was not representative of human exposure to infectious bites when compared to indoor or outdoor HLCs (Fig. 2). Eighty percent of all of the Anopheles samples collected in this survey were captured in outdoor HLCs.
More than half of the specimens caught in indoor and outdoor HLCs, and animal-baited traps were An. farauti 4 (Fig. 3). Fourteen percent of all of the anophelines were collected indoors, with the majority collected outdoors. An. farauti 4 was 54 % of indoor HLCs, 93 % of outdoor HLCs, 56 % of animal-baited traps, and 17 % of backpack-aspirated mosquitoes from vegetation. An. koliensis was the second-most abundant species, representing 45, 6, and 33 % of indoor and outdoor HLCs, and animal baited tents, respectively. An. koliensis was the most abundant species captured in vegetation, representing 62 % of the catch. All other species were captured in very low abundance in each trap type. An. farauti 1 was only collected outdoors, with six specimens captured in HLC and three in vegetation (Fig. 3).
To evaluate trap efficacy for the two most prevalent species in our collections, An. farauti 4 and An. koliensis, multinomial regression was performed to compare these two species captured in indoor HLC, outdoor HLC, and backpack aspiration of surrounding foliage, which were performed every hour from 18:00 to 06:00 hours. The trap effect was not significant for An. koliensis (p = 0.021), and was significant for An. farauti 4 (p ≤ 0.0001). An. farauti 4 was more likely to be found outdoor HLC than indoor HLC, and more likely to be caught in indoor HLC than in backpack aspiration of vegetation. A binomial logistic regression of these same two species was performed to compare their propensity to be captured indoors vs. outdoors. An. farauti 4 was 9.85 times more likely to be found outdoors than An. koliensis. Hourly biting profiles of An. farauti 4 were compared using a Chi square goodness of fit test. More An. farauti 4 were caught indoors before midnight than between the hours of midnight and 6 am (χ2: 15.12, df = 1, p = 0.0001). The number of An. farauti 4 caught indoors was significantly higher between 23 and 24 h when compared to other hours (χ2 = 48.418; df = 1, p < 0.0001) (Additional file 1). The frequency of An. farauti 4 outdoors before midnight was higher than between the hours of midnight and 6 am (χ2: 392.38, df = 1, p = <0.0001).
Four of the species collected in this survey are considered to be zoophilic, An. hinesorum, An. farauti 4, An. koliensis, and An. tessellatus and they were all captured in the animal-baited tents (Fig. 3). Animal-baited tents captured only 3 % of the total anophelines. Two of the four total An. tessellatus samples collected were in the animal-baited tents. Only 2 and 8 % of the total An. farauti 4, and An. koliensis specimens were captured in the animal-baited tents, with the vast majority captured in HLCs, indicating that both of these species have a higher human-biting preference than previously thought in this region. There were village-specific effects of species distribution and biting rates (Additional file 2). Futher sampling must be done to determine location-specific effects on species catch.
Human-biting rates were exceptionally high in Kehiran I village with about 84 bites per person per night outdoors and about 19 indoors. All sporozoite-positive mosquitoes were captured in this village, as were 98.7 % of An. farauti 4 and 51.7 % of An koliensis captured. Four individual specimens from these collections were found to be carrying P. vivax sporozoites, three An. farauti 4, and one An. koliensis (Table 2). In this region, An. farauti and An. punctulatus are considered to be the primary vectors [8, 17, 32]. These two vector species represented a very small portion of the collection, only 0.8 % of the total collection. The majority of specimens collected were An. farauti 4, representing 89 % of the total collection and also found to be sporozoite positive for P. vivax. The three P. vivax-positive An. farauti 4 and one An. koliensis were collected in outdoor HLCs (Table 2), with 0.4 and 1.6 % positivity rates, respectively.
The molecularly confirmed An. farauti 4 specimens were identified morphologically as An. punctulatus 39 % of the time (Table 1). This is an issue for both species identification and the allocation of vector bionomic traits to the appropriate species. In this case, without molecular species identification, the transmission of malaria and early biting behaviour by An. farauti 4 would be attributed to An. farauti s.l. and An. punctulatus. An. opheles farauti 4 is not considered to be a major malaria vector, although it has been shown to have high P. falciparum and P. vivax infection rates during mosquito surveys in PNG [17, 33] and has also been identified to be extremely polymorphic for proboscis characteristics which differentiate members of the An. punctulatus group . Although these species have some overlapping morphological and behavioural traits, they are evolutionarily distinct .
Members of the An. punctulatus group have very similar morphological characteristics, often only being distinguishable by banding patterns on the underside of the proboscis, a key morphological characteristic for this group [8, 9]. Even within these collections in Papua, Indonesia, additional morphological variants of these proboscis types were identified (Table 3; Additional file 3). Detailed morphology of the proboscis of individual specimens revealed morphological variation, even within cryptic species . This high level of morphological variation within species has been observed elsewhere within a single cryptic species in Indonesia, such as An. epiroticus and An. vagus (S Zubaidah, unpublished). The An. farauti 4 specimens in the collections were particularly morphologically variable for this trait, with specimens with each well-known phenotype, and several new phenotypes, represented within this species in the collection (Table 3; Additional file 3) . An. farauti and An. punctulatus had the expected A1 and B phenotypes, respectively, while each other species had unexpected or mixed phenotypes. The assortment of proboscis phenotypes, within even cryptic species, supports the need for molecular tools to distinguish all species in the An. punctulatus group, particularly in areas where these species occur in sympatry.
Results of this survey illustrate the necessity for use of molecular tools to ensure accurate species identification in Indonesian Papua. Morphological identification of these specimens was wrong nearly as often it was correct, highlighting the need for supplementation of morphological keys with molecular information to allow production of better training tools for morphological identification. Morphologically based routine mosquito surveillance is needed and inevitable, but updated keys with periodic integration of molecular identification would minimize inaccuracies. There are an increasing number of species diagnostic tools that are decreasing in cost per sample, including multiplex PCRs, to identify members of the An. punctulatus group. Tools such as these should be incorporated into malaria control programmes in this region to accurately identify vector species and associate them with their specific bionomic traits in order to efficiently target them, particularly in an area with such high transmission rates [6, 25, 35]. Without distinguishing field-collected specimens at a molecular level, the indoor and human-biting behaviours of An. farauti 4 and An. koliensis found in this study would have been attributed to An. farauti s.l. and assumed to be An. farauti.
Although An. farauti 4 and An. koliensis were implicated as malaria vectors in this survey, it remains unclear to what extent this result can be generalized. Samples were collected in 1 month (May 2011) and nearly all of the An. farauti 4 specimens came from a single village. Although this data concludes that this species transmits malaria in Papua province, Indonesia, it is not yet know how widely distributed the species is in Indonesian Papua, or even in Jayapura district. Further studies will be needed to know if the species is seasonally abundant or regularly infected, and information on its larval habitats is lacking. It is unknown whether a change in species abundance of different anophelines may have occurred in Papua subsequent to massive scale-up of malaria control interventions targeting indoor-biting species. Such information would allow the National Malaria Control Programme to make confident inferences as to which interventions might be most effective in the future in this last redoubt of malaria in Indonesia.
Arnott A, Barnadas C, Senn N, Siba P, Mueller I, Reeder JC, et al. High genetic diversity of Plasmodium vivax on the north coast of Papua New Guinea. Am J Trop Med Hyg. 2013;89:188–94.
Elyazar IR, Gething PW, Patil AP, Rogayah H, Sariwati E, Palupi NW, et al. Plasmodium vivax malaria endemicity in Indonesia in 2010. PLoS ONE. 2012;7:e37325.
Gething PW, Elyazar IR, Moyes CL, Smith DL, Battle KE, Guerra CA, et al. A long neglected world malaria map: Plasmodium vivax endemicity in 2010. PLoS Negl Trop Dis. 2012;6:e1814.
Douglas NM, Pontororing GJ, Lampah DA, Yeo TW, Kenangalem E, Poespoprodjo JR, et al. Mortality attributable to Plasmodium vivax malaria: a clinical audit from Papua, Indones. BMC Med. 2014;12:217.
Karyana M, Burdarm L, Yeung S, Kenangalem E, Wariker N, Maristela R, et al. Malaria morbidity in Papua Indonesia, an area with multidrug resistant Plasmodium vivax and Plasmodium falciparum. Malar J. 2008;7:148.
Beebe NW, Russell T, Burkot TR, Cooper RD. Anopheles punctulatus group: evolution, distribution, and control. Annu Rev Entomol. 2015;60:335–50.
Beebe NW, Cooper RD, Morrison DA, Ellis JT. Subset partitioning of the ribosomal DNA small subunit and its effects on the phylogeny of the Anopheles punctulatus group. Insect Mol Biol. 2000;9:515–20.
Cooper RD, Waterson DG, Frances SP, Beebe NW, Sweeney AW. Speciation and distribution of the members of the Anopheles punctulatus (Diptera: Culicidae) group in Papua New Guinea. J Med Entomol. 2002;39:16–27.
O’Connor CT, Soepanto A. Illustrated key to female Anophelines of Indonesia. Jakarta: Directorate of Communicable Disease, Ministry of Health; 1989.
Cooper RD, Waterson DG, Frances SP, Beebe NW, Sweeney AW. The anopheline fauna of Papua New Guinea. J Am Mosq Control Assoc. 2006;22:213–21.
Beebe NW, Cooper RD. Systematics of malaria vectors with particular reference to the Anopheles punctulatus group. Int J Parasitol. 2000;30:1–17.
Sweeney AW. Larval salinity tolerances of the sibling species of Anopheles farauti. J Am Mosq Control Assoc. 1987;3:589–92.
Bangs MJ, Taai K, Howard TM, Cook S, Harbach RE. The mosquito Anopheles (Cellia) oreios sp. n., formerly species 6 of the Australasian Anopheles farauti complex, and a critical review of its biology and relation to disease. Med Vet Entomol. 2015;29:68–81.
Foley DH, Bryan JH. Shared salinity tolerance invalidates a test for the malaria vector Anopheles farauti s.s. on Guadalcanal, Solomon Islands [corrected]. Med Vet Entomol. 2000;14:450–2.
Beebe NW, Bakote’e B, Ellis JT, Cooper RD. Differential ecology of Anopheles punctulatus and three members of the Anopheles farauti complex of mosquitoes on Guadalcanal, Solomon Islands, identified by PCR-RFLP analysis. Med Vet Entomol. 2000;14:308–12.
Charlwood JD, Graves PM, Alpers MP. The ecology of the Anopheles punctulatus group of mosquitoes from Papua New Guinea: a review of recent work. P N G Med J. 1986;29:19–26.
Cooper RD, Waterson DG, Frances SP, Beebe NW, Pluess B, Sweeney AW. Malaria vectors of Papua New Guinea. Int J Parasitol. 2009;39:1495–501.
Benet A, Mai A, Bockarie F, Lagog M, Zimmerman P, Alpers MP, et al. Polymerase chain reaction diagnosis and the changing pattern of vector ecology and malaria transmission dynamics in papua new Guinea. Am J Trop Med Hyg. 2004;71:277–84.
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.
Elyazar IR, Hay SI, Baird JK. Malaria distribution, prevalence, drug resistance and control in Indonesia. Adv Parasitol. 2011;74:41–175.
Hii JL, Smith T, Vounatsou P, Alexander N, Mai A, Ibam E, et al. Area effects of bednet use in a malaria-endemic area in Papua New Guinea. Trans R Soc Trop Med Hyg. 2001;95:7–13.
Sorontou Y, Asih PB, Wanandi SI, Ramelan W, Syafruddin D. Malaria in Jayapura district, Papua province, Indonesia and resistance to sulfadoxine-pyrimethamine. Med J Indones. 2007;16:32–8.
Hii J, Smith T, Mai A, Ibam E, Alpers M. Comparison between anopheline mosquitoes (Diptera: Culicidae) caught using different methods in a malaria endemic area of Papua New Guinea. Bull Entomol Res. 2000;90:211–9.
Logue K, Small ST, Chan ER, Reimer L, Siba PM, Zimmerman PA, et al. Whole-genome sequencing reveals absence of recent gene flow and separate demographic histories for Anopheles punctulatus mosquitoes in Papua New Guinea. Mol Ecol. 2015;24:1263–74.
Elyazar IR, Gething PW, Patil AP, Rogayah H, Kusriastuti R, Wismarini DM, et al. Plasmodium falciparum malaria endemicity in Indonesia in 2010. PLoS ONE. 2011;6:e21315.
Vazquez-Prokopec GM, Galvin WA, Kelly R, Kitron U. A new, cost-effective, battery-powered aspirator for adult mosquito collections. J Med Entomol. 2009;46:1256–9.
Beebe NW, Saul A. Discrimination of all members of the Anopheles punctulatus complex by polymerase chain reaction–restriction fragment length polymorphism analysis. Am J Trop Med Hyg. 1995;53:478–81.
Burkot TR, Zavala F, Gwadz RW, Collins FH, Nussenzweig RS, Roberts DR. Identification of malaria-infected mosquitoes by a two-site enzyme-linked immunosorbent assay. Am J Trop Med Hyg. 1984;33:227–31.
Singh B, Bobogare A, Cox-Singh J, Snounou G, Abdullah MS, Rahman HA. A genus- and species-specific nested polymerase chain reaction malaria detection assay for epidemiologic studies. Am J Trop Med Hyg. 1999;60:687–92.
Beebe NW, Cooper RD, Morrison DA, Ellis JT. A phylogenetic study of the Anopheles punctulatus group of malaria vectors comparing rDNA sequence alignments derived from the mitochondrial and nuclear small ribosomal subunits. Mol Phylogenet Evol. 2000;17:430–6.
Stevenson J, St Laurent B, Lobo NF, Cooke MK, Kahindi SC, Oriango RM, et al. Novel vectors of malaria parasites in the western highlands of Kenya. Emerg Infect Dis. 2012;18:1547–9.
Mueller I, Namuigi P, Kundi J, Ivivi R, Tandrapah T, Bjorge S, et al. Epidemic malaria in the highlands of Papua New Guinea. Am J Trop Med Hyg. 2005;72:554–60.
Charlwood JD, Birley MH, Dagoro H, Paru R, Holmes PR. Assessing survival rates of Anopheles farauti (diptera: Culicidae) from Papua New Guinea. J Anim Ecol. 1985;54:1003–16.
Rozeboom LE, Knight KL. The punctulatus complex of Anopheles (Diptera: Culicidae). J Parasitol. 1946;32:95–131.
Barcus MJ, Basri H, Picarima H, Manyakori C, Sekartuti Elyazar I, et al. Demographic risk factors for severe and fatal vivax and falciparum malaria among hospital admissions in northeastern Indonesian Papua. Am J Trop Med Hyg. 2007;77:984–91.
Bryan JH. Morphological studies on the Anopheles punctulatus Dönitz complex. Trans R Entomol Soc Lond. 1974;125:413–35.
BSTL and NFL carried out the sequence analysis, interpretation of the data, and drafted the manuscript. NFL, DS, SS, and FHC contributed to the conception and design of the study. DS, WAH, SS, PBSA, AS, FL, AB, and NFL participated in the implementation of the study and entomological collections. BSTL, PBSA, DB, JM, and HCM participated in the molecular analysis and sequence data generation. MN and NFL did the statistical analyses. All authors contributed to critical revision of the manuscript and intellectual content. All authors read and approved the final manuscript.
This project has been funded in whole or in part by the Bill and Melinda Gates Foundation under Grant No. 45114. Special thanks to field assistants, collectors and laboratory assistants for specimen processing at the Eijkman Institute of Molecular Biology, Jakarta, Indonesia and the University of Notre Dame. Thanks to Jillian Mueller for help creating Fig. 1 and Catherine Craker for creating the proboscis phenotype additional figure. Thanks to Julie Niedbalski and Katie Cybulski for help and support at the University of Notre Dame and to Rita Kusriastuti and P R Arbani for their support.
The authors declare that they have no competing interests.