# Landscape and Regional Environmental Analysis of the Spatial Distribution of Hantavirus Human Cases in Europe

Caroline Brigitte Zeimes, Sophie Quoilin, Heikki Henttonen, Outi

Lyytikäinen, Olli Vapalahti, Jean-Marc Reynes, Chantal Reusken, Arno N

Swart, Kirsti Vainio, Marika Hjertqvist, et al.

► To cite this version:

Caroline Brigitte Zeimes, Sophie Quoilin, Heikki Henttonen, Outi Lyytikäinen, Olli Vapalahti, et al. Landscape and Regional Environmental Analysis of the Spatial Distribution of Hantavirus Human Cases in Europe. Frontiers in Public Health, 2015, 3, 10.3389/fpubh.2015.00054. pasteur-01834449

HAL Id: pasteur-01834449

https://pasteur.hal.science/pasteur-01834449

Submitted on 10 Jul 2018

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Landscape and regional environmental analysis of the spatial distribution of hantavirus human cases in Europe

Caroline Brigitte Zeimes<fnref n="1" />*, Sophie Quoilin<fnref n="2" />, Heikki Henttonen<fnref n="3" />, Outi Lyytikäinen<fnref n="4" />, Olli Vapalahti<fnref n="5" />, Jean-Marc Reynes<fnref n="6" />, Chantal Reusken<fnref n="7" />,<fnref n="8" />, Arno N. Swart<fnref n="8" />, Kirsti Vainio<fnref n="9" />, Marika Hjertqvist<fnref n="10" /> and Sophie O. Vanwambeke<fnref n="1" />

**Edited by:**
Efstratios Vogiannis, Evangeliki Model
School of Smyrna, Greece

**Reviewed by:**
M. Jahangir Alam, University of
Houston College of Pharmacy, USA
Sofia Kottou, National and
Kapodistrian University of Athens,
Greece

**\*Correspondence:**
Caroline Brigitte Zeimes, Georges
Lemaître Centre for Earth and Climate
Research, Earth and Life Institute,
Université Catholique de Louvain
(UCL), Pace Louis Pasteur 3 bte
L4.03.07 Louvain-la-Neuve 1348,
Belgium
e-mail: caroline.zeimes@uclouvain.be

**Background:** In Europe, the most prevalent hantavirus, Puumala virus, is transmitted by bank voles and causes nephropathia epidemica in human. The European spatial distribution of nephropathia epidemica is investigated here for the first time with a rich set of environmental variables.

**Methods:** The influence of variables at the landscape and regional level is studied through multilevel logistic regression, and further information on their effects across the different European ecoregions is obtained by comparing an overall niche model (boosted regression trees) with regressions by ecoregion.

**Results:** The presence of nephropathia epidemica is likely in populated regions with well-connected forests, more intense vegetation activity, low soil water content, mild summers, and cold winters. In these regions, landscapes with a higher proportion of built-up areas in forest ecotones and lower minimum temperature in winter are expected to be more at risk. Climate and forest connectivity have a stronger effect at the regional level. If variables are staying at their current values, the models predict that nephropathia epidemica may know intensification but should not spread (although southern Sweden, the Norwegian coast, and the Netherlands should be kept under watch).

**Conclusion:** Models indicate that large-scale modeling can lead to a very high predictive power. At large scale, the effect of one variable on disease may follow three response scenarios: the effect may be the same across the entire study area, the effect can change according to the variable value, and the effect can change depending on local specificities. Each of these scenarios impacts large-scale modeling differently.

**Keywords:** hantavirus, environmental modeling, spatial model, multilevel logistic regression, large scale model

## INTRODUCTION

In Europe, Puumala virus (PUUV), the most prevalent hantavirus in Western and Northern Europe (1), is responsible for hemorrhagic fever with renal syndrome, often called nephropathia epidemica (NE) in humans (2, 3). The reservoir is the bank vole (*Myodes glareolus*) (4). Transmission to humans may be direct, by biting, but is mainly indirect, by inhalation of dust of urine and feces of infected bank voles. For example, while cleaning closed and un-aired buildings or handling firewood (5). Transmission between bank voles is also both direct and mainly indirect (6–8).

Nephropathia epidemica cases have been recorded in many European countries. A number of studies have investigated the

environmental factors influencing PUUV distribution, but never more broadly than at the national level [e.g., Ref. (9, 10)]. While the main environmental factors are fairly well known, and hypotheses explaining the broad distribution of PUUV exist (2, 11, 12), no study so far has investigated the distribution of NE cases in humans or in bank voles at the continental scale. This is attempted here for Belgium, Finland, France, the Netherlands, Norway, and Sweden with a rich set of environmental variables related to climate, land use, vegetation, soil, and human distribution. A broad range of conditions and diverse environments are encountered across Europe and various factors identified from the literature may be more relevant for some or other areas. Among climatic factors,

<footnote n="1">¹ Georges Lemaître Centre for Earth and Climate Research, Earth and Life Institute, Université Catholique de Louvain (UCL), Louvain-la-Neuve, Belgium</footnote>

<footnote n="2">² Epidemiology Unit of Infectious Diseases, Scientific Institute of Public Health, Brussels, Belgium</footnote>

<footnote n="3">³ Natural Resources Institute Finland, Vantaa, Finland</footnote>

<footnote n="4">⁴ Infectious Disease Control Unit, Department of Infectious Diseases, National Institute for Health and Welfare, Helsinki, Finland</footnote>

<footnote n="5">⁵ Infectious Disease Control Unit, Department of Virology, University of Helsinki, Helsinki, Finland</footnote>

<footnote n="6">⁶ Unité de Biologie des Infections Virales Emergentes, Centre National de Référence des Hantavirus, Institut Pasteur, Lyon, France</footnote>

<footnote n="7">⁷ Department of Viroscience, ErasmusMC, Rotterdam, Netherlands</footnote>

<footnote n="8">⁸ Centre for Infectious Disease Control, Rijksinstituut voor Volksgezondheid en Milieu (RIVM), Bilthoven, Netherlands</footnote>

<footnote n="9">⁹ Department of Virology, Norwegian Institute of Public Health, Oslo, Norway</footnote>

<footnote n="10">¹⁰ The Swedish Institute for Infectious Disease Control, Stockholm, Sweden</footnote>
higher humidity and cold temperatures favor *ex vivo* survival (8, 9). However, higher winter temperature may relate to more human outdoors activities and a higher exposure to bank vole excreta (13). In Fennoscandia, a thick snow cover shelters bank voles against predators and, in the absence of an ice sheet, preserves food for the bank voles (14–18). Snow also protects the virus from UV light and maintains the preferred stable cold temperature and high humidity for the *ex vivo* virus survival (8). In areas where snow cover is often abundant, higher winter temperatures may also lead to more contacts between human and bank voles migrating near dwellings because of the lack of protective snow cover (19, 20). In areas where the snow cover is generally low, extreme cold winter condition may also force bank voles to migrate near human dwellings (21). In areas where the snow cover is too abundant (more than 120 days/year), there is a lower prey diversity and more pressure on bank voles (14, 16, 17, 22). Land use is also important. Wet environments are preferred over drier environments (23). Forest is the main habitat of the bank vole and broadleaved forests are assumed to be the most preferred type (24). Coniferous mature and moist forests with well-developed undergrowth can also be favorable (25). Contiguous, well-connected forests may favor virus circulation in rodents (26). Buildings within 50 m of forests may lead to an increased human exposure (27). Concerning the vegetation, extended vegetation growing period increases the number of breeding opportunities for bank voles by enhancing the carrying capacity (28). Finally, soil moisture (related among other elements to soil texture) may imply a better *ex vivo* survival (8), and more humans mean higher probability of human infections.

Ecological processes take place at diverse scales, ranging from the planet, continental, regional, landscape, and local ecosystem level (29). Transitions between scales may not follow a gradual change and, in this nested structure, processes encountered at various scales differ. Looking at a broad study area means that environmental factors operating at various scales need to be considered. Our first objective is to understand the level at which variables operate and how they influence PUUV distribution. The effect of the environmental factors is tested jointly at the landscape and regional scale, using multilevel models.

Provided a regional structure is identified, investigating in further detail the regional and local effects of environmental factors are valuable in understanding the spatial distribution of NE cases across Europe. Our second objective is to understand what the effects of environmental variables are across diverse European ecoregions. This is done by using both logistic regressions on ecoregions and overall niche modeling.

# MATERIALS AND METHODS

## DATA

### Disease record collections

Data are collected for Belgium, Finland, France, the Netherlands, Norway, and Sweden. As recording protocols and time span of data availability vary across countries, the study focuses on presence or absence locations. At least 5 years of records are available for the countries included. Presence record is included in the database when spatial information is provided. In France, NE cases from 2003 to 2012 are reported by municipality of exposure (473 records) or by municipality of residence (289 records),

and 192 records had no spatial information (Centre National de Référence des Hantavirus, Institut Pasteur). In Belgium, NE cases from 2000 to 2010 are reported by municipality of residence (Institut scientifique de Santé publique-Wetenschappelijk Instituut Volksgezondheid). In the Netherlands, official notifications for 2008–2012 are provided by municipality of residence [Rijksinstituut voor Volksgezondheid en Milieu (RIVM)]. This dataset is completed with data from a cross-sectional population serological survey of 2929 persons sampled in 2006 (RIVM) and with the RIVM diagnostic database (RIVM is one of three diagnostic centers) for 2007–2011. In Norway, cases reported from 1989 to 2012 are provided by municipality of residence (Folkehelseinstituttet). In southern Sweden, human infections reported for 1997–2012 are collected by municipality of infection (Smittskyddsinstitutet). In northern Sweden, 212 human infections by coordinates of infection are recorded for 1991–1998 (Smittskyddsinstitutet). In Finland, the National Institute of health and Welfare collects the data that are generally available at hospital district level (province level). The coordinates of the administrative unit center were recorded for 2004–2009 (Metsäntutkimuslaitos). For this research, some high-incidence and low-incidence regions are selected.

The map of presences and absences is presented in **Figure 1**. While the resolution is homogeneous across the data (municipality for most countries), the average size of this unit varies greatly. In order to accommodate this, the data are converted into points locations in the following way. The centroid of municipalities between 80 and 150 km² is used as the absence/presence point. Small municipalities (area ≤80 km²) are aggregated in a 10 km resolution grid of points. Larger municipalities (area >150 km²) are represented by points located where the population density is highest (2.5 arc-min resolution, Gridded Population of the World from Center for International Earth Science Information Network). Geographical coordinates are used when available. In Finland, the presence and absence data are aggregated on a 10 km grid of points. In northern Sweden, 300 absences are randomly selected among other dwellings [see Ref. (10)]. There are a total of 1933 presences and 6425 absences.

### Environmental variables

Relevant environmental variables are identified based on a literature review of PUUV ecology in Europe, as summarized in **Table 1**. As this study focuses on the spatial distribution, seasonal or yearly temporal fluctuations are ignored and values are averaged over time. Two ecologically relevant levels are identified for the multi-level analyses: the landscape (individual) and the region (group). For operational purposes, the landscape level is defined as the area covered by a circular radius of 5 km around a presence/absence record. A sample size of at least 50 groups is recommended for a good estimation of regression coefficients in multilevel regressions (30). Fifty regions were defined using Ward clustering (R i386 3.1.0, "Rcmdr" package, "FactoMineR" plug-in) based on the mean annual temperature and the annual temperature range (World-clim, resolution of 1 km, monthly mean from 1950 to 2000). Environmental variables are calculated at the landscape level and at the regional level (ArcGIS 10.1 and FRAGSTATS, version 4).
The regional level is only used in the multilevel model. Specific hypotheses are made for the effect of environmental variable at the landscape and regional levels.

<figure><img src="image_3.png" /><figcaption>FIGURE 1 | Map of presences and absences of human hantavirus cases in Belgium, Finland, France, the Netherlands, Norway, and Sweden.</figcaption></figure>

**Climatic variables.** Climatic variables are assumed to have the same effect at both levels. Higher mean annual precipitations (Worldclim, resolution of 1 km, monthly mean from 1950 to 2000) are assumed to reflect higher air humidity and should increase the probability of presence of NE cases. The mean maximum temperature in the summer (Worldclim, 1 km, monthly mean from 1950 to 2000) is calculated to identify the driest and most unfavorable environments for virus *ex vivo* survival. Summer is defined as the months of June, July, and August. In some areas, this can reflect more human outdoor activities. The mean minimum temperature in the winter (Worldclim, 1 km, monthly mean from 1950 to 2000) is added to our database. Winter is defined as the months of December, January, and February. Lower minimum temperatures in winter are assumed to increase the presence of NE cases (thanks to better virus *ex vivo* survival) as well as the range of minimum temperature suitable for more human outdoor activities. Snow cover is here measured by the annual mean percentage of pixel area covered by snow (MODIS, 0.05°, 2000–2008 monthly mean). Less extensive snow cover (increasing human exposure) and extensive cover (favoring bank voles and virus *ex vivo* survival) should increase the probability of NE cases presence while very extensive cover (reflecting harsh condition) should decrease the probability.

**Land use variables.** We hypothesize that the proportion of forest [Corine 2006 (European Environment Agency, EEA), vector map] indicates the presence of bank voles at the landscape level, and the presence and also the abundance of bank voles at the regional level. The proportions of coniferous, broadleaved, and mixed forest [Corine 2006 (EEA), vector map] are used. It is assumed that broadleaved forest is preferred over coniferous forest but the effect of forest type may differ between the landscape and regional levels. Indeed, if a region is characterized by a large proportion of coniferous forest, the proportion of broadleaved forest could be significant at the landscape level but not at the regional level as, regionally, coniferous forests are used by default by the bank vole. The contiguity index of forest patches [Corine 2006 (EEA), vector map] is used as a proxy to represent contacts between rodents,

<table><caption>Table 1 | Variables and their hypothesized effect (X) on the bank voles' abundance, virus *ex vivo* survival, and human repartition (numbers in brackets indicate references).</caption><thead><tr><th>Variables</th><th>Bank voles</th><th>Virus</th><th>Human</th><th>Resolution</th><th>Units</th><th>Sources</th></tr></thead><tbody><tr><td>Annual precipitation (8, 23)</td><td>X</td><td>X</td><td></td><td>1 km, 1950–2000</td><td>mm</td><td>Worldclim</td></tr><tr><td>Maximum temperature in summer (8)</td><td></td><td>X</td><td></td><td>1 km, 1950–2000</td><td>°C</td><td>Worldclim</td></tr><tr><td>Minimum temperature in winter (8, 9, 13, 19)</td><td>X</td><td>X</td><td>X</td><td>1 km, 1950–2000</td><td>°C</td><td>Worldclim</td></tr><tr><td>Snow cover (8, 14–18, 22)</td><td>X</td><td>X</td><td></td><td>0.05°, 2000–2008</td><td>Area percentage</td><td>MODIS</td></tr><tr><td>Proportion of forest</td><td>X</td><td></td><td></td><td>100 m</td><td>Area percentage</td><td>Corine 2006 (EEA)</td></tr><tr><td>Proportion of coniferous forest (25)</td><td>X</td><td></td><td></td><td>100 m</td><td>Area percentage</td><td>Corine 2006 (EEA)</td></tr><tr><td>Proportion of broadleaved forest (24)</td><td>X</td><td></td><td></td><td>100 m</td><td>Area percentage</td><td>Corine 2006 (EEA)</td></tr><tr><td>Proportion of mixed forest</td><td>X</td><td></td><td></td><td>100 m</td><td>Area percentage</td><td>Corine 2006 (EEA)</td></tr><tr><td>Forest contiguity index (6, 7)</td><td>X</td><td></td><td></td><td>100 m</td><td>None</td><td>Corine 2006 (EEA)</td></tr><tr><td>Built-up areas in forest ecotones (24, 27)</td><td>X</td><td></td><td>X</td><td>100 m</td><td>Area percentage</td><td>Corine 2006 (EEA)</td></tr><tr><td>Enhanced vegetation index (EVI) (28)</td><td>X</td><td></td><td></td><td>0.0083°, 2001–2012</td><td>None</td><td>MODIS</td></tr><tr><td>Number of green days (28)</td><td>X</td><td></td><td></td><td>0.005°, 2006–2010</td><td>Number of days</td><td>MODIS</td></tr><tr><td>Soil water index (SWI) (8, 23)</td><td>X</td><td>X</td><td></td><td>25 km, 2007–2010</td><td>None</td><td>TU-WIEN</td></tr><tr><td>Population proximity index (31)</td><td></td><td></td><td>X</td><td>0.0083°, 2005</td><td>Number of persons</td><td>Environment Research Group Oxford</td></tr></tbody></table>
increasing prevalence in bank voles and so, NE cases presence. At the landscape level, a higher proportion of built-up areas in a buffer of 150 m around forests [Corine 2006 (EEA), vector map] would indicate a higher probability of contact between bank voles and humans. At the regional level, high levels of urbanization may nonetheless limit PUUV transmission.

**Vegetation indices.** The mean enhanced vegetation index (EVI) (MODIS, 0.0083°, 2001–2012 mean) and number of green days per year (MODIS, 0.005°, mean from 2006 to 2010) are included. Green days represent the number of days between the detection of the start of a growth cycle and the end of the cycle, when vegetation greenness decreases. These indices are hypothesized to reflect food availability and the number of litters (more food may reflect smaller territories and more nests) at the landscape level, as well as a longer breeding season at the regional level, and should thus increase the disease presence probability at both levels.

**Soil variable.** Higher soil water index (SWI) (TU-WIEN, 25 km), averaged over 2007–2010, should improve the virus *ex vivo* survival at the landscape level. At the regional level, we assume that a higher soil water index gives information about land use (e.g., presence of marshes and bogs) unfavorable for the bank voles.

**Human distribution variable.** The probability of human NE cases increases with human population density, whether at the landscape or regional level. To reflect the presence of humans, the mean distance weighted population proximity index (Environment Research Group Oxford, 0.0083°, 2005) is considered. The population proximity index represents the population likely to visit a place, taking into account the population in the surrounding environment (31).

# METHODS

## At which levels do the variables operate and how do they influence NE distribution?

The spatial distribution of NE cases in Europe is first modeled using multilevel logistic regression in order to identify significant environmental variables at the landscape and regional level. The intraclass correlation coefficient, a measure of the proportion of the variance found at the regional level (32, 33), is first calculated. If this value is high, a multilevel approach is recommended. A full multilevel model is subsequently fitted (32).

Multilevel regressions allow including intercepts varying between regions (random intercept model) and/or slopes varying between regions (random slopes model). Variables at the landscape and regional level are tested in a model with random intercepts and a model with random intercepts and slopes ("lme4" package in R2.12.0). These models are compared with an empty model and with each other using an ANOVA test based on the Akaike information criterion and the Bayesian information criterion. The final model includes the variables significant in regressions with one explanatory variable and the random intercepts and slopes, which have better Akaike and Bayesian information criterion. To avoid collinearity issues, variables with a variance inflation factor >10 are removed. The area under the curve (AUC) is used to evaluate the predictive power of the model (34). Resulting probabilities of

being in presence of a case and false presences and absences are mapped.

## What are the effects of environmental variables across European ecoregions?

In order to better distinguish the effects of environmental variables across Europe, logistic regressions are fitted for each ecoregion using landscape level variables. Seven ecological regions encountered in the study area are compared (European Topic Centre on Nature Protection and Biodiversity) (Figure 2).

The effect of environmental factors is not necessarily linear, and may be positive in some regions and negative in some other. Because of the broad range of values covered, this is likely to be encountered in our study area. This type of relationship is poorly addressed by classical regression, and a niche model is thus needed to explore the effect of the landscape level variables across the study area.

Boosted regression trees (BRT), which allow to model non-linear response, are used to build an European model including all variables ("gbm" package in R) (35). The AUC is here based on cross-validation on 10 subsets, as overfitting is a known issue with BRT. The relative importance of a variable, a weighted measure of the number of times a variable is used to build the consecutive trees, is calculated. Graphs of the probability of presence of the disease as a function of the value of the variable are produced. Only the global trends should be considered, as local peaks may result from interactions with other variables (e.g., strong interactions between the built-up areas in forest ecotones and the maximum temperature in summer are observed). Boxplots representing the

<figure><img src="image_4.png" /><figcaption>FIGURE 2 | Digital map of European ecological regions for the study area.</figcaption></figure>
distribution of the variable in each ecological region, with colors corresponding to the map in **Figure 2**, are superimposed to this curve in order to jointly examine the global effect as identified by BRT and its subset in the variable range of each ecoregion. Finally, the signs of the significant coefficients (at the level of 0.05) of bivariate logistic regressions by ecoregion are added in front of each boxplot for comparison with the trend modeled by BRT in this ecoregion. An example graph is presented in **Figure 3**.

# RESULTS

## AT WHICH LEVELS DO THE VARIABLES OPERATE AND HOW DO THEY INFLUENCE THE DISTRIBUTION?

The intraclass correlation coefficient shows that 57% of the variance is at the regional level. A multilevel model is therefore advised.

Following the variance inflation factor analysis, snow cover, green days, and the proportion of forest are removed from further analyses to avoid collinearity issues.

In bivariate analyses, all variables at the landscape and regional levels are highly significant, except the maximum temperature in summer and the EVI at the landscape level.

The results of the multilevel logistic regression are presented in **Table 2**. The AUC is $0.895 \pm 0.004$ and indicates a good predictive power.

Random slopes are included for the landscape level variables coniferous forest, broadleaved forest, mixed forest, contiguity of forest, built-up areas in forest ecotones, and population proximity index and, at the regional level, minimum temperature in winter and EVI.

<figure><img src="image_5.png" /><figcaption>**FIGURE 3** | LEFT – Black curve: relative probability of presence of hantavirus human cases according to the evolution of the variable as modeled by BRT. Colored boxplots: distribution of the variable by ecoregions (colors refer to **Figure 2** and the height of the boxplot is proportional to the number of points by ecoregion). + and – signs: sign</figcaption></figure>

of the logistic regression coefficient, if significant, for this variable in this ecoregion. Dashed color lines: boxplot extent (second to third quartile) of an ecoregion on the black curve. RIGHT – Boxplots of the presence and absence distributions of coniferous forest for two ecoregions.

<table><caption>**Table 2** | Multilevel logistic regression (significant at the level of *0.05, **0.01, and ***0.001).</caption><thead><tr><th rowspan="2"></th><th colspan="3">Landscape level</th><th colspan="3">Regional level</th></tr><tr><th>Estimator</th><th>P-value</th><th>Random slope</th><th>Estimator</th><th>P-value</th><th>Random slope</th></tr></thead><tbody><tr><td>Annual precipitations</td><td>0.00</td><td>2.23E-01</td><td></td><td></td><td></td><td></td></tr><tr><td>Maximum temperature in summer</td><td></td><td></td><td></td><td>-0.94**</td><td>2.34E-03</td><td></td></tr><tr><td>Minimum temperature in winter</td><td>-0.15*</td><td>3.34E-02</td><td></td><td>-1.62*</td><td>2.69E-02</td><td>Yes</td></tr><tr><td>Coniferous forest</td><td>-0.00</td><td>2.74E-01</td><td>Yes</td><td>-0.10</td><td>1.75E-01</td><td></td></tr><tr><td>Broadleaved forest</td><td>0.02</td><td>5.99E-02</td><td>Yes</td><td></td><td></td><td></td></tr><tr><td>Mixed forest</td><td>-0.01</td><td>6.97E-01</td><td>Yes</td><td>-0.39</td><td>9.87E-02</td><td></td></tr><tr><td>Contiguity of forest</td><td>0.01</td><td>9.61E-01</td><td>Yes</td><td>37.66*</td><td>1.26E-02</td><td></td></tr><tr><td>Built-up areas in forest ecotones</td><td>0.28***</td><td>1.41E-15</td><td>Yes</td><td></td><td></td><td></td></tr><tr><td>EVI</td><td></td><td></td><td></td><td>0.04*</td><td>2.46E-02</td><td>Yes</td></tr><tr><td>SWI</td><td>0.00</td><td>2.07E-01</td><td></td><td>-0.14*</td><td>2.01E-02</td><td></td></tr><tr><td>Population proximity index</td><td>0.01</td><td>1.42E-01</td><td>Yes</td><td>0.11*</td><td>4.44E-02</td><td></td></tr></tbody></table>
At the landscape level, the coefficient of the variable built-up areas in forest ecotones is significantly positive at the level of 0.001. The coefficient of the minimum temperature in winter is significantly negative at the level of 0.05. At the regional level, the coefficients of forest contiguity, EVI, and population proximity index are significantly positive at the level of 0.05. Also, the maximum temperature in summer is negatively significant at the level of 0.01, and the minimum temperature in winter and SWI at the level of 0.05.

The predicted probabilities and the false presences/absences of the model are mapped in **Figure 4**. The model predicts the highest concentrations of probabilities in northeastern France, southern Belgium, inland Norway, central Sweden, and southern and middle Finland. False presences are mainly located in areas where the diseases is already established and along the Norwegian coast. A few false presences are also predicted in less expected areas: southern and central France and southern Sweden. False absences are more present at the western and southern border of the French disease distribution, in northern Belgium, in the southern Netherlands, in southern and northern Sweden, and in northern Finland.

## WHAT ARE THE EFFECTS OF ENVIRONMENTAL VARIABLES ACROSS EUROPEAN ECOREGIONS?

The multilevel regression indicates that the effect of environmental factors varies across regions and random slopes were included for several variables. Logistic regressions are therefore elaborated for each ecoregion (Figure 2). When significant, forest, mixed forest, forest contiguity, and built-up areas in forest ecotones are

positive; EVI, green days, and population proximity index are positive or negative; annual precipitation is positive except for the Scandinavian montane birch forest and grassland; maximum temperature in summer is negative except for western European broadleaf and sarmatic mixed forest; minimum temperature in winter is negative except for Fennoscandian and Russian taiga; snow cover and the SWI are positive except for Fennoscandian and Russian taiga; coniferous forest is negative except for the northern temperate Atlantic and the Scandinavian montane birch; and broadleaved forest is positive except for northern temperate Atlantic and Fennoscandian and Russian taiga.

The AUC of the boosted regression trees model is 0.921 ± 0.003. The AUC on cross-validation is 0.893 ± 0.003. These indicate a very good predictive power of the model. The variables with a relative importance over 10 are the minimum temperature in winter, the snow cover, and the maximum temperature in summer.

The response curves (Figure 5) of snow cover, broadleaved forest, mixed forest, contiguity of forest, built-up areas in forest ecotones, SWI, and population proximity index show an overall positive trend. The curves of the minimum temperature in winter, forest, coniferous forest, and EVI present an overall negative trend. Green days show a mixed curve. Annual precipitation and maximum temperature in the summer display an increasing then decreasing trend in probabilities.

The boxplots in Figure 5 indicate that variables with a similar width of boxes between regions are the proportion of mixed forest, the built-up areas in forest ecotones, and the population proximity index. Variables with boxplot extents following a gradient from

<figure><img src="image_6.png" /><figcaption>FIGURE 4 | Map of predicted probabilities and false presences/absences of the European multilevel model (threshold for presence = 0.25, when the sensitivity equals the specificity at 84%).</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>FIGURE 5 | Response curves of the variable according to the predicted probabilities of the boosted regression trees, boxplots of the variable per ecoregion, and significant signs of the coefficient of bivariate logistic regressions per ecoregion.</figcaption></figure>
South to North are minimum temperature in winter, snow cover, contiguity of forest, and the SWI. The ecoregional boxplot of EVI and green days shows a strong difference between western European and Fennoscandian ecoregions. Finally, the boxplot extent seems mixed between regions for the annual precipitations, the maximum temperature in summer, forest, coniferous forest, and broadleaved forest.

The general trend by ecoregions can be compared to the sign of the coefficient from the bivariate logistic regressions (Figure 5). The signs of the regression coefficient and the BRT trends are consistent for minimum temperature in winter, snow cover, mixed forest, contiguity of the forest, built-up areas in forest ecotones, green days, and SWI.

# DISCUSSION

## LARGE-SCALE MODELING

Both the multilevel and the BRT models have very good predictive powers. The distribution of zoonoses at a large-scale using a single database built from national databases is successfully modeled. However, large-scale models can obscure the local effect of environmental factors or how they may differ from place to place. With these models, the interpretation of the variables may generate misleading conclusion as wide range of conditions cannot be summarized into one explanation. We identify three scenarios of response to a variable that can be encountered in modeling studies at a large-scale (Figure 6). These scenarios are illustrated by variables from this study, but are relevant to other disease modeling studies.

### Scenario A

The first easiest scenario implies that a variable has the same effect over the complete study area. These variables can be easily included in a large-scale model. In our study, the built-up areas in forest ecotones are assumed to increase the human exposure. This is reflected both by the ecoregional regressions and the global BRT model.

### Scenario B

The second scenario, the most common, is when the variable has a non-linear effect on the probability of the disease. Because of this, the effect of the variable varies from region to region, depending on the value of the variable encountered. For these variables, results and hypotheses from local studies can be used for places with the

same variable extent. They can be used in non-linear large-scale models. For example, looking at the snow cover BRT trend, less extensive snow cover (such as in southern Europe) increases the probability as it relates to more contact between the virus/bank voles and humans; a more extensive snow cover increases the probability even more as that implies a larger bank voles population and a better *ex vivo* virus survival; very extensive snow cover (such as in the extreme north of Europe, in Fennoscandian, and Russian taiga) has a decreasing probability because of the harsh conditions for the bank voles and the increasing pressure due to lower prey diversity. The sign of the ecoregional regression matches these hypotheses as well as the BRT response trend.

### Scenario C

The third scenario, the most difficult to identify, is when the variable does not seem to reply to one hypothesis, whatever the extent of its distributions. Different hypotheses are needed for each region. For these variables, results from local studies cannot be applied and interpolated to other places and their interpretation in large-scale study may lead to misleading conclusion. Taking the coniferous forest as an example, the global BRT trend is negative as it is generally expected, but the sign of the regression coefficient appears significantly positive for the northern temperate Atlantic and the Scandinavian montane birch ecoregions. In the northern temperate Atlantic ecoregion, there are very few forests and it is possible that bank voles are found in coniferous forests, lacking a better habitat. For the Scandinavian montane birch ecoregions, coniferous forests are in valleys while mountains are colonized by birches (36). Soils on the slopes are drier and know important runoff. In those areas, the conditions in broadleaved forest are less favorable for the virus and rodents than in the coniferous forest. This is probably why coniferous forests appear to be positively related to NE presence in these regions. The effect of coniferous is thus not the result of one non-linear process but can be explained by the local particularities of ecoregions. These particularities may result from the level of the other variables or from historical and cultural factors.

## MODEL VARIABLES

So far, the distribution of NE cases has not been studied at the European level, and this study helps to answer whether factors

<figure><img src="image_8.png" /><figcaption>FIGURE 6 | Three scenarios that can happen when modeling variables over a large study area, squared frames represent different subregions of a large study area.</figcaption></figure>
impact the distribution in one region or at a multi-country level [see question raised in Ref. (11)]. Our model highlights the new finding that climatic variables as well as the connectivity of forests act more at the regional level than at the landscape level. The multilevel logistic model shows that NE cases are mostly found in populated regions with connected forests, high EVI, and low SWI. These regions have a low maximum temperature in summer and a low minimum temperature in winter (summers are mild and winters are cold). In these regions, landscapes with a higher proportion of built-up areas in forest ecotones and a lower minimum temperature in winter are expected to be more at risk. Models results are consistent with what has been already hypothesized on NE in Europe (2, 11). Unexpected curves are nonetheless observed for EVI and green days but their relative importance is low, and their effect complex to interpret (37).

By combining ecoregions and a global BRT model, we identified variables for which it is advised to pay attention when hypothesizing effects (Scenario C): precipitations, coniferous forest, broadleaved forests, and population. For the minimum temperature in winter, the snow cover, green days, the SWI, the proportion of mixed forest, the built-up areas in forest ecotones, and the contiguity of forests, our results suggest that they follow one process that can be non-linear.

## MODEL PREDICTIONS

The predicted probabilities map shows that, according to environmental data, fewer cases should be expected in southern, western, central, and the extreme east of France, northern Belgium, the Netherlands, southern and northern Sweden, and northern Finland. The BRT model yields similar predictions. Also, a lot of false presences are found in areas where the disease is already recorded and false absences are located at the edge of these areas. If variables stay at their current values (this model predicts the current situation and not the future), false presences can be interpreted as suitable places to find the disease and where the disease may appear in the future, while false absences can represent unsuitable environments. It seems that the recorded data already covered the potential distribution of NE cases. A spread of the disease should not occur but intensification is possible.

## LIMITATION OF THE MODELS

As we learned here, effects can vary between scales and, even if the landscape scale is interesting, a finer scale may bring further details and reflect more particularities. Exact geographical coordinates are rarely known or disclosed, due to privacy rules but also due to the nature of the hantavirus pathogenesis, and the size of the administrative units may be too big to make proper models. NE cases are not reported in the same way everywhere (between and within countries). Cases can, for example, be identified from antibodies from serological surveys or by the identification of clinical symptoms. A bias may probably result from the difference of these reports. Inclusion of a variable accounting for these differences could solve this issue but is particularly challenging as distinction of the way of records is not always known. Cases are most often reported by residence or where the disease is identified, rather than the true infection location. Therefore, some presences are reported in unexpected places, like in southern Sweden where

PUUV is reportedly absent. These presences were nevertheless included in the database and were all predicted as absences. However, in this area, some false presences are predicted near these false absences. It could indicate that the recorded presences do indicate local transmission, unless these were indeed false predictions. The models presented here are not designed to predict the future distribution of PUUV. False presences indicate places where the confirmation of the absence or presence of cases must be assessed (e.g., through serosurvey among exposed humans). In southern and central France, absence of confirmed case will probably result from the fact that places are too far and isolated from the current distribution to imagine that the disease will reach these places. The other false absences places like southern Sweden, the Netherlands, and the Norwegian coast, even if located in lower probability areas, should be carefully considered in the future for monitoring and surveillance.

## CONCLUSION

This study successfully models the pan-European NE distribution using a compilation of national disease data source and European level environmental databases. Europe covers diverse landscapes and ecoregions and, even if some variables are recurring in the literature for various areas, the underlying hypotheses can be different. This study highlights that these differences can be the result of local specificities or of the non-linearity of the processes. When the effect of a factor varies locally, its inclusion in a large-scale model is compromised.

According to both multilevel logistic and non-linear model, the distribution of PUUV seems limited in space by environmental variables that are identified here. Our model does not predict an important spread of the disease, assuming current conditions, but possibly intensification in places where the disease is already there. Also, some lower probabilities areas such as the South of Sweden, the Netherlands, and the Norwegian inner coast should be further studied.

## AUTHOR CONTRIBUTIONS

CZ and SV contribute to the conception and design of the work, the analysis and interpretation of data, and the drafting of the work. SQ, HH, OL, OV, J-MR, CR, AS, KV, and MH contribute to the design of the work, the acquisition of the data, the critical revision for important intellectual content. All the authors approve the version to be published and agree to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.

## ACKNOWLEDGMENTS

This study was funded by EU grant FP7-261504 EDENext and is cataloged by the EDENext Steering Committee as EDENext310 (http://www.edenext.eu). The contents of this publication are the sole responsibility of the authors and do not necessarily reflect the views of the European Commission.

## SUPPLEMENTAL MATERIALS

The Supplementary Material for this article can be found online at http://journal.frontiersin.org/article/10.3389/fpubh.2015.00054/abstract
One Excel file with three sheets:

1. Table with results for each ecoregion of logistic regression in bivariate, complete, and step models.

2. Table with results for each country of logistic regression in bivariate, complete, and step models.

3. Table of random intercepts and slopes added in the multilevel model.

## REFERENCES

1. Clement J, Lameire N, Keyaerts E, Maes P, Van Ranst M. Hantavirus infections in Europe. *Lancet Infect Dis* (2003) **3**:752–3. doi:10.1016/S1473-3099(03)00827-2

2. Olsson G, Leirs H, Henttonen H. Hantaviruses and their hosts in Europe: reservoirs here and there, but not everywhere? *Vector Borne Zoonotic Dis* (2010) **10**:549–61. doi:10.1089/vbz.2009.0138

3. Heyman P, Ceianu C, Christova I, Tordo N, Beersma M, Alves M, et al. A five-year perspective on the situation of haemorrhagic fever with renal syndrome and status of the hantavirus reservoirs in Europe, 2005–2010. *Euro Surveill* (2011) **16**:15–22.

4. Leduc J. Epidemiology of hemorrhagic-fever viruses. *Rev Infect Dis* (1989) **11**:S730–5. doi:10.1093/clinids/11.Supplement_4.S730

5. Dearing M, Dizney L. Ecology of hantavirus in a changing world. *Ann NY Acad Sci* (2010) **1195**:99–112. doi:10.1111/j.1749-6632.2010.05452.x

6. Vanapeldoorn R, Oostenbrink W, Vanwinden A, Vanderzee F. Effects of habitat fragmentation on the bank vole, *Clethrionomys glareolus*, in an agricultural landscape. *Oikos* (1992) **65**:265–74. doi:10.2307/3545018

7. Kozakiewicz M, Van Apeldoorn R, Bergers P, Gortat T, Kozakiewicz A. Landscape approach to bank vole ecology. *Pol J Ecol* (2000) **48**:149–61.

8. Kallio E, Klingstrom J, Gustafsson E, Manni T, Vaheri A, Henttonen H, et al. Prolonged survival of Puumala hantavirus outside the host: evidence for indirect transmission via the environment. *J Gen Virol* (2006) **87**:2127–34. doi:10.1099/vir.0.81643-0

9. Linard C, Tersago K, Leirs H, Lambin EF. Environmental conditions and Puumala virus transmission in Belgium. *Int J Health Geogr* (2007) **6**:55. doi:10.1186/1476-072X-6-55

10. Zeimes CB, Olsson GE, Ahlm C, Vanwambeke SO. Modelling zoonotic diseases in humans: comparison of methods for hantavirus in Sweden. *Int J Health Geogr* (2012) **11**:39. doi:10.1186/1476-072X-11-39

11. Heyman P, Thoma BR, Marié JL, Cochez C, Essbauer SS. In search for factors that drive hantavirus epidemics. *Front Physiol* (2012) **3**:237. doi:10.3389/fphys.2012.00237

12. Vaheri A, Henttonen H, Voutilainen L, Mustonen J, Sironen T, Vapalahti O. Hantavirus infections in Europe and their impact on public health. *Rev Med Virol* (2012) **23**:35–49. doi:10.1002/rmv.1722

13. Piechotowski I, Brockmann S, Schwarz C, Winter C, Ranft U, Pfaff G. Emergence of hantavirus in south Germany: rodents, climate and human infections. *Parasitol Res* (2008) **103**:S131–7. doi:10.1007/s00436-008-1055-8

14. Hansson L, Henttonen H. Gradients in density variations of small rodents – the importance of latitude and snow cover. *Oecologia* (1985) **67**:394–402. doi:10.1007/bf00384946

15. Hansson L, Henttonen H. Rodent dynamics as community processES. *Trends Ecol Evol* (1988) **3**:195–200. doi:10.1016/0169-5347(88)90006-7

16. Hanski I, Hansson L, Henttonen H. Specialist predators, generalist predators, and the microtine rodent cycle. *J Anim Ecol* (1991) **60**:353–67. doi:10.2307/5465

17. Hanski I, Henttonen H, Korpimaki E, Oksanen L, Turchin P. Small-rodent dynamics and predation. *Ecology* (2001) **82**:1505–20. doi:10.2307/2679796

18. Korslund L, Steen H. Small rodent winter survival: snow conditions limit access to food resources. *J Anim Ecol* (2006) **75**:156–66. doi:10.1111/j.1365-2656.2005.01031.x

19. Evander M, Ahlm C. Milder winters in northern Scandinavia may contribute to larger outbreaks of haemorrhagic fever virus. *Glob Health Action* (2009) **2**:2020. doi:10.3402/gha.v2io.2020

20. Schwarz AC, Ranft U, Piechotowski I, Childs JE, Brockmann SO. Risk factors for human infection with Puumala virus, southwestern Germany. *Emerg Infect Dis* (2009) **15**:1032–9. doi:10.3201/eid1507.081413

21. Faber MS, Ulrich RG, Frank C, Brockmann SO, Pfaff GM, Jacob J, et al. Steep rise in notified hantavirus infections in Germany, April 2010. *Euro Surveill* (2010) **15**:19574.

22. Hansson L. Intraspecific variation in dynamics: small rodents between food and predation in changing landscapes. *Oikos* (1999) **86**:159–69. doi:10.2307/3546581

23. Verhagen R, Leirs H, Tkachenko E, Van Der Groen G. Ecological and epidemiological data on hantavirus in bank vole populations in Belgium. *Arch Virol* (1986) **91**:193–205. doi:10.1007/BF01314280

24. Linard C, Lamarque P, Heyman P, Ducoffre G, Luyasu V, Tersago K, et al. Determinants of the geographic distribution of Puumala virus and Lyme borreliosis infections in Belgium. *Int J Health Geogr* (2007) **6**:20. doi:10.1186/1476-072X-6-15

25. Olsson G, White N, Hjalten J, Ahlm C. Habitat factors associated with bank voles (*Clethrionomys glareolus*) and concomitant hantavirus in northern Sweden. *Vector Borne Zoonotic Dis* (2005) **5**:315–23. doi:10.1089/vbz.2005.5.315

26. Guivier E, Galan M, Chaval Y, Xuéreb A, Ribas Salvador A, Poulle ML, et al. Landscape genetics highlights the role of bank vole metapopulation dynamics in the epidemiology of Puumala hantavirus. *Mol Ecol* (2011) **20**:3569–83. doi:10.1111/j.1365-294X.2011.05199.x

27. Crowcroft NS, Infuso A, Lief D, Le Guenno B, Desenclos JC, Van Loock F, et al. Risk factors for human hantavirus infection: Franco-Belgian collaborative case-control study during 1995-6 epidemic. *Br Med J* (1999) **318**:1737–8. doi:10.1136/bmj.318.7200.1737

28. Barrios JM, Verstraeten WW, Maes P, Clement J, Aerts JM, Haredasht SA, et al. Satellite derived forest phenology and its relation with nephropathia epidemica in Belgium. *Int J Environ Res Public Health* (2010) **7**:2486–500. doi:10.3390/ijerph7062486

29. Forman R. *Land Mosaics: The Ecology of Landscapes and Regions*. Cambridge, IN: Cambridge University Press (1995).

30. Maas CJM, Hox JJ. Sufficient sample sizes for multilevel modeling. *Methodology* (2005) **1**:86–92. doi:10.1027/1614-2241.1.3.86

31. Alexander NS, Wint W. Projected population proximity indices (30km) for 2005, 2030 & 2050. *Open Health Data* (2013) **1:e2**. doi:10.5334/jophd.abd

32. Snijders TAB, Boskers RJ. *Multilevel Analysis: An Introduction to Basic and Advanced Multilevel Modeling*. London: SAGE Publications Ltd (1999).

33. Merlo J, Chaix B, Yang M, Lynch J, Rástam L. A brief conceptual tutorial of multilevel analysis in social epidemiology: linking the statistical concept of clustering to the idea of contextual phenomenon. *J Epidemiol Community Health* (2005) **59**:443–9. doi:10.1136/jech.2004.023929

34. Fawcett T. An introduction to ROC analysis. *Pattern Recognit Lett* (2006) **27**:861–74. doi:10.1016/j.patrec.2005.10.010

35. Elith J, Leathwick J, Hastie T. A working guide to boosted regression trees. *J Anim Ecol* (2008) **77**:802–13. doi:10.1111/j.1365-2656.2008.01390.x

36. Bryn A, Durojeanni P, Hemsing LØ, O'donnell S. A high-resolution GIS null model of potential forest expansion following land use changes in Norway. *Scand J For Res* (2013) **28**:81–98. doi:10.1080/02827581.2012.689005

37. Lhermitte S, Verbesselt J, Verstraeten WW, Coppin P. A comparison of time series similarity measures for classification and change detection of ecosystem dynamics. *Remote Sens Environ* (2011) **115**:3129–52. doi:10.1016/j.rse.2011.06.020

**Conflict of Interest Statement:** The research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Received: 15 January 2015; paper pending published: 10 February 2015; accepted: 17 March 2015; published online: 31 March 2015.

Citation: Zeimes CB, Quoilin S, Henttonen H, Lyytikäinen O, Vapalahti O, Reynes J-M, Reusken C, Swart AN, Vainio K, Hjertqvist M and Vanwambeke SO (2015) *Landscape and regional environmental analysis of the spatial distribution of hantavirus human cases in Europe*. *Front. Public Health* **3**:54. doi: 10.3389/fpubh.2015.00054
This article was submitted to Environmental Health, a section of the journal Frontiers in Public Health.

Copyright © 2015 Zeimes, Quoilin, Henttonen, Lyytikäinen, Vapalahti, Reynes, Reusken, Swart, Vainio, Hjertqvist and Vanwambeke. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.