Modelling the transmission dynamics of severe fever with thrombocytopenia syndrome in Jiangsu Province, China.

BACKGROUND: Severe fever with thrombocytopenia syndrome (SFTS) is an emerging infectious disease that is regionally distributed in Asia, with high fatality. Constructing the transmission model of SFTS could help provide clues for disease control and fill the gap in research on SFTS models.
METHODS: We built an SFTS transmission dynamics model based on the susceptible-exposed-infectious-asymptomatic-recovered (SEIAR) model and the epidemiological characteristics of SFTS in Jiangsu Province. This model was used to evaluate the effect by cutting off different transmission routes and taking different interventions into account, to offer clues for disease prevention and control.
RESULTS: The transmission model fits the reported data well with a minimum R2 value of 0.29 and a maximum value of 0.80, P < 0.05. Meanwhile, cutting off the environmental transmission route had the greatest effect on the prevention and control of SFTS, while isolation and shortening the course of the disease did not have much effect.
CONCLUSIONS: The model we have built can be used to simulate the transmission of SFTS to help inform disease control. It is noteworthy that cutting off the environment-to-humans transmission route in the model had the greatest effect on SFTS prevention and control.

Background

Severe fever with thrombocytopenia syndrome (SFTS) is an emerging infectious disease caused by a novel bunya virus called severe fever with thrombocytopenia syndrome virus (SFTSV) [1]. The clinical manifestations of infection in human are characterized by high fever, a drastic reduction of platelets and leukocytes and hemorrhage [2-5]. The clinical course is mainly divided into three stages: fever, multiple nerve dysfunction and recovery [1, 3, 4, 3, 7].

In previous years, SFTS occurred mainly in the rural areas of central and northeast China, with an estimated 500–1500 cases recorded in 2011–2012. Over time, SFTS has spread to South Korea and Japan, although the related fatality rate of nearly 30% has dropped to 10–12% in recent years [8]. In South Korea, the fatality rate was above 35% among infected patients in 2012–2013, but the incidence rate dropped to 11.52% in 2016 [9]. In Japan, however, the fatality rate is slightly higher than that observed in China (23%), and cases have mainly occurred in the western regions [10]. Currently, the SFTS fatality rate in Asia seems to be decreasing, while the incidence rate is increasing regionally, especially in China. According to the monitoring data from the China Information Network System of Disease Prevention and Control, SFTS incidence and prevalence rates have increased each year from 2010 to 2016, with cases found in 23 provinces [11]. Analyzing the epidemiological characteristics of SFTS, one published study found Henan Province, Hubei Province and Anhui Province to be the main provinces with the largest numbers of cases before 2011 [12]. Including five other provinces, SFTS cases were detected in Jiangsu Province as far back as 2011 [13]. Since then, the number of SFTS cases reported in Jiangsu Province has increased annually, with more than 140 confirmed cases in over 60 locations within the province by 2016 [14].

To fight against this emerging infectious disease, findings from both domestic and foreign research on the etiology and risk factors of the disease are of great importance. As shown by these study findings, SFTS mainly occurred in tick-active months from May to September, and reports of new infections continued until November. The main population affected by SFTS is farmers aged 50 years old and above who live in forested or hilly areas [15-18]. Recent researches suggested that ticks are the main host of SFTS and play a key role in its transmission. Haemaphysalis longicornis, Ixodes nipponensis, Rhipicephalus microplus and Dermacentor sinicus, amongst others, are the susceptible ticks identified, but Haemaphysalis longicornis is the dominant species [1, 19–22]. While SFTS cases occurred from March to November, the disease period peaked in August and September. The seasonal distribution of SFTS cases was also found to be closely related to the seasonal trend of tick density [23, 24]. Some preliminary researches exist to show the process of host animal–tick–human natural transmission routes of SFTS [25]. However, there are very few systematic and definite researches on transmission routes among vector, host and SFTSV-affected populations. Meanwhile, various problems exist in the treatment of SFTS such as unclear pathogenesis, absence of vaccine for prevention and lack of a specific SFTS antiviral therapy. As a result, a low fatality rate cannot be maintained.

SFTS infection also induces a certain amount of economic pressure on affected populations, as a study has shown that patients and their families pay nearly $3158.4 for treatment [2]. This makes SFTS a significant threat to public health, as the treatment amount is almost equal to the disposable income of rural residents in Anhui Province as of 2017 [2]. To date, no research exists on a transmission dynamics model of SFTS, and there is no clear evidence on preventive methods that could help achieve optimal intervention effects.

This study aimed to construct a multi-population and multi-route dynamic model (MMDM) based on the susceptible–exposed–infectious–asymptomatic–recovered (SEIAR) model, the epidemiological characteristics and the transmission features of SFTS in Jiangsu Province.

Methods

Study area

Jiangsu Province is located between 30°45′ and 35°20′ north latitude and 116°18′ to 121°57′ east longitude (Additional file 1: Figure S1). Jiangsu Province is located in the East Asian monsoon climate zone, in the transitional zone between subtropical and warm temperate climate. The total area of Jiangsu Province is 107,200 square kilometers, with an estimated permanent population of 80.7 million as of the end of 2019.

Study design

This study was designed by three steps. First step, the MMDM model of SFTS was established by searching scientific literature and analyzing reported data. Second step, the reported data of SFTS in Jiangsu Province from January 2001 to December 2019 was used to fit the developed model. Third step, the MMDM model was used to simulate the intervention effect of eliminating different routes of transmission. This was achieved by setting the parameters of the model to obtain the best intervention plan, thus providing clues to effective disease control interventions and filling the gap in the research on SFTS models (Fig. 1).

Fig. 1

Flowchart of study design

Establishment of transmission dynamics model

Based on the characteristics of epidemic and the routes of transmission of SFTS in Jiangsu Province and the SEIAR model [14, 26–28], this study established the MMDM to simulate SFTS transmission. This model included four parts, namely human, ticks, host animals and environment, which are represented by subscript 1, 2, 3 and w, respectively (Fig. 2).

Fig. 2

Flowchart of the MMDM without intervention

Human part of the MMDM

In this part, humans were categorized as: susceptible (S1), exposed (E1), asymptomatic (A1), infectious (I1) or recovered/remover (R1). The structure of the MMDM is shown in Fig. 2, and this part of the MMDM was based on the following assumptions:

This model supposed that all human newborns are susceptible to the disease. Under this assumption, the total population is denoted by N, the birth rate is br and the natural mortality rate is dr. According to the literature, the fatality rate of overt infections, f, is 16.1% [5]. Meanwhile, this model took the population-wide mortality rate to approximately replace the mortality rate of the excluded overt infections.

According to the literature, people can be infected with SFTS by coming in contact with the blood of an infected patient [29]. At the same time, the transmissibility of I1 and A1 are different. Therefore, we assumed that the transmissibility of A1 is κ times that of I1 κ. The infection rate coefficients of I1 and A1 are β1 and κβ1, respectively.

We assumed that certain percentage (1 − p) (0 ≤ p ≤ 1) of E1 will transform into I1 after an incubation period, while the other part (p) of E1 will become A1 after a latent period. Hence, at time t, the development speed of E1 to the I1 pathway is proportional to the number of people in the incubation period, and the proportional coefficient is (1 − p) ω1. The development speed of E1 to the A1 pathway is proportional to the number of people in the latent period, and the proportional coefficient is pω1.

The out speed of the population in the I1 and A1 groups is respectively proportional to their population number, and the ratio coefficients are respectively γ and γ'.

The infectors of SFTS will produce antibodies after they are removed, but this study supposed that the R1 will never be infected and that the R1 is the destination of the transmission route between people, because literature on a second infection in the population was not found.

Environment part of the MMDM

In this part, W represent the environmental transmission route in the MMDM, and this part was based on the following assumptions:

1. Although no literature strongly proved that the environment is one of the transmission route, some evidence suggested that the possibility of getting SFTS infection SFTS in the environment was higher than by ticks [30-33]. Hence, we assume that W has the ability to facilitate the spread of SFTS to humans and animals. We set the transmissibility of W to animals and human as β [34].

2. As reported, the blood and excreta of infected human or host animals may diffuse the SFTSV. And as SFTS has commonly occurred extensively in the wild, we assumed that animals have the possibility to infect the environment by their blood or excreta. We set the propagation coefficient j as 10 [35].

Tick part of the MMDM

In this part, the ticks were divided into susceptible ticks (S2), exposed ticks (E2) and asymptomatic ticks (A2). And this part of the MMDM was based on the following assumptions:

According to literature, SFTS mainly spreads by ticks [1, 19–21], and infected ticks can spread the virus by biting host animals or humans [36].

As the long-term storage host, ticks infected with SFTSV can preserve the pathogen and retain the infectious ability of the pathogen without symptoms [36]. So this study set the transmission route among ticks as the susceptible–exposed–asymptomatic (SEA) model.

In transmission among ticks, A2 has transmissibility to infect ticks, and the infection rate coefficient is β2, but our study considered that the SFTS was not diffused among ticks, so we set β2 as 0. And the S2 will become E2 by biting the infected animals and humans.

This study supposed that after one latent period, all the E2 will transform to A2. So at time t, the number of ticks in the E2 to A2 route is ω2E2.

Ticks in the state of latent infection can infect host animals through biting. As the coefficient of infection rate of the virus spreading from ticks to a host animal is β2-3, we denoted the part of host animal infection by a tick’s bite as β2-3S3A2. In addition, ticks can also spread SFTS among humans through direct bites. Since the coefficient of infection rate of the virus spreading from ticks to humans is β2-1, the part of humans getting the virus by tick bites is denoted as β2-1S1A2.

Host animal part of the MMDM

In this part, the animals were divided into susceptible animals (S3), exposed animals (E3) and asymptomatic animals (A3). The assumptions of this part in the MMDM were as follows:

According to reports by one study, the main amplification hosts of arboviruses are small rodents. This meant that arthropods can cause relative viremia by biting small rodents. However, other susceptible arthropods can also be infected and make the virus spread further by biting infected small rodents [37].

For SFTS, ticks infected with SFTSV can cause viremia by biting host animals. At the same time, susceptible ticks can also be infected by biting infected host animals [37-40]. Therefore, we still set the transmission progress of SFTSV as the SEA model in the host animal part of the transmission dynamics model.

Among the transmission of host animals, A3 has transmissibility that can make S3 transform to E3 with the infection rate coefficient β3. At the time t, the number of host animals in the S3 to E3 pathway was β3S3A3.

This study supposed that all the E3 will transform to A3 after one latent period. At the time t, the number of host animals in the E3 to A3 pathway was ω3E3.

The uninfected ticks can be infected after biting the host animal that is in a state of latent infection. The coefficient of host-to-tick transmission of the virus is β3-2, so ticks with a β3-2S2A3 component can acquire the virus by biting an infected host animal.

The differential equation

The above MMDM is expressed by the following differential equations:

The left-hand side of the equation was respectively expressed as the instantaneous change velocity of S, E, I, A and R at time t.

Estimation of parameters

In our MMDM, there were a total of 20 parameters. The values of parameters p,1/ω1,1/γ,1/γ', f,1/ω2, 1/ω3 and j were all according to relative literatures. Among them, p of recessive infection ranged from 3.2% to 5.5% [41], and the mean was 4.3% in this study, i.e., p = 0.043. The 1/ω1 of the incubation period of SFTS in the population ranged from 9 to 14 days [42, 43], and the median was 11 days in this study. The course of the disease 1/γ of overt infection was 2 weeks [44], i.e. 14 days. The research about the course of the disease 1/γ' of latent infection was not found, so we used the same days with the course of the disease of overt infection, i.e., 14 days. The fatality rate f ranged from 11.2% to 30%, and f = 0.16 [5]. The incubation period after tick infection with SFTSV was 1/ ω2 weeks, i.e., 7 days. The incubation period of SFTSV infection in host animals 1/ω3 was 12 days [39]. No valid data or literature support was found for the relative transmissibility coefficient κ of recessive infected persons. In this study, κ = 1 was selected for calculation. The br and dr were obtained from the statistical yearbook of Jiangsu Province, see Table 1.

Table 1

Definition and value of parameters

Parameter	Definition	Value	Unit	Method
Person
β1	Person-to-person transmissibility coefficient	–	(Person*day)−1	Model fitting
κ	Relative transmissibility coefficient of unapparent infection	1	1	–
p	Proportion of latent infection	0.043	1	Ref. [41]
1/ω1	Incubation period	11	Day	Ref. [42]
1/γ	Infectious period of dominant infection	14	Day	Ref. [44]
1/γ'	Latent infection period	14	Day	–
f	Fatality rate	0.16	1	Ref. [5]
br	Birth rate	0.009481	1	Statistical Yearbook of Jiangsu Province
dr	Mortality rate	0.007003	1	Statistical Yearbook of Jiangsu Province
Ticks
β2	Coefficient of transmissibility between ticks	0	(Pieces*day)−1	–
β2-1	Coefficient of tick-to-human transmissibility	8β1	(Person*day)−1	Ref. [34]
β2-3	Coefficient of tick-to-host infection	8β1	(Pieces*day)−1	Ref. [34]
1/ω2	Incubation period of ticks	7	Day	Ref. [26]
Host animal
β3	Coefficient of transmissibility between host animal	2β1	(Pieces*day)−1	Ref. [34]
β3-1	Coefficient of transmissibility rate of host animal to human	2β1	(Pieces*day)−1	Ref. [34]
β3-2	Coefficient of host animal infection to ticks	8β1	(Person*day)−1	Ref. [34]
1/ω3	Host animal incubation period	12	Day	Ref. [39]
j	Rate of host animal discharge to the environment	10	–	Ref. [35]
Environment
βw	Coefficient of transmissibility rate of the environment to humans and host animals	9β1	(Human/animal*day)−1	Ref. [34]
W	Environmental tick density	0.047494291	One km2/person (flag hour)	Jiangsu Province surveillance data

The coefficient of infection β1was generated by model fitting. The Runge–Kutta method of order 4 with tolerance set at 0.001 was used to perform curve fitting of the root mean square deviation between the data and best run so far. In this study, the seasonality of the transmission was considered. According to the MMDM, the seasonality should be dynamic, focusing on β. Therefore, a trigonometric function was adopted and shown as follows:

In the equation, β0, t, α and T refer to the baseline of the transmission relative, time, a constant which adjusts the position of time, and the time span of the season cycle, respectively.

“Knock-out” simulation

“Knock-out” simulation was performed as the reference [45] to quantify cutting different transmission routes. To “knock out” means to cut off the transmission routes between different groups. In the model study, the simulation sets a parameter to 0 and estimates the contribution of the parameter by calculating the number of reduced cases or the total attack rate. For example, in the MMDM, the contribution of environmental transmission β was set to 0, and its effect was reflected by calculating the number of reduced cases.

In this study, the transmissibility coefficient in the MMDM was used to reflect different transmission routes, including transmissibility of the environment to humans or host animals (β), host animals to humans (β3-1), ticks to humans (β2-1), humans to humans (β1), ticks to host animals (β2-3), host animals to ticks (β3-2) and host animals to host animals (β3).

Based on the established SFTS transmissibility model, using 2019 data, a scenario of cutting off various transmission routes was established by setting one of the above infection rates up to seven values or to 0 at the time of simulation. We arranged the matrix between different transmission routes, as shown in Table 2, to perform a more comprehensive re-sectioning of transmission routes. There were seven scenarios at the same time:

Table 2

Combined transmissibility matrix

	β1	β3	βw	β21	β31	β32	β23
β1	1	2β1	9β1	8β1	2β1	8β1	8β1
β3	1/2β3	1	9/2β3	4β3	1β3	4β3	4β3
βw	1/9βw	2/9βw	1	8/9βw	2/9βw	8/9βw	8/9βw
β21	1/8β21	1/4β21	9/8β21	1	1/4β21	1β21	1β21
β31	1/2β31	1β31	9/2β31	4β31	1	4β31	4β31
β32	1/8β32	1/4β32	9/8β32	1β32	1/4β32	1	1β32
β23	1/8β23	1/4β23	9/8β23	1β23	1/4β23	1β23	1

Scenario 1: Only one coefficient β was controlled to be zero.

we made β or β3-1 = 0 or β2-1 = 0 or β1 = 0 or β2-3 = 0 or β3-2 = 0 or β3 = 0. At the same time, we kept the other six β coefficients unchanged.

Scenario 2: Two coefficients β were controlled to be zero.

Then, we combined these seven coefficients β in pairs so that two of them are equal to zero, while keeping the other five coefficients β unchanged.

In the following scenarios, the method adopted was the same as the above two scenarios, that is, the β value of the corresponding number was controlled to be zero, and the remaining β values remain unchanged. Scenario 3: Three coefficients β were controlled to be zero. Scenario 4: Four coefficients β were controlled to be zero. Scenario 5: Five coefficients β were controlled to be zero. Scenario 6: Six coefficients β were controlled to be zero. Scenario 7: Seven coefficients β were controlled to be zero.

The above methods were used to obtain the intervention effect of cutting off different transmission routes of SFTS. The indicators of intervention effect included: cumulative number of cases, total attack rate, peak incidence site and peak incidence. This was to provide guidance for the prevention and control of the occurrence of an SFTS epidemic.

Effectiveness of interventions

In this model, we only considered two interventions, and the model was established as shown in Fig. 3. We were trying to isolate the infections according to the transmission route of SFTS from person to person. In addition, although there are no targeted antiviral drugs to treat SFTS, drug treatment can still reduce the mortality rate and increase the clinical recovery rate. Therefore, considering interventions that can shorten the course of the disease, such as hospitalization or drug treatment, was important [46, 47].

Fig. 3

Flowchart of the MMDM with interventions

For the isolation, Q was the number of people isolated, R was the recovery of isolation, δ was the isolation coefficient and c was the coefficient of shortening the disease course in the model. We set the ϕ as the proportion of isolation which can be calculated by ; we then set ϕ as 98%, 90%, 80%, 70%, 60% and 50% to simulate the solo effect of different isolation ratios.

For the intervention to shorten the disease courses, we used γ to represent this intervention. We set γ as 1/14, 1/12, 1/10, 1/8, 1/6 and 1/4 to simulate the solo effect of shortening different days.

We did not know how many isolated cases and how many days of medication would shorten the course of the disease. Therefore, we set ϕ as 0.5%, 5% and 50%, and set γ as 1/14, 1/12, 1/10, 1/8, 1/6 and 1/4 to simulate the combined effect of SFTS under different intervention intensity.

The indicators of intervention effect included: cumulative number of cases, incidence rate, peak position of incidence and peak size of incidence. Finally, the model provided the optimal way to intervene in the spread of SFTS and provided guidance for the prevention and control of the outbreak of SFTS.

Model simulation and statistical analysis

In this study, Excel was used for data summary and chart drawing, IBM SPSS Statistics (version 21.0.0.0) was used for data statistical analysis, Berkeley Madonna (Version 8.3.18) was used for fitting existing data and models, and epidemic prediction and evaluation were conducted for intervention effects under different transmission routes.

Results

Model fitting

The result of model fitting using 2011–2019 SFTS incidence rates in Jiangsu Province was shown in Fig. 4. By analyzing every month's and year's SFTS incidence rates, we discovered that the incidence rate of SFTS strictly changed seasonally. This meant that in a year, the incidence rate was high in summer and autumn with two peaks, and lowered with no cases in winter. Therefore, the seasonal correction function was added to carry on fitting in the model. Just like the results of our model fitting, the SFTS incidence rate of 2011–2019 in Jiangsu Province was periodically changed. It is often at rest in winter, begins to increase in spring and peaks by July–August. Notably, the peak incidence rate of SFTS in Jiangsu Province has increased to some extent in recent years compared with 2011–2014. The R2 and P values of the model fitting results for 2011–2019 are listed in Table 3. The model fitting results for all years were good, and P was significant (P < 0.05).

Fig. 4

Model fitting results of the model to the reported SFTS incidence data from 2011 to 2019 in Jiangsu Province, China

Table 3

Model fitting result of 2011–2019

Year	R2	P
2011	0.573878	< 0.05
2012	0.299520	< 0.05
2013	0.482155	< 0.05
2014	0.764983	< 0.05
2015	0.408163	< 0.05
2016	0.767604	< 0.05
2017	0.776788	< 0.05
2018	0.801787	< 0.05
2019	0.758567	< 0.05

Model simulation of cutting transmission routes

Table 4 showed the simulation results of truncating different transmission routes based on the SFTS transmission dynamics model in Jiangsu Province in 2019. In general, when the cut transmission route has no effect, the cumulative number of cases was 671, and the peak time was 205 days.

Table 4

Interventional effect of cutting off different transmission routes on SFTS

“Knock-out” scenario*	Cumulative number of cases	Total attack rate	Peak time (days)	Peak attack rate
β1	0	0	0	0
β3	671	0	205	0
βw	0	0	198	0
β21	671	0	205	0
β31	671		205	0
β23	671	0	205	0
β32	671	0	205	0
β3βw	0	0	198	0
β3β21	671	0	205	0
β3β31	671	0	205	0
β3β23	671	0	205	0
β3β32	671	0	205	0
βwβ21	0	0	198	0
βwβ31	0	0	10	0
βwβ23	0	0	198	0
βwβ32	0	0	198	0
β21β23	671	0	205	0
β21β31	671	0	205	0
β21β32	671	0	205	0
β31β23	671	0	205	0
β31β32	671	0	205	0
β23β32	671	0	205	0
β3βwβ21	0	0	198	0
β3βwβ31	0	0	10	0
β3βwβ23	0	0	198	0
β3βwβ32	0	0	198	0
β3β21β31	671	0	205	0
β3β21β32	671	0	205	0
β3β21β23	671	0	205	0
β3β31β23	671	0	205	0
β3β31β32	671	0	205	0
β3β23β32	671	0	205	0
βwβ21β31	0	0	10	0
βwβ21β23	0	0	198	0
βwβ31β23	0	0	10	0
βwβ21β32	0	0	198	0
βwβ31β32	0	0	10	0
βwβ23β32	0	0	198	0
β21β31β23	671	0	205	0
β21β31β32	671	0	205	0
β21β23β32	671	0	205	0
β31β23β32	671	0	205	0
β3βwβ21β31	0	0	7	0
β3βwβ21β23	0	0	198	0
β3βwβ21β32	0	0	198	0
β3βwβ23β32	0	0	198	0
β3βwβ31β23	0	0	10	0
β3βwβ31β32	0	0	10	0
β3β21β31β23	671	0	205	0
β3β21β23β32	671	0	205	0
β3β21β31β32	671	0	205	0
β3β31β23β32	671	0	205	0
βwβ21β31β23	0	0	10	0
βwβ21β31β32	0	0	10	0
βwβ21β23β32	0	0	198	0
βwβ31β23β32	0	0	10	0
β21β31β23β32	671	0	205	0
β3βwβ21β31β23	0	0	10	0
β3βwβ21β31β32	0	0	10	0
β3βwβ21β23β32	0	0	198	0
β3βwβ31β23β32	0	0	10	0
β3β21β31β23β32	671	0	205	0
βwβ21β31β23β32	0	0	10	0
β3βwβ21β31β23β32	0	0	10	0

*“Knock-out” means to simulate cutting off the transmission routes by setting a parameter to 0, and estimates the contribution of the parameter by calculating the number of reduced cases or the total attack rate. The different scenario means which parameter was set as 0

We mainly studied the spread of SFTS in the population. When the transmission route between people β1 was cut off, SFTS cannot be transmitted in the population. With zero cases, the model cannot continue to simulate, and hence, other combinations containing β1 were not shown. When cutting off a single transmission route, it is worth noting that when the environmental transmission β was cut off alone, the number of cases also became zero. None of the other transmission routes cut off alone reduced the number of SFTS infections. When two transmission routes, such as a combination with environmental transmission were cut off, the number of patients was also reduced. When the environmental transmission β and animals-to-humans transmission β31 were cut off together, the peak time of disease was advanced to 10 days. When other combinations that can cause a decrease in SFTS cases were cut off, the peak time was advanced to 198 days. When three transmission routes (such as the combination of environmental transmission β, animals-to-humans transmission β31 and ticks-to-humans transmission β21) were cut off at the same time, the number of cases decreased significantly. And similarly, as above, when β31 and βw were cut off, the peak time was advanced to 10 days. With a cut-off of environmental transmission β and animals-to-ticks transmission β32, ticks-to-animals transmission β23 and ticks-to-humans transmission β21 at the same time, the peak time was also increased to 198 days. When four transmission routes were cut off at the same time (such as β3, β, β21 and β31), the peak time of SFTS was advanced to 7 days. The remaining results were the same as above. When five and six transmission routes were cut off at the same time, SFTS only occurred when environmental transmission was not cut off, and the peak onset time was not advanced. For more details, see Table 4.

Taking different interventions

In the study of separate interventions, the shortening of the disease course was assumed to be no shortening and shortening of 2 days, 4 days, 6 days, 8 days and 10 days; the isolation ratio was assumed to be: 98%, 90%, 80%, 70%, 60% and 50%. When studying the combined effects of the two interventions, the isolation ratios were assumed to be 0.5%, 5% and 50%. It turned out that when a separate intervention was implemented alone, shortening the course of the disease will not limit the spread of fever with thrombocytopenia syndrome. But when the quarantine was implemented, only one case was reduced. When the joint intervention was implemented (i.e., the course of the disease was not shortened and the isolation ratio was 5% and 50%), only one case was reduced. Therefore, these two interventions had little effect on the number of cases of fever with thrombocytopenia syndrome, as shown in Tables 5 and 6.

Table 5

Interventional effect of taking joint interventions

Intervention	Cumulative number of cases	Total attack rate	Peak time (days)	Peak attack rate
γ = 1/14 + ∅ = 0	530	0.00	197	0
γ = 1/14 + ∅ = 0.5%	531	0.00	197	0
γ = 1/14 + ∅ = 5%	530	0.00	197	0
γ = 1/14 + ∅ = 50%	530	0.00	197	0
γ = 1/12 + ∅ = 0	531	0.00	197	0
γ = 1/12 + ∅ = 0.5%	531	0.00	197	0
γ = 1/12 + ∅ = 5%	531	0.00	197	0
γ = 1/12 + ∅ = 50%	531	0.00	197	0
γ = 1/10 + ∅ = 0	531	0.00	197	0
γ = 1/10 + ∅ = 0.5%	531	0.00	197	0
γ = 1/10 + ∅ = 5%	531	0.00	197	0
γ = 1/10 + ∅ = 50%	531	0.00	197	0
γ = 1/8 + ∅ = 0	531	0.00	197	0
γ = 1/8 + ∅ = 0.5%	531	0.00	197	0
γ = 1/8 + ∅ = 5%	531	0.00	197	0
γ = 1/8 + ∅ = 50%	531	0.00	197	0
γ = 1/6 + ∅ = 0	531	0.00	197	0
γ = 1/6 + ∅ = 0.5%	531	0.00	197	0
γ = 1/6 + ∅ = 5%	531	0.00	197	0
γ = 1/6 + ∅ = 50%	531	0.00	197	0
γ = 1/4 + ∅ = 0	531	0.00	197	0
γ = 1/4 + ∅ = 0.5%	531	0.00	197	0
γ = 1/4 + ∅ = 5%	531	0.00	197	0
γ = 1/4 + ∅ = 50%	531	0.00	197	0

aγ represent the intervention of shorten disease course and ϕ represent the proportion of isolation

Table 6

Interventional effect of taking isolation

Isolation ratio	Cumulative number of cases	Cumulative attack rate	Peak time (days)	Peak attack rate
50%	530	0.00	201	0.00
60%	530	0.00	201	0.00
70%	530	0.00	201	0.00
80%	530	0.00	201	0.00
90%	530	0.00	201	0.00
98%	530	0.00	201	0.00

Discussion

SFTS is an emerging infectious disease caused by SFTSV. Recently, the incidence rate has been increasing year by year, and the coverage of the epidemic area has expanded. Jiangsu Province is one of the provinces with the highest incidence of SFTS, and they have conducted stringent SFTS monitoring. This study was established to use a model based on the SFTS monitoring data of Jiangsu Province from 2011 to 2019. The SFTS transmission dynamics model was used to simulate the change in incidence rate by cutting different transmission routes, and taking interventions into consideration and assessing the effect of those interventions.

The effect of knock-out simulation

Previous studies have shown that the number of SFTS patients infected by the environment were more than those who were infected by ticks [19, 30–33]. But most studies only use the density of ticks in the environment and the antibody of SFTSV in ticks to evaluate the hazardous situation of the environment [22, 48–50], so they did not consider the environment as one of the transmission route. As we all know, several factors such as the number of ticks, the dominant phase of ticks, outdoor activities and the change of weather can affect the occurrence of SFTS [51]. There was one study that collected SFTSV from inanimate objects in the room of SFTS patients. Thus, the environment could be considered as one of the transmission routes of SFTS, although a test result of the air sample was negative [52]. Therefore, based on the literature and the 2019 data of SFTS incidents in Jiangsu Province, we established the SFTS transmission dynamics model and used the coefficient of transmissibility to reflect the different routes of transmission. We simulated cutting different transmission routes to show the different effects of interventions. The result of this study showed that in all the combinations of cutting different transmission routes, the cumulative number of SFTS cases and other indicators could be significantly reduced when the transmission route containing the environment was cut off. Moreover, the intervention effect was more obvious when the environmental and host animal transmission routes were cut simultaneously. Those results indicated that the environment is an important route for the transmission of SFTS, and this was consistent with findings of previous research. Besides, this also provided a direction that can be explained as follows. It is possible that people were not exposed to a wild environment, and the chance of contacting ticks or animal excreta, blood, etc. decreased and accounted for the low SFTS incidence rate. In addition, some literatures have shown that SFTS was an environment-tick-host animal-tick cycle and that humans were infected by accidental participation in one of the steps. Therefore, the more transmission ways cut off, the smaller the impact of SFTS was. Interestingly, the effect of cutting off the environment and the transmission route of ticks on SFTS was less than that of cutting off the environment and host animals simultaneously. This might be the main way for people to be infected with SFTS except for tick bites. In previous studies, there were also literatures showing that those infected with SFTS with a clear history of tick bites only accounted for a small proportion (about 22–29% [6, 33, 53]) of all infected individuals. The study also suggested that humans might be infected with SFTS after exposure to animal blood [6]. Although the transmission of SFTS was initiated by ticks, it may be that the two transmission routes of animal and environment are more important for human beings.

The effect of taking interventions

For infectious diseases, there are three links and two factors. The three links are: the source of infection, transmission routes and the population. The two factors include natural factors and social factors. Common measures to control infectious diseases include controlling the source of infection, cutting off the route of transmission and protecting the population. Therefore, after cutting off the route of transmission and protecting the population, there is still a need to control the source of infection. In that regard, we chose to take interventions including taking isolation and shortening the course of the disease. We also assessed the effects of taking isolation and shortening the course of the disease according to the established SFTS transmissibility model. Based on the data of SFTS, the results of the simulation were used to evaluate the effect of the two interventions on SFTS. The results showed that whether implementing isolation alone or shortening the course of the disease, or combining the two interventions, the effect of interventions mentioned above on the incidence of SFTS was limited. And the different isolation ratio had little effect on the cases. This may be because that the incidence of SFTS is low. People living in cities were less likely to be infected with SFTSV as most of the patients are farmers or field workers [19, 32, 33, 51]. Even with humans-to-humans transmission, the greatest chance of occurrence was only among members of the population who come in contact with blood, such as nurses or family members [54]. In addition, secondary SFTS infection leads to milder clinical manifestations compared with primary SFTS. Therefore, secondarily infected patients only experience sudden fever, thrombocytopenia, leukopenia and gastrointestinal symptoms and recover more quickly [55].

In summary, implementing isolation and shortening the course of the disease had little impact on SFTS. This may suggested that prevention and control of the development of SFTS should focus on cutting off the direction of transmission routes. Therefore, according to the results of cutting off the transmission route earlier, interventions to prevent the transmission of SFTS suggest the following several aspects: First, for the environment, a long stay in the wild environment should be avoided. Also, residential areas should be cleaned/sanitized regularly and field working personnel should be equipped with protective clothing and gear such as insect repellents. Secondly, regarding tick transmission, regular repellent work should be done in rural and suburban areas to prevent tick breeding. Third, for animal transmission, the poultry and other domesticated animals should be well protected in rural areas to avoid the livestock being bitten by infected wild animals. At the same time, the livestock should be regularly treated with insecticide, physically examined for ticks and their serum should be tested regularly for SFTSV. Finally, for human-to-human transmission, medical staff and patients’ families should protect themselves when handling SFTS cases and avoid direct contact with the patient’s body or blood. Hospitals and township health centers in high-risk areas should be equipped with corresponding protective gear for use. Additionally, as our research results showed that isolation of SFTS patients or shortening the course of the disease cannot reduce the spread of SFTS, we should reduce the formulation of such plans to avoid wasting resources.

Limitations

There were some limitations in our study: first, the data of SFTS cases collected by the China Information Network System of Disease Prevention and Control. Since the capability and availability of this system was finite, our findings need to be interpreted with caution. Second, this study failed to separate the effects of ticks and host animals on humans in the environment. Third, some parameters of our model were collected from literature that may represent certain error in terms of timeliness.

Conclusions

In this study, the transmission dynamics model of SFTS was constructed based on the epidemiological characteristics of SFTS in Jiangsu Province. The effects of cutting off different transmission routes and taking isolation and shortening the course of the disease on SFTS control were evaluated. Finally, it was determined that cutting off the transmission route from the environment to humans had the greatest effect on SFTS prevention and control.

Additional file 1: Figure S1. The location of Jiangsu Province in China.