- 1Research Center for Veterinary Biosciences and Veterinary Public Health, Faculty of Veterinary Medicine, Chiang Mai University, Chiang Mai, Thailand
- 2Veterinary Public Health and Food Safety Centre for Asia Pacific (VPHCAP), Faculty of Veterinary Medicine, Chiang Mai University, Chiang Mai, Thailand
- 3Department of Veterinary Biosciences and Veterinary Public Health, Faculty of Veterinary Medicine, Chiang Mai University, Chiang Mai, Thailand
- 4The 4th Regional Livestock Office, Department of Livestock Development, Khon Kaen, Thailand
- 5Department of Livestock Development, Bangkok, Thailand
- 6College of Veterinary Science and Medicine, Central Luzon State University, Science City of Muñoz, Nueva Ecija, Philippines
- 7Department of Mathematics, Faculty of Science, Naresuan University, Phitsanulok, Thailand
Introduction: In 2021, Thailand reported the highest incidence of lumpy skin disease (LSD) outbreaks in Asia. In response to the widespread outbreaks in cattle herds, the government's livestock authorities initiated comprehensive intervention measures, encompassing control strategies and a national vaccination program. Yet, the efficacy of these interventions remained unevaluated. This research sought to assess the nationwide intervention's impact on the incidence of new LSD cases through causal impact analysis.
Methods: Data on weekly new LSD cases in Thailand from March to September 2021 was analyzed. The Bayesian structural time series (BSTS) analysis was employed to evaluate the causal relationship between new LSD cases in the pre-intervention phase (prior to the vaccination campaign) and the post-intervention phase (following the vaccination campaign). The assessment involved two distinct scenarios, each determined by the estimated effective intervention dates. In both scenarios, a consistent decline in new LSD cases was observed after the mass vaccination initiative, while other control measures such as the restriction of animal movement, insect control, and the enhancement of the active surveillance approach remained operational throughout the pre-intervention and the post-intervention phases.
Results and discussion: According to the relative effect results obtained from scenario A and B, it was observed that the incidence of LSD cases exhibited reductions of 119% (95% Credible interval [CrI]: −121%, −38%) and 78% (95% CrI: −126, −41%), respectively. The BSTS results underscored the significant influence of these interventions, with a Bayesian one-sided tail-area probability of p < 0.05. This model-based study provides insight into the application of BSTS in evaluating the impact of nationwide LSD vaccination based on the national-level data. The present study is groundbreaking in two respects: it is the first study to quantify the causal effects of a mass vaccination intervention on the LSD outbreak in Thailand, and it stands as the only endeavor of its kind in the Asian context. The insights collected from this study hold potential value for policymakers in Thailand and other countries at risk of LSD outbreaks.
1 Introduction
Lumpy skin disease (LSD) is a significant transboundary disease affecting the bovine population causing substantial economic losses in cattle industries across countries (1). The disease is caused by the lumpy skin disease virus (LSDV) belonging to the Capripoxvirus genus, of the Poxviridae family (2). Its primary symptoms include nodular lesions on the skin and mucous membranes, enlargement of lymph nodes, and elevated body temperature (3, 4). Affected cattle may experience weight loss, decreased milk production, and deteriorated skin quality, all contributing to economic losses. While the morbidity rate of LSD is moderate, its mortality rate remains low (5). The disease primarily spreads through insect vectors (1) but can be effectively controlled and prevented with vaccines (6–8). The primary strategies for managing and eradicating LSD are restricting animal movement and mass vaccination (5, 9). Given its potential for swift spread and significant economic impact, the World Organization for Animal Health (WAOH) classifies LSD as a notifiable disease, necessitating reporting of its occurrences (5). While LSD is endemic in Africa (2), it reached the Middle East in 2012, Europe in 2015, and had outbreaks in Russia between 2017 and 2019. The disease further spread to the South and East Asia between 2019–2020 (3, 10, 11). From 2020 onwards, several Southeast Asian countries reported multiple LSD outbreaks (12–14) with Thailand recording the highest incidence in 2021 (15, 16).
Thailand reported the first LSD outbreak in cattle farming areas in the northeastern region in April 2021 (17). Within a few months, the disease had spread nationwide (18–20). In response, starting from mid-April 2021, authorities in charge of livestock implemented various control measures including regulating cattle movement, disinfecting affected areas, managing insect vectors, treating affected animals, closing live cattle markets, and instituting active surveillance to identify infected premises. Local veterinarians and livestock officers across the country took on the responsibility to control LSD outbreaks in their respective areas to contain the swift dissemination of lumpy skin disease viruses from outbreak areas to other regions (18, 21, 22). Additionally, a public awareness campaign was initiated to spread information about preventive measures (18). This campaign involved in-person meetings between livestock authorities and farmers, along with an online website that regularly updated the status of LSD outbreaks in the country. Despite these efforts, there was a continuous rise in new LSD cases at the national level (23). In June 2021, the government initiated a mass vaccination campaign with the primary goal of controlling nationwide LSD outbreaks (24, 25). After the implementation of mass vaccination, the number of new LSD cases appeared to continuously decrease (18). Nonetheless, there has been no study conducted to assess the impact of this implementation.
Interrupted time series (ITS) analysis is increasingly used to evaluate the effectiveness of interventions, ranging from clinical treatments to broad public health policies (26–29). This analytical method is commonly employed to assess the effectiveness of large-scale population healthcare interventions. For example, ITS has been applied to evaluate large-scale health interventions initiated by the Australian government (27). Similarly, ITS has been instrumental in assessing the causal impact of lockdown measures on COVID-19 in both China (30) and Brazil (31) as well as the influence of social distancing on the COVID-19 epidemic in the United States (32) and Europe (33). Notably, the term “causal impact” is frequently mentioned in research publications that use ITS analysis to determine the effects of specific interventions (34–36).
Various methods of ITS have been proposed such as ordinary least square regression, generalized least square, restricted maximum likelihood, and autoregressive integrated moving average (29, 37–40). Within a Bayesian framework, Bayesian structural time series (BSTS) model is widely used. This quasi-experimental model is designed to assess the impact of an intervention by predicting a counterfactual time series, which represents the anticipated outcome in the absence of the intervention. The counterfactual is determined through an evaluation of the time series data. Ultimately, the assessment of impact involves comparing the observed time series during the intervention period to the counterfactual scenario generated by the model, thereby determining the difference (41–43). This approach has been adopted in multiple studies (35, 43–45). A distinct advantage of BSTS over alternative models is its counterfactual approach, while other models predominantly aim to identify overarching trend changes (46). Thus, this study attempted to determine the impact of the mass vaccination intervention using the BSTS method. In this regard, we estimated the number of new LSD cases under the premise that no mass vaccination campaign took place. These estimates were then compared with the actual figures from scenarios where mass vaccinations were conducted.
While many studies in Asia provide valuable insights into the epidemiology of LSD, there is a limited number of research focusing on large-scale national efforts to control LSD outbreaks in this region. Hence, this study aims to determine the effectiveness of the vaccination campaign implemented to control LSD outbreaks in Thailand using the BSTS approach.
2 Materials and methods
2.1 LSD outbreaks and intervention policy
In 2021, LSD outbreaks were reported in numerous areas across Thailand by the livestock authorities of the Department of Livestock Development (DLD). The outbreaks were confirmed through clinical diagnosis and by testing samples from LSD-affected animals using the PCR technique at the DLD laboratory, as described in previous studies (17–19). Figure 1 illustrates the LSD outbreak areas in Thailand, indicating a widespread occurrence throughout the country. Additionally, the number of new LSD cases from April to November 2021 is illustrated in Figure 2. For additional details, the graph presenting the number of new LSD cases on a daily basis is included in Supplementary Figure S1.
Figure 1. Outbreak areas of lumpy skin disease in Thailand in 2021 including the outbreaks from April to December (A), April to June (B), and July to December (C).
Figure 2. Number of new LSD cases on a weekly basis. The light-yellow line depicts the period of the mass vaccination campaign. Based on two scenarios, the blue line indicates the pre-intervention period, and the green line represents the post-intervention period.
This study utilized the reporting of the number of LSD-affected cattle from the national livestock system organized by DLD. Initially, provincial livestock officers provided the reports, having received information from district-level officers investigating LSD outbreaks on farms. Subsequently, using the data collected by these livestock officers, the total number of new LSD cases was calculated. The dataset utilized in the present study includes the number of new LSD cases on a weekly basis, spanning from April to November 2021.
The mass vaccination campaign was implemented in June 2021. A total of 360,000 doses of live attenuated homologous vaccines were administered to cattle in high-risk areas for LSD outbreaks across the country (18). The effectiveness of the mass vaccination campaign, initiated in early June, couldn't be expected to manifest instantly. Implementing such a nationwide campaign involves multiple considerations to determine its effective date. Firstly, numerous administrative tasks, including paperwork, vaccine transportation, and identification of target farmers, must be completed to distribute the vast quantity of vaccine doses throughout the country. Secondly, time was allocated for the registration and verification of farm owners, as well as the actual vaccination procedures in different regions. Thirdly, post-vaccination, immediate immunity is not conferred to the cattle (7, 47). Research indicates that antibodies typically emerge within 15 days and peak around 30 days after vaccination (5). Taking these factors into consideration, we propose two start dates at which the majority of the vaccinated cattle population would attain sufficient immunity to the disease: 45 days (July 15) and 60 days (July 30) post-June 1, identified as scenario A and scenario B. Under scenario A, the pre-intervention period extended from April 1 to July 15, while the pre-intervention period for scenario B covered April to July 30. The post-intervention periods for scenario A and scenario B were July 15 to November 31 and August 1 to November 31, respectively (Figure 2).
It is crucial to note that control measures, such as restrictions on animal movement, closure of live cattle markets, insect control, disinfection on farms, active surveillance for LSD cases, and heightened public awareness, have been consistently enforced since April 1. These measures were applied throughout both the pre-intervention and post-intervention periods.
2.2 Bayesian structural time series models
The Bayesian methodologies based on STSM, termed as BSTS models has been proposed by Brodersen and Hauser (41). These BSTS models are state-space models tailored for time series data, which can be defined through a set of equations as described in (41, 48, 49):
where and ηt~N(0, Qt) are independent of all unknown data. Equation 1 serves as the observation equation that links the observed data yt to the latent state vector αt. Equation 2 acts as the state equation detailing the progression of the vector αt over time. The term yt denotes the observed value, Zt represents a d-dimensional output vector, and Tt is a d×d transition matrix. Rt is a d×q matrix, while εt signifies the observation error with a white noise variance αt. The term ηt is a q-dimensional system error with a q×q state-diffusion matrix Qt, where q ≤ d.
The BSTS approach is utilized to learn these parameters. Subsequently, the Markov chain Monte Carlo (MCMC) technique in conjunction with a Gibbs sampler, is used for posterior inference (50). The causal effect is estimated by contrasting the counterfactual (predicted) values with the actual LSD case numbers during the post-intervention phase (49). A comprehensive description of the statistical models and graphical procedures can be found in prior studies (41, 51).
The R statistical software, with the “Causallmpact” package (41), was utilized to determine the causal effect and assess the statistical significance of the intervention using the BSTS approach. The CausalImpact package provides an intuitive interface that conceals the typical intricacies linked with Bayesian analysis such as normalization, setting priors and estimation. It seamlessly integrates with the bsts package for the actual model fitting process. In the initial step, the BSTS model predicted the number of new LSD cases under the assumption that the mass vaccination campaign was not implemented (counterfactual). In the second step, the evaluation of the difference between the predicted and actual number of new LSD cases was conducted to quantify the impact of the mass vaccination. Subsequently, point effects, absolute effects, and relative effects, each accompanied by their 95% credible intervals (CrIs), were derived by comparing the predicted and actual death trends across 1,000 MCMC iterations, which is a default value.
The results from the BSTS process in the “CausalImpact” package analysis provide three graphical representations: original, pointwise, and cumulative. The graph consists of actual observed values and the predicted values during the post-intervention period which represents the causal effect of the intervention. The difference the predicted values to the actual values corresponding to an impact of the intervention can be displayed per time-point and cumulative using pointwise and cumulative graphs. Furthermore, in the analysis, the BSTS model makes the assumption that the relationship between interventions and the time series, which was established during the pre-period, continues to be stable throughout the post-period (50).
3 Results
For scenario A (Table 1), the post-intervention period saw an average of 2,941 new LSD cases. In an event without intervention, this number was projected to be 25,970. The difference between the predicted and actual values amounts to −23,029 in terms of absolute effect, indicating a 119% (95% CrI: −121%, −38%) decrease in relative effect. Additionally, the cumulative count of new LSD cases stood at 58,811, whereas without the intervention, it was projected to be 519,401. This translates to a relative decrease of 119%. The likelihood of this effect occurring by mere chance is minimal, with a Bayesian one-sided tail-area probability of p = 0.04, indicating statistical significance.
Table 1. Summary of actual, predicted, absolute and relative effects from a Bayesian structural time series analysis based on scenario A.
For scenario B (Table 2), in the absence of mass vaccination, the predicted number of cases was 24,273. The actual cases were only 2,085 cases, signifying a relative effect of the intervention, which equates to a 78% (95% CrI: −126, −41%) decrease (p = 0.04).
Table 2. Summary of actual, predicted, absolute and relative effects from a Bayesian structural time series analysis based on scenario B.
Figures 3, 4 displays three panels that visually represent the results. The first panel showcases both the actual data and a counterfactual projection for the post-treatment phase. The second panel highlights the pointwise causal effect, evident from the difference between the observed data and the hypothetical forecasts. The third panel captures the cumulative impact of the intervention by aggregating the pointwise effects from the second panel.
Figure 3. Graphical representation based on a causal impact analysis pertaining to the implementation of a nationwide mass vaccination campaign aimed to control the outbreak of lumpy skin disease, as per scenario A. The graph includes the original data, pointwise and cumulative effects of the intervention.
Figure 4. Graphical representation based on the causal impact analysis pertaining to the implementation of nationwide mass vaccination campaign to control the outbreak of lumpy skin disease, as per scenario B. The graph includes the original data, pointwise and cumulative effects of the intervention.
From both scenario A (Figure 3) and scenario B (Figure 4), it is evident that the interventions influenced the number of new LSD cases. The top panel anticipated a consistent number of new LSD cases (represented by the dashed line) if no intervention took place. However, the actual figures (depicted by the solid line) reveal a decline in new LSD cases post-intervention. The middle panel's pointwise causal effect displays negative values, representing the divergence between the actual numbers and the counterfactual predictions. The cumulative influence of the intervention is delineated in the bottom panel.
4 Discussion
The introduction of control measures and a mass vaccination policy in Thailand represented a significant effort to curb the LSD epidemic. Through a causal impact analysis employing the BSTS methodology, this study underscores the efficacy of these interventions in reducing the number of LSD cases in Thailand.
Following the confirmation of the initial LSD outbreak in April (17), a series of control measures were rolled out to limit the spread of the disease. However, between May and June, there was a consistent rise in new LSD cases (18). This surge could be attributed to the vulnerability of naive herds lacking immunity to the disease, given that the LSD vaccine had never been used in Thailand before the discovery of the first LSD outbreak. The rapid spread of LSD could also be linked to challenges in controlling insect vectors (19, 20) and unauthorized animal movements (18).
In the context of an LSD outbreak, the most effective control measure involves widespread cattle immunization coupled with restrictions on animal movement (2). Predominantly, commercially available vaccines for LSD are live attenuated formulations derived from LSDV strains, sheep pox virus (SPPV), or goat pox virus (GTPV) (6, 16). Specifically, live attenuated LSDV vaccines, commonly referred to as Neethling vaccines, have been widely used in cattle across various countries. It is acknowledged that these vaccines exhibit notable efficacy, providing robust protection for cattle against virulent field LSDV strains (6). In line with the government control strategy, mass vaccination of live attenuated homologous LSD vaccines was conducted in multiple areas across Thailand. The mass vaccination was implemented following the ring vaccination strategy, with a focus on administering LSD vaccines to control outbreaks, especially in high-risk areas. This typically included neighboring regions within a 5–50 km radius of the LSD outbreak areas (18, 25). The effectiveness of mass vaccination efforts was notably evident when comparing actual values to projection values from the BSTS model, assuming no vaccination. The observed decrease in the number of cases by more than 75%, as predicted by the two scenarios, serves as compelling evidence of the success of mass vaccination initiatives. It is important to note that our findings align with prior studies highlighting the efficacy of mass vaccination in controlling LSD outbreaks. In South-Eastern Europe, a regional mass vaccination campaign using a homologous LSD vaccine covered over 1.8 million bovines, proving instrumental in managing LSD in that region (9). Similarly, the deployment of Neethling vaccines was recognized as an effective measure against LSD in Israel in 2012 (52) and in the Balkans from 2015 to 2017 (53).
A noteworthy control measure implemented by DLD was the establishment of a reporting system for LSD outbreaks. Local livestock authorities, including district and provincial livestock offices, actively participated by providing daily updates on the situation. Once these reports are reviewed and verified, DLD consolidates the data and publishes it on its website for public access. This online resource offers biweekly and monthly overviews of the LSD outbreak status in Thailand (https://sites.google.com/view/dldlsd/home). Through this mechanism, farmers and other stakeholders are consistently updated on the nation's LSD situation. Moreover, the proactive involvement of local livestock authorities is crucial in managing the LSD outbreak. They diligently identify LSD cases in herds by purposefully visiting farms in their respective localities (18, 54).
Following the nationwide mass vaccination, there was a noticeable decline in new LSD cases in Thailand toward the end of 2021 (22). However, some LSD cases were still reported in 2022 (55). One plausible reason for the persistence of a few new LSD cases could be the non-vaccination of certain cattle, such as the young ones or those in specific herds during the campaign. Additionally, farmers might have faced challenges in procuring vaccines for newborn animals post-campaign. Therefore, we advocate for the continuation of the vaccination process to ensure sustained immunity against LSD in Thai cattle. Given the limited research on the longevity of cattle immunity against LSD post-vaccination in Thailand, annual vaccination might be a prudent approach. This could be funded either by farmers and stakeholders directly or through government subsidies. It's imperative to achieve comprehensive vaccination coverage, and farmers should be motivated to persist with the vaccination regimen for multiple years. This should be overseen by livestock authorities and policymakers to prevent a decline in vaccination practices, as witnessed in Israel, which could potentially lead to a resurgence in LSD outbreaks (56). Strengthening LSD surveillance in cattle herds is also crucial. Livestock authorities should guide farmers on vigilant monitoring for new LSD cases within their herds. In the event of a suspected case, immediate notification to livestock authorities is essential, emphasizing the need for a swift response from these authorities to the farmers' concerns.
This study utilized univariate time series data for BSTS analysis, aligning with methodologies employed in previous studies (40, 43). In this approach, individual characteristics, such as the breed and age of animals, were not integrated into the BSTS models. Therefore, it is important to note that the estimations produced by the BSTS models in this study were solely contingent on the number of new LSD cases presented in the form of time series data. Furthermore, the interpretation of the results, which implies a decrease in the number of new LSD cases during the post-intervention period, should not be exclusively attributed to the mass vaccination campaign. Rather, it should be viewed as reflective of the impact of the mass vaccination campaign, taking into account the sustained implementation of other control measures both before and after the vaccination campaign.
This study comes with certain limitations. The reported number of LSD cases may underestimation of the actual count due to under-reporting issue. This could arise from herd owners not reporting outbreaks, occurrences in remote areas beyond regular surveillance, or certain LSD cases developing after the survey period. However, the impact of this limitation may not be significant, given the proactive approach in detecting LSD outbreaks by livestock authorities and the high level of cooperation from farmers. Furthermore, some farmers might purchase vaccines from private sectors and vaccinate their cattle independently. While the number of farmers engaging in this practice may not be high compared to those receiving government-provided vaccines at no cost, it is a situation that should be addressed.
5 Conclusion
This study represents the first application of the BSTS methodology to evaluate the impact of a government-led mass vaccination campaign against the LSD epidemic, drawing from national-level data. The methodologies employed in this research offer a deeper understanding of how established statistical tools can be used to determine the causal effects of intervention policies on large-scale LSD outbreaks. The findings from this investigation highlight the efficacy of the intervention strategy in mitigating LSD outbreaks. Notably, Thailand stands out as the only Asian country to have implemented a nationwide mass vaccination with extensive coverage. Consequently, the strategies adopted by Thai livestock authorities and stakeholders concerning mass vaccination and control measures could serve as valuable references for other countries at risk to a widespread LSD epidemic.
Data availability statement
The datasets presented in this article are not readily available because the dataset used in this research is sourced from the Department of Livestock Development (DLD) in Thailand. Requests to access the datasets should be directed to foreign@dld.go.th.
Ethics statement
This study is based on secondary data sourced from the Department of Livestock Development, Ministry of Agriculture and Cooperatives, Thailand. As our research exclusively employs pre-existing data and does not involve direct interactions with live animals, there is no requirement for ethical approval related to animal use.
Author contributions
VP: Conceptualization, Data curation, Formal analysis, Funding acquisition, Methodology, Project administration, Software, Supervision, Validation, Visualization, Writing—original draft, Writing—review & editing. OA: Data curation, Investigation, Methodology, Resources, Validation, Visualization, Writing—review & editing. NB: Data curation, Investigation, Methodology, Resources, Validation, Writing—review & editing. CJ: Data curation, Validation, Visualization, Writing—review & editing. RS: Methodology, Validation, Visualization, Writing—original draft, Writing—review & editing. KJ: Conceptualization, Methodology, Project administration, Software, Supervision, Writing—review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. The authors are grateful for research funding from Chiang Mai University (Grant: R66IN00356).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2023.1301546/full#supplementary-material
References
1. Sprygin A, Pestova Y, Wallace D, Tuppurainen E, Kononov A. Transmission of lumpy skin disease virus: a short review. Virus Res. (2019) 269:197637. doi: 10.1016/j.virusres.2019.05.015
2. Tuppurainen E, Oura C. Review: lumpy skin disease: an emerging threat to Europe, the Middle East and Asia. Transbound Emerg Dis. (2012) 59:40–8. doi: 10.1111/j.1865-1682.2011.01242.x
3. Parvin R, Chowdhury EH, Islam MT, Begum JA, Nooruzzaman M, Globig A, et al. Clinical epidemiology, pathology, and molecular investigation of lumpy skin disease outbreaks in bangladesh during 2020–2021 indicate the re-emergence of an old African strain. Viruses. (2022) 14:2529. doi: 10.3390/v14112529
4. Mathewos M, Dulo F, Tanga Z, Sombo M. Clinicopathological and molecular studies on cattle naturally infected with lumpy skin diseases in selected districts of Wolaita Zone, Southern Ethiopia. BMC Vet Res. (2022) 18:297. doi: 10.1186/s12917-022-03403-4
5. Tuppurainen E, Venter EH, Shisler J, Gari G, Mekonnen G, Juleff N, et al. Capripoxvirus diseases: current status and opportunities for control. Transbound Emerg Dis. (2017) 64:729–45. doi: 10.1111/tbed.12444
6. Tuppurainen E, Dietze K, Wolff J, Bergmann H, Beltran-Alcrudo D, Fahrion A, et al. Vaccines and vaccination against lumpy skin disease. Vaccines. (2021) 9:1136. doi: 10.3390/vaccines9101136
7. Haegeman A, De Leeuw I, Mostin L, Van Campe W, Philips W, Elharrak M, et al. Duration of immunity induced after vaccination of cattle with a live attenuated or inactivated lumpy skin disease virus vaccine. Microorganisms. (2023) 11:210. doi: 10.3390/microorganisms11010210
8. Nuvey FS, Arkoazi J, Hattendorf J, Mensah GI, Addo KK, Fink G, et al. Effectiveness and profitability of preventive veterinary interventions in controlling infectious diseases of ruminant livestock in sub-Saharan Africa: a scoping review. BMC Vet Res. (2022) 18:1–19. doi: 10.1186/s12917-022-03428-9
9. Authority EFS, Calistri P, De Clercq K, Gubbins S, Klement E, Stegeman A, et al. Lumpy skin disease epidemiological report IV: data collection and analysis. EFSA J. (2020) 18:e06010. doi: 10.2903/j.efsa.2020.6010
10. Uddin MA, Islam MA, Rahman AA, Rahman MM, Khasruzzaman A, Ward MP, et al. Epidemiological investigation of lumpy skin disease outbreaks in Bangladeshi cattle during 2019–2020. Transbound Emerg Dis. (2022) 69:3397–404. doi: 10.1111/tbed.14696
11. Chouhan CS, Parvin MS, Ali MY, Sadekuzzaman M, Chowdhury MGA, Ehsan MA, et al. Epidemiology and economic impact of lumpy skin disease of cattle in Mymensingh and Gaibandha districts of Bangladesh. Transbound Emerg Dis. (2022) 69:3405–18. doi: 10.1111/tbed.14697
12. Trinh TBN, Nguyen VT, Nguyen TTH, Mai NTA, Le PN, Lai TNH, et al. Molecular and histopathological characterization of lumpy skin disease in cattle in northern Vietnam during the 2020–2021 outbreaks. Arch Virol. (2022) 167:2143–9. doi: 10.1007/s00705-022-05533-4
13. Maw MT, Khin MM, Hadrill D, Meki IK, Settypalli TBK, Kyin MM, et al. First report of lumpy skin disease in Myanmar and molecular analysis of the field virus isolates. Microorganisms. (2022) 10:897. doi: 10.3390/microorganisms10050897
14. Ratyotha K, Prakobwong S, Piratae S. Lumpy skin disease: a newly emerging disease in Southeast Asia. Veter World. (2022) 15:2764. doi: 10.14202/vetworld.2022.2764-2771
15. Anwar A, Na-Lampang K, Preyavichyapugdee N, Punyapornwithaya V. Lumpy skin disease outbreaks in Africa, Europe, and Asia (2005–2022): multiple change point analysis and time series forecast. Viruses. (2022) 14:2203. doi: 10.3390/v14102203
16. Liang Z, Yao K, Wang S, Yin J, Ma X, Yin X, et al. Understanding the research advances on lumpy skin disease: a comprehensive literature review of experimental evidence. Front Microbiol. (2022) 13:1065894. doi: 10.3389/fmicb.2022.1065894
17. Arjkumpa O, Suwannaboon M, Boonrawd M, Punyawan I, Laobannu P, Yantaphan S, et al. First emergence of lumpy skin disease in cattle in Thailand, 2021. Transbound Emerg Dis. (2021) 68:3002–4. doi: 10.1111/tbed.14246
18. Suwankitwat N, Songkasupa T, Boonpornprasert P, Sripipattanakul P, Theerawatanasirikul S, Deemagarn T, et al. Rapid spread and genetic characterisation of a recently emerged recombinant lumpy skin disease virus in Thailand. Veter Sci. (2022) 9:542. doi: 10.3390/vetsci9100542
19. Arjkumpa O, Suwannaboon M, Boonrod M, Punyawan I, Liangchaisiri S, Laobannue P, et al. The first lumpy skin disease outbreak in Thailand 2021: epidemiological features and spatio-temporal analysis. Front Veter Sci. (2022) 8:1580. doi: 10.3389/fvets.2021.799065
20. Punyapornwithaya V, Seesupa S, Phuykhamsingha S, Arjkumpa O, Sansamur C, Jarassaeng C. Spatio-temporal patterns of lumpy skin disease outbreaks in dairy farms in northeastern Thailand. Front Veter Sci. (2022) 9:957306. doi: 10.3389/fvets.2022.957306
21. Singhla T, Boonsri K, Kreausukon K, Modethed W, Pringproa K, Sthitmatee N, et al. Molecular characterization and phylogenetic analysis of lumpy skin disease virus collected from outbreaks in northern Thailand in 2021. Veter Sci. (2022) 9:194. doi: 10.3390/vetsci9040194
22. Moonchai S, Himakalasa A, Rojsiraphisal T, Arjkumpa O, Panyasomboonying P, Kuatako N, et al. Modelling epidemic growth models for lumpy skin disease cases in Thailand using nationwide outbreak data, 2021–2022. Infect Dis Model. (2023) 8:282–93. doi: 10.1016/j.idm.2023.02.004
23. Punyapornwithaya V, Arjkumpa O, Buamithup N, Kuatako N, Klaharn K, Sansamur C, et al. Forecasting of daily new lumpy skin disease cases in Thailand at different stages of the epidemic using fuzzy logic time series, NNAR, and ARIMA methods. Prevent Veter Med. (2023) 217:105964. doi: 10.1016/j.prevetmed.2023.105964
24. Sariya L, Paungpin W, Chaiwattanarungruengpaisan S, Thongdee M, Nakthong C, Jitwongwai A, et al. Molecular detection and characterization of lumpy skin disease viruses from outbreaks in Thailand in 2021. Transbound Emerg Dis. (2022) 69:e2145–52. doi: 10.1111/tbed.14552
25. Premashthira S. Lumpy skin disease (LSD): an update on LSD prevention & control with focus on LSD vaccination. In: Lumpy Skin Disease Update Meeting for Southeast Asia. (2021).
26. Bernal JL, Cummins S, Gasparrini A. Interrupted time series regression for the evaluation of public health interventions: a tutorial. Int J Epidemiol. (2017) 46:348–55. doi: 10.1093/ije/dyw098
27. Schaffer AL, Dobbins TA, Pearson S-A. Interrupted time series analysis using autoregressive integrated moving average (ARIMA) models: a guide for evaluating large-scale health interventions. BMC Med Res Methodol. (2021) 21:1–12. doi: 10.1186/s12874-021-01235-8
28. Gebski V, Ellingson K, Edwards J, Jernigan J, Kleinbaum D. Modelling interrupted time series to evaluate prevention and control of infection in healthcare. Epidemiol Infect. (2012) 140:2131–41. doi: 10.1017/S0950268812000179
29. Turner SL, Karahalios A, Forbes AB, Taljaard M, Grimshaw JM, McKenzie JE. Comparison of six statistical methods for interrupted time series studies: empirical evaluation of 190 published series. BMC Med Res Methodol. (2021) 21:1–19. doi: 10.1186/s12874-021-01306-w
30. Figueiredo AM, Codina AD, Figueiredo D, Saez M, León AC. Impact of lockdown on COVID-19 incidence and mortality in China: an interrupted time series study. Bull World Health Organ. (2020) 6:2471. doi: 10.2471/BLT.20.256701
31. Silva L, Figueiredo Filho D, Fernandes A. The effect of lockdown on the COVID-19 epidemic in Brazil: evidence from an interrupted time series design. Cadernos de Saúde Pública. (2020) 36:e00213920. doi: 10.1590/0102-311x00213920
32. Doucette ML, Tucker A, Auguste ME, Watkins A, Green C, Pereira FE, et al. Initial impact of COVID-19′s stay-at-home order on motor vehicle traffic and crash patterns in Connecticut: an interrupted time series analysis. Injury Prevent. (2021) 27:3–9. doi: 10.1136/injuryprev-2020-043945
33. Vokó Z, Pitter JG. The effect of social distance measures on COVID-19 epidemics in Europe: an interrupted time series analysis. GeroScience. (2020) 42:1075–82. doi: 10.1007/s11357-020-00205-0
34. Navas Thorakkattle M, Farhin S, Khan AA. Forecasting the trends of covid-19 and causal impact of vaccines using bayesian structural time series and arima. Ann Data Sci. (2022) 9:1025–47. doi: 10.1007/s40745-022-00418-4
35. Feroze N. Forecasting the patterns of COVID-19 and causal impacts of lockdown in top five affected countries using Bayesian structural time series models. Chaos, Solit Fractals. (2020) 140:110196. doi: 10.1016/j.chaos.2020.110196
36. Sung H. Causal impacts of the COVID-19 pandemic on daily ridership of public bicycle sharing in Seoul. Sustain Cities Soc. (2023) 89:104344. doi: 10.1016/j.scs.2022.104344
37. Tweed EJ, Mackay DF, Boyd KA, Brown A, Byrne T, Conaglen P, et al. Evaluation of a national smoke-free prisons policy using medication dispensing: an interrupted time-series analysis. Lancet Public Health. (2021) 6:e795–804. doi: 10.1016/S2468-2667(21)00163-8
38. Jandoc R, Burden AM, Mamdani M, Lévesque LE, Cadarette SM. Interrupted time series analysis in drug utilization research is increasing: systematic review and recommendations. J Clin Epidemiol. (2015) 68:950–6. doi: 10.1016/j.jclinepi.2014.12.018
39. Hategeka C, Ruton H, Karamouzian M, Lynd LD, Law MR. Use of interrupted time series methods in the evaluation of health system quality improvement interventions: a methodological systematic review. BMJ Global Health. (2020) 5:e003567. doi: 10.1136/bmjgh-2020-003567
40. Li Y, Liu X, Li X, Xue C, Zhang B, Wang Y. Interruption time series analysis using autoregressive integrated moving average model: evaluating the impact of COVID-19 on the epidemic trend of gonorrhea in China. BMC Public Health. (2023) 23:2073. doi: 10.1186/s12889-023-16953-5
42. Barría-Sandoval C, Ferreira G, Lagos B, Bacigalupo CM. Assessing the effectiveness of quarantine measures during the COVID-19 pandemic in Chile using Bayesian structural time series models. Infect Dis Model. (2022) 7:625–36. doi: 10.1016/j.idm.2022.08.007
43. Hossain MS, Khan JR, Al Mamun SA, Islam MT, Raheem E. Excess mortality during the COVID-19 pandemic (2020–2021) in an urban community of Bangladesh. PLOS Global Public Health. (2023) 3:e0002176. doi: 10.1371/journal.pgph.0002176
44. Zhang Y, Fricker JD. Quantifying the impact of COVID-19 on non-motorized transportation: a Bayesian structural time series model. Transport policy. (2021) 103:11–20. doi: 10.1016/j.tranpol.2021.01.013
45. Pinilla J, Negrín M, González-López-Valcárcel B, Vázquez-Polo F-J. Using a Bayesian structural time–series model to infer the causal impact on cigarette sales of partial and total bans on public smoking. Jahrbücher Nationalökonomie Statistik. (2018) 238:423–39. doi: 10.1515/jbnst-2017-0125
46. Green M, Musi E, Rowe F, Charles D, Pollock FD, Kypridemos C, et al. Identifying how COVID-19-related misinformation reacts to the announcement of the UK national lockdown: an interrupted time-series study. Big Data Soc. (2021) 8:20539517211013869. doi: 10.1177/20539517211013869
47. Bazid A-H, Wasfy M, Fawzy M, Nayel M, Abdelmegeid M, Thabet RY, et al. Emergency vaccination of cattle against lumpy skin disease: Evaluation of safety, efficacy, and potency of MEVAC® LSD vaccine containing Neethling strain. Veter Res Commun. (2022) 47:767–777. doi: 10.1007/s11259-022-10037-2
48. Pontoh R, Zahroh S, Sunengsih N. Causal impact of PSBB and new normal in West Java. J Phys Conf Ser. (2021) 1722:012088. doi: 10.1088/1742-6596/1722/1/012088
49. Whata A, Chimedza C, A. Machine learning evaluation of the effects of South Africa's COVID-19 lockdown measures on population mobility. Mach Learn Knowl Extr. (2021) 3:481–506. doi: 10.3390/make3020025
50. De Vocht F, Tilling K, Pliakas T, Angus C, Egan M, Brennan A, et al. The intervention effect of local alcohol licensing policies on hospital admission and crime: a natural experiment using a novel Bayesian synthetictime-series method. J Epidemiol Commun Health. (2017) 71:912–8. doi: 10.1136/jech-2017-208931
51. Brodersen KH, Gallusser F, Koehler J, Remy N, Scott SL. Inferring causal impact using Bayesian structural time-series models. Ann Appl Stat. (2015) 9:247–274. doi: 10.1214/14-AOAS788
52. Ben-Gera J, Klement E, Khinich E, Stram Y, Shpigel N. Comparison of the efficacy of Neethling lumpy skin disease virus and x10RM65 sheep-pox live attenuated vaccines for the prevention of lumpy skin disease–The results of a randomized controlled field study. Vaccine. (2015) 33:4837–42. doi: 10.1016/j.vaccine.2015.07.071
53. Klement E, Broglia A, Antoniou S-E, Tsiamadis V, Plevraki E, Petrović T, et al. Neethling vaccine proved highly effective in controlling lumpy skin disease epidemics in the Balkans. Prev Vet Med. (2020) 181:104595. doi: 10.1016/j.prevetmed.2018.12.001
54. Modethed W, Singhla T, Boonsri K, Pringproa K, Sthitmatee N, Vinitchaikul P, et al. Identifying the patterns and sizes of the first lumpy skin disease outbreak clusters in Northern Thailand with a high degree of dairy farm aggregation using spatio-temporal models. PLoS ONE. (2023) 18:e0291692. doi: 10.1371/journal.pone.0291692
55. Suwankitwat N, Deemagarn T, Bhakha K, Songkasupa T, Lekcharoensuk P, Arunvipas P. Monitoring of genetic alterations of lumpy skin disease virus in cattle after vaccination in Thailand. J Animal Sci Technol. (2023). doi: 10.5187/jast.2023.e72
Keywords: lumpy skin disease, intervention, mass vaccination, Bayesian structural time series, causal impact, Thailand
Citation: Punyapornwithaya V, Arjkumpa O, Buamithup N, Jainonthee C, Salvador R and Jampachaisri K (2024) The impact of mass vaccination policy and control measures on lumpy skin disease cases in Thailand: insights from a Bayesian structural time series analysis. Front. Vet. Sci. 10:1301546. doi: 10.3389/fvets.2023.1301546
Received: 25 September 2023; Accepted: 11 December 2023;
Published: 05 January 2024.
Edited by:
Heinzpeter Schwermer, Federal Food Safety and Veterinary Office (FSVO), SwitzerlandReviewed by:
Mian Muhammad Awais, Bahauddin Zakariya University, PakistanKashif Hussain, Cholistan University of Veterinary and Animal Sciences, Pakistan
Copyright © 2024 Punyapornwithaya, Arjkumpa, Buamithup, Jainonthee, Salvador and Jampachaisri. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Katechan Jampachaisri, katechanj@nu.ac.th