INTRODUCTION
African swine fever virus (ASFV), a member of the Asfarviridae family, causing a severe haemorrhagic disease in domestic and wild swine is currently circulating in Eastern Europe [1, 2]. Controlling the spread remains challenging due to the lack of effective vaccines and treatments, the potentially long survival of the virus in the environment and the potential role of wild boar in virus spread [1]. African swine fever (ASF) also represents a serious concern for European pig-producing countries due to the cost of stamping-out policies and the economic burden that bans impose on international trade, especially as production and consumption of pork are expected to further increase in future [1, Reference Madec3]. ASFV spreads by direct contact between infectious and susceptible animals [Reference de Carvalho4–Reference Eblé10]. Infection may also occur by indirect contact with, for example, contaminated surfaces of transport vehicles, contaminated clothing of animal workers or through feeding of swill that contains contaminated pig products [Reference Costard11–Reference Mur, Martinez-Lopez and Sanchez-Vizcaino13]. Airborne transmission has been demonstrated over short distances [Reference Wilkinson and Donaldson14]. For each of the potential routes for ASFV transmission between individual pigs, quantitative parameters have to be determined so that control programmes can be targeted based on the relative importance of each pathway.
The basic reproduction number (R 0), defined as the average number of newly infected cases caused by one infectious individual during its infectious period in a susceptible population, is a key quantitative indicator of the transmission potential and therefore very useful for evaluating the likely impact of different components in control strategies [Reference Anderson and May15–Reference Wonham, de-Camino-Beck and Lewis18]. So far, R 0 for the currently circulating ASFV strain in Eastern Europe has been reported in only one study [Reference Gulenkin19]. The authors developed a mathematical model to describe the transmission of ASFV within farms of the Russian Federation and estimated R 0 at 9·8 [95% confidence interval (CI) 3·9–15·6] based on field and experimental data. However, estimates for transmission parameters strongly depend on assumptions regarding the length of the infectious, latent and incubation periods [Reference Diekmann and Heesterbeek20–Reference Keeling and Rohani23]. For example, different R 0 estimates for classical swine fever virus (CSFV) were obtained depending on the diagnostic tool (detection of genome or infectious virus) or the biological marker (viraemia, infectious oropharyngeal or rectal fluids, clinical symptoms) used to determine the time of infection or the duration of infectiousness [Reference Gogin24, Reference de Carvalho25].
The objective of the current study was to estimate quantitative pig-to-pig transmission parameters for the currently circulating ASFV strain in Eastern Europe. For this purpose, we fitted stochastic Susceptible-Exposed-Infectious-Removed (SEIR) models to detailed data available from transmission experiments, in which susceptible pigs were exposed by direct and indirect contact to pigs infected with the ASFV Georgia 2007/1 strain. Furthermore, we assessed the sensitivity of the estimates to assumptions about the time of infection and duration of infectiousness for the infected animals.
Ethical standards
The authors assert that all procedures contributing to this work comply with the ethical standards of the relevant national and institutional guides on the care and use of laboratory animals. All experimental procedures, carried out under UK Home Office Licence number 70/7198, were approved by the Pirbright Ethical Review process. All animals were euthanized for animal welfare reasons.
METHODS
Experimental ASFV dataset
Data were obtained from experimental ASFV transmission in weaner pigs as described in Guinat et al. [Reference Guinat26]. In short, four groups of pigs were allocated to different rooms (A, B, C, D) and randomly selected pigs were inoculated intramuscularly with the highly virulent Georgia 2007/1 ASFV strain (Table 1). The inoculated pigs were in direct contact with susceptible pigs within each pen in all rooms. Rooms B and C were each split into two adjacent pens allowing indirect contact between inoculated and susceptible pigs, most likely through airborne transmission and small amounts of urine and faeces passing under the fence. The partition between pens did not allow nose-to-nose contacts. Blood samples were collected every 2 days and oral and nasal fluid samples daily. Samples were tested for presence of live ASFV and ASFV genome by virus isolation and quantitative real-time polymerase chain reaction (qPCR), respectively. All inoculated animals became infected and contact transmission was observed in all experimental groups. The inoculated pigs were generally viraemic 1–3 days before starting to excrete ASFV genome through the nasal and oral routes. The within-pen and between-pen contact pigs generally started to excrete ASFV genome through nasal and oral fluid samples 1–4 days before showing viraemia. Table 2 summarizes the diagnostic results in blood samples during the study period for all rooms.
Numbers of pigs used and transmission results with the Georgia 2007/1 African swine fever virus strain
Results of African swine fever virus (ASFV) isolation in blood samples from transmission experiments with domestic pigs infected with the Georgia 2007/1 ASFV strain
WP, Within-pen contact pigs; BP, between-pen contact pigs; –, negative virus isolation in blood; +, positive virus isolation in blood.
* Day of euthanasia.
Model assumptions
Three models were developed in this study, each assuming a different latent period (summarized in Table 3). In all models, the animals were assumed to be infectious as soon as ASFV was isolated from blood samples, supported by the fact that ASFV was mainly shed in blood [Reference Guinat26]. Since data on the moment of infection for the contact pigs were lacking, they were considered infected before the onset of infectiousness using an assumed latent period (L, the time from infection to the onset of infectiousness [Reference Vynnycky and White21]) of either 3 days (model 1), 4 days (model 2) or 5 days (model 3). This was based on values observed among the inoculated pigs (Table 2).
Model inputs
* A 2-day latent period had also been considered (data not shown).
The infectious period (T) represents the time between the onset of infectiousness and death or recovery [Reference Vynnycky and White21]. The transmission experiments did not allow T to be estimated because animals were euthanized for welfare reasons, resulting in censored values [± standard deviation (s.d.)] for T (e.g. average of 4·7 (±1·4), 2·4 (±0·7) and 3·0 (±1·2) days observed for the inoculated, within-pen and between-pen contact pigs, respectively [Reference Guinat26]). Therefore, animals were considered infectious either for a time period of 3–6 days (minimum infectious period) usually reported in previous studies or for 3–14 days (maximum infectious period) which has occasionally been observed [Reference Blome, Gabriel and Beer27–Reference Gabriel29]. We assumed T was normally distributed, as is considered to be a biologically realistic distribution for the infectious period of a virus [Reference Keeling and Grenfell30], either with a mean (±s.d.) of 4·5 (±0·75) days or 8·5 (±2·75) days.
The animals were assumed to mix randomly either within or between pens, i.e. to have equal opportunity for contact with each other [Reference Vynnycky and White21], which was considered plausible given the limited room/pen size in the experiments [Reference Guinat26]. The number of new infections per unit of time was assumed to be proportional to the product of the proportion (or frequency) of infectious animals and the number of susceptible animals [Reference Vynnycky and White21] which was considered adequate to describe the transmission dynamics process [Reference Bouma, De Jong and Kimman31, Reference McCallum, Barlow and Hone32].
Model analyses
We considered a standard transmission scenario with I A,t and I B,t , the number of infectious pigs, N A,t and N B,t the total population size, in two pens A and B at time t, respectively. Pen A housed the inoculated pigs and the within-pen contact pigs, while pen B housed the between-pen contact pigs. We used a stochastic SEIR model to estimate the within- and between-pen transmission rate parameters (β w and β b). In the model (Fig. 1), the population was divided into four classes over time t [Reference Vynnycky and White21]: pigs leaving the susceptible compartment (S t ) became infected (E A,t and E B,t ) with probability of infection p t , those leaving the infected compartment became infectious (I A,t and I B,t ) after L days, and those leaving the infectious compartment died (R t ) after T days. The probability p t depends on the number of infectious pigs (I A,t and I B,t ), the total population size (N A,t and N B,t ) and how frequently they make effective contact with each other within and between pens (β w and β b), so that
$$p_t = 1 - {\rm exp}\; \left( { - \left( {\displaystyle{{\beta _{\rm w} I_{{\rm B},t} \;} \over {N_{{\rm B},t}}} + \displaystyle{{\beta _{\rm b} I_{{\rm A},t}} \over {\left( {N_{{\rm A},t} + N_{{\rm B},t}} \right)}}} \right)} \right) [33].$$
Schematic representation of the SEIR model used for estimating the experimental pig-to-pig transmission parameters for African swine fever virus. Considering two adjacent pens A and B at time t, pigs leave the susceptible compartment (S t ) and become infected (E A,t or E B,t ) with probability p i,t . They leave the infected compartment and become infectious (I A,t or I B,t ) after L days. Finally, they leave the infectious compartment and die (R t ) after T days.
The number of new infections at time t (C t ) follows a binomial distribution with the number of trials given by the susceptible population (S t ) and the probability of success given by the probability of infection (p t ). This was used to construct the likelihood for the data and we estimated β w and β b by maximizing this likelihood [Reference Klinkenberg34, Reference Myung35]. In addition, 95% CIs were determined using ±1·96 times the standard error (s.e.) with the s.e. obtained from the inverse of the negative Hessian matrix for the log likelihood. Basic reproduction numbers (R 0w and R 0b) were calculated using the relationship R 0w = β w T and R 0b = β b T [Reference Vynnycky and White21] and the 95% CIs were based on the variance of β w, β b and T. Models were implemented using R statistical software [Reference Myung35]. Models were compared based on Akaike's Information Criterion (AIC) [Reference Burnham and Anderson36].
Vaccination coverage
Using the transmission parameter estimates, we calculated the vaccination coverage required to prevent pig-to-pig disease transmission based on the herd immunity threshold (HIT). HIT represents the minimum proportion of pigs that need to be immune in the population to limit ASFV transmission, given by HIT = 1 – (1/R 0w) [Reference Vynnycky and White21].
Sensitivity analysis
We simulated ASF outbreaks in a pig unit to assess the effects of different assumptions regarding the different latent and infectious period durations. The weaner pigs on this farm unit are assumed to be housed in a ring of six adjacent pens, with 90 pigs in each (i.e. a total of 540 pigs). We simulated 1000 independent outbreaks, initiated by the introduction of one infected pig at time t = 0. For each simulation, β w, β b and T were sampled from zero-truncated normal distributions under the assumption that these parameters are independently distributed. The mean and standard deviation for each distribution was defined so that 95% of the sampled values would fall within their respective 95% CI and L was considered to be either 3, 4 or 5 days. The likelihood surface for β w and β b suggested that the estimates for these parameters were not correlated (data not shown). We extended equation (1) to account for the fact that each pen is adjacent to two others, rather than one. Therefore, the probability that a susceptible pig in pen j becomes infected at time t is given by
$$\eqalign{p_{\,j,t} & = 1 - {\rm exp} \bigg( \! - \bigg( \displaystyle{{\beta _{\rm w} \; I_{\,j,t}} \over {N_{\,j,t}}} + \displaystyle{{\beta _{\rm b} \; I_{\,j - 1,t}} \over {( {N_{\,j - 1,t} + N_{\,j,t}} )}} \cr & \quad + \displaystyle{{\beta _{\rm b} \; I_{\,j + 1,t}} \over {( {N_{\,j + 1,t} + N_{\,j,t}} )}} \bigg) \bigg)\comma }$$
where I j,t, Ij −1,t and I j+1,t are the number of infectious animals in pens j, j – 1 and j + 1 at time t, respectively, and N j,t , Nj −1,t and N j+1,t are the total number of pigs in these pens. From these simulations, we estimated the proportion of outbreaks that fail or result in limited numbers of infected animals, the number of newly infected pigs and the total number of infectious pigs per day. Simulations were implemented with R statistical software [Reference Myung35].
RESULTS
Transmission parameters for ASFV
Assuming blood as a biological marker for infectiousness, the increase in the latent period duration from model 1 to model 3 (i.e. from L = 3 to L = 5 days) resulted in an increase in β w and β b values (Table 4). The highest β w and β b estimates were obtained for model 3 with 0·9 (95% CI 0·4–1·3) and 0·4 (95% CI 0·2–0·7) newly infected pigs per day within and between pens, respectively. The lowest AIC was obtained for model 2, followed by model 3, and model 1 had the highest AIC. A similar pattern was observed for the R 0w and R 0b estimates rather than for the β w and β b estimates. In addition, the increase in the infectious period duration (i.e. from T = 3–6 to T = 3–14 days) also resulted in an increase in R 0w and R 0b values (Table 4). The highest R 0w and R 0b estimates were generated by model 3 showing that infectious pigs would cause up to 7·2 (95% CI 2·1–14·2) and 3·5 (95% CI 1·2–7·0) new infected cases within and between pens, respectively. R 0b estimates were generally lower than R 0w estimates, but not significantly greater than the threshold value of 1 in all models (except for model 3 with the maximum infectious period) considering the lower bound of the 95% CI. Similarly, the HITs reflected the R 0 values, with the highest levels of vaccination required for model 3, i.e. from 74% (95% CI 47–84) to 86% (95% CI 55–93) depending on the infectious period duration.
Maximum likelihood estimates (95% confidence intervals) for experimental pig-to-pig transmission parameters for Georgia 2007/1 African swine fever virus strain
HIT, Herd immunity threshold; AIC, Akaike's Information Criterion.
Simulation of ASFV outbreaks
The decrease in the latent period from model 3 to model 1 (i.e. from L = 5 to L = 3 days) resulted in an increase in the probability that an outbreak does not occur after ASFV is introduced into the farm and in the probability that this outbreak does not lead to the infection of all animals (Table 5). This effect was particularly marked for a shorter infectious period, i.e. T = 3–6 days. Decreasing the infectious period in all models from T = 3–14 days to T = 3–6 days resulted in an increase in the probability that the outbreak goes extinct (i.e. non-occurrence and incomplete population infection). This effect was particularly strong for a shorter latent period of L = 3 days (model 1). As a result, model 3 produced a 7% probability of outbreak failure and 10–17% probability of not all pigs on the farm becoming infected.
Description of outbreaks simulated in a pig unit for Georgia 2007/1 ASFV strain
ASFV, African swine fever virus.
Varying durations of latent and infectious periods resulted in similar epidemic curves (Fig. 2). For example, at the peak of the epidemic curve in all models, the number of newly infected pigs per day remained around 10 (95% CI 0–25) [Fig. 2(a1–c2)].
Median (dots) and 95% confidence intervals (shaded area) of the number of newly infected pigs (red) and of the total number of infectious pigs (blue) per day during simulated outbreaks within a farm unit with Georgia ASFV 2007/1 strain based on three different models. (a1, a2) Model 1 assumed a 3-day latent period. (b1, b2) Model 2 assumed a 4-day latent period. (c1, c2) Model 3 assumed a 5-day latent period. Infectious period duration (days) was represented as a normal distribution (mean ± standard deviation) of either 4·5 ± 0·75 days (a1, b1, c1) or 8·5 ± 2·75 days (a2, b2, c2).
The increase in the latent period duration did not seem to modify the prevalence curve although it slightly increased the outbreak duration in all models, regardless of the infectious period duration (Fig. 2). For example, in model 3, the last infectious pigs were reported around 70 (95% CI 0–90) days compared to model 1 at around 60 (95% CI 0–80) days [Fig. 2(c1 vs. a1)]. The increase in the infectious period duration generally resulted in significantly higher numbers of infectious pigs per day and longer outbreak duration. For example, in model 2, more than 100 (95% CI 0–200) infectious pigs were observed at the peak of the prevalence curve and these animals were reported until 80 (95% CI 0–100) days when using the infectious period of 3–14 days (Fig. 2b2). When decreasing the infectious period to 3–6 days, approximately 50 (95% CI 0–100) infectious pigs remained at the peak of the prevalence curve and no infectious pigs were reported after 65 (95% CI 0–85) days (Fig. 2b1).
DISCUSSION
This study provides quantitative estimates for the dynamics of pig-to-pig transmission for ASFV. Model 2, assuming a 4-day latent period and using presence of live virus in blood as a marker of infectiousness, had the smallest AIC value and was thus the model with the best fit to the data from the transmission experiments. For this model, the estimate for the transmission rate parameter (β) within a pen was 0·6 (95% CI 0·3–1·0) per day, while between pens it was 0·3 (95% CI 0·1–0·5) per day. The lower AIC for model 2 compared to models 1 and 3 suggests that 4 days might be a more appropriate latent period than 3–5 days. In the transmission experiments, this corresponds to the time when ASFV genome was first recovered in oro-nasal excretions from contact animals before they showed viraemia [Reference Guinat26]. This means that the use of qPCR assay in oro-nasal samples from contact animals appears to be an effective tool for early detection of infection, as has previously been suggested for CSFV [Reference de Carvalho25]. In addition, swabs are easier and cheaper to collect as well as requiring a less stressful sample collection procedure for the pigs than taking blood samples. However, it cannot be excluded that ASFV could be excreted through other routes or that contact pigs could be infectious after different days post-exposure that were not tested here. For example, in this study viraemia was considered closely related to efficiency of ASFV transmission, although the detection of live CSFV in oropharyngeal samples was reported to be a better indicator for infectiousness than blood samples [Reference de Carvalho25]. This was not assessed in the present study due to the limited numbers of positive virus titration assay in oral samples.
Results demonstrate that, assuming a mean infectious period of 4·5 days, infectious pigs would infect on average 2·8 (95% CI 1·3–4·8) animals within their pen and 1·4 (95% CI 0·6–2·4) animals between pens. In the transmission studies, acute disease and fatal outcomes were observed for all individual infected pigs but the dynamics of the transmission process were relatively slow (Table 2), resulting in a prolonged epidemic at the herd level. In addition, specific or apparent clinical signs have been reported infrequently [Reference Guinat26–Reference Blome28]. According to field observations, no particular clinical signs of illness or values of pig mortality were observed by animal keepers in most of the ASF cases [1]. Outcomes also showed that ASF outbreaks could go extinct within a pig unit, with 10–17% probability of outbreak failure and 18–45% probability of small-scale transmission. This therefore emphasizes the potential factors that make farmers or animal workers more likely to under-report suspicious ASF cases.
Estimates for the between-pen basic reproduction number (R 0b) indicate that the disease spread among pigs is less efficient between than within pens. This suggests that adjustments to the design of the pen layout within a pig unit, such as better physically isolating neighbouring pens from each other, are likely to reduce pig-to-pig ASFV transmission. This probably also means that ASFV transmission occurs even more slowly between different pig houses or farms than was estimated among pigs in these transmission experiments. Although better pen separation should be feasible in pig farms [Reference Madec3], this will be more challenging in traditional free-ranging or backyard farms that are commonly found in Eastern Europe and where animals are not confined within a fenced area [1]. These findings also imply that animals moved from one pen to another should be assessed for health status under conditions of strict quarantine considering the latent period observed in animals. However, this will also be unrealistic in areas where lack of interest from stakeholders in eradicating the disease or poor compliance to biosecurity regulations are common [1]. In addition, design considerations for pig pens are unlikely to have practical impacts at the regional or national level for disease control strategies.
Our estimates for R 0 are quite different from those reported in previous studies on ASFV transmission: R 0 for the moderately virulent Malta strain was estimated to be 18·0 (95% CI 6·9–46·9) [Reference Sanchez-Vizcaino37], while that for the highly virulent Russia strain was 9·8 (95% CI 3·9–15·6) [Reference Gulenkin19]. This may reflect differences in virus strain, infection and infectiousness markers, pathogenesis, diagnostic tools, experimental conditions and estimation methods. Moreover, assumptions about when pigs became infected and how long they remained infectious influenced our estimates for the pig-to-pig transmission parameters and, hence, any conclusions about ASFV transmission within a farm (i.e. outbreak size and duration) and control (i.e. critical level of vaccination coverage). A study that focused on CSFV transmission found, for example, a significant difference between weaner and slaughter pigs, explained by potential various types of contact or susceptibility to the infection [Reference Klinkenberg34]. Although the clinical course of ASFV does not seem to be age dependent [Reference Blome28], infection transmission dynamics may also be influenced by host characteristics (e.g. sex or breed). Although the effect of housing systems on ASFV transmission remains unclear, the density and the contact structures that vary between backyard and different types of commercial pig farms are also likely to result in different pig-to-pig transmission parameters. This demonstrates the importance of considering the assumptions behind the model when the aim is to use the estimated values to draw inferences about the actual transmission dynamics [Reference Lloyd-Smith38, Reference Woolhouse39]. In the current study, intramuscularly inoculated pigs were used as virus sources for exposure of susceptible pigs. Although this route of inoculation has been shown to efficiently induce infection [Reference Howey40], it is not the natural route of transmission to susceptible pigs. This is likely to impact on our R 0 estimates, particularly due to significant differences in the course of pathogenesis reported between intramuscularly inoculated pigs and pigs infected by contact [Reference Guinat26]. As a result, secondary contact infection should be used in future experiments for generating R 0 estimates that are a better representation of field conditions [Reference Bouwknegt41]. Therefore, more data from the field or from transmission experiments in relation to ASFV infection dynamics in domestic pigs are needed to be able to better understand the behaviour of the virus, especially in the context of different pig production systems.
This study provides insight into quantitative pig-to-pig transmission parameters for the currently circulating ASFV strain in Eastern European countries. The experimental conditions that were used are likely to match those observed in commercial weaner pig units with good sanitary measures and with the possibility of direct and indirect contact between pigs that are housed within several adjacent pens. Findings show a low to moderate transmissibility of ASFV between pigs and that the transmission is influenced by the contact structure between these animals. These estimated transmission parameters are essential to inform future mathematical models developed with the aim of predicting between-farm ASFV spread in Eastern European countries and beyond, which will also help policy-makers to optimize interventions that are available to impede the ASFV spread.
ACKNOWLEDGEMENTS
We are grateful to the African Swine Fever Group, the Mathematical Biology Group and the Transmission Biology Group of The Pirbright Institute for assistance.
The research leading to these results has received funding from the European Union's Seventh Framework Programme (FP7/2007–2013) under grant agreement no. 311 931 (ASFORCE). This work supported by the Biotechnology and Biological Sciences Research Council (BBSRC), the Department for Environment, Food and Rural Affairs (Defra) and the National Center for Foreign Animal and Zoonotic Disease Defense (FAZD).
DECLARATION OF INTEREST
None.
Open access
