Data-driven network modelling of disease transmission using complete population movement data: spread of VTEC O157 in Swedish cattle

The aquaculture disease network model (AquaNet-Mod): A simulation model to evaluate disease spread and controls for the salmonid industry in England and Wales

Infectious disease causes significant mortality in wild and farmed systems, threatening biodiversity, conservation and animal welfare, as well as food security. To mitigate impacts and inform policy, tools such as mathematical models and computer simulations are valuable for predicting the potential spread and impact of disease. This paper describes the development of the Aquaculture Disease Network Model, AquaNet-Mod, and demonstrates its application to evaluating disease epidemics and the efficacy of control, using a Viral Haemorrhagic Septicaemia (VHS) case study. AquaNet-Mod is a data-driven, stochastic, state-transition model. Disease spread can occur via four different mechanisms, i) live fish movement, ii) river based, iii) short distance mechanical and iv) distance independent mechanical. Sites transit between three disease states: susceptible, clinically infected and subclinically infected. Disease spread can be interrupted by the application of disease mitigation measures and controls such as contact tracing, culling, fallowing and surveillance. Results from a VHS case study highlight the potential for VHS to spread to 96% of sites over a 10 year time horizon if no disease controls are applied. Epidemiological impact is significantly reduced when live fish movement restrictions are placed on the most connected sites and further still, when disease controls, representative of current disease control policy in England and Wales, are applied. The importance of specific disease control measures, particularly contact tracing and disease detection rate, are also highlighted. The merit of this model for evaluation of disease spread and the efficacy of controls, in the context of policy, along with potential for further application and development of the model, for example to include economic parameters, is discussed.

A novel method to jointly estimate transmission rate and decay rate parameters in environmental transmission models

In environmental transmission, pathogens transfer from one individual to another via the environment. It is a common transmission mechanism in a wide range of host-pathogen systems. Incorporating environmental transmission in dynamic transmission models is crucial for gauging the effect of interventions, as extrapolating model results to new situations is only valid when the mechanisms are modelled correctly. The challenge in environmental transmission models lies in not jointly identifiable parameters for pathogen shedding, decay, and transmission dynamics. To solve this unidentifiability issue, we present a stochastic environmental transmission model with a novel scaling method for shedding rate parameter and a novel estimation method that distinguishes transmission rate and decay rate parameters. The core of our scaling and estimation method is calculating exposure and relating exposure to infection risks. By scaling shedding rate parameter, we standardize exposure to pathogens contributed by one infectious individual present during one time interval to one. The standardized exposure leads to a standard definition of transmission rate parameter applicable to scenarios with different decay rate parameters. Hence, we unify direct transmission (large decay rate) and environmental transmission in a continuous manner. More importantly, our exposure-based estimation method can correctly estimate back the transmission rate and the decay rate parameters, while the commonly used trajectory-based method failed. The reason is that exposure-based method gives the correct weight to infection data from previous observation periods. The correct estimation from exposure-based method will lead to more reliable predictions of intervention impact. Using the effect of disinfection as an example, we show how incorrectly estimated parameters may lead to incorrect conclusions about the effectiveness of interventions. This illustrates the importance of correct estimation of transmission rate and decay rate parameters for extrapolating environmental transmission models and predicting intervention effects.

A large-scale epidemiological model of BoHV-1 spread in the Irish cattle population to support decision-making in conformity with the European Animal Health Law

We present a new modelling framework to address the evaluation of national control/surveillance programs planned in line with the European Animal Health Law (AHL) for livestock diseases. Our modelling framework is applied to the cattle sector in Ireland where there is need for policy support to design an optimal programme to achieve bovine herpesvirus type 1 (BoHV-1) free status under the AHL. In this contribution, we show how our framework establishes a regional model that is able to mechanistically reproduce the demography, management practices and transport patterns of an entire cattle population without being dependent on continuous livestock registry data. An innovative feature of our model is the inclusion of herd typing, thereby extending these beyond the categories of dairy, beef and mixed herds that are frequently considered in other regional modelling studies. This detailed representation of herd type-specific management facilitates comparative assessment of BoHV-1 eradication strategies targeting different production types with individual strategy protocols. Finally, we apply our model to support current discussions regarding the structure and implementation of a potential national BoHV-1 eradication programme in Ireland.

Research perspectives on animal health in the era of artificial intelligence

Leveraging artificial intelligence (AI) approaches in animal health (AH) makes it possible to address highly complex issues such as those encountered in quantitative and predictive epidemiology, animal/human precision-based medicine, or to study host × pathogen interactions. AI may contribute (i) to diagnosis and disease case detection, (ii) to more reliable predictions and reduced errors, (iii) to representing more realistically complex biological systems and rendering computing codes more readable to non-computer scientists, (iv) to speeding-up decisions and improving accuracy in risk analyses, and (v) to better targeted interventions and anticipated negative effects. In turn, challenges in AH may stimulate AI research due to specificity of AH systems, data, constraints, and analytical objectives. Based on a literature review of scientific papers at the interface between AI and AH covering the period 2009-2019, and interviews with French researchers positioned at this interface, the present study explains the main AH areas where various AI approaches are currently mobilised, how it may contribute to renew AH research issues and remove methodological or conceptual barriers. After presenting the possible obstacles and levers, we propose several recommendations to better grasp the challenge represented by the AH/AI interface. With the development of several recent concepts promoting a global and multisectoral perspective in the field of health, AI should contribute to defract the different disciplines in AH towards more transversal and integrative research.

Modelling spread and surveillance of Mycobacterium avium subsp. paratuberculosis in the Swedish cattle trade network

To monitor a state of disease freedom and to ensure a timely detection of new introductions of disease, surveillance programmes need be evaluated prior to implementation. We present a strategy to evaluate surveillance of Mycobacterium avium subsp. paratuberculosis (MAP) using simulated testing of bulk milk in an infectious disease spread model. MAP is a globally distributed, chronic infectious disease with substantial animal health impact. Designing surveillance for this disease poses specific challenges because methods for surveillance evaluation have focused on estimating surveillance system sensitivity and probability of freedom from disease and do not account for spread of disease or complex and changing population structure over long periods. The aims of the study were to 1. define a model that describes the spread of MAP within and between Swedish herds; 2. define a method for simulation of imperfect diagnostic testing in this framework; 3. to compare surveillance strategies to support surveillance design choices. The results illustrate how this approach can be used to identify differences between the probability of detecting disease in the population based on choices of the number of herds sampled and the use of risk-based or random selection of these herds. The approach was also used to assess surveillance to detect introduction of disease and to detect a very low prevalence endemic state. The use of bulk milk sampling was determined to be an effective method to detect MAP in the population with as few as 500 herds tested per year if the herd-level prevalence was 0.2 %. However, detection of point introductions in the population was unlikely in the 13-year simulation period even if as many as 2000 herds were tested per year. Interestingly, the use of a risk-based selection strategy was found to be a disadvantage to detect MAP given the modelled disease dynamics.

Modelling infectious viral diseases in swine populations: a state of the art

Mathematical modelling is nowadays a pivotal tool for infectious diseases studies, completing regular biological investigations. The rapid growth of computer technology allowed for development of computational tools to address biological issues that could not be unravelled in the past. The global understanding of viral disease dynamics requires to account for all interactions at all levels, from within-host to between-herd, to have all the keys for development of control measures. A literature review was performed to disentangle modelling frameworks according to their major objectives and methodologies. One hundred and seventeen articles published between 1994 and 2020 were found to meet our inclusion criteria, which were defined to target papers representative of studies dealing with models of viral infection dynamics in pigs. A first descriptive analysis, using bibliometric indexes, permitted to identify keywords strongly related to the study scopes. Modelling studies were focused on particular infectious agents, with a shared objective: to better understand the viral dynamics for appropriate control measure adaptation. In a second step, selected papers were analysed to disentangle the modelling structures according to the objectives of the studies. The system representation was highly dependent on the nature of the pathogens. Enzootic viruses, such as swine influenza or porcine reproductive and respiratory syndrome, were generally investigated at the herd scale to analyse the impact of husbandry practices and prophylactic measures on infection dynamics. Epizootic agents (classical swine fever, foot-and-mouth disease or African swine fever viruses) were mostly studied using spatio-temporal simulation tools, to investigate the efficiency of surveillance and control protocols, which are predetermined for regulated diseases. A huge effort was made on model parameterization through the development of specific studies and methodologies insuring the robustness of parameter values to feed simulation tools. Integrative modelling frameworks, from within-host to spatio-temporal models, is clearly on the way. This would allow to capture the complexity of individual biological variabilities and to assess their consequences on the whole system at the population level. This would offer the opportunity to test and evaluate in silico the efficiency of possible control measures targeting specific epidemiological units, from hosts to herds, either individually or through their contact networks. Such decision support tools represent a strength for stakeholders to help mitigating infectious diseases dynamics and limiting economic consequences.

Data-Driven Network Modeling as a Framework to Evaluate the Transmission of Piscine Myocarditis Virus (PMCV) in the Irish Farmed Atlantic Salmon Population and the Impact of Different Mitigation Measures

Cardiomyopathy syndrome (CMS) is a severe cardiac disease of Atlantic salmon caused by the piscine myocarditis virus (PMCV), which was first reported in Ireland in 2012. In this paper, we describe the use of data-driven network modeling as a framework to evaluate the transmission of PMCV in the Irish farmed Atlantic salmon population and the impact of different mitigation measures. Input data included live fish movement data from 2009 to 2017, population dynamics events and the spatial location of the farms. With these inputs, we fitted a network-based stochastic infection spread model. After assumed initial introduction of the agent in 2009, our results indicate that it took 5 years to reach a between-farm prevalence of 100% in late 2014, with older fish being most affected. Local spread accounted for only a small proportion of new infections, being more important for sustained infection in a given area. Spread via movement of subclinically infected fish was most important for explaining the observed countrywide spread of the agent. Of the targeted intervention strategies evaluated, the most effective were those that target those fish farms in Ireland that can be considered the most connected, based on the number of farm-to-farm linkages in a specific time period through outward fish movements. The application of these interventions in a proactive way (before the first reported outbreak of the disease in 2012), assuming an active testing of fish consignments to and from the top 8 ranked farms in terms of outward fish movement, would have yielded the most protection for the Irish salmon farming industry. Using this approach, the between-farm PMCV prevalence never exceeded 20% throughout the simulation time (as opposed to the simulated 100% when no interventions are applied). We argue that the Irish salmon farming industry would benefit from this approach in the future, as it would help in early detection and prevention of the spread of viral agents currently exotic to the country.

A between-herd data-driven stochastic model to explore the spatio-temporal spread of hepatitis E virus in the French pig production network

Hepatitis E virus is a zoonotic pathogen for which pigs are recognized as the major reservoir in industrialised countries. A multiscale model was developed to assess the HEV transmission and persistence pattern in the pig production sector through an integrative approach taking into account within-farm dynamics and animal movements based on actual data. Within-farm dynamics included both demographic and epidemiological processes. Direct contact and environmental transmission routes were considered along with the possible co-infection with immunomodulating viruses (IMVs) known to modify HEV infection dynamics. Movements were limited to 3,017 herds forming the largest community on the swine commercial network in France and data from the national pig movement database were used to build the contact matrix. Between-herd transmission was modelled by coupling within-herd and network dynamics using the SimInf package. Different introduction scenarios were tested as well as a decrease in the prevalence of IMV-infected farms. After introduction of a single infected gilt, the model showed that the transmission pathway as well as the prevalence of HEV-infected pigs at slaughter age were affected by the type of the index farm, the health status of the population and the type of the infected farms. These outcomes could help design HEV control strategies at a territorial scale based on the assessment of the farms' and network's risk.

Bayesian epidemiological modeling over high-resolution network data

Mathematical epidemiological models have a broad use, including both qualitative and quantitative applications. With the increasing availability of data, large-scale quantitative disease spread models can nowadays be formulated. Such models have a great potential, e.g., in risk assessments in public health. Their main challenge is model parameterization given surveillance data, a problem which often limits their practical usage. We offer a solution to this problem by developing a Bayesian methodology suitable to epidemiological models driven by network data. The greatest difficulty in obtaining a concentrated parameter posterior is the quality of surveillance data; disease measurements are often scarce and carry little information about the parameters. The often overlooked problem of the model's identifiability therefore needs to be addressed, and we do so using a hierarchy of increasingly realistic known truth experiments. Our proposed Bayesian approach performs convincingly across all our synthetic tests. From pathogen measurements of shiga toxin-producing Escherichia coli O157 in Swedish cattle, we are able to produce an accurate statistical model of first-principles confronted with data. Within this model we explore the potential of a Bayesian public health framework by assessing the efficiency of disease detection and -intervention scenarios.

Reviewing age-structured epidemiological models of cattle diseases tailored to support management decisions: Guidance for the future

Mechanistic simulation models are being increasingly used as tools to assist with animal health decision-making in the cattle sector. We reviewed scientific literature for studies reporting age-structured cattle management models in application to infectious diseases. Our emphasis was on papers dedicated to support decision making in the field. In this systematic review we considered 1290 manuscripts and identified 76 eligible studies. These are based on 52 individual models from 10 countries addressing 9 different pathogens. We provide an overview of these models and present in detail their theoretical foundations, design paradigms and incorporated processes. We propose a structure of the characteristics of cattle disease models using three main features: [1] biological processes, [2] farming-related processes and [3] pathogen-related processes. It would be of benefit if future cattle disease models were to follow this structure to facilitate science communication and to allow increased model transparency.

Assessing biases in phylodynamic inferences in the presence of super-spreaders

Phylodynamic analyses using pathogen genetic data have become popular for making epidemiological inferences. However, many methods assume that the underlying host population follows homogenous mixing patterns. Nevertheless, in real disease outbreaks, a small number of individuals infect a disproportionately large number of others (super-spreaders). Our objective was to quantify the degree of bias in estimating the epidemic starting date in the presence of super-spreaders using different sample selection strategies. We simulated 100 epidemics of a hypothetical pathogen (fast evolving foot and mouth disease virus-like) over a real livestock movement network allowing the genetic mutations in pathogen sequence. Genetic sequences were sampled serially over the epidemic, which were then used to estimate the epidemic starting date using Extended Bayesian Coalescent Skyline plot (EBSP) and Birth–death skyline plot (BDSKY) models. Our results showed that the degree of bias varies over different epidemic situations, with substantial overestimations on the epidemic duration occurring in some occasions. While the accuracy and precision of BDSKY were deteriorated when a super-spreader generated a larger proportion of secondary cases, those of EBSP were deteriorated when epidemics were shorter. The accuracies of the inference were similar irrespective of whether the analysis used all sampled sequences or only a subset of them, although the former required substantially longer computational times. When phylodynamic analyses need to be performed under a time constraint to inform policy makers, we suggest multiple phylodynamics models to be used simultaneously for a subset of data to ascertain the robustness of inferences.

Neighbourhood contacts and trade movements drive the regional spread of bovine viral diarrhoea virus (BVDV)

To explore the regional spread of endemic pathogens, investigations are required both at within and between population levels. The bovine viral diarrhoea virus (BVDV) is such a pathogen, spreading among cattle herds mainly due to trade movements and neighbourhood contacts, and causing an endemic disease with economic consequences. To assess the contribution of both transmission routes on BVDV regional and local spread, we developed an original epidemiological model combining data-driven and mechanistic approaches, accounting for heterogeneous within-herd dynamics, animal movements and neighbourhood contacts. Extensive simulations were performed over 9 years in an endemic context in a French region with high cattle density. The most uncertain model parameters were calibrated on summary statistics of epidemiological data, highlighting that neighbourhood contacts and within-herd transmission should be high. We showed that neighbourhood contacts and trade movements complementarily contribute to BVDV spread on a regional scale in endemically infected and densely populated areas, leading to intense fade-out/colonization events: neighbourhood contacts generate the vast majority of outbreaks (72%) but mostly in low immunity herds and correlated to a rather short presence of persistently infected animals (P); trade movements generate fewer infections but could affect herds with higher immunity and generate a prolonged presence of P. Both movements and neighbourhood contacts should be considered when designing control or eradication strategies for densely populated region.

Spatio-temporal modelling of verotoxigenic Escherichia coli O157 in cattle in Sweden: exploring options for control

A spatial data-driven stochastic model was developed to explore the spread of verotoxigenic Escherichia coli O157 (VTEC O157) by livestock movements and local transmission among neighbouring holdings in the complete Swedish cattle population. Livestock data were incorporated to model the time-varying contact network between holdings and population demographics. Furthermore, meteorological data with the average temperature at the geographical location of each holding was used to incorporate season. The model was fitted against observed data and extensive numerical experiments were conducted to investigate the model’s response to control strategies aimed at reducing shedding and susceptibility, as well as interventions informed by network measures. The results showed that including local spread and season improved agreement with prevalence studies. Also, control strategies aimed at reducing the average shedding rate were more efficient in reducing the VTEC O157 prevalence than strategies based on network measures. The methodology presented in this study could provide a basis for developing disease surveillance on regional and national scales, where observed data are combined with readily available high-resolution data in simulations to get an overview of potential disease spread in unobserved regions.

Enterohaemorrhagic and other Shiga toxin-producing Escherichia coli (STEC): Where are we now regarding diagnostics and control strategies?

Escherichia coli comprises a highly diverse group of Gram-negative bacteria and is a common member of the intestinal microflora of humans and animals. Generally, such colonization is asymptomatic; however, some E. coli strains have evolved to become pathogenic and thus cause clinical disease in susceptible hosts. One pathotype, the Shiga toxigenic E. coli (STEC) comprising strains expressing a Shiga-like toxin is an important foodborne pathogen. A subset of STEC are the enterohaemorrhagic E. coli (EHEC), which can cause serious human disease, including haemolytic uraemic syndrome (HUS). The diagnosis of EHEC infections and the surveillance of STEC in the food chain and the environment require accurate, cost-effective and timely tests. In this review, we describe and evaluate tests now in routine use, as well as upcoming test technologies for pathogen detection, including loop-mediated isothermal amplification (LAMP) and whole-genome sequencing (WGS). We have considered the need for improved diagnostic tools in current strategies for the control and prevention of these pathogens in humans, the food chain and the environment. We conclude that although significant progress has been made, STEC still remains an important zoonotic issue worldwide. Substantial reductions in the public health burden due to this infection will require a multipronged approach, including ongoing surveillance with high-resolution diagnostic techniques currently being developed and integrated into the routine investigations of public health laboratories. However, additional research requirements may be needed before such high-resolution diagnostic tools can be used to enable the development of appropriate interventions, such as vaccines and decontamination strategies.