INTRODUCTION
Hantaviruses of the family Bunyaviridae, genus Hantavirus are predominantly rodent-borne pathogens that may cause haemorrhagic fever with renal syndrome (HFRS) in Eurasia [Reference Vaheri1]. Humans mainly become infected by inhaling aerolized rodent excreta [Reference Vaheri1, Reference Hardestam2]. Five hantaviruses circulate among rodents in Europe, but the majority of human HFRS cases are caused by Puumala virus (PUUV) [Reference Vaheri1, Reference Olsson, Leirs and Henttonen3]. The bank vole (Myodes glareolus) is the reservoir for PUUV and is a widely spread vole species in Europe although absent in many areas around the Mediterranean [Reference Vaheri1]. HFRS is a reportable disease in many countries in the European Union. However, most cases are reported from Finland (~70% of cases), Sweden and forest-rich regions of Belgium and Germany [Reference Faber4–Reference Clement8].
PUUV infections show pronounced seasonality in Finland with highest incidence rates in late summer and early winter [Reference Makary7, Reference Brummer-Korvenkontio9–Reference Kallio11]. The bank vole undergoes 3- to 4-year population cycles, seemingly driven by specialist predators, which coincide with human HFRS epidemics [Reference Vaheri1]. However, it has been unclear how seasonality and long-term cycles of these epidemics vary regionally and temporally.
The aim of this study was to report recent trends in PUUV infections in Finland and to evaluate in detail how different regions contribute to the general epidemiological pattern in Finland; whether there are regional differences in seasonality and long-term cycles and whether the patterns have changed over time.
METHODS
Since 1995, all Finnish clinical microbiology laboratories have reported serological tests that were positive for PUUV to the National Infectious Disease Register (NIDR). Each notification includes data on age, sex and, as well as date of diagnostic specimen, place of residence and treatment. Multiple notifications of persons with the same date of birth, sex and place of treatment received within a 12-month period were combined as one report; since 2004, this was done by a unique national identity number.
A case was defined as a person who had a positive serological test for PUUV infection which was reported to the NIDR from 1 April 1995 to 31 March 2014. We considered each 12-month period from April until the following March as a ‘seasonal year’, based on previous observation: the lowest number of cases occurred from March to April in each year. Data from the Population Information System for 1995–2014 were used as a denominator to calculate incidence rates. For regional analyses, we divided Finland into six regional states according to the division of Regional State Administrative Agencies (Fig. 1). This division follows the natural geographical division of Finland into Southern, Southwestern, Western & Inland, Eastern, and Northern Finland as well as Lapland.
Fig. 1. Geographical division showing the regional states in Finland.
To estimate how different regions contribute to the general epidemiological pattern and to better understand the long-term cycles of PUUV infections in Finland, we constructed a first-order Poisson autoregression model using a Bayesian approach [Reference Congdon12]. We estimated the parameters under a model-averaging perspective allowing for model uncertainty. We estimated the probability of different long-term cycles being in the model and Gibbs variable selection was performed [Reference Ntzoufras13]. We included cycles from 1·5 years to 4 years with increments of 0·5 years. Cycles with a >50% inclusion probability were considered significant cycles. We used weakly informative priors. To control the correlation of errors, we made a predictive check on overdispersion.
We calculated the average monthly incidence/100 000 person-years with 95% credible intervals in each regional states using the same model. To study changes in seasonality and incidence, we divided the study into three periods: first period (1 April 1995 to 31 March 2001), second period (1 April 2001 to 31 March 2007) and third period (1 April 2007 to 31 March 2014). A negative binomial regression model was used to estimate differences in incidence rate ratios between time periods. The analyses were performed with Stata v. 13 (Stata Corp., USA) and WinBUGS v. 1.4.3 (http://www.mrc-bsu.cam.ac.uk/software/bugs/). Figures were created by Excel (Microsoft Corp., USA).
RESULTS
Descriptive data
A total of 30 942 cases of PUUV infection were identified during the study period (range by seasonal year 559–3779). Cases were mostly (61%) males (median age 48, range 3–92 years). Most cases were notified from Eastern Finland (n = 9445, 31%) and least cases from Southwestern Finland (n = 1534, 5%). The average annual incidence was 31 cases/100 000 person-years, with the highest in Eastern Finland and the lowest in Southwestern Finland (range by regional state 11·1–80·5, range by seasonal year 10·8–71·2).
Changes in incidence over-time
Overall, the model including the quadratic trend had the best fit (based on Akaike's Information Criterion) for describing the long-term trend in Finland. Compared to the first time period (reference), the average monthly incidence increased significantly over the three time periods in Southern and Eastern Finland. In Southwestern, and Western & Inland Finland incidence increased from the first to second period but somewhat decreased or remained stable in the third period (Table 1). In Northern Finland and Lapland, incidence was higher during the second period compared to the first but not significantly different in the third period (Northern Finland) and in fact lower (Lapland) from the first time period.
Table 1. Incidence rate ratios of Puumala virus infections by regional state over the three time periods, Finland, 1995–2014
Values given are incidence rate ratio (95% credible interval).
* First period (1 April 1995 to 31 March 2001); Second period (1 April 2001 to 31 March 2007); Third period (1 April 2007 to 31 March 2014).
Long-term cycles and seasonality
The model fitted well with the actual observations in all the regional states: the model fit is shown in Supplementary Figure 1 for Southern Finland. Based on the posterior probabilities of different long-term cycles, five regional states had 3-year cycles with high posterior probability while a 4-year cycle with 100% probability was observed in Lapland. The 3-year cycle probabilities for different regions were 73% (Southern), 100% (Southwestern), 86% (Western & Inland), 86% (Eastern) and 74% (Northern). Other cycle probabilities in the model were negligible except for a 22% probability for a 3·5-year cycle in Northern Finland. The peak years differed between regions. Incidence peaked in 2008 in Southern, Southwestern, Western & Inland, and Eastern Finland while the peak was observed in 2007 in Northern Finland and Lapland.
The seasonality of incidences in Southern Finland was characterized by two peaks; the first in August and the second in early winter (Fig. 2). The timing of the second peak differed between the studied periods; in the first time period, a weak peak was observed in December but during the second and third periods the incidence peaked in November. In the third period, the average monthly incidence was highest in winter, not in August as in the other periods and incidence remained high until January.
Fig. 2. Average incidence of Puumala hantavirus infection by month in the three time periods, Southern Finland, 1 April 1995–31 March 2014. The vertical bars represent 95% credible intervals.
In Southwestern Finland, a minor peak in October was observed in the first period, which was not evident in the second and third periods (Fig. 3). Incidence peaked in early winter during the third period (December) but later in winter during the second period (January); notably, similar winter peaks were not observed in the first period at all and overall, seasonality with two distinct peaks was most apparent in the third period. In Western & Inland Finland, seasonality remained similar throughout the time periods with minor peaks in August followed by major peaks in early winter. The seasonality in Eastern Finland was distinctly different from Southern, and Western & Inland Finland. Here the incidence peaked in early winter (December) with other seasonal peaks absent (data not shown).
Fig. 3. Average incidence of Puumala hantavirus infection by month in the three time periods, Southwestern Finland, 1 April 1995–31 March 2014. The vertical bars represent 95% credible intervals.
In Northern Finland, there were some differences in the seasonal pattern between the time periods: the major peak in winter was observed earlier (November) in the second period than in the first and third periods (December). Furthermore, an additional increase in average incidence was observed in January during the second period after a significant decrease in December. In Lapland, the seasonal pattern was characterized by the early winter peak in November (Fig. 4). Incidence remained high in December in the first period. In the second period, incidence increased again in January after a drop in December, analogous to Northern Finland.
Fig. 4. Average incidence of Puumala hantavirus infection by month in the three time periods, Lapland, 1 April 1995–31 March 2014. The vertical bars represent 95% credible intervals
DISCUSSION
We report that PUUV infections are not increasing uniformly across Finland, but rather, changes in incidence, seasonality and long-term cycles differ regionally. Incidence of PUUV infection, which is reportedly the highest in the world, increased from the first to second period in all regions albeit subsequently decreased or remained stable in most regions. Marked regional differences were observed in the long-term cycles and seasonality as the cycle lengths increased from south towards north and the pattern of seasonality changed over time in certain areas.
We used ‘regional state’ as the spatial unit in this study. These administrative regions, unlike the hospital district level often used in epidemiological studies in Finland [Reference Makary7, Reference Brummer-Korvenkontio9], correspond better although not perfectly with the northern, middle and southern boreal vegetation zones in which the vole populations follow similar dynamics [Reference Korpela14]. We divided the study period into three 6-year (first and second periods) or 7-year (third period) periods, with prior assumptions of cycle length, in order to better allow analyses for temporal changes. Our results on changes in incidence or seasonality are not expected to change significantly if other time division is used given the data on continuous trend. We acknowledge that the true incidence of PUUV infection is higher as only laboratory-confirmed cases seeking healthcare were included in the numerator in our register-based study.
The increase was highest in Southwestern Finland where incidence had been lowest in the country [Reference Makary7]. The bank-vole cycles may have been weaker during the first period and thus, fewer human PUUV infections were reported. Lack of awareness of PUUV infection and/or diagnostic testing by physicians in the region may also contribute to possible under-reporting. Interestingly, the amplitude of vole cycles has shown an increasing trend in this Southwestern area, at least until 2011 [Reference Korpela14], which may explain the increased incidence of human PUUV infections. Overall, the increased incidence of PUUV infections in sub-Lapland Finland from the beginning towards the end of the study period closely mirrors the observed increase in amplitude and extent of spatial synchrony in bank-vole population cycles across Finland [Reference Olsson, Leirs and Henttonen3].
The long-term cycle lengths of PUUV infections differed between the south and the north shifting from 3 to 4 years, respectively. These cycle lengths corroborate well with the spatio-temporal population dynamics of bank voles, which are to a large degree influenced by predation and climatic conditions [Reference Vaheri1, Reference Olsson, Leirs and Henttonen3, Reference Korpela14]. The pattern in Northern Finland, with a clear 3-year cycle and low but discernible support for 3·5-year cycle probably reflects the transitional position of this area between the northern and southern vegetation zones.
The regional states exhibited differing HFRS seasonality. In urban areas, particularly in Southern Finland, a clear peak in incidence occurred in August as residents from cities are exposed to the virus during their summer vacations around July in more rural areas (visiting summer cottages), as demonstrated in previous studies [Reference Makary7, Reference Clement8]. In general and in the other regions, incidence peaked in November or December when rodents looking for shelter due to the approaching winter normally enter human dwellings. It has been demonstrated previously that the increase in incidence follows vole abundance with a lag of a few months [Reference Vaheri1] although this may vary regionally. Besides the changes in bank-vole populations, aspects of human behaviour, e.g. frequency of human contact with voles and whether protective measures are used, e.g. during handling of firewood or making house repairs, both known risk factors for PUUV infection [Reference Vapalahti15–Reference Van Loock17], also play a role in the epidemiology. Furthermore, smoking, a well-established risk factor for PUUV infection, may play a role in increasing the likelihood of exposure to rodent excreta [Reference Vapalahti15, Reference Clement18]. There may be regional differences in the use of protective measures, frequency of visits and activities undertaken in summer cottages.
Temporal changes in seasonality were evident in many regions. In Southern Finland, seasonality was characterized by a clear summer peak in August in the first time period but over time, the winter peak became more apparent (with high incidence until January) and seasonality in the third period differed from the first. In Southwestern Finland, the seasonality pattern changed from the rather weak seasonal first period towards two rather distinct, albeit relatively weak peaks in the second and third periods. In Lapland, in addition to the significant decrease in incidence in the third period, incidence declined faster from the peak than during the first two periods.
These observed regional changes and differences in human incidence of PUUV infections are probably linked to a complex interaction of temporal and geographical variations in bank-vole population dynamics and climate. A recent study showed that increasing growing-season temperatures have weakened the cyclicity of voles (not specifically bank voles) in cold regions but strengthened cycles in warmer regions of Finland [Reference Korpela14]. How these findings translate to temporal and geographical changes in the incidence of human PUUV infections observed in our study warrants further modelling. The model we developed could be utilized for in-depth analyses including bank-vole populations and climatic factors with time lags in order to better understand regional differences in the epidemiology of PUUV infections. Furthermore, regional differences in human behaviour patterns should be assessed with knowledge, attitudes and practice surveys in order to effectively guide intervention measures.
SUPPLEMENTARY MATERIAL
For supplementary material accompanying this paper visit http://dx.doi.org/10.1017/S0950268816000765.
ACKNOWLEDGEMENTS
We thank Dr Otso Huitu for critically reviewing the manuscript.
DECLARATION OF INTEREST
None.
