Mapping risk of plague in Qinghai-Tibetan Plateau, China
- Quan Qian†1, 2,
- Jian Zhao†3,
- Liqun Fang1,
- Hang Zhou3,
- Wenyi Zhang1,
- Lan Wei1,
- Hong Yang1,
- Wenwu Yin3,
- Wuchun Cao1 and
- Qun Li3Email author
© Qian et al.; licensee BioMed Central Ltd. 2014
Received: 8 April 2014
Accepted: 3 July 2014
Published: 10 July 2014
Qinghai-Tibetan Plateau of China is known to be the plague endemic region where marmot (Marmota himalayana) is the primary host. Human plague cases are relatively low incidence but high mortality, which presents unique surveillance and public health challenges, because early detection through surveillance may not always be feasible and infrequent clinical cases may be misdiagnosed.
Based on plague surveillance data and environmental variables, Maxent was applied to model the presence probability of plague host. 75% occurrence points were randomly selected for training model, and the rest 25% points were used for model test and validation. Maxent model performance was measured as test gain and test AUC. The optimal probability cut-off value was chosen by maximizing training sensitivity and specificity simultaneously.
We used field surveillance data in an ecological niche modeling (ENM) framework to depict spatial distribution of natural foci of plague in Qinghai-Tibetan Plateau. Most human-inhabited areas at risk of exposure to enzootic plague are distributed in the east and south of the Plateau. Elevation, temperature of land surface and normalized difference vegetation index play a large part in determining the distribution of the enzootic plague.
This study provided a more detailed view of spatial pattern of enzootic plague and human-inhabited areas at risk of plague. The maps could help public health authorities decide where to perform plague surveillance and take preventive measures in Qinghai-Tibetan Plateau.
Plague, caused by Yersinia pestis, is a rapidly progressing, highly infectious, and highly feared disease that is likely to be fatal without prompt antibiotic treatment . Plague is maintained among wild rodents in distinct geographic foci showing a serious threat to humans in Asia, America, and Africa . Although the mechanisms by which plague is maintained between epizootic cycles are not well understood, it is generally accepted that the disease cycles between enzootic infections and occasional epizootic outbreaks among susceptible hosts . Human plague cases are relatively low incidence but high mortality, which presents unique surveillance and public health challenges, because early detection through surveillance may not always be feasible and infrequent clinical cases may be misdiagnosed.
In mainland China, natural plague foci of plague were divided into 12 types according to their primary reservoirs, principal vectors, landscapes and genotypes of Y. pestis , which were distributed in 19 of total 31 provinces and autonomous regions., The marmots (Marmota himalayana) foci in Qinghai-Tibetan Plateau was most active and widespread. The primary host, M. himalayana usually hibernates from October and comes out of hibernation from April the next year . Callopsylla dolabris and Oropsylla silantiewi are known as the principal vectors . The epizootics of plague occurs during the whole period when the rodent hosts are active on the ground, and usually reaches the peak in June and July . Human cases occurred almost every year in the enzootic region . In 2009, an outbreak involving 12 pneumonic plague cases occurred in a remote village of Qinghai Province leading to three deaths . In 2010, cases with bubonic plague and pneumonic plague were reported in Gansu Province and Tibet Autonomous Region.
In mainland China, the natural foci of plague had been discovered and described since 1950, which were plotted out based on administrative boundary. This approach is rather rough and enlarged the potential risk area. Natural foci of plague are related to the particular landscapes which are favorable for a high and stable number of rodent reservoirs and flea vectors of Y. pestis [8–10]. Thus, understanding how these host species are geographically distributed in China and relating that information to how effciently each can transmit plague is critical to under-standing the ecology of plague and to recognizing where the greatest threats for people exist in the country.
The “Host Niche Hypothesis” (HNH) postulates that plague distributions are mediated by host distributions, such that the distribution of plague depends on an amalgam of host ranges, and the presence of a particular host species could extend the distributional potential of the pathogen . Thus, by predicting the potential distribution of host animal, the plague risk could be mapped. Topography, vegetation, climate and other environmental factors are thought to influence spatial distibution and temporal dynamics of the host animal [12–18]. Novel spatial modelling methods such as maximum entropy (MAXENT) and the genetic algorithm for rule set production (GARP) require only disease presence data and have been used extensively in the fields of ecology and conservation, to model species distribution and habitat suitability . In this paper, we aimed to predict the potential natural foci of plague in Qinghai-Tibetan Plateau by ecological niche modeling based on environmental parameters derived from remote sensing data and the transmission risk to support our control and prevention.
Data collection and pre-processing
All plague surveillance data used in this study were derived from the plague surveillance database of ‘Chinese Plague Prevention and Control Management Information System’, a national surveillance network for plague maintained by Chinese Center for Disease Control and Prevention. According to the surveillance from 2007 to 2009, 233 infected himalayan marmots were confirmed by bacteriological tests including microscopy, culture, phagelysis, and mouse inoculation, geocoded by globe positioning system (GPS) or detail address. In addition, 51 human plague cases in the study area from 2004-2010 were collected from Chinese Center for Disease Control and Prevention. As the spatial resolution of enrironmental data was 1 km, we removed duplicated records of the same pixel and got 218 spatially unique points with infected marmots with a spatial precision of 5–10 km (~0.05–0.1°).
Environmental variables supposed to be associated with the existence of natural foci of plague were collected. The distribution of marmot was mainly driven by geomorphology, climate and vegetation . Geomorphology variables, as Elevation, slope, aspect, were derived from SRTM30 (Available from: http://www2.jpl.nasa.gov/srtm/). Climate factors were characterized by Annual average daytime land surface temperature (LSTd) and annual average nighttime land surface temperature (LSTn), composited from eight days MODIS LST products (MOD11A2) from 2004-2009 (Available from: https://lpdaac.usgs.gov/lpdaac/products/modis_products_table). The amount and productivity of vegetation was featured by annual maximum normalized difference of vegetation index (NDVI-) from ten days SPOT-VGT S10 products from 2004-2009 (Available from: http://free.vgt.vito.be/). In addition, land use and cover change (LUCC) revealed the proper grassland and effect of human activity, which was derived from GlobCover Land Cover version V2.2 released by European Space Agency (Available from: http://www.esa.int/esaEO/SEMXB7TTGOF_index_0.html). In addition, Chinese population distribution raster map with 1 km2 resolution was used to identify human-inhabited areas with high plague risk by overlayanalysis. This map was obtained from Data Center for Resources and Environmental Sciences Chinese Academy of Sciences.
All data layers were projected in Albers coordinates and generalized to a pixel resolution of ~1 × 1 km for analysis in ENVI 4.5 (ITT Visual Information Solutions, Boulder, CO, USA) and ArcGIS Desktop 9.2 (Environmental Systems Research Institute, Redlands, CA, USA) environment. We also considered the spatial heterogeneity in the whole china using the GeoDetector software [23, 24]. The power of determinant (PD) reflects the degree to which a determinant explains the prevalence of the disease.
Ecological niche modeling
Recently ecological niche modeling (ENM) methods, such as GARP  and Maxent , have been widely applied in species distribution modeling when only presence data is available for prediction [27, 28]. ENM also has been proved to be effective of in applications to questions regarding the geography and ecology of disease transmission [29–31]. In this study, presence-only method Maxent was applied [32, 33]. Maxent is a general-purpose method for characterizing probability distributions from incomplete information based on the principle of maximum entropy .
Maxent outputs the maximum entropy distribution that satisfies a set of environmental constraints. In place of true absences, Maxent uses background points (pseudo-absences) to evaluate commission . We run the model in the support of the “Maxent” Version 3.3.1 k software . Within Maxent processing, 75% occurrence points (164) were randomly selected for training model, and the rest 25% points (54) were used for model test and validation. 11000 ‘pseudo-absence’ points created by random sampling from areas lacking known presences. To measure the relative contribution of each environmental variable to the predictive model, a jackknife manipulation was performed. Receiver operating characteristic (ROC) analysis was used to evaluate the discrimination ability of models, and to determine the optimal probability cut-off value for classifying the risk areas of plague. We measured Maxent model performance as test gain and test AUC (random 25% testing) for all analyses omitting individual variable, and for each individual variable alone. The AUC is considered as an effective indicator of model performance. The larger the AUC, the highest is the sensitivity rate and the lower is the 1-specificity rate. Usually AUC values of 0 · 5–0 · 7 are taken to indicate low accuracy, values of 0 · 7–0 · 9 indicate useful applications and values of > 0 · 9 indicate high accuracy . The optimal probability cut-off value was chosen by maximizing training sensitivity and specificity simultaneously . All pixels with a probability value at least equal to the optimal value were classified as plague risk area, then we used binomial test based on omission of independent test points and predicted area to determine whether a model predicts the test localities significantly better than random .
Summary of ‘Jackknife analysis ’ used to determine importance of each environmental variable
National surveillance network for plague recorded the information of infected marmots as the primary data source of our study. We hypothesis that these geocoded points represented the suitable habitat for plague host animal and fleas that could spread Yersinia pestis. In the niche modeling, the relevant environmental variables reflects the distribution of suitable areas of host animals, rather than plague foci or Yersinia pestis. This plague-positive host animal extent would be smaller than the whole species on the landscape, but we could not tell the exact proportion or overlap between them because of presence data limitation.
In this study, we applied field surveillance data with GIS/RS-based ecological niche modeling approach to increase our understanding of the suitable habitat of enzootic plague and anticipate the distribution of the natural foci in Qinghai-Tibetan Plateau. The reliable absence data are hardly available in enzootic plague surveillance as the plague-enzootic foci are easily missed during field surveys, thus the “presence-only” ecological niche modeling Maxent was preferred. Independent tests indicated ENM can anticipate the distribution of host animal robustly.
Environmental variables, such as topography, climate and vegetation, can play important roles in natural focal diseases by affecting pathogens directly, or by influencing the distribution and abundance of disease hosts and vectors [37, 38]. In this study, all remotely sensed environmental variables appear to contribute positively to the final predictive model. Elevation, LSTd, LSTn and NDVI are likely to be most important factors which influence the distribution of host marmots. The environmental ranges of host marmots were clear revealed in the two-dimensional scatter plots. The environmental conditions of host marmots are proved to be highly diverse. The elevation range of host marmots are various, consistent with Qinghai-Tibetan Plateau . The risk increases up to an elevation of 3177 m and then declines as elevation increase, because in the lower height, forest and thickly grass would be harmful for marmot making doggishness, while in the higher height, there would be lack of food. The marmot preferred temperature between 5°C and 20°C in the daylight. The extremely low temperature below -8.38°C would limit their activities. Extremely low and high temperature of land surface seems negatively affect the distribution of host animals. The temperature of land surface may influence the plague persistence in ‘marmot-flea’ communities in complex ways. We suspected extreme temperature may negatively affect marmot and flea ecology, and blockage of vector fleas . Host marmots were observed and predicted in the sparely vegetated areas, but in our modeling, NDVI variation had no significant limitation in host animal distribution (Figure 3). This result suggested that plague prefer alpine desert/semi-desert grasslands confirmed previously, such as alpine meadow, alpine grassland and alpine shrub [5, 39], while habitats with extravagant vegetation seemed to be unsuitable for enzootic plague. This could be due to that luxuriant vegetation obstructs them to protect from predators, and the high subterranean biomass make them hard to burrow underground .
Maxent outputs the distribution probability of the host marmots, which was segmented by maximizing training sensitivity and specificity simultaneously. The binarization map reveal massive extent of natural foci of plague in Qinghai-Tibetan Plateau, 16.5% (575,355 km2) of the study area was defined as natural foci. Since the maps show where the enzootic plague are more likely to occur and where natural foci probably exist, they can be used to guide surveillance of plague in future.
Generally speaking, human plague cases are directly or indirectly associated with the epizootic plague activity. In Qinghai-Tibetan Plateau, most human plague outbreaks stemmed from contact with infected marmots. We estimated more than 2,000,000 people at the risk of exposure to the enzootic plague. Most risk population are distributed in the east and south of Qinghai-Tibetan Plateau. Our findings can guide the design and spatial targeting of plague prevention and control efforts against human plague infection in Qinghai-Tibetan Plateau. Thereby, limit resource could be appropriately allocated to the a few areas where most human plague cases are most likely to occur.
We got a risk map that shows the full extent of areas of potential natural foci with similar environmental conditions to the observed presence infected host. However, it is important to keep in mind that some unmappable risk factors were not included in our analysis. Natural barriers and other factors may also accidentally obstruct the dispersal of marmots, vector fleas, or Y. pestis into the certain areas where habitats are suitable for persistence of natural foci. In addition, epizootic activity may be quite dynamic in the same plague focus in different periods [15, 18, 42, 43]. Temporal variation of climate and other environmental factors may drive the dynamics of epizootic activity in foci. Thus multivariate time series analysis of plague epizootic is required. Finally, with the implementation of the western China development strategy, the traffic transport between Qinghai-Tibetan Plateau and other parts of China has become much more convenient than before, especially since the completion of Qinghai-Tibet Railway in 2006. There will be more and more people such as tourists and transient workers get into Qinghai-Tibetan Plateau in future. The risk of human exposure to enzootic plague and long-distance transmission will increase consequently. Thus the risk of human plague should be further assessed under the scenario of increasing mobility of the population in Qinghai-Tibetan Plateau.
Maxent was suitable for this 1 km resolution modeling which outputs the distribution probability of the host marmots, its feasibility for other resolution would be tested in the following study. We got a more detailed view of spatial pattern of potential plague natural foci by maximizing training sensitivity and specificity simultaneously. Human-inhabited areas and population at the risk of exposure to enzootic plague could be identified by overlying Chinese population distribution map to the risk map created above. The maps could help public health authorities decide where to perform plague surveillance and take preventive measures in Qinghai-Tibetan Plateau.
This study was supported by the National Science Fund for Distinguished Young Scholars (No. 30725032), the Special Program for Prevention and Control of Infectious Diseases (Nos. 2008ZX10004-012 and 2009ZX10004-720), the Natural Science Foundation of China (No. 30590374), and the Special Fund for Health Sector, People’s Republic of China (No. 201202021). We thank Dr. Peng Gong at State Key Laboratory of Remote Sensing Science, Jointly Sponsored by Institute of Remote Sensing Applications, Chinese Academy of Sciences and Beijing Normal University; Dr. Sake J. de Vlas and Dr. Jan Hendrik Richardus at Department of Public Health of Erasmus Medical Center (MC), University Medical Center Rotterdam, for their comments on draft of the paper. We thank Dr. Feng-ming Hui at College of Global Change and Earth System Science, Beijing Normal University, for valuable suggestions of data analysis.
- Stenseth N, Atshabar B, Begon M, Belmain S, Bertherat E, Carniel E, Gage K, Leirs H, Rahalison L: Plague: past present and future. PLoS Med. 2008, 5: E3-View ArticlePubMedPubMed CentralGoogle Scholar
- Tan J, Liu Y, Shen E, Zhu W, Wang W, Li R, Yang L: Plague Foci of the World. The Atlas of Plague and Its Environment in the People’s Republic of China. 2000, Beijing: Science press, 42-46.Google Scholar
- Gage K, Kosoy M: Natural history of plague: perspectives from more than a century of research. Annu Rev Entomol. 2005, 50: 505-528. 10.1146/annurev.ento.50.071803.130337.View ArticlePubMedGoogle Scholar
- Cong X, Yin W: Plague Overview. Chinese Plague Manual: Prevention, Control and Emergency Response. 2009, Beijing: Beijing University Medical Press, 1-49.Google Scholar
- Tian F: Investigation of the nature focus of Marmota himalayana [in Chinese]. Chin J Zoonoses. 2000, 16: 95-97.Google Scholar
- Yin W: Analysis and Prediction of Plague Epidemic of 2009 [in Chinese]. China National Plague Surveillance Summary of China CDC. 2009, 32-38.Google Scholar
- World Health Organization: Plague in China. http://www.who.int/csr/don/2009_08_11/en/index.html,
- Tikhomirov E: Epidemiology and Distribution of Plague. Plague Manual: Epidemiology, Distribution, Surveillance and Control. http://www.who.int/entity/csr/resources/publications/plague/whocdscsredc992a.pdf,
- Davis S, Trapman P, Leirs H, Begon M, Heesterbeek J: The abundance threshold for plague as a critical percolation phenomenon. Nature. 2008, 454: 634-637. 10.1038/nature07053.View ArticlePubMedGoogle Scholar
- Davis S, Begon M, De Bruyn L, Ageyev V, Klassovskiy N, Pole S, Viljugrein H, Stenseth N, Leirs H: Predictive thresholds for plague in Kazakhstan. Science. 2004, 304: 736-10.1126/science.1095854.View ArticlePubMedGoogle Scholar
- Maher S, Ellis C, Gage K, Enscore R, Peterson A: Range-wide determinants of plague distribution in North America. Am J Trop Med Hyg. 2010, 83 (4): 736-742. 10.4269/ajtmh.2010.10-0042.View ArticlePubMedPubMed CentralGoogle Scholar
- Neerinckx S, Peterson AT, Gulinck H, Deckers J, Leirs H: Geographic distribution and ecological niche of plague in sub-Saharan Africa. Int J Health Geogr. 2008, 7: 54-10.1186/1476-072X-7-54.View ArticlePubMedPubMed CentralGoogle Scholar
- Holt A, Salkeld D, Fritz C, Tucker J, Gong P: Spatial analysis of plague in California: niche modeling predictions of the current distribution and potential response to climate change. Int J Health Geogr. 2009, 8: 38-10.1186/1476-072X-8-38.View ArticlePubMedPubMed CentralGoogle Scholar
- Adjemian J, Girvetz E, Beckett L, Foley J: Analysis of Genetic Algorithm for Rule-Set Production (GARP) modeling approach for predicting distributions of fleas implicated as vectors of plague, Yersinia pestis, in California. J Med Entomol. 2006, 43: 93-103. 10.1603/0022-2585(2006)043[0093:AOGAFR]2.0.CO;2.PubMedGoogle Scholar
- Eisen R, Reynolds P, Ettestad P, Brown T, Enscore R, Biggerstaff B, Cheek J, Bueno R, Targhetta J, Montenieri J, Gage K: Residence-linked human plague in New Mexico: a habitat-suitability model. Am J Trop Med Hyg. 2007, 77: 121-PubMedGoogle Scholar
- Eisen R, Glass G, Eisen L, Cheek J, Enscore R, Ettestad P, Gage K, Ettestad P, Gage K: A spatial model of shared risk for plague and hantavirus pulmonary syndrome in the southwestern United States. Am J Trop Med Hyg. 2007, 77: 999-PubMedGoogle Scholar
- Winters A, Staples J, Ogen-Odoi A, Mead P, Griffith K, Owor N, Babi N, Enscore R, Eisen L, Gage K, Eisen R: Spatial risk models for human plague in the West Nile region of Uganda. Am J Trop Med Hyg. 2009, 80: 1014-PubMedGoogle Scholar
- Stenseth N, Samia N, Viljugrein H, Kausrud K, Begon M, Davis S, Leirs H, Dubyanskiy V, Esper J, Ageyev V, Klassovskiy N, Pole S, Chan K: Plague dynamics are driven by climate variation. Proc Natl Acad Sci U S A. 2006, 103: 13110-10.1073/pnas.0602447103.View ArticlePubMedPubMed CentralGoogle Scholar
- Stevens K, Pfeiffer D: Spatial modelling of disease using data- and knowledge-driven approaches. Spat Spatio-tempor Epidemiol. 2011, 2 (3): 125-133. 10.1016/j.sste.2011.07.007.View ArticleGoogle Scholar
- Wang Z, Chen H, Zheng Y, Wu K: Plague Epizoology. Types and Distribution of Plague Natural Foci in China. 2009, Beijing: Science Press, 264-Google Scholar
- Zhao S: A new proposal of Chinese physical geographical zoning. J Geogr Sci. 1983, 38: 1-10.Google Scholar
- Gao M: Spatial Distribution Prediction and Analysis of Himalayan Marmot Plague Natural Foci in Qinghai-Tibet Plateau. 2011, Beijing China: Graduate University of Chinese Academy SciencesGoogle Scholar
- Wang J, Li X, Christakos G, Liao Y, Zhang T, Gu X, Zheng X: Geographical detectors-based health risk assessment and its application in the neural tube defects study of the Heshun region, China. Int J Geogr Inf Sci. 2010, 24 (1): 107-127. 10.1080/13658810802443457.View ArticleGoogle Scholar
- Wang J, Hu Y: Environmental health risk detection with GeogDetector. Environ Model Software. 2012, 33: 114-115.View ArticleGoogle Scholar
- Pereira R: Desktop Garp. Lawrence, Kansas: University of Kansas Biodiversity Research Center, http://www.nhm.ku.edu/desktopgarp/,
- Phillips S, Anderson R, Schapire R: Maximum entropy modeling of species geographic distributions. Ecol Model. 2006, 190: 231-259. 10.1016/j.ecolmodel.2005.03.026.View ArticleGoogle Scholar
- Raxworthy C, Martinez-Meyer E, Horning N, Nussbaum R, Schneider G, Ortega-Huerta M, Peterson T: Predicting distributions of known and unknown reptile species in Madagascar. Nature. 2003, 426: 837-841. 10.1038/nature02205.View ArticlePubMedGoogle Scholar
- Saatchi S, Buermann W, Ter Steege H, Mori S, Smith T: Modeling distribution of Amazonian tree species and diversity using remote sensing measurements. Remote Sens Environ. 2008, 112: 2000-2017. 10.1016/j.rse.2008.01.008.View ArticleGoogle Scholar
- Peterson A, Sánchez-Cordero V, Beard C, Ramsey J: Ecologic niche modeling and potential reservoirs for Chagas disease, Mexico. Emerg Infect Dis. 2002, 8: 662-667. 10.3201/eid0807.010454.View ArticlePubMedPubMed CentralGoogle Scholar
- Mak S, Klinkenberg B, Bartlett K, Fyfe M: Ecological Niche Modeling of Cryptococcus gattii in British Columbia, Canada. Environ Health Persp. 2009, 118: 653-10.1289/ehp.0901448.View ArticleGoogle Scholar
- Peterson A, Bauer J, Mills J: Ecologic and geographic distribution of filovirus disease. Emerg Infect Dis. 2004, 10: 40-47. 10.3201/eid1001.030125.View ArticlePubMedPubMed CentralGoogle Scholar
- Steven J, Miroslav D, Robert E: A Maximum Entropy Approach to Species Distribution Modeling. 2004, Machine Learning: In Proceedings of the Twenty-First International Conference on, 655-662.Google Scholar
- Jane E, Steven J, Trevor H, Miroslav D, Yung E, Colin J: A statistical explanation of MaxEnt for ecologists. Divers Distrib. 2011, 17: 43-57. 10.1111/j.1472-4642.2010.00725.x.View ArticleGoogle Scholar
- Zhao J, Xu M, Lu S, Cao C: Human settlement evaluation in mountain areas based on remote sensing. GIS and ecological niche modeling. J Mount Sci. 2013, 10 (3): 378-387. 10.1007/s11629-013-2413-2.View ArticleGoogle Scholar
- Guisan A, Zimmermann N: Predictive habitat distribution models in ecology. Ecol Model. 2000, 135: 147-186. 10.1016/S0304-3800(00)00354-9.View ArticleGoogle Scholar
- Anderson R, Gómez-Laverde M, Peterson A: Geographical distributions of spiny pocket mice in South America: insights from predictive models. Global Ecol Biogeogr. 2002, 11: 131-141. 10.1046/j.1466-822X.2002.00275.x.View ArticleGoogle Scholar
- Menne B, Davies C, Daniel M, Kriz B, Gratz N: Environmental changes and vector and rodent borne diseases in Europe. Epidemiology. 2004, 15: S111-View ArticleGoogle Scholar
- Walther G, Post E, Convey P, Menzel A, Parmesan C, Beebee T, Fromentin J, Guldberg O, Bairlein F: Ecological responses to recent climate change. Nature. 2002, 416: 389-395. 10.1038/416389a.View ArticlePubMedGoogle Scholar
- Tan J, Liu Y, Shen E, Zhu W, Wang W, Ribang L, Linsheng Y: Maps of Plague Natural Foci and Their Environments. The Atlas of Plague and Its Environment in the People’s Republic of China. 2000, Beijing: Science pressGoogle Scholar
- Perry R, Fetherston J: Yersinia pestis–etiologic agent of plague. Clin Microbiol Rev. 1997, 10: 35-66.PubMedPubMed CentralGoogle Scholar
- Shi S, Liu F, Yan X: Habitat selection by Marmota himalayana in the eastern Qilian mountains [in Chinese]. J Gansu Agr Univ. 2008, 43: 28-Google Scholar
- World Health Organization: Plague. http://www.who.int/csr/disease/plague/en/,
- Parmenter R, Yadav E, Parmenter C, Ettestad P, Gage K: Incidence of plague associated with increased winter-spring precipitation in New Mexico. Am J Trop Med Hyg. 1999, 61: 814-PubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2334/14/382/prepub
This article is published under license to BioMed Central Ltd. 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.