Modelling H5N1 in Bangladesh across spatial scales: Model complexity and zoonotic transmission risk

EPINEST, an agent-based model to simulate epidemic dynamics in large-scale poultry production and distribution networks

The rapid intensification of poultry production raises important concerns about the associated risks of zoonotic infections. Here, we introduce EPINEST (EPIdemic NEtwork Simulation in poultry Transportation systems): an agent-based modelling framework designed to simulate pathogen transmission within realistic poultry production and distribution networks. We provide example applications to broiler production in Bangladesh, but the modular structure of the model allows for easy parameterization to suit specific countries and system configurations. Moreover, the framework enables the replication of a wide range of eco-epidemiological scenarios by incorporating diverse pathogen life-history traits, modes of transmission and interactions between multiple strains and/or pathogens. EPINEST was developed in the context of an interdisciplinary multi-centre study conducted in Bangladesh, India, Vietnam and Sri Lanka, and will facilitate the investigation of the spreading patterns of various health hazards such as avian influenza, Campylobacter, Salmonella and antimicrobial resistance in these countries. Furthermore, this modelling framework holds potential for broader application in veterinary epidemiology and One Health research, extending its relevance beyond poultry to encompass other livestock species and disease systems.

Mathematical modeling of the infectious spread and outbreak dynamics of avian influenza with seasonality transmission for chicken farms

A compartmental model with a time-varying contact rate, the seasonality effect, and its corresponding nonautonomous model are investigated. The model is developed based on the six compartments: susceptible, latent, infected, asymptomatic, treated, and recovered individuals. We determine the effective reproduction number for this nonautonomous system, and analytic discussion shows that at least one positive periodic solution exists for R0>1. The model is simulated using the RK-45 numerical method, and the parameter values for the model are taken from the available literature. From the numerical results, we observe that the degree of seasonality and vaccine efficacy significantly impact the amplitude of the epidemic curve. The latent-infected phase plane shows that periodic solutions exhibit a period-doubling bifurcation as the amplitude of seasonality increases. Finally, the model outcome was compared with the actual field data and found to be consistent.

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.

An exploration of within-herd dynamics of a transboundary livestock disease: A foot and mouth disease case study

Transboundary livestock diseases are a high priority for policy makers because of the serious economic burdens associated with infection. In order to make well informed preparedness and response plans, policy makers often utilize mathematical models to understand possible outcomes of different control strategies and outbreak scenarios. Many of these models focus on the transmission between herds and the overall trajectory of the outbreak. While the course of infection within herds has not been the focus of the majority of models, a thorough understanding of within-herd dynamics can provide valuable insight into a disease system by providing information on herd-level biological properties of the infection, which can be used to inform decision making in both endemic and outbreak settings and to inform larger between-herd models. In this study, we develop three stochastic simulation models to study within-herd foot and mouth disease dynamics and the implications of different empirical data-based assumptions about the timing of the onset of infectiousness and clinical signs. We also study the influence of herd size and the proportion of the herd that is initially infected on the outcome of the infection. We find that increasing herd size increases the duration of infectiousness and that the size of the herd plays a more significant role in determining this duration than the number of initially infected cattle in that herd. We also find that the assumptions made regarding the onset of infectiousness and clinical signs, which are based on contradictory empirical findings, can result in the predictions about when infection would be detectable differing by several days. Therefore, the disease progression used to characterize the course of infection in a single bovine host could have significant implications for determining when herds can be detected and subsequently controlled; the timing of which could influence the overall predicted trajectory of outbreaks.

Environmental factors and spatiotemporal distribution characteristics of the global outbreaks of the highly pathogenic avian influenza H5N1

The spread of highly pathogenic avian influenza H5N1 has posed a major threat to global public health. Understanding the spatiotemporal outbreak characteristics and environmental factors of H5N1 outbreaks is of great significance for the establishment of effective prevention and control systems. The time and location of H5N1 outbreaks in poultry and wild birds officially confirmed by the World Organization for Animal Health from 2005 to 2019 were collected. Spatial autocorrelation analysis and multidistance spatial agglomeration analysis methods were used to analyze the global outbreak sites of H5N1. Combined with remote sensing data, the correlation between H5N1 outbreaks and environmental factors was analyzed using binary logistic regression methods. We analyzed the correlation between the H5N1 outbreak and environmental factors and finally made a risk prediction for the global H5N1 outbreaks. The results show that the peak of the H5N1 outbreaks occurs in winter and spring. H5N1 outbreaks exhibit aggregation, and a weak aggregation phenomenon is noted on the scale close to 5000 km. Water distance, road distance, railway distance, wind speed, leaf area index (LAI), and specific humidity were protective factors for the outbreak of H5N1, and the odds ratio (OR) were 0.985, 0.989, 0.995, 0.717, 0.832, and 0.935, respectively. Temperature was a risk factor with an OR of 1.073. The significance of these ORs was greater than 95%. The global risk prediction map was obtained. Given that the novel coronavirus (COVID-19) is spreading globally, the methods and results of this study can provide a reference for studying the spread of COVID-19.

Risk Areas for Influenza A(H5) Environmental Contamination in Live Bird Markets, Dhaka, Bangladesh

We evaluated the presence of influenza A(H5) virus environmental contamination in live bird markets (LBMs) in Dhaka, Bangladesh. By using Bernoulli generalized linear models and multinomial logistic regression models, we quantified LBM-level factors associated with market work zone–specific influenza A(H5) virus contamination patterns. Results showed higher environmental contamination in LBMs that have wholesale and retail operations compared with retail-only markets (relative risk 0.69, 95% 0.51–0.93; p = 0.012) and in March compared with January (relative risk 2.07, 95% CI 1.44–2.96; p<0.001). Influenza A(H5) environmental contamination remains a public health problem in most LBMs in Dhaka, which underscores the need to implement enhanced biosecurity interventions in LBMs in Bangladesh.

Trade-off between local transmission and long-range dispersal drives infectious disease outbreak size in spatially structured populations

Transmission of infectious diseases between immobile hosts (e.g., plants, farms) is strongly dependent on the spatial distribution of hosts and the distance-dependent probability of transmission. As the interplay between these factors is poorly understood, we use spatial process and transmission modelling to investigate how epidemic size is shaped by host clustering and spatial range of transmission. We find that for a given degree of clustering and individual-level infectivity, the probability that an epidemic occurs after an introduction is generally higher if transmission is predominantly local. However, local transmission also impedes transfer of the infection to new clusters. A consequence is that the total number of infections is maximal if the range of transmission is intermediate. In highly clustered populations, the infection dynamics is strongly determined by the probability of transmission between clusters of hosts, whereby local clusters act as multiplier of infection. We show that in such populations, a metapopulation model sometimes provides a good approximation of the total epidemic size, using probabilities of local extinction, the final size of infections in local clusters, and probabilities of cluster-to-cluster transmission. As a real-world example we analyse the case of avian influenza transmission between poultry farms in the Netherlands.

Transmission of infectious diseases between immobile hosts depends on the transmission characteristics of the infection and on the spatial distribution of hosts. Examples include infectious diseases of plants that are spread by wind or via vectors (e.g., Asiatic citrus canker spread between citrus trees), diseases that are transmitted between local host populations (e.g., sylvatic plague transmitted between rodents living in burrows), diseases of production animals that are spread between farms (e.g., avian influenza in poultry transmitted from farm to farm). We use spatial transmission modelling to investigate how the total number of infections over the course of an epidemic is determined by host clustering and spatial range of transmission. We find that for a given degree of clustering and infectivity of hosts, the number of infections is maximal if the spatial range of transmission is intermediate. In highly clustered populations we show that epidemic size can be approximated by a metapopulation model, illustrating that in such populations the transmission dynamics is dominated by transmission between clusters of hosts.

Time-aggregated mobile phone mobility data are sufficient for modelling influenza spread: the case of Bangladesh

Human mobility plays a major role in the spatial dissemination of infectious diseases. We develop a spatio-temporal stochastic model for influenza-like disease spread based on estimates of human mobility. The model is informed by mobile phone mobility data collected in Bangladesh. We compare predictions of models informed by daily mobility data (reference) with that of models informed by time-averaged mobility data, and mobility model approximations. We find that the gravity model overestimates the spatial synchrony, while the radiation model underestimates the spatial synchrony. Using time-averaged mobility resulted in spatial spreading patterns comparable to the daily mobility model. We fit the model to 2014–2017 influenza data from sentinel hospitals in Bangladesh, using a sequential version of approximate Bayesian computation. We find a good agreement between our estimated model and the case data. We estimate transmissibility and regional spread of influenza in Bangladesh, which are useful for policy planning. Time-averaged mobility appears to be a good proxy for human mobility when modelling infectious diseases. This motivates a more general use of the time-averaged mobility, with important implications for future studies and outbreak control. Moreover, time-averaged mobility is subject to less privacy concerns than daily mobility, containing less temporal information on individual movements.

Approximate Bayesian Computation for infectious disease modelling

Approximate Bayesian Computation (ABC) techniques are a suite of model fitting methods which can be implemented without a using likelihood function. In order to use ABC in a time-efficient manner users must make several design decisions including how to code the ABC algorithm and the type of ABC algorithm to use. Furthermore, ABC relies on a number of user defined choices which can greatly effect the accuracy of estimation. Having a clear understanding of these factors in reducing computation time and improving accuracy allows users to make more informed decisions when planning analyses. In this paper, we present an introduction to ABC with a focus of application to infectious disease models. We present a tutorial on coding practice for ABC in R and three case studies to illustrate the application of ABC to infectious disease models.

The epidemic potential of avian influenza A (H7N9) virus in humans in mainland China: A two-stage risk analysis

Background

From 2013 to 2017, more than one thousand avian influenza A (H7N9) confirmed cases with hundreds of deaths were reported in mainland China. To identify priorities for epidemic prevention and control, a risk assessing framework for subnational variations is needed to define the epidemic potential of A (H7N9).

Methods

We established a consolidated two-stage framework that outlined the potential epidemic of H7N9 in humans: The Stage 1, index-case potential, used a Boosted Regression Trees model to assess population at risk due to spillover from poultry; the Stage 2, epidemic potential, synthesized the variables upon a framework of the Index for Risk Management to measure epidemic potential based on the probability of hazards and exposure, the vulnerability and coping capacity.

Results

Provinces in southern and eastern China, especially Jiangsu, Zhejiang, Guangzhou, have high index-case potential of human infected with A (H7N9), while northern coastal provinces and municipalities with low morbidity, i.e. Tianjin and Liaoning, have an increasing risk of A (H7N9) infection. Provinces in central China are likely to have high potential of epidemic due to the high vulnerability and the lack of coping capacity.

Conclusions

This study provides a unified risk assessment of A (H7N9) to detect the two-stage heterogeneity of epidemic potential among different provinces in mainland China, allowing proactively evaluate health preparedness at subnational levels to improve surveillance, diagnostic capabilities, and health promotion.

The role of case proximity in transmission of visceral leishmaniasis in a highly endemic village in Bangladesh

BACKGROUND

Visceral leishmaniasis (VL) is characterised by a high degree of spatial clustering at all scales, and this feature remains even with successful control measures. VL is targeted for elimination as a public health problem in the Indian subcontinent by 2020, and incidence has been falling rapidly since 2011. Current control is based on early diagnosis and treatment of clinical cases, and blanket indoor residual spraying of insecticide (IRS) in endemic villages to kill the sandfly vectors. Spatially targeting active case detection and/or IRS to higher risk areas would greatly reduce costs of control, but its effectiveness as a control strategy is unknown. The effectiveness depends on two key unknowns: how quickly transmission risk decreases with distance from a VL case and how much asymptomatically infected individuals contribute to transmission.

METHODOLOGY/PRINCIPAL FINDINGS

To estimate these key parameters, a spatiotemporal transmission model for VL was developed and fitted to geo-located epidemiological data on 2494 individuals from a highly endemic village in Mymensingh, Bangladesh. A Bayesian inference framework that could account for the unknown infection times of the VL cases, and missing symptom onset and recovery times, was developed to perform the parameter estimation. The parameter estimates obtained suggest that, in a highly endemic setting, VL risk decreases relatively quickly with distance from a case-halving within 90m-and that VL cases contribute significantly more to transmission than asymptomatic individuals.

CONCLUSIONS/SIGNIFICANCE

These results suggest that spatially-targeted interventions may be effective for limiting transmission. However, the extent to which spatial transmission patterns and the asymptomatic contribution vary with VL endemicity and over time is uncertain. In any event, interventions would need to be performed promptly and in a large radius (≥300m) around a new case to reduce transmission risk.

Dynamics of the 2004 avian influenza H5N1 outbreak in Thailand: The role of duck farming, sequential model fitting and control

The Highly Pathogenic Avian Influenza (HPAI) subtype H5N1 virus persists in many countries and has been circulating in poultry, wild birds. In addition, the virus has emerged in other species and frequent zoonotic spillover events indicate that there remains a significant risk to human health. It is crucial to understand the dynamics of the disease in the poultry industry to develop a more comprehensive knowledge of the risks of transmission and to establish a better distribution of resources when implementing control. In this paper, we develop a set of mathematical models that simulate the spread of HPAI H5N1 in the poultry industry in Thailand, utilising data from the 2004 epidemic. The model that incorporates the intensity of duck farming when assessing transmision risk provides the best fit to the spatiotemporal characteristics of the observed outbreak, implying that intensive duck farming drives transmission of HPAI in Thailand. We also extend our models using a sequential model fitting approach to explore the ability of the models to be used in “real time” during novel disease outbreaks. We conclude that, whilst predictions of epidemic size are estimated poorly in the early stages of disease outbreaks, the model can infer the preferred control policy that should be deployed to minimise the impact of the disease.

The impact of surveillance and control on highly pathogenic avian influenza outbreaks in poultry in Dhaka division, Bangladesh

In Bangladesh, the poultry industry is an economically and socially important sector, but it is persistently threatened by the effects of H5N1 highly pathogenic avian influenza. Thus, identifying the optimal control policy in response to an emerging disease outbreak is a key challenge for policy-makers. To inform this aim, a common approach is to carry out simulation studies comparing plausible strategies, while accounting for known capacity restrictions. In this study we perform simulations of a previously developed H5N1 influenza transmission model framework, fitted to two separate historical outbreaks, to assess specific control objectives related to the burden or duration of H5N1 outbreaks among poultry farms in the Dhaka division of Bangladesh. In particular, we explore the optimal implementation of ring culling, ring vaccination and active surveillance measures when presuming disease transmission predominately occurs from premises-to-premises, versus a setting requiring the inclusion of external factors. Additionally, we determine the sensitivity of the management actions under consideration to differing levels of capacity constraints and outbreaks with disparate transmission dynamics. While we find that reactive culling and vaccination policies should pay close attention to these factors to ensure intervention targeting is optimised, across multiple settings the top performing control action amongst those under consideration were targeted proactive surveillance schemes. Our findings may advise the type of control measure, plus its intensity, that could potentially be applied in the event of a developing outbreak of H5N1 amongst originally H5N1 virus-free commercially-reared poultry in the Dhaka division of Bangladesh.

Ongoing circulation of avian influenza H5N1 viruses in poultry pose a global public health risk and cause extensive damage to the livestock industry. One of several countries in South Asia gravely affected is Bangladesh, where the poultry industry is an economically and socially important sector. Identifying the optimal control response in anticipation of further outbreaks is therefore a key challenge for policy-makers. This study tested a series of culling, vaccination and active surveillance intervention actions, assessing specific control objectives related to the burden or duration of H5N1 outbreaks among commercial poultry farms in the Dhaka division. This assessment was achieved through performing computational simulations of a previously developed H5N1 influenza transmission mathematical model. The findings of this assessment indicate that the impact of reactive culling and vaccination control policies are dependent upon transmission characteristics, control objectives and availability of resources to enact the control action, whereas proactive surveillance schemes significantly outperform reactive surveillance procedures irrespective of these conditions.