Skip to main content

Vertical stratification of adult mosquitoes (Diptera: Culicidae) within a tropical rainforest in Sabah, Malaysia

Abstract

Background

Malaria cases caused by Plasmodium knowlesi, a simian parasite naturally found in long-tailed and pig-tailed macaques, are increasing rapidly in Sabah, Malaysia. One hypothesis is that this increase is associated with changes in land use. A study was carried out to identify the anopheline vectors present in different forest types and to observe the human landing behaviour of mosquitoes.

Methods

Mosquito collections were carried out using human landing catches at ground and canopy levels in the Tawau Division of Sabah. Collections were conducted along an anthropogenic disturbance gradient (primary forest, lightly logged virgin jungle reserve and salvage logged forest) between 18:00 and 22:00 h.

Results

Anopheles balabacensis, a vector of P. knowlesi, was the predominant species in all collection areas, accounting for 70 % of the total catch, with a peak landing time of 18:30–20:00 h. Anopheles balabacensis had a preference for landing on humans at ground level compared to the canopy (p < 0.0001). A greater abundance of mosquitoes were landing in the logged forest compared to the primary forest (p < 0.0001). There was no difference between mosquito abundance in the logged forest and lightly logged forest (p = 0.554). A higher evening temperature (p < 0.0001) and rainfall (p < 0.0001) significantly decreased mosquito abundance during collection nights.

Conclusions

This study demonstrates the potential ability of An. balabacensis to transmit P. knowlesi between canopy-dwelling simian hosts and ground-dwelling humans, and that forest disturbance increases the abundance of this disease vector. These results, in combination with regional patterns of land use change, may partly explain the rapid rise in P. knowlesi cases in Sabah. This study provides essential data on anthropophily for the principal vector of P. knowlesi which is important for the planning of vector control strategies.

Background

Malaria still remains a public health problem throughout tropical and sub-tropical regions of the world, with an estimated 214 million new cases and 438,000 deaths in 2015 [1]. Four Plasmodium species are recognised as causing human malaria; Plasmodium falciparum, Plasmodium malariae, Plasmodium ovale and Plasmodium vivax, but recently a fifth species, Plasmodium knowlesi, has been recognized as causing symptomatic malaria in humans [2–4]. Plasmodium knowlesi, transmitted by the forest-dwelling Anopheles from the Leucosphyrus group, is an emerging cause for zoonotic human malaria in Southeast Asia [3, 5–7]. Malaysia has had a successful malaria control programme, aimed to eliminate malaria by 2020, with marked reductions in reported cases of P. falciparum and P. vivax, but there has been a recent increase in P. knowlesi cases [8–11]. Plasmodium knowlesi is now the most common cause of malaria in the Malaysian state of Sabah, accounting for 62 % of all malaria incidences in 2013 and presenting a threat to malaria elimination [8].

It has been proposed that land use change, including deforestation, forest fragmentation and agricultural practices, has increased the incidence of P. knowlesi by increasing the encroachment of humans into previously forested areas, allowing a higher interaction between vectors and human and macaque hosts [9, 12]. The increase in P. knowlesi cases may also be underestimated due to misdiagnosis during microscopic examination [13, 14]. Microscopy of stained blood smears allows differentiation between species, but frequent misdiagnosis occurs in areas containing P. falciparum, P. vivax and P. knowlesi [14].

Monkeys, particularly the long-tailed macaque (Macaca fascicularis) and the pig-tailed macaque (Macaca nemestrina) found in Southeast Asia, are the two main natural hosts of P. knowlesi [15]. A study in Sabah showed nearly all patients with P. knowlesi malaria had a recent history of forest or forest-edge exposure, and had seen a monkey in the preceding month [16]. Most members of the Leucosphyrus group, from the genus Anopheles, feed primarily on monkeys in the canopy and are capable of transmitting various Plasmodium species [17]. Anopheles balabacensis is the predominant vector of human malaria in Sabah [18, 19], and has also been incriminated as a P. knowlesi vector [20, 21].

Most primates are arboreal. Although some species of chimpanzees, baboons and macaques rest and feed at ground level during the day, primates almost always sleep in the canopy during the night [22]. As Anopheles species generally bite between 6 p.m. and 6 a.m., primate roosting sites are potentially a key location for disease transmission between primate hosts. It is hypothesized that vectors are biting humans at ground level, but if given the opportunity, will bite at canopy level. This study expects that key vector species should be present in the disturbed forest habitats where people come into contact with monkeys, but it is unknown which species should expect to be present in primary forest.

Since transmission is increasing in Sabah, it is important to identify the P. knowlesi vectors present and understand their biting behaviour within forest habitats. While other vector assessments are ongoing in the Interior, West Coast, Kudat and Sandakan Divisions in Sabah, the Tawau Division has not been studied. This study was conducted to determine the vertical distribution of mosquitoes and their biting preference in Sabah, Malaysia.

Methods

Study site

The study was conducted in the Tawau Division of Sabah, Malaysia. Three areas were selected along a forest disturbance gradient; primary lowland dipterocarp rainforest (PF), virgin jungle reserve (VJR) and twice-logged disturbed dipterocarp rainforest (LF) with forest disturbance quantified using on-the-ground forest plots [23]. Primary forest survey points were selected in the vicinity of Danum Valley Field Centre (4°58′N, 117°42′E), located within the Danum Valley Conservation Area. This area consists of 43,800 ha of protected dipterocarp rainforest [24]. Virgin jungle reserve (4°40′N, 117°32′E) and logged forest survey points (4°41′N, 117°34′E) were selected within the Benta Wawasan oil palm plantation. The 45,601 ha area is a mixture of twice-logged rainforest, virgin jungle reserve, acacia and oil palm. The VJR, of 2200 ha, has been logged around the edge, but never logged in the steep interior [25]. Survey points were selected 500–1000 m from the VJR edge in locations that had undergone light logging. Logged forest survey points were in selectively twice-logged forest, logged during 1970s, 1990s–2000s resulting in the cumulative removal of ~180 m3 ha−1 of timber [26], and currently being further disturbed by additional salvage logging activity in surrounding areas. Further details of the project area are given by Ewers et al. [25].

Three survey points, with a minimum separation distance of 500 m, were selected in each area. One tree was selected at each point based on its accessibility into the canopy, low density of epiphytes and height. Visual tree assessments were carried out to make sure every tree was safe to climb. The trees selected ranged from a height of 15 m in the logged forest to 30 m in the virgin jungle reserve and primary forest. Survey points were a subset of those designated as part of the central sampling design of the ‘Stability of Altered Forest Ecosystems (SAFE) Project’, a large-scale fragmentation experiment, which is investigating the long-term effects of forest fragmentation [25].

Data collection

All data collection was carried out from April to July 2014, and mosquitoes were collected using human landing catches at ground and canopy height between 18:00 and 22:00 h. Four consecutive nights of collections were carried out in PF and VJR, and five nights in LF, using a rotation of collectors. Access was gained into the canopy using line insertion to high branches [27], followed by the double rope climbing technique taught by Canopy Access Limited [28]. Canopy samples were collected at a height of two-thirds the average canopy height at that location (10–20 m). Ground and canopy collections were conducted simultaneously. The average canopy height surrounding each selected tree was calculated using a laser rangefinder.

The collectors, with the aid of a red torch-light, aspirated mosquitoes off their own legs. Collected mosquitoes were placed into cups covered with a net cloth, and a new cup was used during every half an hour of collection. Mosquitoes were taken back to the field laboratory to be killed and sorted into individual tubes with silica gel. All mosquitoes were identified morphologically using keys [29–35].

Meteorological data

Air temperature (°C) and relative humidity (%) were measured at each site using microclimate dataloggers at the base and in the canopy of each tree during sampling. Nightly rainfall (mm) data were obtained from the nearest rain gauges located at SAFE Project and Danum Valley Field Centre for the duration of the field survey (March to July 2014; range 1.0–3.5 km from sample sites). In addition, lunar illumination (%), cloud cover and unusual climatic events (e.g. strong winds) were recorded every half an hour by the collectors.

Data analysis

Analyses were performed using R version 3.1.1 [36]. Simpson and Shannon diversity indices were calculated in each area using the vegan function ‘diversity’ (R package vegan, ‘diversity’) [37]. Species accumulation curves were calculated in each area using the vegan function ‘specaccum’ (R package vegan, ‘specaccum’) [37]. To estimate the number of undetected species and add them to the observed richness, true richness estimators were used. These estimators included the Chao species estimator (Chao 1) and Abundance Coverage Estimator (ACE). Chao 1 and ACE were used to estimate the extrapolated species richness in each area using ‘chao1’ and ‘ACE’ functions in the R package ‘fossil’ (R package fossil, ‘chao1’, ‘ACE’) [38].

The effect of canopy height (a dichotomous variable representing ground or canopy level) and forest disturbance (PF, VJR and LF) on mosquito abundance was analysed using a generalized linear mixed-effect model (R package lme4) [39], using day and site as random factors, with Poisson error distribution. Chi squared tests were used to compare the relative abundance of vector and non-vector species in each area, and between ground and canopy level. Differences in community composition at ground and canopy height were explored using Detrended Correspondence Analysis (DCA, vegan function ‘decorana’) [37]. The number of mosquitoes landing per night was used as a measure of relative abundance. Significant differences were tested for in community composition using a linear model with the first DCA axis as the response variable against canopy height and forest disturbance.

Results

Mosquito abundance

A total of 807 mosquitoes were collected from 39 human landing catch nights, consisting of 743 (92.1 %) anophelines, and 64 (7.9 %) culicines. A total of 555 (68.8 %) mosquitoes from 21 species were found at ground level in comparison to 252 (31.2 %) mosquitoes from 10 species at canopy level. Anopheles balabacensis was the predominant species at ground (62.4–87.5 % of all individuals) and canopy level (43.4–73.5 %) at each collection site. A full list of species is given in Table 1.

Table 1 Mosquitoes collected from different collection sites in the district of Tawau, Sabah, Malaysia

The number of species collected, Shannon index, Simpson index, Chao1 and ACE varied across forest disturbance and between ground and canopy level (Tables 2, 3). The Chao 1 and ACE predicted a higher number of species in the logged forest than the primary forest, and a higher number of species at ground level (Tables 2, 3). These patterns were also seen with the number of species collected, Shannon index and Simpson index. The species accumulation curves did not reach an asymptote at ground and canopy level, indicating not all species of mosquitoes had been collected (see Additional file 1).

Table 2 Mean species richness and diversity indices (±SE) of mosquito communities, collected at ground level
Table 3 Mean species richness and diversity indices (±SE) of mosquito communities, collected at canopy level

Effect of height and forest disturbance on mosquito abundance

Mosquito abundance in the canopy was significantly lower than at ground level (χ2 = 81.89, df = 1, p < 0.0001) (Table 4; Fig. 1a). Logged forest had a higher abundance than virgin jungle reserve or primary forest (χ2 = 10.94, df = 2, p < 0.05) (Fig. 1a). Similar patterns were seen with abundance of An. balabacensis (Table 4; Fig. 1b). Rainfall and a higher evening temperature decreased mosquito abundance (χ2 = 14.55, df = 1, p < 0.0001), whereas moonlight increased the abundance (χ2 = 4.49, df = 1, p < 0.05) (Table 4). Collector identity had no effect on mosquito abundance (χ2 = 3.871, df = 2, p = 0.144). Peak biting of An. balabacensis during the collection period was observed between 19:00 and 20:00 h in logged forest and virgin jungle reserve (Fig. 2).

Table 4 Effects of parameters on mosquito abundance
Fig. 1
figure 1

Effects of collection height on the human landing rate across a forest disturbance gradient. Effects of collection height on the human landing rate (number of mosquitoes per night per bait) across a forest disturbance gradient: primary forest (PF), lightly logged virgin jungle reserve (VJR), and twice-logged forest (LF). a Total abundance of all species combined, b Abundance of the most common species, Anopheles balabacensis, alone. Error bars show ± SE of the mean

Fig. 2
figure 2

Hourly number of Anopheles balabacensis landing per person per night. Hourly number of Anopheles balabacensis landing per person per night, at ground and canopy level, across an anthropogenic disturbance gradient from a Primary forest, b Virgin jungle reserve and c Logged forest. Error bars show ± SE of the mean

Community composition

The DCA plot indicated that the community composition of the canopy was significantly different to ground level collection and across forest disturbance, with the community from the two height strata being strongly separated on the first axis (F1, 60 = 24.72, p < 0.0001, Fig. 3). The communities in the three categories of forest disturbance were not separated along the first axis (F2, 60 = 0.92, p = 0.4), but the interaction between height and forest disturbance was significant (F2, 60 = 8.37, p < 0.001). The first two axes accounted for 66.9 % of the total variance. Anopheles balabacensis was prevalent at both ground and canopy level, but the ground level community also included species such as Anopheles barbirostris, Armigeres jugraensis, Culex vishnui and the Anopheles Aitkenii group, which were not present in the canopy. The species turnover that causes the difference in community composition between ground and canopy was driven by species that have not been reported to cause disease transmission.

Fig. 3
figure 3

Detrended correspondence analysis for adult mosquito abundance at ground level and in the canopy. Detrended correspondence analysis (DCA) plot showing the major axes of variation for adult mosquito abundance at ground level and in the canopy of a tropical rainforest. The two axes represent linear summaries of the variation in the species numbers and areas

Discussion

Studying the anthropophily of simian malaria vectors in the canopy of tropical forests is essential because the hosts are primarily arboreal. In order to develop, sustain or adapt a good control programme, it is important to monitor mosquito populations as well as their hosts and host-seeking preference, distribution and behaviour. Although previous studies in Southeast Asia have used monkey-baited traps at different canopy heights [12, 40, 41], this is the first study to attempt human landing catches, using this method, in the canopy. This study found that there was a higher abundance and human landing rate of mosquitoes at ground level, where people tend to be, than in the canopy where the simian hosts reside. This trend was driven by An. balabacensis, a key malaria vector in Sabah, and highlights the potential importance of this species in transmitting Plasmodium species from simian to human hosts.

Anopheles balabacensis was the most abundant mosquito in all sampled areas, accounting for 70 % of all collected species. Anopheles balabacensis is considered the most important vector of human malaria parasites on Banggi Island and mainland Sabah, Malaysia [18, 19, 30, 42]. In Sabah, An. balabacensis was found to be mainly exophagic, but could also be endophagic and exophilic [30, 43]. These behaviours cause An. balabacensis to be an effective vector of P. knowlesi from human to primate hosts. There were also two distinct subpopulations, one more zoophilic and one more anthropophilic [44–46]. Anopheles balabacensis occurs in forested areas, and readily bites human and monkey hosts, making it an ideal vector of simian malaria [20, 47].

Currently insecticide-treated bed nets and indoor residual spraying are the two main control methods in Malaysia [48]. This study showed An. balabacensis bites as early as 18:00 h in the Tawau Division. Other studies have shown the species bites as early as 18:00 h in recent years in comparison to late night biters in previous decades [44, 49–51]. Given that An. balabacensis is early evening biting, highly anthropophilic, exophagic and exophilic, current control methods are not sufficient to break the transmission cycle of P. knowlesi [52]. In Vietnam and Cambodia, long-lasting insecticidal hammocks (LLIH) were shown to reduce malaria incidence and prevalence in forested areas, and may prove to be an additional effective tool in reduction of malaria in Malaysia [52–54]. The use of repellents have been used for malaria control, but need to be tested in forest and plantations areas.

This study found a different community composition of mosquitoes in the canopy to that at ground level. Different mosquito species have particular flight distributions, with certain species flying and feeding close to ground, some species showing a preference for higher canopy heights, while others show a random distribution [55, 56]. The percentage biting at different canopy heights can be affected by microclimate conditions, such as relative humidity, temperature, wind speed and rainfall [57, 58], but may also change according to time of day [59].

Moonlight appeared to have a significant impact on mosquito activity, with human landing rates increasing on bright nights. Although some studies have shown moonlight increases relative abundance of biting vectors [60–64], others have shown a decrease [65–68] or no effect at all [69]. Collection bias was reduced in this study by collecting in each area under different phases of the moon.

This study also showed how forest disturbance affected mosquito abundance, species richness and human landing rates. Vector abundance was greater in the lightly modified virgin jungle reserve and heavily modified logged forest than in the unmodified primary forest. These results may be explained by the availability of larval breeding sites. Wheel tracks in logged areas due to logging activities can provide breeding sites for a range of mosquito species, whereas wheel tracks are not present within primary forests or virgin jungle reserves [31]. Species richness, estimated by the Chao1 index and ACE, differed across forest disturbance and height, with logged forest and ground level having a higher species richness than primary forest and canopy.

Conclusions

This study has given an overview of mosquito species found in the Tawau Division, including the anthropophily of P. knowlesi vectors at canopy and ground levels. Anopheles balabacensis was the predominant species found in primary forest, virgin jungle reserve and logged forest with a preference for landing on humans at ground level. As An. balabacensis is a vector of human and simian malaria, these findings will be useful for the planning of control strategies of malaria vectors.

Abbreviations

ACE:

abundance coverage estimator

DCA:

detrended correspondence analysis

LF:

logged forest

LLIH:

long-lasting insecticidal hammocks

PF:

primary forest

VJR:

virgin jungle reserve

References

  1. WHO. World malaria report. Geneva: World Health Organization; 2015 (2015).

    Google Scholar 

  2. White NJ. Plasmodium knowlesi: The fifth human malaria parasite. Clin Infect Dis. 2008;46:172–3.

    Article  CAS  PubMed  Google Scholar 

  3. Cox-Singh J, Davis TME, Lee K-S, Shamsul SSG, Matusop A, Ratnam S, et al. Plasmodium knowlesi malaria in humans is widely distributed and potentially life threatening. Clin Infect Dis. 2008;46:165–71.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Singh B, Kim Sung L, 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  PubMed  Google Scholar 

  5. Kantele A, Jokiranta TS. Review of cases with the emerging fifth human malaria parasite, Plasmodium knowlesi. Clin Infect Dis. 2011;52:1356–62.

    Article  PubMed  Google Scholar 

  6. Antinori S, Galimberti L, Milazzo L, Corbellino M. Plasmodium knowlesi: the emerging zoonotic malaria parasite. Acta Trop. 2013;125:191–201.

    Article  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  PubMed  PubMed Central  Google Scholar 

  9. William T, Rahman HA, Jelip J, Ibrahim MY, Menon J, Grigg MJ, et al. Increasing incidence of Plasmodium knowlesi malaria following control of P. falciparum and P. vivax malaria in Sabah, Malaysia. PLoS Negl Trop Dis. 2013;7:e2026.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Rajahram GS, Barber BE, William T, Menon J, Anstey NM, Yeo TW. Deaths due to Plasmodium knowlesi malaria in Sabah, Malaysia: association with reporting as Plasmodium malariae and delayed parenteral artesunate. Malar J. 2012;11:284.

    Article  PubMed  PubMed Central  Google Scholar 

  11. 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  PubMed  PubMed Central  Google Scholar 

  12. 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  PubMed  PubMed Central  Google Scholar 

  13. Singh B, Daneshvar C. Human infections and detection of Plasmodium knowlesi. Clin Microbiol Rev. 2013;26:165–84.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Barber BE, William T, Grigg MJ, Yeo TW, Anstey NM. Limitations of microscopy to differentiate Plasmodium species in a region co-endemic for Plasmodium falciparum, Plasmodium vivax and Plasmodium knowlesi. Malar J. 2013;12:8.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Garnham PCC. Malaria parasites and other haemosporidia. Oxford: Blackwell Scientific Publications Ltd.; 1966.

    Google Scholar 

  16. 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  PubMed  Google Scholar 

  17. Sallum MAM, Peyton EL, Wilkerson RC. Six new species of the Anopheles Leucosphyrus group, reinterpretation of An. elegans and vector implications. Med Vet Entomol. 2005;19:158–99.

    Article  CAS  PubMed  Google Scholar 

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

  19. Khoon CC. Status of malaria vectors in Malaysia. Southeast Asian J Trop Med Publ Health. 1985;16:133–8.

    CAS  Google Scholar 

  20. Vythilingam I. Plasmodium knowlesi in humans: a review on the role of its vectors in Malaysia. Trop Biomed. 2010;27:1–12.

    PubMed  Google Scholar 

  21. Collins WE, Contacos PG, Guinn EG. Studies on the transmission of simian malarias II. Transmission of the H strain of Plasmodium knowlesi by Anopheles balabacensis balabacensis. J Parasitol. 1967;53:841–4.

    Article  CAS  PubMed  Google Scholar 

  22. Anderson JR. Sleep-related behavioural adaptations in free-ranging anthropoid primates. Sleep Med Rev. 2000;4:355–73.

    Article  PubMed  Google Scholar 

  23. Pfeifer M, Lefebvre V, Turner E, Cusack J, Khoo M, Chey VK, et al. Deadwood biomass: an underestimated carbon stock in degraded tropical forests? Environ Res Lett. 2015;10:044019.

    Article  Google Scholar 

  24. Marsh CW, Greer AG. Forest land-use in Sabah, Malaysia: an introduction to Danum Valley. Philos Trans R Soc Lond B Biol Sci. 1992;335:331–9.

    Article  Google Scholar 

  25. Ewers RM, Didham RK, Fahrig L, Ferraz G, Hector A, Holt RD, et al. A large-scale forest fragmentation experiment: the stability of altered forest ecosystems project. Philos Trans R Soc Lond B Biol Sci. 2011;366:3292–302.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Struebig MJ, Turner A, Giles E, Lasmana F, Tollington S, Bernard H, et al. Quantifying the biodiversity value of repeatedly logged rainforests: gradient and comparative approaches from Borneo. Adv Ecol Res. 2013;48:183–224.

    Article  Google Scholar 

  27. Ellwood MDF, Foster WA. Line insertion techniques for the study of high forest canopies. Selbyana. 2001;22:97–102.

    Google Scholar 

  28. Canopy Access Ltd. http://www.canopyaccess.co.uk (2016). Accessed 24 Jun 2016.

  29. Reid JA. Anopheline mosquitoes of Malaya and Borneo. Stud Inst Med Res Malaysia. 1968;31:1–520.

    Google Scholar 

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

  31. Rattanarithikul R, Harrison BA, Panthusiri P, Coleman RE. Illustrated keys to the mosquitoes of Thailand. I. Background; geographic distribution; lists of genera, subgenera, and species; and a key to the genera. Southeast Asian J Trop Med Publ Health. 2005;36:1–80.

    Google Scholar 

  32. Rattanarithikul R, Harbach RE, Harrison BA, Panthusiri P, Jones JW, Coleman RE. Illustrated keys to the mosquitoes of Thailand. II. Genera Culex and Lutzia. Southeast Asian J Trop Med Publ Health. 2005;36:1–97.

    Google Scholar 

  33. Rattanarithikul R, Harrison BA, Panthusiri P, Peyton EL, Coleman RE. Illustrated keys to the mosquitoes of Thailand. III. Genera Aedeomyia, Ficalbia, Mimomyia, Hodgesia, Coquillettidia, Mansonia, and Uranotaenia. Southeast Asian J Trop Med Publ Health. 2006;37:1–85.

    Google Scholar 

  34. Rattanarithikul R, Harrison BA, Harbach RE, Panthusiri P, Coleman RE. Illustrated keys to the mosquitoes of Thailand. IV. Anopheles. Southeast Asian J Trop Med Publ Health. 2006;37:1–128.

    Google Scholar 

  35. Rattanarithikul R, Harbach RE, Harrison BA, Panthusiri P, Coleman RE, Richardson JH. Illustrated keys to the mosquitoes of Thailand. VI. Tribe Aedini. Southeast Asian J Trop Med Publ Health. 2010;41:1–225.

    Google Scholar 

  36. R Core Development Team. R: a language and environment for statistical computing. Austria; 2014.

  37. Oksanen J, Blanchet FG, Kindt R, Legendre P, Minchin PR, O’Hara RB, et al. Vegan: community ecology package. R package version 2.2-1. http://www.CRAN.R-project.org/package=vegan. 2015.

  38. Vavrek MJ. Fossil: palaeoecological and palaeogeographical analysis tools. R package version 0.3.7. https://www.cran.r-project.org/web/packages/fossil/fossil.pdf. 2015.

  39. Bates D, Maechler M, Bolker B, Walker S. lme4: linear mixed-effects models using Eigen and S4. R package version 1.1–7. http://www.CRAN.Rproject.org/package=lme4. 2015.

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

    Article  PubMed  PubMed Central  Google Scholar 

  41. Tan CH, Vythilingam I, Matusop A, Chan ST, Singh B. Bionomics of Anopheles latens in Kapit, Sarawak, Malaysian Borneo in relation to the transmission of zoonotic simian malaria parasite Plasmodium knowlesi. Malar J. 2008;7:52.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Hii JLK, Kan S, Vun YS, Chin KF, Tambakau S, Chan MK, 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.

    CAS  PubMed  Google Scholar 

  43. 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  PubMed  PubMed Central  Google Scholar 

  44. Hii JLK. Evidence for the existence of genetic variability in the tendency of Anopheles balabacensis to rest in houses and to bite man. Southeast Asian J Trop Med Publ Health. 1985;16:173–82.

    CAS  Google Scholar 

  45. Hii JLK, Birley MH, Sang VY. Estimation of survival rate and oviposition interval of Anopheles balabacensis mosquitoes from mark-recapture experiments in Sabah, Malaysia. Med Vet Entomol. 1990;4:135–40.

    Article  CAS  PubMed  Google Scholar 

  46. Hii JLK, Chew M, Sang VY, Munstermann LE, Tan SG, Panyim S, et al. Population genetic analysis of host seeking and resting behaviors in the malaria vector, Anopheles balabacensis (Diptera: Culicidae). J Med Entomol. 1991;28:675–84.

    Article  CAS  PubMed  Google Scholar 

  47. Eyles DE, Wharton RH, Cheong WH, Warren M. Studies on malaria and Anopheles balabacensis in Cambodia. Bull World Health Organ. 1964;30:7–21.

    CAS  PubMed  PubMed Central  Google Scholar 

  48. Rundi C. Malaria elimination in Malaysia. Third Annual APMEN technical and business meeting. 9–12 May 2011. Kota Kinabalu; 2011. http://www.apmen.org/apmen-iii-meeting-proceedings/. Accessed 28 Oct 2015.

  49. 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  PubMed  Google Scholar 

  50. Rohani A, Lokman Hakim S, Hassan AR, Chan ST, Ong YF, Abdullah AG, et al. Bionomics of Anopheles balabacensis Baisas, the principal malaria vector in Ranau, Sabah. Trop Biomed. 1999;16:31–8.

    Google Scholar 

  51. 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  PubMed  PubMed Central  Google Scholar 

  52. Manguin S. Anopheles mosquitoes—new insights into malaria vectors. Rijeka: InTech; 2013.

    Book  Google Scholar 

  53. Thang ND, Erhart A, Speybroeck N, Xa NX, Thanh NN, Van Ky P, et al. Long-lasting insecticidal hammocks for controlling forest malaria: a community-based trial in a rural area of central Vietnam. PLoS One. 2009;4:e7369.

    Article  PubMed  PubMed Central  Google Scholar 

  54. Sochantha T, Van Bortel W, Savonnaroth S, Marcotty T, Speybroeck N, Coosemans M. Personal protection by long-lasting insecticidal hammocks against the bites of forest malaria vectors. Trop Med Int Health. 2010;15:336–41.

    Article  CAS  PubMed  Google Scholar 

  55. Bates M. Observations on the distribution of diurnal mosquitoes in a tropical forest. Ecology. 1944;25:159–70.

    Article  Google Scholar 

  56. Clements AN. The biology of mosquitoes. Vol. 2. In: Sensory reception and behaviour. Wallingford: CABI Publishing; 1999.

  57. Bates M. Observations on climate and seasonal distribution of mosquitoes in Eastern Colombia. J Anim Ecol. 1945;14:17–25.

    Article  Google Scholar 

  58. Clements AN. The biology of mosquitoes. Vol. 1. In: Development, nutrition and reproduction. Wallingford: CABI Publishing; 1992.

  59. Silver JB. Mosquito ecology: field sampling methods. 3rd ed. London: Springer; 2008.

    Book  Google Scholar 

  60. Bidlingmayer WL. The effect of moonlight on the flight activity of mosquitoes. Ecology. 1964;45:87–94.

    Article  Google Scholar 

  61. Kampango A, Cuamba N, Charlwood JD. Does moonlight influence the biting behaviour of Anopheles funestus? Med Vet Entomol. 2011;25:240–6.

    Article  CAS  PubMed  Google Scholar 

  62. Charlwood JD, Paru R, Dagaro H, Lagog M. Influence of moonlight and gonotrophic age on biting activity of Anopheles farauti (Diptera, Culicidae) from Papua New Guinea. Bull Entomol Res. 1986;23:132–5.

    CAS  Google Scholar 

  63. Birley MH, Charlwood JD. The effect of moonlight and other factors on the ovipositon cycle of malaria vectors in Madang, Papua New Guinea. An Trop Med Parasitol. 1989;83:415–22.

    CAS  Google Scholar 

  64. Chadee DD. Indoor and outdoor host-seeking rhythms of Anopheles bellator (Diptera: Culicidae) in Trinidad, West Indies. J Med Entomol. 1992;29:567–9.

    Article  CAS  PubMed  Google Scholar 

  65. Miller TA, Stryker RG, Wilkinson RN, Esah S. The influence of moonlight and other environmental factors on the abundance of certain mosquito species in light-trap collections in Thailand. J Med Entomol. 1970;7:555–61.

    Article  CAS  PubMed  Google Scholar 

  66. Davies JB. Moonlight and the biting activity of Culex (Melanoconion) portesi Senevet & Abonnenc and C. (M.) taeniopus D. & K. (Diptera, Culicidae) in Trinidad forests. Bull Entomol Res. 1975;65:81–96.

    Article  Google Scholar 

  67. Rubio-Palis Y. Influence of moonlight on light trap catches of the malaria vector Anopheles nuneztovari in Venezuela. J Am Mosq Control Assoc. 1992;8:178–80.

    CAS  PubMed  Google Scholar 

  68. Souza NA, Andrade-Coelho CA, Silva VC, Peixoto AA, Rangel EF. Moonlight and blood-feeding behaviour of Lutzomyia intermedia and Lutzomyia whitmani (Diptera: Psychodidae: Phlebotominae) vectors of American cutaneous leishmaniasis in Brazil. Mem Inst Oswaldo Cruz. 2005;100:39–42.

    Article  PubMed  Google Scholar 

  69. Singh N, Mishra AK, Curtis CF, Sharma VP. Influence of moonlight on light trap catches of the malaria vector Anopheles culicifacies (Diptera: Culicidae) in India. Bull Entomol Res. 1996;86:475–9.

    Article  Google Scholar 

Download references

Authors’ contributions

HLB conducted the fieldwork, analysis and wrote the manuscript. HLB, RME, IV, CD and JDM were involved in the design of the study. RME and JDM supervised the study. All authors read and approved the final manuscript.

Acknowledgements

For research permission, we thank the Sabah Biodiversity Council, Royal Society Southeast Asia Rainforest Research Project, Maliau Basin Management Committee, SAFE Project, Danum Valley Management Committee, Benta Wawasan, National Institute of Health, Institute for Medical Research, Medical Research and Ethics Committee and Ministry of Health. Thanks to SAFE Project coordinators (M. Khoo, R. Gray) and research assistants. We would like to thank the Natural Environmental Research Council (NE/I527902/1) and the Sime Darby Foundation for funding this work. This paper represents a contribution to Imperial College’s Grand Challenges in Ecosystems and the Environment initiative.

Competing interests

The authors declare that they have no competing interests.

Availability of data and material

Data is available in the SAFE Project database (http://www.safe.earthcape.com).

Ethics approval and consent to participate

This project was approved by the Medical Research and Ethics Committee (NMRR-12-689-12521), Ministry of Health, Malaysia and Imperial College London Research Ethics Committee (ICREC-12-5-6), UK.

Funding

This research was funded by Natural Environmental Research Council (NE/I527902/1).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Hayley L. Brant.

Additional file

12936_2016_1416_MOESM1_ESM.pdf

Additional file 1. Species accumulation curves for human landing catch sampling in primary forest at (a) Ground and (b) Canopy level, virgin jungle reserve at (c) Ground and (d) Canopy level, and logged forest at (e) Ground and (f) Canopy level. Shaded area indicates 95 % confidence intervals.

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

Brant, H.L., Ewers, R.M., Vythilingam, I. et al. Vertical stratification of adult mosquitoes (Diptera: Culicidae) within a tropical rainforest in Sabah, Malaysia. Malar J 15, 370 (2016). https://doi.org/10.1186/s12936-016-1416-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12936-016-1416-1

Keywords