Life-table studies revealed significant effects of deforestation on the development and survivorship of Anopheles minimus larvae
© The Author(s). 2016
Received: 16 February 2016
Accepted: 27 May 2016
Published: 6 June 2016
Many developing countries are experiencing rapid ecological changes such as deforestation and shifting agricultural practices. These environmental changes may have an important consequence on malaria due to their impact on vector survival and reproduction. Despite intensive deforestation and malaria transmission in the China-Myanmar border area, the impact of deforestation on malaria vectors in the border area is unknown.
We conducted life table studies on Anopheles minimus larvae to determine the pupation rate and development time in microcosms under deforested, banana plantation, and forested environments.
The pupation rate of An. minimus was 3.8 % in the forested environment. It was significantly increased to 12.5 % in banana plantations and to 52.5 % in the deforested area. Deforestation reduced larval-to-pupal development time by 1.9–3.3 days. Food supplementation to aquatic habitats in forested environments and banana plantations significantly increased larval survival rate to a similar level as in the deforested environment.
Deforestation enhanced the survival and development of An. minimus larvae, a major malaria vector in the China-Myanmar border area. Experimental determination of the life table parameters on mosquito larvae under a variety of environmental conditions is valuable to model malaria transmission dynamics and impact by climate and environmental changes.
Many developing countries are experiencing rapid ecological changes such as deforestation and shifting agricultural practices [1–4]. These ecological changes may have an important consequence on vector-borne disease transmission due to their impact on vector survival and reproduction [5–8]. In southeast Asia, including China and Myanmar, deforestation and cultivation of cash crops represent the most important environmental changes in rural areas [9, 10]. Many forests were depleted through illegal logging, agricultural clearing, and land development for housing and hydroelectric projects. Deforestation has led to major changes in the environment and subsequently may affect the ecology of malaria vectors. Deforestation may provide more favorable conditions for the larval development of anopheline species. For example, enhanced larval survivorship and faster larval-to-pupal development were found in Anopheles gambiae and An. arabiensis larvae in habitats from open-canopy areas in comparison to shaded habitats [11, 12]. It is possible that deforestation may be detrimental to the survival of other mosquito species. For example, An. dirus larvae breed in well shaded, temporary ground pools, or slow-moving water [13, 14]. These breeding habitats normally occur in the deep forest, making An. dirus a common vector at the forest edges . Consequently, environmental disturbance such as deforestation can have positive or negative impacts on malaria vector populations, depending on the ecology of the vector species.
The malaria elimination goal set by China, Thailand and other countries in the Greater Mekong subregion faces a number of challenges [16, 17]. These challenges include parasite re-introduction by migratory human populations from the endemic neighboring country Myanmar. For example, among malaria patients over 15 years old in Yunnan, imported cases constituted 22 % of total cases in men and 13 % in women, and Myanmar was the predominant source of infection . The second is intense malaria transmission in the remote border areas where impoverished minority ethnic groups are concentrated and resources for malaria control are very limited . In the border area between China and Myanmar, illegal logging of the tropical rainforest and forest clearing for cash crops have been particularly severe [20, 21]. There is no information on how these environmental changes alter the ecology of malaria vectors in the border area. However, this knowledge is critical to malaria control in the border area where malaria burden is the highest and residents are most prone to malaria epidemics.
The objective of this study is to determine the impact of deforestation and agricultural land development on the larval ecology of An. minimus, a key malaria vector species in China-Myanmar and Thai-Myanmar border areas. Knowledge of the response to anopheline vectors to environmental changes will provide a better understanding of malaria transmission dynamics, which is crucial for implementing effective malaria control strategies.
Blood-engorged Anopheles adult mosquitoes were collected from the deforested area and separated to species by morphology. F1 eggs were collected and allowed to hatch in the insectary. The female parents were collected for molecular identification using the multiplex polymerase chain reaction (PCR) method [22, 23]. Because An. minimus was the predominant vector species in the field collection, it was the focus of the present study. We constructed microcosms using plastic bins of 24 cm in diameter and 27 cm deep (Fig. 1). Two liters of rainwater were added to the plastic bin and allowed to acclimate for 1 day. To prevent other insects from invading the microcosms or other mosquitoes from laying eggs, the microcosms were placed in an insect-proof 61 × 61 × 61 cm3 BugDorm tent (BioQuip, Rancho Dominguez, CA, USA) (Fig. 1). All sides of the BugDorm tent were made of clear polyester netting materials, so that sunlight was not blocked. Microcosms mimicked natural habitats, and were homogeneous. The homogeneous feature of microcosms had the advantage over the natural habitats which were highly variable in habitat size, larval food conditions (e.g. organic matters), vegetation coverage and predators.
Fifty newly hatched An. minimus larvae were introduced to each microcosm. Each day the number of surviving larvae was counted and their stage of development recorded. Pupae were counted and removed daily. Water levels in the microcosm were checked daily, and were maintained by adding rainwater if needed. Water temperature was measured using the HOBO data loggers placed in the microcosm, 1 cm below the water surface. There were 8 microcosms in each of the three land use and land cover types (deforested, banana field and forest area). Microcosms in the forested environment may not receive sufficient sunlight for microbial growth and thus mosquito larvae failed to develop due to a lack of necessary nutrients. To test this hypothesis, 12 microcosms (4 replicates per land use type) were supplemented with Tetramin fish food daily, and daily larval mosquito survivorship was measured.
The pupation rate of An. minimus larvae was calculated as the proportion of 1st instar larvae that developed into pupae. Mean larval-to-pupal development time was calculated. The t-test and analysis of variance (ANOVA) were used to determine the effect of land use and land cover types and larval food supplementation on pupation rate and larval development time wherever appropriate. Kaplan-Meier survival analysis was performed to determine the effects on larval survivorship. The log-rank test was used to test the statistical significance. Stage-specific larval development time and mortality rate were calculated. All data analyses were performed using STATISTICA 10.0 (StatSoft, Tulsa, USA).
Field-collected mosquito species
A total of 218 blood-engorged Anopheles adult mosquitoes were collected from indoors and outdoors (pig and cow shelter) using aspirators. Among these, 94 % were An. minimus and the remaining mosquitoes were An. maculatus. The PCR method confirmed morphological identification.
Effects of land use and land cover on An. minimus larvae
Developmental stage-specific survivorship of immature An. minimus under three land use and land cover environments under natural condition and larval food supplementation
Development time (day)
Cumulative survival rate
Stage mortality rate
Development time (day)
Cumulative survival rate
Stage mortality rate
Effects of food supplementation
With food supplementation, pupation rate increased significantly for all three land use and land cover types (Fig. 2a). A particularly large change was found in banana fields and forested environments where a 3.7 fold (t (1) = 2.92, P = 0.02) and 15.4 fold (t (1) = 9.26, P < 0.0001) increase was found in comparison to the natural conditions. The pupation rate did not vary significantly among the three land use and land cover types when larval food was added to the microcosms, suggesting that sunlight contributed significantly to the process of larval food production in the habitats. Stage-specific mortality analysis found a dramatic reduction in the mortality rate for all 4 larval instars (Table 1). The average larval-to-pupal development pupation time did not vary significantly among the three land cover types. Comparing the treatment without food supplementation, food supplementation reduced larval development time by 3.1 days (t (1) = 1.66, P = 0.17) in the deforested area, by 3.3 days in the banana field (t (1) = 2.74, P = 0.05) and by 1.9 days (t (1) = 1.25, P = 0.34) in the forested environment (Fig. 2b). Kaplan-Meier survivor analysis found similar mortality rate across time for the three land use and land cover types (Fig. 3).
Effects of land use and land cover on habitat temperature
The present study found that in comparison to larval habitats in the forested area, the pupation rate of An. minimus was about 4-fold higher in banana plantation fields and 13-fold higher in deforested area, and the larval-to-pupal development time was shortened by 1.9–3.3 days. The microcosms were covered by insect-proof materials and no predator had been observed in the microcosms. Because the experiments used microcosms made of the same rain water and same substrate, and the only difference was the environment in which the microcosms were placed (forested, banana field and deforested); the difference in pupation rate and development time of mosquito larvae among the three environments was entirely due to land use and land cover differences. This finding was consistent with the study in the Peruvian Amazon where An. darlingi larvae were more commonly found in open sunlit aquatic habitats . Tuno et al. demonstrated that the survivorship of An. gambiae larvae was reduced from 55 to 57 % in habitats fully exposed to sunlight to 1–2 % in habitats with full forest canopy coverage (forest habitats) in western Kenya highlands . Whether the observation that deforestation facilitates the survival of the immature stage of anopheline mosquitoes is a common phenomenon across various malaria vector species and across multiple continents warrants further investigation.
The mechanisms for reduced larval survivorship in microcosms placed in the forested environment and banana plantation with high canopy coverage are unknown. Land use and land cover may modify the temperature and food source availability for mosquito larvae. We observed an average of 0.4 °C lower temperature in microcosms in the forested environment than deforested environment although the minimum and maximum temperature in forested and sunlit deforested microcosms ranged between 25.5–30.2 °C, within the optimal temperature for mosquito development. Lower pupation rate in habitats from the forest environment could lead to large early-instar larval mortality, subsequently leading to lower pupation rate as observed in Culex quinquefasciatus . The larval density used in the present study was comparable to the density in natural habitats. Field sampling in over 100 habitats of various habitat types found that the mean density of 3rd and 4th instar Anopheles larvae ranged from 1.5–2 per dipper (300 ml volume), or 10–13 larvae per 2 l water (unpublished data). In the present experiments we observed ~20–40 % of the larvae developed into 3rd or 4th instar (or 8–20 larvae in 2 l microcosms). Further, the same larval density used in all experiments enabled us to compare the environmental impact on larval development and survivorship. Another hypothesis is the lack of food for mosquito larvae in the shaded microcosms. This hypothesis was supported by food supplementation to shaded microcosms and resulting recuperation of mosquito pupation rate to a level comparable to the deforested environment. Additionally, faster larval-to-pupal development was found in all three land use and land cover types. The potential food source of anopheline larvae may include bacteria, fungi, diatom, detritus, organic matter, and others. Using stable isotopes of carbon and nitrogen, Gilbreath et al. demonstrated inter-specific resource partitioning between Culex quinquefasciatus and An. gambiae larvae in natural habitats in western Kenya , but resource utilization of bacteria, algae, organic matter or others by Anopheles larvae and their nutritional contribution to mosquito larval biomass are virtually unknown. The abundance and structure of microbes such as algae and photosynthetic cyanobacteria in aquatic habitats may have changed in response to shading . Information on the contribution of specific microbes to larval mosquito development and survival will significantly enhance our understanding of larval ecology.
Our finding on the impact of environmental changes on larval development of malaria vectors has important implications on understanding malaria transmission dynamics in the China-Myanmar border area. First, An. minimus larvae can successfully develop into pupae and adults and complete its life cycle in the forest despite of much reduced pupation rate. If aquatic habitats are abundant in the forest, a considerable number of adult mosquitoes may be produced. Because An. minimus is highly anthropophilic and highly susceptible to malaria parasites [28, 29], this presents a significant risk to forest workers and hunters. Secondly, significantly enhanced larval survivorship and pupation rate in An. minimus mosquitoes in the deforested and cultivated areas render the local residents highly prone to malaria outbreaks because more human-mosquito contact is expected in such environments. Therefore, close malaria surveillance for local residents living near deforested and cultivated areas is needed to prevent outbreaks. Thirdly, climate warming is expected to affect the ecology of Anopheles mosquitoes and their malaria transmission potential. It is imperative to use various tools, including mathematical modeling to examine the impact of climate and environmental changes on vector and malaria transmission dynamics . Larval survivorship and development time are important parameters for modeling transmission dynamics. The present study provided valuable information on these important parameters.
This study has several limitations. First, we examined the response of one of the most abundant malaria vector species to environmental changes. There are multiple malaria vectors in the area, including An. maculatus, An. culicifacies, An. dirus and others . It is possible that other vector species may respond to environmental changes differently. It would be particularly interesting to study the mechanisms and larval ecology of An. dirus that was found exhibiting reduced abundance in deforested areas. Second, we did not examine the effect of deforestation on the survival and reproduction of adult mosquitoes. Enhanced survivorship and blood feeding frequency in adult mosquitoes by deforestation through the effects on the microclimatic condition of the mosquitoes were found in other Anopheles species , and it is possible that such effects are also manifested in An. minimus adults. Third, the study was conducted in microcosms. The findings from the present study need validating in natural habitats.
The present study shows that deforestation significantly enhanced the pupation rate and shortened the larval-to-pupal development time of An. minimus larvae, a major malaria vector in southeast Asia. The results of the present study on the life table parameters of mosquitoes under a variety of environmental conditions is valuable to model malaria transmission dynamics and the impact by climate and environmental changes in southeast Asia.
ANOVA, analysis of variance; PCR, polymerase chain reaction.
The authors wish to thank all field staff for their assistance with the field work and Amruta Dixit for her editorial assistance and Ming-chieh Lee for his assistance in study site map preparation. We thank two anonymous reviewers for their constructive comments.
This study was supported by grants from the National Institutes of Health (U19 AI089672 and D43 TW009527). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Availability of data and materials
The data sets supporting the results are included within the article and its Additional file 1.
XMW, GZ and GY conceived and designed the research; GZ, DZ, XLW, YW and ZY participated in the field work, collected the data and samples. XMW and GZ did the data analysis. XMW and GZ drafted the manuscript, GY edited the manuscript. All authors were involved in the interpretation and discussion of the results and provided comments. All authors read and approved the final version of the manuscript.
The authors declare that they have no competing interests.
Ethics approval and consent to participate
This study did not involve vertebrate animals or human subjects. All necessary permits required for this field study were obtained from Kunming Medical University, China.
Open AccessThis 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.
- Abdullah SA, Hezri AA. From forest landscape to agricultural landscape in the developing tropical country of Malaysia: pattern, process, and their significance on policy. Environ Manag. 2008;42(5):907–17.View ArticleGoogle Scholar
- Chua KB, Chua BH, Wang CW. Anthropogenic deforestation, El Nino and the emergence of Nipah virus in Malaysia. Malays J Pathol. 2002;24(1):15–21.PubMedGoogle Scholar
- Garedew E, Sandewall M, Soderberg U, Campbell BM. Land-use and land-cover dynamics in the central rift valley of Ethiopia. Environ Manag. 2009;44(4):683–94.View ArticleGoogle Scholar
- Rudel TK. The national determinants of deforestation in sub-Saharan Africa. Philos Trans R Soc Lond B Biol Sci. 2013;368(1625):20120405.View ArticlePubMedPubMed CentralGoogle Scholar
- Afrane YA, Little TJ, Lawson BW, Githeko AK, Yan G. Deforestation and vectorial capacity of Anopheles gambiae Giles mosquitoes in malaria transmission, Kenya. Emerg Infect Dis. 2008;14(10):1533–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Vittor AY, Pan W, Gilman RH, Tielsch J, Glass G, Shields T, et al. Linking deforestation to malaria in the Amazon: characterization of the breeding habitat of the principal malaria vector, Anopheles darlingi. Am J Trop Med Hyg. 2009;81(1):5–12.PubMedPubMed CentralGoogle Scholar
- Walsh JF, Molyneux DH, Birley MH. Deforestation: effects on vector-borne disease. Parasitology. 1993;106(Suppl):S55–75.View ArticlePubMedGoogle Scholar
- Chaves LF, Koenraadt CJ. Climate change and highland malaria: fresh air for a hot debate. Q Rev Biol. 2010;85(1):27–55.View ArticlePubMedGoogle Scholar
- Hillman AL, Yu JQ, Abbott MB, Cooke CA, Bain DJ, Steinman BA. Rapid environmental change during dynastic transitions in Yunnan Province, China. Quaternary Sci Rev. 2014;98:24–32.View ArticleGoogle Scholar
- Webb EL, Jachowski NRA, Phelps J, Friess DA, Than MM, Ziegler AD. Deforestation in the Ayeyarwady Delta and the conservation implications of an internationally-engaged Myanmar. Global Environ Chang. 2014;24:321–33.View ArticleGoogle Scholar
- Munga S, Minakawa N, Zhou G, Mushinzimana E, Barrack OO, Githeko AK, et al. Association between land cover and habitat productivity of malaria vectors in western Kenyan highlands. Am J Trop Med Hyg. 2006;74(1):69–75.PubMedGoogle Scholar
- Afrane YA, Zhou G, Lawson BW, Githeko AK, Yan G. Life-table analysis of Anopheles arabiensis in western Kenya highlands: effects of land covers on larval and adult survivorship. Am J Trop Med Hyg. 2007;77(4):660–6.PubMedGoogle Scholar
- Rattanarithikul R, Green CA, Panyim S, Noigamol C, Chanaimongkol S, Mahapibul P. Larval habitats of malaria vectors and other Anopheles mosquitoes around a transmission focus in northwestern Thailand. J Am Mosq Control Assoc. 1995;11(4):428–33.PubMedGoogle Scholar
- Rosenberg R, Maheswary NP. Forest malaria in Bangladesh. II. Transmission by Anopheles dirus. Am J Trop Hyg. 1982;31(2):183–91.Google Scholar
- Obsomer V, Defourny P, Coosemans M. The Anopheles dirus complex: spatial distribution and environmental drivers. Malar J. 2007;6:26.View ArticlePubMedPubMed CentralGoogle Scholar
- Hsiang MS, Gosling RD. Striding toward malaria elimination in China. Am J Trop Med Hyg. 2015;93(2):203–4.View ArticlePubMedPubMed CentralGoogle Scholar
- Maude RJ, Lubell Y, Socheat D, Yeung S, Saralamba S, Pongtavornpinyo W, et al. The role of mathematical modelling in guiding the science and economics of malaria elimination. Int Health. 2010;2(4):239–46.View ArticlePubMedPubMed CentralGoogle Scholar
- Lin H, Lu L, Tian L, Zhou S, Wu H, Bi Y, et al. Spatial and temporal distribution of falciparum malaria in China. Malar J. 2009;8:130.View ArticlePubMedPubMed CentralGoogle Scholar
- Cui L, Yan G, Sattabongkot J, Cao Y, Chen B, Chen X, et al. Malaria in the Greater Mekong Subregion: heterogeneity and complexity. Acta Trop. 2012;121(3):227–39.View ArticlePubMedGoogle Scholar
- Achard F, Eva HD, Stibig HJ, Mayaux P, Gallego J, Richards T, et al. Determination of deforestation rates of the world's humid tropical forests. Science. 2002;297(5583):999–1002.View ArticlePubMedGoogle Scholar
- Gibbs HK, Ruesch AS, Achard F, Clayton MK, Holmgren P, Ramankutty N, et al. Tropical forests were the primary sources of new agricultural land in the 1980s and 1990s. Proc Natl Acad Sci U S A. 2010;107(38):16732–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang Y, Zhong D, Cui L, Lee MC, Yang Z, Yan G, et al. Population dynamics and community structure of Anopheles mosquitoes along the China-Myanmar border. Parasit Vectors. 2015;8:445.View ArticlePubMedPubMed CentralGoogle Scholar
- Yu G, Yan G, Zhang N, Zhong D, Wang Y, He Z, et al. The Anopheles community and the role of Anopheles minimus on malaria transmission on the China-Myanmar border. Parasit Vectors. 2013;6(1):264.View ArticlePubMedPubMed CentralGoogle Scholar
- Tuno N, Okeka W, Minakawa N, Takagi M, Yan G. Survivorship of Anopheles gambiae sensu stricto (Diptera: Culicidae) larvae in western Kenya highland forest. J Med Entomol. 2005;42(3):270–7.View ArticlePubMedGoogle Scholar
- Chaves LF, Keogh CL, Vazquez-Prokopec GM, Kitron UD. Combined sewage overflow enhances oviposition of Culex quinquefasciatus (Diptera: Culicidae) in urban areas. J Med Entomol. 2009;46(2):220–6.View ArticlePubMedGoogle Scholar
- Gilbreath 3rd TM, Kweka EJ, Afrane YA, Githeko AK, Yan G. Evaluating larval mosquito resource partitioning in western Kenya using stable isotopes of carbon and nitrogen. Parasit Vectors. 2013;6:353.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang Y, Gilbreath 3rd TM, Kukutla P, Yan G, Xu J. Dynamic gut microbiome across life history of the malaria mosquito Anopheles gambiae in Kenya. PLoS One. 2011;6(9):e24767.View ArticlePubMedPubMed CentralGoogle Scholar
- Polsomboon S, Poolprasert P, Suwonkerd W, Bangs MJ, Tanasinchayakul S, Akratanakul P, et al. Biting patterns of Anopheles minimus complex (Diptera: Culicidae) in experimental huts treated with DDT and deltamethrin. J Vector Ecol. 2008;33(2):285–92.View ArticlePubMedGoogle Scholar
- Tisgratog R, Tananchai C, Juntarajumnong W, Tuntakom S, Bangs MJ, Corbel V, et al. Host feeding patterns and preference of Anopheles minimus (Diptera: Culicidae) in a malaria endemic area of western Thailand: baseline site description. Parasit Vectors. 2012;5:114.View ArticlePubMedPubMed CentralGoogle Scholar
- Reiner Jr RC, Perkins TA, Barker CM, Niu T, Chaves LF, Ellis AM, et al. A systematic review of mathematical models of mosquito-borne pathogen transmission: 1970-2010. J R Soc Interface. 2013;10(81):20120921.View ArticlePubMedPubMed CentralGoogle Scholar
- Afrane YA, Lawson BW, Githeko AK, Yan G. Effects of microclimatic changes caused by land use and land cover on duration of gonotrophic cycles of Anopheles gambiae (Diptera: Culicidae) in western Kenya highlands. J Med Entomol. 2005;42(6):974–80.View ArticlePubMedGoogle Scholar