Search

Abstract 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.

Abstract We analysed the interplay between palmiped farm density and the vulnerability of the poultry production system to highly pathogenic avian influenza (HPAI) H5N8. To do so, we used a spatially-explicit transmission model, which was calibrated to reproduce the observed spatio-temporal distribution of outbreaks in France during the 2016–2017 epidemic of HPAI. Six scenarios were investigated, in which the density of palmiped farms was decreased in the municipalities with the highest palmiped farm density. For each of the six scenarios, we first calculated the spatial distribution of the basic reproduction number (R 0), i.e. the expected number of farms a particular farm would be likely to infect, should all other farms be susceptible. We also ran in silico simulations of the adjusted model for each scenario to estimate epidemic sizes and time-varying effective reproduction numbers. We showed that reducing palmiped farm density in the densest municipalities decreased substantially the size of the areas with high R 0 values (> 1.5). In silico simulations suggested that reducing palmiped farm density, even slightly, in the densest municipalities was expected to decrease substantially the number of affected poultry farms and therefore provide benefits to the poultry sector as a whole. However, they also suggest that it would not have been sufficient, even in combination with the intervention measures implemented during the 2016–2017 epidemic, to completely prevent the virus from spreading. Therefore, the effectiveness of alternative structural preventive approaches now needs to be assessed, including flock size reduction and targeted vaccination.

Modeling of infectious diseases at the livestock-wildlife interface is a unique subset of mathematical modeling with many innate challenges. To ascertain the characteristics of the models used in these scenarios, a scoping review of the scientific literature was conducted. Fifty-six studies qualified for inclusion. Only 14 diseases at this interface have benefited from the utility of mathematical modeling, despite a far greater number of shared diseases. The most represented species combinations were cattle and badgers (for bovine tuberculosis, 14), and pigs and wild boar [for African (8) and classical (3) swine fever, and foot-and-mouth and disease (1)]. Assessing control strategies was the overwhelming primary research objective (27), with most studies examining control strategies applied to wildlife hosts and the effect on domestic hosts (10) or both wild and domestic hosts (5). In spatially-explicit models, while livestock species can often be represented through explicit and identifiable location data (such as farm, herd, or pasture locations), wildlife locations are often inferred using habitat suitability as a proxy. Though there are innate assumptions that may not be fully accurate when using habitat suitability to represent wildlife presence, especially for wildlife the parsimony principle plays a large role in modeling diseases at this interface, where parameters are difficult to document or require a high level of data for inference. Explaining observed transmission dynamics was another common model objective, though the relative contribution of involved species to epizootic propagation was only ascertained in a few models. More direct evidence of disease spill-over, as can be obtained through genomic approaches based on pathogen sequences, could be a useful complement to further inform such modeling. As computational and programmatic capabilities advance, the resolution of the models and data used in these models will likely be able to increase as well, with a potential goal being the linking of modern complex ecological models with the depth of dynamics responsible for pathogen transmission. Controlling diseases at this interface is a critical step toward improving both livestock and wildlife health, and mechanistic models are becoming increasingly used to explore the strategies needed to confront these diseases.

Abstract In late 2015, an epizootic of Highly Pathogenic Avian Influenza (H5Nx) was registered in Southwestern France, including more than 70 outbreaks in commercial poultry flocks. Phylogenetic analyses suggested local emergence of H5 viruses which differed from A/goose/Guangdong/1/1996 clade 2.3.4.4b lineage and shared a unique polybasic cleavage site in their hemagglutinin protein. The present work provides an overview of the pathobiological picture associated with this epizootic in naturally infected chickens, guinea fowls and ducks. Upon necropsy examination, selected tissues were sampled for histopathology, immunohistochemistry and quantitative Real Time Polymerase Chain Reaction. In Galliformes, HPAIVs infection manifested as severe acute systemic vasculitis and parenchymal necrosis and was associated with endothelial expression of viral antigen. In ducks, lesions were mild and infrequent, with sparse antigenic detection in respiratory and digestive mucosae and leukocytes. Tissue quantifications of viral antigen and RNA were higher in chickens and guinea fowls compared to duck. Subsequently, recombinant HA (rHA) was generated from a H5 HPAIV isolated from an infected duck to investigate its glycan-binding affinity for avian mucosae. Glycan-binding analysis revealed strong affinity of rHA for 3’Sialyl-LacNAc and low affinity for Sialyl-LewisX, consistent with a duck-adapted virus similar to A/Duck/Mongolia/54/2001 (H5N2). K222R and S227R mutations on rHA sequence shifted affinity towards Sialyl-LewisX and led to an increased affinity for chicken mucosa, confirming the involvement of these two mutations in the glycan-binding specificity of the HA. Interestingly, the rHA glycan binding pattern of guinea fowl appeared intermediate between duck and chicken. The present study presents a unique pathobiological description of the H5 HPAIVs outbreaks that occurred in 2015–2016 in Southwestern France.

Diagnostic tests are used to classify individual animals’ infection statuses. However, validating test performance in wild animals without gold standard tests is extremely challenging, and the issue is further complicated in chronic conditions where measured immune parameters vary over time. Here, we demonstrate the value of combining evidence from different diagnostic approaches to aid interpretation in the absence of gold standards, large sample sizes, and controlled environments. Over a two-year period, we sampled 268 free-living meerkats (Suricata suricatta) longitudinally for Mycobacterium suricattae (a causative agent of tuberculosis), using three ante-mortem diagnostic tests based on mycobacterial culture, and antigen-specific humoral and cell-mediated immune responses, interpreting results both independently and in combination. Post-mortem cultures confirmed M. suricattae infection in 22 animals, which had prior ante-mortem information, 59% (13/22) of which were test-positive on a parallel test interpretation (PTI) of the three ante-mortem diagnostic assays (95% confidence interval: 37–79%). A similar ability to detect infection, 65.7% (95% credible interval: 42.7–84.7%), was estimated using a Bayesian approach to examine PTI. Strong evidence was found for a near doubling of the hazard of death (Hazard Ratio 1.75, CI: 1.14–2.67, p = 0.01), associated with a positive PTI result, thus demonstrating that these test results are related to disease outcomes. For individual tests, small sample sizes led to wide confidence intervals, but replication of conclusions, using different methods, increased our confidence in these results. This study demonstrates that combining multiple methodologies to evaluate diagnostic tests in free-ranging wildlife populations can be a useful approach for exploiting such valuable datasets.

Risk factors are key epidemiological concepts that are used to explain disease distributions. Identifying disease risk factors is generally done by comparing the characteristics of diseased and non-diseased populations. However, imperfect disease detectability generates disease observations that do not necessarily represent accurately the true disease situation. In this study, we conducted an extensive simulation exercise to emphasize the impact of imperfect disease detection on the outcomes of logistic models when case reports are aggregated at a larger scale (e.g., diseased animals aggregated at farm level). We used a probabilistic framework to simulate both the disease distribution in herds and imperfect detectability of the infected animals in these herds. These simulations show that, under logistic models, true herd-level risk factors are generally correctly identified but their associated odds ratio are heavily underestimated as soon as the sensitivity of the detection is less than one. If the detectability of infected animals is not only imperfect but also heterogeneous between herds, the variables associated with the detection heterogeneity are likely to be incorrectly identified as risk factors. This probability of type I error increases with increasing heterogeneity of the detectability, and with decreasing sensitivity. Finally, the simulations highlighted that, when count data is available (e.g., number of infected animals in herds), they should not be reduced to a presence/absence dataset at the herd level (e.g., presence or not of at least one infected animal) but rather modeled directly using zero-inflated count models which are shown to be much less sensitive to imperfect detectability issues. In light of these simulations, we revisited the analysis of the French bovine abortion surveillance data, which has already been shown to be characterized by imperfect and heterogeneous abortion detectability. As expected, we found substantial differences between the quantitative outputs of the logistic model and those of the zero-inflated Poisson model. We conclude by strongly recommending that efforts should be made to account for, or at the very least discuss, imperfect disease detectability when assessing associations between putative risk factors and observed disease distributions, and advocate the use of zero-inflated count models if count data is available.

In September 2019, African swine fever (ASF) was reported in South Korea for the first time. Since then, more than 651 ASF cases in wild boars and 14 farm outbreaks have been notified in the country. Despite the efforts to eradicate ASF among wild boar populations, the number of reported ASF-positive wild boar carcasses have increased recently. The purpose of this study was to characterize the spatial distribution of ASF-positive wild boar carcasses to identify the risk factors associated with the presence and number of ASF-positive wild boar carcasses in the affected areas. Because surveillance efforts have substantially increased in early 2020, we divided the study into two periods (2 October 2019 to 19 January 2020, and 19 January to 28 April 2020) based on the number of reported cases and aggregated the number of reported ASF-positive carcasses into a regular grid of hexagons of 3-km diameter. To account for imperfect detection of positive carcasses, we adjusted spatial zero-inflated Poisson regression models to the number of ASF-positive wild boar carcasses per hexagon. During the first study period, proximity to North Korea was identified as the major risk factor for the presence of African swine fever virus. In addition, there were more positive carcasses reported in affected hexagons with high habitat suitability for wild boars, low heat load index (HLI), and high human density. During the second study period, proximity to an ASF-positive carcass reported during the first period was the only significant risk factor for the presence of ASF-positive carcasses. Additionally, low HLI and elevation were associated with an increased number of ASF-positive carcasses reported in the affected hexagons. Although the proportion of ASF-affected hexagons increased from 0.06 (95% credible interval (CrI): 0.05–0.07) to 0.09 (95% CrI: 0.08–0.10), the probability of reporting at least one positive carcass in ASF-affected hexagons increased from 0.49 (95% CrI: 0.41–0.57) to 0.73 (95% CrI: 0.66–0.81) between the two study periods. These results can be used to further advance risk-based surveillance strategies in the Republic of Korea.

Highly pathogenic avian influenza (HPAI) H5N1 virus has been circulating in Vietnam since 2003, whilst outbreaks of HPAI H5N6 virus are more recent, having only been reported since 2014. Although the spatial distribution of H5N1 outbreaks and risk factors for virus occurrence has been extensively studied, there have been no comparative studies for H5N6. Data collected through active surveillance of Vietnamese live bird markets (LBMs) between 2011 and 2015 were used to explore and compare the spatiotemporal distributions of H5N1- and H5N6-positive LBMs. Conditional autoregressive models were developed to quantify spatiotemporal associations between agroecological factors and the two HPAI strains using the same set of predictor variables. Unlike H5N1, which exhibited a strong north–south divide, with repeated occurrence in the extreme south of a cluster of high-risk provinces, H5N6 was homogeneously distributed throughout Vietnam. Similarly, different agroecological factors were associated with each strain. Sample collection in the months of January and February and higher average maximum temperature were associated with higher likelihood of H5N1-positive market-day status. The likelihood of market days being positive for H5N6 increased with decreased river density, and with successive Rounds of data collection. This study highlights marked differences in spatial patterns and risk factors for H5N1 and H5N6 in Vietnam, suggesting the need for tailored surveillance and control approaches.