Spatial characterization of colonies of the flying fox bat, a carrier of Nipah Virus in Thailand
© Thanapongtharm et al.; licensee Biomed Central. 2015
Received: 8 September 2014
Accepted: 10 March 2015
Published: 28 March 2015
A major reservoir of Nipah virus is believed to be the flying fox genus Pteropus, a fruit bat distributed across many of the world’s tropical and sub-tropical areas. The emergence of the virus and its zoonotic transmission to livestock and humans have been linked to losses in the bat’s habitat. Nipah has been identified in a number of indigenous flying fox populations in Thailand. While no evidence of infection in domestic pigs or people has been found to date, pig farming is an active agricultural sector in Thailand and therefore could be a potential pathway for zoonotic disease transmission from the bat reservoirs. The disease, then, represents a potential zoonotic risk. To characterize the spatial habitat of flying fox populations along Thailand’s Central Plain, and to map potential contact zones between flying fox habitats, pig farms and human settlements, we conducted field observation, remote sensing, and ecological niche modeling to characterize flying fox colonies and their ecological neighborhoods. A Potential Surface Analysis was applied to map contact zones among local epizootic actors.
Flying fox colonies are found mainly on Thailand’s Central Plain, particularly in locations surrounded by bodies of water, vegetation, and safe havens such as Buddhist temples. High-risk areas for Nipah zoonosis in pigs include the agricultural ring around the Bangkok metropolitan region where the density of pig farms is high.
Passive and active surveillance programs should be prioritized around Bangkok, particularly on farms with low biosecurity, close to water, and/or on which orchards are concomitantly grown. Integration of human and animal health surveillance should be pursued in these same areas. Such proactive planning would help conserve flying fox colonies and should help prevent zoonotic transmission of Nipah and other pathogens.
KeywordsFlying foxes Nipah Species distribution model Ensemble modeling Potential surface analysis
Habitat loss is the greatest threat to wildlife and biodiversity. The loss and fragmentation of wildlife habitats can lead to increasing contact among wildlife, domestic animals, and people, potentially leading to the emergence and spread of zoonotic diseases . The Nipah virus (NiV) is one such pathogen. The novel RNA paramyxovirus (genus Henipavirus), closely related to Hendra virus, is named after the village Sungai Nipah in the State of Negeri Sembilan, Malaysia from which the virus was first isolated from a human patient in 1998 . In humans, NiV causes Nipah virus infection, presenting a range of clinical outcomes, from asymptomatic infection to acute respiratory syndrome and fatal encephalitis .
Investigations of the origins of NiV identified the flying fox genus Pteropus to be a major reservoir [4,5]. Subclinical infections have been found in flying fox populations in Malaysia, Cambodia, Thailand and Madagascar [4-8]. Flying foxes are mammals, members of the Pteropididae or fruit bat family, and are the largest of all bats . They are found throughout tropical and sub-tropical Asia and Australia and on islands of the Indian Ocean and the western Pacific . Pteropididae play a crucial role in rainforest ecosystems . They pollinate flowers and disperse seeds as they forage on the nectar and pollen of plants and on the fruits of rainforest trees and vines . In Thailand, flying foxes are protected by the Wildlife Preservation and Protection Act, B.E. 2535 (1992), which forbids hunting protected wild animals and protects wildlife sanctuaries. A better understanding of the flying fox and its habitat preferences and dispersal would be a useful contribution to its conservation in Thailand. In addition, such an investigation should help efforts in better preventing potential disease transmission.
Work outside Thailand shows that in response to losses in its natural foraging areas, the adaptive Pteropus have turned to foraging in orchards, including those grown on pig farms where the NiV it carries are intermittently passed to pigs via urine or the contamination of partially-eaten fruit [4,5,11]. Investigation showed the virus to subsequently spill over from pigs to other animals and humans via respiratory droplets or close contact [2,12]. Pig farmers and workers exposed to respiratory illness and encephalitis in pigs were the first group of humans infected with the virus . In 1999, abattoir workers in Singapore developed Nipah virus encephalitis . Investigation showed direct contact with live pigs imported from Malaysia appeared to be the most important risk factor for those infections . In contrast, a retrospective study of human cases in Bangladesh in 1999, the consumption of raw date palm sap proved one of the main risk factors of infection [16-18]. The result suggests NiV may have passed directly from bats to humans without an amplification host, as was apparently the case in Malaysia [11,12]. Human-to-human transmission was observed in several outbreaks in Bangladesh and India [18-20].
The situation of Nipah virus infection in Thailand showed that there has been no evidence of the viruses in domestic animals but they have been found in wildlife. Thailand’s National Institute of Animal Health (NIAH), the Department of Livestock Development (DLD)’s central laboratory, conducted a retrospective study of all specimens of swine interstitial pneumonia submitted during 1998 to 2001 using immunohisto-chemistry (IHC) technique . All samples reported negative for NiV. Since 2002, The DLD has conducted a sero-surveillance of 4,000 – 5,000 samples of pig per year by using Modified ELISA technique. The pig blood samples have been collected in high pig density areas and bordering area of Thailand and Malaysia (south). Simultaneously, the veterinarians of the DLD have conducted clinical surveillance by investigating any suspected cases of NiV, they can consider to collect samples submitting for laboratory confirmation  but NiV has never been found so far . On the other hand, the Molecular Biology Laboratory for Neurological Diseases, Chulalongkorn University Hospital conducted surveillance for NiV antibody by using enzyme immunoassay and for NiV by using the duplex reverse transcription–polymerase chain reaction (RT-PCR) in Thailand’s bat population during 2002–2004. The results showed 82 of 1,304 positives to NiV antibody and the tests for NiV presence in the urine and saliva of 12 bat species produced positives for 3 species of fruit bats (P. hypomelanus, P. vampyrus, and P. lylei) and 1 species of insect-eating bat (Hipposideros larvatus) with being a probable accidental case . In only one species of flying fox (P. lylei) was NiV found in both saliva and urine. A longitudinal study subsequently conducted on P. lylei populations between 2005 to 2007 in Thailand showed that 2 NiV strains previously identified circulating in Malaysia and Bangladesh were found in the bat’s urine . The study also highlighted a seasonal pattern with peaks between April and June, when viral RNA could be detected in urine. This seasonal pattern was associated with the observed fluctuation of population numbers, as May corresponds to the time of the year when young bats fledge .
The objectives of the present study were threefold. First, we aimed to describe the characteristics of the flying fox colonies and their neighborhoods in the central plain of Thailand (including central and eastern Thailand) from field observations, remote sensing (RS), and geographic information systems (GIS) data. Second, we aimed to predict the potential distribution of flying foxes in the study area using species distribution models (SDM). Finally, we aimed to map the areas where the three key elements of NiV ecology coincide, specifically flying fox habitat, human population, and pig farms, with the aim of informing NiV surveillance on the central plain of Thailand.
Characteristics of bat colonies and their vicinities
The GIS layer of the colonies was overlaid on other layers, including of bodies of water, human population density, elevation, and land cover. The vector map of permanent bodies of water was provided by the Ministry of Transportation. A human population density raster map at 100 m resolution was obtained from the Worldpop project . We used the SRTM elevation database with 90 m spatial resolution produced by NASA .
A land cover map was developed using LANDSAT images with the Exelis VIS ENVI image processing software. Eleven scenes of the LANDSAT 7 Enhanced Thematic Mapper Plus (LANDSAT 7-ETM+) were used to cover the 23 provinces of the study area (path/row = 131/49-50, 130/49-51,129/49-51,128/49-51). The LANDSAT 7 ETM+ sensor has six optical spectral bands at 30 m spatial resolution and one panchromatic band at 15 m spatial resolution and a 16 day revisit cycle. We searched the LANDSAT image archive at the United States Geological Survey EROS Data Center (http://glovis.usgs.gov) and downloaded images with low cloud cover acquired in January 2014. All images were mosaicked and the minimum distance technique supervised classification method  was used to classify images into 4 land cover types most-related to flying fox habitat, including forest, irrigated vegetation, settlement/rainfed vegetation, and bodies of water. The regions of interest (ROIs) were built as classification training sets using ground truth data, 2D scatter plots, visible composition images, and spectral profiles. We evaluated the accuracy of the classes with 100 points per class of additional ground truth data and high-resolution data from Google Earth images for accuracy. Overall accuracy was 93%, with 98% accuracy for forest, 95% for irrigated vegetation, 91% for settlement/rainfed vegetation, and 86% for bodies of water, results considered acceptable and sufficient for the analysis.
For each bat colony polygon, we estimated summary descriptive statistics on the environment, geography, and anthropogenic variables. Specifically, we estimated the area of each colony, the distance from each colony to its nearest neighbor (another colony), the distance from each colony to the nearest body of water, the distance from each colony to the nearest temple, the average elevation within the colony, the average human population density, the proportion of irrigated vegetation land cover in a 10 km buffer around the colony, and the mean normalized difference vegetation index (NDVI) within the colony acquired from LANDSAT images. The vector map of Buddhist temples was provided by the Ministry of Transportation.
Species distribution models
In this study, we used ensemble modeling (EM) (or consensus methods or ensemble forecasting) with the ‘dismo’ and ‘raster’ packages in R, which combines the predictions from several different statistical modeling techniques into a single prediction. Species distribution models (SDM) were initially used to map the ecological suitability for flying fox colonies across the study area. SDM can be used to predict the geographical distribution of species as a function of a series of spatial variables, as they relate species distribution data (occurrence or abundance in known locations) to information on the environmental and/or spatial characteristics of those locations . They have been widely used both for describing patterns and making predictions across terrestrial, freshwater and marine ecosystems [30,31]. Flying fox colonies can occasionally move, and such modeling should allow inferring other areas to which colonies might move, even those from which they are presently absent. The variables used to build the models were selected according to field observations and the results of the descriptive analysis. This would allow, for example, to map areas where colonies are not present at the time of the study, but where the colonies may move in the future if they are too disturbed, or if their current habitat became degraded. The seven different SDM methods used in analyses include: Bioclim, Domain, generalized linear model (GLM), generalized additive model (GAM), maximum entropy model (Maxent), boosted regression tree (BRT), and random forests (RF). The Bioclim and Domain are presence-only modeling methods. Bioclim characterizes the occurrences that are located within the environmental hyper-space occupied by a species, whereas Domain is a distance-based method that assesses new locations in terms of their environmental similarity to locations of presence . The GLM and GAM are presence-absence models based on the regression framework. The GLM is a generalization of ordinary least squares regression using maximum likelihood allowing the linear model to be related to the response variable via a logit link function. The GAM is an extension of the GLM, where the linear predictor is the sum of smoothing functions. It is more flexible and much as machine learning methods can fit very complex functions .
Maxent, BRT, and RF are machine learning methods using presence-absence data. Maxent, sometimes misleadingly referred to as presence-only methods, actually does require the use of background data . It estimates species’ distributions by finding the distribution of maximum entropy (i.e. closest to uniform) subject to the constraint that the expected value of each environmental variable (or its transform and/or interactions) under the distribution matches its empirical average . BRT combines the strengths of two algorithms, regression trees and boosting, creating a single best model from a large numbers of relatively simple models, each formed by a regression tree . RF combines tree predictors such that each tree depends on the values of a random vector sampled independently and with the same distribution for all trees in the forest . When compared with other methods, RF shows a very high accuracy, an ability to model complex interactions among predictors, and the flexibility to perform several types of statistical analysis . The predictions of the seven SDM were then combined into a single ensemble model prediction by weighting each prediction by the performance of its source model, a procedure called ensemble modeling (EM) recognized as producing significantly more robust predictions than all the single models alone [38-42].
For three reasons, all our models were subject to 10 bootstraps. First, there was a very low proportion of positive samples in our data set, which can introduce bias into the logistic regression analysis framework. So for each trial we bootstrapped a different set of pseudo-absences [33,43]. Second, the bootstrapping also aims at preventing over-fitting. That is, we aim at avoiding modeling the noise rather than the main pattern in the data by assembling across a population of models trained with different subsets of data. Third, the pseudo-absences were distributed within a given distance of the presence sites. We wanted to bootstrap through different distance values. The 10 sets of absence data were randomly selected from the background and from 6–15 kilometers beyond the presence sites. Then each set was randomly selected again and divided into two parts equally: a model set used to train the model and a test set used to evaluate the models. These were then used as weights in combining the methods. Nine times the number of positives was randomly selected at each bootstrap to maintain 10% of the positive values of the outcome variable. This 10% ratio was chosen because previous studies compared the various prevalences across models and reported that GAM was not influenced by prevalence, whereas the accuracy increased up to an asymptote when the number of presences reached one tenth of the number of absences for GLM, BRT, and RF . All predictor variables were simultaneously tested in the models.
The performance of the models was evaluated using the area under the curve (AUC) of the receiver operating characteristics (ROC) plots. AUC is a quantitative measure of the overall fit of models that varies from 0.5 (chance event) to 1.0 (perfect fit) . Although AUC was recently criticized as an absolute measure of goodness of fit by many authors, it remains valuable in comparing the performances of several models tested on the same data set .
Mapping the risk area of NiV
Where S represents a summation of scores, W 1 − W n , the weight of each factor according to its importance, and R 1 − R n the rating score of each variable, which corresponds to its scaling into bins. For these maps we assumed two potential scenarios for human infection: 1) humans are directly infected by the virus from the bats, and 2) humans are infected through a pig intermediate host. The overlay corresponding to the first scenario was hence based on three factors: the flying fox distribution map, the distance to the flying fox colonies, and the house density in the sub-district. For the second scenario, the pig farm density map was added to the first three factors. For mapping the overlay of factors important to pig infection, we used three factors: the flying fox distribution map, distance to the flying fox colonies, and the pig farm density at the sub-district level.
Scores given for a Potential Surface Analysis (PSA)
Flying fox distribution zones (probability during 0 to1)
Distance to the flying foxes colonies (km)
Pig farm density (farm per km 2 )
House density (house per km 2 )
0.83x10−3 - 0.83x10−2
0.13 – 907
0.84x10−2 - 0.14
0.15 - 0.27
0.28 - 0.87
This study was approved by the Research Committee of the Bureau of Disease Control and Veterinary Services, Department of Livestock Development (Permit Number: 0601/1325).
Characteristics of the trees roosted by flying foxes
8 and 12
Bamboo is generally found interspersed in many other types of forest and as a pioneer species. It is a fast growing species that easily colonizes disturbed forest sites, both natural and man-made. As such, and due to the logging excesses in Thailand in the past, many bamboo forests have become established in man-made disturbed sites .
17,18, 19, and 21
Trees in mangrove forest are evergreen species with a very dense forest floor. The roots of the trees are both for anchoring it in the soil and for breathing. This type of forest is found close to the seat of the mouth of major rivers where the sea washes ashore. The important tree species include the Kongklang (Rhizophora spp.), Prasak. The plants grow on the forest floor include the various types of sea grasses [73,74].
2,5, and 11
The rubber (Dipterocarpus alatus) is a tropical forest tree, of dense evergreen or mixed dense forests, common in Thailand, Cambodia, Laos and Vietnam. It is a medium-sized to fairly large tree of up to 40 m tall (sometimes more), bole tall, straight, cylindrical, branchless up to 20 m, up to 150 cm in diameter. Leaves narrowly ovate to ovate to elliptical-oblong, 9–25 cm x 3.5-15 cm, base cuneate to rounded, apex acute or shortly indistinctly acuminate .
1,3,4,6,7,9,10,13,14, 15,16, 20 and 22
The various type of trees (mostly are in the Buddhist temple) composed of rubber (Dipterocarpus spp), ficus tree (Ficus spp), bohhi tree (Ficus religiosa L.), bengal almond (Terminalia catappa L.), rain tree (Samanea saman), neem plant (Azadirachta indica), tamarind (Tamarindus indica Linn), sal tree (Shorea robusta Roxb.), bamboo (Bambusa sp.), coconut (Cocos nucifera Linn.), and others . Most of trees are a medium-sized to fairly large tree of up to 40 m tall, which mostly found in the tropical evergreen forest which is distributed in all areas of Thailand. They are concentrated in pockets of high moisture such as valleys and close to water sources such as rivers, streams and mountains. A common characteristic of tropical evergreen forest is the appearance of a lush green vegetation all year round.
Descriptive statistics of the flying foxes’ colonies and vicinity
Size (m 2 )
Distance from each colony
Human density (people/km 2 )
Land cover in 10 km. of radius (m 2 )
To nearest neighbor (km)
To nearest water (m)
To nearest temple (m)
Bodies of water
Our field observations indicated that flying foxes choose a variety of tree types, especially members of evergreen forests, for roosting, even if the trees no longer protect the bats from sunlight. The observations are supported by remote sensing, showing a normalized difference vegetation index acquired from January 2014 LANDSAT imagery with relatively low values at some of the roosting sites. This suggests flying foxes prefer evergreen forests to protect themselves from the sunlight while roosting, but also show a tolerance to trees damaged by bat urine and roosting . We found roosting sites in relatively safe places, including Buddhist temples, islands, the Wildlife Conservation Park, and private lands, as marked by SDMs that included distance to a temple as an important predictor. Even as flying foxes are protected by the Wildlife Preservation and Protection Act, B.E. 2535 (1992), they are still threatened by human hunting, efforts to protect fruit orchards, and informal efforts at disease prevention. Most of the population is Buddhist (>90%) and would largely refrain from threatening animals in the vicinity of temples. Human density appears to correspond positively with roosting sites for temple communities but is negatively associated for the greatest densities in and around urban areas. Some bat colonies are located in private lands and study informants indicated landlords and/or the people in the local community around these sites had tried to protect the bats against hunters. Finally, the other colonies were located in isolated areas such as islands, riverside, and at seaside that are hard to reach by hunters. The Wildlife Conservation Park is closed off as a unit of wildlife conservation. Therefore, all bat colonies, across a variety of locales, were protected from hunters for an array of reasons, including cultural practices, ownership, local sentiment, and remoteness.
The distribution of flying fox colonies is dynamic and changes are observed over time. In 2004–2014, new colonies were observed and a few colonies moved away from their previously known sites . Apart from disturbances caused by hunters, other factors may trigger colony migration. Disturbance by visitors or tourists is assumed to have caused colony no. 21 to move away from its old site to a new isolated site seaside. The roosting trees may have been damaged or killed by the flying foxes themselves by way of their urine and/or roosting . Colony no. 17 moved away from a site with numerous destroyed mangrove forest trees to an adjacent area. Competition with other species using the same habitat could also play an important role as former bat colonies sites were observed colonized by large bird populations. Finally, colonies may move if their sizes increase beyond the capacity of a roosting site, if the foraging areas are reduced or too impacted by urban development, or in relation to the mating season . Colony mobility supports the concept of mapping potentially suitable sites. Even should these sites be presently empty, they may be occupied in the near future.
The distance to bodies of water was found to be an important factor, both in the field and through statistical analysis. Rainho and Palmeirim  made similar observations in two cave-dwelling species (Rhinolophus mehelyi and Miniopterus schreibersii), for which proximity to a source of drinking water was an important factor. The Department of Environment of the Australian government also reported that flying foxes sites were usually found close to water . Several studies indicated that bats lose a significant amount of water while they are roosting, especially under conditions of low relative humidity and high temperature [51-53]. Furthermore, lactating females need more frequent drinking than non-reproductive females . Flying foxes may also need water for cooling down. Welbergen et al.  reported temperatures exceeding 42°C in January 2002 in New South Wales, Australia, causing the deaths of thousands of flying foxes from hyperthermia. The high temperature may lead flying foxes to dip their bellies into water to cool down . The maximum temperatures in central Thailand in most months are above 30°C, with temperatures of 40°C commonly recorded in April . As some roosting trees fail to protect bats from sunlight, the availability of nearby water may help those populations to resist the worst of the heat during the hottest months. Informants living nearby bat colonies suggested flying foxes may use bodies of water as landmarks for foraging. They reported flying foxes frequently flying along the river when they depart their roosting sites in the evening and when flying back along the river in the morning. Using water bodies as foraging landmarks was reported in insect-eating bats. For example, in the little brown bat, water bodies have been shown to be used as landmarks to help foraging on patches of insects found in abundance above rivers, streams, ponds, or lakes . The association that we found between flying foxes and areas located in the lowland central plain, which is surrounded by vegetation, may also simply correspond to the extensive irrigation that allows greater vegetation than elsewhere, as observed in Australia .
Further studies, focusing on the distribution, ecology, behaviors, and disease status of flying foxes should be conducted in Thailand in the central region but also elsewhere. Although the foraging plants and some of the environmental factors associated with flying fox colonies have been reported in other countries, a follow-up should be pursued in Thailand and for its singular ecologies [9,59]. Such data would be useful for conserving flying fox populations and in disease prevention and surveillance. Flying fox movements, heat relief, water usage, and other behaviors should be more fully characterized as they are likely to have impacts upon transmission patterns. For instance, during the mating season, large aggregations of individuals migrating from different sites are observed, and, as documented in Arctic waterfowl, could potentially contribute to the spread of pathogens across bat and other populations [49,60].
The SDM maps converged with the observations discussed above, showing highly suitable areas for flying foxes mainly located along riversides, in river basins in the central plain, and in areas of moderate human population density. The number of known occurrences in our study was low (n = 22) and many studies note that small sample sizes can significantly reduce the predictive potential of models [31,61-63]. Several methods have been proposed to deal with the problem [33,64-66]. While some methods are more effective at predicting species’ distributions than others, no modeling method has proven to be the best in all situations. The ensemble modeling approach used in this study appears as a way to limit the potential influence of one particular modeling method, which was found to provide good results in previous studies [38,67,68]. However, we recognize that one of the limitations of our study may be the low sample size. One option for improvement might be to pool locations from wider areas and across countries, to have a larger sample size and sets of environmental conditions.
Even more challenging than mapping the suitability for a colony is to map the suitability for NiV infection. Generally, identifying risk factors associated with the spatial distribution of disease relies on disease distribution data that are used to quantify the effect of a set of explanatory variables on the spatial distribution of a particular disease outcome . The outcome variable can be a count of disease events in a unit area or more simply a binary response indicating the presence or absence of disease at a given location. Each outcome can be used to map other areas sharing similar risk factors . However, such an approach was not possible for NiV in Thailand since no case of NiV infection in human or pig has yet been found . What we do have outside an etiological agent, in this case NiV, are susceptible hosts (bats, pigs and humans) and environments that connect hosts and the potential agent. By PSA we mapped areas where the virus’s documented reservoirs potentially coincide. The approach has not been used in epidemiological study but may be useful in the absence of disease data, as a way to spatialize disease surveillance and regionally plan livestock production. Even though it has not been used in epidemiological study and is not based on a formal statistical model, it remained useful in the present case of a disease that is absent (and hence provides no data to train a model) as a way to integrate different factors in a risk map that can inform further planning and disease surveillance in a context of very limited knowledge. A limitation of the approach is, however, the somewhat arbitrary choices on weights that are made along the process, that are defined in a more explicit and thorough way in using MCDA approaches. Ultimately, the spatial validity of both approaches could only be formally evaluated in retrospect, if NiV infection were eventually identified in the country.
Broad-scale delineation of areas where three potential host types—bat, pig and human—are present could improve NiV surveillance strategy . Indeed, in a context of limited financial support for animal disease surveillance systems, a more optimal use of resources could be implemented if active surveillance is targeted at higher-risk farms or areas . One approach could circle around developing passive and active surveillance programs on pig farms of predicted risk, for example, with particular focus on farms of low biosecurity, nearby bodies of water, and/or hosting orchards as additional risk factors . The surveillance program should be integrated with those for other diseases to reduce cost and manpower. Simultaneously, such surveillance efforts could be reinforced with enhanced communication on good farm management practices and public awareness campaigns.
In addition, preventing direct transmission of NiV from bats to humans could be adapted to the characteristic habitats identified in this study. For instance, it is apparent that flying foxes on the central plain of Thailand are found in particular conditions in spatial (e.g., distance to water, vegetation) and social terms (e.g., undisturbed environment and community). An active surveillance program could be conducted on the people who live closely to flying fox colonies. A new colony detected in 2011 (no. 17) is surrounded by commercial orchards, in particular coconut trees . Testing NiV in a fresh coconut-palm sugar, usually produced by leaving a container on the trees overnight, may be useful for a focal study. In Bangladesh, sap harvesters were encouraged to use bamboo skirts on their trees to prevent contacts between fruit bats and raw date palm sap. Authorities educated locals to avoid drinking raw date palm sap or eat partially eaten fruit, and these efforts could be adapted for Thailand . Finally, the central plain of Thailand is an area with intense farming activities, including pig husbandry, reflecting strongly the convergence across multiple risk models here. Surveillance programs in pigs and humans should be integrated to mutually increase their effectiveness.
Part of this work was supported through the NIH NIAID grant (1R01AI101028-01A1). We thank the staff of the Department of Livestock Development (DLD) for conducting the surveillance on the animals; National Institute of Animal Health (NIAH) and Regional Veterinary Research and Development Centers for evaluating the NiV; Ministry of Transportation for geodata; Department of Provincial Administration, Ministry of Interior for population data; and Department of National Parks, Wildlife and Plant Conservation, Ministry of Resources and Environment for location data for the fruit bats. We also thank colleagues at Lutte biologique et Ecologie spatiale (LUBIES), ULB, Belgium, for assistance and suggestions.
- Suzán G, Marcé E, Giermakowski JT, Armién B, Pascale J, Mills J, et al. The effect of habitat fragmentation and species diversity loss on hantavirus Prevalence in Panama. Ann N Y Acad Sci. 2008;1149:80–3.View ArticlePubMedGoogle Scholar
- Mohd Nor MN, Gan CH, Ong BL. Nipah virus infection of pigs in peninsular Malaysia. Rev Sci Tech Int Off Epizoot. 2000;19:160–5.Google Scholar
- WHO | Nipah Virus (NiV) Infection [http://www.who.int/csr/disease/nipah/en/]
- Chua KB, Lek Koh C, Hooi PS, Wee KF, Khong JH, Chua BH, et al. Isolation of Nipah virus from Malaysian Island flying-foxes. Microbes Infect. 2002;4:145–51.View ArticlePubMedGoogle Scholar
- Mohd Yob J, Hume F, Azmin Mohd R, Christopher M, Van Der Heide B, Paul R, et al. Nipah virus infection in bats (order Chiroptera) in peninsular Malaysia. Emergin Infect Dis. 2001;7:439–41.View ArticleGoogle Scholar
- Reynes J-M, Counor D, Ong S, Faure C, Seng V, Molia S, et al. Nipah virus in Lyle’s flying foxes, Cambodia. Emerg Infect Dis. 2005;11:1042–7.View ArticlePubMed CentralPubMedGoogle Scholar
- Wacharapluesadee S, Lumlertdacha B, Boongird K, Wanghongsa S, Chanhome L, Rollin P, et al. Bat Nipah Virus, Thailand. Emerg Infect Dis. 2005;11:1949–51.View ArticlePubMed CentralPubMedGoogle Scholar
- Iehlé C, Razafitrimo G, Razainirina J, Andriaholinirina N, Goodman SM, Faure C, et al. Henipavirus and tioman virus antibodies in pteropodid bats, Madagascar. Emerg Infect Dis. 2007;13:159–61.View ArticlePubMed CentralPubMedGoogle Scholar
- DEPI - Flying-foxes [http://www.depi.vic.gov.au/environment-and-wildlife/wildlife/flying-foxes]
- Flying-foxes and national environmental law [http://www.environment.gov.au/node/16394]
- Chua K, Chua B, Wang C. Anthropogenic deforestation, El Niño and the emergence of Nipah virus in Malaysia. Malays J Pathol. 2001;24(1):15–21.Google Scholar
- Tan K-S, Tan C-T, Goh K-J. Epidemiological aspects of Nipah virus infection. Neurol J Southeast Asia. 1999;4:77–81.Google Scholar
- Chua KB, Goh KJ, Wong KT, Kamarulzaman A, Tan PSK, Ksiazek TG, et al. Fatal encephalitis due to Nipah virus among pig-farmers in Malaysia. Lancet. 1999;354:1257–9.View ArticlePubMedGoogle Scholar
- Paton NI, Leo YS, Zaki SR, Auchus AP, Lee KE, Ling AE, et al. Outbreak of Nipah-virus infection among abattoir workers in Singapore. Lancet. 1999;354:1253–6.View ArticlePubMedGoogle Scholar
- Chew MHL, Arguin PM, Shay DK, Goh K-T, Rollin PE, Shieh W-J, et al. Risk factors for Nipah virus infection among abattoir workers in Singapore. J Infect Dis. 2000;181:1760–3.View ArticlePubMedGoogle Scholar
- Hsu VP, Hossain MJ, Parashar UD, Ali MM, Ksiazek TG, Kuzmin I, et al. Nipah virus encephalitis reemergence, Bangladesh. Emerg Infect Dis. 2004;10:2082–7.View ArticlePubMed CentralPubMedGoogle Scholar
- Luby S, Rahman M, Hossain M, Blum L, Husain M, Gurley E, et al. Foodborne transmission of Nipah virus, Bangladesh. Emerg Infect Dis. 2006;12:1888–94.View ArticlePubMed CentralPubMedGoogle Scholar
- Luby SP, Gurley ES, Hossain MJ. Transmission of human infection with Nipah virus. Clin Infect Dis. 2009;49:1743–8.View ArticlePubMed CentralPubMedGoogle Scholar
- Luby SP, Hossain MJ, Gurley ES, Ahmed B-N, Banu S, Khan SU, et al. Recurrent zoonotic transmission of Nipah virus into humans, Bangladesh, 2001–2007. Emerg Infect Dis. 2009;15:1229–35.View ArticlePubMed CentralPubMedGoogle Scholar
- Chadha MS, Comer JA, Lowe L, Rota PA, Rollin PE, Bellini WJ, et al. Nipah virus-associated encephalitis outbreak, Siliguri, India. Emerg Infect Dis. 2006;12:235–40.View ArticlePubMed CentralPubMedGoogle Scholar
- Pathchimasiri T, Kalpravidh W, Damrongwatanapokin S, Chantamaneechote T, Daniels P, Buranathai C: Immunohistochemistry Investigation of Nipah Virus : A Retrospective study in Thailand. In The 11th International Symposium of the World Association of Veterinary Laboratory Diagnosticials and OIE Seminar on Biotechnology. Thai Association of Veterinary Laboratory Diagnosticians; 2003:44–45.Google Scholar
- Department of Livestock Development. Animal health in Thailand 2011. Bangkok, Thailand: Department of Livestock Development; 2012.Google Scholar
- OIE World Animal Health Information System [http://www.oie.int/wahis_2/public/wahid.php/Countryinformation/Animalsituation]
- Wacharapluesadee S, Boongird K, Wanghongsa S, Ratanasetyuth N, Supavonwong P, Saengsen D, et al. A longitudinal study of the prevalence of Nipah virus in Pteropus lylei bats in Thailand: evidence for seasonal preference in disease transmission. Vector-Borne Zoonotic Dis. 2010;10:183–90.View ArticlePubMedGoogle Scholar
- Boonkird K, Wanghongsa S. On the population number and distribution fo flying foxes (Pteropus lylei) in central plain. Wildl Yearb. 2004;5:89–100.Google Scholar
- Sedwisai P, Changbunjong T, Chamsai T, Yongyuttawichai P, Sangkachai N, Weluwanarak T, et al. The distribution of flying fox (Pteropus spp.) in the central region of Thailand. J Appl Anim Sci. 2011;4:21–9.Google Scholar
- Gaughan AE, Stevens FR, Linard C, Jia P, Tatem AJ. High resolution population distribution maps for Southeast Asia in 2010 and 2015. PLoS One. 2013;8:e55882.View ArticlePubMed CentralPubMedGoogle Scholar
- CGIAR-CSI SRTM 90 m DEM Digital Elevation Database [http://srtm.csi.cgiar.org/]
- Classification References (Using ENVI) | Exelis VIS Docs Center [http://www.exelisvis.com/docs/ClassificationReferences.html]
- Elith J, Leathwick JR. Species distribution models: ecological explanation and prediction across space and time. Annu Rev Ecol Evol Syst. 2009;40:677–97.View ArticleGoogle Scholar
- McPherson JM, Jetz W, Rogers DJ. The effects of species’ range sizes on the accuracy of distribution models: ecological phenomenon or statistical artefact?: Species’ range and distribution model accuracy. J Appl Ecol. 2004;41:811–23.View ArticleGoogle Scholar
- Elith J, Graham CH, Anderson RP, Dudík M, Ferrier S, Guisan A, et al. Novel methods improve prediction of species’ distributions from occurrence data. Ecography. 2006;29:129–51.View ArticleGoogle Scholar
- Barbet-Massin M, Jiguet F, Albert CH, Thuiller W. Selecting pseudo-absences for species distribution models: how, where and how many? Methods Ecol Evol. 2012;3:327–38.View ArticleGoogle Scholar
- Phillips SJ, Anderson RP, Schapire RE. Maximum entropy modeling of species geographic distributions. Ecol Model. 2006;190:231–59.View ArticleGoogle Scholar
- Elith J, Leathwick JR, Hastie T. A working guide to boosted regression trees. J Anim Ecol. 2008;77:802–13.View ArticlePubMedGoogle Scholar
- Breiman L. Random forests. In: Machine learning, vol. 45. The Netherlands: Kluwer Academic Publishers; 2001. p. 5–32.Google Scholar
- Cutler DR, Edwards TC, Beard KH, Cutler A, Hess KT, Gibson J, et al. Random forests for classificaiton in ecology. Ecology. 2007;88:2783–92.View ArticlePubMedGoogle Scholar
- Marmion M, Parviainen M, Luoto M, Heikkinen RK, Thuiller W. Evaluation of consensus methods in predictive species distribution modelling. Divers Distrib. 2009;15:59–69.View ArticleGoogle Scholar
- Araujo M, New M. Ensemble forecasting of species distributions. Trends Ecol Evol. 2007;22:42–7.View ArticlePubMedGoogle Scholar
- Robert J. Hijmans, Jane Elith: Species distribution modeling with R. 2013.Google Scholar
- Engler R, Waser LT, Zimmermann NE, Schaub M, Berdos S, Ginzler C, et al. Combining ensemble modeling and remote sensing for mapping individual tree species at high spatial resolution. For Ecol Manag. 2013;310:64–73.View ArticleGoogle Scholar
- Thuiller W, Lafourcade B, Engler R, Araújo MB. BIOMOD - a platform for ensemble forecasting of species distributions. Ecography. 2009;32:369–73.View ArticleGoogle Scholar
- Augustin NH, Mugglestone MA, Buckland ST. An autologistic model for the spatial distribution of wildlife. J Appl Ecol. 1996;33:339–47.View ArticleGoogle Scholar
- Fielding AH, Bell JF. A review of methods for the assessment of prediction errors in conservation presence/absence models. Environ Conserv. 1997;24:38–49.View ArticleGoogle Scholar
- Stevens KB, Gilbert M, Pfeiffer DU. Modeling habitat suitability for occurrence of highly pathogenic avian influenza virus H5N1 in domestic poultry in Asia: a spatial multicriteria decision analysis approach. Spat Spatio-Temporal Epidemiol. 2013;4:1–14.View ArticleGoogle Scholar
- Ano T, Hormwichian R, Jitrapinate N, Compliew S, Kangrang A. The estimation of drought risk area using potential surface analysis technique. UBU Eng J. 2013;6:13–21.Google Scholar
- Udomsap T, Iamtrakul P. Accessibility improvement for district’s urban diversity: case study of Rachadamnoen klong avenue, Bangkok. J Soc Transp Traffic Stud JSTS. 2011;2:1–17.Google Scholar
- Nakya S, Leopairojna SK, Rangsiraksa L. Use of satellite data and potential surface analysis for Urban expansion of Hua Hin Municipality, Prachuap Khiri Khan Province, Thailand. In: 31st Asian conference on remote sensing 2010. Volume 1. Hanoi, Vietnam: Asian Association on Remote Sensing; 2010. p. 293–302.Google Scholar
- Eby P. Seasonal movements of grey-headed flying-foxes, Pteropus poliocephalus (Chiroptera : Pteropodidae), from two maternity camps in northern New South Wales. Wildl Res. 1991;18:547.View ArticleGoogle Scholar
- Rainho A, Palmeirim JM. The importance of distance to resources in the spatial modelling of bat foraging habitat. PLoS One. 2011;6:e19227.View ArticlePubMed CentralPubMedGoogle Scholar
- Webb PI. The comparative ecophysiology of water balance in microchiropteran. 67;1995:203–218Google Scholar
- Webb PI, Speakman JR, Racey PA. Evaporative water loss in two sympatric species of vespertilionid bat, Plecotus auritus and Myotis daubentoni: relation to foraging mode and implications for roost site selection. J Zool. 2009;235:269–78.View ArticleGoogle Scholar
- Studier EH, O’Farrell MJ. Biology of Myotis thysanodes and M. lucifugus (Chiroptera: Vespertilionidae)—III. Metabolism, heart rate, breathing rate, evaporative water loss and general energetics. Comp Biochem Physiol A Physiol. 1976;54:423–32.View ArticleGoogle Scholar
- Adams RA, Hayes MA. Water availability and successful lactation by bats as related to climate change in arid regions of western North America. J Anim Ecol. 2008;77:1115–21.View ArticlePubMedGoogle Scholar
- Welbergen JA, Klose SM, Markus N, Eby P. Climate change and the effects of temperature extremes on Australian flying-foxes. Proc R Soc B Biol Sci. 2008;275:419–25.View ArticleGoogle Scholar
- Killer climate: tens of thousands of flying foxes dead in a day [http://www.brisbanetimes.com.au/queensland/killer-climate-tens-of-thousands-of-flying-foxes-dead-in-a-day-20140225-33drr.html]
- Thai Meteorological Department [http://www.tmd.go.th/en/]
- Clare EL, Barber BR, Sweeney BW, Hebert PDN, Fenton MB. Eating local: influences of habitat on the diet of little brown bats (Myotis lucifugus): molecular detection of variation in diet. Mol Ecol. 2011;20:1772–80.View ArticlePubMedGoogle Scholar
- Flying foxes | NSW Environment & Heritage [http://www.environment.nsw.gov.au/animals/flyingfoxes.htm]
- Olsen B. Global patterns of influenza A virus in wild birds. Science. 2006;312:384–8.View ArticlePubMedGoogle Scholar
- Stockwell DR, Peterson AT. Effects of sample size on accuracy of species distribution models. Ecol Model. 2002;148:1–13.View ArticleGoogle Scholar
- Chen H, Wood MD, Linstead C, Maltby E. Uncertainty analysis in a GIS-based multi-criteria analysis tool for river catchment management. Environ Model Softw. 2011;26:395–405.View ArticleGoogle Scholar
- Hirzel A, Guisan A. Which is the optimal sampling strategy for habitat suitability modelling. Ecol Model. 2002;157:331–41.View ArticleGoogle Scholar
- Segurado P, Araújo MB. An evaluation of methods for modelling species distributions: methods for modelling species distributions. J Biogeogr. 2004;31:1555–68.View ArticleGoogle Scholar
- Engler R, Guisan A, Rechsteiner L. An improved approach for predicting the distribution of rare and endangered species from occurrence and pseudo-absence data. J Appl Ecol. 2004;41:263–74.View ArticleGoogle Scholar
- Hanberry BB, He HS, Palik BJ. Pseudoabsence generation strategies for species distribution models. PLoS One. 2012;7:e44486.View ArticlePubMed CentralPubMedGoogle Scholar
- Lauzeral C, Grenouillet G, Brosse S. Dealing with noisy absences to optimize species distribution models: an iterative ensemble modelling approach. PLoS One. 2012;7:e49508.View ArticlePubMed CentralPubMedGoogle Scholar
- Grenouillet G, Buisson L, Casajus N, Lek S. Ensemble modelling of species distribution: the effects of geographical and environmental ranges. Ecography. 2011;34:9–17.View ArticleGoogle Scholar
- Graham A, Atkinson P, Danson F. Spatial analysis for epidemiology. Acta Trop. 2004;91:219–25.View ArticlePubMedGoogle Scholar
- Katharina DS, Regula G, Hernandez J, Knopf L, Fuchs K, Morris RS, et al. Concepts for risk-based surveillance in the field of veterinary medicine and veterinary public health: Review of current approaches. BMC Health Serv Res 2006:6–20Google Scholar
- Khan SU, Gurley ES, Hossain MJ, Nahar N, Sharker MAY, Luby SP. A randomized controlled trial of interventions to impede date palm sap contamination by bats to prevent Nipah virus transmission in Bangladesh. PLoS One. 2012;7:e42689.View ArticlePubMed CentralPubMedGoogle Scholar
- Forest Types in Thailand [http://wildlifethailand.com]
- Mangrove Forest Habitat in Tropical Thailand [http://www.kohphrathong.com/thailand_mangrove.html]
- Asia-Pacific Forestry Sector Outlook Study: Country Report - Thailand [http://www.fao.org/docrep/003/x2649e/X2649E03.htm]
- Biodiversity & Expert Database [http://www.bedo.or.th/lcdb/default.aspx]
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. 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.