http://orcid.org/0000-0003-0376-325X
Figures
Abstract
Estimating the tick-borne encephalitis (TBE) infection risk under substantial uncertainties of the vector abundance, environmental condition and human-tick interaction is important for evidence-informed public health intervention strategies. Estimating this risk is computationally challenging since the data we observe, i.e., the human incidence of TBE, is only the final outcome of the tick-host transmission and tick-human contact processes. The challenge also increases since the complex TBE virus (TBEV) transmission cycle involves the non-systemic route of transmission between co-feeding ticks. Here, we describe the hidden Markov transition process, using a novel TBEV transmission-human case reporting cascade model that couples the susceptible-infected compartmental model describing the TBEV transmission dynamics among ticks, animal hosts and humans, with the stochastic observation process of human TBE reporting given infection. By fitting human incidence data in Hungary to the transmission model, we estimate key parameters relevant to the tick-host interaction and tick-human transmission. We then use the parametrized cascade model to assess the transmission potential of TBEV in the enzootic cycle with respect to the climate change, and to evaluate the contribution of non-systemic transmission. We show that the TBEV transmission potential in the enzootic cycle has been increasing along with the increased temperature though the TBE human incidence has dropped since 1990s, emphasizing the importance of persistent public health interventions. By demonstrating that non-systemic transmission pathway is a significant factor in the transmission of TBEV in Hungary, we conclude that the risk of TBE infection will be highly underestimated if the non-systemic transmission route is neglected in the risk assessment.
Citation: Nah K, Magpantay FMG, Bede-Fazekas Á, Röst G, Trájer AJ, Wu X, et al. (2019) Assessing systemic and non-systemic transmission risk of tick-borne encephalitis virus in Hungary. PLoS ONE 14(6): e0217206. https://doi.org/10.1371/journal.pone.0217206
Editor: Sam R. Telford III, Tufts University Cummings School of Veterinary Medicine, UNITED STATES
Received: October 23, 2018; Accepted: May 7, 2019; Published: June 4, 2019
Copyright: © 2019 Nah et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the manuscript and its Supporting Information files.
Funding: This research is part of the project entitled ‘`Estimating risk of tick-borne encephalitis with changes in climate, habitation and recreational activities" which received funding from GlaxoSmithKline SA (HO-17-18014 (K.N., X.Z., J.W.) https://www.gsk.com. GlaxoSmithKline Biologicals SA was provided the opportunity to review a preliminary version of this manuscript for factual accuracy but the authors are solely responsible for final content and interpretation. The project was also supported by Nemzetgazdasági Minisztérium (English: Ministry for National Economy, Hungary) (GINOP-2.3.2-15-2016-00019 [Á.B., A.T.] http://www.kormany.hu/en/ministry-for-national-economy), National Natural Science Foundation of China (11701072 [X.Z.]), 11501358 (X.W.) http://www.nsfc.gov.cn/english/site_1/index.html and Natural Sciences and Engineering Research Council of Canada (105588-2011 [J.W.]) http://www.nserc-crsng.gc.ca/. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: This research is part of the project entitled ‘`Estimating risk of tick-borne encephalitis with changes in climate, habitation and recreational activities" which received funding from GlaxoSmithKline SA (HO-17-18014 [K.N., X.Z., J.W.] https://www.gsk.com). GlaxoSmithKline Biologicals SA was provided the opportunity to review a preliminary version of this manuscript for factual accuracy but the authors are solely responsible for final content and interpretation. This does not alter our adherence to PLOS ONE policies on sharing data and materials.
Introduction
Tick-borne encephalitis (TBE), an arboviral infection of the central nervous system caused by bites from infected ticks, has been a major public health concern for much of Central Europe. In certain parts of Central Europe (e.g. in the Czech Republic) the occurrence of TBE is increasing since climate change pushes the vector to the higher altitudes [1, 2]. The risk of human infection depends on the prevalence of the virus spread in the tick-host ecological cycle, and on the tick-human contact that is effective for the disease transmission. Evaluating the human infection needs accurate estimation of the tick population dynamics which is highly regulated by the environmental and climatic conditions [3–8]; the effectiveness of pathogen transmission per tick-host contact, and the human case reporting given infection probability which is quite stochastic [9]. All of these involve substantial uncertainties in ecological, epidemiological and environmental parameters.
In the affected countries of the European Union, TBE incidence data is available and indicates the disease trend. For example, in Hungary TBE cases have been reported to National Database of Epidemiological Surveillance System since 1977. Between 1977 to 1996, the average annual incidence was approximately 2.7 per 100,000, and this incidence has shown dramatic decrease since 1997. Multiple factors may have contributed to this decrease, including under-reporting followed by decreased serological examination, or public vaccination in early 1990s [10, 11]. Whether or not the decrease in the incidence reflects the actual decrease in human TBE infection level, it seems probable that TBE infection level has not been decreased in the ecological cycle. The studies on the seropositivity of animals in 1960 to 1970 and 2005 support this argument [12, 13].
One of our goals is to develop a mathematical model to infer the TBEV transmission level in the ecological cycle. In particular, we aim to use a mathematical model fitting the surveillance data to determine whether there was significant increase in the basic reproduction number of TBE, and whether there was significant decrease in the expected TBE cases between 1980 to 2015.
There may be substantial amount of unreported cases, since many TBE infections go undiagnosed or unreported [10, 14]. Therefore, estimation of the TBE infection risk under significant uncertainties of the tick-host interaction and tick-human transmission process is challenging since the data we observe, i.e., human incidence of TBE, is only the final outcome of the tick-host transmission, tick-human contact and the case reporting upon infection. To meet this challenge, here we develop a coupled TBEV transmission-human report cascade model, that consists of two important components: the TBEV transmission dynamics among ticks, animal hosts and humans, and the stochastic process of TBE reporting given human TBE infection. We will use this coupled cascade model to facilitate fitting the human TBE incidence data to the model in generating good estimation of key parameters including seasonal human-tick encounter rates and human case reporting probability.
TBEV transmission involve two important routes, systemic and non-systemic transmission. Transovarial transmission is possible, but it only occurs at a low frequency [15]. Systemic transmission of TBEV involves infected ticks and vertebrate hosts. The tick Ixodes ricinus is the vector of the European TBEV subtype [16]. The Ixodes ticks undergo complex developmental cycle involving egg, larva, nymph and adult stages and the full life cycle takes average of 3 years [17]. Ticks take blood meals from hosts to develop from one stage to the next stage. At each blood meal, ticks are integrated into the epidemiological chain of the virus transmission. Systemic transmission of TBEV has the following cycle: hosts acquire infection from infected nymphs. The infected hosts pass the virus to feeding larvae. As the infected larvae develop into nymphs, the nymphs can again transmit virus into a new host. Several mathematical models, see for example [16, 18–21], have been developed and analyzed to examine the ecological or epidemiological factors that govern the abundance of Ixodes ricinus ticks or TBE infections. However, only few studies have focused on the significance of the non-systemic route of TBEV transmission [22–24], through which a susceptible vector can acquire the infection by co-feeding with infected vectors on the same host [25–29] even when the pathogen has not established within the host for systemic transmission. It was observed that TBE in Europe may be mainly maintained by non-systemic transmission between co-feeding ticks [30], and an early modeling study [23] shows that the basic reproduction number (R0) of TBE without non-systemic route of transmission is estimated to be less than 1, which means that TBE would not persist without non-systemic route of transmission. In comparison, the work of [31] shows that systemic transmission cycle alone can sustain the transmission of TBEV in a natural focus. Both systemic and non-systemic transmissions, as well as the vector abundance, are heavily influenced by the climatic and environmental conditions, which are characterized by uncertainty and seasonal variations. Some recently developed stage-structured tick population and tick-borne diseases (Lyme disease in particular) models [20, 32, 33] have used these environmental condition data. However, because of uncertainties in some of the parameters in the aforementioned models, and in the consideration of both systemic and non-systematic transmission with the co-feeding transmission efficacy unsettled, we will need not only refine existing model structure but also tune and estimate some of the model parameters by fitting a disease transmission model to time series incidence data [34, 35].
Data fitting is an important issue in TBE modeling, since human is only a dead-end host of TBEV transmission [36], and can be infected by bites of both infected nymphs and female ticks, although rarely, consumption of unpasteurized milk from infected animals can also cause TBE infection [37]. In addition, two-thirds of human TBEV infections are either sub-clinical or asymptomatic [38], suggesting that there may be more infected cases than cases reported, therefore the probability of reporting given infection needs to be estimated in order to accurately evaluate the TBE risk in the population.
In our developed TBEV transmission-human report cascade model, we couple the TBEV (systemic and non-systemic) transmission dynamics among ticks and animal hosts with the stochastic process of TBE reporting given human TBE infection. We will use this coupled system and the TBE incidence data, together with climate data to estimate such key epidemiological parameters as seasonal human-tick encounter rate, case reporting probability, tick-host contact for disease transmission, and tick-human transmission. With these key parameters and probabilities estimated, we can use our coupled system to estimate the basic reproduction numbers of TBE in Hungary, evaluate the transmission potential of TBEV in the enzootic cycle along with the climate change, and assess the significance of non-systemic pathway in the transmission of TBEV in Hungary.
Materials and methods
We introduce a tick-host and tick-human transmission cascade model with two integrated parts to describe the dynamical relationship between the epidemic in the tick-host population and the subsequent tick-to-human infections. It is a hidden Markov Model with a transition process characterized by a deterministic TBEV (tick-host) transmission model among tick and host populations and an observation process described by a stochastic model of human TBE reporting given infection.
TBEV transmission model among ticks and hosts
Our TBEV transmission model describes the seasonal transmission of a pathogen among ticks and hosts with a system of ordinary differential equations with periodic coefficients, modified from some earlier work of modeling Lyme disease dynamics [32, 33] but adapted for some unique features of TBEV dynamics including in particular co-feeding transmission. We will first stratify the tick population and host population into susceptibles and infecteds, depending on the state of infection; we further stratify the tick population by their physiological stages.
Stages and development.
Ticks in questing stages move to engorged stages upon successful host-feeding. Nymphal and larval Ixodes ricinus feed on small-to-medium-sized mammals such as rodents, while adult ticks feed on large-sized mammals [39, 40]. Engorged adults are reproducible and lay eggs while engorged larvae and nymphs move into the next stages after completing the maturation process. As these reproduction and development processes involve a few key biological activities which are highly relevant to the host density and environment conditions, we need to further stratify the tick population within each physiological stage. Namely, we need to consider the following activity states of tick population: questing larvae (Lq), engorged larvae (Le), questing nymphs (Nq), engorged nymphs (Ne), questing adults (Aq), engorged adults (Ae) and eggs (E).
Eggs develop into the questing larval stage with the developmental rate del(t). Questing larvae attach to hosts with rate αl(t). Among them, only the proportion (fl) who survive the feeding stage move to the engorged larval stage. Once the engorged larvae completes maturation they are accounted as questing nymphs. The same process is repeated from questing nymphs to engorged nymphs and from questing adults to engorged adults. Then, the female engorged adults lay eggs with oviposition rate dpop(t). A parameter η refers to the proportion of engorged female adults to engorged adults. We use Ricker function to describe the birth rate of eggs, with the parameter p being the maximal egg-laying rate and the parameter ω representing the degree of density dependence in fecundity [41]. Each stages have distinct mortalities (μe, μql, μel, μqn, μen, μqa, μea).
Seasonal dependence.
Developmental rates and the host-attaching rates have strong seasonal dependence and those parameters are thus time-dependent and functions of the independent time variable t. Other parameters showing less seasonal dependence will be assumed to be constant. Since the tick questing activity depends on a climatic condition [42], we decompose the host-attaching rates (αl(t), αn(t) and αa(t)) into the proportion of actively searching ticks among questing ticks at time t (pl(t), pn(t) and pa(t)) (and here and in what follows, subindices l, n and a are indicating these parameters are relevant to the particular stages of larvae, nymphs and adults respectively) and the host-finding rate of those actively searching ticks (λl, λn and λa). That is,
Tick population growth (ecological) model.
We can now formulate a system of ordinary differential equations for the tick population growth (ecological) model:
(1)
while the list of variables appear in Table 1, and the full list of the parameters appears in Table 2. The number of egg-laying adults is implicitly accounted in the model, see S1 Appendix for the details.
Transmission diagram and systemic transmission.
Fig 1 illustrates some major epidemiological processes that must be incorporated into our transmission dynamics model. Transovarial infection (transmission from infected females to their eggs) for TBEV occurs only rarely, so we will not consider this vertical transmission route, and hence there is no infection by questing larvae stage [43]. Since we ignore the vertical transmission, there is no need to divide the engorged adult population into susceptible and infecteds. In the systemic transmission route, questing ticks acquire viruses by feeding infected hosts, and infected ticks can also pass the virus to susceptible hosts during feeding. Note that the hosts for adult ticks can also be infected and pass the virus via untreated dairy products, however, TBE is mainly transmitted by tick bite and thus not incorporated explicitly in our model system [44]. Questing adults are further stratified into susceptibles and infecteds since questing adults involve in the human infection. In what follows, the standard incidence is used to describe the systemic transmission of TBEV through the contact of a susceptible host (tick) with an infected tick (host).
Dashed lines depict infection routes. Tick population is stratified by the stages: questing larvae (Lq), engorged larvae (Le), questing nymphs (Nq), engorged nymphs (Ne), questing adults (Aq), engorged adults (Ae) and eggs (E). The additional subscripts s and i refer to susceptible and infecteds, respectively. Host population is divided into susceptible (Hs) and infecteds (Hi). The variables are described in Table 1.
Co-feeding transmission.
Our main focus is on how to model the transmission of TBEV between co-feeding ticks [26, 45]. Let δ(Nqi(t), H) be the probability of a susceptible feeding tick being infected by co-feeding nymphs during which the ticks are co-feeding a host. Let c be the probability of an infected nymph to induce non-systemic infection to the co-feeding susceptible tick. We assume that the events of each infected nymph triggering non-systemic infection to co-feeding larvae are independent. Then, the probability of a larvae to get non-systemic infection while it is co-feeding with n number of infected nymphs is
(2)
The number of feeding infected nymphs at time t is
where Tf is the average duration of feeding and
is the mortality of feeding infected nymphs. Considering that the duration of feeding is relatively short, we estimate the average number of feeding infected nymphs per host at time t with
From the above assumptions, we have the following non-linear co-feeding probability formulation:
Non-systemic transmission can also happen when the host is immune to the infection [46]. Therefore, the force of infection for the questing larvae through non-systemic transmission route is
This formulation captures the fact that non-systemic transmission is only possible when both infective ticks and susceptible ticks are actively questing [7].
By the above discussions, the tick population growth model (1) can be generalized to the following TBEV transmission dynamics among the relevant hosts and ticks:
(3)
The list of variables appear in Table 1, and the full list of the parameters appears in Table 2.
Because of the limited data on the host population, we remove explicit dependence on the parameter H from the model system by normalizing other variables with H. The normalized system is given in S1 Appendix.
TBE human infection and reporting model.
Most I. ricinus tick bites to human are from ticks in nymphal and adult stages [47]. In our model, we assume that human can be infected with TBEV by the bites of infected nymphs and infected adult ticks. The number of newly infected cases between time t − Δ and t is given by
(4)
where Λn(t) and Λa(t) are the temperature-dependent human-attaching rate of nymphs and adults,
adapting the temperature related human outdoor activity in Hungary [48]. Here, the parameter α representing a degree of temperature dependency on the human-attaching rate of nymphs will be estimated from the surveillance data, and κ (also a parameter to be estimated) is the relative ratio between nymphs and adults for the attachments to the human.
We denote this number of newly infected cases between time t − Δ and t by It, with Δ being the reporting period. We assume that Ct, the number of cases reported during reported period Δ follows a negative binomial distribution with mean ρIt and variance , where ρ is the reporting probability and τ stands for the overdispersion parameter. Again, these parameters will be estimated using the parameter identification procedures described below, when other parameters, listed in Tables 3 and 4 are used. The transmission efficacy of non-systemic transmission from a single infected nymph (c) is estimated from the equation 1 − (1 − c)2 = 0.65 in (2), where 0.65 is the probability of transmission of TBEV when co-feeding with two infected ticks [46]. Indeed, it is empirically observed that non-systemic transmission of TBEV in non-viraemic host is less efficient than transmission in viraemic host [25], as we see from our estimates, c < βhl. According to [6], the questing nymph activity is observed when the daily maximal temperature was over the 7°C. We assume that the questing nymph is inactive in a temperature below 7°C.
Data fitting and results
We use weekly mean temperature data between 1901 and 2015 obtained from a weather station in Szombathely (coordinates: 47.20N, 16.65E, 201.0m) in Vas County [54]. Vas region is one of the highly endemic area in Hungary [55]. Weekly human TBE incidence data from 1998 to 2008 is obtained from National Epidemiological Center of Hungary. The clinically diagnosed cases confirmed by laboratory serological test are counted as the reported TBE case.
With given parameter values in Tables 3 and 4 and temperature data, we fit human reporting model (4) with TBE incidence data between 1998 and 2008. The climate data between 1901-1998 was used for the initialization of the model, as explained in detail in S1 Appendix. By the maximum likelihood estimation, we estimate the unknown parameters: probabilities of successful feeding (fl, fn, fa), host-attaching rate of actively questing ticks (λl, λn, λa), degree of density dependent fecundity (ω), host recovery rate (γ), mortality of hosts (b), relative ratio between nymphs and adults for the human attachments (κ), the degree of temperature dependency on the human-attaching rate of nymphs (α), the minimum temperature for the activity of questing larvae (ml) and the reporting probability (ρ). The maximum likelihood estimation was performed by the R package pomp, using trajectory matching [56, 57]. Using Latin hypercube sampling, we choose 104 number of initial set of parameters and compared the likelihoods at each parameters set. We repeat the sampling process at the parameter sets with the maximum likelihoods. The convergence is checked by computing the likelihood profiles (S1 Appendix) over each of the unknown parameters.
Fig 2 shows a model fit and the 1998–2008 TBE case reports in Hungary. The sample simulation of the case report model at the maximal likelihood estimation depicts the binomial curve, which is known to be the characteristic of the TBE incidence in Hungary as well as many other European countries [58]. With the set of parameter values listed in Table 5, we obtain the parameterized model for TBEV transmission. All these values are in the 95% interval that yields the highest likelihood values. The model is then validated by the comparison of the model estimators with the data obtained in the observational studies.
Red line shows the monthly accumulation of the reported case data and skyblue line shows the sample simulation of the case report model at the maximal likelihood estimation. Fixed parameter values and estimated parameter values used for simulation are given in Tables 3 and 5.
Prevalence of TBE in small mammals
According to the recent field study of [59, 60], 3.7% to 20.5% of rodents collected in Hungary during 2010-2013 were seropositive. In comparison, the annual average of infected rodents calculated from our parameterized model is 18.6% of the total rodents. For this estimation, we used the annual cycles of temperature data in Szombathely from 2010 to 2013 [54].
Reporting probability of human TBE infecteds
More than two-thirds of TBE cases are known to be asymptomatic [61]. In the observation of [62], 30–50% of the studied tick-borne disease cases remembered a tick bite prior to the disease onset. For these reasons, the reported case of TBE is expected to be lower than the actual number of infected case. In comparison, the reporting probability (ρ) through our data fitting is estimated to be 0.64.
Stage-dependent frequency of attachment to human
Estimated ratio of human-attaching rate for nymphs to the human-attaching rate for adults (κ) is less than 1, indicating shorter human questing time for adult ticks to nymphs.
Duration of host infectivity
Estimated recovery rate of hosts is 0.1 leading to the average duration of the host infectivity as 10 days, while the experimentally-observed duration of host infectivity to ticks is known to be 2 to 3 days [26].
Probability of surviving the feeding stage
From our model, the probability of successful feeding for host-attached larvae (fl) and nymphs (fn) are estimated to be 0.19 and 0.96, respectively. According to the experimental study on the infestation of Ixodes ricinus, 90 to 99 percent of nymphs and 46 to 72 percent of larvae has successfully fed the mice in the laboratory setting [53].
Temperature dependence of the questing activity for larvae
In our model fitting, minimum temperature for the coincidence of host availability and the activity of questing larvae (ml) is estimated to be 17°C. This coincides with the average temperature of Szombathely in May, when the larvae starts to be detected in the field [63]. As ml also depends on the host abundance, this estimated is larger than the minimum temperature for the larval ticks to become active. We also note that larvae of Ixodes ricinus are observed to have the normal activity between 15°C-27°C [64].
Basic reproduction number
The basic reproduction number (R0) is a universally-recognized metric of the capacity of a parasite or a pathogen to reproduce given particular environmental conditions. We can use our parametrized model to evaluate the transmission potential of TBEV in the enzootic cycle by calculating basic reproduction number (R0) of TBE infection.
The calculation is detailed in S1 Appendix. Fig 3 shows the estimated R0 and the growing degree days in Szombathely. Between year of 1980 and 2015, R0 of TBE is estimated to be between 1.13-1.98. The estimated values are similar to the ones obtained from a modeling study in [23] and the study in a nearby region Borska nizina in Slovakia [7], where R0 ranged between 0.85 to 3.27. Sensitivity of R0 with respect to the estimated parameters are provided in S1 Appendix.
Solid blue line depicts 3-year average of GD4, sum of the daily mean temperatures exceeding 4°C. Dashed red line shows the basic reproduction numbers calculated at every 3 years from 1980 to 2016. Parameter values used for simulation are listed in Table 3, together with the estimated parameters at the best fit in Fig 2.
Trend
To see the trend of the transmission potential of TBEV in the enzootic cycle with respect to the change of temperatures, we computed R0 by assuming that (3) is a periodic system with a period of three years taking account of the average tick’s life cycle. We observe the increasing trend in both the growing degree days and the basic reproduction number. Our result clearly shows increase of the transmission potential of TBEV in the enzootic cycle along with the increased temperature between 1980 to 2013. This is in contrast with the observation that TBE human incidence rate has dropped in mid-1990s, see Fig 4. This discrepancy between transmission risk of TBEV in the ecological tick-host cycle and the TBE human case reporting indicates that other public health interventions have been effective in preventing human infection from a large pool of infected tick, unless the reporting rate of TBE has dropped since mid-1990s.
Blue bars represent TBE incidence (per 100000 people) in Hungary between 1980-2014 [38, 65, 66], with missing data in 2011. White bar shows the expected human case report based on the fitting result. Red line is the estimated R0. Parameter values used for computing R0 is in Table 3, together with the estimated parameters at the best fit in Fig 2.
Further quantification of co-feeding transmission
In order to study the significance of the non-systemic transmission route in the transmission of TBEV, we compared the estimated R0 of TBEV transmission with R0 of the system which excludes the non-systemic transmission, see Fig 5. As a result, values of R0 for the system excluding non-systemic transmission lied between 0.84-1.34 in year 1980 to 2015, which are 25%-33% less than the values of R0 for the system including both systemic and non-systemic transmission route of TBEV (1.13-1.98). The significance of co-feeding transmission estimated in our study is less than the estimates of [25], by which the non-systemic co-feeding pathway is estimated to be 60% greater degree of amplification of TBEV compared with the systemic pathway. Nonetheless, our analysis confirms that co-feeding transmission route is very significant, since the value of R0 for the system excluding the non-systemic transmission is near the threshold value 1 that determines the sustainability of the disease.
Blue and red lines show the basic reproduction numbers, including and excluding non-systemic transmission, respectively. Parameter values used for simulation are listed in Table 3, together with the estimated parameters at the best fit in Fig 2.
Discussion
TBE has been a major public health concern for much of the European countries. The risk of human infection depends on the prevalence of the virus spread in the tick-host ecological cycle, and on the tick-human contact that is effective for the disease transmission. Evaluating the human infection needs accurate estimation of the tick-host transmission, the tick population dynamics highly regulated by the environmental and climatic conditions, and the human case reporting given infection probability which is quite stochastic, all involving substantial uncertainties in ecological, epidemiological and environmental parameters. Here we developed a deterministic TBEV transmission model, coupled with a stochastic human case reporting given infection formulation, so that the relatively limited TBE incidence data and relatively complete climate data in the considered region (Hungary) can be used in conjunction with this cascade coupled model to infer key parameters aforementioned.
TBE incidence in Hungary has shown dramatic decrease in 1997. The reason may be a result of under-reporting followed by decreased serological examination or a result of public vaccination in early 1990s [10, 11]. Whether or not the decrease in TBE incidence reflects the actual decrease in human TBE infection level, it seems probable that TBEV transmission level has not been decreased in ecological cycle. The studies on the seropositivity of animals in 1960 to 1970 and 2005 support the argument [12, 13] that the prevalence of the TBEV has not decreased in the hosts. Our results is in line with this argument. The model simulations do not show significant decrease in the basic reproduction number of TBEV or the expected TBE cases between 1980 to 2015. One factor which has contributed to the decline of the incidence is the use of vaccines which protects human from infection upon tick-human contact. The annual TBE cases in a neighboring country Austria has also decreased, while the incidence in the unvaccinated population remains to be unchanged compared with the incidence in the pre-vaccination era [67]. In comparison, there is no vaccine available against Lyme disease transmitted by Ixodes ticks and there has been an increasing trend of Lyme disease case in Hungary since 1998 [68]. Our simulations suggest instead that the basic reproduction number in the tick-host transmission cycle has been increasing along with the increasing temperature. In addition, in 2007, the massive use of pesticides to control ticks has stopped due to regulations. This may have further increased the transmission potential of TBEV in the enzootic cycle by allowing better tick survival [69, 70].
We have also observed the significance of non-systemic transmission route on the TBEV transmission. This observation is consistent with the result from the early modeling study [23]. To our knowledge, this is the first modeling study which incorporates explicitly both seasonal dependence and the non-systemic transmission pathway in a single compartmental model setting. By considering seasonal dependence together with the non-systemic transmission pathway, we can examine the seasonal factors in the non-systemic transmission of the virus. For example, non-systemic transmission is only possible when both infective ticks and susceptible ticks are actively questing and the questing activity of ticks shows strong seasonal dependence. It should be mentioned that in our formulation of non-systemic transmission, we assumed that feeding ticks are equally distributed in all hosts. It has been observed that few hosts are attached with the most of the co-feeding ticks [25]. Also, the spatiotemporal distance between co-feeding ticks, which affects the transmissibility, is also neglected in our model formulation. Modification of the model is required to study the effect on disease dynamics of tick distribution over hosts, and the effect of spatiotemporal proximity of feeding ticks on hosts.
We quantified tick-to-host and tick-to-human infections in Hungary. The developed tick-borne disease transmission model and the methodologies we designed is expected to be adopted in the future studies assessing immunization programs for TBE in Hungary.
Supporting information
S1 Appendix. Interpretation of parameter p, normalized TBE transmission model, R0 computation, statistical inference.
https://doi.org/10.1371/journal.pone.0217206.s001
(PDF)
Acknowledgments
The author wishes to thank Viktor Zöldi and Katalin Krisztalovics (National Epidemiology Centre of Hungary, Department of Parasitology (Budapest)) for sharing epidemiological data.
References
- 1. Daniel M, Danielova V, Kříž B, Jirsa A, Nožička J. Shift of the tick Ixodes ricinus and tick-borne encephalitis to higher altitudes in central Europe. European Journal of Clinical Microbiology & Infectious Diseases 2003 May 25;22(5):327–8.
- 2. Zeman P, Benes C. A tick-borne encephalitis ceiling in Central Europe has moved upwards during the last 30 years: possible impact of global warming? International Journal of Medical Microbiology 2004 Apr 2;293:48. pmid:15146984
- 3. Campbell JA. Life history and development of the sheep tick ixodes ricinus linnaeus in Scotland, under natural and controlled conditions. PhD thesis, University of Edinburgh 1948.
- 4. Greenfield BP. Environmental parameters affecting tick (Ixodes ricinus) distribution during the summer season in Richmond Park, London. Bioscience Horizons 2011 May 4;4(2):140–8.
- 5. Mansfield KL, Johnson N, Phipps LP, Stephenson JR, Fooks AR, Solomon T. Tick-borne encephalitis virus–a review of an emerging zoonosis. Journal of General Virology 2009 Aug 1;90(8):1781–94. pmid:19420159
- 6. Perret JL, Guigoz E, Rais O, Gern L. Influence of saturation deficit and temperature on Ixodes ricinus tick questing activity in a Lyme borreliosis-endemic area (Switzerland). Parasitology Research 2000 Jun 1;86(7):554–7. pmid:10935905
- 7. Randolph SE, Miklisova D, Lysy J, Rogers DJ, Labuda M. Incidence from coincidence: patterns of tick infestations on rodents facilitate transmission of tick-borne encephalitis virus. Parasitology 1999 Feb;118(2):177–86. pmid:10028532
- 8. Randolph SE, Green RM, Hoodless AN, Peacey MF. An empirical quantitative framework for the seasonal population dynamics of the tick Ixodes Ricinus. International Journal for Parasitology 2002 Jul 1;32(8):979–89. pmid:12076627
- 9. Randolph SE. The shifting landscape of tick-borne zoonoses: tick-borne encephalitis and Lyme borreliosis in Europe. Philosophical Transactions of the Royal Society of London B: Biological Sciences 2001 Jul 29;356(1411):1045–56. pmid:11516382
- 10. Kollaritsch H, Chmelík V, Dontsenko I, Grzeszczuk A, Kondrusik M, Usonis V, et al. The current perspective on tick-borne encephalitis awareness and prevention in six Central and Eastern European countries: report from a meeting of experts convened to discuss TBE in their region. Vaccine 2011 Jun 20;29(28):4556–64. pmid:21549781
- 11. Ferenczi E, Rácz G, Faludi G, Czeglédi A, Mezey I, Berencsi G. Natural foci of classical and emerging viral zoonoses in Hungary. Nato Science Series I Life And Behavioral Sciences 2005 Sep 27;370:43.
- 12. Šikutová S, Hornok S, Hubálek Z, Doležálková I, Juřicová Z, Rudolf I. Serological survey of domestic animals for tick-borne encephalitis and Bhanja viruses in northeastern Hungary. Veterinary Microbiology 2009 Mar 30;135(3):267–71. pmid:19036537
- 13. Molnár E. Occurrence of tick-borne encephalitis and other arboviruses in Hungary. Geographia Medica 1982;12:78–120. pmid:6129178
- 14. Zöldi V, Juhász A, Nagy C, Papp Z, Egyed L. Tick-borne encephalitis and Lyme disease in Hungary: the epidemiological situation between 1998 and 2008. Vector-Borne and Zoonotic Diseases 2013 Apr 1;13(4):256–65. pmid:23428090
- 15. Pettersson JH, Golovljova I, Vene S, Jaenson TG. Prevalence of tick-borne encephalitis virus in Ixodes ricinus ticks in northern Europe with particular reference to Southern Sweden. Parasites & Vectors 2014 Dec;7(1):102.
- 16.
Donoso-Mantke O, Karan LS, Ruzek D. Tick-borne encephalitis virus: a general overview. Flavivirus Encephalitis. InTech 2011.
- 17.
Sonenshine DE, Mather TN. Ecological dynamics of tick-borne zoonoses. Oxford University Press 1994 Oct 20.
- 18. Bolzoni LU, Rosà RO, Cagnacci FR, Rizzoli AN. Effect of deer density on tick infestation of rodents and the hazard of tick-borne encephalitis. II: population and infection models. International Journal for Parasitology 2012 Apr 1;42(4):373–81. pmid:22429768
- 19. Ghosh M, Pugliese A. Seasonal population dynamics of ticks, and its influence on infection transmission: a semi-discrete approach. Bulletin of mathematical biology 2004 Nov 1;66(6):1659–84. pmid:15522350
- 20. Hancock PA, Brackley R, Palmer SC. Modelling the effect of temperature variation on the seasonal dynamics of Ixodes ricinus tick populations. International Journal for Parasitology 2011 Apr 1;41(5):513–22. pmid:21295037
- 21. Nonaka E, Ebel GD, Wearing HJ. Persistence of pathogens with short infectious periods in seasonal tick populations: the relative importance of three transmission routes. PLoS One 2010 Jul 23;5(7):e11745. pmid:20668521
- 22. Rosa R, Pugliese A. Effects of tick population dynamics and host densities on the persistence of tick-borne infections. Mathematical Biosciences 2007 Jul 1;208(1):216–40. pmid:17125804
- 23. Hartemink NA, Randolph SE, Davis SA, Heesterbeek JA. The basic reproduction number for complex disease systems: Defining R0 for tick-borne infections. The American Naturalist 2008 Apr 3;171(6):743–54. pmid:18462128
- 24. Zhang X, Wu X, Wu J. Critical contact rate for vector-host-pathogen oscillation involving co-feeding and diapause. Journal of Biological Systems 2017 Dec;25(4):657–75.
- 25. Randolph SE, Gern L, Nuttall PA. Co-feeding ticks: epidemiological significance for tick-borne pathogen transmission. Parasitology Today 1996 Dec 1;12(12):472–9. pmid:15275266
- 26. Randolph SE. Transmission of tick-borne pathogens between co-feeding ticks: Milan Labuda’s enduring paradigm. Ticks and Tick-borne Diseases 2011 Dec 1;2(4):179–82. pmid:22108009
- 27. Mogl C, de Mendonca PG, Harsch A, Heyl J. Potential for cofeeding transmission of tick-borne pathogens in common voles and roe deer-supportive molecular evidence from field samples. Acta Zoologica Bulgarica 2011 Dec 1;63(3):301–6.
- 28. Ogden NH, Nuttall PA, Randolph SE. Natural Lyme disease cycles maintained via sheep by co-feeding ticks. Parasitology 1997 Dec;115(6):591–9. pmid:9488870
- 29. Voordouw MJ. Co-feeding transmission in Lyme disease pathogens. Parasitology 2015 Feb;142(2):290–302. pmid:25295405
- 30. Jaenson TG, Hjertqvist M, Bergström T, Lundkvist Å. Why is tick-borne encephalitis increasing? A review of the key factors causing the increasing incidence of human TBE in Sweden. Parasites & Vectors 2012 Dec;5(1):184.
- 31. Pretzmann G, Radda A, Loew J. Studies of a natural focus of the spring-summer meningoencephalitis (tick-borne Encephalitis) in Lower Austria. 5. Additional investigations of the virus circulation in a natural focus. Zentralblatt fur Bakteriologie, Parasitenkunde, Infektionskrankheiten und Hygiene. 1. Abt. Medizinisch-hygienische Bakteriologie, Virusforschung und Parasitologie. Originale. 1964 Dec;194(4):431. pmid:5890645
- 32. Ogden NH, Bigras-Poulin M, O’callaghan CJ, Barker IK, Lindsay LR, Maarouf A, et al. A dynamic population model to investigate effects of climate on geographic range and seasonality of the tick Ixodes scapularis. International Journal for Parasitology 2005 Apr 1;35(4):375–89. pmid:15777914
- 33. Wu X, Duvvuri VR, Lou Y, Ogden NH, Pelcat Y, Wu J. Developing a temperature-driven map of the basic reproductive number of the emerging tick vector of Lyme disease Ixodes scapularis in Canada. Journal of Theoretical Biology 2013 Feb 21;319:50–61 pmid:23206385
- 34. King AA, de Cellès MD, Magpantay FM, Rohani P. Avoidable errors in the modelling of outbreaks of emerging pathogens, with special reference to Ebola. Proceedings of the Royal Society B: Biological Sciences 2015 May 7;282(1806):20150347. pmid:25833863
- 35. Magpantay FM, de Cellès MD, Rohani P, King AA. Pertussis immunity and epidemiology: mode and duration of vaccine-induced immunity. Parasitology 2016 Jun;143(7):835–49. pmid:26337864
- 36.
European Centre for Disease Prevention and Control. Factsheet about tick-borne encephalitis (TBE). https://ecdc.europa.eu/en/tick-borne-encephalitis/facts/factsheet (Accessed: 2018-10-09)
- 37. Balogh Z, Ferenczi E, Szeles K, Stefanoff P, Gut W, Szomor KN, et al. Tick-borne encephalitis outbreak in Hungary due to consumption of raw goat milk. Journal of Virological Methods 2010 Feb 1;163(2):481–5. pmid:19836419
- 38.
European Centre for Disease Prevention and Control. Epidemiological situation of tick-borne encephalitis in the European Union and European Free Trade Association countries. Stockholm: ECDC 2012.
- 39.
European Centre for Disease Prevention and Control. Ixodes ricinus—Factsheet for experts. https://ecdc.europa.eu/en/disease-vectors/facts/tick-factsheets/ixodes-ricinus (Accessed: 2018-06-11)
- 40. Matuschka FR, Fischer P, Musgrave K, Richter D, Spielman A. Hosts on which nymphal Ixodes Ricinus most abundantly feed. The American Journal of Tropical Medicine and Hygiene 1991 Jan 1;44(1):100–7. pmid:1996733
- 41. Ricker WE. Stock and recruitment. Journal of the Fisheries Board of Canada 1954 May 1;11(5):559–623.
- 42. Qviller L, Grøva L, Viljugrein H, Klingen I, Mysterud A. Temporal pattern of questing tick Ixodes ricinus density at differing elevations in the coastal region of western Norway. Parasites & Vectors 2014 Dec;7(1):179.
- 43. Rehacek J. Transovarial transmission of tick-borne encephalitis virus by ticks. Acta Virologica 1962;6(3):220–6. pmid:14038646
- 44. Brockmann SO, Oehme R, Buckenmaier T, Beer M, Jeffery-Smith A, Spannenkrebs M, et al. A cluster of two human cases of tick-borne encephalitis (TBE) transmitted by unpasteurised goat milk and cheese in Germany, May 2016. Eurosurveillance 2018 Apr 12;23(15):17–00336.
- 45. Labuda M, Danielova V, Jones LD, Nuttall PA. Amplification of tick-borne encephalitis virus infection during co-feeding of ticks. Medical and Veterinary Entomology 1993 Oct;7(4):339–42. pmid:8268488
- 46. Khasnatinov MA, Tuplin A, Gritsun DJ, Slovak M, Kazimirova M, Lickova M, et al. Tick-borne encephalitis virus structural proteins are the primary viral determinants of non-viraemic transmission between ticks whereas non-structural proteins affect cytotoxicity. PloS One 2016 Jun 24;11(6):e0158105. pmid:27341437
- 47. Wilhelmsson P, Lindblom P, Fryland L, Nyman D, Jaenson TG, Forsberg P, et al. Ixodes ricinus ticks removed from humans in Northern Europe: seasonal pattern of infestation, attachment sites and duration of feeding. Parasites & Vectors 2013 Dec;6(1):362.
- 48. Trájer A, Bede-Fazekas Á, Hufnagel L, Bobvos J, Páldy A. The paradox of the binomial Ixodes ricinus activity and the observed unimodal Lyme borreliosis season in Hungary. International Journal of Environmental Health Research 2014 May 4;24(3):226–45. pmid:23859420
- 49. Dobson AD, Finnie TJ, Randolph SE. A modified matrix model to describe the seasonal population ecology of the European tick Ixodes ricinus. Journal of Applied Ecology 2011 Aug 1;48(4):1017–28.
- 50. Gray JS The fecundity of Ixodes ricinus (L.)(Acarina: Ixodidae) and the mortality of its developmental stages under field conditions. Bulletin of Entomological Research 1981 Sep;71(3):533–42.
- 51. Zolotov PE. Sex ratio in ixodid tick collections from Leningrad Province. Parazitologiia 1981;15(5):469–71. pmid:6457273
- 52. Harrison A, Bennett NC. The importance of the aggregation of ticks on small mammal hosts for the establishment and persistence of tick-borne pathogens: an investigation using the R0 model. Parasitology 2012 Oct;139(12):1605–13. pmid:23036641
- 53. Dusbabek F, Borský I, Jelinek F, Uhlíř J. Immunosuppression and feeding success of Ixodes ricinus nymphs on BALB/c mice. Medical and Veterinary Entomology 1995 Apr;9(2):133–40.
- 54.
KNMI Climate Explorer Daily SZOMBATHELY mean temperature https://climexp.knmi.nl/ecatemp.cgi?id=someone@somewhere&WMO=2042&STATION=SZOMBATHELY&extraargs (Accessed: 2017-10-17)
- 55. Zavadska D, Anca I, André F, Bakir M, Chlibek R, Čižman M, et. al. Recommendations for tick-borne encephalitis vaccination from the Central European Vaccination Awareness Group (CEVAG). Human Vaccines & Immunotherapeutics 2013 Feb 1;9(2):362–74.
- 56. King AA, Nguyen D, Ionides EL. Statistical inference for partially observed Markov processes via the R package pomp. Journal of Statistical Software 2016;69(12).
- 57.
R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing 2018.
- 58. Daniel M, Danielová V, Fialová A, Malý M, Kříž B, Nuttall PA. Increased relative risk of tick-borne encephalitis in warmer weather. Frontiers in Cellular and Infection Microbiology 2018 Mar 22;8:90. pmid:29623261
- 59. Zöldi V, Papp T, Reiczigel J, Egyed L. Bank voles show high seropositivity rates in a natural TBEV focus in Hungary. Infectious Diseases 2015 Mar 4;47(3):178–81. pmid:25549698
- 60. Zöldi V, Papp T, Rigó K, Farkas J, Egyed L. A 4-Year Study of a Natural Tick-Borne Encephalitis Virus Focus in Hungary, 2010–2013. EcoHealth 2015 Mar 1;12(1):174–82. pmid:25227182
- 61. Bogovic P, Strle F. Tick-borne encephalitis: a review of epidemiology, clinical characteristics, and management. World Journal of Clinical Cases 2015 May 16;3(5):430. pmid:25984517
- 62. Eisen L, Eisen RJ. Critical evaluation of the linkage between tick-based risk measures and the occurrence of Lyme disease cases. Journal of Medical Entomology 2016 Jun 21;53(5):1050–62.
- 63. Egyed L, Élő P, Sréter-Lancz Z, Széll Z, Balogh Z, Sréter T. Seasonal activity and tick-borne pathogen infection rates of Ixodes ricinus ticks in Hungary. Ticks and Tick-borne Diseases 2012 Apr 1;3(2):90–4. pmid:22445929
- 64.
Sixl W, Nosek J. Einfluss von. Temperatur und Feuchtigkeit auf das Verhalten von Ixodes ricinus. Dermacentor marginatus und Haemaphysalis inermis 1971.
- 65. Randolph SE. Tick-borne encephalitis incidence in Central and Eastern Europe: consequences of political transition. Microbes and Infection 2008 Mar 1;10(3):209–16. pmid:18316221
- 66.
European Centre for Disease Prevention and Control. Annual epidemiological report for 2015. Tick-borne encephalitis. European Centre for Disease Prevention and Control 2018.
- 67. Heinz FX, Stiasny K, Holzmann H, Kundi M, Sixl W, Wenk M, Kainz W, Essl A, Kunz C. Emergence of tick-borne encephalitis in new endemic areas in Austria: 42 years of surveillance. Eurosurveillance 2015 Apr 2;20(13):21077. pmid:25860391
- 68. Trájer A, Bobvos J, Páldy A, Krisztalovics K. Association between incidence of Lyme disease and spring-early summer season temperature changes in Hungary-1998-2010. Annals of Agricultural and Environmental Medicine 2013;20(2). pmid:23772569
- 69. Medlock JM, Hansford KM, Bormane A, Derdakova M, Estrada-Peña A, George JC, et al. Driving forces for changes in geographical distribution of Ixodes ricinus ticks in Europe. Parasites & Vectors 2013 Dec;6(1):1.
- 70. National Center of Epidemiology, Hungary. A kullancsok elleni védekezésről. Epidemiológiai Információs Hetilap 2009 May 4;16.