Ecology and geography of avian influenza (HPAI H5N1) transmission in the Middle East and northeastern Africa
© Williams and Peterson; licensee BioMed Central Ltd. 2009
Received: 17 March 2009
Accepted: 20 July 2009
Published: 20 July 2009
The emerging highly pathogenic avian influenza strain H5N1 ("HPAI-H5N1") has spread broadly in the past decade, and is now the focus of considerable concern. We tested the hypothesis that spatial distributions of HPAI-H5N1 cases are related consistently and predictably to coarse-scale environmental features in the Middle East and northeastern Africa.
We used ecological niche models to relate virus occurrences to 8 km resolution digital data layers summarizing parameters of monthly surface reflectance and landform. Predictive challenges included a variety of spatial stratification schemes in which models were challenged to predict case distributions in broadly unsampled areas.
In almost all tests, HPAI-H5N1 cases were indeed occurring under predictable sets of environmental conditions, generally predicted absent from areas with low NDVI values and minimal seasonal variation, and present in areas with a broad range of and appreciable seasonal variation in NDVI values. Although we documented significant predictive ability of our models, even between our study region and West Africa, case occurrences in the Arabian Peninsula appear to follow a distinct environmental regime.
Overall, we documented a variable environmental "fingerprint" for areas suitable for HPAI-H5N1 transmission.
Highly pathogenic avian influenza of the strain H5N1 (hereafter "HPAI-H5N1") has received considerable attention as an emerging virus with human pandemic potential [1, 2] since it was first shown to be the cause of human morbidity and mortality in Hong Kong in 1997 . To date, however, its most serious impacts have been on domestic poultry: millions of domestic birds have been killed by HPAI-H5N1 infection, and >230 million domestic birds have been culled to contain the spread of the virus . In contrast to the dramatic publicity, relatively few human cases are confirmed: at the time of writing, 385 human HPAI-H5N1 cases had been documented, of which 243 (63.1%) were fatal , from 60 countries . Human cases however, may eventually prove to be significantly underreported, reducing case-fatality rates.
Until Spring 2005, HPAI-H5N1 was restricted to East and Southeast Asia . Between May and June 2005, however, >6000 birds of 8 wild waterbird species were found dead at Qinghai Lake, in central China: HPAI-H5N1 was detected in 15 birds of 6 wild species , some migratory, fueling fears of broader spread . This event apparently marked a turning point in the spread of the virus: by early 2006, it had been detected widely across South Asia, Western Europe, and parts of Africa . However, whether this rapid spread resulted from accelerated dispersal or from improved surveillance detecting existing infections remains debatable .
The first Middle Eastern detection of HPAI-H5N1 was in Turkey in October 2005, in a flock of "backyard" turkeys (seeional file 1). Further detections followed in 7 Balkan countries (Bosnia-Herzegovina, Bulgaria, Croatia, Greece, Romania, Serbia and Montenegro, and Slovenia; November 2005 – March 2006), more broadly in the Middle East (Egypt, Iraq, Iran, Israel, Jordan, Kuwait, Palestinian Territories; November 2005 – March 2006), and the Caucasus (Azerbaijan and Georgia; January – February 2006) by March 2006. The virus was detected in Sudan and Djibouti in April 2006, and in Saudi Arabia in March 2007 . Countries in the region yet to record cases include the richest (Bahrain, Qatar, United Arab Emirates) and the poorest (Eritrea, Somalia, Yemen).
Most of the 9 tests conducted as part of this study indicated that independent test points coincided with ENM predictions significantly better than random expectations (see Additional file 2), although not without exceptions. In other words, in general, models based on known HPAI-H5N1 occurrences were able to anticipate spatial distributions of independent samples of HPAI-H5N1 based on their environmental attributes. The details of these test results follow.
Predictivity across study region
Single testing regions
Single predictor regions
Predictions of independent points across landscapes based on single training regions were less successful (see Additional file 2; Figure 3). Indeed, only 2 of 4 models showed any predictive ability. Predictions from northeastern Africa to the rest of the region were significant at 8 of 10 thresholds, and projections from Levant-Iran to the rest of the region were significant at 5 of 10 thresholds. Projections based on models trained in the Arabian Peninsula and Balkans-Caucasus showed no significant ability when challenged to predict occurrences in the remaining regions. Once again, visually, the Arabian Peninsula models performed particularly poorly (Figure 3).
Partial ROC analyses
The partial ROC analyses (see Additional file 2) were largely consistent with the cumulative binomial probability results (see Additional file 2). According to these tests, all single-testing-region predictions were successful (i.e., P ≤ 0.001) while 2 of 4 single predictor regions (Levant-Iran, northeastern Africa) were significantly better than random (P ≤ 0.005). The partial ROC evaluation of the overall prediction of the ProMed data was similarly significant (P < 0.01)
Our results are generally consistent with earlier predictions of the ecological niche of HPAI-H5N1 in West Africa . Most Middle Eastern and northeastern African models predicted suitable areas for HPAI-H5N1 transmission in human-habitable areas, such as the Nile Valley, the Levant, the Fertile Crescent, and the savannas of southern Sudan. The major difference between the two sets of models is that most Middle Eastern and northeastern African models predicted suitability in montane areas (Caucasus, Ethiopian Highlands, northern and western Iran, and Turkey), whereas the West African models focused prediction of suitable areas in lowlands. These models agree most clearly in implicating areas with greatest seasonal variation as representing high HPAI-H5N1 risk.
It is interesting that Arabian models should produce predictions so inconsistent with those from the rest of the study area (see Additional file 2): for example, the distribution of Arabian Peninsula occurrences could not be predicted with any confidence by models trained in the remainder of the region, and conversely, Arabian Peninsula points were unable to predict occurrences across the Balkans, Caucasus, Levant, Iran, or northeastern Africa successfully. Three major HPAI-H5N1 outbreaks occurred in the Arabian Peninsula: in Kuwait, Ar-Riyad (city), and southern Ar-Riyad (province), none of which is predicted strongly by models trained elsewhere (Figures 2 and 3). Given the rather extreme arid conditions in the region, the Arabian Peninsula seems a harsh environment for both poultry and poultry diseases. We suspect that Arabian HPAI-H5N1 occurs chiefly or only in human-subsidized habitats that would permit poultry to be raised: indeed, 26 of 30 reported Saudi Arabian cases were detected in commercial farms containing thousands to hundreds of thousands of poultry . Perhaps, Arabian occurrence points reflect something other than the "ecological niche" of HPAI-H5N1 in the subregion; for example, they may reflect principally the conditions under which poultry can be raised, albeit with considerable subsidy of water and shade, irrespective of disease distributions. We should add, though we suspect that such is not the case, the total lack of predictivity in the Arabian Peninsula raises the more troubling possibility that the correspondence between NDVI and disease occurrence in the rest of the region is coincidental. It is possible that HPAI H5N1 distribution is not driven by factors correlated with NDVI seasonality, but by something that cannot be detected in the remotely sensed landscape.
Gilbert et al.  mapped the geographic distribution of suitable conditions for HPAI-H5N1 across Southeast Asia, finding close associations between free-grazing domestic ducks in rice paddies and HPAI-H5N1 cases. This result suggests that transmission risk could be mapped successfully in Southeast Asia, where duck production and rice cultivation are both extensive and intertwined, and that duck production may be an important driver of HPAI-H5N1 persistence. The authors stated that large numbers of Anatidae concentrate in the Nile Delta, and that the Hadejia Jama'are river system of Nigeria is also an important area of duck production. FAO reports a combined domestic duck and goose population of 18.3 million for Egypt in 2004 , presumably concentrated in the Nile Delta and Valley (along with virtually the entire human population and all productive agricultural land), joined in winter by large flocks (several hundreds of thousands ) of wild aquatic birds. Figures are unavailable for domestic Anatidae in Nigeria, although numbers of undifferentiated "exotic poultry" (ducks, geese, turkeys, guinea-fowl, ostriches, etc.) in the 5 states bordering Hadejia Jama'are were around 7.5 million birds in 2003 . Egypt and Nigeria both produce substantial rice crops (on 613 000 and 2 725 000 ha of land, respectively) .
Although total area under rice cultivation and total Anatidae populations are far higher in East Asia than in Egypt, the ratio of domestic Anatidae to area of rice production is considerably higher than in Thailand and Vietnam (see Additional file 3), and about the same as that found in China. If grazing of domestic Anatidae in rice paddies does play an important role in driving HPAI-H5N1 persistence and if duck-raising in the Middle East parallels that in East Asia, we might, expect persistence in China, Egypt, and Iran, all countries with higher duck-to-rice production area ratios than Thailand (Additional file 3). On the other hand, cases of HPAI-H5N1 have been numerous and widespread in Turkey, despite low numbers of Anatidae and little rice cultivation, suggesting that duck grazing in rice paddies is not the only factor in HPAI-H5N1 transmission and persistence. Free mingling of backyard poultry and wild birds has been identified as a risk factor for HPAI-H5N1 transmission [17, 18]. In Egypt, most domestic Anatidae are considered to be backyard (64% of ducks and "all" geese), whereas the majority of chickens (63%) are produced in commercial operations, apparently typifying the poultry industry of North Africa and the Middle East .
HPAI H5N1 detection data used for the development of these models are dominated by transmission among flocks of several poultry species. Given that detection data are so variable in terms of species composition (i.e., taxa, and number of taxa affected), husbandry method (high biosecurity, backyard, etc), origin (home-hatched, purchased, native-hatched, imported legally or illegally), and domestication, it is hard to define mechanisms driving transmission. We do not, however, find that our models are simply reproducing the spatial distributions of poultry flocks. Several ecologically-biased elements in the HPAI-H5N1 transmission cycle could explain the predictivity we detected: introduction of HPAI-H5N1 by migratory birds [20, 21], transmission among poultry flocks [22, 23], areas important for importation of poultry or hobby birds (legal or illegal) , or even transportation routes (e.g., roads, rivers). Inconsistencies in predictions based on HPAI-H5N1 occurrences from different subregions suggest that certain of these factors may have greater importance in some subregions than in others. In the Middle East, at least, we observe coincidence between human populations and HPAI-H5N1 cases, although, of course, this observation may simply point to the fact that influenza surveillance is more intensive in populated areas.
Ecological Niche Models
The ecological niche models (ENMs) developed in this study are based on the idea that organisms and other biological phenomena (including disease transmission) have distinct ecological requirements that determine their occurrences in time and space . In general, disease applications of ENM balance between focusing on individual species in the transmission system and using the integration of the whole system as a "black box" determining transmission to some species or biological phenomenon of interest [26, 27]. In this contribution, given the as-yet poorly characterized avian reservoir of HPAI-H5N1, we focus on all cases of HPAI-H5N1, effectively treating the transmission system as a black box. We thus attempt to model the transmission of a single pathogen based on its appearance in a multi-species system (i.e., the subset of animals in which HPAI H5N1 has been detected), in this case, dominated by distributions of domestic birds. In this sense, we deviate somewhat from the classical ENM approaches, which are based on single-species occurrence-environment correspondence. ENMs have been developed via diverse methodological approaches [28–31]; however, the method most frequently applied to questions of disease transmission has been the Genetic Algorithm for Rule-set Prediction (GARP), an evolutionary-computing approach [32, 33].
This study was based on HPAI-H5N1 animal case-occurrence data for 2005–2008 from the Middle East and northeastern Africa. Data were drawn from the World Organisation for Animal Health (OIE) , consisting of 610 unique locations, including isolations from wild birds, zoo birds, commercial poultry, and backyard poultry (Figure 1). This survey of occurrences includes birds assumed to be raised under strict biosecurity control, as well as birds raised with none; it similarly includes birds raised in strictly monospecific farms, multispecies assemblies mingling freely with wild birds (and other fauna), and even pets in a children's kindergarten. The database is composed of detections in at least 18 species of birds, although reporting standards are not consistent, so all too frequently information about hosts is either vague or absent. Most detections occurred in anthropogenic habitats. Our study area included 25 countries and one territory, ranging from Greece to the northwest, Somalia to the southwest, Georgia to the north, and Iran to the east. We assembled a complementary set of 17 unique and non-overlapping human cases occurrences from the archives of the International Society for Infectious Disease (ProMed Avian Influenza archive)  from the region (Figure 2) with which to test model predictions. All textual descriptions of occurrence localities were converted to geographic coordinates accurate to the nearest 0.01° using the GeoNet Names Server http://gnswww.nga.mil/geonames/GNS/index.jsp, Alexandria Digital Library Gazetteer http://middleware.alexandria.ucsb.edu/client/gaz/adl/index.jsp, and other sources .
We based ENM development on the 610 OIE localities for which geographic coordinates were provided with a precision of at least 0.01°; duplicate localities (i.e., multiple occurrences in the same 8 × 8 km grid square) were discarded. Geographic coordinates in the OIE data set were drawn from global positioning system recordings for the point of detection of HPAI-H5N1 cases . They thus specify the spatial position of HPAI-H5N1 occurrences, and probably represent the coarse-scale ecological conditions under which HPAI-H5N1 transmission occurs.
Given that the spatial pattern of H5N1 outbreaks has been on rather fine spatial scales, our previous experience with niche modeling and H5N1 outbreaks indicates that spatial resolutions on the order of 1–10 km are necessary, making use of climate-based data layers impractical. Environmental data sets included 12 monthly composite remotely-sensed data layers for Nov 1999 – Oct 2000, each summarizing maximum Normalized Difference Vegetation Index (NDVI; native spatial resolution 8 × 8 km) values ; although not exactly coincident with occurrence data temporally, these data provided an exemplar year of landscape variation in greenness. As NDVI is derived from reflectance in the visible and near-infrared domains, and as such is sensitive to photosynthetic activity and is closely correlated with photosynthetic mass , the NDVI time series used here summarizes aspects of land cover and vegetation phenology across the region. A year 2001 MODIS-based vegetation continuous fields dataset summarizing percent tree cover was also used (native spatial resolution 500 m) . Finally, we also included 3 data sets summarizing aspects of topography: slope, aspect, and compound topographic index (which summarizes tendency to pool water), from the U.S. Geological Survey's Hydro-1K data set (native resolution 1 km) . We deliberately excluded data on elevation from the study to avoid confusion caused by indirect variables. Climate data were not included in these analyses for lack of sufficiently high-resolution data sets across the region.
The GARP algorithm
The Genetic Algorithm for Rule-set Prediction (GARP) has been applied widely to questions of disease transmission [26, 39], and its predictive ability has been tested under diverse circumstances [30, 40, 41]. Although GARP has seen criticism in some comparative studies , more recent studies have indicated considerably better performance [42, 43] and some artifactual causation of previous results . As such, we used GARP for ENM development.
In general, we developed tests based on spatially stratified subsets of available occurrence information set aside prior to model development. Of occurrence data actually input into GARP, the program divides occurrence data randomly into three subsets: training data (25%; for rule development), intrinsic testing data (25%; for evaluation of rules) and extrinsic testing data (50%; for evaluation of model quality, see below). Spatial predictions of presence versus absence can include two types of error: false negatives (areas of actual presence predicted absent) and false positives (areas of actual absence predicted present)  – rule performance in each of these dimensions is evaluated via the intrinsic testing data set. Changes in predictive accuracy from one iteration to the next are used to evaluate whether particular rules should be incorporated into the model or not, and the algorithm runs either 1000 iterations or until convergence . The final rule-set is then used to query the environmental data sets across the study region to identify areas fitting the rule set predictions to produce a hypothesis of the potential geographic distribution of the species .
Since GARP processing includes several random-walk components, each replicate model produces distinct results, representing alternative solutions to the optimization challenge. Following best-practices approaches , we developed 100 replicates of each model. We filtered these replicates based on their error characteristics, retaining the 20 with lowest false negative rates (= percentage of independent testing points falling in areas not predicted to be suitable), and then retained the 10 (of the 20) closest to the median of proportional area predicted present, an index of false-positive error rates . A consensus of these 'best subset' models was then developed by summing values for each pixel in the map to produce final predictions of potential distributions with 11 thresholds (integers from 0 to 10).
Modeling and testing approach
This study focuses on the question of whether HPAI-H5N1 transmission in the Middle East and northeastern Africa occurs under a consistent and predictable set of environmental conditions. As such, we developed a series of tests of model predictivity; in each case, models were developed and predictions tested using spatially independent suites of occurrence data. Model tests were based on 4 spatial subsets of the Middle Eastern and northeastern African occurrence data (Figure 1): Arabian Peninsula (Bahrain, Kuwait, Oman, Qatar, Saudi Arabia, United Arab Emirates, Yemen; N = 31), Balkans-Caucasus (Armenia, Azerbaijan, Cyprus, Georgia, Greece, Turkey; N = 175), Levant-Iran (Iran, Iraq, Israel, Lebanon, Palestinian Territories, Syria; N = 18), and northeastern Africa (Djibouti, Egypt, Eritrea, Ethiopia, Somalia, Sudan; N = 386).
The basic design of testing included three schemes for subdividing available occurrence data, as follows:
1. Single testing regions: We combined each possible set of 3 subregional occurrence datasets to develop ENMs that were tested with the fourth subregion. Total 4 tests.
2. Single predictor regions: Occurrence data for each subregion were used to develop predictive models that were projected to the rest of the region for testing (e.g., Arabian Peninsula data points used to build predictions for the combination of Levant-Iran, northeastern Africa, and Balkans-Caucasus). Total 4 tests.
3. Predictivity across study region: We developed ENM predictions based on all OIE veterinary cases in the region, and tested its prediction based on coincidence of predictions with the 17 independent ProMed human cases. One test.
The customary approaches to spatial model validation (e.g. simple receiver operating characteristic, kappa statistics) are not applicable to situations in which presence-only data are the only information available [45, 46]. As such, we validated models using two approaches. First, we calculated binomial probabilities that observed coincidence of predictions and independent test data is no better than random, with the probability of k successes in N trials depending on p, the probability of success in any one trial; we estimated p as the proportion of the testing area predicted present, and k as the number of the N testing points successfully predicted by the model prediction . Binomial probabilities were calculated for each of the 10 thresholds representing predictions of presence (1 = broad, 10 = narrow), in each case testing whether predictivity was better than expected by chance.
Second, we followed Phillips et al.  in modifying receiver operating characteristic curves (ROCs) so as not to depend on absence data. We calculated the area under the curve (AUC) of the ROC, a statistical technique that has become a dominant tool in evaluating the accuracy of models predicting distributions of species 16. However, when comparing two ROCs, AUC systematically undervalues models that do not provide predictions across the entire spectrum of proportional areas in the study area (such as GARP, the modeling approach used here) . In addition current ROC approaches inappropriately weight the two error components (omission and commission) equally. Accordingly, we use a modification of ROC that remedies these problems: partial-area ROC approaches that evaluate only over the spectrum of the prediction and that allow for differential weighting of the two error components .
We carried out partial ROC analyses  for each model, all based solely on independent testing points not used to train the models in areas from distinct regions(s) to which models were projected. AUCs were limited to the proportional areas over which models actually made predictions, and only omission errors of <5% were considered (i.e., E = 5% ). We calculated partial AUCs using a program based on the trapezoid method  kindly developed by N. Barve (in prep.), and present our ROC results as the ratio of the model AUC to the null expectation ("AUC ratio") . Bootstrapping manipulations to permit evaluation of statistical significance of AUCs (as compared with null expectations) were achieved by resampling 50% of the test points with replacement 1000 times from the overall pool of testing data; one-tailed significance of differences in AUC (i.e. elevation above the line of null expectation) was assessed by counting the number of bootstrap replicates with AUC ratios <1.
Thanks to Monica Papeş and Yoshinori Nakazawa for expert assistance with GIS and satellite imagery. Thanks to Wildlife Conservation Society's Global Avian Influenza Network for Surveillance and the U.S. Geological Survey's National Biological Information Infrastructure for funding this work.
- Guan Y, Poon LLM, Cheung CY, Ellis TM, Lim W, Lipatov AS, Chan KH, Sturm-Ramirez KM, Cheung CL, Leung YHC: H5N1 influenza: a protean pandemic threat. Proceedings of the National Academy of Sciences, USA. 2004, 101: 8156-8161. 10.1073/pnas.0402443101.View ArticleGoogle Scholar
- Webby RJ, Webster RG: Are we ready for pandemic influenza?. Science. 2003, 302: 1519-1522. 10.1126/science.1090350.PubMedView ArticleGoogle Scholar
- de Jong JC, Claas ECJ, Osterhaus A, Webster RG, Lim WL: A pandemic warning?. Nature. 1997, 389: 554-10.1038/39218.PubMedView ArticleGoogle Scholar
- Whitworth D, Newman SH, Mundkur T, Harris P, Eds: Wild birds and avian influenza: an introduction to applied field research and disease sampling techniques. FAO Animal Production and Health Manual, Rome. 2007, 5
- Cumulative number of confirmed human cases of avian influenza A(H5N1) reported to WHO.http://www.who.int/csr/disease/avian_influenza/country/cases_table_2008_06_19/en/index.html
- H5N1 avian influenza: Timeline of major events.http://www.who.int/csr/disease/avian_influenza/Timeline_09_02_23.pdf
- Chen HX, Shen HG, Li XL, Zhou JY, Hou YQ, Guo JQ, Hu JQ: Seroprevalance and identification of influenza A virus infection from migratory wild waterfowl in China (2004–2005). Journal of Veterinary Medicine Series B-Infectious Diseases and Veterinary Public Health. 2006, 53: 166-170.View ArticleGoogle Scholar
- Liu J, Xiao H, Lei F, Zhu Q, Qin K, Zhang XW, Zhang XL, Zhao D, Wang G, Feng Y: Highly pathogenic H5N1 influenza virus infection in migratory birds. Science. 2005, 309: 1206-10.1126/science.1115273.PubMedView ArticleGoogle Scholar
- Peterson AT, Benz BW, Papeş M: Highly pathogenic H5N1 avian influenza: entry pathways into North America via bird migration. PLOS One. 2007, 2: e261-10.1371/journal.pone.0000261.PubMedPubMed CentralView ArticleGoogle Scholar
- Williams RAJ, Fasina FO, Peterson AT: Predictable ecology and geography of avian influenza (H5N1) transmission in Nigeria and West Africa. Transactions of the Royal Society of Tropical Medicine and Hygiene. 2008, 102: 471-479. 10.1016/j.trstmh.2008.01.016.PubMedView ArticleGoogle Scholar
- Update on highly pathogenic avian influenza in animals (Type H5 and H7).http://www.oie.int/downld/avian%20influenza/A_AI-Asia.htm
- Gilbert M, Xiao X, Pfeiffer DU, Epprecht M, Boles S, Czarnecki C, Chaitaweesub P, Kalpravidh W, Minh PQ, Otte MJ: Mapping H5N1 highly pathogenic avian influenza risk in Southeast Asia. Proceedings of the National Academy of Sciences, USA. 2008, 105: 4769-4774. 10.1073/pnas.0710581105.View ArticleGoogle Scholar
- GLiPHA, Global Livestock Production and Health Atlas.http://www.fao.org/ag/aga/glipha/index.jsp
- Hughes RH, Hughes JS: A directory of African wetlands. 1992, Gland, Switzerland and Cambridge, UK IUCNGoogle Scholar
- Adene DF, Oguntade AE: The structure and importance of the commercial and village based poultry industry in Nigeria. 2006, Rome, Italy: FAOGoogle Scholar
- Rice around the World.http://www.irri.org/science/ricestat/
- Chotpitayasunondh T, Ungchusak K, Hanshaoworakul W, Chunsuthiwat S, Sawanpanyalert P, Kijphati R, Lochindarat S, Srisan P, Suwan P, Osotthanakorn Y: Human disease from influenza A (H5N1), Thailand, 2004. Emerging Infectious Diseases. 2005, 11: 201-209.PubMedPubMed CentralView ArticleGoogle Scholar
- de Benedictis P, Joannis TM, Lombin LH, Shittu I, Beato MS, Rebonato V, Cattoli G, Capua I: Field and laboratory findings of the first incursion of the Asian H5N1 highly pathogenic avian influenza virus in Africa. Avian Pathology. 2007, 36: 115-117. 10.1080/03079450601161406.PubMedView ArticleGoogle Scholar
- Taha FA: The poultry sector in middle-income countries and its feed requirements: the case of Egypt. Agriculture and Trade Reports. 2003, US Department of Agriculture, Washington, DC, 47-Google Scholar
- Chen HL, Li YB, Li ZJ, Shi JZ, Shinya K, Deng GH, Qi QL, Tian GB, Fan SF, Zhao HD: Properties and dissemination of H5N1 viruses isolated during an influenza outbreak in migratory waterfowl in western China. Journal of Virology. 2006, 80: 5976-5983. 10.1128/JVI.00110-06.PubMedPubMed CentralView ArticleGoogle Scholar
- Olsen B, Munster VJ, Wallensten A, Waldenstrom J, Osterhaus A, Fouchier RAM: Global patterns of influenza A virus in wild birds. Science. 2006, 312: 384-388. 10.1126/science.1122438.PubMedView ArticleGoogle Scholar
- Gilbert M, Xiao XM, Domenech J, Lubroth J, Martin V, Slingenbergh J: Anatidae migration in the western palearctic and spread of highly pathogenic avian influenza H5N1 virus. Emerging Infectious Diseases. 2006, 12: 1650-1656.PubMedPubMed CentralView ArticleGoogle Scholar
- Yasue M, Feare CJ, Bennun L, Fiedler W: The epidemiology of H5N1 avian influenza in wild birds: why we need better ecological data. BioScience. 2006, 56: 923-929. 10.1641/0006-3568(2006)56[923:TEOHAI]2.0.CO;2.View ArticleGoogle Scholar
- Kilpatrick AM, Chmura AA, Gibbons DW, Fleischer RC, Marra PP, Daszak P: Predicting the global spread of H5N1 avian influenza. Proceedings of the National Academy of Sciences, USA. 2006, 103: 19368-19373. 10.1073/pnas.0609227103.View ArticleGoogle Scholar
- Soberón J, Peterson AT: Interpretation of models of fundamental ecological niches and species' distributional areas. Biodiversity Informatics. 2005, 2: 1-10.View ArticleGoogle Scholar
- Peterson AT: Ecological niche modeling and spatial patterns of disease transmission. Emerging Infectious Diseases. 2007, 12: 1822-1826.View ArticleGoogle Scholar
- McCormack JE, Peterson AT, Bonaccorso E, Smith TB: Speciation in the highlands of Mexico: genetic and phenotypic divergence in the Mexican jay (Aphelocoma ultramarina). Molecular Ecology. 2008, 17: 2505-2521. 10.1111/j.1365-294X.2008.03776.x.PubMedView ArticleGoogle Scholar
- Austin MP, Nicholls AO, Margules CR: Measurement of the realized qualitative niche: environmental niches of five Eucalyptus species. Ecological Monographs. 1990, 60: 161-177. 10.2307/1943043.View ArticleGoogle Scholar
- Carpenter G, Gillison AN, Winter J: DOMAIN: a flexible modeling procedure for mapping potential distributions of animals and plants. Biodiversity and Conservation. 1993, 2: 667-680. 10.1007/BF00051966.View ArticleGoogle Scholar
- Elith J, Graham CH, Anderson RP, Dudik M, Ferrier S, Guisan A, Hijmans RJ, Huettmann F, Leathwick JR, Lehmann A: Novel methods improve prediction of species' distributions from occurrence data. Ecography. 2006, 29: 129-151. 10.1111/j.2006.0906-7590.04596.x.View ArticleGoogle Scholar
- Pearson RG, Dawson TP, Berry PM, Harrison PA: SPECIES: a spatial evaluation of climate impact on the envelope of species. Ecological Modelling. 2002, 154: 289-300. 10.1016/S0304-3800(02)00056-X.View ArticleGoogle Scholar
- Stockwell DRB, Noble IR: Induction of sets of rules from animal distribution data: a robust and informative method of analysis. Mathematics and Computers in Simulation. 1992, 33: 385-390. 10.1016/0378-4754(92)90126-2.View ArticleGoogle Scholar
- Stockwell DRB, Peters DP: The GARP modelling system: problems and solutions to automated spatial prediction. International Journal of Geographical Information Science. 1999, 13: 143-158. 10.1080/136588199241391.View ArticleGoogle Scholar
- Avian Influenza Archive.http://www.promedmail.org
- Rand McNally: New Millenium World Atlas Deluxe. 1998, Skokie, Illinois, USA: Rand McNally & CoGoogle Scholar
- Tucker CJ: Red and photographic infrared linear combinations for monitoring vegetation. Remote Sensing of Environment. 1979, 8: 127-150. 10.1016/0034-4257(79)90013-0.View ArticleGoogle Scholar
- Hansen M, DeFries R, Townshend JR, Carroll M, Dimiceli C, Sohlberg R: MOD44B: Vegetation Continuous Fields Collection 3, Version 3.0.0. 2003, University of Maryland, College ParkGoogle Scholar
- HYDRO1k Elevation Derivative Database.http://edc.usgs.gov/products/elevation/gtopo30/hydro/index.html
- Costa J, Peterson AT, Beard CB: Ecologic niche modeling and differentiation of populations of Triatoma brasiliensis neiva, 1911, the most important Chagas' disease vector in northeastern Brazil (Hemiptera, Reduviidae, Triatominae). American Journal of Tropical Medicine and Hygiene. 2002, 67: 516-520.PubMedGoogle Scholar
- Anderson RP, Lew D, Peterson AT: Evaluating predictive models of species' distributions: criteria for selecting optimal models. Ecological Modelling. 2003, 162: 211-232. 10.1016/S0304-3800(02)00349-6.View ArticleGoogle Scholar
- Peterson AT, Papeş M, Kluza DA: Predicting the potential invasive distributions of four alien plant species in North America. Weed Science. 2003, 51: 863-868. 10.1614/P2002-081.View ArticleGoogle Scholar
- McNyset K, Blackburn J: Does GARP really fail miserably? A response to Stockman et al. (2006). Diversity and Distributions. 2006, 12: 782-786. 10.1111/j.1472-4642.2006.00281.x.View ArticleGoogle Scholar
- Tsoar A, Allouche O, Steinitz O, Rotem D, Kadmon R: A comparative evaluation of presence-only methods for modelling species distribution. Diversity and Distributions. 2007, 13: 397-405. 10.1111/j.1472-4642.2007.00346.x.View ArticleGoogle Scholar
- Peterson AT, Papeş M, Eaton M: Transferability and model evaluation in ecological niche modeling: a comparison of GARP and Maxent. Ecography. 2007, 30: 550-560.View ArticleGoogle Scholar
- Fielding AH, Bell JF: A review of methods for the assessment of prediction errors in conservation presence/absence models. Environmental Conservation. 1997, 24: 38-49. 10.1017/S0376892997000088.View ArticleGoogle Scholar
- Manel S, Dias JM, Buckton ST, Ormerod SJ: Alternative methods for predicting species distribution: an illustration with Himalayan river birds. Journal of Applied Ecology. 1999, 36: 734-747. 10.1046/j.1365-2664.1999.00440.x.View ArticleGoogle Scholar
- Phillips SJ, Anderson RP, Schapire RE: Maximum entropy modeling of species geographic distributions. Ecological Modelling. 2006, 190: 231-259. 10.1016/j.ecolmodel.2005.03.026.View ArticleGoogle Scholar
- Peterson AT, Papeş M, Soberon J: Rethinking receiver operating characteristic analysis applications in ecological niche modeling. Ecological Modelling. 2008, 213: 63-72. 10.1016/j.ecolmodel.2007.11.008.View ArticleGoogle Scholar
- Burden RL, Faires JD: Numerical Analysis. 2005, Belmont, California: Thomson Books, 8Google Scholar
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/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.