Determinants of tick-borne encephalitis in counties of southern Germany, 2001-2008
International Journal of Health Geographics volume 9, Article number: 42 (2010)
Tick-borne encephalitis (TBE) virus can cause severe symptoms in humans. The incidence of this vector-borne pathogen in humans is characterised by spatial and temporal heterogeneity. To explain the variation in reported human TBE cases per county in southern Germany, we designed a time-lagged, spatially-explicit model that incorporates ecological, environmental, and climatic factors.
We fitted a logistic regression model to the annual counts of reported human TBE cases in each of 140 counties over an eight year period. The model controlled for spatial autocorrelation and unexplained temporal variation. The occurrence of human TBE was found to be positively correlated with the proportions of broad-leafed, mixed and coniferous forest cover. An index of forest fragmentation was negatively correlated with TBE incidence, suggesting that infection risk is higher in fragmented landscapes. The results contradict previous evidence regarding the relevance of a specific spring-time temperature regime for TBE epidemiology. Hunting bag data of roe deer (Capreolus capreolus) in the previous year was positively correlated with human TBE incidence, and hunting bag density of red fox (Vulpes vulpes) and red deer (Cervus elaphus) in the previous year were negatively correlated with human TBE incidence.
Our approach suggests that a combination of landscape and climatic variables as well as host-species dynamics influence TBE infection risk in humans. The model was unable to explain some of the temporal variation, specifically the high counts in 2005 and 2006. Factors such as the exposure of humans to infected ticks and forest rodent population dynamics, for which we have no data, are likely to be explanatory factors. Such information is required to identify the determinants of TBE more reliably. Having records of TBE infection sites at a finer scale would also be necessary.
Tick-borne encephalitis (TBE) is the most important flavivirus infection of the central nervous system in Europe and Russia. The annual number of cases is estimated to be as high as 10,000 in Russia and about 3,000 in European countries [1–5]. Severe TBE infections caused by European virus strains typically take a biphasic course: After a short incubation period (usually 7-14 days, with extremes of 4-28 days), the first (viraemic) phase presents as an uncharacteristic influenza-like illness lasting 2-4 days (range 1-8 days) with fever, malaise, headache, myalgia, gastrointestinal symptoms, leukocytopenia, thrombocytopenia and elevated liver enzymes, and is often followed by a symptom-free interval of about one week (range 1-33 days). The second phase of TBE occurs in 20-30% of infected patients and is marked by four clinical features of different severity (meningitis, meningoencephalitis, meningoencephalomyelitis or meningoencephaloradiculitis) and the appearance of specific antibodies in the serum and cerebrospinal fluid. This is usually the time when patients with high fever and severe headache seek medical advice. The fatality rate in adult patients is less than 2%. However, severe courses of TBE infection with higher mortality and long-lasting sequelae, often affecting the patient's quality of life, are correlated with increased age [6–9].
As with most zoonotic diseases , TBE incidence in humans is characterised by considerable temporal (Figure 1) and spatial (Figure 2) heterogeneity . The main determinant of infection risk is the density of infected ticks, i.e. the product of pathogen prevalence in the ticks and tick density . Ixodes ricinus, the main vector of TBE virus (TBEV), has three distinct stages, larval, nymphal and adult ticks  whereas nymphs are responsible for the majority of tick bites in humans . The prevalence of TBEV in (nymphal) ticks depends on a combination of factors; the basic reproduction number (R0) of TBEV, which crucially depends on the occurrence of non-viraemic virus transmission between ticks co-feeding, especially on forest rodents [15–17]. Furthermore, the persistence of TBEV depends on a threshold value between the density of competent (e.g. yellow-necked mice Apodemus flavicollis) and incompetent hosts (e.g. red deer Cervus elaphus) . Possibly resulting from the species-specific reservoir competence, peaks in rodent populations in a given year are positively correlated with TBE incidence in humans in the succeeding year . Co-feeding of larval and nymphal ticks on rodents critically depends on seasonal activity synchrony of these immature tick stages. Since Randolph et al.'s study , which related larval and nymphal synchrony to a certain temperature regime in autumn, subsequent research attention has been directed towards a specific temperature regime during spring time  causing the seasonal synchrony of immature ticks. Tick density in a given landscape is primarily determined by the availability of suitable forest habitat . The findings of Allan et al. , relating to Lyme borreliosis risk, suggest that forest fragmentation might affect epidemiological risk. Tick density is further influenced by the abundance of host species such as roe (Capreolus capreolus) and red (Cervus elaphus) deer [23, 24] which feed large numbers ticks in all stages of development [25, 26]. Since estimates of wildlife population densities do not exist at the appropriate temporal and spatial resolution, we used hunting bag statistics as proxy for density. We also used hunting bag data of red fox (Vulpes vulpes), which was found to be positively correlated (with a time lag of one year) with TBE incidence in humans in Sweden . The explanatory variables investigated in this paper include hunting bag data for red deer, roe deer and red fox, land cover, spring warming increase, and an index of forest fragmentation. Specifically we apply statistical tests, which take account of spatial correlation, to assess whether these factors are related to TBE incidence in humans.
Annual symptomatic TBE infections in humans registered by the patients' place of permanent residence (county level) for the period 2001-2008 were obtained from the data base of the Robert Koch-Institut (SurvStat, http://www3.rki.de/SurvStat, 11/02/2009). We included the human population size for each county and year (provided by the federal statistical bureau) in our model. For a given county in a given year we modelled p, the proportion of reported clinical TBE cases y out of n inhabitants, assuming that TBE incidence in humans and infection risk are closely correlated.
A model for the TBE-count data
Since our dependent variable consists of counts of successes (y) and failures (n-y) in n trials, we used a generalized linear model (GLM) with binomial error distribution, i.e. a logistic regression model . The statistical package R  was used to estimate parameters of the model and to compute the test statistics. We model the probability that a person in a given county, and in a given year, is infected as a function of the explanatory variables.
Forest type and forest cover
Information on forest type and forest cover for each county was obtained from the CORINE dataset with a spatial resolution of 100 meters . We used the most recent available land cover information (year 2000) for all years, assuming that forest type and cover did not change substantially over this relatively short time period. For each county, we computed the proportion of coniferous, broad-leafed and mixed forest area.
Given the relative coarse resolution of the CORINE dataset, we used the 'largest patch index' (largest forest patch/total forest area) to approximate the level of forest fragmentation for each county .
We extracted mean diurnal monthly interpolated temperatures from 30 arc s resolution climate surface maps based on the 1950-2000 period . These data describe the spatial variation of temperatures which seems reasonable for our research question, given that temporal variation of temperatures is unlikely to explain the temporal variation of TBE incidence [33, 34]. Data pre-processing included a validity check of the temperature data and a re-projection to the WGS84 cartographic system in GIS ArcView version 9 (ESRI, Redlands, CA, USA). Mean interpolated temperatures (in °C) were computed for each county, for the months of January, February, March and April. From these we computed the temperature increase in spring as the mean temperature increase between February and April corrected for the mean temperature of January .
Hunting bag data
We obtained annual hunting statistics of game species from wildlife authorities of both states (Baden-Württemberg and Bavaria) for the period 2000/2001-2007/2008. To deal with the unequal sizes of the 140 counties, we standardised the hunting bag counts by dividing them by area available for wildlife, the latter being approximated by the area covered by agriculture and forest. For Bavaria, hunting bag data on roe deer were not available for the period 2000/2001. Considering the minor temporal variation in roe deer hunting bags, we used the roe deer bag data of the period 2001/2002 also for the missing period. In Bavaria, roe deer hunting is based upon a three years management plan but hunting bags are approximately evenly realised among years. In Baden-Württemberg, roe deer hunting is based on yearly management plans. Management plans for carnivores are not in place.
Unexplained spatial and temporal variation
Spatial autocorrelation  was incorporated by using the proportion of TBE cases per inhabitant in the neighbouring counties as an explanatory variable in the model. To control for unexplained temporal variation, we included each year as fixed factor (with year 2001 as reference) in the model. Summary statistics of dependent and explanatory variables are provided in Table 1.
Before fitting a full model, we tested the effect of single covariates on the probability of TBE infections in humans after correcting for unexplained spatial and temporal variation (Table 2). We tested for multicollinearity between the explanatory variables (Table 3) by computing the condition number. For model checking the deviance residuals were computed and displayed using GRASS GIS .
After correcting for unexplained spatial and temporal variation, five of the eight candidate variables were significantly (p < 0.05) related with TBE incidence in humans (Table 2). The candidate variables showed rather weak correlations, resulting in a condition number of 6.5, which indicates that multicollinearity is not of concern . The final model was fitted using all candidate variables.
The full model suggested that all candidate variables were significantly associated with TBE infection risk in humans (Table 4). The proportion of forest cover was positively correlated with probability of TBE infection; the effect of coniferous forest cover was stronger than that of broad-leafed and mixed forest cover (Figure 3a-c). When tested on its own, the proportion of broad-leafed forest cover was not significant (Table 2).
TBE risk was significantly negatively associated with the forest fragmentation index. This suggests that, taking all other variables as fixed, TBE risk in humans would decline if forest cover were more continuous (Figure 3d). Admittedly, we did not detect a significant relationship between the forest fragmentation index and TBE risk if tested without other explanatory variables (Table 2).
Wildlife population densities (measured indirectly as described above) were significantly correlated with TBE risk in humans in the following year (Figure 3e-g). The correlation with roe deer density was positive but those with red fox and red deer were negative. Except in the case of red deer, all coefficients in the full model have the same sign as those in the corresponding restricted models. TBE risk was significantly correlated with temperature increase during the spring months (Figure 3h) both in the restricted and in the full model (Tables 2, 4). Unexpectedly, the correlation is negative.
There is strong evidence of spatial autocorrelation (Figure 3i) and, unfortunately, also of unexplained temporal heterogeneity. The incidence of TBE infections varied substantially from year to year (see Figure 1) peaking in 2005 and 2006. Of course we would have wished that the model could explain these peaks in terms of the explanatory variables considered here. As is evident in Table 4 we were unable to account for this temporal heterogeneity with the covariates available to us.
With the above limitations in mind, the model (Null deviance = 3580.0 on 1119 degrees of freedom; residual deviance = 2494.8 on 1101 degrees of freedom) explained a large fraction of the observed variance (Cox-Snell-R2 = 0.96, Nagelkerke's R2 = 0.65). A histogram of the deviance residuals of the model (Figure 4) indicates that the distribution of the residuals is somewhat skewed, which is due to the discrete nature of the response; about 40% of the observed counts are equal to zero, and almost 75% of them are two or less. The spatial distribution of the residuals (Figure 5) varied considerably from year to year. There was a slight tendency to underestimate the TBE incidence in counties with few reported cases (south central part of the study area) and to overestimate the incidence in counties with many reported cases (western, northern and eastern regions of the study area) (Figures 2 and 5).
The model explaining TBE incidence in southern Germany contains eight variables (three related to forest cover, one related to forest fragmentation, three related to wildlife species population dynamics and one related to climatic conditions during spring time) plus one variable each to account for unobserved spatial and temporal variation. With one exception (spring warming) the results are broadly in line with previous research, highlighting the range of factors influencing TBE transmission dynamics and thus infection potential in humans [33, 34, 38]. The fact that the factor "year" had to be included in the model constitutes a failure of the model. Other factors, especially those taking account of human activities, would need to be considered to explain the remaining and substantial temporal heterogeneity. The admittedly incomplete model does, however, consider the influence of several explanatory variables simultaneously, rather than one variable or a few variables.
The proposed model suggests that the proportion of coniferous forests had a stronger effect on TBE risk for humans compared to other forest types. Results from field studies provide equivocal results concerning the epidemiological risk [densities of (infected) ticks] according to forest type, but generally broad-leafed forest is considered as a suitable habitat for ticks. In North American forests, black-legged ticks (Ixodes scapularis) are more abundant in deciduous than in coniferous woodlands . In Scotland, tick densities were highest in coniferous forests compared to deciduous woodland and pastures . In Denmark, no differences in Ixodes ricinus tick densities between spruce (conifer) forests and non spruce forests were found . It is, however, unlikely that forest type itself accounts directly for tick density patterns. Potentially, forest type indirectly affects tick densities by providing appropriate moisture and host species for ticks or by sustaining more TBE competent reservoir species relative to incompetent ones [18, 23, 33], or forest types are unevenly frequented by humans.
In addition to forest extent, we hypothesized that forest fragmentation would enhance TBE risk either by increasing density of infected ticks  or by enhancing contact rates between humans and ticks. The effect of forest fragmentation was unclear. Despite showing weak correlations with other variables (Table 3), which indicates that it is not confounded with other variables, the sign of the coefficient in the simple model was the opposite of that in the full model. Yet, removing this variable from the full model did not alter the signs of other coefficients.
Hunting bag densities of three game species were found to be significantly associated with TBE risk in humans. Since hunting bags are subject to management plans (deer) and/or influenced by hunting effort (red fox), these observations should be regarded cautiously . Yet, "important ecological questions simply have to be addressed on the right scale - which often means an uncomfortable large scale - even if that means a certain degree of imprecision" . There is growing empirical evidence that roe deer populations increase the epidemiological risk of ticks , Lyme disease risk  and TBE risk [, , this study]. Regarding the role of red deer in determining TBE risk, our results are inconclusive (changing signs simple model vs. full model) but do not necessarily contradict findings from Italy, where a non-significant effect was found between red deer and TBE incidence in humans .
The negative relationship found in this study between red fox hunting bags and human TBE cases in the succeeding year is in contrast to a study in TBE endemic areas in Sweden, where this relationship was found to be positive . These contradictory findings will need to be clarified in future studies on the specific role of foxes in the epidemiology of TBEV.
Unfortunately, reliable data on forest rodent dynamics are not available at the appropriate temporal and spatial scale, and could thus not enter our model. Further attempts should be made to test the hypothesis that high rodent (especially Apodemus flavicollis) densities translate into high nymphal densities in the following year  and hence to elevate TBE incidence in humans . Such a finding might be of use as an effective early warning system for public health.
Unexpectedly, the direction of the regression coefficient of the spring warming variable is not in line with the claimed importance of simultaneous activity of larval and nymphal ticks for TBE maintenance . That theory suggests that fast warming in spring is required to allow nymphs (critical temperature ~7°C) and larvae (critical temperature ~10°C) to feed synchronously  and to transmit pathogens among the tick population . However, it should be emphasised that spring warming has so far been used only to explain the distribution of TBE foci and not the variable incidence.
As in most epidemiological studies, we based our analysis on the assumption that TBE incidence is strongly correlated with the epidemiological risk of TBE (i.e. the density of infected nymphs). Although the model provides a reasonable fit, some potentially important factors were not investigated due to lack of data. Included here are differing virulence of TBE, variable exposure of humans to ticks, differing immune responses, demography and TBE immunisation coverage. E.g. Kimmig et al.  found discrepancies in spatial distribution of clinical TBE cases and prevalence of TBE antibodies in sera of persons at risk. In addition, demographic analyses indicate that men above the age of 35 years are disproportional represented (Figure 6) among the persons with clinical TBE symptoms, which might be related to higher exposure risk to ticks and/or to differing immune responses . Furthermore, locations associated with the counts available to us are the patients' place of permanent residence, and not the place of infection. Although it is likely that most infections occur in the infected person's home county, exceptions in this respect are clearly a potential source of bias. Due to data availability this analysis is based on a rather coarse scale [cf. 42]. Counties do not necessarily represent homogenous units with respect to environmental and ecological conditions. Generally forests, animal communities and interpolated temperatures are very heterogeneous within a given county. Thus any analysis based on values that have been averaged over variable conditions cannot be expected to yield precise results. Consequential we strongly recommend recording the suspected places of TBE infection as accurately as possible i.e. at the resolution of the forest patch.
Our modelling approach focussed on ten biotic and abiotic factors possibly influencing the epidemiological risk of TBE. Except for the influence of spring warming, the results are in line with previously published findings. The presented approach might therefore be useful for predicting at least the potential direction, and approximate magnitude of TBE risk, as a function of ecological change (land cover change, wildlife population dynamics). To make further progress, a higher spatial resolution is needed in order to make more reliable interferences about the relationships between eco-environmental factors and TBE incidence. Moreover, neither this analysis nor the in-depth analyses of yearly temperature variability [33, 34, 48] could entirely explain the temporal dynamics of tick abundance and TBE in the past decade. This underlines previous notions that besides the epidemiological risk, the contact rate probability between infected ticks and humans is of crucial importance and should be incorporated in models for tick-borne diseases [11, 43, 48].
Pugliese A, Beltramo T, Torre D: Emerging and re-emerging viral infections in Europe. Cell Biochem Funct. 2007, 25 (1): 1-13. 10.1002/cbf.1342.
Günther G, Haglund M: Tick-borne encephalopathies: epidemiology, diagnosis, treatment and prevention. CNS Drugs. 2005, 19 (12): 1009-32. 10.2165/00023210-200519120-00004.
Charrel RN, Attoui H, Butenko AM, Clegg JC, Deubel V, Frolova TV, Gould EA, Gritsun TS, Heinz FX, Labuda M, Lashkevich VA, Loktev V, Lundkvist A, Lvov , Mandl CW, Niedrig M, Papa A, Petrov VS, Plyusin A, Randolph S, Süss J, Zlobin VL, de Lamballerie X: Tick-borne virus diseases of human interest in Europe. Clin Microbiol Infect. 2004, 10 (12): 1040-55. 10.1111/j.1469-0691.2004.01022.x.
Gritsun TS, Lashkevich VA, Gould EA: Tick-borne encephalitis. Antiviral Res. 2003, 57 (1-2): 129-46. 10.1016/S0166-3542(02)00206-1.
Süss J: Tick-borne encephalitis in Europe and beyond - the epidemiological situation as of 2007. Euro Surveill. 2008, 13 (26): pii=18916
Holzmann H: Diagnosis of tick-borne encephalitis. Vaccine. 2003, 21 (Suppl 1): S36-40. 10.1016/S0264-410X(02)00819-8.
Kunze U, Baumhackl U, Bretschneider R, Chmelik V, Grubeck-Loebenstein B, Haglund M, Heinz F, Kaiser R, Kimmig P, Kunz C, Kunze M, Mickiene A, Misic-Majerus L, Randolph S, Rieke B, Stefanoff P, Süss J, Wimmer R, International Scientific Working Group on Tick-borne encephalitis: The Golden Agers and Tick-borne encephalitis: Conference report and position paper of the International Scientific Working Group on Tick-borne encephalitis. Wiener Med Wochenschr. 2005, 155 (11-12): 289-94. 10.1007/s10354-005-0178-0.
Kaiser R: Tick-borne encephalitis (TBE) in Germany and clinical course of the disease. Int J Med Microbiol. 2002, 291 (Suppl 33): 58-61. 10.1016/S1438-4221(02)80012-1.
Mickiene A, Laiskonis A, Günther G, Vene S, Lundkvist A, Lindquist L: Tickborne encephalitis in an area of high endemicity in Lithuania: disease severity and long-term prognosis. Clin Infect Dis. 2002, 35 (6): 650-658. 10.1086/342059.
Jones E, Patel NG, Levy MA, Storeygard A, Balk D, Gittleman JL, Daszak P: Global trends in emerging infectious diseases. Nature. 2008, 451: 990-994. 10.1038/nature06536.
Randolph SE, Sumilo D: Tick-borne encephalitis in Europe: dynamics of changing risk. Emerging Pests and Vector-borne Diseases in Europe. Edited by: Takken W, Knols B. 2007, Wageningen, University Publishers, 187-206.
Ostfeld RS, Canham CD, Oggenfuss K, Winchcombe RJ, Keesing F: Climate, deer, rodents, and acorns as determinants of variation in Lyme-disease risk. PLoS Biol. 2006, 4 (6): e145-10.1371/journal.pbio.0040145.
Randolph SE: Tick ecology: processes and patterns behind the epidemiological risk posed by ixodid ticks as vectors. Parasitology. 2004, 129 (Suppl 1): S37-S65. 10.1017/S0031182004004925.
Vassalo M, Perez-eid C: Comparative behavior of different life cycle stages of Ixodes ricinus (Acari: Ixodidae) to human-produced stimuli. J Med Entomol. 2002, 39: 234-236. 10.1603/0022-2585-39.1.234.
Labuda M, Nuttal PA, Kozuch O, Eleckova E, Zuffova E, Sabo A: Non-viraemic transmission of tick-borne encephalitis virus: a mechanism for arbovirus survival in nature. Experientia. 1993, 49: 802-805. 10.1007/BF01923553.
Randolph SE, Gern L, Nutall PA: Co-feeding ticks: Epidemiological significance for tick-borne pathogen transmission. Parasitol Today. 1996, 12 (12): 472-479. 10.1016/S0169-4758(96)10072-7.
Hartemink NA, Randolph SE, Davis SA, Heesterbeek : The basic reproduction number for complex disease systems: defining R0 for tick-borne infections. Am Nat. 2008, 171 (8): 743-754. 10.1086/587530.
Rosà R, Pugliese A: Effects of tick population dynamics and host densities on the persistence of tick-borne infections. Math Biosci. 2007, 208: 216-240. 10.1016/j.mbs.2006.10.002.
Zeman P, Benes C: A tick-borne encephalitis ceiling has moved upwards during the last 30 years: possible impact of global warming?. Int J Med Microbiol. 2004, 293 (Suppl 37): 48-54.
Randolph SE, Green RM, Peacey MF, Rogers DJ: Seasonal synchrony: the key to tick-borne encephalitis foci identified by satellite data. Parasitology. 2000, 121: 15-23. 10.1017/S0031182099006083.
Racz GR, Ban E, Ferenczi E, Berencsi G: A simple spatial model to explain the distribution of human tick-borne encephalitis cases in Hungary. Vector-Borne Zoonot. 2006, 6: 369-378. 10.1089/vbz.2006.6.369.
Allan BF, Keesing F, Ostfeld RS: Effect of forest fragmentation on lyme disease risk. Conserv Biol. 2003, 17 (1): 267-272. 10.1046/j.1523-1739.2003.01260.x.
Jensen PM, Hansen H, Frandsen F: Spatial risk assessment for Lyme borreliosis in Denmark. Scand J Infect Dis. 2000, 32: 545-550. 10.1080/003655400458857.
Gilbert L: Altitudinal patterns of tick and host abundance: a potential role for climate change in regulating tick-borne diseases?. Oecologia. 2010, 162: 217-225. 10.1007/s00442-009-1430-x.
Kiffner C, Lödige C, Alings M, Vor T, Rühe F: Abundance estimation of Ixodes ticks (Acari: Ixodidae) on roe deer (Capreolus capreolus). Exp Appl Acarol. 2010, 52: 73-84. 10.1007/s10493-010-9341-4.
Vor T, Kiffner C, Niedrig M, Hagedorn P, Rühe F: Tick burden on European roe deer (Capreolus capreolus L.). Exp Appl Acarol. 2010, 51: 405-417. 10.1007/s10493-010-9337-0.
Haemig PD, Lithner S, Sjöstedt de Luna S, Lundkvist Å, Waldenström J, Hansson L, Arneborn M, Olsen B: Red fox and tick-borne encephalitis (TBE) in humans: Can predators influence public health?. Scand J Infect Dis. 2008, 40: 527-532. 10.1080/00365540701805446.
Dobson AJ, Barnett AG: An introduction to generalized linear models. 2008, Boca Raton, Chapman & Hall/CRC
R Development Core Team: R: A language and environment for statistical computing. 2008, R Foundation for Statistical Computing, Vienna, Austria,http://www.R-project.org
Anonymous: CORINE Land Cover. 2004, Federal Environmental Agency, German Aeospace Center
Haines-Young R, Chopping M: Quantifying landscape structure: a review of landscape indices and their application to forested landscapes. Prog Phys Rev. 1996, 20: 418-445. 10.1177/030913339602000403.
Hijmans RJ, Cameron SE, Parra JL, Jones PG, Jarvis A: Very high resolution interpolated climate surfaces for global land areas. Int J Climat. 2005, 25: 1965-1978. 10.1002/joc.1276.
Rizzoli A, Hauffe HC, Tagliapietra V, Neteler M, Rosà R: Forest structure and roe deer abundance predict tick-borne encephalitis risk in Italy. PloS ONE. 2009, 4: e4336-10.1371/journal.pone.0004336.
Sumilo D, Asokliene L, Bormane A, Vasilenko V, Golovljova I, Randolph SE: Climate change cannot explain the upsurge of tick-borne encephalitis in the Baltics. PLoS ONE. 2007, 6: e500-10.1371/journal.pone.0000500.
Wimberly MC, Baer AD, Yabsley MJ: Enhanced spatial model for predicting the geographic distributions of tick-borne pathogens. Int J Health Geogr. 2008, 7: 15-10.1186/1476-072X-7-15.
GRASS Development Team: Geographic Resources Analysis Support System (GRASS) Software,l 2010. Open Source Geospatial Foundation,http://grass.osgeo.org
Mason CH, Perreault WD: Collinearity power and interpretation of multiple regression analysis. J Marketing Res. 1991, 28: 268-280. 10.2307/3172863.
Estrada-Pena A: Tick-borne pathogens, transmission rates and climate change. Front Biosci. 2009, 14: 2674-2780. 10.2741/3405.
Ginsberg HS, Zhioua E, Mitra S, Fischer J, Buckley PA, Verret F, Underwood HB, Buckley FG: Woodland type and spatial distribution of nymphal Ixodes scapularis (Acari : Ixodidae). Environ Entomol. 2004, 33 (5): 1266-1273. 10.1603/0046-225X-33.5.1266.
Walker AR, Alberdi MP, Urquhart KA, Rose H: Risk factors in habitat of the tick Ixodes ricinus influencing human exposure to Ehrlichia phagocytophila bacteria. Med Vet Entomol. 2001, 15: 40-49. 10.1046/j.1365-2915.2001.00271.x.
Ranta E, Lindström J, Lindén H, Helle P: How reliable are harvesting data for analyses of spatio-temporal population dynamics?. Oikos. 2008, 117: 1461-1468. 10.1111/j.0030-1299.2008.16879.x.
May RM: Crash test for the real. Nature. 1999, 398: 371-372. 10.1038/18769.
Linard C, Lamarque P, Heyman P, Ducoffre G, Luyasu V, Tersago K, Vanwambeke SO, Lambin EF: Determinants of the geographic distribution of Puumala virus and Lyme borreliosis infections in Belgium. Int J Health Geogr. 2007, 6: 15-10.1186/1476-072X-6-15.
Hudson PJ, Rizzoli A, Rosà R, Chemini C, Jones LD, Gould EA: Tick-borne encephalitis virus in northern Italy: molecular analysis, relationships with density and seasonal dynamics of Ixodes ricinus. Med Vet Entomol. 2001, 15: 304-313. 10.1046/j.0269-283x.2001.00317.x.
Rosà R, Pugliese A, Ghosh M, Perkins SE, Rizzoli A: Temporal variation of Ixodes ricinus intensity on the rodent host Apodemus flavicollis in relation to local climate and host dynamics. Vector-Borne Zoonotic Dis. 2007, 285-295. 10.1089/vbz.2006.0607.
Kimmig P, Oehme R, Backe H: Epidemiologie der Frühsommer-Meningoenzephalitis (FSME) und Lyme-Borreliose in Südwestdeutschland. Ellipse. 1998, 14 (4): 95-105.
Sheldon BC, Verhulst S: Ecological immunology: costly parasite defences and trade-offs in evolutionary ecology. TREE. 1996, 11: 317-321.
Randolph SE, Asokliene L, Avsic-Zupanc T, Bormane A, Burri C, Gern L, Golovljova I, Hubalek Z, Knap N, Kondrusik M, Kupca A, Pejcoch M, Vasilenko V, Zygutiene M: Variable spikes in tick-borne encephalitis incidence in 2006 independent of variable tick abundance but related to weather. Parasites & Vectors. 2008, 1: 44-
This study was funded by a research grant from the Federal Ministry of Education and Research (BMBF grant 1363120: Emerging arthropod-borne-viral infections in Germany: Pathogenesis, diagnostics and surveillance).The Bavarian State Ministry of Alimentation, Agriculture and Forestry and the Wildlife Research Unit of the State Baden-Württemberg kindly provided hunting bag data. Data on human population were kindly supplied by the German Federal Statistical Office. We sincerely thank F. Herschleb for assistance in digitalising hunting bag data and H. Reinecke for generating the map of the residuals. We thank four anonymous reviewers whose comments greatly improved previous drafts of this manuscript.
The authors declare that they have no competing interests.
CK designed the experiment, acquired the data, performed the statistical analysis, interpreted the results and wrote the manuscript. WZ designed the experiment, performed the statistical analysis, contributed to the interpretation of the results and the writing of the MS. PS performed the GIS analyses. FR acquired the data, contributed to the interpretation of the results and the writing of the MS. TV, PH and MN contributed to the interpretation of the results and the writing of the MS. All authors read and approved the final manuscript.
Authors’ original submitted files for images
About this article
Cite this article
Kiffner, C., Zucchini, W., Schomaker, P. et al. Determinants of tick-borne encephalitis in counties of southern Germany, 2001-2008. Int J Health Geogr 9, 42 (2010). https://doi.org/10.1186/1476-072X-9-42