Estimating spatial and temporal variations of the reproduction number for highly pathogenic avian influenza H5N1 epidemic in Thailand

A systematic review of mechanistic models used to study avian influenza virus transmission and control

The global spread of avian influenza A viruses in domestic birds is causing increasing socioeconomic devastation. Various mechanistic models have been developed to better understand avian influenza transmission and evaluate the effectiveness of control measures in mitigating the socioeconomic losses caused by these viruses. However, the results of models of avian influenza transmission and control have not yet been subject to a comprehensive review. Such a review could help inform policy makers and guide future modeling work. To help fill this gap, we conducted a systematic review of the mechanistic models that have been applied to field outbreaks. Our three objectives were to: (1) describe the type of models and their epidemiological context, (2) list estimates of commonly used parameters of low pathogenicity and highly pathogenic avian influenza transmission, and (3) review the characteristics of avian influenza transmission and the efficacy of control strategies according to the mechanistic models. We reviewed a total of 46 articles. Of these, 26 articles estimated parameters by fitting the model to data, one evaluated the effectiveness of control strategies, and 19 did both. Values of the between-individual reproduction number ranged widely: from 2.18 to 86 for highly pathogenic avian influenza viruses, and from 4.7 to 45.9 for low pathogenicity avian influenza viruses, depending on epidemiological settings, virus subtypes and host species. Other parameters, such as the durations of the latent and infectious periods, were often taken from the literature, limiting the models' potential insights. Concerning control strategies, many models evaluated culling (n = 15), while vaccination received less attention (n = 6). According to the articles reviewed, optimal control strategies varied between virus subtypes and local conditions, and depended on the overall objective of the intervention. For instance, vaccination was optimal when the objective was to limit the overall number of culled flocks. In contrast, pre-emptive culling was preferred for reducing the size and duration of an epidemic. Early implementation consistently improved the overall efficacy of interventions, highlighting the need for effective surveillance and epidemic preparedness.

A review of estimated transmission parameters for the spread of avian influenza viruses

Avian influenza poses an increasing problem in Europe and around the world. Simulation models are a useful tool to predict the spatiotemporal risk of avian influenza spread and evaluate appropriate control actions. To develop realistic simulation models, valid transmission parameters are critical. Here, we reviewed published estimates of the basic reproduction number (R₀ ), the latent period and the infectious period by virus type, pathogenicity, species, study type and poultry flock unit. We found a large variation in the parameter estimates, with highest R₀ estimates for H5N1 and H7N3 compared with other types; for low pathogenic avian influenza compared with high pathogenic avian influenza types; for ducks compared with other species; for estimates from field studies compared with experimental studies; and for within-flock estimates compared with between-flock estimates. Simulation models should reflect this observed variation so as to produce more reliable outputs and support decision-making. How to incorporate this information into simulation models remains a challenge.

The effectiveness of travel restriction measures in alleviating the COVID-19 epidemic: evidence from Shenzhen, China

With the expansion of the global novel coronavirus disease (COVID-19) pandemic, unprecedented interventions have been widely implemented in many countries, including China. In view of this scenario, this research aims to explore the effectiveness of population mobility restriction in alleviating epidemic transmission during different stages of the outbreak. Taking Shenzhen, a city with a large immigrant population in China, as a case study, the real-time reproduction number of COVID-19 is estimated by statistical methods to represent the dynamic spatiotemporal transmission pattern of COVID-19. Furthermore, migration data between Shenzhen and other provinces are collected to investigate the impact of nationwide population flow on near-real-time dynamic reproductive numbers. The results show that traffic flow control between populated cities has an inhibitory effect on urban transmission, but this effect is not significant in the late stage of the epidemic spread in China. This finding implies that the government should limit international and domestic population movement starting from the very early stage of the outbreak. This work confirms the effectiveness of travel restriction measures in the face of COVID-19 in China and provides new insight for densely populated cities in imposing intervention measures at various stages of the transmission cycle.

Estimation of the Basic Reproduction Numbers of the Subtypes H5N1, H5N8, and H5N6 During the Highly Pathogenic Avian Influenza Epidemic Spread Between Farms

It is important to understand pathogen transmissibility in a population to establish an effective disease prevention policy. The basic reproduction number (R 0) is an epidemiologic parameter for understanding the characterization of disease and its dynamics in a population. We aimed to estimate the R 0 of the highly pathogenic avian influenza (HPAI) subtypes H5N1, H5N8, and H5N6, which were associated with nine outbreaks in Korea between 2003 and 2018, to understand the epidemic transmission of each subtype. According to HPAI outbreak reports of the Animal and Plant Quarantine Agency, we estimated the generation time by calculating the time of infection between confirmed HPAI-positive farms. We constructed exponential growth and maximum likelihood (ML) models to estimate the basic reproduction number, which assumes the number of secondary cases infected by the index case. The Kruskal-Wallis test was used to analyze the epidemic statistics between subtypes. The estimated generation time of H5N1, H5N8, and H5N6 were 4.80 days [95% confidence interval (CI) 4.23-5.38] days, 7.58 (95% CI 6.63-8.46), and 5.09 days (95% CI 4.44-5.74), respectively. A pairwise comparison showed that the generation time of H5N8 was significantly longer than that of the subtype H5N1 (P = 0.04). Based on the ML model, R 0 was estimated as 1.69 (95% CI 1.48-2.39) for subtype H5N1, 1.60 (95%CI 0.97-2.23) for subtype H5N8, and 1.49 (95%CI 0.94-2.04) for subtype H5N6. We concluded that R 0 estimates may be associated with the poultry product system, climate, species specificity based on the HPAI virus subtype, and prevention policy. This study provides an insight on the transmission and dynamics patterns of various subtypes of HPAI occurring worldwide. Furthermore, the results are useful as scientific evidence for establishing a disease control policy.

Evaluation of strategies using simulation model to control a potential outbreak of highly pathogenic avian influenza among poultry farms in Central Luzon, Philippines

The Philippines confirmed its first epidemic of Highly Pathogenic Avian Influenza (HPAI) on August 11, 2017. It ended in November of 2017. Despite the successful management of the epidemic, reemergence is a continuous threat. The aim of this study was to conduct a mathematical model to assess the spatial transmission of HPAI among poultry farms in Central Luzon. Different control strategies and the current government protocol of 1 km radius pre-emptive culling (PEC) from infected farms were evaluated. The alternative strategies include 0.5km PEC, 1.5km PEC, 2 km PEC, 2.5 km PEC, and 3 km PEC, no pre-emptive culling (NPEC). The NPEC scenario was further modeled with a time of government notification set at 24hours, 48 hours, and 72 hours after the detection. Disease spread scenarios under each strategy were generated using an SEIR (susceptible-exposed-infectious-removed) stochastic model. A spatial transmission kernel was calculated and used to represent all potential routes of infection between farms. We assumed that the latent period occurs between 1–2 days, disease detection at 5–7 days post-infection, notification of authorities at 5–7 days post-detection and start of culling at 1–3 days post notification. The epidemic scenarios were compared based on the number of infected farms, the total number of culled farms, and the duration of the epidemic. Our results revealed that the current protocol is the most appropriate option compared with the other alternative interventions considered among farms with reproductive ratio (R_i) > 1. Shortening the culling radius to 0.5 km increased the duration of the epidemic. Further increase in the PEC zone decreased the duration of the epidemic but may not justify the increased number of farms to be culled. Nonetheless, the no-pre-emptive culling (NPEC) strategy can be an effective alternative to the current protocol if farm managers inform the government immediately within 24 hours of observation of the presence of HPAI in their farms. Moreover, if notification is made on days 1–3 after the detection, the scale and length of the outbreak have been significantly reduced. In conclusion, this study provided a comparison of various control measures for confronting the spread of HPAI infection using the simulation model. Policy makers can use this information to enhance the effectiveness of the current control strategy.

Age-Dependent Lethality in Ducks Caused by Highly Pathogenic H5N6 Avian Influenza Virus

Ducks show notably higher resistance to highly pathogenic avian influenza viruses as compared to chickens. Here, we studied the age-dependent susceptibility in ducks to the infections caused by highly pathogenic avian influenza viruses. We intranasally infected ducks aged 1, 2, 4, and 8 weeks with highly pathogenic H5N6 avian influenza viruses isolated in South Korea in 2016. All the 1-and 2-week-old ducks died after infection, 20% of 3-week-old ducks died, and from the ducks aged 4 and 8 weeks, all of them survived. We performed microarray analysis and quantitative real-time PCR using total RNA isolated from the lungs of infected 2- and 4-week-old ducks to determine the mechanism underlying the age-dependent susceptibility to highly pathogenic avian influenza virus. Limited genes were found to be differentially expressed between the lungs of 2- and 4-week-old ducks. Cell damage-related genes, such as CIDEA and ND2, and the immune response-related gene NR4A3 were notably induced in the lungs of infected 2-week-old ducks compared to those in the lungs of infected 4-week-old ducks.

Comparison of methods to Estimate Basic Reproduction Number (R 0) of influenza, Using Canada 2009 and 2017-18 A (H1N1) Data

Background

The basic reproduction number (R ₀) has a key role in epidemics and can be utilized for preventing epidemics. In this study, different methods are used for estimating R ₀'s and their vaccination coverage to find the formula with the best performance.

Materials and Methods

We estimated R ₀ for cumulative cases count data from April 18 to July 6, 2009 and 35-2017 to 34-2018 weeks in Canada: maximum likelihood (ML), exponential growth rate (EG), time-dependent reproduction numbers (TD), attack rate (AR), gamma-distributed generation time (GT), and the final size of the epidemic. Gamma distribution with mean and standard deviation 3.6 ± 1.4 is used as GT.

Results

The AR method obtained a R _{0 (}95% confidence interval [CI]) value of 1.116 (1.1163, 1.1165) and an EG (95%CI) value of 1.46 (1.41, 1.52). The R ₀ (95%CI) estimate was 1.42 (1.27, 1.57) for the obtained ML, 1.71 (1.12, 2.03) for the obtained TD, 1.49 (1.0, 1.97) for the gamma-distributed GT, and 1.00 (0.91, 1.09) for the final size of the epidemic. The minimum and maximum vaccination coverage were related to AR and TD methods, respectively, where the TD method has minimum mean squared error (MSE). Finally, the R ₀ (95%CI) for 2018 data was 1.52 (1.11, 1.94) by TD method, and vaccination coverage was estimated as 34.2%.

Conclusion

For the purposes of our study, the estimation of TD was the most useful tool for computing the R ₀, because it has the minimum MSE. The estimation R ₀ > 1 indicating that the epidemic has occurred. Thus, it is required to vaccinate at least 41.5% to prevent and control the next epidemic.

The implementation of cattle market closure strategies to mitigate the foot-and-mouth disease epidemics: A contact modeling approach

Foot-and-mouth disease (FMD) is one of the most endemic diseases in livestock worldwide. The disease occurrence generally results in a huge economic impact. The virus may distribute across countries or even continents along the contact network of animal movements. The present study, therefore, aimed to explore a cattle movement network originated in Tak, a Thailand-Myanmar bordered province and to demonstrate how FMDV spread among the nodes of market, source and destination. Subsequently, we examined the effectiveness of market closure intervention. The market-market (M-M) network was constructed to highlight the inter-market connections and the FMDV was modeled to spread along the trade chain. Four market closure scenarios based on rapidness and duration of implementation were examined. Our results indicate that two of the three major markets located in the province were highly connected and a strongly connected component was identified. The intra-provincial animal movements, which were currently overlooked, should be moved into sights as most of the high-risk sources for FMD epidemics were recognized in a close proximity to the cattle markets. Simultaneously, remote destinations across the country were identified. The inter-provincial animal movement control must be strengthened once FMD outbreak is notified. Based on our simulations, closing markets with low inter-market connectivity may not prevent the spread of FMDV. A selective market closure strategy targeting highly connected markets together with cattle trader tracking system was an alternative approach. However, socio-economic consequences regarding this intervention must be considered.

Estimating the between-farm transmission rates for highly pathogenic avian influenza subtype H5N1 epidemics in Bangladesh between 2007 and 2013

Highly Pathogenic Avian Influenza (HPAI) is classified by the World Organization for Animal Health as one of the notifiable diseases. Its occurrence is associated with severe socio-economic impacts and is also zoonotic. Bangladesh HPAI epidemic data for the period between 2007 and 2013 were obtained and split into epidemic waves based on the time lag between outbreaks. By assuming the number of newly infected farms to be binomially distributed, we fit a Generalized Linear Model to the data to estimate between-farm transmission rates (β). These parameters are then used together with the calculated infectious periods to estimate the respective basic reproduction numbers (R₀ ). The change in β and R₀ with time during the course of each epidemic wave was explored. Finally, sensitivity analyses of the effects of reducing the delay in detecting infection on a farm as well as extended infectiousness of a farm beyond the day of culling were assessed. The point estimates obtained for β ranged from 0.08 (95% CI: 0.06-0.10) to 0.11 (95% CI: 0.08-0.20) per infectious farm per day while R₀ ranged from 0.85 (95% CI: 0.77-1.02) to 0.96 (95% CI: 0.72-1.20). Sensitivity analyses reveal that the estimates are quite robust to changes in the assumptions about the day in reporting infection and extended infectiousness. In the analysis allowing for time-varying transmission parameters, the rising and declining phases observed in the epidemic data were synchronized with the moments when R₀ was greater and less than one, respectively. From an epidemiological perspective, the consistency of these estimates and their magnitude (R₀  ≈ 1) indicate that the effectiveness of the deployed control measures was largely invariant between epidemic waves and the trend of the time-varying R₀ supports the hypothesis of sustained farm-to-farm transmission that is possibly initiated by a few unique introductions.

Quantitative effects of a declaration of a state of emergency on foot-and-mouth disease

OBJECTIVES

The law in Japan requires the declaration of a state of emergency and implementation of countermeasures for an epidemic of a new infectious disease. However, because a state of emergency has never been declared in Japan, its effects remain unknown. The required countermeasures are similar to those implemented in the foot-and-mouth disease epidemic in Miyazaki in 2010. This study aimed to quantitatively estimate the effect of the declaration in 2010 and investigate the nature of the epidemic based on the day on which the declaration took effect.

METHODS

Only publicly available data were used. Data for farms in the most affected town were analyzed. A modified susceptible-infected-recovered model was used to estimate the effect and for the simulation. Another model was used to estimate the effective reproduction number.

RESULTS

After the declaration, the intra-bovine transmission rate decreased by 18.1 %, and there were few days when the effective reproduction number was >1.0. A few weeks delay in the declaration significantly increased the possibility of epidemic, number of farms at peak, and final infection scale.

CONCLUSIONS

Based on the substantial decrease in the transmission rate after the declaration of a state of emergency in 2010, a future declaration will have a similar effect for a new infectious disease even though a direct extrapolation is not valid. Although a declaration should be carefully considered owing to the potential socioeconomic effects, it is essential to prepare for the implementation given that a delay of only a few weeks should be acceptable.

Quantitative transmission characteristics of different H5 low pathogenic avian influenza viruses in Muscovy ducks

EU annual serosurveillance programs show that domestic duck flocks have the highest seroprevalence of H5 antibodies, demonstrating the circulation of notifiable avian influenza virus (AIV) according to OIE, likely low pathogenic (LP). Therefore, transmission characteristics of LPAIV within these flocks can help to understand virus circulation and possible risk of propagation. This study aimed at estimating transmission parameters of four H5 LPAIV (three field strains from French poultry and decoy ducks, and one clonal reverse-genetics strain derived from one of the former), using a SIR model to analyze data from experimental infections in SPF Muscovy ducks. The design was set up to accommodate rearing on wood shavings with a low density of 1.6 ducks/m(2): 10 inoculated ducks were housed together with 15 contact-exposed ducks. Infection was monitored by RNA detection on oropharyngeal and cloacal swabs using real-time RT-PCR with a cutoff corresponding to 2-7 EID50. Depending on the strain, the basic reproduction number (R0) varied from 5.5 to 42.7, confirming LPAIV could easily be transmitted to susceptible Muscovy ducks. The lowest R0 estimate was obtained for a H5N3 field strain, due to lower values of transmission rate and duration of infectious period, whereas reverse-genetics derived H5N1 strain had the highest R0. Frequency and intensity of clinical signs were also variable between strains, but apparently not associated with longer infectious periods. Further comparisons of quantitative transmission parameters may help to identify relevant viral genetic markers for early detection of potentially more virulent strains during surveillance of LPAIV.

Transmission rate and reproductive number of the H5N1 highly pathogenic avian influenza virus during the December 2005-July 2008 epidemic in Nigeria

We quantified the between-village transmission rate, β (the rate of transmission of H5N1 HPAI virus per effective contact), and the reproductive number, Re (the average number of outbreaks caused by one infectious village during its entire infectious period), of H5N1 highly pathogenic avian influenza (HPAI) virus in Nigeria using outbreak data collected between December 2005 and July 2008. We classified the outbreaks into two phases to assess the effectiveness of the control measures implemented. Phase 1 (December 2005-October 2006) represents the period when the Federal Government of Nigeria managed the HPAI surveillance and response measures, while Phase 2 (November 2006-July 2008) represents the time during which the Nigeria Avian Influenza Control and Human Pandemic Preparedness project (NAICP), funded by a World Bank credit of US$ 50 million, had taken over the management of most of the interventions. We used a total of 204 outbreaks from 176 villages that occurred in 78 local government areas of 25 states. The compartmental susceptible-infectious model was used as the analytical tool. Means and 95% percentile confidence intervals were obtained using bootstrapping techniques. The overall mean β (assuming a duration of infectiousness, T, of 12 days) was 0.07/day (95% percentile confidence interval: 0.06-0.09). The first and second phases of the epidemic had comparable β estimates of 0.06/day (0.04-0.09) and 0.08/day (0.06-0.10), respectively. The Re of the virus associated with these β and T estimates was 0.9 (0.7-1.1); the first and second phases of the epidemic had Re of 0.84 (0.5-1.2) and 0.9 (0.6-1.2), respectively. We conclude that the intervention measures implemented in the second phase of the epidemic had comparable effects to those implemented during the first phase and that the Re of the epidemic was low, indicating that the Nigeria H5N1 HPAI epidemic was unstable.