Next Article in Journal
Characterization of Regulatory T Cells in Patients Infected by Leishmania Infantum
Next Article in Special Issue
Retrospective Modeling of the Omicron Epidemic in Shanghai, China: Exploring the Timing and Performance of Control Measures
Previous Article in Journal
Vertical Transfer of Humoral Immunity against Nipah Virus: A Novel Evidence from Bangladesh
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Dynamic Compartmental Model to Explore the Optimal Strategy of Varicella Vaccination: An Epidemiological Study in Jiangsu Province, China

1
Department of Expanded Program on Immunization, Jiangsu Provincial Center for Disease Control and Prevention, Nanjing 210009, China
2
Department of Biomedical Engineering and Imaging Medicine, Army Medical University, Chongqing 400038, China
3
The First Clinical Medical College of Nanjing Medical University, Nanjing 210029, China
4
Department of Epidemiology and Biostatistics, School of Public Health, Nanjing Medical University, Nanjing 211166, China
5
Department of Mathematics and Statistics, Huaiyin Normal University, Huaian 223300, China
6
Department of Mathematics and Statistics, Southwest University, Chongqing 400715, China
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Trop. Med. Infect. Dis. 2023, 8(1), 17; https://doi.org/10.3390/tropicalmed8010017
Submission received: 12 November 2022 / Revised: 11 December 2022 / Accepted: 23 December 2022 / Published: 27 December 2022
(This article belongs to the Special Issue Response Strategies for Emerging Infectious Diseases)

Abstract

:
Varicella (chickenpox) is highly contagious among children and frequently breaks out in schools. In this study, we developed a dynamic compartment model to explore the optimal schedule for varicella vaccination in Jiangsu Province, China. A susceptible-infected-recovered (SIR) model was proposed to simulate the transmission of varicella in different age groups. The basic reproduction number was computed by the kinetic model, and the impact of three prevention factors was assessed through the global sensitivity analysis. Finally, the effect of various vaccination scenarios was qualitatively evaluated by numerical simulation. The estimated basic reproduction number was 1.831 ± 0.078, and the greatest contributor was the 5–10 year-old group (0.747 ± 0.042, 40.80%). Sensitivity analysis indicated that there was a strong negative correlation between the second dose vaccination coverage rate and basic reproduction number. In addition, we qualitatively found that the incidence would significantly decrease as the second dose vaccine coverage expands. The results suggest that two-dose varicella vaccination should be mandatory, and the optimal age of second dose vaccination is the 5–10 year-old group. Optimal vaccination time, wide vaccine coverage along with other measures, could enhance the effectiveness of prevention and control of varicella in China.

1. Background

Varicella is a highly contagious disease caused by varicella-zoster virus (VZV) [1,2]. VZV infection arises with the inhalation of aerosols containing pathogens or direct body contact. In the absence of varicella vaccination, varicella affects almost everyone by mid-adulthood [3]. The average annual incidence of varicella increased by 30% from 2005 to 2015 [4], and the number of people with varicella in China continues to rise. A seroepidemiological study in Suzhou, China, showed that about half of Chinese children aged 4–6 years were at high risk of contracting chickenpox virus [5]. In China, the direct and indirect damage, burden of disease and cost of medical services due to varicella is becoming increasingly severe [6,7,8]. All these pose a high economic burden on the public health system [9,10].
Universal varicella vaccination (UVV) has proven to be the most effective way to control the spread of VZV [1,2]. In the United States, every child is required to receive a one-dose varicella vaccination (VarV) after 1995 (the dose doubled in 2006), causing a dramatic decrease of the morbidity, hospitalization rate and mortality [2,11]. This immunization strategy has also been implemented in Germany, Australia, and Canada [12]. However, despite that varicella vaccine is recommended by the World Health Organization (WHO) [13], some developing countries have not incorporated it in national programs [14,15]. At present, some provinces and cities in China have started to implement the two-dose vaccination strategy, such as Beijing, Tianjin, Shanghai and Guangdong, which also use the above vaccination strategy.
In Jiangsu Province, the number of varicella cases has increased every year, even exponentially after 2014. Varicella emergency frequently occurs, with 128 outbreaks in 2017 and 219 outbreaks in 2018. Currently, only Suzhou (a city in Jiangsu) has implemented two-dose varicella vaccination free of charge, and the others have only given one or no dose. Many studies have shown that one-dose of varicella vaccine was not enough to prevent school outbreaks [16,17]. A two-dose vaccination strategy should be mandatory in China. In addition, an optimal schedule of vaccination should be designed based on many other factors. A modelling study from Italy demonstrated that coverage, efficacy, number of doses or dosing interval are also critical success factors for vaccination [18].
It is important to establish efficient strategies to control the spread of varicella in a highly susceptible population. The major objective of this research is to estimate the optimal strategy for varicella vaccination. First, we discussed the optimal time of second dose inoculation. Second, we analyzed the impact of other parameters on the varicella epidemic, which are vaccination coverage rate, vaccination efficacy and contact rate. Furthermore, simulations were conducted to investigate the effect of the two-dose vaccination coverage rate and related strategies and measures on varicella transmission.

2. Methods

2.1. Data Source

Every month from 2005 to 2018, we collected the reported varicella cases, either clinically diagnosed or lab-confirmed from Jiangsu Provincial Center for Disease Control and Prevention. All cases were voluntarily reported by local doctors.

2.2. Population

According to demographic information from the Statistics Bureau of Jiangsu Province, the population was divided into seven age groups, and k = 1 , 2 , 3 , 4 , 5 , 6 , 7 was used to represent the >1 year (0–1), 1 ≤ age < 3 (1–3), 3 ≤ age < 5 (3–5), 5 ≤ age < 10 (5–10), 10 ≤ age < 15 (10–15), 15 ≤ age < 20 (15–20), >20, respectively. The number of new births per year was acquired from the Statistics Bureau [19]. For simplicity, the difference in disease-related mortality was ignored, which was taken as 7‰ from Statistics Bureau in all age groups [19].
According to the primary epidemiological properties of varicella, we divided the total population of each age group at time t into three compartments: susceptible ( S ( k , t ) ), infected ( I ( k , t ) ) and recovered or immune ( R ( k , t ) ). Considering that vaccination in the susceptible may fail, we further divided the susceptible of each age groups into susceptible without vaccination ( S 0 ( k , t ) ), the susceptible with first dose vaccination but failed ( S 1 ( k , t ) ) and the susceptible with second dose vaccinations but all failed ( S 2 ( k , t ) ).

2.3. Vaccination Strategy

Children below ten years old (0–1, 1–3, 3–5 and 5–10 age groups) in any condition can be vaccinated with the first or second dose. We assumed that the first dose vaccination can only be given to individuals susceptible but not yet vaccinated ( S 0 ( k , t ) ), and the second dose vaccination to those susceptible with failed primary vaccination ( S 1 ( k , t ) ). Note that for individuals who were susceptible (due to primary failure) but whose age was less than one year old ( S 1 ( 1 , t ) ) they were not given the second dose, because the guidelines recommend that the minimum interval between the two doses must be greater than 12 months [20].
The effective rates of the first and second doses were presented as constants ε 1 and ε 2 respectively, and the corresponding vaccination coverage rates in the susceptible people at kth age group as η 1 ( k ) and η 2 ( k ) respectively. We obtained the vaccination outcomes in each age group (Table 1). Here, V 0 f ( k , t ) and V 0 s ( k , t ) were defined as individuals with failed and successful first dose vaccination in the kth age group, respectively. Similarly, V 1 f ( k , t ) and V 1 s ( k , t ) were the individuals with failed and successful second dose vaccination in the kth age group, respectively.

2.4. Model Structure

According to the above mentioned population division and vaccination strategy, as well as the epidemic characteristics of varicella, we set up the compartments in each age group (Figure 1). For the kth age group, the flowchart describes how individuals can move among the five mutually exclusive states. Considering the interpersonal behaviors, we roughly assumed that every susceptible person may contact all infected persons and become infected.
Based on the flowchart in the graphical abstract image, the following deterministic SIR dynamic model with age structure can be proposed to mimic the varicella epidemic in Jiangsu province, China.
d S 0 ( k , t ) d t = F 0 ( k 1 , t ) V 0 f ( k , t ) V 0 s ( k , t ) β ( k ) S 0 ( k , t ) k = 1 7 I ( k , t ) F 0 ( k , t ) d S 0 ( k , t ) , d S 1 ( k , t ) d t = F 1 ( k 1 , t ) + V 0 f ( k , t ) V 1 f ( k , t ) V 1 s ( k , t ) β ( k ) S 1 ( k , t ) k = 1 7 I ( k , t ) F 1 ( k , t ) d S 1 ( k , t ) , d S 2 ( k , t ) d t = F 2 ( k 1 , t ) + V 1 f ( k , t ) β ( k ) S 2 ( k , t ) k = 1 7 I ( k , t ) F 2 ( k , t ) d S 2 ( k , t ) , d I ( k , t ) d t = β ( k ) S 0 ( k , t ) + S 1 ( k , t ) + S 2 ( k , t ) k = 1 7 I ( k , t ) c I ( k , t ) d I ( k , t ) , d R ( k , t ) d t = F R ( k 1 , t ) + V 0 s ( k , t ) + V 1 s ( k , t ) + c I ( k , t ) F R ( k , t ) d R ( k , t ) ,
For simplicity, (1) is only a sub-model of the kth age group, and the full model of varicella transmission is given in Supplementary File S1. Here F i ( k 1 , t ) represents the transfer rate from (k − 1)th age group to the kth age group of the susceptible population who have received ith vaccination at time t , i = 0 , 1 , 2 and k = 2 , 3 , , 7 . F R ( k 1 , t ) represents the transfer rate from (k − 1)th age group to the kth age group of the recovered population.
Since the infectious period of varicella is only one month, we assumed that all the infected people recovered at the same age of infection, so they do not flow between different age groups, and the transfer due to age growth only appeared in the first, second, third and last equations in (1). The value of parameters and initial conditions were detailed in Supplementary File S2.

2.5. Simulations

Based on the reported varicella cases in all age groups and the proposed dynamic model, an adaptive Metropolis-Hastings algorithm was used to carry out the Markov Chain Monte Carlo (MCMC) procedure to estimate the undetermined parameters in model (Supplementary File S3 for details). The algorithm runs for 2 × 106 iterations with a burn-in of 106 iterations. The Geweke convergence diagnostic method was employed to assess the convergence [21]. The median values of the estimated parameters and their corresponding first and third quartiles (Q1–Q3) were listed in Table S3, and their distribution were depicted in Figure S1.

2.6. Basic Reproduction Number

Basic reproduction number, R0, is a key parameter in infectious epidemiology, which indicates the expected number of secondary infections produced by an infected case. If R0 < 1, this disease cannot invade the population and will eventually disappear, otherwise, if R0 > 1, the disease will continue to spread and become endemic [22]. See Supplementary File S4 for the detailed calculation process.

2.7. Global Sensitivity Analysis

Because the basic reproduction number represents the epidemic intensity of the disease, we used it as a dependent variable in Latin Hypercube Sampling (LHS). Partial rank correlation coefficients (PRCCs) were used to examine the global sensitivity of the parameters in the model. Clearly, expanding vaccination coverage, reducing the contact frequency with infected population, and improving vaccination efficacy are all potentially feasible strategies for better preventing and controlling varicella epidemics. The effect of these three parameters on the basic regeneration number is discussed further in Supplementary File S5.

2.8. Vaccination Scenarios

In this study, in addition to the effective contact rate, it has been found that increasing the coverage rate of vaccination is more conducive to disease control than increasing the efficiency of vaccination. Thus, in order to simulate the effect of the second vaccination since 2018, we simulated the long-term trend of varicella epidemic in the next 50 years according to the following two scenarios: (1) increase the coverage of the second vaccination while keeping the contact rate unchanged in the model; (2) increase the coverage of the second vaccination when the contact rate decreases in the model, i.e., assuming that the contact rate is 80% of its original baseline value.

3. Results

3.1. Varicella Data and Underdetermined Parameters

The annual varicella cases in Jiangsu Province from 2005 to 2018 (red circles) and the best-fit simulation of full model with the fitted parameters (Figure 1). According to the goodness-of-fit test for nonlinear fitting [23], our model achieved 0.8211 in the overall goodness, which indicates that the model results were in good agreement with the reported data. We observed that the reported varicella cases in all age groups have been growing exponentially from 2005 to 2018 in Jiangsu Province.

3.2. Basic Reproduction Numbers

The overall basic reproduction number can be decomposed as R 0 = k = 1 7 R 0 ( k ) , where R 0 ( k ) was the contribution in each age group (See Supplementary File S4 for detailed procedures). From Figure 2A, we observed that the overall basic reproduction number was 1.831 ± 0.078, indicating that varicella has become an endemic disease in Jiangsu Province. Furthermore, from high to low, the contribution of each age group to the epidemic was 5–10-years-old (0.747 ± 0.042, 40.80%), 10–15-years-old (0.401 ± 0.035, 21.90%), 3–5-years-old (0.219 ± 0.012, 11.98%), more than 20-years-old (0.213 ± 0.027, 11.66%), 15–20-years-old (0.128 ± 0.020, 6.98%), 1–3-years-old (0.085 ± 0.005, 4.62%), 0–1-year-old (0.038 ± 0.003, 2.06%).

3.3. Sensitivity Analysis

Taking the basic reproduction number as an output variable, we calculated the PRCCs of the three factors, including vaccination coverage rate, vaccination efficacy and effective contact rate. From Figure 2B, we found that almost all parameters were sensitive to the basic reproduction number except η 1 ( 1 ) (first dose vaccination coverage rate in the first age group) and β ( 1 ) (the effective contact rate in the first age group), which can explain the low percentage of basic reproduction number in the total in the first age group (Figure 2A). For the vaccination coverage related parameters, the simulation results showed that the PRCC values of the second dose vaccination coverage were always greater than those of the first dose, which indicates that the second immunization is critical to control the varicella epidemic. The urgent task is to find ways to improve the second vaccination coverage rate. For effective contact rate parameters, we found that the order of PRCC values of each age group was almost as same as that of their corresponding basic reproduction number (Figure 2). The age groups with the largest PRCC values were 5–10-years-old, 10–15-years-old and 15–20-years-old, respectively, and the possible explanation is that these age groups have long-time social activities, such as attending school [24]. The parameters of both vaccinations have a very significant negative impact on the prevalence of varicella.

3.4. Association between Basic Reproduction Number and the Three Prevention Factors

According to the basic reproduction number in Supplementary File S4, we simulated the critical basic reproduction number surface (Figure 3A) and the corresponding contour maps (Figure 3B–D) with these parameters and their different combinations.

3.5. Impact of Second Dose Vaccination Coverage Rate on the Varicella Transmission

Based on the methods described earlier, Figure 4 shows the qualitative prediction of a varicella epidemic if we increase second dose vaccination coverage under two scenarios. The results indicated that increasing the coverage rate is always beneficial to control the varicella epidemic, which can not only reduce the basic reproduction number, but also reduce the peak of the epidemic or make the peak earlier. Based on Figure 4, it is not surprising that reducing the contact rate at the same time (Figure 4B) can better control the epidemic than keeping the contact rate unchanged (Figure 4A).

4. Discussion

In this study, we developed a dynamic compartmental model to explore the optimal age for second dose varicella vaccination and benefits of increasing the coverage rate. Our model fitted well with annually reported varicella cases from 2005 to 2018 in Jiangsu Province, China. Based on the SIR model, the estimated basic reproduction number was 1.831 ± 0.078, indicating that varicella is endemic in Jiangsu Province. The WHO recommends that in countries where varicella is a major public health issue, varicella vaccination should be introduced into their routine immunization programs [24]. To date, the free second dose of varicella vaccination is still not fully promoted in Jiangsu Province.
To our knowledge, this study is the first to use a mathematic model to evaluate the optimal age for the varicella vaccination strategy in Jiangsu Province. We concluded that the optimal age for second dose varicella vaccination is 5–10 years old, with more evident results illustrated in Figure 2. For the two-dose vaccination schedule there is continuous debate about whether a short-term (months) or longer term (years) interval should be kept. Some hold that a short interval between two doses would lead to a lower risk of varicella breakthrough infections, but with less durable protection over time, and a long interval would produce durable immunity, but children might face a higher risk of breakthrough infections in their early life [25]. A recent study [25] used an agent-based-model to compare the effectives of the two vaccination schedules: (1) first dose at 12 months, second dose at 18 months, and (2) first dose at 12 months, second dose at 4–6 years. In that study, the model predicted that the second dose at 18 months may be more effective than that at 4–6 years, but the difference between the two schedules was insignificant. Therefore, policymakers should consider other factors to determine the appropriate varicella vaccine schedule. Considering factors like security, simplicity and persistence of immunity, the second dose should be given to an older group. Since varicella frequently breaks out in the kindergarten and primary schools [26], it is reasonable to come to this conclusion.
Several other factors also determine the effectiveness of UVV. Sensitivity analysis confirmed a stronger negative correlation between the coverage rate and basic reproduction number, especially for the second dose (Figure 4B). However, model simulation revealed that the second dose coverage rate was far lower than the first dose (about 5% vs. 65%). The WHO has already recommended that vaccines should cover above 80% of the population [27]. High vaccination coverage not only effectively reduces the incidence of breakthrough varicella, but also generates a herd protection effect [28]. We found that when the second dose coverage rate increased by 7.347 times, the varicella transmission could be terminated, at the same time, with contract rate reduced by 80%, and varicella cases could be reduced 3.751 times (Figure 4. Holl et al. conducted a model-based study, suggesting that increasing the coverage rate can compensate for less effective vaccination strategies, such as a longer interval between two doses and lower second dose efficacy [18]. Similarly, the combination of these measures can control varicella more effectively.
Our kinetic model provides an effective method to study the effect of the second dose of vaccine on the spread of the disease. Dynamic models are widely used to generate optimal policies for controlling the spread of infectious diseases [29]. Using mathematical models, we can analyze big data or make predictions based on a limited data in epidemiology.
There are several limitations in our study. First, we only considered the varicella cases from Jiangsu Province, thus the generalization of the conclusion was limited to some extent. Secondly, this paper only studied the effect of vaccination strategies on varicella, ignoring the effect on herpes zoster. This is because the study was not designed to prove the causality between herpes zoster and varicella vaccination. Third, many factors affecting the implementation of vaccination strategy should be taken into consideration, including duration of vaccination effectiveness, catch-up program, breakthrough varicella, and financial cost.

5. Conclusions

Increasing vaccine coverage can improve the prevention and control of varicella. The two-dose varicella vaccination should be compulsory, and the optimal age of two-dose vaccination is 5–10 years in Jiangsu Province, China.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/tropicalmed8010017/s1, Supplementary File S1: Full model description; Supplementary File S2: Parameters and initial conditions; Supplementary File S3: Parameter estimation; Supplementary File S4: Basic reproduction number; Supplementary File S5: Global sensitivity analysis.

Author Contributions

X.S., C.D. and K.W. (Kai Wang) curated and provided data and drafted the manuscript. Y.L., C.W. and X.J. screened the data. Y.Y. and Z.D. performed laboratory tests. W.W., Z.L. and Z.W. provided guidance for epidemiological analysis throughout the study. F.T., K.W. (Kaifa Wang) and Z.P. adjusted the research framework, provided financial support, and revised the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (Grant number 82073673, 12171396), the National S&T Major Project Foundation of China (2018ZX10715002) and the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD).

Institutional Review Board Statement

Our research have been performed in accordance with the Declaration of Helsinki and have been approved by the ethics committee of Nanjing Medical University (“F”, “CH”, “Nanjing Med U”, “FWA00001501”, “NANJING”, 11/21/2004). In this paper an age-structured susceptible-infected-recovered (SIR) model was proposed to simulate the transmission of varicella based on reported varicella cases of Jiangsu Province from 2005–2018. Ethical review and approval were waived for this study due to all data are from public data and do not involve personal privacy data, and there was no intervention treatment. The data used in this study was anonymised before the analyses. Further ethics statement was not needed based on the internal regulations.

Informed Consent Statement

Patient consent was waived due to all data are from public data and do not involve personal privacy data, and there was no intervention treatment. The data used in this study was anonymised before the analyses. Further ethics statement was not needed based on the internal regulations.

Data Availability Statement

All data generated or analysed during this study are included in this published article. The data that support the findings of this study are available from the corresponding author upon reasonable request.

Acknowledgments

The authors thank the study stuff in Jiangsu Provincial Center for Disease Control and Prevention for data management within the Infectious disease surveillance system and thank all participants enrolled in this study.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

Susceptible-infected-recovered (SIR); Varicella-zoster virus (VZV); herpes zoster (HZ); universal varicella vaccination (UVV); varicella vaccination (VarV); the World Health Organization (WHO); Latin Hypercube Sampling (LHS). Partial rank correlation coefficients (PRCCs).

References

  1. Gershon, A.A. Is chickenpox so bad, what do we know about immunity to varicella zoster virus, and what does it tell us about the future? J. Infect. 2017, 74, S27–S33. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Wolfson, L.J.; Daniels, V.J.; Altland, A.; Black, W.; Huang, W.; Ou, W. The Impact of Varicella Vaccination on the Incidence of Varicella and Herpes Zoster in the United States: Updated Evidence from Observational Databases, 1991–2016. Clin. Infect. Dis. 2020, 70, 955–1002. [Google Scholar] [CrossRef] [PubMed]
  3. Xu, Y.; Liu, Y.; Zhang, X.; Zhang, X.; Du, J.; Cai, Y.; Wang, J.; Che, X.; Gu, W.; Jiang, W.; et al. Epidemiology of varicella and effectiveness of varicella vaccine in Hangzhou, China, 2019. Hum. Vaccine Immunother. 2021, 17, 211–216. [Google Scholar] [CrossRef] [PubMed]
  4. Siu, H.T.; Li, J.C.; Wang, M.; Liu, Y.M.; Yin, D.P. Varicella epidemiology in China, 2005–2015. Chin. J. Vaccines Immun. 2019, 2, 155–159. [Google Scholar]
  5. Luan, L.; Shen, X.; Qiu, J.; Jing, Y.; Zhang, J.; Wang, J.; Zhang, J.; Dong, C. Seroprevalence and molecular characteristics of varicella-zoster virus infection in Chinese children. BMC Infect. Dis. 2019, 19, 643. [Google Scholar] [CrossRef] [Green Version]
  6. Pan, X.; Shu, M.; Ma, R.; Fang, T.; Dong, H.; Sun, Y.; Xu, G. Varicella breakthrough infection and effectiveness of 2-dose varicella vaccine in China. Vaccine 2018, 36, 5665–5670. [Google Scholar] [CrossRef]
  7. Britton, P.N.; Dale, R.C.; Blyth, C.C.; Clark, J.E.; Crawford, N.; Marshall, H.; Elliott, E.J.; Macartney, K.; Booy, R.; Jones, C.A. Causes and Clinical Features of Childhood Encephalitis: A Multicenter, Prospective Cohort Study. Clin. Infect. Dis. 2019, 70, 2517–2526. [Google Scholar] [CrossRef]
  8. Gowin, E.; Wysocki, J.; Michalak, M. Don’t forget how severe varicella can be—Complications of varicella in children in a defined Polish population. Int. J. Infect. Dis. 2013, 17, e485–e489. [Google Scholar] [CrossRef] [Green Version]
  9. Riera-Montes, M.; Bollaerts, K.; Heininger, U.; Hens, N.; Gabutti, G.; Gil, A.; Nozad, B.; Mirinaviciute, G.; Flem, E.; Souverain, A.; et al. Estimation of the burden of varicella in Europe before the introduction of universal childhood immunization. BMC Infect. Dis. 2017, 17, 353. [Google Scholar] [CrossRef] [Green Version]
  10. Gater, A.; Uhart, M.; McCool, R.; Préaud, E. The humanistic, economic and societal burden of herpes zoster in Europe: A critical review. BMC Public Health 2015, 15, 193. [Google Scholar] [CrossRef] [Green Version]
  11. Baxter, R.; Ray, P.; Tran, T.N.; Black, S.; Shinefield, H.R.; Coplan, P.M.; Lewis, E.; Fireman, B.; Saddier, P. Long-term effectiveness of varicella vaccine: A 14-Year, prospective cohort study. Pediatrics 2013, 131, e1389–e1396. [Google Scholar] [CrossRef] [PubMed]
  12. Yin, M.; Xu, X.; Liang, Y.; Ni, J. Effectiveness, immunogenicity and safety of one vs. two-dose varicella vaccinationa meta-analysis. Expert Rev. Vaccines 2018, 17, 351–362. [Google Scholar] [CrossRef] [PubMed]
  13. Barrenechea, G.G.; Bastos, L.S. Evaluation of impact of one dose varicella vaccine on the incidence of chickenpox in Argentina. Vaccine 2020, 38, 330–335. [Google Scholar] [CrossRef] [PubMed]
  14. Hoshi, S.L.; Kondo, M.; Okubo, I. Cost-effectiveness of varicella vaccine against herpes zoster and post-herpetic neuralgia for elderly in Japan. Vaccine 2017, 35, 3264–3271. [Google Scholar] [CrossRef] [PubMed]
  15. Suo, L.; Lu, L.; Wang, Q.; Yang, F.; Wang, X.; Pang, X.; Marin, M.; Wang, C. Varicella outbreak in a highly-vaccinated school population in Beijing, China during the voluntary two-dose era. Vaccine 2017, 35, 4368–4373. [Google Scholar] [CrossRef] [Green Version]
  16. Lopez, A.S.; Guris, D.; Zimmerman, L.; Gladden, L.; Moore, T.; Haselow, D.T.; Loparev, V.N.; Schmid, D.S.; Jumaan, A.O.; Snow, S.L. One dose of varicella vaccine does not prevent school outbreaks: Is it time for a second dose? Pediatrics 2006, 117, e1070–e1077. [Google Scholar] [CrossRef] [Green Version]
  17. Perella, D.; Wang, C.; Civen, R.; Viner, K.; Kuguru, K.; Daskalaki, I.; Schmid, D.S.; Lopez, A.S.; Tseng, H.F.; Newbern, E.C.; et al. Varicella Vaccine Effectiveness in Preventing Community Transmission in the 2-Dose Era. Pediatrics 2016, 137, e20152802. [Google Scholar] [CrossRef] [Green Version]
  18. Holl, K.; Sauboin, C.; Amodio, E.; Bonanni, P.; Gabutti, G. Coverage, efficacy or dosing interval: Which factor predominantly influences the impact of routine childhood vaccination for the prevention of varicella? A model-based study for Italy. BMC Public Health 2016, 16, 1103. [Google Scholar] [CrossRef] [Green Version]
  19. Statistics Bureau of Jiangsu Province. Jiangsu Statistical Yearbook 2019. Available online: http://tj.jiangsu.gov.cn/2019/nj03/nj0303.htm (accessed on 6 September 2019).
  20. Gabutti, G.; Bolognesi, N.; Sandri, F.; Florescu, C.; Stefanati, A. Varicella zoster virus vaccines: An update. ImmunoTargets Ther. 2019, 8, 15–28. [Google Scholar] [CrossRef] [Green Version]
  21. Haario, H.; Laine, M.; Mira, A.; Saksman, E. DRAM: Efficient adaptive MCMC. Stat. Comput. 2006, 16, 339–354. [Google Scholar] [CrossRef]
  22. Van den Driessche, P.; Watmough, J. Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. Math. Biosci. 2002, 180, 29–48. [Google Scholar] [CrossRef] [PubMed]
  23. Zhang, S. Approach on the fitting optimization index of curve regression. Chin. J. Health Stat. 2002, 19, 9–11. [Google Scholar]
  24. Lu, L.; Suo, L.; Li, J.; Zhai, L.; Zheng, Q.; Pang, X.; Bialek, S.R.; Wang, C. A varicella outbreak in a school with high one-dose vaccination coverage, Beijing, China. Vaccine 2012, 30, 5094–5098. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Rafferty, E.R.S.; McDonald, W.; Osgood, N.D.; Qian, W.; Doroshenko, A. Seeking the optimal schedule for chickenpox vaccination in Canada: Using an agent-based model to explore the impact of dose timing, coverage and waning of immunity on disease outcomes. Vaccine 2020, 38, 521–529. [Google Scholar] [CrossRef] [PubMed]
  26. Wu, Q.; Zhang, L.; Liu, J.; Chen, Y.; Zhou, Q.; Xuan, Z.; Hu, Y. Factors Influencing the Size and Duration of School Varicella Outbreaks in Shanghai, China, in the Era of Voluntary Single-Dose Vaccination. Open Forum Infect. Dis. 2019, 6, ofz410. [Google Scholar] [CrossRef] [Green Version]
  27. Wutzler, P.; Bonanni, P.; Burgess, M.; Gershon, A.; Sáfadi, M.A.; Casabona, G. Varicella vaccination—The global experience. Expert Rev. Vaccines 2017, 16, 833–843. [Google Scholar] [CrossRef] [Green Version]
  28. Drolet, M.; Bénard, É.; Boily, M.-C.; Ali, H.; Baandrup, L.; Bauer, H.; Beddows, S.; Brisson, J.; Brotherton, J.M.L.; Cummings, T.; et al. Population-level impact and herd effects following human papillomavirus vaccination programmes: A systematic review and meta-analysis. Lancet Infect. Dis. 2015, 15, 565–580. [Google Scholar] [CrossRef] [Green Version]
  29. Cardim, L.L.; Pinho, S.T.R.; Teixeira, M.G.; Costa, M.C.N.; Esteva, M.L.; Ferreira, C.P. Heterogeneities in dengue spatial-temporal transmission in Brazilian cities and its influence on the optimal age of vaccination. BMC Public Health 2019, 19, 155. [Google Scholar] [CrossRef]
Figure 1. Illustration of the fitting result of the full model from 2005 to 2018 underestimated parameters. The blue solid line and azure shade represent the model-predicted cases and 95% confidence interval (CI). Red circles denote survey-reported cases. (A) is the fitting result of total infected population, and (BH) is the fitting result of each age group.
Figure 1. Illustration of the fitting result of the full model from 2005 to 2018 underestimated parameters. The blue solid line and azure shade represent the model-predicted cases and 95% confidence interval (CI). Red circles denote survey-reported cases. (A) is the fitting result of total infected population, and (BH) is the fitting result of each age group.
Tropicalmed 08 00017 g001
Figure 2. (A) Estimated basic reproduction number based on the full model. (B) Partial rank correlation coefficients illustrating the sensitivity of basic reproduction number on three types of parameters, named as vaccination coverage rate, vaccination efficacy and effective contact rate, which are shown by blue, green and orange bars, respectively. The red dashed line in (A) denotes the critical threshold. The blue and red lines in (B) represent moderate and strong correlation thresholds of PRCC values, respectively.
Figure 2. (A) Estimated basic reproduction number based on the full model. (B) Partial rank correlation coefficients illustrating the sensitivity of basic reproduction number on three types of parameters, named as vaccination coverage rate, vaccination efficacy and effective contact rate, which are shown by blue, green and orange bars, respectively. The red dashed line in (A) denotes the critical threshold. The blue and red lines in (B) represent moderate and strong correlation thresholds of PRCC values, respectively.
Tropicalmed 08 00017 g002
Figure 3. (A) The surface of R0 = 1 in terms of parameters related to vaccination coverage rate, vaccination efficacy and contact rate. The surface above the plane represents R0 < 1 and otherwise R0 > 1. The contours of R0 versus vaccination coverage rate and vaccination efficacy (B), versus vaccination coverage rate and contact rate (C), versus the contact rate and vaccination efficacy (D).
Figure 3. (A) The surface of R0 = 1 in terms of parameters related to vaccination coverage rate, vaccination efficacy and contact rate. The surface above the plane represents R0 < 1 and otherwise R0 > 1. The contours of R0 versus vaccination coverage rate and vaccination efficacy (B), versus vaccination coverage rate and contact rate (C), versus the contact rate and vaccination efficacy (D).
Tropicalmed 08 00017 g003
Figure 4. The impact of increasing second dose vaccination coverage rate staring in 2018 on the varicella transmission with (A) unchanged contact rate and (B) reduced contact rate (80% baseline). η2 represents all parameters related to second dose vaccination. β represents all parameters of contact rate.
Figure 4. The impact of increasing second dose vaccination coverage rate staring in 2018 on the varicella transmission with (A) unchanged contact rate and (B) reduced contact rate (80% baseline). η2 represents all parameters related to second dose vaccination. β represents all parameters of contact rate.
Tropicalmed 08 00017 g004
Table 1. Vaccination outcomes in each age group.
Table 1. Vaccination outcomes in each age group.
Age Group V 0 f ( k , t ) V 0 s ( k , t ) V 1 f ( k , t ) V 1 s ( k , t )
0–1 year η 1 ( 1 ) ( 1 ε 1 ) S 0 ( 1 , t ) η 1 ( 1 ) ε 1 S 0 ( 1 , t ) NANA
1–3 years η 1 ( 2 ) ( 1 ε 1 ) S 0 ( 2 , t ) η 1 ( 2 ) ε 1 S 0 ( 2 , t ) η 2 ( 2 ) ( 1 ε 2 ) S 1 ( 2 , t ) η 2 ( 2 ) ε 2 S 1 ( 2 , t )
3–5 years η 1 ( 3 ) ( 1 ε 1 ) S 0 ( 3 , t ) η 1 ( 3 ) ε 1 S 0 ( 3 , t ) η 2 ( 3 ) ( 1 ε 2 ) S 1 ( 3 , t ) η 2 ( 3 ) ε 2 S 1 ( 3 , t )
5–10 years η 1 ( 4 ) ( 1 ε 1 ) S 0 ( 4 , t ) η 1 ( 4 ) ε 1 S 0 ( 4 , t ) η 2 ( 4 ) ( 1 ε 2 ) S 1 ( 4 , t ) η 2 ( 4 ) ε 2 S 1 ( 4 , t )
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Sun, X.; Dai, C.; Wang, K.; Liu, Y.; Jin, X.; Wang, C.; Yin, Y.; Ding, Z.; Lu, Z.; Wang, W.; et al. A Dynamic Compartmental Model to Explore the Optimal Strategy of Varicella Vaccination: An Epidemiological Study in Jiangsu Province, China. Trop. Med. Infect. Dis. 2023, 8, 17. https://doi.org/10.3390/tropicalmed8010017

AMA Style

Sun X, Dai C, Wang K, Liu Y, Jin X, Wang C, Yin Y, Ding Z, Lu Z, Wang W, et al. A Dynamic Compartmental Model to Explore the Optimal Strategy of Varicella Vaccination: An Epidemiological Study in Jiangsu Province, China. Tropical Medicine and Infectious Disease. 2023; 8(1):17. https://doi.org/10.3390/tropicalmed8010017

Chicago/Turabian Style

Sun, Xiang, Chenxi Dai, Kai Wang, Yuanbao Liu, Xinye Jin, Congyue Wang, Yi Yin, Zhongxing Ding, Zhenzhen Lu, Weiming Wang, and et al. 2023. "A Dynamic Compartmental Model to Explore the Optimal Strategy of Varicella Vaccination: An Epidemiological Study in Jiangsu Province, China" Tropical Medicine and Infectious Disease 8, no. 1: 17. https://doi.org/10.3390/tropicalmed8010017

APA Style

Sun, X., Dai, C., Wang, K., Liu, Y., Jin, X., Wang, C., Yin, Y., Ding, Z., Lu, Z., Wang, W., Wang, Z., Tang, F., Wang, K., & Peng, Z. (2023). A Dynamic Compartmental Model to Explore the Optimal Strategy of Varicella Vaccination: An Epidemiological Study in Jiangsu Province, China. Tropical Medicine and Infectious Disease, 8(1), 17. https://doi.org/10.3390/tropicalmed8010017

Article Metrics

Back to TopTop