Skip directly to site content Skip directly to page options Skip directly to A-Z link Skip directly to A-Z link Skip directly to A-Z link
Volume 23, Number 8—August 2017

Density-Dependent Prevalence of Francisella tularensis in Fluctuating Vole Populations, Northwestern Spain

Article Metrics
citations of this article
EID Journal Metrics on Scopus
Ruth Rodríguez-Pastor, Raquel Escudero, Dolors Vidal, François Mougeot, Beatriz Arroyo, Xavier Lambin, Ave Maria Vila-Coro, Isabel Rodríguez-Moreno, Pedro Anda, and Juan J. Luque-LarenaComments to Author 
Author affiliations: Universidad de Valladolid, Palencia, Spain (R. Rodríguez-Pastor, J.J. Luque-Larena); Instituto Universitario de Investigación en Gestión Forestal Sostenible, Palencia (R. Rodríguez-Pastor, J.J. Luque-Larena); Instituto de Salud Carlos III, Madrid, Spain (R. Escudero, A.M. Vila-Coro, I. Rodríguez-Moreno, P. Anda); Universidad de Castilla-La Mancha, Ciudad Real, Spain (D. Vidal); Instituto de Investigación en Recorsosos Cinégeticos del Consejo Superior de Investigaciones Científicas, Ciudad Real (F. Mougeot, B. Arroyo); University of Aberdeen, Aberdeen, Scotland, UK (X. Lambin)

Cite This Article


Tularemia in humans in northwestern Spain is associated with increases in vole populations. Prevalence of infection with Francisella tularensis in common voles increased to 33% during a vole population fluctuation. This finding confirms that voles are spillover agents for zoonotic outbreaks. Ecologic interactions associated with tularemia prevention should be considered.

Emerging infectious zoonotic diseases are increasing worldwide, and most zoonoses are linked to wildlife (1,2). Thus, quantifying disease prevalence in potential wildlife hosts is critical to understanding the outbreak dynamics of zoonoses (3). Tularemia, which is caused by Francisella tularensis, is a problematic zoonotic disease worldwide, but its ecology remains poorly understood. This pathogen is classified by the US Centers for Disease Control and Prevention as a class A bio-threat agent (F. tularensis subsp. holarctica in Europe) because only a few bacteria are needed to induce tularemia in humans or susceptible animal species (>250 hosts described) (4). However, the relative epidemiologic roles (i.e., reservoir, spillover, and amplification agents) for different hosts are uncertain.

A major hotspot for tularemia in Europe is northwestern Spain (Castilla and León region), where the largest recent outbreaks of the disease have been recorded (>1,000 officially confirmed human cases during 1997–1998 and 2007–2008) (5). In intensive farmlands in Europe, rodents and lagomorphs are the main putative mammalian hosts (5,6), but most studies addressing the epidemiologic roles of these species have been correlative or used opportunistic sampling.

A recent study suggests that common voles (Microtus arvalis) are a key agent for human tularemia in northwestern Spain because of a spatial and temporal coincidence between human tularemia cases and increases in number of voles (5). Voles periodically fluctuate in density and can reach high numbers during specific periods in farming areas (5). Dead voles infected with F. tularensis subsp. holarctica have been reported in northwestern Spain during massive decreases in vole populations (7). If, as hypothesized (5), common voles are a key amplifying and spillover agent for tularemia in intensive farming areas in northwestern Spain, we should expect an increased prevalence of tularemia in voles as their numbers increase. Thus, it is crucial to empirically evaluate whether such a density-dependent pattern occurs in natural populations.

We obtained samples from live voles periodically collected during population increases (2013–2015) in northwestern Spain. Our goal was to determine how prevalence of F. tularensis in common voles varies with population density.

The Study

We complied with all necessary licenses and permits for conducting this study. During 2013–2015, a common vole population fluctuation, which peaked in 2014, was observed in agricultural areas of Castilla and León, Spain (8). This increase in vole numbers was moderate (in terms of peak density) compared with previous increases when tularemia outbreaks among humans were reported (1997–1998 and 2007–2008) by the National Network of Epidemiologic Surveillance of Spain (5,8). In 2014, no outbreak of tularemia was reported. However, there was a higher-than-average number of identified cases of tularemia (n = 95) among humans in the area (the regional average is 3 [range 0–11] cases/y, excluding outbreak years) (5).

To monitor vole abundance during the complete population fluctuation, we sampled 80 km2 of farmland in Palencia Province, Spain (42°1′N, 4°42′W), where human tularemia cases have been reported (5,8). We live-trapped voles seasonally (every 4 months) during March 2013–March 2015. Our vole trapping effort was constant (840 traps set for 24 h/seasonal sampling), and our sampling design was spatially stratified (we obtained random samples from 8 alfalfas fields, 8 grain fields, and 8 fallow fields at each seasonal sampling). Vole abundance was estimated as the number of captures/100 traps/24 h in each season. Trapping was extractive, and animals were brought alive to our laboratory in rodent cages provided with food, water, and bedding immediately after their capture. Voles were euthanized by using CO2. Carcasses were individually frozen at −30°C.

We extracted DNA from a homogenized mixture of liver and spleen (≈25 mg). DNA was extracted by using standard procedures (QIAamp DNA Mini Kit; QIAGEN, Valencia, CA, USA). A phylogenetically informative region of the lipoprotein A (lpnA) gene (231 bp) was amplified by conventional PCR and hybridized with specific probes by reverse-line blotting as described (9). We tested positive samples by using a real-time multitarget TaqMan PCR and tul4 and ISFtu2 assays (10). Negative controls for PCR and DNA extraction were included in each group of samples processed. We used R 3.2.4 software ( for statistical analyses.


Thumbnail of Vole abundance and tularemia prevalence, northwestern Spain. A) Temporal variations in vole abundance (no. captures/100 traps/24 h) and tularemia prevalence. Spain. Four voles were tested in March 2013, 15 in July 2013, 32 in November 2013, 63 in March 2014, 102 in July 2014, 19 in November 2014, and 8 in March 2015. B) Relationship between tularemia prevalence and vole abundance. Histograms show number of positive (top) or negative (bottom) voles sampled at each level of vole densi

Figure. Vole abundance and tularemia prevalence, northwestern Spain. A) Temporal variations in vole abundance (no. captures/100 traps/24 h) and tularemia prevalence. Four voles were tested in March 2013, 15 in July 2013,...

We tested 243 live voles and found an average prevalence of F. tularensis of 20.16%. Prevalence greatly varied between samplings (range 0%–33%) and was strongly related to vole abundance (generalized linear model, χ2 = 21.64, df = 1, p<0.001) with a direct and positive density-dependent association (Figure). The predicted odds of tularemia infection increased by 1.037 (95% CI 1.021–1.056) when vole density increased by +1 captured vole/100 traps/24 h (range during the study 1–60 voles/100 traps/24 h). During the vole population peak in July 2014, a total of 34 (33%) of 102 sampled live voles were infected with F. tularensis.


We report a direct and positive density-dependent association between prevalence of F. tularensis in common voles and their abundance in agricultural landscapes. These findings are consistent with vole-to-vole transmission and amplification of the bacterium as vole density increases. Voles experimentally infected with F. tularensis die within a few days after rapid acute infection and generally show high bacterial loads in organs (11). Thus, transmission between voles might involve direct contact, cannibalism, or contamination of the environment. In our study, all voles tested were alive and free of obvious signs of disease when captured, which implied that prevalence could be higher than what we estimated if moribund voles were less trappable and underrepresented in trapped animals.

The role that exogenous sources (i.e., other animals, environmental sources) might play in modulating infection prevalence among vole populations still needs to be clarified. Notwithstanding and irrespective of the precise mechanism(s) of transmission, our results support the hypothesis that exponential growth of common vole populations is crucial for amplification of tularemia transmission in farmlands, and that increases in vole populations are linked to periodic emergence of human cases of tularemia in Spain (5). Vole density can reach >1,000 voles/hectare (i.e., >300 tularemia-infected voles/hectare) during outbreaks, potentially leading to contamination of the environment and other wildlife, including harvestable species, such as crayfish and hares, which have higher contact rates with humans than voles (5).

Tularemia is probably not completely enzootic in vole populations because we did not detect F. tularensis at low densities of voles, which suggests involvement of animal or environmental reservoirs. A key unknown facet of the ecologic cycle of F. tularensis is where does it persist between epizootic periods (5,6). There is no evidence for F. tularensis replication in arthropods. However, ticks might be a reservoir of this pathogen because they have life-long infections. Thus, mammalian populations are probably needed to amplify tularemia in the environment (11).

Characteristic spatial and social behaviors of voles, including increased contact rates, aggression, and wounding, during massive population increases readily account for amplification of disease transmission rates and spread (5,12). Although reservoir and vector hosts of F. tularensis at variable densities can play major roles in the ecologic cycle of tularemia in different ecosystems, there appears to be a common link between tularemia outbreaks and rodent population fluctuations across Europe (5,6,11,13).

Common voles are useful for surveillance of tularemia, and strategic prevention programs should incorporate their temporal fluctuations in planned preventive actions. Because vole numbers seem to modulate the risk for disease exposure in humans, monitoring vole population dynamics can help anticipate and increase awareness of the risk for tularemia in rural areas of Spain.

Ms. Rodríguez-Pastor is a biologist and doctoral candidate in the Department of Agroforestry Sciences, University of Valladolid, Palencia, Spain. Her primary research interests are ecology of zoonoses and disease dynamics in animal populations.



J.J.L.L., F.M., and R.R.P. held official licenses for trapping wildlife in Spain. Capture permits were provided by the Dirección General del Medio Natural, Junta de Castilla y León.

This study was supported by projects ECOVOLE (grant CGL2012-35348), ECOTULA (grant CGL2015-66962-C2-1-R), and RESERTULA (grant CLG2015-66962-C2-2-R), which were funded by the Ministerio de Economía y Competitividad MINECO/FEDER, Spain. R.R.P. was supported by a PhD studentship from the University of Valladolid (co-funded by Banco Santander).



  1. Karesh  WB, Dobson  A, Lloyd-Smith  JO, Lubroth  J, Dixon  MA, Bennett  M, et al. Ecology of zoonoses: natural and unnatural histories. Lancet. 2012;380:193645. DOIPubMedGoogle Scholar
  2. Han  BA, Schmidt  JP, Bowden  SE, Drake  JM. Rodent reservoirs of future zoonotic diseases. Proc Natl Acad Sci U S A. 2015;112:703944. DOIPubMedGoogle Scholar
  3. Walton  L, Marion  G, Davidson  RS, White  PCL, Smith  LA, Gavier-Widen  D, et al. The ecology of wildlife disease surveillance: demographic and prevalence fluctuations undermine surveillance. J Appl Ecol. 2016;53:14609. DOIGoogle Scholar
  4. World Health Organization. Guidelines on tularaemia. Geneva: The Organization; 2007 [cited 2017 May 3].
  5. Luque-Larena  JJ, Mougeot  F, Roig  DV, Lambin  X, Rodríguez-Pastor  R, Rodríguez-Valín  E, et al. Tularemia outbreaks and common vole (Microtus arvalis) irruptive population dynamics in northwestern Spain, 1997–2014. Vector Borne Zoonotic Dis. 2015;15:56870. DOIPubMedGoogle Scholar
  6. Gyuranecz  M, Reiczigel  J, Krisztalovics  K, Monse  L, Szabóné  GK, Szilágyi  A, et al. Factors influencing emergence of tularemia, Hungary, 1984-2010. Emerg Infect Dis. 2012;18:137981. DOIPubMedGoogle Scholar
  7. Vidal  D, Alzaga  V, Luque-Larena  JJ, Mateo  R, Arroyo  L, Viñuela  J. Possible interaction between a rodenticide treatment and a pathogen in common vole (Microtus arvalis) during a population peak. Sci Total Environ. 2009;408:26771. DOIPubMedGoogle Scholar
  8. Luque-Larena  JJ, Mougeot  F, Viñuela  J, Jareño  D, Arroyo  L, Lambin  X, et al. Recent large-scale range expansion and outbreaks of the common vole (Microtus arvalis) in NW Spain. Basic Appl Ecol. 2013;14:43241. DOIGoogle Scholar
  9. Escudero  R, Toledo  A, Gil  H, Kovácsová  K, Rodríguez-Vargas  M, Jado  I, et al. Molecular method for discrimination between Francisella tularensis and Francisella-like endosymbionts. J Clin Microbiol. 2008;46:313943. DOIPubMedGoogle Scholar
  10. Versage  JL, Severin  DD, Chu  MC, Petersen  JM. Development of a multitarget real-time TaqMan PCR assay for enhanced detection of Francisella tularensis in complex specimens. J Clin Microbiol. 2003;41:54929. DOIPubMedGoogle Scholar
  11. Rossow  H, Forbes  KM, Tarkka  E, Kinnunen  PM, Hemmilä  H, Huitu  O, et al. Experimental Infection of voles with Francisella tularensis indicates their amplification role in tularemia outbreaks. PLoS One. 2014;9:e108864. DOIPubMedGoogle Scholar
  12. Ostfeld  RS, Mills  JN. Social behavior, demography, and rodent-borne pathogens. In: Wolff JO, Sherman PW, editors. Rodent societies: an ecological and evolutionary perspective. Chicago: University of Chicago Press; 2007. p. 478–86.
  13. Rossow  H, Ollgren  J, Hytönen  J, Rissanen  H, Huitu  O, Henttonen  H, et al. Incidence and seroprevalence of tularaemia in Finland, 1995 to 2013: regional epidemics with cyclic pattern. Euro Surveill. 2015;20:21209. DOIPubMedGoogle Scholar




Cite This Article

DOI: 10.3201/eid2308.161194

Table of Contents – Volume 23, Number 8—August 2017

EID Search Options
presentation_01 Advanced Article Search – Search articles by author and/or keyword.
presentation_01 Articles by Country Search – Search articles by the topic country.
presentation_01 Article Type Search – Search articles by article type and issue.



Please use the form below to submit correspondence to the authors or contact them at the following address:

Juan J. Luque-Larena, Departamento de Ciencias Agroforestales, Escuela Técnica Superior de Ingenierías Agrarias, Universidad de Valladolid, Avenida de Madrid 44, Palencia, E-34004, Spain

Send To

10000 character(s) remaining.


Page created: July 17, 2017
Page updated: July 17, 2017
Page reviewed: July 17, 2017
The conclusions, findings, and opinions expressed by authors contributing to this journal do not necessarily reflect the official position of the U.S. Department of Health and Human Services, the Public Health Service, the Centers for Disease Control and Prevention, or the authors' affiliated institutions. Use of trade names is for identification only and does not imply endorsement by any of the groups named above.