Skip to main content

Effect of different habitat types on abundance and biting times of Anopheles balabacensis Baisas (Diptera: Culicidae) in Kudat district of Sabah, Malaysia

Abstract

Background

We investigated the effect of five common habitat types on the diversity and abundance of Anopheles spp. and on the biting rate and time of Anopheles balabacensis (currently the only known vector for Plasmodium knowlesi in Sabah) at Paradason village, Kudat, Sabah. The habitats were forest edge, playground area, longhouse, oil palm plantation and shrub-bushes area. Sampling of Anopheles was done monthly using the human landing catch method in all habitat types for 14 months (October 2013 to December 2014, excluding June 2014). The Anopheles species were morphologically identified and subjected to PCR assay for the detection of Plasmodium parasites. Generalised linear mixed models (GLMM) were applied to test the variation in abundance and biting rates of An. balabacensis in different habitat types.

Results

A total of 1599 Anopheles specimens were collected in the village, of which about 90% were An. balabacensis. Anopheles balabacensis was present throughout the year and was the dominant Anopheles species in all habitat types. The shrub bushes habitat had the highest Anopheles species diversity while forest edge had the greatest number of Anopheles individuals caught. GLMM analysis indicated that An. balabacensis abundance was not affected by the type of habitats, and it was more active during the early and late night compared to predawn and dawn. PCR assay showed that 1.61% of the tested An. balabacensis were positive for malaria parasites, most of which were caught in oil palm estates and infected with one to two Plasmodium species.

Conclusions

The identification of infected vectors in a range of habitats, including agricultural and farming areas, illustrates the potential for humans to be exposed to P. knowlesi outside forested areas. This finding contributes to a growing body of evidence implicating environmental changes due to deforestation, expansion of agricultural and farming areas, and development of human settlements near to forest fringes in the emergence of P. knowlesi in Sabah.

Background

Borneo Island, the third largest island in the world, is famous for its diverse and unique flora and fauna. Unfortunately, Borneo Island has been experiencing forest losses recently. From 1975 to 2015, approximately 18.5 million hectares (Mha) of old forests were cleared for industrial plantations [1]. In 2015, an estimated 7.9 Mha of the total Borneo landmass (73.7 Mha) were planted with oil palm. Changes in the environment due to deforestation, mainly for development, logging and agricultural plantations, are also changing the ecology of animals on the island, including anopheline ecology [2]. With continued deforestation across Borneo [3], resulting fragmentation of forest areas may increase the spillover risk of zoonotic, vector-borne diseases due to increased contact between humans and disease reservoirs and vectors [4].

Sabah, the second largest state in Malaysia, situated in northern Borneo Island with an estimated area of 73,620 km2 has recently experienced tremendous forest loss. Prior to the 1980s, about 80% of Sabah was covered with mainly dipterocarp forests [5]. However, by the mid-1980s, only 60% of Sabah was covered with forest, a loss of 20% in a short time due to logging and agriculture activities. By 1995, about 90% of the primary forests in Sabah had been lost, with the area reduced from around 2.8 to 0.3 Mha [5].

One of the current major mosquito borne-diseases in Sabah is zoonotic knowlesi malaria caused by Plasmodium knowlesi transmitted from long-tailed and pig-tailed macaques by Anopheles mosquitoes [6]. In 2016, 93% of the knowlesi malaria cases in Malaysia (1263/1357) were reported from Malaysian Borneo [7]. Deforestation and land alteration have been thought as the main reasons behind the increase of knowlesi malaria cases by creating new habitat types suitable for the Anopheles vector [8].

In Southeast Asia, Anopheles species are present in high diversity, including species complexes that occur in sympatry but have different behavior [9, 10]. The high diversity of these vectors and the complexity of the species could lead to misidentification which may influence the success of vector control programmes [11, 12]. In East Malaysia, An. balabacensis and An. latens have been confirmed as the main vectors for P. knowlesi in Sabah and Sarawak, respectively [13, 14]. In Peninsular Malaysia, several species of Anopheles belonging to the Leucosphyrus group had been identified as the natural vectors for P. knowlesi [15]. Nevertheless, other Anopheles sp. (e.g. An. kochi) can act as vectors for P. knowlesi but further studies are needed to establish their role in transmission [16].

In Asia, P. knowlesi was probably the first simian malaria confirmed to cause widespread human infection after a group of villagers had been found infected naturally in Kapit, Sarawak, East Malaysia in 2004 [17]. This species was previously misidentified as P. falciparum or P. malariae [18] until molecular techniques were introduced. Since then, P. knowlesi cases have also been reported in Sabah [19,20,21] and Peninsular Malaysia [22,23,24]. Similarly, many other countries in Southeast Asia have also reported cases, e.g. Kalimantan, Indonesia [25, 26]; Sumatra, Indonesia [27]; Brunei [28]; Thailand [29, 30]; Vietnam [31, 32]; Singapore [16, 33]; Phillippines [34]; Myanmar [35]; Cambodia [36]; and recently Laos [37].

In this study, we conducted an entomological investigation to determine the diversity and abundance of Anopheles species, especially An. balabacensis, in the five types of habitats that are now commonly found in the rural areas of Sabah. Anopheles balabacensis has been known as the dominant species in Kudat Division, mainly in Banggi Island located off the northern coast of Sabah near Marudu Bay [38]. Current changes in the landscape of Kudat district of Kudat Division have provided new habitats for An. balabacensis closer to human settlements, have altered its distribution, increased their opportunity of feeding on humans and infecting them with malaria parasites.

Methods

Sampling site

Kudat Division of Sabah is the smallest division, and covers the districts of Kudat, Pitas, Kota Marudu, Balambangan, Banggi Island and several other islands. The climate in Sabah is tropical with an average temperature of 32 °C for lowland areas and 21 °C for highland areas, typically hot and sunny the whole year. Rainfall is common throughout the year and influenced by northeast (November to March) and southwest monsoons (May to September).

This study was carried out in Kudat district (Fig. 1) which receives a high volume of rain during the northeast monsoon but less rain during the southwest monsoon. The population of Kudat district was about 85,300 in July 2010 with the majority of Rungus ethnicity [39]. They live in several hamlets each consisting of several houses or longhouses. The villagers practice swidden farming, growing rice, corn and vegetables for their own consumption. Some families might have their own small plantation of coconut, rubber or oil palm but many work as laborers in oil palm plantations.

Fig. 1
figure 1

Map showing the location of Kudat district, location of Paradson village and the five types of habitat (or sampling sites) labeled as FE, PG, LH, OP and BU. Abbreviations: FE, forest edge between forest and deforested area; PG, playground area surrounded with small trees including bamboo, coconut and rubber trees; LH, area outside a longhouse; OP, an oil palm plantation; BU, shrub-bushes area with small and tall grasses near to a coconut plantation

Sampling of Anopheles mosquitoes was conducted in Paradason village [40] (6°46′11.75″N, 116°47′3.23″E) which had 28 households and 122 residents in October 2013 census (Fig. 1). Prior to and within the study period five knowlesi malaria cases had been reported from this village. Long-tailed macaques were also sighted in the surrounding areas. The village has an area of about 1.98 km2 comprising forest, plantation, farm and human settlement areas. Five types of habitats were selected: a forest edge between secondary forest and deforested area (FE); a playground area (or an open space area) surrounded with small trees including bamboo, coconut and rubber trees where the village youths gather during evening to play football (PG); an area outside a longhouse (LH); a farm which was converted into an oil palm plantation (OP); and a shrub-bushes area with small and tall grasses near to a coconut plantation (BU) (Fig. 2).

Fig. 2
figure 2

The habitat types where mosquitoes were sampled are: a FE, forest edge between forest and deforested area; b PG, playground area surrounded with small trees including bamboo, coconut and rubber trees; c LH, area outside a longhouse; d OP, an oil palm plantation; and e BU, shrub-bushes area with small and tall grasses near to a coconut plantation

Anopheles mosquitoes were collected by the human landing catch (HLC) method, with two persons in a team, which involved exposing lower part of their legs to attract mosquitoes. Mosquitoes that landed on their legs were captured using a plastic specimen tube (2 cm diameter × 6 cm length). Each captured mosquito was placed in a separate specimen tube which was labeled with the time of capture. The mosquitoes were sampled from 18:00–06:00 h (12 h). Sampling was done once a month for 14 months starting in October 2013 and ending in December 2014 with no sampling in June 2014. Meteorological data such as temperature, humidity and rainfall were obtained from Kudat meteorological station which was located 17 km from Paradason village.

Morphological identification of Anopheles mosquitoes

The Anopheles mosquitoes were killed by keeping them in a freezer (− 20 °C) and each individual was identified based on morphological characteristics using published keys [41,42,43,44]. Anopheles specimens were mounted onto Nu poly strip using ultra-thin micro-headless pin and identified under a compound microscope. Each individual mosquito was kept separately in a clean microfuge tube and transported to the Faculty of Medicine & Health Sciences, Universiti Malaysia, Sabah, for further processing.

Species composition, diversity and evenness

The diversity index of Anopheles mosquitoes in each habitat type was determined using the Shannon–Wiener index, \(H\) [45]. This index is commonly used to characterize species diversity in a community which accounts for both abundance and evenness of the species present. Evenness assumes a value between 0 and 1 with 1 being complete evenness. The index is calculated based on the formula, \(H = - \mathop \sum \nolimits_{i = 1}^{S} P_{i} \left( {{\text{ln }}P_{i} } \right)\) where \(S\) is the number of species and \(P_{\text{i}}\) is the proportion of individuals of ith species.

Extraction of total DNA from Anopheles mosquitoes

Each Anopheles individual was placed separately in a clean mortar and the tissue was homogenized using a sterile pestle. The total DNA was extracted from each individual using the DTAB-CTAB method [46] with some modifications. First, 600 µl of DTAB solution was added into the mortar and the tissue was ground using pestle until homogenized. Then, the homogenized tissue was transferred into a clean 1.5 ml microcentrifuge tube and incubated at 68 °C for 30 min. Subsequently, 700 µl of chloroform was added into the microcentrifuge tube which was inverted ten times to mix the contents and centrifuged at 13,000× rpm for 5 min. Then, 400 µl of the upper aqueous layer was carefully transferred into a new clean 1.5 ml microcentrifuge tube and mixed with 900 µl of sterile dH2O and 100 µl CTAB solution by gently inverting the microcentrifuge tube several times. The tube was then allowed it to sit at room temperature for 5 min and spun at 13,000× rpm for 10 min. The supernatant was discarded and the DNA pellet was re-suspended in 300 µl of 1.2 M NaCl solution. Total DNA was precipitated by adding 750 µl of absolute ethanol and centrifuged at 13,000× rpm for 5 min. Then, the supernatant was discarded and the DNA was washed with 500 µl of 70% ethanol and centrifuged at 13,000× rpm for 2 min. The DNA pellet was air-dried at 45 °C for 10 min, re-suspended in 30 µl of Tris-EDTA (pH8.0) buffer and stored at − 30 °C until use.

Detection of malaria parasites

The presence of malaria parasites in the mosquitoes was detected using nested PCR assay by targeting the small subunit ribosomal RNA (SSU rRNA) gene of Plasmodium. The PCR primer pair rPLU1 and rPLU5 were used in first PCR reaction, while rPLU3 and rPLU4 were used in the second reaction [47]. For internal control, another nested PCR reaction was performed separately at the same time to amplify the cytochrome c oxidase subunit 2 (cox2) gene of Anopheles [48]. When a specimen was found positive for malaria parasites, further identification of Plasmodium sp. was conducted using nine species-specific primers (Additional file 1: Table S1). For positive control, known positive DNA of each Plasmodium species was obtained from the Parasitology and Medical Entomology Laboratory, Faculty of Medicine and Health Sciences, Universiti Malaysia Sabah and tested with the nine species-specific PCR primers (Additional file 2: Figure S1).

Both PCR genus-specific and species-specific reactions were performed in a final volume of 25.0 µl. Reaction components were prepared by mixing 5.0 µl of 5× PCR buffer (Promega, Madison, WI, USA), 0.5 µl of (10 mM) dNTPs (Promega), 3.0 µl of (25 mM) MgCl2, 1.0 µl of (10 µM) forward and reverse primers, 0.3 µl of (5.0 U/µl) Taq DNA polymerase (Promega), 2.0 µl of DNA template and sterile dH2O up to 25.0 μl final volume. After completion of the first PCR, 2.0 µl of the PCR products was used as a DNA template in the second PCR. The reaction was carried out using a thermal cycler (T100™ Thermal Cycler, Bio-Rad, California, CA, USA) with the following conditions: initial denaturation at 95 °C for 5 min; 35 cycles of denaturation at 94 °C for 1 min, annealing for 1 min and extension at 72 °C for 1 min; final extension at 72 °C for 5 min. The annealing temperature used was based on the optimal temperature for each set of primers (Additional file 1: Table S1). PCR products were analyzed on 1.5% agarose gel electrophoresis stained with RedSafe™ nucleic acid staining solution (iNtRON Biotechnology, Jungwon-Seongnam, South Korea) and visualized using UV transilluminator.

Analysis of relationship between An. balabacensis and meteorological data

The correlation between the number of An. balabacensis caught per person per night with temperature, relative humidity and rainfall from Kudat meteorological station was examined using Pearson’s linear correlation coefficient (r) and checked for significance. The test was conducted using SPSS v.18 for Windows 10.

Analysis of abundance of An. balabacensis

Statistical analysis of mosquito abundance was conducted using R programming software for statistical analysis (v.3.2.2). Generalised linear mixed models (GLMM) were constructed to test the variation in abundance of An. balabacensis related to month of collection, types of habitat and time of biting. This model included both fixed and random predictor variables in the analysis. Time of biting was categorized into four periods for the total catch of An. balabacensis: 18:00–21:00 h (early night); 21:00–00:00 h (late night); 0:00–03:00 h (predawn); and 03:00–06:00 h (dawn). In the analysis, the type of habitat and the time of biting were considered as fixed effects. Month was fitted alternatively as fixed (to predict monthly value) and random effect (to test differences between types of habitat and time of biting while controlling variation in seasonality). In order to determine the best model, both negative binomial and Poisson distributions including zero inflation were tested in the model. Tukey’s post-hoc test was used to compare predicted means between the fixed effects and the interaction between the two fixed effects.

Results

Anopheles species composition, diversity and evenness

A total of ten Anopheles species were recorded from Paradason village (Table 1). The highest number of species was recorded from shrub bushes (BU, 10 species) and the least was from oil palm plantation (OP, 4 species). Anopheles balabacensis was the dominant species in all types of habitat with the total number caught reaching 90% (1437/1599). The highest number of Anopheles was recorded from the forest edge (FE, 21.58%) followed by oil palm (OP, 20.39%) and shrub bushes (BU, 20.39%), longhouse (LH, 19.07%) and playground (PG, 18.57%). Overall, the diversity of Anopheles mosquito in Paradason village is considered low as indicated by the evenness value of 0.2.

Table 1 Total numbers of Anopheles species caught from each type of habitat from October 2013 to December 2014

Abundance of An. balabacensis by month and habitat type

A total of 1437 individuals of An. balabacensis were caught for the 70 sampling nights pooling all the habitats at Paradason village, equivalent to 10.3 ± 0.47 bites per man per night. The fluctuations in the monthly number of An. balabacensis caught in the five types of habitat showed similar trends (Fig. 3). The GLMM predicted the highest mean monthly number of An. balabacensis for February (15.0 ± 3.06 bites per man per night) and the lowest for August (5.40 ± 1.52 bites per man per night). Fitting the model with a Poisson distribution and comparing the predicted mean between all months showed that there was a significant difference only between December 2013 and August 2014 (12.0 ± 1.55 bites per man per night versus 5.40 ± 1.52 bites per man per night, Z = 3.44, P = 0.0372) and between February 2014 and August 2014 (15.0 ± 3.06 bites per man per night versus 5.40 ± 1.52 bites per man per night, Z = − 3.59, P = 0.0231).

Fig. 3
figure 3

Number of An. balabacensis caught in the five types of habitat in each month from October 2013 to December 2014. FE, forest edge between forest and deforested area; PG, playground area surrounded with small trees including bamboo, coconut and rubber trees; LH, area outside a longhouse; OP, an oil palm plantation; and BU, shrub-bushes area with small and tall grasses near to a coconut plantation

The mean bites per man per night for each habitat type predicted by GLMM were 10.7 ± 1.67, 10.1 ± 1.58, 9.8 ± 1.03, 9.3 ± 1.50 and 9.0 ± 1.50 for forest edge, oil palm estates, bushes, longhouse and playground, respectively. However, using a Poisson distribution in the model, no difference was detected in the number caught between the different habitat types.

During the study period, the range of the mean monthly temperature, relative humidity and rainfall in Kudat district were 25.9–28.4 °C, 80.4–87.5% RH and 0–27.8 mm, respectively (Fig. 4). Analysis of the number of An. balabacensis caught on each sampling night against the environmental factors of the day showed significant correlation between the numbers of An. balabacensis caught with temperature (r = − 0.243, P = 0.043), but not with relative humidity (r = 0.220, P = 0.067) and rainfall (r = 0.204, P = 0.090) (Fig. 5).

Fig. 4
figure 4

Monthly mean temperature, relative humidity and rainfall and the number of An. balabacensis collected in Kudat district from October 2013 to December 2014

Fig. 5
figure 5

Plot of An. balabacensis (bites/man/night) caught on each sampling night (n = 70) against the environmental factors: mean daily temperature; mean daily relative humidity; and mean daily rainfall

Biting time of An. balabacensis

Overall, biting times of An. balabacensis from different types of habitat was similar with a peak biting time of 19:00–20:00 h in PG, LH, OP and BU and 20:00–21:00 h in FE (Fig. 6). More An. balabacensis were caught within the first three hours after dusk (18:00–21:00 h) with the number decreasing gradually until dawn. Only one mosquito bite per person or less in each hour after midnight was recorded. Fitting the model with a Poisson distribution, there were significant differences between the catch at early night versus predawn, early night versus dawn, late night versus predawn, late night versus dawn, and predawn versus dawn (Table 2).

Fig. 6
figure 6

Mean number of bites/man/night (and 95% standard error bars) of An. balabacensis recorded in each hour for the various types of habitat

Table 2 Mean biting rate for each time period of catch as predicted by GLMM

Analysis of interaction between habitat type and time-period showed that only three interactions were significant, namely FE × early night versus FE × dawn, BU × early night versus BU × dawn, and BU × late night versus BU× dawn (Additional file 3: Table S2).

Presence of malaria parasites in Anopheles specimens

A total of 1586 Anopheles specimens were tested for malaria parasites using a PCR assay, and only 23 (1.45%) were found positive. These were all An. balabacensis, constituting 1.6% or 23/1425 of the species collected. Most infected specimens were caught in the oil palm plantation and shrub bushes, and during early and late night. These An. balabacensis were infected with one or two simian malaria species with P. inui being the most common species, while P. knowlesi was found in two mosquitoes only (Table 3; representative figures of agarose gels are shown in Additional file 4: Figure S2).

Table 3 Number of Plasmodium-infected An. balabacensis caught from each habitat type, the recorded biting time and the Plasmodium species as detected using PCR assay

Discussion

We conducted a survey at Paradason village, in Kudat district, Sabah to investigate if the diversity and abundance of adult Anopheles is influenced by the habitat types that are generally found in the rural areas of Sabah, namely forest edge, playground, longhouse, oil palm plantation and shrub bushes.

We collected ten Anopheles species from all the habitats, among which An. balabacensis was the dominant species in all types of habitat. The dominance of this species in Kudat Division was also noted previously [14], although our study in Tajau Laut village located further north of Kudat indicated that An. balabacensis and An. tessellatus were almost equally common [48]. However, differences in the dominant Anopheles species have been reported in different locations across the wider state of Sabah, possibly due to heterogeneities in the surrounding environments. For example, in the northern area of Ranau district, An. maculatus was recorded as the dominant species, followed by An. balabacensis and An. donaldi [49], while in Kinabatangan district An. donaldi was the dominant species replacing An. balabacensis, possibly due to the development and malaria control activities conducted there [50]. Our recent data from Paus village in southern part of Ranau district showed that the dominant species is An. donaldi followed by An. balabacensis and An. barbumbrosus (unpublished data). Further south in Keningau district, in Keritan Ulu village, An. maculatus, An. balabacensis and An. barbumbrosus were almost equally abundant (Additional file 5: Table S3). Overall, the diversity of Anopheles species in Paradason village is considered low as indicated by the evenness index due to overwhelming abundance of An. balabacensis in all types of habitat.

Anopheles balabacensis has been known as the dominant and common Anopheles species in Kudat Division of Sabah since 1980s [38] and as an efficient human malaria vector in Sabah since 1950s [51, 52]. Now, almost 40 years later, this species still ranks high among the dominant Anopheles species in Kudat district despite significant land use and land cover changes, due to development for human settlements, infrastructure, farming and plantations [8, 53]. Our findings suggest that An. balabacensis is adaptable and able to thrive in the man-made environment, with relatively high densities in plantation and farm areas. This is evident by the presence of An. balabacensis breeding sites observed in oil palm plantation near to human settlements. Similar inference was drawn for An. cracens in Peninsular Malaysia [54] and An. dirus (s.l.) in Thailand [55], both of which have adapted and colonized the forest edge and plantations after deforestation.

The shrub-bush habitat recorded the highest number of Anopheles species. This could be a result of shrub and bushes providing ample refuge areas for resting before and after blood meals [56], and shade for the breeding sites in both permanent and temporary water bodies. In the Paradason area, this bush habitat was observed in close proximity to human settlements, farms, plantations and deforested areas, which could be a source of infective Anopheles. Furthermore, a troop of long-tailed macaques was also observed in study site especially in the morning (06:00–07:00 h) and evening (16:00–18:00 h) coming from a nearby forest area into a coconut farm to forage for food. These monkeys were positive for P. knowlesi [57], just as those in other parts of Malaysia which recorded 50–97% positive [58,59,60].

Although forest activities have previously been associated with increased human risk of P. knowlesi within this area [61], GLMM analysis indicated that there was no difference in number of An. balabacensis caught at various habitat types, indicating individuals have a risk of being exposed to an infected mosquito across habitat types. Although the HLC was conducted on different nights at each habitat type which may introduce bias between sampling, An. balabacensis is known to have a short flight distance with an average maximum flight distance of less than 500 meters [62]. This suggests An. balabacensis is likely to stay in close proximity to breeding sites, as observed with An. farauti (s.l.) [10] and conducting HLCs on different nights between habitat types is not likely to substantially impact results of this study.

Data on the biting times indicate that An. balabacensis starts to be actively biting as soon as it gets dark, especially at the forest edge or bush areas that have been deforested, highlighting the need for vector control measures targeting these times when people are unlikely to be sleeping under a bednet. A previous study [63] also showed more An. balabacensis bite humans in logged forest as compared to primary and virgin reserve forests. This study did not identify a significant difference in abundance of An. balabacensis between early and late nights; however, human exposure may be higher earlier in the night when people are more likely to visit these environments. A previous study in Kudat demonstrated the risk of getting bitten by An. balabacensis was higher for those who stay outside the house after dusk (18:00–00:00 h) as compared to those who stay inside [64]. The number of Anopheles mosquitoes seeking blood meals at night might be dependent on natural preferences, host availability and environmental conditions such as temperature, relative humidity, rainfall, moonlight and wind speed [65,66,67]. Anopheles balabacensis in this study area showed a significant negative correlation with temperature but did not show a significant correlation with the increment of relative humidity and rainfall. The lack of correlation between the mosquito catch and rainfall could be due to the localized variations in rainfall and microclimate as data was only available from the Kudat meteorological station, located 17 kilometers from the sampling site at Paradason village.

Only 1.45% of all Anopheles specimens were found infected with simian Plasmodium and of these only An. balabacensis was positive. However, since whole mosquitoes were processed, infectious (oocyst positive) and infective (salivary glands positive) individuals cannot be differentiated. Most of the infections were by single species, P. inui (47.8%), followed by P. cynomolgi (17.4%). In addition, the two mosquitoes infected with P. knowlesi were also co-infected either with P. inui or P. cynomolgi. These two specimens were caught at the forest edge and shrub bushes habitats, respectively, where a troop of long-tailed macaques had been observed.

Anopheles balabacensis is a vector for P. knowlesi as well as for other simian malaria parasites [14, 64]. Furthermore, we could not detect any human malaria parasites by PCR in the Anopheles specimens although P. falciparum and P. vivax cases have been reported in Kudat district [19, 68]. This contrasts with findings in Vietnam in which An. dirus (s.l.), a primary vector for P. knowlesi, had co-infections with P. vivax and P. falciparum [31, 69]. The majority of the positive An. balabacensis was infected with P. inui, another species of simian malaria parasite that has a high possibility of becoming zoonotic. Although there has not yet been a report of naturally acquired P. inui in humans, experimentally it has been shown that humans are susceptible to this species and will develop fever if infected [70]. Another simian Plasmodium that is of great concern is P. cynomolgi as this has been reported infecting humans under natural conditions [71, 72]. Furthermore, the real number of P. cynomolgi infection in humans might be underestimated as this species could be misidentified as P. vivax by microscopy [71].

Although there was no difference in vector densities, molecular results indicated that a higher proportion (65.2%) of infected An. balabacensis was caught earlier in the evening between 18:00 h and midnight. Although these numbers are limited, this suggests people may be at a higher risk of getting infected with malaria during this period, especially when they are still outside their houses.

Conclusions

The dominant Anopheles species in Kudat district is An. balabacensis which could be found equally abundant in different types of habitats, suggesting people are at a risk of exposure to infected vectors across a range of environments. Further work is needed to characterise the breeding sites within these different habitats. While limited numbers of infected An. balabacensis were caught, the presence of infected vectors in oil palm estates demonstrates the need for vector control measures targeting these habitats. Furthermore, mixed infections with P. knowlesi and other zoonotic malaria species illustrates the potential for spillover into human populations and the need for further surveillance. While deforestation and agricultural expansion can create economic opportunities, these changes may also increase spatial overlap between human, macaque and vector populations in new habitat types, presenting a challenge for malaria elimination.

Availability of data and materials

The datasets supporting the conclusions of this article are included within the article and its additional files.

Abbreviations

PCR:

polymerase chain reaction

GLMM:

generalised linear mixed model

MHa:

million hectares

HLC:

human landing catch

DNA:

deoxyribonucleic acid

DTAB:

dodecyltrimethylammonium bromide

CTAB:

cetyltrimethylammonium bromide

SSU rRNA:

small subunit ribosomal RNA

cox2:

cytochrome c oxidase subunit II

References

  1. Gaveau DLA, Sheil D, Husnayaen A, Salim MA, Arjasakusuma S, Ancrenaz M, et al. Rapid conversions and avoided deforestation: examining four decades of industrial plantation expansion in Borneo. Sci Rep. 2016;6:32017.

    Article  CAS  Google Scholar 

  2. Yasuoka J, Levins R. Impact of deforestation and agricultural development on anopheline ecology and malaria epidemiology. Am J Trop Med Hyg. 2007;76:450–60.

    Article  Google Scholar 

  3. Gaveau DLA, Sloan S, Molidena E, Yaen H, Sheil D, Abram NK, et al. Four decades of forest persistence, clearance and logging on Borneo. PLoS ONE. 2014;9:e101654.

    Article  Google Scholar 

  4. Imai N, White MT, Ghani AC, Drakeley CJ. Transmission and control of Plasmodium knowlesi: a mathematical modelling study. PLoS Negl Trop Dis. 2014;8:e2978.

    Article  Google Scholar 

  5. Gunggut H, Saufi K, Zaaba Z, Liu MSM. Where have all the forests gone Deforestation in land below the wind. Procedia Soc Behav Sci. 2014;153:363–9.

    Article  Google Scholar 

  6. Millar SB, Cox-Singh J. Human infections with Plasmodium knowlesi—zoonotic malaria. Clin Microbiol Infect. 2015;21:640–8.

    Article  CAS  Google Scholar 

  7. World Health Organization Malaria Policy Advisory Committee. Outcomes from the evidence review group on Plasmodium knowlesi. 2017. http://www.who.int/malaria/mpac/mpac-mar2017-plasmodium-knowlesi-presentation. Accessed 7 Apr 2018.

  8. Fornace KM, Abidin TR, Alexander N, Brock P, Grigg MJ, Murphy A, et al. Association between landscape factors and spatial patterns of Plasmodium knowlesi infections in Sabah, Malaysia. Emerg Infect Dis. 2016;22:201–8.

    Article  CAS  Google Scholar 

  9. Morgan K, Somboon P, Walton C. Understanding Anopheles diversity in Southeast Asia. Rijeka: InTech; 2013. p. 327–55.

    Google Scholar 

  10. Sinka ME, Bangs MJ, Manguin S, Chareonviriyaphap T, Patil AP, Temperley WH, et al. The dominant Anopheles vectors of human malaria in the Asia-Pacific region: occurrence data, distribution maps and bionomic précis. Parasit Vectors. 2011;4:89.

    Article  Google Scholar 

  11. Bortel WV, Harbach RE, Trung HD, Roelants P, Backeljau T, Coosemans M. Confirmation of Anopheles varuna in Vietnam, previously misidentified and mistargeted as the malaria vector Anopheles minimus. Am J Trop Med Hyg. 2001;65:729–32.

    Article  Google Scholar 

  12. Saeung A. Anopheles (Diptera: Culicidae) species complex in Thailand: identification, distribution, bionomics and malaria-vector importance. Int J Parasitol Res. 2012;4:75–82.

    Article  Google Scholar 

  13. Vythilingam I, Tan CH, Asmad M, Chan ST, Lee KS, Singh B. Natural transmission of Plasmodium knowlesi to humans by Anopheles latens in Sarawak, Malaysia. Trans R Soc Trop Med Hyg. 2006;100:1087–8.

    Article  CAS  Google Scholar 

  14. Wong ML, Chua TH, Leong CS, Khaw LT, Fornace K, Wan-Sulaiman WY, et al. Seasonal and spatial dynamics of the primary vector of Plasmodium knowlesi within a major transmission focus in Sabah, Malaysia. PLoS Negl Trop Dis. 2015;9:e0004135.

    Article  Google Scholar 

  15. Vythilingam I, Wong ML, Wan-Yussof WS. Current status of Plasmodium knowlesi vectors: a public health concern? Parasitology. 2018;145:32–40.

    Article  CAS  Google Scholar 

  16. Jeslyn WPS, Huat TC, Vernon L, Irene LMZ, Sung LK, Jarrod LP, et al. Molecular epidemiological investigation of Plasmodium knowlesi in humans and macaques in Singapore. Vector Borne Zoonotic Dis. 2011;11:131–5.

    Article  Google Scholar 

  17. Singh B, Sung LK, Matusop A, Radhakrishnan A, Shamsul SSG, Cox-Singh J, et al. A large focus of naturally acquired Plasmodium knowlesi infections in human beings. Lancet. 2004;363:1017–24.

    Article  Google Scholar 

  18. Lee KS, Cox-Singh J, Singh B. Morphological features and differential counts of Plasmodium knowlesi parasites in naturally acquired human infections. Malar J. 2009;8:73.

    Article  Google Scholar 

  19. Barber BE, William T, Dhararah P, Anderios F, Grigg M, Yeo TW, et al. Epidemiology of Plasmodium knowlesi malaria in north-east Sabah, Malaysia: family clusters and wide age distributions. Malar J. 2012;11:401.

    Article  Google Scholar 

  20. Barber BE, William T, Grigg MJ, Menon J, Auburn S, Marfurt J, et al. A prospective comparative study of knowlesi, falciparum, and vivax malaria in Sabah, Malaysia: high proportion with severe disease from Plasmodium knowlesi and Plasmodium vivax but no mortality with early referral and artesunate therapy. Clin Infect Dis. 2013;56:383–97.

    Article  CAS  Google Scholar 

  21. Goh XT, Lim YA, Vythilingam I, Chew CH, Lee PC, Ngui R, et al. Increased detection of Plasmodium knowlesi in Sandakan division, Sabah as revealed by PlasmoNexâ„¢. Malar J. 2013;12:264.

    Article  Google Scholar 

  22. Jiram AI, Hisam S, Reuben H, Husin SZ, Roslan A, Ismail WRW. Submicroscopic evidence of the simian malaria parasite, Plasmodium knowlesi, in an orang asli community. SE Asian J Trop Med Public Health. 2016;47:591–9.

    Google Scholar 

  23. Lee WC, Chin PW, Lau YL, Chin LC, Fong MY, Yap CJ, et al. Hyperparasitaemic human Plasmodium knowlesi infection with atypical morphology in peninsular Malaysia. Malar J. 2013;12:88.

    Article  Google Scholar 

  24. Yusof R, Lau YL, Mahmud R, Fong MY, Jelip J, Ngian HU, et al. High proportion of knowlesi malaria in recent malaria cases in Malaysia. Malar J. 2014;13:168.

    Article  Google Scholar 

  25. Figtree M, Lee R, Bain L, Kennedy T, Mackertich S, Urban M, et al. Plasmodium knowlesi in human, Indonesian Borneo. Emerging. Infect Dis. 2010;16:672–4.

    Article  Google Scholar 

  26. Setiadi W, Sudoyo H, Trimarsanto H, Sihite BA, Saragih RJ, Juliawaty R, et al. A zoonotic human infection with simian malaria, Plasmodium knowlesi, in central Kalimantan. Indonesia. Malar J. 2016;15:218.

    Article  Google Scholar 

  27. Lubis IND, Wijaya H, Lubis M, Lubis CP, Divis PCS, Beshir KB, et al. Contribution of Plasmodium knowlesi to multispecies human malaria infections in North Sumatera, Indonesia. J Infect Dis. 2017;215:1148–55.

    Article  CAS  Google Scholar 

  28. Ramaswami A, Pisharam JK, Aung H, Ghazala K, Maboud K, Chong VH, et al. Co-incidental Plasmodium knowlesi and mucormycosis infections presenting with acute kidney injury and lower gastrointestinal bleeding. Am J Case Rep. 2013;14:103–5.

    Article  Google Scholar 

  29. Jongwutiwes S, Putaporntip C, Iwasaki T, Sata T, Kanbara H. Naturally acquired Plasmodium knowlesi malaria in human, Thailand. Emerg Infect Dis. 2004;10:2211–3.

    Article  Google Scholar 

  30. Jongwutiwes S, Buppan P, Kosuvin R, Seethamchai S, Pattanawong U, Sirichaisinthop J, et al. Plasmodium knowlesi malaria in humans and macaques, Thailand. Emerg Infect Dis. 2011;17:1799–806.

    Article  Google Scholar 

  31. Marchand RP, Culleton R, Maeno Y, Quang NT, Nakazawa S. Co-infections of Plasmodium knowlesi, P. falciparum, and P. vivax among humans and Anopheles dirus mosquitoes, Southern Vietnam. Emerg Infect Dis. 2011;17:1232–8.

    Article  Google Scholar 

  32. Van den Eede P, Van HN, Overmeir CV, Vythilingam I, Duc TN, Hung LX, et al. Human Plasmodium knowlesi infections in young children in central Vietnam. Malar J. 2009;8:249.

    Article  Google Scholar 

  33. Ng OT, Ooi EE, Lee CC, Lee PJ, Ng LC, Wong PS, et al. Naturally acquired human Plasmodium knowlesi infection, Singapore. Emerg Infect Dis. 2008;14:814–6.

    Article  CAS  Google Scholar 

  34. Luchavez J, Espino F, Curameng P, Espina R, Bell D, Chiodini P, et al. Human infections with Plasmodium knowlesi, the Philippines. Emerg Infect Dis. 2008;14:811–3.

    Article  Google Scholar 

  35. Jiang N, Chang Q, Sun X, Lu H, Yin J, Zhang Z, et al. Co-infections with Plasmodium knowlesi and other malaria parasites, Myanmar. Emerg Infect Dis. 2010;16:1476–8.

    Article  Google Scholar 

  36. Khim N, Siv S, Kim S, Mueller T, Fleischmann E, Singh B, et al. Plasmodium knowlesi infection in humans, Cambodia, 2007–2010. Emerg Infect Dis. 2011;17:1900–2.

    Article  Google Scholar 

  37. Iwagami M, Nakatsu M, Khattignavong P, Soundala P, Lorphachan L, Keomalaphet S, et al. First case of human infection with Plasmodium knowlesi in Laos. PLoS Negl Trop Dis. 2018;12:e0006244.

    Article  Google Scholar 

  38. Hii JLK, Kan S, Vun YS, Chin KF, Tambakau S, Chan MKC, et al. Transmission dynamics and estimates of malaria vectorial capacity for Anopheles balabacensis and An. flavirostris (Diptera: Culicidae) on Banggi island, Sabah, Malaysia. Ann Trop Med Parasitol. 1988;82:91–101.

    Article  CAS  Google Scholar 

  39. Sabah State Department of Health. 2016 Annual Report of Unit for Control of Vector Borne Diseases (RKPBV), Division of Public Health, Sabah, Malaysia; 2016 (Jabatan Kesihatan Negeri Sabah. Laporan Tahunan RKPBV 2016. Unit Rancangan Kawalan Penyakit Bawaan Vektor, Bahagian Kesihatan Awam, JKNS; 2016).

  40. Fornace KM, Herman LS, Abidin TR, Chua TH, Daim S, Lorenzo PJ, et al. Exposure and infection to Plasmodium knowlesi in case study communities in northern Sabah, Malaysia and Palawan, The Philippines. PLoS Negl Trop Dis. 2018;12:e0006432.

    Article  Google Scholar 

  41. Cagampang-Ramos A, Darsie Jr RF. Illustrated keys to the Anopheles mosquitoes of the Philippine Islands. Malaria Eradication Training Center, Technical Report 70-1, USAF Fifth Epidemiological Flight. Manila: PACAF; 1970.

  42. Rattanarithikul R, Harrison B, Harbach RE, Panthusiri P, Coleman RE. Illustrated keys to the mosquitoes of Thailand IV. Anopheles. Southeast Asian J Trop Med Public Health. 2006;37(2):1–128.

    PubMed  Google Scholar 

  43. Reid JA. Anopheline mosquitoes of Malaya and Borneo. Kuala Lumpur: Government of Malaysia; 1968.

    Google Scholar 

  44. Sallum MAM, Peyton EL, Harrison BA, Wilkerson RC. Revision of the Leucosphyrus group of Anopheles (Cellia) (Diptera, Culicidae). Rev Bras Entomol. 2005;49:1–152.

    Article  Google Scholar 

  45. Shannon CE, Weaver W. The mathematical theory of communication. Urbana: University of Illinois Press; 1963.

    Google Scholar 

  46. Phillips A, Simon C. Simple, efficient, and nondestructive DNA extraction protocols for arthropods. Ann Entomol Soc Am. 1995;88:281–3.

    Article  CAS  Google Scholar 

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

    Article  CAS  Google Scholar 

  48. Hawkes F, Manin BO, Ng SH, Torr SJ, Drakeley C, Chua TH, Ferguson HM. Evaluation of electric nets as means to sample mosquito vectors host seeking on humans and primates. Parasit Vectors. 2017;10:338.

    Article  Google Scholar 

  49. Rohani A, Chan ST, Abdullah AG, Tanrang H, Lee HL. Species composition of mosquito fauna in Ranau, Sabah, Malaysia. Trop Biomed. 2008;25:232–6.

    CAS  PubMed  Google Scholar 

  50. Vythilingam I, Chan ST, Shanmugratnam C, Tanrang H, Chooi KH. The impact of development and malaria control activities on its vectors in the Kinabatangan area of Sabah, East Malaysia. Acta Trop. 2005;96:24–30.

    Article  CAS  Google Scholar 

  51. Hii JLK, Vun YS. A study of dispersal, survival and adult population estimates of the malaria vector, Anopheles balabacensis Baisas (Diptera: Culicidae) in Sabah, Malaysia. Trop Biomed. 1985;2:121–31.

    Google Scholar 

  52. McArthur. The control of Anopheles leucospyrus. Malaria conference for the Western Pacific and Southeast Asia Regions, Taipei, 15–27 November 1954.

  53. Bryan JE, Shearman PL, Asner GP, Knapp DE, Aoro G, Lokes B. Extreme differences in forest degradation in Borneo: comparing practices in Sarawak, Sabah, and Brunei. PLoS ONE. 2013;8:e69679.

    Article  CAS  Google Scholar 

  54. Jiram AI, Vythilingam I, NoorAzian YM, Yusof YM, Azahari AH, Fong MY. Entomologic investigation of Plasmodium knowlesi vectors in Kuala Lipis, Pahang, Malaysia. Malar J. 2012;11:213.

    Article  Google Scholar 

  55. Rosenberg R, Andre RG, Somchit L. Highly efficient dry season transmission of malaria in Thailand. Trans R Soc Trop Med Hyg. 1990;84:22–8.

    Article  CAS  Google Scholar 

  56. Dewald JR, Fuller DO, Müller GC, Beier JC. A novel method for mapping village-scale outdoor resting microhabitats of the primary African malaria vector, Anopheles gambiae. Malar J. 2016;15:489.

    Article  Google Scholar 

  57. Stark DJ, Fornace KM, Brock PM, Abidin TR, Gilhooly L, Jalius C, et al. Long-tailed macaque response to deforestation in a Plasmodium knowlesi-endemic area. EcoHealth. 2019. https://doi.org/10.1007/s10393-019-01403-9.

    Article  PubMed  Google Scholar 

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

    Article  Google Scholar 

  59. Lee K, Divis PCS, Zakaria SK, Matusop A, Julin RA, Conway DJ, et al. Plasmodium knowlesi: reservoir hosts and tracking the emergence in humans and macaques. PLoS Pathog. 2011;7:e1002015.

    Article  CAS  Google Scholar 

  60. Vythilingam I, NoorAzian YM, Huat TC, Jiram AI, Yusri YM, Azahari AH, et al. Plasmodium knowlesi in humans, macaques and mosquitoes in peninsular Malaysia. Parasit Vectors. 2008;1:26.

    Article  Google Scholar 

  61. Grigg MJ, Cox J, William T, Jelip J, Fornace KM, Brock PM, et al. Individual-level factors associated with the risk of acquiring human Plasmodium knowlesi malaria in Malaysia: a case-control study. Lancet Planet Health. 2017;1:e97–104.

    Article  Google Scholar 

  62. Verdonschot PFM, Besse-Lototskaya AA. Flight distance of mosquitoes (Culicidae): a metadata analysis to support the management of barrier zones around rewetted and newly constructed wetlands. Limnologica. 2014;45:69–79.

    Article  Google Scholar 

  63. Brant HL, Ewers RM, Vythilingam I, Drakeley C, Benedick S. Vertical stratification of adult mosquitoes (Diptera: Culicidae) within a tropical rainforest in Sabah, Malaysia. Malar J. 2016;15:370.

    Article  Google Scholar 

  64. Manin BO, Ferguson HM, Vythilingam I, Fornace K, William T, Torr SJ, et al. Investigating the contribution of peri-domestic transmission to risk of zoonotic malaria infection in humans. PLoS Negl Trop Dis. 2016;10:e0005064.

    Article  Google Scholar 

  65. Basseri H, Raeisi A, Khakha MJ, Pakarai A, Abdolghafar H. Seasonal abundance and host-feeding patterns of anopheline vectors in malaria endemic area of Iran. J Parasitol Res. 2010. https://doi.org/10.1155/2010/671291.

    Article  PubMed  PubMed Central  Google Scholar 

  66. Freire MG, Schweigmann NJ, Svagelj WS, Loetti MV, Jensen O, Burroni NE. Relationship between environmental conditions and host-seeking activity of Ochlerotatus albifasciatus (Diptera: Culicidae) in an agroecosystem and in an urban area in Chubut, central Patagonia, Argentina. J Nat Hist. 2016;50:1369–80.

    Article  Google Scholar 

  67. Stone C, Chitnis N, Gross K. Environmental influences on mosquito foraging and integrated vector management can delay the evolution of behavioral resistance. Evol Appl. 2015. https://doi.org/10.1111/eva.12354.

    Article  PubMed  PubMed Central  Google Scholar 

  68. William T, Jelip J, Menon MY, Anderios F, Mohammad R, Mohammad TAAM, et al. Changing epidemiology of malaria in Sabah, Malaysia: increasing incidence of Plasmodium knowlesi. Malar J. 2014;13:390.

    Article  Google Scholar 

  69. Nakazawa S, Marchand RP, Quang NT, Culleton R, Manh ND, Maeno Y. Anopheles dirus co-infection with human and monkey malaria parasites in Vietnam. Int J Parasitol. 2009;39:1533–7.

    Article  Google Scholar 

  70. Coatney GR, Chin W, Contacos PG, King HK. Plasmodium inui, a quartan-type malaria parasite of old world monkeys transmissible to man. J Parasitol. 1966;52:660–3.

    Article  CAS  Google Scholar 

  71. Imwong M, Madmanee W, Suwannasin K, Kunasol C, Peto TJ, Tripura R, et al. Asymptomatic natural human infections with the simian malaria parasites Plasmodium cynomolgi and Plasmodium knowlesi. J Infect Dis. 2018. https://doi.org/10.1093/infdis/jiy519.

    Article  PubMed Central  Google Scholar 

  72. Ta TH, Hisam S, Lanza M, Jiram AI, Ismail N, Rubio JM. First case of naturally acquired human infection with Plasmodium cynomolgi. Malar J. 2014;13:68.

    Article  Google Scholar 

Download references

Acknowledgments

The authors would like to acknowledge the Universiti Malaysia Sabah for all the research facilities provided, Mr Fazreen and Mr Nemran for helping in the fieldwork. We would also like to acknowledge the community of Paradason for their cooperation and warm hospitality during this study.

Funding

This study was funded by the Medical Research Council, Natural Environment Research Council, Economic and Social Research Council, and Biotechnology and Biosciences Research Council through the Environmental and Social Ecology of Human Infectious Diseases Initiative, grant no. G1100796.

Author information

Authors and Affiliations

Authors

Contributions

THC, IV, KF and CJD designed the conceptualization of the study. CJD acquired the funding. THC, BOM, IV and CJD designed the experiments. BOM and THC implemented the experiments and collected the data. BOM identified the specimens and carried out molecular analysis. THC and BOM carried out statistical analysis. BOM and THC drafted the manuscript. All authors contributed to drafting the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Tock H. Chua.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the National Medical Research Register of the Malaysian Ministry of Health (NMRR, Reference No. NMRR-12-786-13048). Consent to carry out mosquito collection was obtained from the village council or the village headman and the land owners. All volunteers who carried out mosquito collections signed informed consent forms and were provided with antimalarial prophylaxis during the study period.

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 files

Additional file 1: Table S1.

PCR primers used for detecting Plasmodium spp. in Anopheles specimens.

Additional file 2: Figure S1.

Amplification of known Plasmodium DNA using Plasmodium genus and species-specific PCR primers. Lanes M: 1.0 kb DNA ladder (Promega); Lane a: genus Plasmodium, 240 bp (rPLU3 + rPLU4); Lane b: P. coatneyi, 504 bp (PctF1 + PctR1); Lane c: P. inui, 479 bp (PinF2 + INAR3); Lane d: P. fieldi, 421 bp (PfldF1 + PfldR2); Lane e: P. cynomolgi, 137 bp (CY2F + CY4R); Lane f: P. knowlesi, 424 bp (PkF1140 + PkR1550); Lane g: P. falciparum, 370 bp (NewPLFshort + FARshort); Lane h: P. vivax, 476 bp (NewPLFshort + VIRshort); Lane i: P. malariae, 241 bp (NewPLFshort + MARshort); Lane j: P. ovale, 407 bp (NewPLFshort + OVRshort); Lane k: negative control (no DNA template).

Additional file 3: Table S2.

Mean interaction statistic between habitat type and time period as predicted by GLMM. Habitat types: FE: forest edge; PG: playground; LH: longhouse; OP: oil palm plantation; BU: shrub bushes. Time periods are: early night (18:00–21:00 h); late night (21:00–00:00 h); predawn (00:00–03:00 h); dawn (03:00–06:00 h).

Additional file 4: Figure S2.

Detection and identification of Plasmodium in Anopheles. Top: Anopheles specimens PD171-PD183 using Plasmodium genus-specific PCR primers. Internal control targeting cox2 gene of Anopheles was used. The numbers represent replicate 1 and 2. Lanes M: 1.0 kb DNA ladder; J: infected Anopheles with Plasmodium; K: non-infected Anopheles; L: negative control (no DNA template). Bottom: Anopheles specimen PD178 using nine species-specific PCR primers. Lanes M: 1.0 kb DNA ladder; a: genus Plasmodium; b: P. coatneyi; c: P. inui; d: P. fieldi; e: P. cynomolgi; f: P. knowlesi; g: P. falciparum; h: P. vivax; i: P. malariae; j: P. ovale; k: negative control.

Additional file 5: Table S3.

Total Anopheles specimens collected at study sites in Paus, Ranau and Keritan Ulu, Keningau.

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

Chua, T.H., Manin, B.O., Vythilingam, I. et al. Effect of different habitat types on abundance and biting times of Anopheles balabacensis Baisas (Diptera: Culicidae) in Kudat district of Sabah, Malaysia. Parasites Vectors 12, 364 (2019). https://doi.org/10.1186/s13071-019-3627-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13071-019-3627-0

Keywords