Skip to main content

Feasibility of controlling hepatitis E in Jiangsu Province, China: a modelling study



Hepatitis E, an acute zoonotic disease caused by the hepatitis E virus (HEV), has a relatively high burden in developing countries. The current research model on hepatitis E mainly uses experimental animal models (such as pigs, chickens, and rabbits) to explain the transmission of HEV. Few studies have developed a multi-host and multi-route transmission dynamic model (MHMRTDM) to explore the transmission feature of HEV. Hence, this study aimed to explore its transmission and evaluate the effectiveness of intervention using the dataset of Jiangsu Province.


We developed a dataset comprising all reported HEV cases in Jiangsu Province from 2005 to 2018. The MHMRTDM was developed according to the natural history of HEV cases among humans and pigs and the multi-transmission routes such as person-to-person, pig-to-person, and environment-to-person. We estimated the key parameter of the transmission using the principle of least root mean square to fit the curve of the MHMRTDM to the reported data. We developed models with single or combined countermeasures to assess the effectiveness of interventions, which include vaccination, shortening the infectious period, and cutting transmission routes. The indicator, total attack rate (TAR), was adopted to assess the effectiveness.


From 2005 to 2018, 44 923 hepatitis E cases were reported in Jiangsu Province. The model fits the data well (R2 = 0.655, P < 0.001). The incidence of the disease in Jiangsu Province and its cities peaks are around March; however, transmissibility of the disease peaks in December and January. The model showed that the most effective intervention was interrupting the pig-to-person route during the incidence trough of September, thereby reducing the TAR by 98.11%, followed by vaccination (reducing the TAR by 76.25% when the vaccination coefficient is 100%) and shortening the infectious period (reducing the TAR by 50.05% when the infectious period is shortened to 15 days).


HEV could be controlled by interrupting the pig-to-person route, shortening the infectious period, and vaccination. Among these interventions, the most effective was interrupting the pig-to-person route.

Graphic Abstract


Hepatitis E is an acute zoonotic disease caused by hepatitis E virus (HEV). The virus is mainly transmitted through the faecal-oral route upon ingestion of contaminated water and food [1]. A small percentage is transmitted through vertical and extraintestinal routes [2,3,4]. Among the four main genotypes of HEV, genotype 1 (HEV-1) and genotype 2 (HEV-2) are found in water and the external environment and only infect humans [5], whereas genotype 3 (HEV-3) and genotype 4 (HEV-4) can simultaneously infect humans and animals, such as pigs [6]. Consumption of undercooked pork is one of the main causes of infection, and the rate of HEV infection caused by direct contact with pigs is high [7,8,9,10]. A Swedish study on the positive rates of HEV in Sweden showed that the positive rates for HEV IgG were 13% and 9.3%, respectively, in pig farmers and the control group, respectively [11]. Currently, the major HEV genotype in China has changed from HEV-1 to HEV-4 [12]. Therefore, the main causes of infection in China are contact with infected pigs and humans, as well as environmental factors.

Globally, acute viral hepatitis has become a major public health problem and is now considered to be comparable with the three major infectious diseases (HIV/AIDS, malaria, and tuberculosis) [13]. According to the estimates of studies on the global burden of hepatitis E, approximately 20.1 million people are infected with HEV yearly and 3.4 million symptomatic cases are reported annually [14]. Many developing countries, including those in Asia (such as India, China, and Myanmar), Africa (Sudan, Somalia, Uganda), and North America (Mexico), show a high incidence of HEV infection [15, 16]. In China, the highest incidence is concentrated in the northwest and east regions, such as in Jiangsu Province [17, 18]. Therefore, the analysis of hepatitis E in Jiangsu Province is of great significance.

The current research model on hepatitis E mainly uses experimental animal models (such as pigs, chickens, and rabbits) to explain the transmission of HEV [19, 20]. Mathematical models [Discrete Poisson model, Auto Regressive Integrated Moving Average (ARIMA) model, and logistic model] can analyse the influencing factors of the disease and predict its developing trend [21,22,23]. A study used the susceptible-infectious-recovered (SIR) model to explain the person-to-person transmission of HEV [24]. Another study included the effect of water/food (W) to the SIR model to explore the environment-to-person route and simulate interventions for hepatitis E [25]. However, no study has developed a multi-host and multi-route transmission dynamic model (MHMRTDM) to explore the transmission features of HEV.

In this study, the MHMRTDM was used to simulate the dataset of Jiangsu Province, China. The main factors affecting the incidence and transmissibility of HEV were also investigated. Several measures were employed to evaluate the effects of interventions, including cutting off the routes of transmission, reducing the infectious period, and vaccination.


Research setting

Jiangsu Province is located on the eastern coast of China, with numerous lakes and complex water systems. It covers 107 200 km2 and administers over 13 prefecture-level cities. According to the Statistical Yearbook of Jiangsu Province, the resident population of Jiangsu Province at the end of 2019 was 80.7 million, making it the most densely populated province in China. Jiangsu Province has a transitional climate from a temperate to a subtropical zone, with mild weather, moderate rainfall, and four distinct seasons.

Data collection

This study collected information regarding hepatitis E cases from 13 cities in Jiangsu Province from 1 January 2005 to 31 December 2018. Data were provided by the Jiangsu Provincial Center for Disease Control and Prevention (CDC) and a hepatitis E epidemic dataset, containing data on almost 45 000 cases (variables included number, gender, age, and date of onset), was established. Demographic data, including the total number of people, birth rate, and natural death rate, were obtained by consulting the Jiangsu Statistical Yearbook 2019.

The MHMRTDM model without intervention

This study investigated three routes of HEV transmission: person-to-person, environment-to-person, and pig-to-person. These are the known routes through which HEV-susceptible people can be infected (Fig. 1).

Fig. 1

Study design of three transmission routes

In the pig-to-person route, healthy pigs that drink contaminated water or come into contact with infected pigs will become infected. Infected pigs excrete faeces that then contaminate water sources, forming a cycle of pig infection. When the infected pigs are slaughtered and processed for human consumption, those working in slaughter facilities or consuming undercooked pork become infected. In the environment-to-person route, contaminated water contaminates seafood, which is then consumed by people. Contaminated water may also enter the drinking water system directly or indirectly and then be consumed by people. In the person-to-person route, infected people can infect healthy people through direct contact. Infected people can become healthy again through self-healing or after treatment.

In the MHMRTDM, we adopted the subscript i and p to represent individuals and pigs, respectively. Person-to-person transmission includes the following: susceptible person (Si), exposed person (Ei), infectious person (Ii), asymptomatic person (Ai), and recovered person (Ri). Environment-to-person transmission is mainly water (W) transmission. The pig-to-person transmission route includes the following: susceptible pig (Sp), exposed pig (Ep), infectious pig (Ip), and dead pig (Dp). (Model 1 in Fig. 2). The definitions of each compartment are shown in Table 1.

Fig. 2

Establishing the transmission dynamics of the MHMRTDM model of hepatitis E

Table 1 Variables definition table of the MHMRTDM

The model was based on the following assumptions:

  • Considering the natural birth and death factors of the population and assuming the disease will not spread vertically, all new individuals born in various groups are susceptible. The population birth rate is bri, and the natural death rate is dri.

  • A susceptible person is infected through the “person-to-person”, “environment-to-person”, and “pig-to-person” routes and the infection rates are β, βw, and βp′, respectively. At the same time, it was assumed that the infection rate of an asymptomatic person is k times that of an infected person, where 0 ≤ k ≤ 1.

  • The incubation period of an exposed person (E) is 1/ω. In the model, the changed rate from Ei to Ai was assumed to be equal to that of Ii. Meanwhile, p was defined as the proportion of asymptomatic people. The rate (per day) at which an exposed person became infected was (1 − p)ωEi (0 ≤ p ≤ 1), and the rate at which they acquired an asymptomatic infection was pωEi.

  • Infected and asymptomatic people will recover after an infection period of 1/γ and 1/γ′, respectively. According to Chinese hepatitis E epidemic data released by the Ministry of Health in 2019, hepatitis E has a low case fatality rate. We therefore did not include the case fatality rate in the model for this study. The antibodies produced by a recovered person can only last for a certain period of time (antibody maintenance period, 1/f); thereafter, Ri becomes Si again.

  • An asymptomatic person (Ai), infected person (Ii), and pig (Ip) can all excrete HEV, and it is assumed that the coefficient of HEV excretion is μ, μ′, and μ″, respectively. Previous studies showed that infected people present with diarrhea, whereas infected pigs remain asymptomatic. We assumed that the ability to detoxify while presenting symptoms is twice as great as the ability to detoxify while remaining asymptoms; thus, μ = ′ = 2μ″.

  • The HEV will cease to exist after a certain period in the environment, and the survival time is 1/ε.

  • In this study, we assigned the pig birth rate as brp and the pig natural mortality rate as drp.

  • A susceptible pig is infected through "pig-to-pig" or "environment-to-pig" routes, and the infection rates are βp and βw′, respectively.

  • The incubation period of the exposed pig (Ep) is 1/ωp, and the time from pig infection to slaughter is 1/\(\tau\). Most people are infected with HEV after ingesting undercooked pork. A slaughtered infected pig becomes meat (pork) after a certain period of time depending on the rate of slaughtering.

  • It was assumed that the number of pigs in the Jiangsu Province market is a dynamic balance, and the birth rate and death rates of pigs are the equal. Owing to infection, disease, recovery, and other reasons, several groups of people and pigs change dynamically at a certain speed, and the virus in the environment also changes dynamically over time.

The equations of the model containing the three transmission routes are as follows:

$$\frac{d{S}_{i}}{dt}=b{r}_{i}{N}_{i}+f{R}_{i}-{dr}_{i}{S}_{i}-{\beta S}_{i}\left({I}_{i}+{kA}_{i}\right)-{\beta }_{w}{S}_{i}W-{\beta }_{p}^{{\text{'}}}{S}_{i}{D}_{p}$$
$$\frac{d{E}_{i}}{dt}=\beta {S}_{i}\left({I}_{i}+k{A}_{i}\right)+{\beta }_{w}{S}_{i}W+{\beta }_{p}^{\text{'}}{S}_{i}{D}_{p}-(1-p)\omega {E}_{i}-p{\omega }^{\text{'}}{E}_{i}-{dr}_{i}{E}_{i}$$
$$\frac{d{I}_{i}}{dt}=(1-p)\omega {E}_{i}-(\gamma +{dr}_{i}){I}_{i}$$
$$\frac{d{A}_{i}}{dt}=p{\omega }^{{\text{'}}}{E}_{i}-({\gamma }^{{\text{'}}}+{dr}_{i}){A}_{i}$$
$$\frac{d{R}_{i}}{dt}=\gamma {I}_{i}+{\gamma }^{{{'}}}{A}_{i}-(f+{dr}_{i}){R}_{i}$$
$$\frac{dW}{dt}=\mu {I}_{i}+{\mu}^{{{'}}}{A}_{i}+{\mu }^{{{''}}}{I}_{p}-\varepsilon W$$
$$\frac{d{S}_{p}}{dt}=b{r}_{p}{N}_{p}-{\beta }_{w}^{{\text{'}}}{S}_{p}W-{\beta }_{p}{I}_{p}{S}_{p}$$
$$\frac{d{E}_{p}}{dt}={\beta }_{w}^{{\text{'}}}{S}_{p}W+{\beta }_{p}{I}_{p}{S}_{p}-{\omega }_{p}{E}_{p}$$
$$\frac{d{I}_{p}}{dt}={\omega }_{p}{E}_{p}-{\tau I}_{p}$$
$$\frac{d{D}_{p}}{dt}={\tau I}_{p}-{dr}_{p}{D}_{p}$$

To consistently and easily calculate the dimensions across all groups of people and modes and reduce the number of parameters, the variables and parameters of the model were normalised, so that si = Si/Ni, ei = Ei/Ni, ii = Ii/Ni, ai = Ai/Ni, ri = Ri/Ni, w = εW/(μN + μ″Np), b = βNi, bW = μβW Ni/ε, μ = μ′, sp = Sp/Np, ep = Ep/Np, ip = Ip/Np, dp = Dp/Np, x = Np/Ni, σ = bw′w + bpip, bp = βpNp, bw = μ″βW′Np/ε.

The normalised equations are as follows:

$$\frac{{de}_{i}}{dt}=b{s}_{i}\left({i}_{i}+k{a}_{i}\right)+{b}_{w}{s}_{i}w+{b}_{p}^{{\text{'}}}{s}_{i}{d}_{p}-(1-p)\omega {e}_{i}-p{\omega }^{{\text{'}}}{e}_{i}-{dr}_{i}{e}_{i}$$
$$\frac{{di}_{i}}{dt}=(1-p)\omega {e}_{i}-\gamma {i}_{i}-{dr}_{i}{i}_{i}$$
$$\frac{{da}_{i}}{dt}=p{\omega }^{{\text{'}}}{e}_{i}-{\gamma }^{{\text{'}}}{a}_{i}-{dr}_{i}{a}_{i}$$
$$\frac{{dr_{i} }}{{dt}} = \gamma ^{\prime}a_{i} + ~\gamma i_{i} - fr_{i} - dr_{i} r_{i}$$
$$\frac{dw}{dt}=\varepsilon ({i}_{i}+{a}_{i}+2{i}_{p}-w)$$
$$\frac{{ds}_{p}}{dt}=x{br}_{p}-{\sigma s}_{p}$$
$$\frac{{de}_{p}}{dt}=\sigma {s}_{p}-{\omega }_{p}{e}_{p}$$
$$\frac{{di}_{p}}{dt}={\omega }_{p}{e}_{p}-\tau {i}_{p}$$
$$\frac{{dd}_{p}}{dt}=\tau {i}_{p}-{dr}_{p}{d}_{p}$$

We built a new model (Model 2) to estimate the vaccination effects (Fig. S1 in Additional file 2) by adding the following assumptions to Model 1:

  • Since the hepatitis E vaccine in China is not compulsory, we assumed that the vaccination coefficient of the vaccine is δ. V and G may also be infected by the people, environment, and pigs, and the infection rates are β, βw, and βp′, respectively.

  • The time from the absence of antibody to the production of antibody after vaccination was assumed to be 1/φ, and the protection rate of the vaccine was assumed to be 1-λ.

The equations of Model 2 are as follows:

$$\frac{d{S}_{i}}{dt}=b{r}_{i}{N}_{i}+f{R}_{i}-{dr}_{i}{S}_{i}-{\beta S}_{i}\left({I}_{i}+{kA}_{i}\right)-{\beta }_{w}{S}_{i}W-{\beta }_{p}^{{\text{'}}}{S}_{i}{D}_{p}-\delta {S}_{i}$$
$$\frac{d{E}_{i}}{dt}=\beta {S}_{i}\left({I}_{i}+k{A}_{i}\right)+{\beta }_{w}{S}_{i}W+{\beta }_{p}^{{\text{'}}}{S}_{i}{D}_{p}+{\beta }_{p}V{D}_{p}+\beta V\left({I}_{i}+k{A}_{i}\right)+{\beta }_{w}VW+(1-\lambda ){(\beta }_{p}G{D}_{p}+\beta G\left({I}_{i}+k{A}_{i}\right)+{\beta }_{w}GW)-(1-p)\omega {E}_{i}-p{\omega }^{{\text{'}}}{E}_{i}-{dr}_{i}{E}_{i}$$
$$\frac{d{I}_{i}}{dt}=(1-p)\omega {E}_{i}-(\gamma +{dr}_{i}){I}_{i}$$
$$\frac{d{A}_{i}}{dt}=p{\omega }^{{\text{'}}}{E}_{i}-({\gamma }^{{\text{'}}}+{dr}_{i}){A}_{i}$$
$$\frac{d{R}_{i}}{dt}=\gamma {I}_{i}+{\gamma }^{{\text{'}}}{A}_{i}-(f+dr){R}_{i}$$
$$\frac{dW}{dt}=\mu {I}_{i}+{\mu }^{{{'}}}{A}_{i}+{\mu }^{{{''}}}{I}_{p}-\varepsilon W$$
$$\frac{dV}{dt}=\delta {S}_{i}-{\beta }_{p}V{D}_{p}-\beta V\left({I}_{i}+k{A}_{i}\right)-{\beta }_{w}VW-{dr}_{i}V-\varphi V$$
$$\frac{{dG}}{{dt}} = \varphi V - dr_{i} G - \left( {1 - \lambda } \right)~\left( {\beta _{p} GD_{p} + \beta G\left( {I_{i} + kA_{i} } \right) + \beta _{w} GW} \right)$$
$$\frac{d{S}_{p}}{dt}=b{r}_{p}{N}_{p}-{\beta }_{w}^{{{'}}}{S}_{p}W-{\beta }_{p}{I}_{p}{S}_{p}$$
$$\frac{d{E}_{p}}{dt}={\beta }_{w}^{{\text{'}}}{S}_{p}W+{\beta }_{p}{I}_{p}{S}_{p}-{\omega }_{p}{E}_{p}$$
$$\frac{d{I}_{p}}{dt}={\omega }_{p}{E}_{p}-{\tau I}_{p}$$
$$\frac{d{D}_{p}}{dt}={\tau I}_{p}-{dr}_{p}{D}_{p}$$

Simulating the effectiveness of interventions

This study simulates three interventions: cutting off transmission routes, shortening the infectious period, and vaccination. We forecasted the incidence of infection for 2019–2023 based on the reported incidence data in Jiangsu Province during 2005–2018 and started to implement the intervention forecast for the 2018–2023 period. In addition, to investigate the difference in the effect of interventions at provincial and municipal levels, we selected three cities in Jiangsu Province with high, moderate, and low annual average incidence rates for the same operation: Zhenjiang City, Yancheng City, and Wuxi City. We compared the differences under the following scenarios:

  1. (a)

    Scenario 1: Cutting off the transmission routes

    1. (1)

      If all patients were isolated and we set the isolation time to March, June, September, and December 2018, then b = 0.

    2. (2)

      If we blocked the route of external transmission (contaminated food and water resources) and set the intervention time points to March, June, September, and December 2018, we only set bw = 0.

    3. (3)

      If the pig-to-person transmission was cut off and we set the cut off time points to March, June, September, and December 2018. At the same time, we set bp = 0.

  2. (b)

    Scenario 2: Shortening the infectious period

When an infected patient receives treatment promptly after symptoms appear, the time of treatment is shortened. The recovery period 1/γ and 1/γ′ will also be shortened. In this study, the recovery/infectious period was shortened from 30 days to 27, 24, 21, 18, and 15 days, respectively. We adjusted γ and γ′ after December 2018 to 1.111, 1.25, 1.429, 1.667, and 2. The other parameter settings remained unchanged.

  1. (iii)

    Scenario 3: Vaccinating susceptible people

The vaccine intervention model has two more warehouses than the original model; we used V and G to denote those who were vaccinated but not resistant to infectious diseases and those who are resistant after vaccination, respectively. We then set the monthly vaccination ratio of the vaccine using a monthly vaccination coefficient, δ. We set δ to 20%, 40%, 60%, 80%, and 100%, respectively. All other parameters remained unchanged.

Parameter estimation

The parameters of the population in the MHMRTDM used in this study were bri, dri, b, bW, p, ωi, γ, γ′, ε, and f (Table 2). According to the 2006 Statistical Yearbook of each city in Jiangsu Province, this study used month as the basic time unit of the birth rate and natural death rate; thus, bri and dri were taken as 1/12 of each city’s annual birth rate and death rates.

Table 2 Parameter definitions and values

According to previous studies [25], during an hepatitis E epidemic, the proportion of asymptomatic infections is p = 0.9, the reciprocal of the average disease course is γ = 1, and the reciprocal of the incubation period is ωi = 0.75.

An individual’s HEV IgG “antibody maintenance time” after infection with HEV exceeds 14 years [29]; thus, f = 1/(14 × 12) = 0.00595. Based on the literature, the survival time of HEV in water and the environment is 1/ε, about 3 weeks; thus, thus we set ε = 1.3333 [28].

The parameters of the pig groups are brp, drp, bp, ωp, and τ. According to the literature [23], the time from birth to slaughter of pigs is approximately 5 months; this was used to estimate the monthly birth rate and natural mortality of pigs, namely brp = 0.2 and drp = 0.2. The time from birth to infection of pigs is approximately 2 months, and the monthly slaughter rate of infected pigs is τ = 0.333. According to the literature, the expected value of pig HEV infectivity is between 2.68 × 10–2 and 3.45 × 10–2 per day [26]. In the previously mentioned study, the rate of infection of pigs was 3 × 10–2 per day, which corresponds to 0.9 per month (bp = 0.9/month). According to the "China Rural Statistical Yearbook 2019" and the main data of the sixth national census in 2010, the ratio of the total number of pigs to the total population density is x = 0.333.

The curve fitting time of this study began in 2005. For models with intervention measures, the values of the parameters at the beginning of the intervention are the simulated output values of the model at the corresponding time without intervention. The simulated time for vaccination was 2019–2023, and we simulated with inoculation coefficients of five scenarios, which were set to 20%, 40%, 60%, 80%, and 100%.

Seasonality estimation

Hepatitis E has obvious seasonality and periodicity; thus, it was necessary to correct the model for these factors. Based on a previous study [31], the seasonal and periodic correction results of the model can be obtained using cosine function correction \(\beta\).

$$\beta \left( t \right) = \beta _{0} \left[ {{{1~ + ~\rho \cos \left( {~2\pi t} \right)} \mathord{\left/ {\vphantom {{1~ + ~\rho \cos \left( {~2\pi t} \right)} T}} \right. \kern-\nulldelimiterspace} T}} \right]$$

where β(t) represents β at t time and β0 represents β at the initial time of simulation, with a period of 1 year. In this study, month is the basic time unit of hepatitis E incidence; thus, it is necessary to introduce the parameter seasonal cycle T = 12. It was also observed that the incidence of hepatitis E in the non-epidemic season in Jiangsu Province was close to zero, indicating that it was the lowest epidemic month of the disease, that is, when T = 6 k (k is odd), β (t) = β(0) (1-ρ) = 0, and c = 1. Because the sine function can be regarded as the translation of the cosine function, the sine expression was used in this study. Considering that the monthly incidence period of hepatitis E in Jiangsu Province is different from that in Changsha, we introduced the new parameter Δt, under which the correction effect of Δt can be divided in two parts, that is, the conversion from cosine function to sine function and the seasonal difference between Changsha and Jiangsu Provinces. The equation used in this study was as follows:

$$b~ = b_{0} \left[ {1~ + ~\sin \left( {\frac{{2\pi \left( {t + \Delta t} \right)}}{{12}}} \right)} \right]$$
$$b_{w} ~ = b_{{w0}} \left[ {1~ + ~\sin \left( {\frac{{2\pi \left( {t + \Delta t} \right)}}{{12}}} \right)} \right]$$

Among them, Δt is set and adjusted in sections according to the principle of best fit.

Evaluating the effectiveness of interventions

Measures to assess the effectiveness of interventions were used in this study with a total attack rate (TAR). The formula is as follows:

$$TAR = \frac{{TN}}{N} \times 100\%$$

In the equation, TN and N refer to the total number of new cases and the total population number, respectively.

Simulation method and statistical analysis

In this study, Microsoft Excel (2020 version; Microsoft Corp., Redmond, WA) was used for the entry and management of related data and related mapping, and Berkeley Madonna 8.3.18 software (developed by Robert Macey and George Oster of the University of California at Berkeley, CA, USA) was used for modelling, while the Runge–Kutta method of order 4 with tolerance set at 0.001 was used to perform curve fitting. When curve fit is in progress, Berkeley Madonna displays the root mean square deviation between the data and best run so far [32, 33]. The IBM SPSS Statistics for Windows, (version 23.0; JBM Corp., Armonk, NY) curve estimation function was used to obtain R2 to describe the goodness of fit.

Sensitivity analysis

Sensitivity analysis was conducted for 14 parameters, which were divided into 42 values according to their range. In this study, sensitivity analysis was performed for the 2004–2018 date of Jiangsu Province. In general, the TAR of infectious diseases has different rise and decline trends. Each parameter was divided into three segments and calibrated with a minimum, moderate, and maximum or greater segment when it did not have a maximum.


Epidemiological characteristics

A total of 44 923 cases of hepatitis E were reported in Jiangsu Province, China, from 2005 to 2018 (average annual incidence of 4.12 per 100 000 people). The annual incidence of hepatitis E in Jiangsu Province peaked in 2007 (4.35 per 100 000 people), 2011 (5.21 per 100 000 people), and 2013 (4.94 per 100 000 people). The lowest incidence was 2.96 per 100 000 people in 2005. The three cities with the highest average annual incidence in Jiangsu Province (Fig. 3) were Zhenjiang City (7.90 per 100 000 people), Nantong City (6.96 per 100 000 people), and Yangzhou City (5.72 per 100 000 people). The three cities with the lowest average annual incidence were Wuxi City (1.32 per 100 000 people), Suzhou City (2.24 per 100 000 people), and Changzhou City (2.32 per 100 000 people). The annual incidence rates in Lianyungang City (4.37 per 100 000 people) and Taizhou City (4.06 per 100 000 people) are increasing each year, while the annual incidence rates in Changzhou City (2.32 per 100 000 people) and Nanjing City (3.07 per 100 000 people) are decreasing each year. The overall annual incidence rates in other regions showed an upward trend first, followed by a downward trend. The peak incidence rate of hepatitis E in Jiangsu Province and other cities was around March; however, the peak transmissibility was in December and January (Table 3). The results of curve fitting showed that the MHMRTDM fits the data well (Fig. 4). Jiangsu Province (R2 = 0.655, P < 0.001) and the first three cities that fit well are Huai’an City (R2 = 0.632, P < 0.001), Nantong City (R2 = 0.619, P < 0.001), and Changzhou City (R2 = 0.596, P < 0.001).

Fig. 3

Incidence map of each city in Jiangsu Province from 2005 to 2018

Table 3 Peak time and transmission capacity of different transmission routes of hepatitis E in Jiangsu Province
Fig. 4

Simulated situation of hepatitis E incidence in Jiangsu Province and other cities during 2005–2018. A Jiangsu Province; B Suqian City; C Huaian City; D Changzhou City; E Nantong City; F Yangzhou City; G Lianyungang City; H Nanjing City; I Wuxi City; J Xuzhou City; K Suzhou City; L Taizhou City; M Yancheng City; N Zhenjiang City

Effectiveness of cutting off transmission routes

As shown in Fig. 5, the results of cutting off the person-to-person and environment-to-person routes in Jiangsu Province are not obvious regarding prevention and control. Cutting off the pig-to-person route, however, can effectively control the spread of the diseases. If the pig-to-person transmission route was cut off in March 2018, the TAR from March 2018 to December 2023 was predicted to be 1.398 per 100 000 people, which was lower than the predicted TAR of 17.437 per 100 000 people (reduction of 92.0%). If the transmission route was cut off in June 2018, the TAR from June 2018 to December 2023 was 0.504 per 100 000 people, which is lower than the predicted TAR of 16.310 per 100 000 people (reduction of 96.9%). If the transmission route was cut off in September 2018, the TAR from September 2018 to December 2023 was 0.301 per 100 000 people, which was lower than the predicted TAR of 15.895 per 100 000 people (reduction of 98.1%). If the transmission route was cut off in December 2018, the TAR between December 2018 and December 2023 was 1.109 per 100 000 people, which is lower than the predicted TAR of 15.508 per 100 000 people (reduction of 92.9%). The effectiveness of cutting off transmission routes in Zhenjiang City, Yancheng City, and Wuxi City was shown are shown in Additional file 1: Table S1 and Additional file 2: Figs. S2–S4.

Fig. 5

Effects of cutting the transmission routes on the incidence of hepatitis E

Effectiveness of shortening the infectious period

The incidence of hepatitis E in Jiangsu Province was reduced when the infectious period was shortened (Fig. 6). When the predicted infection period was 27, 24, 21, 18, and 15 days without intervention from 2019 to 2023, the TAR values (and reduction proportions) per 100 000 people were 13.721 (9.8%), 12.195 (19.9%), 10.666 (29.9%), 9.130 (40.0%), and 7.601 (50.1%), respectively. The effectiveness of shortening the infectious period in Zhenjiang City, Yancheng City, and Wuxi City are shown in the Additional file 1: Table S1 and Additional file 2: Figs. S5–S7.

Fig. 6

Simulation of shortening the incidence of hepatitis E infection

Effectiveness of assessing the effectiveness of vaccination

As shown in Fig. 7, the simulation of general population immunisation in Jiangsu Province decreased the TAR after vaccination and the vaccine effect increased with the vaccination coefficient. Within five years of the simulated 2019–2023 period, when the vaccination coefficient was 0.2, the TAR was 4.496/100 000, showing a decrease of 70.5% compared with that without intervention. Similarly, vaccination coefficients of 0.4, 0.6, 0.8, and 1 yielded TAR values per 100 000 people (and decreasing proportions) of 3.932 (74.2%), 3.754 (75.3%), 3.667 (75.9%), and 3.614 (76.3%), respectively. The effectiveness of vaccination in Zhenjiang City, Yancheng City, and Wuxi City are shown in the Additional file 1: Table S1 and Additional file 2: Figs. S8–S10.

Fig. 7

Effects of different vaccination coefficients on the incidence of hepatitis E infection

Sensitivity analysis

In this study, all parameters executed in the model are consistent with the range of reference values. The MHMRTDM was sensitive to p, γ, and x, but not sensitive to other parameters. (Additional file 2: Fig. S11).


The World Health Organization has formulated a plan called “Hepatitis Elimination in 2030”, but owing to its complex transmission, only few mathematical model studies have focused on HEV transmission, such as the SIR model and the Susceptible-Exposed-Symptomatic-Asymptomatic-Recovered/Removed-Water (SEIARW) model [24, 25]. Therefore, this study first clarified the transmission mechanism of how the virus transmits through multiple hosts and the environment ("person-to-person", "environment-to-person", and "pig-to-person") and then explored the influences of meteorological factors on the disease incidence. These factors are important for the prevention and control of hepatitis E.

According to the determination coefficient (R2) of the goodness of fit test, the model fits the reported incidence in different cities in Jiangsu Province well, indicating that the model is suitable for this study. At the same time, the effectiveness of the model fitting is consistent with the multi-pathway HEV SEIARW transmission model and its application in Changsha city [25].

In recent years, the annual incidence of hepatitis E in China has been rising. From 2004 to 2017, the annual incidence of hepatitis E in China increased from 1.27/100 000 to 2.10/100 000[34]. The average annual incidence of hepatitis E in Jiangsu Province from 2005 to 2018 was 4.12/100 000, which was significantly higher than the national level. However, from 2013 to 2018, the incidence of hepatitis E in Jiangsu Province displayed a downward trend, possibly owing to the strengthened prevention and control measures of the local CDC or due to vaccination. In 2011, the world’s first hepatitis E vaccine (HEV 239, Xiamen Innovax Biotech, Xiamen, China) was licensed in China. Volunteers were recruited in Jiangsu Province for the vaccine study, which was a large clinical test involving 120 000 participants [35].

At the same time, our results indicated that the incidence of hepatitis E varies from region to region. According to the Statistical Yearbook of Jiangsu Province, the gross domestic product rankings of Zhenjiang City and Yangzhou City are relatively low, while those of Suzhou City, Wuxi City, and Changzhou City rankings are high Economic development is likely to be an important reason for the uneven distribution of disease burden in Jiangsu Province. Additionally, climatic environment, population, urban environment, and human factors may represent additional reasons for the uneven disease burden distribution in Jiangsu Province.

The peak time of the disease in Jiangsu Province and its cities is in February and March each year, which is consistent with a previous study [36]. This period coincides with the Chinese Lunar New Year [22]. Perhaps because there are more holidays during the Chinese New Year than usual and people eat out and socialise more, the number of people going out for social activities increases. The greater demand for food will therefore also increase the purchase of poultry, livestock (such as pigs), and seafood; thus, these two factors probably increase transmission rates. Concurrently, the months of February and March are relatively cold, which may lower people’s immunity and consequently increase their susceptibility to the virus. Therefore, health education regarding food quality control and personal hygiene should be strengthened during this period. However, we found that transmissibility usually peaked in December and January, 3 months before the peak incidence of hepatitis E; this also coincides with the time pigs grow up and become infected with the virus, which can then be transmitted to people.

The peak of the disease in each city of Jiangsu Province is mainly in March, while the peak in August is very small (not included when fitting the model). However, the peak of rainfall is in August, which is consistent with the second incidence peak of the incidence, indicating that the incidence of the disease in Jiangsu Province is mainly related to the pig-to-person transmission route, but not closely related to the environment-to-person transmission route.

According to the model’s prediction, without intervention, the incidence will show continual cyclical changes and the epidemic will continue to develop without a spontaneous end. Intervention effects at provincial and municipal levels in Jiangsu Province were very similar. When the person-to-person and environment-to-person transmission routes were cut off, the incidence hardly decreased; conversely, cutting off the pig-to-person transmission route has achieved significant results. This indicates that controlling the contact between pigs and contaminated water as well as the quality of pork and corresponding meat products in the market can effectively control the transmission of HEV. Residents should ensure that they only consume pork meat that has been properly and thoroughly cooked, and this precaution can result in improved eating habits of the residents. The gross domestic product ranking of Jiangsu Province is among the highest in China (Statistical Yearbook of China 2019), and many related studies have indicated that HEV is mainly contracted from food in areas with better economic development [7,8,9].

Furthermore, the results showed that shortening the period of infection can lower the incidence of HEV, representing another effective intervention and suggesting that we should aim for early detection, early diagnosis, and early treatment for controlling HEV transmission. To achieve this, it is necessary to improve the detection rate and treatment of hepatitis E. In addition, a study has shown that the hepatitis E vaccine has a good effect on controlling the disease [29], and our study also showed that the higher the vaccination coefficient, the more the incidence rate decreases. When the vaccination coefficient was 60%, the incidence has already decreased even more significantly. When the vaccination rate reached 80%, the disease was effectively under control. However, this level of vaccination coefficient is very difficult to achieve because (1) only those 16 years and older are eligible for vaccination. However, according to the Statistical Yearbook of Jiangsu Province in 2018, 14.7% of the population was aged 0–14 years; and (2) the vaccination is voluntary and is self-limited for most people; thus, the vaccination rate will likely remain low. However, given that hepatitis E has a mortality rate of 1–3% in people aged 14–40 years [37] and up to 30% in pregnant women [38], and an even worse effect on those with underlying hepatitis [39], vaccination is necessary.

This study has some limitations. We used secondary information, mainly from the Statistical Yearbook of Jiangsu Province, and the key parameters p, γ, and x used in our model were derived from references instead of first-hand data. This may have influenced the accuracy of the model. The model should be further adjusted according to specific parameters for each city. In addition, further studies that model a larger variety of hepatitis E interventions and explore different variable factors, such as meteorological data, would be useful.


In Jiangsu Province, the main route of HEV transmission is from pigs to humans in Jiangsu Province. Effective prevention measures for hepatitis E are needed to control the transmission from pigs to humans, strengthen the management of contact between pigs and contaminated water, and improve the quality control of market pork and related meat products. In addition, effective treatment and hepatitis E vaccination can also effectively prevent the spread of this disease.

Availability of data and materials

Data supporting the conclusions of this article are included within the article.



Hepatitis E virus


Genotype 1 of hepatitis E virus


Genotype 2 of hepatitis E virus


Genotype 3 of hepatitis E virus


Genotype 4 of hepatitis E virus


Human immunodeficiency virus


Acquired immune deficiency syndrome


Autoregressive integrated moving average (model)


Susceptible–infectious–recovered (model)


Multi-host and multi-route transmission dynamic model


Immunoglobin G


Total attack rate

R 2 :

Coefficient of determination


Susceptible–exposed–symptomatic–asymptomatic–recovered/removed–water (model)


  1. 1.

    Schlosser J, Eiden M, Vina-Rodriguez A, Fast C, Dremsek P, Lange E, et al. Natural and experimental hepatitis E virus genotype 3-infection in European wild boar is transmissible to domestic pigs. Vet Res. 2014;45:121.

    Article  Google Scholar 

  2. 2.

    Khuroo MS, Kamili S, Jameel S. Vertical transmission of hepatitis E virus. Lancet. 1995;345:1025–6.

    CAS  Article  Google Scholar 

  3. 3.

    Kumar RM, Uduman S, Rana S, Kochiyil JK, Usmani A, Thomas L. Sero-prevalence and mother-to-infant transmission of hepatitis E virus among pregnant women in the United Arab Emirates. Eur J Obstet Gynecol Reprod Biol. 2001;100:9–15.

    CAS  Article  Google Scholar 

  4. 4.

    Arankalle VA, Chobe LP. Hepatitis E virus: can it be transmitted parenterally? J Viral Hepat. 1999;6:161–4.

    CAS  Article  Google Scholar 

  5. 5.

    Kamar N, Bendall R, Legrand-Abravanel F, Xia NS, Ijaz S, Izopet J, et al. Hepatitis E. Lancet. 2012;379:2477–88.

    Article  Google Scholar 

  6. 6.

    Melgaço JG, Gardinali NR, de Mello VDM, Leal M, Lewis-Ximenez LL, Pinto MA. Hepatitis E: update on prevention and control. Biomed Res Int. 2018;2018:5769201.

    PubMed  PubMed Central  Google Scholar 

  7. 7.

    Chan MCW, Kwok K, Hung TN, Chan PKS. Molecular epidemiology and strain comparison between hepatitis E viruses in human sera and pig livers during 2014 to 2016 in Hong Kong. J Clin Microbiol. 2017;55:1408–15.

    Article  Google Scholar 

  8. 8.

    Colson P, Borentain P, Queyriaux B, Kaba M, Moal V, Gallian P, et al. Pig liver sausage as a source of hepatitis E virus transmission to humans. J Infect Dis. 2010;202:825–34.

    Article  Google Scholar 

  9. 9.

    Matsuda H, Okada K, Takahashi K, Mishiro S. Severe hepatitis E virus infection after ingestion of uncooked liver from a wild boar. J Infect Dis. 2003;188:944.

    Article  Google Scholar 

  10. 10.

    Krumbholz A, Mohn U, Lange J, Motz M, Wenzel JJ, Jilg W, et al. Prevalence of hepatitis E virus-specific antibodies in humans with occupational exposure to pigs. Med Microbiol Immunol. 2012;201:239–44.

    CAS  Article  Google Scholar 

  11. 11.

    Olsen B, Axelsson-Olsson D, Thelin A, Weiland O. Unexpected high prevalence of IgG-antibodies to hepatitis E virus in Swedish pig farmers and controls. Scand J Infect Dis. 2006;38:55–8.

    Article  Google Scholar 

  12. 12.

    Chen C, Zhang SY, Zhang DD, Li XY, Zhang YL, Li WX, et al. Clinical features of acute hepatitis E super-infections on chronic hepatitis B. World J Gastroenterol. 2016;22:10388–97.

    CAS  Article  Google Scholar 

  13. 13.

    Locarnini S, Chen DS, Shibuya K. No more excuses: viral hepatitis can be eliminated. Lancet. 2016;387:1703–4.

    Article  Google Scholar 

  14. 14.

    WHO. Global Hepatitis Report. World Health Organization. 2017. Accessed 1 Mar 2021.

  15. 15.

    Khuroo MS, Khuroo MS. Hepatitis E: an emerging global disease—from discovery towards control and cure. J Viral Hepat. 2016;23:68–79.

    CAS  Article  Google Scholar 

  16. 16.

    Fierro NA, Realpe M, Meraz-Medina T, Roman S, Panduro A. Hepatitis E virus: an ancient hidden enemy in Latin America. World J Gastroenterol. 2016;22:2271–83.

    CAS  Article  Google Scholar 

  17. 17.

    Huang RT, Li DR, Wei J, Huang XR, Yuan XT, Tian X. Isolation and identification of hepatitis E virus in Xinjiang, China. J Gen Virol. 1992;73:1143–8.

    Article  Google Scholar 

  18. 18.

    Yang S, Wu J, Ding C, Cui Y, Zhou Y, Li Y, et al. Epidemiological features of and changes in incidence of infectious diseases in China in the first decade after the SARS outbreak: an observational trend study. Lancet Infect Dis. 2017;17:716–25.

    Article  Google Scholar 

  19. 19.

    Yugo DM, Cossaboom CM, Meng XJ. Naturally occurring animal models of human hepatitis E virus infection. Ilar J. 2014;55:187–99.

    CAS  Article  Google Scholar 

  20. 20.

    Corneillie L, Banda DH, Meuleman P. Animal models for hepatitis E virus. Viruses. 2019;11:564.

    CAS  Article  Google Scholar 

  21. 21.

    Ren H, Li J, Yuan ZA, Hu JY, Yu Y, Lu YH. The development of a combined mathematical model to forecast the incidence of hepatitis E in Shanghai, China. BMC Infect Dis. 2013;13:421.

    Article  Google Scholar 

  22. 22.

    Liu K, Cai J, Wang S, Wu Z, Li L, Jiang T, et al. Identification of distribution characteristics and epidemic trends of hepatitis E in Zhejiang Province, China from 2007 to 2012. Sci Rep. 2016;6:25407.

    CAS  Article  Google Scholar 

  23. 23.

    Satou K, Nishiura H. Transmission dynamics of hepatitis E among swine: potential impact upon human infection. BMC Vet Res. 2007;3:9.

    Article  Google Scholar 

  24. 24.

    Backer JA, Berto A, McCreary C, Martelli F, van der Poel WH. Transmission dynamics of hepatitis E virus in pigs: estimation from field data and effect of vaccination. Epidemics. 2012;4:86–92.

    CAS  Article  Google Scholar 

  25. 25.

    Tianmu C, Qiong H, Aichun T. Development of a multiple transmission pathways dynamic model of hepatitis E and its application dynamics in Changsha. Chin J Health Stat. 2014;2:20.

    Google Scholar 

  26. 26.

    Satou K, Nishiura H. Transmission dynamics of hepatitis E among swine: potential impact upon human infection. BMC Vet Res. 2007;3:9.

  27. 27.

    Zhuang H, Cao XY, Liu CB, Wang GM. Epidemiology of hepatitis E in China. Gastroenterol Jpn. 1991;26(Suppl 3):135–8.

    Article  Google Scholar 

  28. 28.

    Parashar D, Khalkar P, Arankalle VA. Survival of hepatitis A and E viruses in soil samples. Clin Microbiol Infect. 2011;17:E1-4.

    CAS  Article  Google Scholar 

  29. 29.

    Zhang J, Zhao Q, Xia N. Prophylactic hepatitis E vaccine. Adv Exp Med Biol. 2016;948:223–46.

    Article  Google Scholar 

  30. 30.

    Zhang J, Zhang XF, Huang SJ, Wu T, Hu YM, Wang ZZ, et al. Long-term efficacy of a hepatitis E vaccine. N Engl J Med. 2015;372:914–22.

    CAS  Article  Google Scholar 

  31. 31.

    Su YY, Huang SJ, Guo M, Zhao J, Yu H, He WG, et al. Persistence of antibodies acquired by natural hepatitis E virus infection and effects of vaccination. Clin Microbiol Infect. 2017;23(336):e331–4.

    Google Scholar 

  32. 32.

    Zhao ZY, Zhu YZ, Xu JW, Hu SX, Hu QQ, Lei Z, et al. A five-compartment model of age-specific transmissibility of SARS-CoV-2. Infect Dis Poverty. 2020;9:117.

    Article  Google Scholar 

  33. 33.

    Zhao ZY, Chen Q, Zhao B, Hannah MN, Wang N, Wang YX, et al. Relative transmissibility of shigellosis among male and female individuals: a modeling study in Hubei Province. China. Infect Dis Poverty. 2020;9:39.

    Article  Google Scholar 

  34. 34.

    Sun XJ, Zhang GM, Zheng H, Miao N, Wang HQ, Yin ZD, et al. Epidemiological analysis of viral hepatitis E in China, 2004–2017. Zhonghua Yu Fang Yi Xue Za Zhi. 2019;53:382–7.

    CAS  PubMed  Google Scholar 

  35. 35.

    Jiangsu Net, China. The world's first HEV vaccine was born in China. 2019. Accessed 1 Mar 2021.

  36. 36.

    Tan J, Chen Y, Wang L, Chan TC, Amer S, Xu X, et al. Acute sporadic hepatitis E in the Zhejiang coastal area of China: a 14-year hospital-based surveillance study. Virol J. 2019;16:16.

    Article  Google Scholar 

  37. 37.

    Khuroo MS. Discovery of hepatitis E: the epidemic non-A, non-B hepatitis 30 years down the memory lane. Virus Res. 2011;161:3–14.

    CAS  Article  Google Scholar 

  38. 38.

    Khuroo MS, Kamili S. Aetiology, clinical course and outcome of sporadic acute viral hepatitis in pregnancy. J Viral Hepat. 2003;10:61–9.

    CAS  Article  Google Scholar 

  39. 39.

    Hamid SS, Atiq M, Shehzad F, Yasmeen A, Nissa T, Salam A, et al. Hepatitis E virus superinfection in patients with chronic liver disease. Hepatology. 2002;36:474–8.

    Article  Google Scholar 

Download references


The authors thank all study participants for providing the data and field investigators for collecting the data.


This study was partly supported by the Bill & Melinda Gates Foundation (INV-005834), the Science and Technology Program of Fujian Province (No: 2020Y0002), the Xiamen New Coronavirus Prevention and Control Emergency Tackling Special Topic Program (No: 3502Z2020YJ03).

Author information




Conceptualisation: T-mC, MY, X-qC, Z-yZ, and P-hL; Methodology: T-mC, MY, Z-yZ, X-qC, and P-hL, JR, S-nL, X-fZ, J-wX and Y-zZ; Data analysis: W-dL, YH, QY, S-yM and Y-hS; Writing—original draft: T-mC, MY, Z-yZ, P-hL, YW, X-cL, LL, BD, CL and J-fH; Writing—review & editing: T-lY, Z-yL, W-kL and B-hZ; Data collection: X-fZ, W-dL. All authors have read and agreed to the final version of the manuscript.

Corresponding authors

Correspondence to Yan-Hua Su, Xue-Feng Zhang or Tian-Mu Chen.

Ethics declarations

Ethics approval and consent to participate

The effort towards disease control is part of the Centre for Disease Control and Prevention’s routine responsibility in Jiangsu Province, China. Therefore, institutional review and informed consent were not required for this study. All data analysed were anonymised.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Supplementary Information

Additional file 1: Table S1.

Evaluates the effectiveness of intervention in 3 cities of Jiangsu Province, China.

Additional file 2: Figure S1.

Establishing the transmission dynamics of the vaccination intervention model of hepatitis E. Figure S2. Effects of cutting the transmission routes on the incidence of hepatitis E of Zhenjiang City. Figure S3. Effects of cutting the transmission routes on the incidence of hepatitis E of Yancheng City. Figure S4. Effects of cutting the transmission routes on the incidence of hepatitis E of Wuxi City. Figure S5. Simulation of shortening the incidence of hepatitis E infection of Zhenjiang City. Figure S6. Simulation of shortening the incidence of hepatitis E infection of Yancheng City. Figure S7. Simulation of shortening the incidence of hepatitis E infection of Wuxi City. Figure S8. Effects of different vaccination coefficients on the incidence of hepatitis E infection of Zhenjiang City. Figure S9. Effects of different vaccination coefficients on the incidence of hepatitis E infection of Yancheng City. Figure S10. Effects of different vaccination coefficients on the incidence of hepatitis E infection of Wuxi City. Figure S11. The sensitivity analysis of parameter.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Yang, M., Cheng, XQ., Zhao, ZY. et al. Feasibility of controlling hepatitis E in Jiangsu Province, China: a modelling study. Infect Dis Poverty 10, 91 (2021).

Download citation


  • Hepatitis E
  • Transmission dynamic model
  • Meteorological factor
  • Intervention
  • Transmissibility