Skip to main content

Diversity of Anopheles species and zoonotic malaria vector of the Buton Utara Wildlife Sanctuary, Southeast Sulawesi, Indonesia



The recent deforestation for agricultural, mining, and human re-settlement has significantly reduced the habitat of many non-human primates (NHPs) in Indonesia and intensifies interaction between the NHPs and humans and thus opening the possibility of pathogen spill-over. The emergence of zoonotic malaria, such as Plasmodium knowlesi, poses an immense threat to the current malaria control and elimination that aims for the global elimination of malaria by 2030. As malaria in humans and NHPs is transmitted by the female Anopheles mosquito, malaria vector control is very important to mitigate the spill-over of the malaria parasite to humans. The present study aims to explore the Anopheles species diversity in human settlements adjacent to the wildlife sanctuary forest in Buton Utara Regency, Southeast Sulawesi, Indonesia, and identify the species that potentially transmit the pathogen from monkey to human in the area.


Mosquito surveillance was conducted using larval and adult collection, and the collected mosquitoes were identified morphologically and molecularly using the barcoding markers, cytochrome oxidase subunit I (COI), and internal transcribed species 2 (ITS2) genes. Plasmodium sporozoite carriage was conducted on mosquitoes collected through human landing catch (HLC) and human-baited double net trap (HDNT).


The results revealed several Anopheles species, such as Anopheles flavirostris (16.6%), Anopheles sulawesi (3.3%), Anopheles maculatus (3.3%), Anopheles koliensis (1.2%), and Anopheles vagus (0.4%). Molecular analysis of the sporozoite carriage using the primate-specific malaria primers identified An. sulawesi, a member of the Leucosphyrus group, carrying Plasmodium inui sporozoite.


This study indicates that the transmission of zoonotic malaria in the area is possible and alerts to the need for mitigation efforts through a locally-tailored vector control intervention and NHPs habitat conservation.


Malaria still poses a major public health problem in 77 countries around the globe with the majority of death occurring in Africa. According to the global malaria report, there were 247 million malaria cases and 619,000 deaths in 2021 [1]. Despite significant progress in malaria incidence reduction within the last two decades in the Southeast Asian region, an immense emergence of zoonotic malaria poses a new challenge in areas where human malaria has been successfully eliminated such as Singapore, Malaysia, Indonesia, and other Southeast Asian countries [2]. Malaria is caused by a protozoan parasite of the genus Plasmodium and is transmitted to humans and other non-human primates (NHPs) by female Anopheles mosquitoes. Approximately 537 species of Anopheles are distributed in different parts of the world [3], but only about 30 species play a major role as malaria vectors [1, 4] as they prefer human blood meals to nurture their eggs [5].

Climate and environmental changes affect the presence of disease vectors that transmit many of the vector-borne diseases to humans and NHPs [6, 7]. Ecosystems are influenced by many factors, many of which are controlled by climate [8]. Disease vectors require certain ecosystems for their survival and reproduction [9]. The potential for malaria transmission is associated with water bodies acting as larval sites for malaria vectors, such as puddles, rivers, lakes, ponds, ground pools, and animal footpaths [10]. The wildlife sanctuary in Buton Utara Regency provides a suitable site or environment for the mosquito, including Anopheles to live as it also has some non-human primates, such as macaque and tarsius [11].

The archipelago of Indonesia is home for an approx. 80 Anopheles species, 26 of which have been incriminated as malaria vectors. Furthermore, ten species of Anopheles have been confirmed by molecular identification as malaria vectors, such as: Anopheles vagus, Anopheles barbirostris, Anopheles kochi, Anopheles nigerrimus, Anopheles tessellatus, Anopheles maculatus, Anopheles flavirostris, Anopheles aconitus, Anopheles karwari, and Anopheles peditaeniatus. Anopheles flavirostris is a species attracted to large mammals, including buffalos, cows, and people [12,13,14,15].

To achieve the goal of malaria elimination by 2030, efforts to control malaria through vector control are very important. Across Indonesia, there is huge variability in the diversity and distribution of malaria vectors; therefore, localized knowledge of malaria vectors is crucial to establishing suitable vector intervention strategies [16]. Studies conducted in various sites in Southeast Asia to date identified several vectors of zoonotic malaria, such as Anopheles dirus, Anopheles balabacensis and Anopheles leucosphyrus [17]. These mosquitoes are all classified as forest mosquitoes.

This study aims to identify the diversity of Anopheles spp. as well as the vectors of primate malaria in the Buton Utara Wildlife Sanctuary (BUWS) that potentially transmit zoonotic malaria to the adjacent human settlement.


Study area

The BUWS is located on Buton island precisely the area of Buton Utara Regency. It has four division resort areas: Resort I, Resort II, Resort III, and Resort Muna. Study areas were chosen based on the adjacent settlement area to a forest. The climate in this area is typically tropical, with a drier season proceeding from April to September and a wetter season from October to March. Mosquitoes were collected in larval and adult stages in Resort III of the wildlife sanctuary forest and adjacent human settlements in the villages: Labuan bajo, Lasiwa, Labajaya, and Laeya (Fig. 1). Local inhabitants in the four villages are subsistence farmers that grow corn, sweet potato, cashew nuts, and others. Cattle and chicken are reared, and some people enter the forest to look for rattan shoots and other edible flora.

Fig. 1
figure 1

Resort III of the wildlife sanctuary forest and adjacent human settlements in the villages: Labuan bajo, Lasiwa, Labajaya, and Laeya

Mosquito collections

Larval collection

The larval collection was conducted through water bodies surveillance in the forest adjacent to or within 500 m of the human settlement in 4 villages. The larval survey was conducted from the temporal house for crop protection in a forest. Villagers built this house to keep their farm in the forest, from pests such as wild boar and macaques. All potential aquatic larval habitats near the households and up to 500 m into the forest were sampled for mosquito larvae. First, different habitats were observed to determine if mosquito larvae were present. The habitats were rain pools, tree hollows, wells, springs, stream margins, and ditches. The habitats that were positive for larvae were included as sampling sites. Each sampling site was given a permanent number for repeated sampling. The sampling sites were georeferenced by Global Positioning System (GPS), and the coordinates (latitude and longitude), habitat type, and vegetation coverage were recorded. Distance from the nearest human settlement was measured by a tape measure. The water depth in each aquatic habitat was measured with a meter ruler at three spots and then the average was calculated. Salinity and pH were recorded using a Horiba Salt Laquatwin and Horiba pH Laquatwin. All physicochemical parameters were determined on-site only during the first sampling event.

Quantification of Anopheles larval density (D) was calculated as follows [18,19,20]:

$$D \left(per\;dipper\right)=\frac{total\;number\;of\;Anopheles\;larvae}{total\;number\;of\;dips}$$

Adult mosquitoes collection

Adult mosquito collections were conducted for 8 consecutive nights from the 15th of September until the 15th of October 2020. One each sequential night, a new area of Labuan resort was sampled. Mosquitoes were collected hourly from 7 p.m. to 7 a.m. On each night of sampling, mosquitoes were collected with one replicate of each of the following traps: indoor and outdoor human landing catch (HLC), human-baited double net trap (HDNT), animal-baited tent (ABT), and a CDC miniature light trap (CDC-LT). Although livestock for the ABTs was not universally available, this trap was only used in Labajaya (Hamlet 2) area. The biting mosquitoes were collected every hour from 7 p.m. to 7 a.m. by either the volunteer or entomology staff.

Human landing catch (HLC)

The HLC was performed by two male adult volunteers, one indoors and one outdoors. Each collector sat on a chair with the legs exposed from foot to knee and captured mosquitoes as soon as they land on the exposed legs before they commenced feeding using a flashlight and mouth aspirator [21, 22]. Each hour’s collection was kept separately in labelled paper cups. A supervisor was assigned to coordinate the activities and watch volunteers not fall asleep during the collection nights. All collectors were provided with anti-malaria prophylaxis to avoid a risk of contracting malaria during the collection period. Mosquitoes were identified as species the next morning. The human biting rate was calculated separately for indoor and outdoor catches. Man biting rate, mean hour density, and entomological inoculation rate will be calculated whenever applicable.

Human-baited double net trap (HDNT)

The HDNT in this study consisted of two box nets (inner and outer nets) with a roof made of canvas and set up in the study area except for Labuan Bajo (one place in Resort III). The inner net (97 cm high × 200 cm long × 100 cm wide) fully protects a human volunteer who rests on a mattress. The outer net (100 cm high × 250 cm long × 150 cm wide) is hung over the inner net and raised 30 cm off the ground. Mosquitoes attracted to human bait are collected every hour by research staff. The HDNT is an exposure-free tool since the lured mosquitoes are prevented to bite volunteers by the inner net. Outdoor mosquito sampling using the HDNT was conducted from 7 p.m. to 7 a.m. during each collection night.

Animal-baited tent (ABT)

An adult cow was loosely tied to a stake in the ground at the center of a large tent (Coleman 13 × 15 ft. screened canopy tent). ABT was conducted in once per 1 night. Mosquitoes attracted to bite and resting on the tent’s interior walls were collected hourly for 5 min using mouth aspirators.

CDC miniature light traps (CDC-LT)

A conventional CDC miniature light trap was also set outdoors at about 2 m from each of the selected trees in the forested areas, at the height of 1.5 m from the ground. The tree is between the breeding site and the house and where the CDC-LT was hung. The traps were set with a white LED light and an octanol lure (Octanol Emplura ® Merck). The captured mosquitoes were collected the following morning at 7 a.m. CDC-LT was run one unit per area.

Mosquito identification and molecular analyses

Genomic DNA was isolated from individual mosquito heads and thorax using the cetyltrimethylammonium bromide (CTAB)-based DNA extraction method. The ribosomal DNA internal transcribed spacer region (rDNA ITS2) was PCR-amplified using ITS2A and ITS2B primers [23,24,25] to identify Anopheles species. Two set primers (LCO-HCO and UEA 9.2-UEA 10.2) targeting mitochondrial DNA subunit 1 region were used in addition to confirming Leucosphyrus group identification [26, 27]. PCR products were visualized on a 2% agarose gel and purified by mixing 8 µl of the PCR product with 2U of exonuclease 1 (USB Corporation, Cleveland, OH), 1U of shrimp alkaline phosphatase (USB Corporation), and 1.8 µl of ddH20. This clean-up mixture was incubated at 37 °C for 15 min and then at 80 °C for 15 min to inactivate the enzymes. PCR products were sequenced and clean ITS2 sequences for each specimen were blasted using BLASTn against the NCBI GenBank database to confirm molecular species identification compared to voucher and published sequences.

Anopheles sporozoite carriage

Extracted DNA from each individual mosquito head and thorax was also used to test for Plasmodium infection using a nested PCR to amplify either a portion of the Plasmodium mitochondrial cytochrome oxidase subunit I gene (COI) or Plasmodium rDNA [28]. The Sanger sequences of these PCR-positive amplicons were assigned to Plasmodium species by comparing them to known Plasmodium voucher rDNA or COI sequences in the NCBI database.

Data analysis

Quantitative data was captured, managed, and analysed using Microsoft office excel basic functions and open-source software, RStudio version 2023.03.0 + 386 on R version 4.1.0 [29, 30]. Statistical analysis using Chi-Square and Pearson tests was used to correlate between the type of breeding site with the presence and density of Anopheles larvae. Tukey honest significant differences were used to check for multiple pairwise comparisons of larval densities. Statistical analysis for different methods used the Kruskall–Wallis test.


Mosquito larval sites

Surveys on water bodies at the four localities in the wildlife sanctuary identified 24 larval sites, from which 11 were positive for Anopheles larvae (Table 1; Fig. 2). The larval sites that contained Anopheles larvae were rain pools, springs, ditches, wells, and stream margins (Fig. 2). Ditches and river stream were the most frequently used larval sites, in which 88% contained Anopheles. Anopheles coexisted with non-Anopheles in ditches was 75%. Statistical analysis using the Chi-Square test showed no significant correlation between the type of breeding sites and the existence of Anopheles sp. larvae (P > 0.05) at 95% CI and between the type of breeding sites and Anopheles sp. larval density (P > 0.05) at 95% CI. The pH of the larval sites ranged from 8 to 10 whereas the salinity was 0 ppm, indicating that all had fresh water.

Table 1 Aquatic larval sites utilized by mosquito larvae in Buton Utara Wildlife Sanctuary
Fig. 2
figure 2

Types of aquatic larval habitats: a rain pool, b tree hollow, c well, d spring, e stream margin, and f ditch

Adult mosquito collection

The only location with cows and other livestock was Labajaya (Hamlet 2); therefore, all adult mosquito sampling methods were compared there. In the other three localities, the ABT method was not performed. A total of 246 adult mosquitoes were collected, primarily through HLC and HDNT, with the vast majority of the mosquitoes from genus Anopheles, Culex, and Aedes. The Anopheles and non-Anopheles collected at 4 localities for total of 7 collection nights is shown in Table 2. Of the 4 methods used, HLC and HDNT yielded highest mosquito catches. The CDC-LT and ABT were initially installed near the human settlement. This study found that the average mosquito number collected by HLC, HBDNT, and CDC-LT were 3.83, 12.8, and 2.57, respectively. The mosquito obtained from each method were statistically different (P < 0.05) at 95% CI. The density and diversity of mosquitoes collected at Labajaya were higher than in other areas.

Table 2 Total numbers of adult mosquitoes collected for localities at Resort BUTUR III from morphological identification

Mosquito species identification

Based on adult morphology, anophelines were roughly identified to 8 species of Anopheles, with An. flavirostris, Anopheles hackery, An. sulawesi as the most frequent species found (Table 2). Other species such as An. maculatus, An. vagus, Anopheles vanus, Anopheles subpictus, and Anopheles schuffneri were found mainly in the Labajaya area where human settlement is more densely built and more forest part was converted to agricultural land.

Molecular identification of the mosquitoes

Of the total 8 representative morphological species were PCR amplified and DNA sequenced using ITS2 fragment of the rDNA as a target, DNA sequence alignment revealed 5 distinct species, such as An. flavirostris, An. sulawesi, An. maculatus, An. subpictus, and An. vagus (Table 3). For the An. sulawesi, as the ITS2 fragment did not result in a conclusive species. Therefore, the COI fragment of mtDNA was amplified using 2 sets of primers LCO-HCO and UEA 9.2-UEA 10.2 as previously described. DNA sequence alignment of the UEA fragment of mtDNA successfully fished out the DNA sequence of mtDNA of An. sulawesi (GenBank ON908464; ON955533; DQ897967), with DNA sequence similarity of 97.01%, indicating that the species examined was An. sulawesi. Subsequently, the newly found ITS2 fragment of rDNA and the barcoding COI fragment was submitted to the GenBank as molecular markers for species identification (GenBank OP020395).

Table 3 Molecular identification using Internal Transcribed Spacer II gene (ITS2) at Resort BUTUR III

Vector incrimination for zoonotic malaria in BUWS

To determine the potential vectors of zoonotic malaria in BUWS, DNA was extracted from head-thoracal parts of all Anopheles species collected through HLC, and HDNT was amplified using rPLU primers in a nested PCR. Two mosquito samples were found to be positive for Plasmodium. DNA sequencing of the amplicon showed 2 samples of An. sulawesi carried DNA of Plasmodium inui. Based on the HLC and HDNT methods performed in 6 collection nights, the man-biting rate (MBR) and the entomological inoculation rate (EIR) of the An. sulawesi were calculated and as shown in Table 4. The finding confirmed that An. sulawesi is a vector for primate malaria and could potentially transmit zoonotic malaria as it bites humans during HLC and HDNT.

Table 4 Anopheles biting behavior in BUWS


Entomologic surveys at the BUWS revealed a high species diversity of Anopheles and non-Anopheles mosquitoes The findings indicate the potential spill-over of mosquito-borne pathogens from the NHPs to humans in the BUWS. In the context of zoonotic malaria, the previous finding of a high prevalence of Plasmodium spp. infection among the Macaque combined with relatively high Anopheles species diversity reflects a possible ransmission from NHPs to humans [31]. In support of this notion, An. sulawesi was found [32] to carry P. inui sporozoite during HLC with a relatively high sporozoite positivity rate and EIR. Although screening and treatment in the human population did not find any zoonotic malaria cases, the identification of one human malaria case of mixed Plasmodium falciparum and Plasmodium vivax may explain this phenomenon [33]. To date, zoonotic malaria cases have been mostly reported in areas where human malaria cases have been eliminated [34]. Other evidence also reveals that sporadic zoonotic malaria cases were reported on forest workers or forest goers [35].

The Anopheles mosquito vector that transmits zoonotic malaria in South and Southeast Asia is primarily of the leucosphyrus group, which is commonly named as forest Anopheles [36, 37]. The finding of An. sulawesi as a vector for primate malaria in BUWS further confirms the evidence as it is one of the members of the Hackery species complex of the Leucosphyrus group. Several other species of non-Leucosphyrus group, such as Anopheles letifer, An. kochi and An. sundaicus contributed less role in the transmission. In Southeast Asia, simian Plasmodium species are mainly transmitted by mosquitoes in the An. leucosphyrus and An. dirus complexes [17]. Whereas in BUWS particularly, An. sulawesi shows a relatively high man-biting rate and entomological inoculation rate for P. inui. Plasmodium knowlesi is currently the most common species of primate malaria that cause zoonotic infection. Other species, namely P. inui, Plasmodium cynomolgi, Plasmodium simiovale and Plasmodium coatneyi have been rarely reported. As the primate malaria species found in BUWS have all been reported to cause zoonotic infection in other areas, regular surveillance to monitor human infection is mandatory, particularly to human settlers at the fringe of BUWS. To prevent and reduce the potential for zoonotic malaria, the main efforts to reduce contact between mosquitoes and humans must be considered. Other efforts, especially through the conservation of the NHPs habitat and thus the forest in the BUWS should be sufficient to support food resources for their life. Therefore, illegal conversion of the forest into agricultural land by the human settlers in the area should be strictly prohibited as the BUWS has been assigned as a wildlife sanctuary by the government.

The environmental change caused by massive deforestation may result in climate change that affects the ecology of malaria vectors [38]. Various environmental conditions are valuable to model malaria transmission dynamics and the impact of climate and environmental change in Southeast Asia [39, 40]. In the context of BUWS, forest clearing for the agricultural land need to be controlled to prevent further impact on NHPs habitat and environment as the condition is certainly favorable for zoonotic malaria transmission.


The BUWS exhibits a high species diversity of mosquitoes, particularly Anopheles spp. The primate malaria prevalence among the endemic NHPs and the forest clearing for agricultural land may increase the pathogen spill-over to humans such as zoonotic malaria. Of the 5 Anopheles found, An. sulawesi carried P. inui and could potentially transmit primate malaria to human settlers. Knowledge of vector bionomics and behaviour is important for establishing zoonotic malaria mitigation efforts through vector intervention, NHPs habitat conservation, and regular malaria surveillance on human settlers.

Availability of data and materials

All relevant data are within the manuscript.


  1. WHO. World malaria report. Geneva: World Health Organization; 2022.

    Google Scholar 

  2. Lempang MEP, Dewayanti FK, Syahrani L, Permana DH, Malaka R, Asih PBS, et al. Primate malaria: an emerging challenge of zoonotic malaria in Indonesia. One Health. 2022;14: 100389.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Harbach RE. The phylogeny and classification of Anopheles. In: Manguin S, editor. Anopheles mosquitoes: new insights into malaria vectors. London: IntechOpen; 2013. p. 3–55.

    Google Scholar 

  4. WHO. World malaria report. Geneva: World Health Organization; 2019.

    Google Scholar 

  5. Marques J, Cardoso JCR, Felix RC, Santana RAG, Guerra MGB, Power D, et al. Fresh-blood-free diet for rearing malaria mosquito vectors. Sci Rep. 2018;8:17807.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Githeko AK, Lindsay SW, Confalonieri UE, Patz JA. Climate change and vector-borne diseases: a regional analysis. Bull World Health Organ. 2000;79:1136–47.

    Google Scholar 

  7. Dhimal M, Ahrens B, Kuch U. Species composition, seasonal occurrence, habitat preference and altitudinal distribution of malaria and other disease vectors in eastern Nepal. Parasit Vectors. 2014;7:540.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Montoya JM, Raffaelli D. Climate change, biotic interactions and ecosystem services. Philos Trans R Soc Lond B Biol Sci. 2010;365:2013–8.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Campbell-Lendrum D, Manga L, Bagayoko M, Sommerfeld J. Climate change and vector-borne diseases: what are the implications for public health research and policy? Philos Trans R Soc Lond B Biol Sci. 2015;370:20130552.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Minakawa N, Dida GO, Sonya GO, Futami K, Njenga SM. Malaria vectors in Lake Victoria and adjacent habitats in western Kenya. PLoS ONE. 2012;7:e32725.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Supriatna J, Shekelle M, Fuad HAH, Winbarni N, Dwiyahreni AA, Farid M, et al. Deforestation on the Indonesian island of Sulawesi and the loss of primate habitat. Glob Ecol Conserv. 2020;24:e01205.

    Article  Google Scholar 

  12. Mahdalena V, Wurisastuti T. Review of Anopheles species and its malaria vector status in East Nusa Tenggara, Papua and West Papua provinces (Gambaran distribusi spesies anopheles dan perannya sebagai vektor malaria di Provinsi Nusa Tenggara Timur, Papua dan Papua Barat). Spirakel. 2020;12:46–59.

    Article  Google Scholar 

  13. Wigati RA, Mardiana, Arianti Y, Mujiono. Incrimination of An. vagus Donitz, 1902 mosquito (Diptera: Culicidae) as a malaria vector in Kokap sub-district, Kulon Progo district, Yogyakarta. Ind J Health Sci. 2006;19:503–16.

    Google Scholar 

  14. Elyzar IRF, Sinka ME, Gething PW, Tarmizi SN, Surya A, Kusriastuti R, et al. The distribution and bionomics of Anopheles malaria vector mosquitoes in Indonesia. Adv Parasitol. 2013;83:173–266.

    Article  Google Scholar 

  15. Sugiarto, Kesumawati HU, Soviana S, Hakim L. Confirmation of Anopheles peditaeniatus and Anopheles sundaicus as malaria vectors (Diptera: Cullicidae) in Sungai Nyamuk village, Sebatik island Utara Kalimantan, Indonesia using enzyme-linked immunosorbent assay. J Med Entomol. 2016;53:1422–32.

    Article  Google Scholar 

  16. Rosenberg R, Lindsey NP, Ficher M, Gregory CJ, Hinckley AF, Mead PS, et al. Vital signs: trends in reported vector borne diseases cases-United States and Territories, 2004–2016. Morb Mortal Wkly Rep. 2018;67:496–501.

    Article  Google Scholar 

  17. van de Straat B, Sebayang B, Grigg MJ, Staunton K, Garjito TA, Vythilingam I, et al. Zoonotic malaria transmission and land use change in Southeast Asia: what is known about the vectors. Malar J. 2022;21:109.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Wang H, Wang Y, Cheng P, Wang H, Wang H, Liu H, et al. The larval density of mosquitoes (Diptera: Culicidae) in Jiaxiang County, Shandong Province, China: influence of bacterial diversity, richness, and physicochemical factors. Front Ecol Evol. 2021;9:82.

    Google Scholar 

  19. Villarreal-Treviño C, Penilla-Navarro RP, Vázquez-Martínez MG, Moo-Llanes DA, Ríos-Delgado JC, Fernández-Salas I, et al. Larval habitat characterization of Anopheles darlingi from its northernmost geographical distribution in Chiapas, Mexico. Malar J. 2015;14:517.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Mattah PAD, Futagbi G, Amekudzi LK, Mattah MM, De Souza DK, Kartey-Attipoe WD, et al. Diversity in breeding sites and distribution of Anopheles mosquitoes in selected urban areas of southern Ghana. Parasit Vectors. 2017;10:25.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Kenea O, Balkew M, Tekie H, Gebre-Michael T, Deressa W, Loha E, et al. Comparison of two adult mosquito sampling methods with human landing catches in south-central Ethiopia. Malar J. 2017;16:30.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Wong J, Bayoh N, Olang G, Killeen GF, Hamel MJ, Vulule JM, et al. Standardizing operational vector sampling techniques for measuring malaria transmission intensity: evaluation of six mosquito collection methods in western Kenya. Malar J. 2013;12: 143.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Zomuanpuiia R, Ringnghetia L, Brindhaa S, Gurusubramanianb G, Senthil KN. ITS2 characterization and Anopheles species identification of the subgenus Cellia. Acta Trop. 2013;125:309–19.

    Article  Google Scholar 

  24. Sum JE, Lee WC, Amir A, Braima KA, Jeffery J, Abdul-Aziz NM, et al. Phylogenetic study of six species of Anopheles mosquitoes in Peninsular Malaysia based on inter-transcribed spacer region 2 (ITS2) of ribosomal DNA. Parasit Vectors. 2014;7:309.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Walton C, Somboon P, O’Loughlin SM, Zhang S, Harbach RE, Linton YM, et al. Genetic diversity and molecular identification of mosquito species in the An. maculatus group using ITS2 region of rDNA. Infect Genet Evol. 2007;7:93–102.

    Article  CAS  PubMed  Google Scholar 

  26. Sallum MA, Foster PG, Li C, Sithiprasasna R, Wilkerson RC. Phylogeny of the leucosphyrus group of Anopheles (cellia) (Diptera: Culicidae) based on mitochondrial gene sequences. Ann Entomol Soc Am. 2007;100:27–35.

    Article  CAS  Google Scholar 

  27. 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.

    Article  CAS  PubMed  Google Scholar 

  28. 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 epidemiology studies. Am J Trop Med Hyg. 1999;60:687–92.

    Article  CAS  PubMed  Google Scholar 

  29. Posit team. RStudio: integrated development environment for R. Boston: Posit Software, PBC; 2023.

  30. R Core Team. R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing. 2021.

  31. Handayani S, Dewi RM. Plasmodium cynomolgi: Ancaman Zoonotik Malaria Kedua Setelah Plasmodium knowlesi di Era Eliminasi Malaria di Indonesia. J Biotek Medisiana Indonesia. 2019;8:105–13.

    Google Scholar 

  32. Warren M, Wharton RH. The vectors of simian malaria: identity, biology, and geographical distribution. J Parasitol. 1963;49:892–904.

    Article  CAS  PubMed  Google Scholar 

  33. Ramasamy R. Zoonotic malaria—global overview and research and policy needs. Front Public Health. 2014;2:123.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Akter R, Vythilingam I, Khaw LT, Qvist R, Lim YA, Sitam FT, et al. Simian malaria in wild macaques: first report from Hulu Selangor district, Selanngor Malaysia. Malar J. 2015;14:386.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Fahmi F, Pasaribu AP, Theodora M, Wangdi K. Spatial analysis to evaluate risk of malaria in Northern Sumatera, Indonesia. Malar J. 2022;21:241.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Maeno Y, Quang NT, Culleton R, Kawai S, Masuda G, Nakazawa S, et al. Humans frequently exposed to a range of non-human primate malaria parasite species through the bites of Anopheles dirus mosquitoes in South-central Vietnam. Parasit Vectors. 2015;8:376.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Maeno Y. Molecular epidemiology of mosquitoes for the transmission of forest malaria in South-central Vietnam. Trop Med Health. 2017;45:27.

    Article  PubMed  PubMed Central  Google Scholar 

  38. MacDonald AJ, Mordecai EA. Amazon deforestation drives malaria transmission, and malaria burden reduces forest clearing. Proc Natl Acad Sci USA. 2019;116:22212–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Wang X, Zhou G, Zhong D, Wang X, Ying W, Yan Z, et al. Life table studies revealed significant effect of deforestation on the development and survivorship of Anopheles minimus larvae. Parasit Vectors. 2016;9:323.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Roiz D, Ruiz S, Soriguer R, Figuerola J. Climatic effects on mosquito abundance in Mediterranean wetlands. Parasit Vectors. 2014;7: 333.

    Article  PubMed  PubMed Central  Google Scholar 

Download references


The authors are grateful to the Directorate General of Higher Education, Ministry of Education, Culture, Research and Technology, Republic of Indonesia for providing the “Pendidikan Magister menuju Doktor untuk Sarjana Unggul (PMDSU) Program” and express special thanks to all staff at the laboratory of malaria and vector resistance research group, the Eijkman Research Center for Molecular Biology, National Research and Innovation Agency (BRIN), Cibinong, Indonesia, for their support.


This project was funded by a research grant from the University of Hasanuddin.

Author information

Authors and Affiliations



DS and PA concepted the study and drafted the original draft. ML, DP, PA, and DS responsible to formal analysis and writing review. ML, DP, PA, SW, FD, IR, LP, WS, RW, LW, and DS responsible for investigation, contributed in data curation, and design methodology. ML, PA, and DS have the same contribution as main contributor. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Din Syafruddin.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Ethics Committee for Medical Research, Faculty of Medicine, University of Hasanuddin, Makassar, Indonesia, No. 56/UN4.

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.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Lempang, M.E.P., Permana, D.H., Asih, P.B.S. et al. Diversity of Anopheles species and zoonotic malaria vector of the Buton Utara Wildlife Sanctuary, Southeast Sulawesi, Indonesia. Malar J 22, 221 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: