Mathematical modelling of Echinococcus multilocularis abundance in foxes in Zurich, Switzerland
Parasites & Vectors volume 10, Article number: 21 (2017)
In Europe, the red fox (Vulpes vulpes) is the main definitive host of Echinococcus multilocularis, the aetiological agent of a severe disease in humans called alveolar echinococcosis. The distribution of this zoonotic parasite among the fox population is remarkably aggregated with few heavily infected animals harbouring much of the parasite burdens and being responsible for most of the environmental parasitic egg contamination. Important research questions explored were: (i) spatial differences in parasite infection pressure related to the level of urbanization; (ii) temporal differences in parasite infection pressure in relation to time of the year; (iii) is herd immunity or an age-dependent infection pressure responsible for the observed parasite abundance; (iv) assuming E. multilocularis infection is a clumped process, how many parasites results from a regular infection insult.
By developing and comparing different transmission models we characterised the spatio-temporal variation of the infection pressure, in terms of numbers of parasites that foxes acquired after exposure per unit time, in foxes in Zurich (Switzerland). These included the variations in infection pressure with age of fox and season and the possible regulating effect of herd immunity on parasite abundance.
The model fitting best to the observed data supported the existence of spatial and seasonal differences in infection pressure and the absence of parasite-induced host immunity. The periodic infection pressure had different amplitudes across urbanization zones with higher peaks during autumn and winter. In addition, the model indicated the existence of variations in infection pressure among age groups in foxes from the periurban zone.
These heterogeneities in infection exposure have strong implications for the implementation of targeted control interventions to lower the intensity of environmental contamination with parasite eggs and, ultimately, the infection risk to humans.
Echinococcus multilocularis is a zoonotic cestode present in large parts of the northern hemisphere. The parasite is sustained by a wildlife cycle with carnivores (mainly foxes) as definitive hosts and small mammals (mainly rodents) as intermediate hosts . However, domestic dogs are believed to be an infection source for humans in Asia [2, 3]. Humans are accidental hosts that become infected through the ingestion of parasitic eggs excreted through the faeces of the infected canids . The metacestode stage of this tapeworm causes chronic life-threatening alveolar echinococcosis (AE), which can have a high economic impact in highly endemic resource-poor settings . The geographic distribution of E. multilocularis seems to be expanding and it is considered an emerging/re-emerging pathogen in many countries [5–8]. In Europe, high prevalences (23.9–57.3%) of E. multilocularis have been frequently reported in the red fox population (Vulpes vulpes) [9–11], which is increasingly colonising urban areas [12, 13]. In Zurich (Switzerland), the abundant availability of anthropogenic food seems to have contributed to the gradual increase of the urban fox population . Moreover, the establishment of an E. multilocularis transmission cycle in the urban and periphery of Zurich is well documented [11, 15–17] as conditions appear to sustain high densities of foxes that support an active parasite life-cycle. These findings, along with the increasing incidence of human AE  have raised public health concerns and the demand to implement disease control strategies [15, 19].
Variations between individuals in their exposure and susceptibility to parasite infective stages result in the aggregated distribution of parasites within their hosts . The distribution of E. multilocularis in foxes is also characteristically aggregated with most animals carrying low numbers of parasites whereas a few harbour thousands of them. The risk of developing human AE depends, among other factors, on the amount of infective eggs found in the environment and their accessibility to humans . Due to the parasite aggregation, the degree of egg contamination in the environment depends greatly on a few highly infected animals [11, 15, 22]. However, as eggs can survive in the environment for some time, there may also be some contribution from less heavily infected foxes. Information on prevalence in foxes has been often used to characterize the infection risk for AE as its estimation is more reliable and straightforward than other epidemiological parameters . However, there is not a clear correspondence between prevalence rates and parasite abundance in the fox population . Hence, there are major limitations when using prevalence in foxes to describe the epidemiology of E. multilocularis infection . The determination of parasite abundance in animal hosts can provide valuable information to optimize parasite control strategies. For instance, if there is evidence of spatial heterogeneities in E. multilocularis infection pressure, anthelmintic baits can be distributed in areas where superinfected animals are predicted to be in order to reduce more efficiently the environmental contamination of eggs and ultimately, human infection. A key epidemiological parameter to predict parasite abundance in the animal host is the infection pressure. The parasite burden in the definitive host depends on the number of infectious stages encountered per infection insult, meaning the number of viable protoscoleces contained in the hydatid cysts that the intermediate host carries. This study complements the results reported on the force of infection by a study on the mathematical modelling of E. multilocularis infection in foxes in Zurich . There, the force of infection is defined as the number of exposures per unit time regardless of the quantity of parasites to which a fox is exposed .
The infection pressure cannot be estimated through direct observation; hence, we use mathematical models that allow inference on processes relevant to transmission as well as their quantification, in conjunction with field data. Besides the specific research question we want to address and the identification and incorporation of the epidemiological knowledge available, the selection of an appropriate model will depend on its ability to represent the available field data. The data for the present study consisted of parasite counts found in necropsied foxes collected in three different spatial zones within the municipality of Zurich . Several studies have been carried out on E. multilocularis transmission in foxes in Switzerland providing an extensive prior knowledge for model construction and hypothesis formulation. Previous studies of E. multilocularis in Switzerland have shown that transmission dynamics in animal hosts are influenced by multiple interrelated factors that contribute to its spread [11, 17, 26–28]. Decreasing parasite prevalences along with the increasing level of urbanization have been reported in foxes in the two largest cities of Switzerland [11, 27, 28]. Special attention was brought to the intermediate areas between the rural and urban habitats where the proportion of E. multilocularis coproantigen-positive fox faeces was higher compared to the urban area . These areas are believed to be heavily contaminated by infective eggs, and thus may represent hot-spots for human infection . In addition, there is evidence of seasonal variation in parasite abundance in Swiss foxes, which has been found to be related with the age of the host [11, 17, 26]. In addition, juvenile foxes of less than one-year-old have frequently been reported bearing higher infection rates and parasite burdens [11, 17, 26, 27, 29]. The study quantifying the force of infection in E. multilocularis in foxes in Zurich, defined as the number of fox exposures to parasite infection (insults) per unit time, reported spatial and seasonal variations in incidence of exposure . However, it did not address parasite abundance, which is important for the transmission dynamics. Here, we adapted existing transmission models describing the number of parasites depending on host age [30, 31] to estimate the spatio-temporal variation of the infection pressure. We aim to address further specific research questions: (i) are there spatial differences in parasite infection pressure related to the level of urbanization; (ii) are there temporal differences in parasite infection pressure in relation to time of the year; (iii) is herd immunity or an age-dependent infection pressure responsible for the observed parasite abundance; (iv) assuming E. multilocularis infection is a clumped process, how many parasites results from a regular infection insult.
The data used for this study was sourced from the necropsies of red foxes collected from January 1996 to April 2000 within the political community of Zurich as part of the Integrated Fox Project and before the implementation of an anthelmintic baiting study [11, 17]. The age of each fox was determined in years through dental examination  assuming all cubs were born on the first of April, as described previously . In this study, we used the dates when the foxes were collected to estimate the approximate age in days of each fox. Foxes less than 1-year-old were classified as juveniles . Each animal was further classified as coming from the periurban, border or urban zone, depending on where it was collected. The characteristics of each spatial zone have already been described in detail by Hegglin et al. . The periurban zone refers to the external ring delimiting the city of Zurich and which mainly comprises a green belt of forests, fields, pastures, and meadows. The border area refers to the intermediate ring that contains residential areas, allotment gardens, cemeteries, sports fields, and public parks. The urban area refers to the center of the city encompassing much of the built-up zone.
For the analysis, we used E. multilocularis intestinal counts from 531 foxes aged up to 9 years old. Thus, the parasite biomass is the total number of parasites recovered from the foxes. The median age was less than 1-year-old in all zones. The group had an overall prevalence of 41.4% and a median abundance of 0 parasites (95% central range 0–10,488 parasites). All the data is provided in Additional file 1.
Age-based abundance model
The association between parasite burden and age in foxes [11, 26, 27] suggested the use of an age-stratified SIR model originally developed by Roberts et al. . It stratifies the host population into compartments that represent their infection and immune status at a particular age and the transition between states can be described by a set of ordinary differential equations. A schematic representation of the model is given in Fig. 1.
where S is the proportion of susceptibles, t is the age of the host, γ is the rate of loss of immunity to parasites by foxes, a is the rate of acquisition of immunity, h is the infection pressure in number of parasites per year, M is the parasite abundance and μ is the parasite death rate (1/μ is the parasite life expectancy). The infection pressure, in the present report, is defined as the number of adult worms that would develop in the definitive host after parasite exposure in the absence of density-dependent constraints.
By adapting these equations, we attempt to answer our questions on E. multilocularis infection pressure and build models describing different plausible scenarios that might explain the parasite abundance observed in the foxes. As a result, a series of models assessing the existence of spatio-temporal and age-dependent variations in infection pressure were developed. The model parameters and descriptions are summarized in Tables 1 and 3, respectively.
Models assessing spatial differences in the infection pressure
The study area was divided into three spatial zones, periurban, border, and central urban covering 20%, 41%, and 39% of the study area respectively. Three different scenarios were considered: (i) the study area comprised just one spatial zone; (ii) the study area comprised two different spatial zones, the periurban and the suburban which includes the border and urban zones and; (iii) the study area comprised three different spatial zones, periurban, border and urban. The border and urban area were merged into one zone in the second scenario to consider the possibility of no spatial differences in infection pressure between both areas, as they are quite alike. These scenarios for different spatial zones were analyzed by using either 1 model with a single value for the infection pressure h, 2 different values for h for 2 zones and 3 values of h for 3 zones (in equations 1–4). Likewise, it was analyzed if there were potentially 1, 2, or 3 baseline infection pressure β 0 between the zones (see equations 3 and 4)
Models assessing time-dependent infection pressure
The evaluation of the time dependence of the infection pressure using host age as proxy for time resulted in three different functions (equation 3), where β 0 represented the baseline number of parasites in a year and β the amplitude by which this baseline could vary according to a linear or periodic relationship. The models accounted for different baseline infection pressures and amplitudes for each of the three urbanization zones: periurban (β p ), border (β b ) and urban (β u ). A log link function was implemented to ensure positive estimates of β 0 and β, as previously described .
where the infection pressure (h) at age t is given by the amplitude (β) by which the baseline number of parasites (β 0 ) varies, or does not, in a year.
Models assessing age-dependent infection pressure
The models that assume the existence of an age-dependent infection pressure include a parameter (ϕ) representing the reduction in the number of parasites that foxes acquired after exposure which is proportional to the increase of host age. Thus, the infection pressure at age t where there is both periodic infection pressure and a decrease with age is given by equation 4:
where h is the infection pressure in number of parasites per year, β 0 is the baseline number of parasites in a year, β is the amplitude by which this baseline varies periodically and ϕ is the decrease in the number of parasites related with fox age.
Echinococcus multilocularis follows a highly aggregated distribution within the animal hosts, thus, we used the negative binomial likelihood function to fit age-based abundance models to the observed data (equation 5).
where the probability of the number of parasites (s) for each sample (Z) at age (t) is given by the mean number of parasites (M) predicted by the model, where Γ represents the gamma distribution and k is the negative binomial constant of aggregation. The values for the aggregation constants for each spatial zone were estimated from the observed data using the glm.nb function from the MASS package in R . In addition, we assumed a common negative binomial constant for all age groups, as it has been previously reported . We explored variable aggregation between the zones by assigning different values of k to each zone.
The life expectancy of the parasite (1/μ) was estimated from the model fit, allowing μ to be data driven. This was compared to an estimate of μ of 8.6 from the data presented in Kapel et al. . Equations (1) and (2) including any variation in h over time, as described by equation (3), were numerically integrated using the ode function in the deSolve package in R .
Based on this probability model, a likelihood function was computed stating the probability to observe the data, given the model. The transformed, negative log-likelihood (NLL) function was minimized using the optim function of the statistical package in R . All R code is provided in Additional file 2.
For model comparison and selection we followed the method described in Rüegg et al. . The NLL of each competing model was plotted against the number of parameters. This method provided a visual aid to identify the best fitting models for each number of parameters. The selected models were then compared in pairs in increasing order of complexity, starting with the simplest model (M1). The difference of NLL between each pair of models was tested against an empirical probability distribution of the null hypothesis that the simpler model provides a better fit to the data. To this end, 500 populations were simulated from the simpler model. For each population the two competing models were fitted and the difference in NLL was computed. The NLL difference estimated from the data was then compared to the 95%-quantile of this distribution. Therefore, the more complex model would give a better fit just by chance in less than 5% of the cases (α = 0.05).
Key epidemiological parameters were quantified from the best fitting model and confidence intervals were estimated by bootstrapping. We generated 1,000 data sets and estimated the parameter values by resampling with replacement from the originally data set. That is creating 1,000 data sets, each of 531 data points being the size of the original sample of data. For the confidence interval, we reported the 2.5th and 97.5th percentiles of the bootstrap samples. For the confidence bands of the most parsimonious model, these 1,000 data sets were used to generate predicted abundances at each time point to then report their 2.5th and 97.5th percentiles.
Number of parasites per infectious insult
Using the results of the model quantifying the force of infection with the same data set , we estimated the numbers of parasites per infectious insult acquired by foxes in the periurban and urban zone at times of highest and lowest infection pressure, by using the simple equation:
Complete analysis of the data was performed using the open source software in R .
The exploratory analysis of the data showed that foxes aged up to 3 years old, which represented 86% of the total samples, accounted for 88% of all infected animals and harboured 94% of the total parasite biomass. The study data encompassed 531 observations categorised by age of the host (juveniles, n = 309; adults between 1 and 2 years, n = 99; adults between 2 and 3 years, n = 50; and adults over 3 years, n = 73), type of urbanization zone (periurban, n = 185; border, n = 200; and urban, n = 146) and season when the fox was collected (spring, n = 31; summer, n = 123; autumn, n = 113; and winter, n = 264). The seasons were defined in groups of three months: spring (March to May), summer (June to August), autumn (September to November) and winter (December to February).
The parasite counts varied widely between observations with no parasites in 59% of the foxes, 21% foxes found with 1–99 worms, 17% foxes with 100–9,999 worms, and 3% of them with more than 10,000 worms. The proportion of parasite loads found in the foxes by type of urbanization zone, fox age and season are displayed in Table 2.
The distribution of E. multilocularis in foxes was highly aggregated (overall negative binomial constant k = 0.05).
Transmission models comparing the possibility of acquired immunity or changes in infection pressure were compared to explore the hypotheses whether parasite induced immunity, seasonality, spatial differences and host age may be contributing to the observed pattern of parasite abundance in the foxes. A total of 20 models describing different scenarios for parasite transmission were compared based on their goodness-of-fit to the data and the number of parameters used, as it is illustrated in Fig. 2. The best fitting model was M20 (Table 3). Thus, models M1-M5 which described no spatial variation in transmission and models M6-M12 in which there were 2 spatial zones of transmission had a poorer fit to the data generally than models M13 to M20 where there were three spatial zones. Of these latter models those with a periodic infection pressure (M15-M20) described the data better than a non-periodic infection pressure (M13-M14). M20 and M19, with a decreasing abundance only in the periurban zone described the data better than M15 and M16 with either no decrease in infection pressure in any zone or a decrease in all 3 zones. M20 where the lower abundance in old foxes in the periurban zone is best explained by decreasing infection pressure in old fox gives a better description of the data than M19 where is hypothesizes it is due to parasite induced immunity. The difference between M17 and M20 is fixing the life expectancy of the parasite to that experimentally observed (M20) rather than using the data.
The best-fitting model, M20, assumed different parasite burdens in foxes from the periurban, border, and urban zones. The estimations of the negative binomial constants indicated variability in the degree of aggregation of the parasites (k) between the three zones (k peri = 0.1, k border = 0.02, k urban = 0.05). The model also considered age-dependent infection pressure, but only resulted in a better fit for such a model in the periurban area (i.e. old foxes had a lower exposure rate). In addition, the model suggested the existence of a sinusoidal infection pressure that varied with time with higher peaks during autumn and winter in foxes in all spatial zones, even though this seasonality was highest and most marked in the periurban zone. However, the baseline number of parasites (β 0 ) was found to be similar among the three zones and thus a single β 0 was applied to all zones. Finally, the model did not find evidence of parasite-induced immunity. Table 4 gives the maximum likelihood estimates (MLE) of the five parameters estimated by the model. Thus, the infection pressure, as described by equation (4), can be estimated at any time point (t) in any spatial zone. For example, a 10-month-old fox in the periurban zone has β 0 = 8.5, βp = 2.6 and ϕ = 0.5 (t is in years, so in this case = 0.833). Thus ln[(h(t))] = 8.5 + 2.6*sin(2*π*0.833) -0.5*0.833 = 10.3. Taking the exponent gives an infection pressure (or exposure) of 30,880 parasites per year at that time point. Likewise a 10-month-old fox in the urban zone has β 0 = 8.5 and βu = 1.2 (ϕ = 0 as in model M20) and hence an infection pressure of 13,911 parasites per year.
A graphical representation of the seasonal variation of the infection pressure on each urbanization zone can be found in Fig. 3. The model gives predictions of the infection pressure greater than zero, even for newborn foxes, due to the baseline parameter (β 0 ). However, the fox cubs are not exposed to infection during their lactation period (c.4 weeks) thus foxes less than two months were assigned no parasites to their model predictions.
The most parsimonious model (M20) therefore indicated that there were spatial variations in infection pressure, with the periurban area having the highest value of h. The infection pressure varied throughout the year in all three spatial zones, with the highest infection pressure occurring in the winter months. Variations in abundance with age that were most notable in the periurban zone were better explained by an age-related decrease in infection pressure rather than prevention of reinfection by immunity resulting from an earlier exposure.
Parasites per infectious exposure
In periurban foxes, the maximum infection pressure occurred in winter and varied between 36,000 parasites in year 1 (1st winter), 22,000 in year 2 (2nd winter), and 13,000 in year 3 (3rd winter). This results in an approximate mean of 24,000 parasites per fox over the three winters. Lewis et al. , using the same data set reported around 9.5 infectious insults per year in winter. Therefore, about 2,500 parasites result from a single infectious insult in periurban foxes during winter. Likewise, in summer periurban foxes are exposed to an average of 230 parasites per year derived from 2.3 insults or 100 parasites per insult.
In urban foxes, we predict an infection pressure during winter of approximately 15,000 parasites per year from 1.8 insults or 8,300 parasites per infection event. In summer infections there are 1,500 parasites per year from 0.5 insults or approximately 3,000 parasites per infection event.
Model 20 was the best fitting model, so the data gives support to the hypotheses that there are: (i) spatial differences in parasite infection pressure among the three zones; (ii) temporal differences in parasite infection pressure in relation to time of the year and; (iii) there are infection pressure variations across different age groups only in the periurban area. These findings are consistent with some of the often interrelated and frequently reported risk factors in EM infection in foxes . Nevertheless, some of the model implications are not in line with previous research. These findings are discussed further below in detail.
First, the model describes spatial differences in infection pressure across urbanization zones. Urban resident foxes in Zurich have been found to display small home ranges (c.25 ha) and they pursuit their daily activities mainly within this area, although some movement of foxes among urbanization zones also occurs [41, 42]. The level of urbanization of their limited territories determines the number of rodents and foxes and their predator-prey interactions, influencing ultimately parasite transmission . Therefore, the model hypothesis of an existing heterogeneous distribution of infected foxes within the city is consistent with numerous studies that found an association between infection status and type of urbanization zone . Even though this association has been often linked to other risk factors such as season [11, 16] and fox age . Most of these studies also reported a gradual decrease in parasite prevalence from the rural areas and the periphery of the cities towards the more urbanised zones [27–29]. Similarly, Lewis et al.  estimated a higher number of infection exposures in periurban foxes (maximum of 9.35 to 9.7 insults/year) compared to urban foxes (maximum of 1.6 to 2 insults/year) in Zurich. Foxes from the outside and transition areas of the cities seemed to prey more on rodents and hence, be more exposed to parasite infection [17, 29]. This comes as a result of the presence of high densities of suitable intermediate hosts bearing high parasite prevalences in the outskirt of the cities [16, 17]. In contrast, urban foxes rely more on anthropogenic food for their diet, being highly abundant and accessible in the city center . Similarly, our model estimated the highest amplitude to be in the periurban zone, implying that the highest infection pressure was borne by the periurban foxes. However, in our study this is just applicable to juvenile foxes since M20 describes the infection pressure in the periurban area is age-dependent resulting in adult foxes being exposed to less number of parasites per infection insult than their young. Consequently, only the periurban juveniles presented the highest infection pressure across all areas. In fact, the model predictions referring to adult foxes suggested the urban foxes harboured the highest infection pressure among zones. This is a surprising finding since it would be expected that animals living in the edge of the cities would prey more on rodents and thus be more exposed to infection as previously discussed. The model suggests that periurban adults are being infected more frequently on average, but that each infection event results in fewer parasites than a typical infection event in urban foxes. In absence of host immune responses to infection it may indicate that infected rodents that are consumed by urban foxes have greater numbers of protoscolices than those consumed by periurban foxes, even though urban foxes are consuming fewer infected rodents. This hypothesis differs to what has been reported in terms of parasite infection in city foxes . A potential explanation might be that some super infected foxes collected in the urban area were in fact dispersal foxes whose home range encompassed mainly the border area but they died in the urban area during an excursion looking for feeding or mating opportunities. The occurrence of so-called floating individuals with larger home ranges has been previously recorded in Zurich . These foxes are commonly young males that expand their range during the mating season (autumn and winter) . In this case, seasonal variations in the spatial behaviour of foxes could explain the higher amplitude in the infection pressure found in the urban area. Alternatively, it could be due to the short history of urban colonization of foxes, suggesting that the transmission cycle is not yet equilibrated, showing typical flickering in transiting complex systems . In any case, there is an increasing individual risk of developing AE mainly in areas where high densities of humans and urban foxes intersect [15, 24], which it is not the case of the city centre. The existence of a high infection pressure in the periphery of the cities and in the transition areas and adjacent spatial zones are still the main cause of concern in terms of AE transmission risk.
Secondly, the model also accounts for a sinusoidal infection pressure throughout the year with peaks during the cold months of autumn and winter. In Zurich, higher infection rates have been previously recorded during winter, in association with the host age and the city zone where the fox was retrieved [11, 16]. Likewise, Lewis et al.  found a periodic force of infection with an annual minimum of 0.27–1.27 parasite insults and a maximum of 6.87–7.05 parasite insults per year in foxes collected in Zurich. Evidence of seasonal variation in prevalence in foxes has been frequently reported in other locations [26, 45, 46]. In fact, seasonal fluctuations in temperature and precipitation are often proposed as infection determinants because of their influence on hosts’ numbers and parasite survival in the environment . The low temperatures and humidity favour the survival of E. multilocularis eggs in the environment, potentially contributing to the occurrence of higher infection rates during the cold and rainy seasons . In addition, the influence of climatic changes on the abundance and age-structure of vole populations [48–50] may have an impact on fox predation on voles  and consequently in the degree of parasite infection in foxes. Some studies have found that foxes showed higher predation on voles in autumn compared to spring coinciding with prey availability [11, 29]. Whereas, it has been reported a correlation between low day temperatures and higher infection rates in rodents [49, 50].
Thirdly, the model suggests the existence of an age-dependent infection pressure was found only in the periurban area. Periurban foxes prey more on rodents compare to other urbanised zones , thus if there is any effect of fox age in parasite exposure it is more likely to be evident in this area. In Zurich, higher worm burdens have been recorded in younger foxes compared to adults . In addition, seasonal variations in prevalence were more marked in the juvenile animals in the same Swiss city . Despite being a young fox has been repeatedly reported as an infection determinant for E. multilocularis , the underlying cause remains unclear. Potential reasons for decreasing parasite abundance in the periurban adult fox may include predatory behaviour or diet preferences. Juveniles might have a higher proportion of rodents in their diet, as they are abundant and easy to prey whereas adult foxes might have better access to more difficult prey or have more experience finding food from anthropogenic sources. Alternatively, inexperienced juveniles might be inclined to prey on infected voles if parasite infection adversely affects the intermediate host . Nevertheless, although variation in feeding behaviour across age groups of foxes has not been demonstrated, it has been hypothesized that juvenile foxes may have best access to voles . The dietary response of red foxes is complex when abundant alternative resources are available such as anthropogenic food and multiple intermediate host species . The diet of urban foxes has a dominance of scavenged meat and other scavenged and cultivated fruit and crops with more than half of an average stomach content being anthropogenic. The proportion of scavenged food recovered from foxes' stomachs increases in foxes found in the city center compared to the periurban area [14, 17]. Such a variation in fox dietary preference correlates with the spatial variations in infection pressure reported in the best fitting model M20.
Other proposed explanation for age-related differences in burdens of parasites considers the existence of a developing immunological response after repeated infection [11, 26, 54]. There are previous studies in highly endemic regions for E. granulosus which document a negative correlation of parasite abundance with age in dogs which would be predicted if parasite-induced host immunity limited infection [31, 55, 56]. In addition, experimental infections have shown evidence of parasite-specific antibody responses in animal hosts although it remains unclear their effect on parasite infection . Nonetheless, previous models assuming presence of acquired host immunity did not give a better fit to E. multilocularis data in dogs or foxes [25, 55]. Furthermore, studies in Kyrgyzstan and Lithuania failed to demonstrate a decrease in E. multilocularis abundance with increasing fox age [58, 59]. Moreover, in our study the parameter values on which the models incorporating immunity converged indicated a very high value of γ – the rate of loss of immunity. This would indicate that the duration of immunity following exposure would only be a matter of weeks at most, and require conditions of extremely high infection pressure to be maintained. Even if immunity were present its effect on parasite abundance would be negligible with this SIR model. Thus, the better fit to data given by the models without immunity or the potential very high rates of loss of immunity, if present, are evidence that definitive host immunity is not regulating the parasite population in this system.
Experimental studies where foxes were artificially infected have reported a pre-patent period of 29–33 days  and a patent period of up to three months . In the study of Kapel et al.  it took approximately 42 days for foxes experimentally infected with 20,000 protoscolices to reduce their worm load to 50%. We used this measure to calculate the parasite death rate in models M12 and M20 as the data itself was not able to define this parameter well (Table 3).
We have attempted to quantify the infection pressure of E. multilocularis in foxes in Zurich to gain a better insight on parasite epidemiology through hypothesis testing using a relatively simple transmission model. The modelling of the E. multilocularis infection is potentially a complex task since the dynamics of parasite transmission are influenced by a wide range of interrelated factors such as, hosts’ population densities, predator-prey interactions, landscape characteristics, climate conditions and human-related activities . Additionally, the modelling of parasite abundance brings extra challenges due to the extreme aggregation of parasites within their hosts. This intense aggregation produces a high degree of uncertainty to the model predictions. This is reflected in the wide confidence intervals related to the model predictions. Such wide confidence intervals could have been narrowed if the data set had had more data points. However, despite this the major findings of the study are robust as the most parsimonious model had an improved likelihood (or statistical fit) compared to the competing models (representing competing hypotheses). Furthermore, the basic model was first proposed by Roberts et al.  and has since been used on several data sets [31, 56, 62] and has proved robust, despite its relative simplicity. The present study introduced potential seasonal variations in infection pressure by allowing the parameter h to have a sinusoidal relationship with age (and hence time). Other data sets analyzed with this model were taken at 1 time point and thus could not analysed in such a way. Thus, in terms of the hypotheses we tested the model appears to have validity and robustness, although it would further support our findings if or when another similar data set becomes available to confirm this.
Nevertheless, models are conceived to be a simplified representation of the highly complex processes in nature and provide a useful tool to assess different hypotheses. All models are wrong, but the question remains as to how wrong the model must be to lose its usefulness . Given the assumptions in the model, the data suggests spatial and age-related variations in infection pressure to foxes. Therefore, even considering all limitations, the model offers a practical platform to improve knowledge on parasite epidemiology and to allow the quantification of epidemiological parameters that cannot be measure directly in the field, such as the infection pressure. Some of the implications derived from the model concurred with previous epidemiological knowledge on E. multilocularis infection, such as the existence of spatial heterogeneities, seasonal fluctuations and age-related differences. Alternatively, other conclusions diverged from previous reports, such as finding that the highest number of parasites developing in the fox after infection exposure occurs in the urban area. However, it cannot be ruled out the possibility that the few urban foxes found harbouring high loads of parasites might have become infected in the neighboring area previously to their incursion into the urban zone. The model also challenged the hypothesis that parasite-induced host immunity may play a role in the transmission dynamics of E. multilocularis. Using the models described we found no convincing evidence that this may be the case. The decrease in abundance in foxes, only observed in the periurban zone, is better explained by a decrease in infection pressure in older foxes, although the differences in model predictions are quite subtle.
In conclusion, the model gives a picture of the overdispersed infection pressure borne by foxes in Zurich, highlighting the potentially large contribution of young periurban foxes and foxes from the outside perimeter of urban areas towards environmental contamination. Previous studies have proved the efficacy of the use of anthelmintic baiting to reduce the environmental contamination with parasitic eggs [64–66]. Similarly, in Zurich the placement of monthly baits along the urban periphery has been shown to successfully decrease the amount of coproantigen-positive fox faeces and reduce infection rates in intermediate hosts (A. terrestris) in baits areas . However, temporal anthelmintic interventions mostly failed to achieve permanent parasite elimination [64, 68]. Hence, there is a need to ensure the optimisation of potential long-term baiting campaigns . Model results suggest that a reduction in parasite biomass in Zurich foxes could be more effectively achieved if baiting strategies were to be intensified in the periphery of the city and during the autumn and winter months. The quantification of the temporal-spatial variation of the number of parasites in foxes can help to optimise the designing of targeted bait programmes aiming to reduce the level of environmental contamination and ultimately, infection risk in humans.
maximum likelihood estimate
Eckert J, Deplazes P. Biological, epidemiological, and clinical aspects of echinococcosis, a zoonosis of increasing concern. Clin Microbiol Rev. 2004;17:107–35.
Craig P, Giraudoux P, Shi D, Bartholomot B, Barnish G, Delattre P, et al. An epidemiological and ecological study of human alveolar echinococcosis transmission in south Gansu. China Acta Trop. 2000;77:167–77.
Budke CM, Campos-Ponce M, Qian W, Torgerson PR. A canine purgation study and risk factor analysis for echinococcosis in a high endemic region of the Tibetan plateau. Vet Parasitol. 2005;127:43–9.
Torgerson PR, Keller K, Magnotta M, Ragland N. The global burden of alveolar echinococcosis. PloS Neg Trop Dis. 2010;4:722.
König A, Romig T, Thoma D, Kellermann K. Drastic increase in the prevalence of Echinococcus multilocularis in foxes (Vulpes vulpes) in southern Bavaria, Germany. Eur J Wildl Res. 2005;51:277–82.
Osterman Lind E, Juremalm M, Christensson D, Widgren S, Hallgren G, Agren EO, et al. First detection of Echinococcus multilocularis in Sweden, February to March 2011. Euro Surveill. 2011;16(14). http://www.eurosurveillance.org/ViewArticle.aspx?ArticleId=19836.
Combes B, Comte S, Raton V, Raoul F, Boué F, Umhang G, et al. Westward spread of Echinococcus multilocularis in foxes, France, 2005–2010. Emerg Infect Dis. 2012;18(12):2059–62.
Sikó SB, Deplazes P, Ceica C, Tivadar CS, Bogolin I, Popescu S, Cozma V. Echinococcus multilocularis in southeastern Europe (Romania). Parasitol Res. 2011;108:1093–7.
Bruzinskaite R, Marcinkute A, Strupas K, Sokolovas V, Deplazes P, Mathis A. Alveolar echinococcosis, Lithuania. Emerg Infect Dis. 2007;13:1618–9.
Malczewski A, Gawor J, Malczewska M. Infection of red foxes (Vulpes vulpes) with Echinococcus multilocularis during the years 2001–2004 in Poland. Parasitol Res. 2008;103:501–5.
Hofer S, Gloor S, Muller U, Mathis A, Hegglin D, Deplazes P. High prevalence of Echinococcus multilocularis in urban red foxes (Vulpes vulpes) and voles (Arvicola terrestris) in the city of Zurich, Switzerland. Parasitology. 2000;120:135–42.
Gloor S, Bontadina F, Hegglin D, Deplazes P, Breitenmoser U. The rise of urban fox populations in Switzerland. Mamm Biol. 2001;66:155–64.
Harris S, Rayner JA. Discriminant analysis of the current distribution of urban foxes (Vulpes vulpes) in Britain. J Anim Ecol. 1986;55:605–11.
Contesse P, Hegglin D, Gloor S, Bontadina F, Deplazes P. The diet of urban foxes (Vulpes vulpes) and the availability of anthropogenic food in the city of Zurich. Switzerland Mamm Biol. 2004;69:81–95.
Deplazes P, Hegglin D, Gloor S, Romig T. Wilderness in the city: the urbanization of Echinococcus multilocularis. Trends Parasitol. 2004;20:77–84.
Stieger C, Hegglin D, Schwarzenbach G, Mathis A, Deplazes P. Spatial and temporal aspects of urban transmission of Echinococcus multilocularis. Parasitology. 2002;124:631–40.
Hegglin D, Bontadina F, Contesse P, Gloor S, Deplazes P. Plasticity of predation behaviour as a putative driving force for parasite life-cycle dynamics: the case of urban foxes and Echinococcus multilocularis tapeworm. Funct Ecol. 2007;21:552–60.
Schweiger A, Ammann RW, Candinas D, Clavien P-A, Eckert J, Gottstein B, et al. Human alveolar echinococcosis after fox population increase, Switzerland. Emerg Infect Dis. 2007;13:878–82.
Hegglin D, Deplazes P. Control of Echinococcus multilocularis: Strategies, feasibility and cost–benefit analyses. Int J Parasitol. 2013;43:327–37.
Wilson K, Bjornstad ON, Dobson AP, Merler S, Poglayen G, Randolph SE, et al. Heterogeneities in macroparasite infections: patterns and processes. In: Hudson PJ, Rizzoli A, Grenfell BT, Heesterbeek H, Dobson AP, editors. The ecology of wildlife diseases. Oxford: Oxford University Press; 2002. p. 6–44.
Eckert J, Deplazes P, Kern P. Alveolar echinococcosis (Echinococcus multilocularis) and neotropical forms of echinococcosis (Echinococcus vogeli and Echinococcus oligarthrus). In: Palmer SR, Soulsby L, Torgerson PR, Brown DWG, editors. Oxford Textbook of Zoonoses Biology, Clinical Practice, and Public Health Control. Oxford: Oxford University Press; 2011. p. 669–99.
Raoul F, Deplazes P, Nonaka N, Piarroux R, Vuitton DA, Giraudoux P. Assessment of the epidemiological status of Echinococcus multilocularis in foxes in France using ELISA coprotests on fox faeces collected in the field. Int J Parasitol. 2001;31:1579–88.
Romig T. Spread of Echinococcus multilocularis in Europe? In: Craig P, Pawlowski Z, editors. Cestode Zoonoses: Echinococcosis and Cysticercosis. Amsterdam: Ios Press; 2002. p. 65–80.
König A, Romig T. Fox tapeworm Echinococcus multilocularis, an underestimated threat: a model for estimating risk of contact. Wildl Biol. 2010;16:258–66.
Lewis FI, Otero-Abad B, Hegglin D, Deplazes P, Torgerson PR. Dynamics of the force of infection: insights from Echinococcus multilocularis infection in foxes. PLoS Negl Trop Dis. 2014;8(3), e2731.
Brossard M, Andreutti C, Siegenthaler M. Infection of red foxes with Echinococcus multilocularis in western Switzerland. J Helminthol. 2007;81:369–76.
Fischer C, Reperant LA, Weber JM, Hegglin D, Deplazes P. Echinococcus multlocularis infections of rural, residential and urban foxes (Vulpes vulpes) in the canton of Geneva, Switzerland. Parasite. 2005;12:339–46.
Reperant L, Hegglin D, Fischer C, Kohler L, Weber JM, Deplazes P. Influence of urbanization on the epidemiology of intestinal helminths of the red fox (Vulpes vulpes) in Geneva, Switzerland. Parasitol Res. 2007;101:605–11.
Robardet E, Giraudoux P, Caillot C, Boue F, Cliquet F, Augot D, et al. Infection of foxes by Echinococcocus multilocularis in urban and suburban areas of Nancy, France: influence of feeding habits and environment. Parasite. 2008;15:77–85.
Roberts MG, Lawson JR, Gemmell MA. Population dynamics in echinococcosis and cysticercosis: mathematical model of the life-cycle of Echinococcus granulosus. Parasitology. 1986;92:621–41.
Torgerson PR, Shaikenov BS, Rysmukhambetova AT, Ussenbayev AE, Abdybekova AM, Burtisurnov KK. Modelling the transmission dynamics of Echinococcus granulosus in dogs in rural Kazakhstan. Parasitology. 2003;126:417–24.
Kappeler A. Die orale immunisierung von Füchsen gegen Tollwut in der Schweiz. Bern: University of Bern; 1991.
Wandeler AI. Altersbestimmung bei Füchsen. Revue Suisse Zoologie. 1976;83:956–63.
Harris S. Age determination in the red fox (Vulpes vulpes): an evaluation of technique efficiency as applied to a sample of suburban foxes. J Zool. 1978;184:91–117.
Venables WN, Ripley BD. Modern Applied Statistics with S. Fourth Edition. New York: Springer; 2002. ISBN 0-387-95457-0.
Kapel C, Torgerson PR, Thompson RCA, Deplazes P. Reproductive potential of Echinococcus multilocularis in experimentally infected foxes, dogs, raccoon dogs and cats. Int J Parasitol. 2006;36:79–86.
Soetaert K, Petzoldt T, Setzer RW. Solving differential equations in R: package deSolve. J Stat Softw. 2010;33:1–25.
Core Team R. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2015. http://www.R-project.org/.
Rüegg SR, Heinzmann D, Barbour AD, Torgerson PR. Estimation of the transmission dynamics of Theileria equi and Babesia caballi in horses. Parasitology. 2008;135:555–65.
Otero-Abad B, Torgerson PR. A systematic review of the epidemiology of echinococcosis in domestic and wild animals. PloS Neg Trop Dis. 2013;7(6), e2249.
Gloor S. The rise of urban foxes (Vulpes vulpes) in Switzerland and ecological and parasitological aspects of a population in the recently colonised city of Zürich [Dissertation]. Zurich: University of Zurich; 2002.
Baker P, Funk SM, Harris S, White PCL. Flexible spatial organization of urban foxes, Vulpes vulpes, before and during an outbreak of sarcoptic mange. Anim Behav. 2000;59:127–46.
Cavallini P. Variation in the social system of the red fox. Ethol Ecol Evol. 1996;8:323–42.
Scheffer M, Bascompte J, Brock WA, Brovkin V, Carpenter S, Dakos V, et al. Early-warning signals for critical transitions. Nature. 2009;461:53–9.
Kritsky DC, Leiby PD. Studies on sylvatic echinococcosis. V. Factors influencing prevalence of Echinococcus multilocularis Leuckart 1863, in red foxes from North Dakota, 1965-1972. J Parasitol. 1978;64:625–34.
Hanosset R, Saegerman C, Adant S, Massart L, Losson B. Echinococcus multilocularis in Belgium: prevalence in red foxes (Vulpes vulpes) and in different species of potential intermediate hosts. Vet Parasitol. 2008;151:212–7.
Veit P, Bilger B, Schad V, Schafer J, Frank W, Lucius R. Influence of environmental factors on the infectivity of Echinococcus multilocularis eggs. Parasitology. 1995;110:79–86.
Stien A, Voutilainen L, Haukisalmi V, Fuglei E, Mork T, Yoccoz NG, et al. Intestinal parasites of the Arctic fox in relation to the abundance and distribution of intermediate hosts. Parasitology. 2010;137:149–57.
Wang Q, Raoul F, Budke C, Craig PS, Xiao YF, Vuitton DA, et al. Grass height and transmission ecology of Echinococcus multilocularis in Tibetan communities, China. Chin Med J. 2010;123:61–7.
Burlet P, Deplazes P, Hegglin D. Age, season and spatio-temporal factors affecting the prevalence of Echinococcus multilocularis and Taenia taeniaeformis in Arvicola terrestris. Parasit Vectors. 2011;19:4–6.
Saitoh T, Takahashi K. The role of vole populations in prevalence of the parasite (Echinococcus multilocularis) in foxes. Res Pop Ecol. 1998;40:97–105.
Vervaeke M, Davies S, Leirs H, Verhagen R. Implications of increased susceptibility to predation for managing the sylvatic cycle of Echinococcus multilocularis. Parasitology. 2006;132:893–901.
Raoul F, Hegglin D, Giraudoux P. Trophic ecology, behaviour and host population dynamics in Echinococcus multilocularis transmission. Vet Parasitol. 2015;213:162–71.
Tackmann K, Loschner U, Mix H, Staubach C, Thulke HH, Conraths FJ. Spatial distribution patterns of Echinococcus multilocularis (Leuckart 1863) (Cestoda: Cyclophyllidea: Taeniidae) among red foxes in an endemic focus in Brandenburg, Germany. Epidemiol Infect. 1998;120:101–9.
Lahmar S, Kilani M, Torgerson PR. Frequency distributions of Echinococcus granulosus and other helminths in stray dogs in Tunisia. Ann Trop Med Parasitol. 2001;95:69–76.
Budke CM, Jiamin Q, Craig PS, Torgerson PR. Modeling the transmission of Echinococcus granulosus and Echinococcus multilocularis in dogs for a high endemic region of the Tibetan plateau. Int J Parasitol. 2005;35:163–70.
Kato N, Nonaka N, Oku Y, Kamiya M. Modified cellular immune responses in dogs infected with Echinococcus multilocularis. Parasitol Res. 2005;95:339–45.
Ziadinov I, Deplazes P, Mathis A, Mutunova B, Abdykerimov K, Nurgaziev R, Torgerson PR. Frequency distribution of Echinococcus multilocularis and other helminths of foxes in Kyrgyzstan. Vet Parasitol. 2010;171:286–92.
Bružinskaitė-Schmidhalter R, Sarkūnas M, Malakauskas A, Mathis A, Torgerson PR, Deplazes P. Helminths of red foxes (Vulpes vulpes) and raccoon dogs (Nyctereutes procyonoides) in Lithuania. Parasitology. 2012;139:120–7.
Nonaka N, Iida M, Yagi K, Ito T, Ooi HK, Oku Y, Kamiya M. Time course of coproantigen excretion in Echinococcus multilocularis infections in foxes and an alternative definitive host, golden hamsters. Int J Parisitol. 1996;26(11):1271–8.
Matsumoto J, Yagi K. Experimental studies on Echinococcus multilocularis in Japan, focusing on biohazardous stages of the parasite. Exp Parasitol. 2008;119:534–41.
Torgerson PR. Canid immunity to Echinococcus spp: impact on transmission. Parasit Immunol. 2006;21:295–303.
Box GEP, Draper NR. Empirical Model Building and Response Surfaces. New York: John Wiley & Sons; 1987.
Schelling U, Frank W, Will R, Romig T, Lucius R. Chemotherapy with praziquantel has the potential to reduce the prevalence of Echinococcus multilocularis in wild foxes (Vulpes vulpes). Ann Trop Med Parasitol. 1997;91:179–86.
Tackmann K, Löschner U, Mix H, Staubach C, Thulke HH, Ziller M, Conraths FJ. A field study to control Echinococcus multilocularis-infections of the red fox (Vulpes vulpes) in an endemic focus. Epidemiol Infect. 2001;127:577–87.
Tsukada H, Hamazaki K, Ganzorig S, Iwaki T, Konno K, Lagapa JT, et al. Potential remedy against Echinococcus multilocularis in wild red foxes using baits with anthelmintic distributed around fox breeding dens in Hokkaido, Japan. Parasitology. 2002;125:119–29.
Hegglin D, Ward PI, Deplazes P. Anthelmintic baiting of foxes against urban contamination with Echinococcus multilocularis. Emerg Infect Dis. 2003;9:1266–72.
Hegglin D, Deplazes P. Control strategy for Echinococcus multilocularis. Emerg Infect Dis. 2008;14:1626–8.
The data set has been kindly provided by the Institute of Parasitology at the Vetsuisse Faculty of Zurich. The authors would like to also acknowledge the valuable contribution of Dr. Fraser I. Lewis.
This study has been funded by the Swiss National Science Foundation (SNSF) CR3313_132482/1. Data collection carried out in the context of the Integrated Fox Project, supported by the Swiss Federal Veterinary Office, Berne (Projekt No. 1.99.03) and the Swiss Federal Office for Education and Science, Berne (EchinoRisk Projekt QLK2-CT-2001-001995/BBW No. 00.0586-2).
Availability of data and materials
The datasets are available within the article and its Additional files.
BOA and PRT conceived and designed the study. DH and PD undertook associated studies which supplied the data. BOA, SRR and PRT undertook the analysis. All authors contributed to drafting of the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
About this article
Cite this article
Otero-Abad, B., Rüegg, S.R., Hegglin, D. et al. Mathematical modelling of Echinococcus multilocularis abundance in foxes in Zurich, Switzerland. Parasites Vectors 10, 21 (2017). https://doi.org/10.1186/s13071-016-1951-1