INTRODUCTION
Many zoonoses are listed as recent emerging or re-emerging infectious diseases [Reference Woolhouse and Gowtage-Sequeria1]. Often it is the synergistic relationship between anthropogenic and ecological parameters which leads to an increased number of reported infections. Recently, 15 emerging endemic zoonoses with potential human health risks were reported within Europe during the period 2000–2006 and most of these zoonoses were related to changes in ecological variables, e.g. increasing host population abundance, climate change or host migration [Reference Vorou, Papavassiliou and Tsiodras2, Reference Mills3]. The rodent-borne hantaviruses (genus Hantavirus, family Bunyaviridae) are part of this list, partly due to the abrupt increase in nephropathia epidemica (NE) incidence during 2005 in Belgium, Germany, France, Luxembourg and The Netherlands [Reference Mailles4, Reference Heyman5]. NE, caused by Puumala virus (PUUV) infection, is a widespread zoonosis throughout Europe. The disease is marked by a range of symptoms going from abrupt onset of fever and headache towards vomiting and acute renal failure [Reference Vapalahti6]. The distribution of NE is mostly related to the occurrence of the predominant reservoir host Myodes glareolus (bank vole). The virus is transmitted among bank voles through virus particles shed in their excretions. Transmission occurs directly through aggressive or sexual contacts among voles or indirectly by inhaling virus aerosols originated from contaminated urine or faeces. The latter is also assumed to be the main transmission route to accidental hosts, e.g. humans, who are, as far as we know, dead-end hosts. If humans get infected with PUUV, it can take up to several weeks to develop symptoms due to PUUV infection [Reference Vapalahti6].
In the present study, we focus on NE in Belgium. A distinct difference in NE incidence is present throughout this small country with rare occurrence of NE in the north, yet high incidences in south Belgium [Reference Heyman7, Reference Linard8]. Additionally, there is a lot of annual variation between ‘epidemic’ and ‘non-epidemic’ years. Until 2000 a 3-year epidemic cycle was described, after which it abruptly changed into a 2-year cycle that currently remains [Reference Heyman5, Reference Heyman7]. Furthermore, in 2005, 365 human cases were reported, this was the highest incidence ever recorded after the 1996 outbreak (224 cases). Previous studies have shown that environmental conditions which favour the indirect transmission route, combined with specific human activities, could play a major role in explaining the spatial distribution patterns of both bank vole and human infection [Reference Linard8–Reference Linard10]. Temporal variation in PUUV infection risk, however, was rarely assessed. When it comes to actual risk of infection as previously described [Reference Davis, Calvet and Leirs11], it is hypothesized that it is the abundance of infected voles at a given time that remains predominantly crucial. This seems to hold as increased densities of infected bank voles have been related to outbreak years on several occasions [Reference Olsson9, Reference Escutenaire12]. A recent model, however, suggests importance of the abundance of ‘newly’ infected voles [Reference Sauvage, Langlais and Pontier13], due to the higher rate of infectious virus shedding in their excretions [Reference Gavrilovskaya14]. Thus the model states that only extreme changes in bank vole demography can affect PUUV infection risk towards humans. Rodent dynamics are influenced by many aspects in the environment. Unlike the situation in Northern Europe where high bank vole abundance is mostly explained by density-dependent effects on predators and vole maturation [Reference Yoccoz15], major peaks in bank vole abundance in the West–Central European region have almost solely been related to mast years in tree seed production [Reference Jensen16–Reference Crespin19]. High seed production affects bank vole abundance in several ways. It can improve winter survival, elongate the bank vole breeding period and even induce winter breeding [Reference Pucek17–Reference Crespin19]. Reproduction of bank voles is regulated by a density-dependent process involving social interaction and occupancy of female breeding territories. In case of high resource availability, female territories become smaller. Hence, during high tree seed production, new unoccupied territories are created in autumn. As a consequence, young female subadults, who would normally delay breeding until the next spring, will still be able to acquire their own breeding territory and initiate reproduction in autumn. Due to the synergistic effect of higher winter survival, extension of the breeding period and a high proportion of breeding females, the population density can remain high until the next spring [Reference Pucek17–Reference Bujalska and Saitoh21]. It is generally believed that rodent peaks are caused by massive production of seeds by native oak or beech and relate to subsequent NE outbreaks in Belgium [Reference Heyman7, Reference Escutenaire12, Reference Clement22], although this was never thoroughly investigated. Furthermore, if tree seeding is indeed the driving force behind NE outbreaks in Belgium, the ultimate explanations for the recent patterns in NE epidemiology must be sought in those factors that trigger tree seed production. Forestry studies indicate that, in case of sufficient tree resources, it is mainly climate that affects tree fructification. For Fagus sylvatica, Quercus robur and Q. petraea there is a period of dormancy between flower bud initiation (early summer/summer) and flowering (pollination)/seeding (respectively spring and autumn in the following year). Both flower bud induction and flowering (pollination) can be affected by local weather conditions and there is no overall consensus on the dominant factors due to inter-species differences [Reference Johnson, Shifley and Rogers23–Reference Piovesan and Adams25]. Only one factor seems to be consistent through most of the literature and that is summer conditions. Especially for beech there is a general agreement that warm or dry summer conditions during flower bud formation is a predominant factor for high seed production in the next year [Reference Piovesan and Adams25, Reference Övergaard, Gemmel and Karlsson26]. Climate conditions during high bank vole abundances, however, can also affect human behaviour and transmission to humans as humid and colder conditions could lead to better survival of the virus in the surroundings, hence, increasing the probability of transmission between voles and to humans [Reference Linard8, Reference Kallio27].
The aim of this study is first to evaluate the relationship between NE numbers and categories of tree seed production in Belgium. Second, to test whether known climatic factors that trigger high seed production can be directly related to human NE incidence or if climatic conditions during the period of human infection cause a higher number of NE cases.
METHODS
NE data
NE data were provided by the Scientific Institute of Public Health (IPH) and the Reference Laboratory for hantavirus infection. The surveillance of human hantavirus infections in Belgium is carried out through a network of sentinel laboratories, coordinated by the IPH. These sentinel laboratories send positive samples to the Reference Laboratory for confirmation. Reliable annual NE case numbers are available from 1994 until 2007.
Tree seed production categories
Categories of seed production of beech and native oak species (Q. robur, Q. petraea) were provided by the Tree Seed Centre of the Ministry of the Walloon Region. Since 1995, this centre has collected seeds for economic purposes. Due to their access to independent information sources, they are able to make reliable global assessments of tree seed production throughout southern Belgium. These assessments are based on tree flowering, pollination and seed production evaluations in the field, which are related to information from the forest service that gives the potential of seed collection in their selected stands. Tree seed production for each tree species is divided into four categories. Very good years (the species is fruiting throughout the Walloon territory and practically all trees are bearing seed in high quantities), Good years (the species is fruiting throughout the territory, but the trees are bearing much less seed and some trees do not fruit), Moderate years (there is a reduced number of trees bearing seeds and sometimes only located in a portion of the territory) and Low years (years without fructification in significant quantities).
Climate data
The meteorogical data were provided by the Royal Meteorological Institute from Belgium (Ukkel). Average monthly and seasonal temperature and precipitation values were calculated from daily weather data. Meteorological seasons were used: winter (December, January, February), spring (March, April, May), summer (June, July, August), autumn (September, October, November).
Statistical analyses
Generalized Linear Models (GLM) in SAS® version 9.1 [28] were used to investigate the effect of beech and native oak tree seed categories on the number of NE cases in the following year. Log-linear models with log link function and both a Poisson and negative binomial distribution were compared using the likelihood ratio model comparison method [Reference Agresti29]. The model with a negative binomial error distribution fitted our data significantly better (P<0·001). Pairwise comparisons were assessed using least square means.
Similar GLM analyses were used to evaluate the relationship between annual NE occurrence in year t and ecologically relevant monthly and seasonal values of average daily temperature and precipitation from summer of year t−3 until autumn of year t. Monthly data were only used for those months or combinations that have been shown to be ecologically relevant in previous studies; June–July (year t−2), June–July (year t−3), July (year t−2), April (year t−1), May (year t−1), June (year t−2), June (year t−3), July (year t−2), July (year t−3), September (year t−2) [Reference Pucek17, Reference Johnson, Shifley and Rogers23–Reference Övergaard, Gemmel and Karlsson26].
To take into account first-order autocorrelation, the effect of the preceding year t−1 number of NE cases was added to the model and independence among explanatory variables was tested using Spearman correlations. All variables were modelled in log-linear models with log link function and both a Poisson and negative binomial distribution. Again, likelihood ratio model comparison showed that models with a negative binomial response distribution fitted best (P<0·001). Overall results however did not differ significantly. Final model selection was based on AIC (Akaike Information Criterion) statistics [Reference Agresti29].
RESULTS
Variation in NE case numbers and tree seed production
Figure 1 shows the number of human PUUV cases (NE) from 1995 to 2007, demonstrating an abrupt increase in the number of NE cases in 2005. This pattern of increased human infection did not cease in the following years, 160 cases were confirmed in 2006 and again 290 in 2007. Remarkably, as prior to 2005, each epidemic year (>100 NE cases annually) was followed by at least one non-epidemic year with <60 annual NE cases.
Fig. 1. Annual number of nephropathia epidemica (NE) cases and the category of seed production in the whole southern Belgian territory of beech, F. sylvatica (white boxes) and native oak, Q. robur and Q. petraea (grey boxes) in the respective years. Categories of fructification are ordered from 0 to 3 (0=low, 1=moderate, 2=good, 3=very good).
The categorical evaluation of annual beach and native oak fruiting in southern Belgium is indicated in Figure 1. We can observe that each peak year in NE incidence is preceded by a year with high seed production of at least native oak, beech or both. Separate analyses show a significant main effect of both beech (F 2,9=9·85, P=0·007) and oak fruiting (F 3,8=12·79, P=0·005). Very good or good seed production resulted in a significantly higher number of NE cases in the following year compared to low seed production for both beech (0–2: χ12=7·67, P=0·006; 0–3: χ12=12·37, P=<0·001) and native oak (0–1: χ12=9·68, P=0·002; 0–2: χ12=13·11, P=<0·001; 0–3: χ12=7·21, P=<0·007). No significant differences were found among the higher seed production categories.
Climate
Only a limited number of significant explanatory variables were found. All showed a positive relationship to NE case numbers and most of the parameters are related to the summer season before tree seeding, i.e. 2 years before NE occurrence; summer temperature year t−2 (χ12=8·01, P=0·005); June–July temperature year t−2 (χ12=8·43, P=0·004); June temperature year t−2 (χ12=4·30, P=0·038). All of these parameters were correlated (P<0·01). Moreover, June temperature year t−3 was significant (χ12=4·76, P=0·029). Additionally, we also found a significant effect of temperature during autumn, the year preceding human case occurrence; autumn temperature year t−1 (χ12=5·35, P=0·021). The number of NE cases in year t−1 did not affect the NE cases in year t significantly throughout the time-series (χ12=0·01, P=0·936). All combinations of significant independent parameters were made and the model with the lowest AIC value included summer temperature year t−2 and autumn temperature year t−1 (ΔAIC>13) (Table 1, Figs 2, 3].
Fig. 2. Contour plot of the annual number of NE cases against the average summer temperature (°C) 2 years before and average autumn temperature (°C) 1 year before nephropathia epidemica (NE) infection. Contours are based on the selected model explained in Table 1. Intervals of raw NE data are given in the key.
Fig. 3. Scatterplot of the observed vs. the predicted number of NE cases (year) by the selected model with average summer temperature (year −2) and average autumn temperature (year −1). Model details are given in Table 1.
Table 1. Parameter estimates (log scale) of the best fitting model of annual nephropathia epidemica cases (year) and climate factors; average summer temperature year −2 (°C) and average autumn temperature year −1 (°C)
DISCUSSION
Our analyses show that regional homogeneous high seed production of both beech and native oak are closely related to increased numbers of NE cases in Belgium. These results, together with the known relationships between both bank vole abundance/high seed production and bank vole abundance/PUUV infection risk [Reference Olsson9, Reference Escutenaire12, Reference Jensen16, Reference Verhagen, Leirs and Verheyen18, Reference Clotfelter30], make us confident in stating that, if either beech, native oak or both show a very high seed production during autumn, this can be regarded as an early warning tool for public health workers. Centres like the Tree Seed Centre of the Walloon Region can therefore play a key role in alerting the IPH during the autumn seasons preceding NE outbreaks.
Furthermore, we gained insight in the ultimate processes changing NE epidemiology in Belgium and the rest of Western Europe. These processes seem to be directly related to climate. Although climate variables typically have a subtle or no measurable effect on disease systems, hantavirus infection has been shown to be very sensitive to climate variability [Reference Tamerius31]. There is the well known example of Sin nombre virus infection in the United States and its close link to increased precipitation patterns [Reference Engelthaler32, Reference Hjelle and Glass33]. Even though it remains under discussion, a higher vegetation growth due to humid conditions causing increased resource availability would lead to higher deer mouse densities and a consequent higher infection risk. This process is known as the Trophic Cascade Hypothesis [Reference Yates34, Reference Glass35]. Our system seems highly similar although other parameters link climate and human infection patterns. We found a significant effect of average autumn temperature 1 year before human infection occurrence. This finding is difficult to relate to tree seeding, during this period seeds fall from the trees. Higher temperatures from summer to autumn, however, have been related to increased bank vole reproduction. The underlying hypothesis here would be the importance of green biomass and body condition. During the normal reproductive season, herbaceous plants constitute a high percentage of the bank vole diet. Normally green foliage decreases at the beginning of autumn, higher temperatures, however, can postpone this effect [Reference Pucek17]. Furthermore, higher temperatures mean a lower use of resources to maintain the high metabolic rate of bank voles. The positive effect of average summer temperature in the 2 years preceding NE occurrence is directly related to the period of flower bud initiation in oak and beech. This is the one factor on which most studies agree that it is essential for good mast crops [Reference Övergaard, Gemmel and Karlsson26]. Therefore, higher temperatures during summer constitute an indirect link to higher PUUV infection risk and higher numbers of NE cases. Shorter mast year intervals and increasing mast crops have been observed in recent years [Reference Övergaard, Gemmel and Karlsson26]. The changing epidemiological cycles can thus be linked to these shorter intervals between masting and shorter intervals between more than average warm summers.
Finally, the warm summer temperatures could not only explain NE occurrence on a temporal scale, but also on a spatial scale. During 2003 an extreme climatic anomaly was registered in Europe. There was a summer heatwave with average summer temperatures exceeding the 1961–1990 mean by up to 5 s.d. and the summer might have been the warmest since 1540 [Reference Zaitchik36]. This anomaly was most pronounced in West–Central Europe and affected local deciduous forests in Germany and part of France and Belgium [Reference Reichstein37]. Based on the findings in our study and observations, we therefore conclude that the extreme PUUV infection outbreak in 2005, which was most reported in those countries where the focus of the climate anomaly was found [Reference Heyman5, Reference Abu Sin38], is a result of this heatwave occurring 2 years earlier in summer 2003. Based on the average temperature data of summer 2006 (18·97°C) and autumn 2007 (13·90°C), we again predict a very high number of NE cases in 2008; this seems to hold as during the first 5 weeks of 2008 already 30 NE cases have been reported in Belgium (IPH).
High seed production by beech and native oak and elevated average summer temperatures followed by warmer autumn conditions in the next year should be considered as early warning indicators for NE outbreaks in Western Europe. The relationship with higher temperatures is especially valuable. Since future climate change scenarios predict higher temperatures in Europe, we suggest regarding PUUV as an increasing health threat.
ACKNOWLEDGEMENTS
We thank Stefan Van Dongen and Laurent Crespin from the University of Antwerp for help with the statistics. We thank two anonymous reviewers for their comments. K. Tersago works as an Academic Assistant at the University of Antwerp. This research was partially funded by EU grant GOCE-2003-010284 EDEN and the paper is catalogued by the EDEN Steering Committee as EDEN0097 (http://www.eden-fp6project.net/).
DECLARATION OF INTEREST
None.
