Determinants of the geographic distribution of Puumala virus and Lyme borreliosis infections in Belgium
© Linard et al; licensee BioMed Central Ltd. 2007
Received: 09 March 2007
Accepted: 02 May 2007
Published: 02 May 2007
Vector-borne and zoonotic diseases generally display clear spatial patterns due to different space-dependent factors. Land cover and land use influence disease transmission by controlling both the spatial distribution of vectors or hosts, and the probability of contact with susceptible human populations. The objective of this study was to combine environmental and socio-economic factors to explain the spatial distribution of two emerging human diseases in Belgium, Puumala virus (PUUV) and Lyme borreliosis. Municipalities were taken as units of analysis.
Negative binomial regressions including a correction for spatial endogeneity show that the spatial distribution of PUUV and Lyme borreliosis infections are associated with a combination of factors linked to the vector and host populations, to human behaviours, and to landscape attributes. Both diseases are associated with the presence of forests, which are the preferred habitat for vector or host populations. The PUUV infection risk is higher in remote forest areas, where the level of urbanisation is low, and among low-income populations. The Lyme borreliosis transmission risk is higher in mixed landscapes with forests and spatially dispersed houses, mostly in wealthy peri-urban areas. The spatial dependence resulting from a combination of endogenous and exogenous processes could be accounted for in the model on PUUV but not for Lyme borreliosis.
A large part of the spatial variation in disease risk can be explained by environmental and socio-economic factors. The two diseases not only are most prevalent in different regions but also affect different groups of people. Combining these two criteria may increase the efficiency of information campaigns through appropriate targeting.
Numerous human infectious diseases, particularly vector-borne and zoonotic diseases, are thought to emerge or re-emerge in some parts of the world due to environmental changes. The spatial distribution of two emerging human diseases were investigated and compared in Belgium: Puumala virus (PUUV) and Lyme borreliosis. Since the former is a zoonotic rodent-borne infection (i.e. transmissible from vertebrate animals to humans) and the latter a tick-borne infection, zoonotic and vector-borne at the same time, both strongly depend on the natural environment.
In Europe, hantaviruses (family Bunyaviridae) are causative agents of human zoonoses called haemorrhagic fever with renal syndrome (HFRS). Hantaviruses are transmitted via rodent excretions, unlike the other viruses from the Bunyaviridae family that are transmitted by arthropod vectors . In western and northern Europe, Puumala virus (PUUV) causes a mild form of HFRS: nephropathia epidemica. PUUV is associated predominantly with the bank vole (Clethrionomys glareolus) [2–5]. Ecological factors such as favourable habitat, food supply, and climatic conditions cause interannual fluctuations in populations of bank voles.
Lyme borreliosis is an infection caused by a spirochete named Borrelia burgdorferi. This zoonosis is a tick-borne disease present in the northern hemisphere where the main vector species in European countries is Ixodes ricinus [6, 7]. The complexity of disease dynamics is due to the life cycle of ticks, composed of three developmental stages (larvae, nymph, and adult). Infected ticks can transmit spirochetes to the host every time they take a blood meal, which occurs once per stage. Each of the components of the system (spirochete, tick vector, and hosts) is influenced by external biotic and abiotic factors . Consequently, human activity can directly or indirectly influence the disease ecology by its impact on landscapes and host populations.
Vector-borne and zoonotic diseases generally display clear spatial patterns that can be due to different space-dependent factors: (i) the limited pathogen dispersal, (ii) the spatial distribution of vectors, hosts and reservoirs, and (iii) the human exposure to the infectious agent. Ostfeld et al. (2005) define spatial epidemiology as "the study of spatial variation in disease risk or incidence" . Land cover – defined by the attributes of the earth's land surface and immediate subsurface, including biota, soil, topography, surface and groundwater, and human structures – and land use – defined by the purposes for which humans exploit the land cover – influence directly the spatial variation in disease risk by controlling the different agents of the disease transmission process (pathogen, vectors and hosts, including humans) and the degree of contact between these agents .
Geographical information systems (GIS) and remote sensing are appropriate tools to describe the spatial distribution of infectious diseases and predict disease risk [16–20]. These tools have already been used to explain or predict tick-borne diseases [21–24] and rodent-borne diseases [25, 26]. These previous studies were aimed at relating spatial data on land cover and climate to the ecology of vectors or hosts.
The objective of this study was to combine environmental and socio-economic factors to explain the geographic distribution of disease incidence of Puumala virus and Lyme borreliosis in Belgium. We used multivariate statistical analysis to explain the number of PUUV and Lyme borreliosis infections per municipality according to variables related to land cover, land use, settlement characteristics, socio-economic variables, risk behaviours, and host population. We used the disease incidence as proxy for disease risk because: (i) data on human infections are better inventoried than data on hosts (in the case of PUUV) or vectors (in the case of Lyme borreliosis), and (ii) human infection results from the presence of vectors/hosts and the pathogen but also from people's exposure to vectors/hosts . Human exposure is generally not taken into account in studies of vector or host ecology (but see ).
Data and methods
The surveillance of Lyme borreliosis and hantavirus infections is carried out through the Scientific Institute of Public Health (IPH) sentinel laboratory network. As the diagnosis of these infections is complex and needs to be confirmed by a specialised laboratory, sentinel laboratories are asked to send positive samples to the Reference Laboratories for confirmation. These Reference Laboratories (e.g., the Reference Laboratory for Vector-borne diseases in Brussels) send monthly data on confirmed cases to the IPH and duplicated records are eliminated. The IPH collects data per postal code, which is a subdivision of the Belgian municipality. As socio-economic census data are only publicly available at the level of municipalities, we simply aggregated data from postal codes to municipalities and conducted the study at that level.
The number of human infections was used as dependent variable and the population number was added as offset variable in the model. Therefore, the model explains the incidence rate per municipality, as it was done in , under the assumption that incidence and disease transmission risk are highly correlated . For Puumala virus, we considered all the human infections that were diagnosed in Belgium between 1994 and 2004. For Lyme borreliosis infections, we only consider the period 1998–2004 given data availability.
The data collection method was slightly different for the two diseases. While PUUV patients were georeferenced according to their residential municipality in the IPH database, Lyme borreliosis patients were georeferenced according to the municipality where they thought the contamination took place. Since, in most cases, Lyme contamination by a tick took place close to one's house, we assumed that the municipality of contamination that was recorded is the municipality of residence. We thus considered municipalities of residence for both diseases. This assumption allows combining socio-economic and environmental characteristics of municipalities in explaining disease occurrence.
Summary statistics of dependent and independent variables
total number of PUUV infections during the 1994–2004 period
total number of Lyme borreliosis infections during the 1998–2004 period
Land cover and land use
proportion of forest area (%)
proportion of broad-leaved forest area (%)
urbanisation level from a morphological point of view (5 categories)
proportion of people living in a separated house (%)
proportion of people living in a twinned house (%)
proportion of people living in an adjoining house (%)
proportion of people living in an apartment (%)
proportion of people living in a caravan or country cottage (%)
average income per 1000 inhabitants in 2002
average population for the 1994–2004 period
average population for the 1998–2004 period
proportion of people with a hunting licence for Flemish or Walloon forests during the hunting year 2004–2005 (%)
density of roe deer (heads/km2)
Land cover and land use
Land cover determines the spatial distribution of vectors, hosts and animal reservoirs according to their habitat preferences. Broad-leaved forests, being considered as the favourite habitat for bank voles and ticks alike are expected to have an impact on disease risk. The CORINE Land Cover 2000 database, with a spatial resolution of 100 meters, was produced between 1999 and 2001 on the basis of satellite image interpretation . The proportion of municipality area covered by forests (categories 3.1.1, 3.1.2 and 3.1.3 from CORINE Land Cover 2000) and broad-leaved forests (category 3.1.1 from CORINE Land Cover 2000) were extracted.
Land use determines human presence in a given place (settlement, leisure, agriculture,...) at particular times of the year and of the day. In built-up areas, the presence of vectors and hosts is very constrained ecologically and disease risk is expected to be low. A higher disease risk is expected in landscape mosaics of human-dominated and natural covers, where contacts between humans and vectors or hosts are more frequent . To test the influence of the level of urbanisation, we used a classification based on the 1991 survey of the Directorate-general Statistics Belgium that classifies municipalities in 5 categories according to their morphological urbanisation level: central municipalities of main towns; high, medium and low level of morphological urbanisation; and rural municipalities . The morphological urbanisation level of municipalities was determined using the population density and the proportion of built-up areas.
Settlement density and types are expected to have an impact on the likelihood of contacts between humans and vectors or hosts. People living at the vicinity of forests are more likely to be infected because they are more often in direct contact with soil and wood potentially infected by PUUV. For some PUUV human cases in France, the only risk factor found was the proximity between habitation and forest . Outdoor activities, even in the garden, also increase the probability of being bitten by a tick. The 2001 socio-economic survey of the Directorate-general Statistics Belgium provides data on population distribution according to residential characteristics . Habitations are divided in 5 categories: separated houses, twinned houses, adjoining houses, apartments, and mobile homes or country cottages. We used the proportion of people living in these different settlement types for each municipality.
Certain professional or leisure activities that are dominant in certain municipalities could be reflected in the socio-professional level. For example, municipalities with many farmers and forest workers, such as loggers or foresters, which are more exposed to disease risk, generally correspond to lower mean incomes. By contrast, peri-urban zones, characterized by large gardens surrounded by woods generally correspond to higher incomes. A higher income may also enhance the propensity to consult a doctor and therefore the probability to detect the disease. Since 2000, the Directorate-general Statistics Belgium performs an annual survey on Belgian household budgets. We used data on average taxable income per tax return form in each municipality from 2002, the most recent publicly available data . Data on the number of inhabitants per municipality are also available for 1993 and 2000 to 2004. Data were interpolated for missing years to estimate the average number of inhabitants for the period 1994–2004 for PUUV and 1998–2004 for Lyme borreliosis.
Specific human behaviours are expected to increase risk for the two diseases. Risk activities for PUUV include working in forest, farming, animal trapping, woodcutting or reopening rodent infested buildings [2, 10, 32], and are mainly linked to professional activities . Lyme borreliosis incidence is more linked to peridomestic and recreational activities such as gardening, hiking or other outdoor activities. Few data on these behaviours are available. A detailed analysis of risk behaviours would require data at the individual level rather than municipalities, as it was done in [2, 5]. In this study, we included only the proportion of people who were holding a hunting licence for Flemish or Walloon forests during the hunting year 2004–2005. These data were provided at the municipality level by institutions that deliver licences to hunt either in Walloon forests (data provided by the "Nature and forest" Division of the Walloon Region, Hunting and Fishing direction, 2005) or in Flemish forests (data provided by district police stations, Nature and Forest agency of the Flanders Region, 2006) and were associated with the municipality of residence. We assumed that people have only one licence, to hunt either in Walloon or Flemish forests.
The roe deer, the main host of Ixodes ricinus – although not a reservoir – plays an important role of maintenance and co-feeding for ticks [35, 36]. Roe deer abundance should therefore be related to tick abundance and Lyme borreliosis transmission risk, as it was verified in Denmark . The absence of deer could also amplify tick abundance at small spatial scales . Censuses on the number of roe deer (estimation of the number of roe deer on feet) are compiled by regional institutions in Belgium. Data were provided by: (i) the "Instituut voor Bosbouw and Wildbehher" (IBW), "Wildbeheer" department, "afschot reewild" for Flanders, (ii) the "grand gibier" statistics from the "Nature and forest" Division of the Walloon Region, and (iii) the Brussels Environment (IBGE-BIM) for Brussels . Data on roe deer abundance are available at the province level in Flanders and at the forest district level in Wallonia and Brussels. These data were spatially disaggregated at the level of municipalities to be compatible with the scale of the dependent variable. For Wallonia, we first overlaid forest district and municipality boundaries in a GIS. Deer populations were allocated to each polygon resulting from the intersection between forest districts and municipalities according to the proportion of forested area of the forest district belonging to the municipality. The number of deer per municipality was then obtained by aggregation of the different polygons in a municipality. For Flanders, as municipalities are fully included in provinces, we simply allocated a deer population to municipalities in proportion to the forest area of each municipality. For Brussels, we only had an estimate of a density of 3 roe deer for 100 hectares of forest . Using this estimate and the proportion of forest area in each municipality of Brussels, we computed the expected number of deer per municipality. For the three regions, we considered the average number of deer over the 1998–2004 period. This variable was used only in the model for Lyme borreliosis.
Statistical distributions analysis
When k = 0, the negative binomial distribution is reduced to a Poisson distribution.
Likelihood ratio test and Dean's tests for overdispersion
Likelihood Ratio Test
Dean's P_B test
Dean's P'_B test
Generalized linear regression models were computed using the R statistical software. All variables were relative to the population size Ni of municipalities and the logarithm of Ni was added as an offset variable in the models. The distribution parameter k is estimated by the maximum likelihood method, as the other parameters of the negative binomial regression model. The Akaike's information criterion (AIC), which takes the number of estimated parameters into account, was used to select the covariates to retain in the model. Variables were tested for collinearity and confounding factors.
Spatial autocorrelation appears when a spatial dependence due to neighbouring effects (i.e., spatial externalities) exists between spatial entities. Spatial dependence usually results from a combination of exogenous and endogenous processes. In an exogenous spatial process, the identified spatial pattern is generated by independent factors, whereas endogenous spatial dependence is an inherent property of the variable of interest . Spatial epidemiological models are particularly prone to spatial autocorrelation since disease transmission requires contacts between infected and susceptible agents. The infection rate in a locality is therefore influenced by the infection rates in neighbouring spatial entities. The spatial diffusion of epidemic diseases and the influence on infection of many unobserved factors also cause spatial autocorrelation.
where Yi is the expected value of the dependent variable y for the municipality i, xin are the dependent variables with their associated regression coefficients β n , σε i is the disturbance term and is the linear prediction of the dependent variable in neighbouring spatial entities j. This model was implemented in two stages, following the method of Mallar (1977) . Firstly, a separate negative binomial regression model, with the number of human infections in the surrounding municipalities as dependent variable, provided linear predictions of Y Ej in each municipality. Covariates included the same variables as for the non-spatial model computed for both the central and its neighbouring municipalities. Interaction variables were also included to maximize the performance of the model. In the second stage, these linear predictions were included among the explanatory variables in the negative binomial regression.
Parameter estimates of significant variables using negative binomial regressions on PUUV infections (1994–2004)
Linear predictors in neigbourhood municipalities
Degrees of freedom
2 × log-likelihood
Parameter estimates of significant variables using negative binomial regressions on Lyme borreliosis infections (1998–2004)
Linear predictors in neigbourhood municipalities
Degrees of freedom
2 × log-likelihood
The variograms of residuals for these non-spatial regressions lead to a different diagnostic for the two diseases (Figure 4). For PUUV, the envelope contains all the points of the empirical variogram. This suggests that the spatial dependence in PUUV incidence is mainly due to an exogenous process that is associated with one or several explanatory variables included in the model, which have their own spatial structure. The proportion of broad-leaved forests, one of the significant independent variables, explains a large part of this spatial dependency in PUUV incidence. By contrast, for Lyme borreliosis, the residuals of the regression model display a significant spatial dependence. This could be caused by unobserved explanatory variables, i.e. an external process whose spatial structure is not represented in the model. It could also be due to an endogenous effect if the incidence of Lyme borreliosis in one location influences its incidence in neighbouring locations.
To test for spatial endogeneity in the Lyme borreliosis model, a linear predictor of the dependent variable in neighbouring municipalities was estimated. For comparative purpose, the same approach was applied to the PUUV model, even though the spatial dependency was already reduced in the non-spatial model for this disease. Both spatial models have a high explanatory power and very significant parameters (Tables 3 and 4). For PUUV, the residual deviance value is closer to DF compared to the non-spatial model, with a ratio that reaches 0.83. For Lyme, the gap between the residual deviance and DF is similar to the non-spatial model (with a ratio of 1.1). Each model includes the same explanatory variables as the non-spatial models. All these covariates – except for roe deer density for Lyme – remain significant and with the same signs. The linear predictors of the dependent variable in neighbouring municipalities are very significant in both models.
The spatial model on PUUV is composed of four variables, each with the expected sign, in addition to the spatial dependency term that is also significant. Two variables are related to the landscape (proportion of the area of the municipality occupied by broad-leaved forests and urbanisation), one to a socio-economic characteristic of the population (income), and another one to a risk behaviour of humans (hunting). The PUUV incidence rate is higher where the proportion of broad-leaved forests is high and where urbanisation is low. Broad-leaved forests are the favourite habitat for bank voles and are thus positively related to the disease. Moreover, bank voles use small areas and only rarely leave forests. Therefore, humans generally catch the disease in the ecological habitat of bank voles, in broad-leaved forests. A high level of urbanisation thus limits PUUV transmission. Income is negatively correlated with the disease incidence because PUUV particularly affects forest workers and farmers who are in frequent contact with bank voles or infected soil in forests or at their edge. These socio-professional categories are generally characterised by a low income. Hunters also spend time in forests and are thus particularly vulnerable to PUUV infection.
The spatial model for Lyme borreliosis infections has four explanatory variables, each of them with a positive sign, in addition to the spatial dependency term that is also highly significant. One variable is related to the host population (roe deer), two variables represent landscape attributes that favour Lyme borreliosis transmission (proportion of the area of the municipality occupied by forest and proportion of people living in separated houses), and one characterises the human population by its socio-professional level. The model confirms the hypothesis that roe deer density has an impact on spatial variations in Lyme disease transmission risk in Belgium, probably by favouring tick reproduction and co-feeding. Lyme borreliosis incidence rate is positively linked to forest cover but also to the proportion of people living in separated houses. This suggests that heterogeneous landscapes with a fragmented land use mixing forests and houses are more at risk. In areas with a high proportion of separated houses, people are more likely to have gardens and thus spend time outdoors. The presence of forests nearby further favours human-vector contacts. Peri-urbanisation could thus be one of the major causes of the recent increase in Lyme borreliosis infections. In Belgium, it largely affects Walloon Brabant and the province of Antwerp, two wealthy peri-urban areas. The positive relationship between Lyme borreliosis incidence rate and mean income supports this hypothesis since higher incomes are likely to be associated with large gardens in wooded areas and thus a larger interface between settlements and forests in peri-urban areas. Higher incomes – especially when they are associated with education level – could also increase the probability of being aware of the disease and of consulting a doctor, thus increasing the likelihood of disease detection.
The spatial distribution of PUUV and Lyme borreliosis infections are thus associated with similar factors, in particular the presence of forests. For both diseases, it is linked to the vector or host habitat preferences. Other variables influence both diseases but in opposite directions, such as mean income and urbanisation. This reflects the fact that PUUV incidence is more associated with professional outdoor activities (in low-income, rural and forested areas), whereas Lyme borreliosis is rather associated with recreational and peridomestic outdoor activities (in high income, peri-urban areas with isolated houses and forests).
Most importantly, both diseases are associated with a combination of factors linked to the vector and host populations, to human behaviours, and to landscape attributes. The latter includes both land cover (forests), as a function of vector/host ecology, and land use that determines human presence in critical habitats, the frequency of contacts between vectors and human populations, and thus human exposure to infection.
Residuals of the spatial regression model of Lyme disease still display strong spatial dependence, which suggests that one or more external factors with spatial structure are missing. However, several factors could influence and skew the spatial effect in the Lyme model. Firstly, the roe deer density data were only available at the province level in Flanders and at the forest district level in the Walloon region and Brussels. This variable was thus spatially disaggregated at the level of municipalities. Removing this variable did not reduce the spatial autocorrelation. However, representing more ecologically-relevant boundaries that separate independent deer populations (e.g., major roads) could reduce spatial autocorrelation compared to administrative boundaries . Secondly, data on Lyme borreliosis incidence are based on the contamination place, whereas data on PUUV infections are based on the residential municipality. The assumption that the municipality of contamination corresponds to that of residence for Lyme disease could create a bias, especially for fine-scale analyses. This bias could increase spatial dependence effects or simply add randomly distributed errors . Thirdly, information campaigns organised for doctors in certain regions could have led to a better diagnostic of Lyme borreliosis infections. High positive residuals are found for example around the municipality of Ottignies-Louvain-la-Neuve in Walloon Brabant, where an important information centre on Lyme disease – Research Group and Information on tick-borne diseases (RILY) – is located. This is highly speculative however since this municipality is also a high-income, peri-urban area with several small forests.
The objective of this study was to explain the spatial distribution of a zoonotic (Puumala virus) and a vector-borne (Lyme borreliosis) disease in Belgium. The combination of environmental and socio-economic factors allowed explaining the major part of the spatial variation in disease risk. Besides the influence of land cover on the spatial distribution of vectors and hosts, land use also determines the prevalence of vector-borne and zoonotic diseases by controlling the degree of exposure of people to contact with vectors or hosts. PUUV infection risk is more prevalent in remote, low-income areas with large forests whereas the probability of Lyme borreliosis infection is higher in areas with a large interface between settlements and forests in wealthy peri-urban areas. This study showed the importance of adding variables related to human behaviour and land use in addition to land cover and habitat variables in epidemiological models. Spatial dependence diagnostics in error terms revealed a combination of endogenous and exogenous processes causing spatial clustering of incidence for both diseases. The spatial dependence could be accounted for in the model of PUUV incidence but not for Lyme borreliosis.
The results of this study have important implications to define priority areas for public health policies aimed at preventing diseases by informing the public and promoting the use of protective measures. For example, the Institute for Public Health produces information sheets on the diseases for distribution to the public. Each disease is not only most prevalent in different regions but also affects different groups of people. Combining these two criteria may increase the efficiency of information campaigns through appropriate targeting.
This research was partially funded by EU grant GOCE-2003-010284 EDEN and the paper is catalogued by the EDEN Steering Committee as EDEN0043 http://www.eden-fp6project.net/. The contents of this publication are the responsibility of the authors and do not necessarily reflect the views of the European Commission. We are grateful to Belgian institutions that provided data on hunting licences and roe deer censuses: the "Nature and forest" Division of the Walloon Region, the "Nature and forest" Agency of the Flanders Region and the "Instituut voor Bosbouw en Wildbeheer" (IBW).
- Krüger DH, Ulrich R, Lundkvist A: Hantavirus infections and their prevention. Microbes Infect. 2001, 3: 1129-1144.PubMedView ArticleGoogle Scholar
- Van Loock F, Thomas I, Clement J, Ghoos S, Colson P: A Case-Control Study after a Hantavirus Infection Outbreak in the South of Belgium: Who Is at Risk?. Clin Infect Dis. 1999, 28 (4): 834-839.PubMedView ArticleGoogle Scholar
- Escutenaire S, Chalon P, Verhagen R, Heyman P, Thomas I, Karelle-Bui L, Avsic-Zupanc T, Lundkvist A, Plyusnin A, Pastoret PP: Spatial and temporal dynamics of Puumala hantavirus infection in red bank vole (Clethrionomys glareolus) populations in Belgium. Virus Res. 2000, 67: 91-107.PubMedView ArticleGoogle Scholar
- Heyman P, Van Mele R, De Jaegere F, Klingström J, Vandenvelde C, Lundkvist A, Rozenfeld F, Zizi M: Distribution of hantavirus foci in Belgium. Acta Trop. 2002, 84: 183-188.PubMedView ArticleGoogle Scholar
- Olsson GE, Dalerum F, Hörnfeldt B, Elgh F, Palo TR, Juto P, Ahlm C: Human Hantavirus Infections, Sweden. Emerg Infect Dis. 2003, 9 (11): 1395-1401.PubMedPubMed CentralView ArticleGoogle Scholar
- Gern L, Humair PF: Ecology of Borrelia burgdorferi sensu lato in Europe. Lyme borreliosis: biology, epidemiology and control. Edited by: Gray JS, Kahl O, Lane RS, Stanek G. 2002, New York: CABI Publishing, 149-174.View ArticleGoogle Scholar
- Hengge UR, Tannapfel A, Tyring SK, Erbel R, Arendt G, Ruzicka T: Lyme borreliosis. Lancet Infect Dis. 2003, 3 (8): 489-500.PubMedView ArticleGoogle Scholar
- Kahl O, Gern L, Eisen L, Lane RS: Ecological Research on Borrelia burgdorferi sensu lato: Terminology and Some Methodological Pitfalls. Lyme borreliosis: biology, epidemiology and control. Edited by: Gray JS, Kahl O, Lane RS, Stanek G. 2002, New York: CABI Publishing, 29-46.View ArticleGoogle Scholar
- Heyman P, Vervoort T, Escutenaire S, Degrave E, Konings J, Vandenvelde C, Verhagen R: Incidence of hantavirus infections in Belgium. Virus Res. 2001, 77: 71-80.PubMedView ArticleGoogle Scholar
- Penalba C, Galempoix JM, Lanoux P: Epidémiologie des infections à hantavirus en France. Med Mal Infect. 2001, 31 (Suppl 2): 272-284.View ArticleGoogle Scholar
- Vapalahti O, Mustonen J, Lundkvist A, Henttonen H, Plyusnin A, Vaheri A: Hantavirus infections in Europe. Lancet Infect Dis. 2003, 3: 653-661.PubMedView ArticleGoogle Scholar
- Randolph SE: The shifting landscape of tick-borne zoonoses: tick-borne encephalitis and Lyme borreliosis in Europe. Philos Trans R Soc Lond B Biol Sci. 2001, 356: 1045-1056.PubMedPubMed CentralView ArticleGoogle Scholar
- Mailles A, Sin MA, Ducoffre G, Heyman P, Koch J, Zeller H: Larger than usual increase in cases of hantavirus infections in Belgium, France and Germany, June 2005. Eurosurveillance. 2005, 10 (7):
- Ostfeld RS, Glass GE, Keesing F: Spatial epidemiology: an emerging (or re-emerging) discipline. Trends Ecol Evol. 2005, 20 (6): 328-336.PubMedView ArticleGoogle Scholar
- Ward SE, Brown RD: A framework for incorporating the prevention of Lyme disease transmission into the landscape planning and design process. Landsc Urban Plann. 2004, 66 (2): 91-106.View ArticleGoogle Scholar
- Kitron U: Landscape Ecology and Epidemiology of Vector-Borne Diseases: Tools for Spatial Analysis. J Med Entomol. 1998, 35 (4): 435-445.PubMedView ArticleGoogle Scholar
- Brooker S, Michael E: The Potential of Geographical Information Systems and Remote Sensing in the Epidemiology and Control of Human Helminth Infections. Remote Sensing and Geographical Information Systems in Epidemiology. Edited by: Hay SI, Randolph SE, Rogers DJ. 2000, Oxford: Academic Press, 245-270. [Baker J, Muller R, Rollinson D (Series Editors): Advances in Parasitology, vol 47.]View ArticleGoogle Scholar
- Curran PJ, Atkinson PM, Foody GM, Milton EJ: Linking Remote Sensing, Land Cover and Disease. Remote Sensing and Geographical Information Systems in Epidemiology. Edited by: Hay SI, Randolph SE, Rogers DJ. 2000, Oxford: Academic Press, 37-78. [Baker J, Muller R, Rollinson D (Series Editors): Advances in Parasitology, vol 47.]View ArticleGoogle Scholar
- Hay SI, Omumbo JA, Craig MH, Snow RW: Earth Observation, Geographic Information Systems and Plasmodium falciparum Malaria in Sub-Saharan Africa. Remote Sensing and Geographical Information Systems in Epidemiology. Edited by: Hay SI, Randolph SE, Rogers DJ. 2000, Oxford: Academic Press, 174-215. [Baker J, Muller R, Rollinson D (Series Editors): Advances in Parasitology, vol 47.]Google Scholar
- Rogers DJ: Satellites, Space, Time and the African Trypanosomiases. Remote Sensing and Geographical Information Systems in Epidemiology. Edited by: Hay SI, Randolph SE, Rogers DJ. 2000, Oxford: Academic Press, 129-171. [Baker J, Muller R, Rollinson D (Series Editors): Advances in Parasitology, vol 47.]View ArticleGoogle Scholar
- Dister SW, Fish D, Bros SM, Frank DH, Wood BL: Landscape characterization of peridomestic risk for lyme disease using satellite imagery. Am J Trop Med Hyg. 1997, 57 (6): 687-692.PubMedGoogle Scholar
- Daniel M, Kolàr J, Zeman P, Pavelka K, Sádlo J: Predictive map of Ixodes ricinus high-incidence habitats and a tick-borne encephalitis risk assessment using satellite data. Exp Appl Acarol. 1998, 22: 417-433.PubMedView ArticleGoogle Scholar
- Randolph SE: Ticks and Tick-borne Disease Systems in Space and from Space. Remote Sensing and Geographical Information Systems in Epidemiology. Edited by: Hay SI, Randolph SE, Rogers DJ. 2000, Oxford: Academic Press, 217-243. [Baker J, Muller R, Rollinson D (Series Editors): Advances in Parasitology, vol 47.]View ArticleGoogle Scholar
- Ogden NH, Barker IK, Beauchamp G, Brazeau S, Charron DF, Maarouf A, Morshed MG, O'Callaghan CJ, Thompson RA, Waltner-Toews D, Waltner-Toews M, Lindsay LR: Investigation of ground level and remote-sensed data for habitat classification and prediction of survival of Ixodes scapularis in habitats of southeastern Canada. J Med Entomol. 2006, 43 (2): 403-414.PubMedView ArticleGoogle Scholar
- Glass GE, Cheek JE, Patz JA, Shields TM, Doyle TJ, Thoroughman DA, Hunt DK, Enscore RE, Gage KL, Irland C, Peters CJ, Bryan R: Using Remotely Sensed Data To Identify Areas at Risk For Hantavirus Pulmonary Syndrome. Emerg Infect Dis. 2000, 6 (3): 238-247.PubMedPubMed CentralView ArticleGoogle Scholar
- Goodin DG, Koch DE, Owen RD, Chu Y, Hutchinson JS, Jonsson CB: Land cover associated with hantavirus presence in Paraguay. Glob Ecol Biogeogr. 2006, 15: 519-527.View ArticleGoogle Scholar
- Horobik V, Keesing F, Ostfeld RS: Abundance and Borrelia burgdorferi-infection prevalence of nymphal Ixodes scapularis ticks along forest-field edges. EcoHealth. 2006, 3 (4): 262-268.View ArticleGoogle Scholar
- Jackson LE, Hilborn ED, Thomas JC: Towards landscape design guidelines for reducing Lyme disease risk. Int J Epidemiol. 2006, 35 (2): 315-322.PubMedView ArticleGoogle Scholar
- Rácz GR, Bán E, Ferenczi E, Berencsi G: A Simple Spatial Model to Explain the Distribution of Human Tick-Borne Encephalitis Cases in Hungary. Vector Borne Zoonotic Dis. 2006, 6 (4): 369-378.PubMedView ArticleGoogle Scholar
- The CORINE Land Cover 2000 database. [http://terrestrial.eionet.europa.eu/CLC2000]
- Halleux JM, Derwael F, Mérenne-Schoumaker B: Typologie des communes selon leur degré d'urbanisation. Monographie n°11A du Recensement Général de la Population et des Logements au 1er mars 1991. Edited by: INS, SSTC. 1998, Bruxelles, 149-182.Google Scholar
- Mailles A, Vaillant V, Haeghebaert S, Fradet MR, Capek I, Zeller H: Recrudescence des infections à Hantavirus, France 2003. Med Mal Infect. 2005, 35 (2): 68-72.PubMedView ArticleGoogle Scholar
- Directorate-general Statistics Belgium: Enquête socio-économique 2001. 2001, [http://statbel.fgov.be/census/]Google Scholar
- Directorate-general Statistics Belgium: Enquête sur le budget des ménages. [http://statbel.fgov.be]
- Jaenson TGT, Tälleklint L: Incompetence of Roe Deer as reservoirs of the Lyme Borreliosis Spirochete. J Med Entomol. 1992, 29 (5): 813-818.PubMedView ArticleGoogle Scholar
- Randolph SE: Faune, climat et politique: Causes possibles aux récents accroissements des zoonoses à tiques. Ann Inst Pasteur (Paris). 2003, 16: 109-120.Google Scholar
- Jensen PM, Hansen H, Frandsen F: Spatial Risk Assessment for Lyme Borreliosis in Denmark. Scand J Infect Dis. 2000, 32 (5): 545-550.PubMedView ArticleGoogle Scholar
- Perkins SE, Cattadori IM, Tagliapietra V, Rizzoli AP, Hudson PJ: Localized deer absence leads to tick amplification. Ecology. 2006, 87 (8): 1981-1986.PubMedView ArticleGoogle Scholar
- Institut Bruxellois pour la Gestion de l'Environnement (IBGE): Le chevreuil, thematic leaflet. 1997, BruxellesGoogle Scholar
- Dean CB: Testing for Overdispersion in Poisson and Binomial Regression Models. J Amer Statistical Assoc. 1992, 87 (418): 451-457.View ArticleGoogle Scholar
- Fortin MJ, Dale M: Spatial Analysis: a Guide for Ecologists. 2005, Cambridge: Cambridge University PressGoogle Scholar
- Robinson TP: Spatial Statistics and Geographical Information Systems in Epidemiology and Public Health. Remote Sensing and Geographical Information Systems in Epidemiology. Edited by: Hay SI, Randolph SE, Rogers DJ. 2000, Oxford: Academic Press, 81-128. [Baker J, Muller R, Rollinson D (Series Editors): Advances in Parasitology, vol 47.]View ArticleGoogle Scholar
- Ribeiro PJ, Diggle PJ: geoR: a package for geostatistical analysis. R-News. 2001, 1 (2): 15-18.Google Scholar
- Ribeiro PJ, Christensen OF, Diggle PJ: geoR and geoRglm: Software for Model-Based Geostatistics. 3rd International Workshop on Distributed Statistical Computing (DSC 2003): 20–22 March 2003; Vienna. Edited by: Hornik K, Leisch F, Zeileis A. 2003, ViennaGoogle Scholar
- Anselin L: Under the hood: Issues in the specification and interpretation of spatial regression models. Agr Econ. 2002, 27 (3): 247-267.View ArticleGoogle Scholar
- Mallar CD: The estimation of simultaneous probability models. Econometrica. 1977, 45 (7): 1717-1722.View ArticleGoogle Scholar
- Jackson LE, Levine JF, Hilborn ED: A comparison of analysis units for associating Lyme disease with forest-edge habitat. Community Ecology. 2006, 7 (2): 189-197.View ArticleGoogle 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.