This is an open-access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work, first published in JMIR Medical Informatics, is properly cited. The complete bibliographic information, a link to the original publication on http://medinform.jmir.org/, as well as this copyright and license information must be included.

Rabies is an acute infectious disease of the central nervous system caused by the rabies virus. The mortality rate of rabies is almost 100%. For some countries with poor sanitation, the spread of rabies among dogs is very serious.

The objective of this paper was to study the ecological transmission mode of rabies to make theoretical contributions to the suppression of rabies in China.

A mathematical model of the transmission mode of rabies was constructed using relevant data from the literature and officially published figures in China. Using this model, we fitted the data of the number of patients with rabies and predicted the future number of patients with rabies. In addition, we studied the effectiveness of different rabies suppression measures.

The results of the study indicated that the number of people infected with rabies will rise in the first stage, and then decrease. The model forecasted that in about 10 years, the number of rabies cases will be controlled within a relatively stable range. According to the prediction results of the model reported in this paper, the number of rabies cases will eventually plateau at approximately 500 people every year. Relatively effective rabies suppression measures include controlling the birth rate of domestic and wild dogs as well as increasing the level of rabies immunity in domestic dogs.

The basic reproductive number of rabies in China is still greater than 1. That is, China currently has insufficient measures to control rabies. The research on the transmission mode of rabies and control measures in this paper can provide theoretical support for rabies control in China.

Since ancient times, infectious diseases have brought much distress to human beings, and recorded deaths from infectious diseases abound. For example, the Black Death that killed a quarter of the European population during the medieval period was transmitted by fleas and black rats. In the mid-14th century, a plague killed another 25 million people [

Rabies is a zoonotic acute infectious disease of the central nervous system caused by the rabies virus. As patients with rabies have the clinical manifestation of being afraid of drinking water, this disease was previously called hydrophobia. However, animals with rabies do not have this characteristic [

Humans have been aware of the disease rabies since the time of Ancient Babylonians. At that time, there were no effective rabies control methods. Due to the pain and high fatality rate of rabies, many people committed suicide after being bitten by dogs [

In China, the vaccination rate of dogs is relatively low, especially in rural areas, where free-range dogs are rarely vaccinated. As such, the death toll from rabies is significant. Over the past two decades, the number of people who died of rabies in China first increased and then decreased. The number of deaths from rabies rose rapidly in the 1990s, reaching a peak in 2007. After the unremitting efforts of relevant personnel in China, the number of deaths from rabies began to decrease; deaths had decreased by about 76% as of 2013 [

Due to the highly infectious nature of the rabies virus, it is being studied by researchers worldwide. Simulation research is a research method that can transform complex biological phenomena into mathematical problems; such mathematical methods may enable better logical reasoning and disease transmission analysis. Some researchers have established mathematical models of rabies transmission from dogs to other animals and studied the transmission mode of rabies in animals to ascertain the optimal time for vaccine control [

To study the transmission mode of rabies and determine effective measures for suppressing rabies, we first constructed a mathematical model for the transmission mode of rabies. To construct the mathematical model, data from the literature and officially published data in China was used. The model fits the number of patients with rabies and predicts the number of patients with rabies in the future. In addition, we studied the effectiveness of different rabies suppression measures. This paper will provide theoretical support for rabies suppression in China.

Rabies transmission models of humans, wild dogs, and domestic dogs were established. Dogs and humans were divided into 4 categories: susceptibility, latency, infection, and recovery. S_{a}, E_{a}, I_{a}, and R_{a} represent the 4 categories of susceptibility, latency, infection, and recovery in wild dogs. S_{b}, E_{b}, I_{b}, and R_{b} indicate the 4 categories of susceptibility, latency, infection, and recovery in domestic dogs. S_{c}, E_{c}, I_{c}, and R_{c} represent the 4 categories of susceptibility, latency, infection, and recovery in humans. Before the simulation study, it was hypothesized that only dogs that have been in contact with an infected person would become infected and dogs in the incubation period have both vaccination and autoimmunity. The transmission model of rabies is shown in

According to this, the corresponding mathematical model was established. The mathematical models for transmission among dogs are shown in Equations 1-8.

The mathematical models for human transmission are shown in Equations 9-12.

Corresponding to the actual situation, the parameters in Equations 1-12 are nonnegative.

For both dog groups, A indicates the birth rate and p indicates the ratio of the two dog groups that are classified as susceptible in the incubation period of rabies but that have not developed disease.

_{a}S_{a}I_{a}_{a}^{’}S_{a}I_{b}_{b}S_{b}I_{b}_{b}^{’}S_{b}I_{a}_{a}_{a}’

For the human population, _{c} indicates the failure rate of vaccination in the human population. _{c}S_{c}I_{a}_{c}^{’}S_{c}I_{b}

Rabies transmission mode.

Model dynamics is used to study the relationship between rabies suppression measures and the number of patients with rabies. For dog groups, it is obtained as follows.

According to the actual situation, the number of groups in this mathematical model is always greater than 0. Then, the equation is obtained as shown below.

Thus, α can obtain the positive invariant set X of the model as follows.

Similarly, the positive invariant set X of the human situation can be obtained as follows:

The rabies-free equilibrium point

In dog groups, the BRN is

There is currently no fixed database of rabies in China. Therefore, the data in this paper were obtained from previously published literature and reports. The number of patients with rabies used for this paper corresponds to the data reported by the Public Health Scientific Data Center and the National Health and Family Planning Commission. Taking the year as the unit, the data shown in

Data list of parameters used for Equations 1-12.

Parameter | Symbol | Value |

Average number of wild dogs born each year | A_{a} |
2.47×10^{5} |

Infection rate of susceptible dogs by wild dogs in unit time | β_{a} |
2.91×10^{-6} |

Infection rate of susceptible dogs by domestic dogs per unit time | β_{b} |
2.20×10^{-7} |

Infection rate of susceptible human by wild dogs in unit time | β_{c} |
3.39×10^{-12} |

Killing rate of wild dogs | c_{a} |
0.06 |

Migration rate from latent wild dogs to infected dogs | σ_{a} |
0.35 |

Failure rate of vaccination in dogs | λ | 0.5 |

Average number of domestic dog births per year | A_{b} |
1.80×10^{6} |

Infection rate of susceptible wild dogs by domestic dogs in unit time | β_{a’} |
6.01×10^{-7} |

Infection rate of susceptible domestic dogs by wild dogs in unit time | β_{b’} |
5.01×10^{-7} |

Infection rate of susceptible humans by domestic dogs in unit time | β_{c’} |
6.79×10^{-12} |

Killing rate of domestic dogs | c_{a’} |
0.4 |

Immunity rate of wild dogs | γ_{a} |
0.05 |

Migration rate of humans from incubation period to infection period | σ_{c} |
0.33 |

Birth rate in the human population | H | 1.49×10^{7} |

The ratio of wild dogs and domestic dogs during the incubation period that did not have rabies outbreak and recovered to susceptible persons | p_{a} |
0.35 |

The ratio of domestic dogs in incubation period who did not have rabies outbreak and recovered to susceptible population | p_{b} |
0.37 |

The ratio of people in incubation period who did not have rabies outbreak and recovered to susceptible population | p_{c} |
0.33 |

Natural mortality of wild dogs | μ_{a} |
0.24 |

Disease-caused mortality of wild dogs | α_{a} |
1 |

Failure rate of vaccination in people | λ_{c} |
1 |

According to the number of patients with rabies in China in the past 10 years, the rabies disease data in the model were fitted. The initial values of each parameter are detailed below.

The initial values of susceptible individuals in the wild dog group, the domestic dog group, and the human population were 2×10^{6}, 3×10^{7}, and 1.29×10^{9}, respectively. The initial values of the latent individuals in the wild dog group, the domestic dog group, and the human population were 7×10^{4}, 2×10^{5}, and 400, respectively. The initial values of infected individuals in the wild dog group, the domestic dog group, and the human population were 2×10^{4}, 5×10^{4}, and 158, respectively. The initial values of recovered individuals in the wild dog group, the domestic dog group, and the human population were 1×10^{5}, 5×10^{6}, and 2×10^{5}, respectively.

Fitting of the model image was performed. During the image fitting process, other unknown parameters were obtained.

It can be seen from the above figure that the number of people with rabies rises at first and then decreases. From the fitted data, the number of patients with rabies peaked around 2005. As people's awareness of rabies increases, more people will be vaccinated after being injured by a dog, which will reduce the incidence of rabies and the number of cases. Additionally, due to the development of modern medicine, the effectiveness of the rabies vaccine is gradually increasing, which greatly reduces the number of patients with rabies. Compared to the number of rabies cases in 2016, the number of infected patients is undergoing a continuous decline. In about 10 years, the number of rabies cases is forecasted to be controlled within a relatively stable range. According to the prediction results of this model, the number of rabies cases will eventually be controlled to about 500 people every year. According to the values of the parameters in the model, the BRN is estimated to be 1.069.

Fitting and prediction data of the number of infected patients.

The effect of the birth number of wild (A_{a}) and domestic dogs (A_{b}) on the BRN is shown in

It can be seen from the figure above that when more wild dogs are killed, fewer wild dogs are born. Therefore, the value of the BRN can be reduced to less than 1 by killing wild dogs. A reduction in the BRN means that the number of people with rabies will also decrease. Therefore, by killing wild dogs, the incidence of rabies can be effectively reduced. Due to the Chinese people's love of dogs and the tradition of having dogs as pets, the number of domestic dogs is far higher than the number of wild dogs. Therefore, the management of domestic dogs in China is also challenging. As can be seen from the figure above, by controlling the birth number of domestic dogs, the value of the BRN can also be effectively controlled. Thus, the management of domestic dogs in China needs to be strengthened to a certain extent, especially the birth number of domestic dogs, which requires strict control.

The effect of the degree of immunity of wild (γ_{a}) and domestic (γ_{b}) dogs on the value of the BRN is shown in

It can be seen in

The values in

The effect of the birth number of wild (A_{a}) and domestic dogs (A_{b}) on the R0.

The effect of the degree of immunity of wild (γ_{a}) and domestic (γ_{b}) dogs on R0.

Rabies is caused by the rabies virus and can cause damage to the human central nervous system, with a fatality rate of 100%. Among infectious diseases, rabies has the highest fatality rate in the world [

In the process of data fitting, through the observation and calculation of each parameter, the BRN of rabies in China was determined to be about 1.069. This figure indicates that the current rabies response is not ideal and cannot effectively prevent the occurrence of rabies. In this paper, when analyzing the effects of each parameter on the BRN, it was found that the higher the number of wild dogs killed, the smaller the number of wild dogs born. Therefore, the value of the BRN can be reduced to less than 1 by killing wild dogs. Controlling the number of births of domestic dogs can also effectively control the value of the BRN. Therefore, the management of domestic dogs in China needs to be strengthened, especially the number of births of domestic dogs, which requires strict control. Although an increase in the degree of immunity of wild dogs can effectively reduce the BRN, the effect of this parameter on the BRN is small and cannot reduce the BRN to less than 1. Improving the degree of immunity of domestic dogs can reduce the value of the BRN to less than 1. The research results in this paper are similar to the results of research conducted by others on how rabies suppression measures affect the number of patients with rabies [

Due to the large number of patients with rabies in China and the severity of this condition, this paper investigated the transmission mode of rabies and rabies suppression measures. By counting the number of patients with rabies in China for many years and using software to fit the data, our model forecast an increase in the number of patients with rabies in China in the next few decades. In addition, this paper studied the effect of various measures on reducing the number of patients with rabies. Ultimately, controlling the birth rate of domestic dogs and wild dogs as well as increasing the degree of immunity in domestic dogs are all relatively effective rabies suppression measures. The BRN of rabies in China is still greater than 1. Therefore, in terms of rabies control, China still needs more research, policy formulation, and grassroots implementation.

Although the research in this paper obtained relatively considerable results, there are still some limitations: (1) The number of wild dogs in China is relatively high, and effective sterilization measures have not been implemented for wild dogs. Therefore, this problem was not considered during the construction of the model used in this paper. (2) The research in this paper is still in the theoretical stage, and further studies are required to determine whether it is effective in practice. Therefore, future research will be based on the results of this paper, and will observe whether the results are correct and feasible in practice.

basic reproductive number

None declared.