# The potential spread of highly pathogenic avian influenza virus via dynamic contacts between poultry premises in Great Britain

- Jennifer E Dent
^{1, 5}Email author, - Istvan Z Kiss
^{3}, - Rowland R Kao
^{4}and - Mark Arnold
^{2}

**7**:59

**DOI: **10.1186/1746-6148-7-59

© Dent et al; licensee BioMed Central Ltd. 2011

**Received: **4 November 2010

**Accepted: **13 October 2011

**Published: **13 October 2011

## Abstract

### Background

Highly pathogenic avian influenza (HPAI) viruses have had devastating effects on poultry industries worldwide, and there is concern about the potential for HPAI outbreaks in the poultry industry in Great Britain (GB). Critical to the potential for HPAI to spread between poultry premises are the connections made between farms by movements related to human activity. Movement records of catching teams and slaughterhouse vehicles were obtained from a large catching company, and these data were used in a simulation model of HPAI spread between farms serviced by the catching company, and surrounding (geographic) areas. The spread of HPAI through real-time movements was modelled, with the addition of spread via company personnel and local transmission.

### Results

The model predicted that although large outbreaks are rare, they may occur, with long distances between infected premises. Final outbreak size was most sensitive to the probability of spread via slaughterhouse-linked movements whereas the probability of onward spread beyond an index premises was most sensitive to the frequency of company personnel movements.

### Conclusions

Results obtained from this study show that, whilst there is the possibility that HPAI virus will jump from one cluster of farms to another, movements made by catching teams connected fewer poultry premises in an outbreak situation than slaughterhouses and company personnel. The potential connection of a large number of infected farms, however, highlights the importance of retaining up-to-date data on poultry premises so that control measures can be effectively prioritised in an outbreak situation.

## Background

For a wide range of epidemic infections, contact structures can be used to describe the potential transmission of infection in a population [1–4]. The validity of such models, however, depends on the parameterisation of the contact structures analysed. The existence of the animal movement licensing scheme and cattle tracing system in Great Britain (GB), for which the identification and movements of cattle, sheep, goats, pigs, deer and horses, must be recorded [5], allows for the reconstruction and analysis of the network of contacts in order to predict the spread of infectious disease across these industries. This is not the case for the poultry industry where, before 2005, there was no national register of poultry farms. Motivated by numerous outbreaks of H5N1 highly pathogenic avian influenza (HPAI) across the world and the occurrence of several incursions of avian influenza viruses (AIV) in GB [6, 7], information was collected by the British government on poultry farm locations, and on the frequency and types of movements between farms.

In GB, the poultry industry can be divided into the primary breeding sector and the production sector. The biosecurity levels in the primary breeding sector are considered to be consistently high, making the probability of introduction of pathogens into this sector extremely low. In the production sector, birds are purchased from a primary breeding company when they are one day old. Birds then remain on specialist rearing farms until approximately eighteen weeks of age before they are moved to production farms or to hatcheries. Before meat birds enter the food chain, a catching company may be brought in to assist in the catching of birds to be sent to slaughter. Some catching companies may operate on multiple independently owned farms, and some farms may not use a catching company at all, choosing to send birds directly to the slaughterhouse. Vehicles used to transport birds between farms and slaughterhouses are often owned by the slaughterhouse and therefore may act as a link between different production farms. Partly due to the increase in the number and types of movements made on to and off production farms and partly due to increased exposure of birds to the environment in the production sector, it is here where diseases such as AIV have the opportunity to enter a farm, rendering the production sector the focus of this study.

It has already been shown that AIVs have the potential to be spread to a large number of poultry premises via movement of humans and fomites [8, 9]. Company personnel, feed lorries, egg collectors, slaughterhouse and catching company personnel (and equipment) have been identified as contact mechanisms between farms, over which disease may transmit [10]. It has also been shown [10] that up to 97%, 42% and 11% of premises associated with slaughterhouses, catching companies and multi-site companies, respectively, are connected. However, this result was reached on the assumption that all premises using the same slaughterhouse, catching company or belonging to the same multi-site company are all potentially connected, with all links being undirected. While most models assume that all potentially infectious connections are always "active" (though see also [1]), in practice there are other factors that will limit the dissemination of disease across the commercial poultry industry. First, over the time that a premises might be expected to be able to transfer infection to other premises, the number of actual connections will depend on the frequency at which the contacts are made. Furthermore, there are likely to be important distance constraints on how far people, vehicles and livestock will travel between premises. Therefore the range over which infection is likely to travel via these means will be limited (although there is currently no maximum journey time for poultry [11]). Also, for catching companies and company personnel, it is possible that there are regional divisions within the company, e.g. geographical sub-divisions within multi-site companies and area-based teams for catching companies, and this has not previously been considered in models of HPAI transmission in GB. Therefore, obtaining detailed data on catching company and slaughterhouse movements for simulation modelling of HPAI in GB has been highlighted as necessary for realistic modelling of HPAI transmission [10]. The objective of this study was to collect detailed movement data from a large catching company and to further explore the potential for HPAI to spread via the routes highlighted.

Based on these new data, an individual farm-based transmission model was developed where nodes are poultry premises with links representing potential transmission routes between premises. Although the static approach (assuming all links between farms are potentially active) that we have previously adopted [10] was appropriate in the absence of detailed link data, the nature of these newly collated data enables us to use a dynamic network model, as is standard practice in disease modelling [2, 12], to determine if the frequency and size (distance travelled) of movements between farms is of such a nature as to reduce the concern for a large HPAI epidemic in GB. In this updated model, the presence of links between premises is drawn from the collected movement data. This model also incorporates link directionality allowing for a more realistic and accurate model. Here, the potential for an HPAI epidemic is determined by considering the results of the individual-based transmission model.

## Methods

### Data sources

Population data on all commercial poultry premises housing 50 or more birds, and within 15 km of farms visited by teams belonging to the catching company, were taken in November 2007 from an extract of the GB Poultry Register (GBPR), provided by the Department for environment, food and rural affairs (Defra). The addition of these premises allowed us to consider disease transmission that can occur between premises that are in close spatial proximity, enabling us to consider how likely infection is to jump from the network of premises for which we have movement data to another network of premises for which no movement data were available for analysis. The addition of these premises brought the number of premises studied up to 10,692. The total number of premises recorded in the GBPR extract studied was 24,075. Although the total number of premises serviced by the catching company is small (approximately 2% of all premises in the GBPR), by comparing the number of premises serviced by the catching company with the number of premises known to be associated with the largest catching companies in England and Wales, we estimate that between 30% and 50% of premises serviced by a major contracted catching company are accounted for in this dataset [13 and unpublished data]. No data were available for small, independent catching companies.

Premises in the catching company database were matched to premises in the GBPR so that the model could be parameterised using location, premises population size and species data. Temporal aspects were accounted for by preserving the order in which poultry premises were visited by individual catching teams on a daily basis. Additional information on which premises belong to integrated companies was obtained from a sample survey of integrated companies.

*pers. comm*.] to inform the model of the likely frequency of movements of personnel between farms, based on species type, farm size and distance between farms belonging to the same company. The probability of companies using staff on multiple farms was estimated as well as the estimated frequency of visits of veterinary officials and area managers to poultry premises (Table 1).

Model parameters.

Parameter | Value | Source |
---|---|---|

Incubation period | Up to 1 day (shedding after 8 hours, death after one day) | |

Survival of virus on seed premises | Up to15 days | [16] |

Probability of staff working on multiple premises (Farm size) | 0.45 (< 50,000 birds) 0.1 (50,000 to 200,000 birds) 0 (> 200,000 birds) | P. McMullin ( |

Distance travelled by staff between premises | Up to 35 km | P McMullin ( |

Frequency of vet visits | Every 50 days | P. McMullin ( |

Frequency of manager visits | Every 10 days (non-layer farms) Every 50 days (layer farms) | P. McMullin ( |

Probability catching team catches species (multi species farms only) | 0.7 Chicken 0.12 Turkey 0.16 Duck/goose 0.01 Other | Calculated from catching company data, where species type available |

Probability of catching team, slaughterhouse and owner transmission | 0 to 0.2, in steps of 0.01, with additional parameter at 0.001 added. | N/A |

Time to detection | 2 to 6 days (15 days later for ducks/geese) |

### Descriptive analysis

A descriptive analysis of the collected data is given in Additional File 1. This analysis was used to identify any trends in the data that could have important implications for disease transmission. In particular we considered how far catching teams and slaughterhouse vehicles travelled between premises, as it was shown in [10] that this had a large impact on the potential size of an HPAI outbreak. The frequency of catching team movements to premises, dependent on farm size, was also considered in order to determine if premises size should be recorded as an output of the simulation model.

### Simulation model (see also Additional File 2)

A stochastic simulation model at the farm level was developed where farms were classed as susceptible, infected, detected or culled. HPAI could be transmitted between premises in close spatial proximity or through contact via catching teams, slaughterhouse vehicles or personnel movements within an integrated (multi-site) company. A random number generator chooses a premises in the network to infect and a random date of infection (within the 886 days covered by the catching company data set). If the seed farm was not visited within 15 days of this time point, assumed to be the maximum time that HPAI would survive in the farm environment, then transmission will be limited to local spread, often resulting in little or no onward transmission from the seed premises. When transmission beyond the seed premises did occur, outbreaks were allowed to run their course.

We assumed a time step of one day in the model, so that for each day of the simulation, once a premises had become infected, we assumed silent spread up to the time of detection. Detection and culling dates were set within the model at time of infection and were dependent on whether the infected premises was in a protection zone (PZ), a surveillance zone (SZ), or neither, as described below.

On detection of notifiable HPAI in poultry, 10 km SZs and 3 km PZs are typically set up around infected premises. In the model, we assumed no transmission within the PZ/SZ via the normal movement of catching companies or slaughterhouse equipment since all movements in those zones would be monitored. Therefore, spread would only continue within the PZ/SZ via local spread. Time to detection within an infected flock was assumed to occur between 2 and 6 days (mean at 4 days) after infection [14–18], with the mean reduced to 3 days for flocks within a PZ or SZ. Culling occurs at the end of an epidemic day. Using the time taken to cull birds in the most recent outbreak of HPAI H5N1 in GB [19], we assumed that culling was completed on infected farms within 3 days of detection. One hundred simulations were run for a range of parameter values (Table 1). Up to 1,000 simulations were run for a subset of parameter combinations with no qualitative change in the results. For each simulation, the times at which each farm had become infected and the times at which their state changed from infected to detected and from detected to culled was recorded.

#### Catching team and slaughterhouse movements

When a movement occurred between infected and susceptible farms, infection was spread between farms with probability relating to the type of movement made. Where multiple species were held on one farm, we assumed that catching teams catch, on average, one species per visit with probability defined in Table 1. This species became infected with probability of transmission via catching team. If a susceptible farm had multiple links from an infected premises, then each link was treated independently, and the probability of infection therefore increased. We assumed that catching team and slaughterhouse vehicles do not remain infectious overnight, as we assume that effective decontamination procedures were in place. A premises was designated as infected if one or more species on the premises was infected.

#### Company (personnel) movements

We assume that spread of infection between farms belonging to the same multi-site, integrated company could occur either via the movements of area managers or veterinary officials between premises, or via staff working on multiple farms. Movements of veterinary officials, of area managers and of company personnel were simulated on a per day basis by using farm size and distance between farms belonging to the same multi-site company to identify using parameters in Table 1 if a) two premises are to be visited on the same day by the same person and, b) if the link made between premises will result in transmission of disease. The model assumes independence between days so that the last time a premises was visited is not accounted for.

#### Local (spatial) spread (see Additional Files 2 and 3)

Based on expert opinion, we assume spatial (primarily airborne) spread in GB is likely to occur with small probability (p ≤ 0.01) and only for distances up to a maximum 0.5 km [D. Alexander, R. Irvine *pers. comm*.], according to the density kernel that is given in Equation 1.

*For dist < 0.5 km:*

*Else p*(*trans|dist*) = 0

A sensitivity analysis of the model to the assumption that airborne spread could occur is given in Additional File 3.

### Analysis of simulation output

We first consider the proportion of outbreaks that resulted in onward spread beyond the seed case. Here, the results follow a linear trend and, as the outcome is a binary variable (essentially secondary spread, or no secondary spread) dependent on explanatory variables that can be categorised into multiple levels, the analysis lends itself to a logistic regression. Thus a binary logistic regression was done (using Minitab v16) on the proportion of outbreaks that resulted in onward spread beyond the seed case. We next consider final epidemic size. In order to determine how the different types of transmission affect the epidemic size, two logistic regression models were fitted (Minitab v16). In the first, the binary response variable describes whether a small (< 25 premises) epidemic occurs or not. In the second, the binary response variable describes whether a large (> 65 premises) epidemic occurs or not. In both cases, the explanatory variables are the simulated transmission probabilities for AIV transmission via catching company, slaughterhouse- and owner-related movements.

## Results

- 1.
When the frequency of movements is not accounted for, slaughterhouse related movements connected 94% of premises, catching team movements connected 76% of premises and owner movements 11% of premises that are associated with the catching company. However, when time is considered, catching teams connect only 2 premises per day and slaughterhouses an average of 3 premises per day.

- 2.
Contrary to expectations, the data presented show that premises do use multiple slaughterhouses and are associated with multiple catching teams within the same company (consecutively over the time period studied). There is no overlap between poultry companies (i.e. poultry premises are associated with only one poultry company).

- 3.
There is an increase in frequency of visits to larger premises, implying that these premises will be at a higher risk of infection should infection be transmitted by catching team or slaughterhouse vehicles/equipment.

- 4.
Slaughterhouse vehicles and catching teams travel long distances between premises, with 72% of movements between premises exceeding 10 km in length.

### Simulation modelling

One hundred simulations were run for a total of 10,648 scenarios. Each scenario represents a different combination of transmission values for each of the transmission routes studied. In particular, each scenario was created by ranging parameters from 0 to 0.2, in a step-wise fashion, such that each parameter took on one of 22 possible values within this range, giving rise to *22*^{
3
} *= 10,648* scenarios.

### Proportion of positive epidemics spread beyond the index case

The aims of the simulation model are to determine if a large outbreak of AIV is possible in the poultry industry in GB and, if so, what might cause a large outbreak to occur. One way of answering the first question is to consider how often infection spreads beyond the seed premises. That is to ask "how many simulated outbreaks result in secondary spread?"

for *p*_{
j
} probability of infection for via link type *j*.

Infection resulted in secondary spread (beyond the seed premises) in up to 35% of scenarios. The simulation that gave the maximum number of cases that spread beyond the seed premises was from the following scenario: catching team (cc) = 0.04, company personnel (owner) = 0.19 and slaughterhouse (sh) = 0.13. This suggests that high probabilities of transmission are not necessary in all three potential transmission routes for infection to (relatively) frequently spread beyond the index case.

*p = 0.001*. Movements related to slaughterhouse vehicles are significant in the probability that an outbreak will result in onward spread only when the probability of transmission is high enough (here, the model predicts a rate of > 0.06 for a significant effect to be seen). Further, the odds ratios also tell us that, under the assumptions made, as the rate of transmission increases for owner movements in particular, the effect on the probability of secondary spread is increasingly large, with the odds ratio rising to 2.11 (2.03, 2.19 (95% CIs)) for a transmission rate of 0.14 compared to zero. This suggests that the probability of secondary spread beyond the seed premises is not uniformly affected by transmission rates across the different link types. This is driven by the characteristics of the networks over which disease can spread.

Binary logistic regression, with odds ratios calculated for the probability of secondary spread versus catching company transmission rates

Transmission rate | Odds Ratio | Lower 95% CI | Upper 95% CI | p-value |
---|---|---|---|---|

0.001 | 0.97 | 0.93 | 1 | 0.064 |

0.01 | 1 | 0.97 | 1.04 | 0.813 |

0.02 | 0.95 | 0.92 | 0.98 | 0.005 |

0.03 | 0.97 | 0.94 | 1.01 | 0.14 |

0.04 | 1.02 | 0.98 | 1.06 | 0.303 |

0.05 | 0.96 | 0.93 | 0.99 | 0.024 |

0.06 | 0.96 | 0.93 | 1 | 0.043 |

0.07 | 0.99 | 0.95 | 1.02 | 0.451 |

0.08 | 0.96 | 0.93 | 1 | 0.038 |

0.09 | 0.95 | 0.92 | 0.99 | 0.011 |

0.1 | 0.97 | 0.94 | 1.01 | 0.124 |

0.11 | 0.98 | 0.94 | 1.01 | 0.23 |

0.12 | 0.99 | 0.96 | 1.03 | 0.677 |

0.13 | 1 | 0.97 | 1.04 | 0.906 |

0.14 | 0.98 | 0.94 | 1.01 | 0.187 |

0.15 | 0.98 | 0.95 | 1.02 | 0.267 |

0.16 | 0.96 | 0.92 | 0.99 | 0.018 |

0.17 | 1 | 0.96 | 1.03 | 0.871 |

0.18 | 0.98 | 0.95 | 1.02 | 0.313 |

0.19 | 0.96 | 0.93 | 1 | 0.047 |

0.2 | 0.98 | 0.95 | 1.02 | 0.359 |

Binary logistic regression, with odds ratios calculated for the probability of secondary spread versus owner transmission rates

Transmission rate | Odds Ratio | Lower 95% CI | Upper 95% CI | p-value |
---|---|---|---|---|

0.001 | 0.98 | 0.94 | 1.03 | 0.488 |

0.01 | 1.09 | 1.04 | 1.14 | 0 |

0.02 | 1.19 | 1.14 | 1.24 | 0 |

0.03 | 1.27 | 1.22 | 1.33 | 0 |

0.04 | 1.38 | 1.33 | 1.44 | 0 |

0.05 | 1.42 | 1.37 | 1.48 | 0 |

0.06 | 1.47 | 1.41 | 1.53 | 0 |

0.07 | 1.53 | 1.47 | 1.59 | 0 |

0.08 | 1.66 | 1.59 | 1.72 | 0 |

0.09 | 1.74 | 1.67 | 1.81 | 0 |

0.1 | 1.75 | 1.69 | 1.83 | 0 |

0.11 | 1.86 | 1.79 | 1.94 | 0 |

0.12 | 1.97 | 1.89 | 2.04 | 0 |

0.13 | 1.94 | 1.87 | 2.02 | 0 |

0.14 | 2.11 | 2.03 | 2.19 | 0 |

0.15 | 2.09 | 2.01 | 2.17 | 0 |

0.16 | 2.19 | 2.11 | 2.27 | 0 |

0.17 | 2.24 | 2.16 | 2.33 | 0 |

0.18 | 2.33 | 2.24 | 2.42 | 0 |

0.19 | 2.38 | 2.29 | 2.47 | 0 |

0.2 | 2.38 | 2.29 | 2.47 | 0 |

Binary logistic regression, with odds ratios calculated for the probability of secondary spread versus slaughterhouse transmission rates

Transmission rate | Odds Ratio | Lower 95% CI | Upper 95% CI | p-value |
---|---|---|---|---|

0.001 | 0.97 | 0.93 | 1.01 | 0.093 |

0.01 | 0.99 | 0.95 | 1.02 | 0.435 |

0.02 | 0.99 | 0.95 | 1.03 | 0.597 |

0.03 | 1 | 0.96 | 1.03 | 0.861 |

0.04 | 1 | 0.96 | 1.03 | 0.824 |

0.05 | 1 | 0.96 | 1.04 | 0.927 |

0.06 | 1.04 | 1 | 1.08 | 0.036 |

0.07 | 1.03 | 1 | 1.07 | 0.068 |

0.08 | 1.04 | 1.01 | 1.08 | 0.019 |

0.09 | 1.04 | 1.01 | 1.08 | 0.019 |

0.1 | 1.05 | 1.01 | 1.09 | 0.012 |

0.11 | 1.07 | 1.03 | 1.11 | 0 |

0.12 | 1.05 | 1.01 | 1.09 | 0.008 |

0.13 | 1.09 | 1.05 | 1.13 | 0 |

0.14 | 1.08 | 1.04 | 1.12 | 0 |

0.15 | 1.09 | 1.06 | 1.13 | 0 |

0.16 | 1.08 | 1.04 | 1.12 | 0 |

0.17 | 1.08 | 1.04 | 1.12 | 0 |

0.18 | 1.09 | 1.06 | 1.13 | 0 |

0.19 | 1.1 | 1.06 | 1.14 | 0 |

0.2 | 1.1 | 1.06 | 1.14 | 0 |

Binary Logistic regression: secondary spread versus transmission rates for interaction between transmission routes at different levels of transmission.

Predictor | Coefficient | SE | Odds Ratio | Lower 95% CI | Upper 95% CI | p-value |
---|---|---|---|---|---|---|

Constant | -2.16587 | 0.071836 | 0 | |||

owncat | ||||||

1 | 0.004215 | 0.076562 | 1 | 0.86 | 1.17 | 0.956 |

2 | 0.415397 | 0.075335 | 1.51 | 1.31 | 1.76 | 0 |

3 | 0.631829 | 0.074869 | 1.88 | 1.62 | 2.18 | 0 |

shcat | ||||||

1 | -0.22299 | 0.077486 | 0.8 | 0.69 | 0.93 | 0.004 |

2 | -0.12265 | 0.077053 | 0.88 | 0.76 | 1.03 | 0.111 |

3 | -0.0784 | 0.076875 | 0.92 | 0.8 | 1.07 | 0.308 |

owncat*shcat | ||||||

1*1 | 0.248773 | 0.082484 | 1.28 | 1.09 | 1.51 | 0.003 |

1*2 | 0.227419 | 0.082041 | 1.26 | 1.07 | 1.47 | 0.006 |

1*3 | 0.229411 | 0.081855 | 1.26 | 1.07 | 1.48 | 0.005 |

2*1 | 0.207304 | 0.0812 | 1.23 | 1.05 | 1.44 | 0.011 |

2*2 | 0.154659 | 0.080772 | 1.17 | 1 | 1.37 | 0.056 |

2*3 | 0.140458 | 0.080593 | 1.15 | 0.98 | 1.35 | 0.081 |

3*1 | 0.228509 | 0.0807 | 1.26 | 1.07 | 1.47 | 0.005 |

3*2 | 0.174293 | 0.080272 | 1.19 | 1.02 | 1.39 | 0.03 |

3*3 | 0.15525 | 0.080095 | 1.17 | 1 | 1.37 | 0.053 |

### Epidemic size

In this part of the analysis, only epidemics that result in spread beyond the seed premises are considered. This accounts for approximately 15% of all simulation results.

Effect of seed premises on outbreak size.

Seed premises size | Number unique premises in small epidemics (seed) | Number unique premises in large epidemics (seed) | Proportion outbreaks resulting in large epidemics |
---|---|---|---|

Small (≤ 100,000 birds) | 35 | 20 | 0.57 |

Medium (100,000 - 200,000 birds) | 59 | 17 | 0.29 |

Large (> 200,000 birds) | 141 | 37 | 0.26 |

In order to determine how the different types of transmission affect the epidemic size, two logistic regression models were fitted. In the first, the binary response variable describes whether a small (< 25 premises) epidemic occurs or not. In the second, the binary response variable describes whether a large (> 65 premises) epidemic occurs or not. In both cases, the explanatory variables are the simulated transmission probabilities for AIV transmission via catching company, slaughterhouse- and owner-related movements. The results are shown in Additional File 4 Table S2 to S7.

For small epidemics, Additional File 4 Tables S2 to S4 show that catching company movements have a significant influence on the results for a range of probability values between 0.02 and 0.16. Interestingly, when these results are significant (the odds ratio confidence intervals do not contain zero), the odds ratios show that the probability of a small epidemic decreases (the odds ratios are less than 1) with an increase in catching company transmission rates, when compared to zero. This suggests that an increase in catching company transmission might result in a higher proportion of epidemics being larger. Additional File 4 Table S3 shows that transmission via owner movements is significant at all levels above *p = 0.001*. Above this value, the odds ratios are all larger than one, implying that increasing the rate of transmission results in the likelihood of a small epidemic occurring to increase. For slaughterhouses (Additional File 4 Table S4), significant results are obtained for transmission rates > 0.05. The strength of the significance does not increase in proportion with the increase in transmission, with all transmission rates > 0.12 having an odds ratio value of between 1.07 and 1.10. These results therefore suggest that the most influential parameter for the probability of a small epidemic to occur is transmission via owner movements.

Additional File 4 Table S5 and S6 show that, contrary to expectations, neither catching company nor owner movements play a significant role in the probability that an outbreak will result in a large epidemic. For large epidemics, the most influential predictor is the transmission rate via slaughterhouse linked movemernts (Additional File 4 Table S7). Analysis of the odds ratios for slaughterhouse transmission versus large epidemics shows that this transmission route is only influential if it is high enough (above 0.12). However, when it is high enough, the upper 95% limits (for the odds ratios) show that an increase from zero transmission to a higher transmission rate will result in a large epidemic being up to 28 times more likely. This is a very striking result with heavy implications on resources, for example, in the event of an outbreak. It is therefore essential to determine the true probability of transmission via this route.

### Spatial spread

While the majority of outbreaks did not result in further onward transmission from the index premises, outbreaks could potentially cover up to 20% of the population (for a range of parameter values < 0.2) for which network data were available, covering distances of up to 730 km between premises (see Additional File 5, Figure S1). Although for the largest distances to be covered (> 700 km), at least one transmission parameter must be as high as 0.12, occasionally, distances over 600 km between premises are reached for transmission parameters between zero and 0.07. This has important implications for the availability of control resources as the number of premises in the SZs will be greater if dissemination of virus is geographically widespread and therefore potentially involving a larger number of local disease control centres. Large epidemics invariably resulted in widespread geographical dissemination of virus.

In this model, infection can only be spread into premises that are not serviced by the catching company, by spatial transmission of disease, to premises within 500 m of infected premises. This results in infection of premises that are potentially connected to different sub-networks (via other catching companies, slaughterhouses or poultry companies for example) in less than 1% of the simulations run. However, we have seen that transmission via sh-linked movements is an important factor in determining final epidemic size, and slaughterhouses that are included in the network studied may also be used by poultry premises not included here. This implies that if this route is important, infection may leak into other sub-networks of the industry much more frequently.

## Discussion

Despite the extent of data previously available on the British poultry industry, the detailed contact structures within the poultry industry in GB have only been poorly understood. Previous studies have been able to identify potential contact structures but assumptions have had to be made on the frequency and patterns of movements between farms [10, 15, 20]. Whilst it is important to acknowledge that the models presented here rely heavily on expert opinion (which is arguably a drawback of such a modelling approach), in the absence of outbreak data for AIV in GB, this cannot be avoided. For this reason, we have considered many scenarios by varying parameter values and by combining expert opinion with real time movement data from a large catching company, we have been able to adopt a similar approach to that used in [1] to investigate the potential spread of AIV in the GB poultry industry.

The results presented here show that restrictions on the frequency of movements can have an important role in determining disease spread risk. In particular, connections via slaughterhouses can connect a large number of premises over a large geographical area, important in the potential for virus dissemination. Spread via slaughterhouse-linked movements is most prominent when partial flock depopulation is being undertaken at a farm, as this action results in more premises being visited in one day and potential infection of birds that remain on the farm. This is also an important output for the control of diseases other than HPAI, such as *Salmonella* or Campylobacter spp., where the slaughterhouse is a more likely reservoir for pathogens [21]. We note here that, whilst slaughterhouses and catching teams are separated in this study, in some cases one might group the two transmission routes together under the assumption that any movement that arises due to a catching team visiting a farm is considered a 'catching company' movement. The results are likely to be sensitive to such an assumption and thus it is important not to misinterpret them. However, the principles used in this study remain valid for the potential transmission of diseases spread by the faeco-oral route, such as Campylobacter spp. and *Salmonella*, as well as different strains of HPAI. The model is well-suited to investigating diseases where expert opinion does not have to be so heavily relied upon for model parameterisation (as expert opinion adds uncertainty to results, resulting in one only being able to answer 'what-if' scenarios, in a situation the model assumptions may affect interpretation of results). In Campylobacter research, for example, one would expect the model results to be quite different due to difference in epidemiological characteristics of Campylobacter spp. compared to HPAI. With a perceived higher prevalence of the pathogen, we would expect the results of the model presented here, applied to Campylobacter spp. to show that catching company movements are likely to have a bigger effect on the spread of disease between farms.

Despite the relatively heavy use of expert opinion to estimate model parameters in this study - in particular for the frequency of movements made by company personnel, we can use the model presented here to hypothesise about the importance of different types of potentially infectious links between poultry premises and we can conclude from these results that, where slaughterhouses can act as a reservoir for pathogens, the spread via this route should be minimized. This can be achieved through additional bio-security measures, such as thorough cleaning of the crates and vehicles that carry the birds, for example.

The results that catching team movements have little effect both on the probability of an outbreak resulting in onward spread beyond the seed premises and on the probability of a large epidemic occurring are important results, as they suggest that the number of farms that a catching team visits during the infectious period of the virus is too low to link a high number of farms, in GB, during an epidemic. For pathogens that can survive for longer periods in the environment or that are more prevalent than HPAI (such as Campylobacter spp.), the number of farms that can be linked by catching team movements will be (potentially significantly) higher. However, while extensive and therefore of value, the data used here correspond to only one (large) catching company that is made up of a 68 distinct catching teams. As each farm may be visited by one or more of the catching teams, there are no distinct regional divisions apparent within this company as was initially expected. Further, these data do not consider further spread once other networks (e.g. connected by slaughterhouses and catching companies) contain infected premises.

Although all three transmission routes were positive when a large proportion of (simulated) outbreaks resulted in spread beyond the seed premises, the fitting of a regression models suggests that only company personnel movements significantly influence the probability that infection will spread beyond the seed premises. This highlights the importance of obtaining more accurate estimates on the frequency of movements of company personnel and the probability of transmission via this route.

There was a significant interaction effect for the owner*slaughterhouse interaction on the proportion of outbreaks that result in onward spread. However, the combinations of potential transmission of disease via catching company and company personnel movements, or slaughterhouse-linked and catching company movements have little effect on the proportion of outbreaks that result in onward spread, particularly compared to the individual owner effect. This can be explained by the frequency of movements relative to premises size (Additional File 1 Figures S6 to S9), such that the increased frequency of catching company movements in particular (and also, but less so for slaughterhouse-linked movements), to larger premises is not high enough to force these potential transmission routes to have a large effect on the proportion of outbreaks that result in spread beyond the seed premises, compared to transmission via owner movements. Having highlighted owner movements as important in previous studies [10] and given that they can have a large effect on the number of outbreaks resulting in an epidemic, it is recommended that data collection is expanded to include movement data from an integrated company, furthering our ability to provide more robust estimates of epidemic size and likelihood.

The results show that there is a "jump" from epidemics of size lower than 23 infected farms (< 5% of premises), to epidemics containing more that 65 infected farms (~20% of premises). This is in line with results published by [22], who report that a predictor of the need to intensify control efforts in GB is whether an outbreak exceeds 20 infected premises. The results follow the pattern of epidemic outbreak sizes (at least qualitatively) as expected for any stochastic epidemic model, with epidemics either going extinct early, or growing to reach a substantial proportion of the population. Whilst this result, which represents a threshold for the basic reproduction number, *R*_{
0
} , will be affected by the structure of the networks, investigating network structure alone is not enough to fully investigate the effect of *R*_{
0
} . To do this, one would need to understand the effect of the individual transmission rates on the probability of a large outbreak.

When comparing the results for small epidemics against those for large epidemics, two factors that differ significantly between the two categories are worth noting: the effect of the probability of transmission via slaughterhouse movements and seed premises size. Large epidemics are up to 28 times more likely for higher levels of slaughterhouse transmission (compared to zero), implying that the characteristics of the network of slaughterhouse links are maintained even when a time component and control measures are added, resulting in connectivity between a higher proportion of premises via this route than via any other route. This result confirms that slaughterhouses are an important factor in this model. The size of seed premises plays a role here as there is an increase in frequency of catching team and slaughterhouse visits to larger premises (Additional File 1 Figure S10). This results in large outbreaks being more likely to occur, as a result of infection in a large seed premises. It is reiterated however that this does not imply that infection seeded in large premises will always result in a large outbreak. Nevertheless, this result does suggests that if premises are to be prioritised during contact tracing, there will be some benefit to targeting large premises ahead of smaller ones in a epidemic situation. Further investigation into all premises included in these epidemics to identify whether the same premises are included in the large epidemics is highlighted here as an area for further research. This will also identify premises that might be considered particularly high risk.

We note that all slaughterhouses that appear in the movement data analysed are recorded as slaughtering birds from farms that are not visited by the catching company studied. This implies that the network of premises studied is not closed; with up to 131 additional farms sending birds to the same slaughterhouse (unpublished data), the possibility of disease spreading into other sub-networks within the industry is potentially high. It is therefore very important to ensure the data held on slaughterhouses and their customers is both complete and up to date. This will enable better prioritisation of the potentially large number of premises that could undergo surveillance in an outbreak situation.

Our results show that the distribution of poultry premises in GB is not dense enough for airborne transmission of AIV contribute significantly to between premises spread amongst premises recorded in the GBPR, so long as the distance for airborne transmission is less than 500 m. This has not been the case in past outbreaks in other countries, such as the Netherlands and Italy, where local spread is likely to have played a role in the transmission of disease from one farm to another. Should a virus strain that can easily transmit via airborne transmission be modelled, then local spread may result in spread between premises that have no other direct connections. For other virus strains, this could have a large impact on the proportion of outbreaks resulting in spread beyond the seed premises and the maximum epidemic size. This implies that there is possible scope to reduce the size of the 10 km SZs, freeing resources for use elsewhere. This could be explored further by using network data currently available to explore how large a SZ should be, taking into account resource constraints and simulating over a range of assumptions regarding transmission rates. The mean number of premises affected by an epidemic may be dependent not only on the underlying epidemiological parameters, but also on the total resources available. Resource constraints were not included in this model but the model could be adapted to aid future work in this area, important for exploring optimal resource allocation in order to provide the most efficient detection of AIV and the curtailing of the outbreak.

## Conclusions

Previous work has shown that large proportions of the poultry industry are potentially connected by catching companies and by slaughterhouse [10]. However, such analyses did not take into account the restriction in the number of interaction events that could occur over the course of a typical infectious period. Including these effects, such as via the explicit spatio-temporal simulations explored in this study, shows that such restrictions can have an important role in determining disease spread risk.

In line with previous work [10], we have shown that slaughterhouses connect the highest number of premises in the poultry industry. Furthermore, the potentially high frequency of company personnel between farms renders this type of movement more important in the beginning of an outbreak. Contrary to expectations, however, the frequency of movements of catching teams between premises is not high enough to connect large numbers of premises, reducing the potential for a large outbreak spread via this route, in GB. The size of seed premises played a role in final epidemic size suggesting that there will be some benefit to targeting large premises ahead of smaller ones in an epidemic situation. The ability of the virus to jump from one sector of the industry to another highlights the importance of keeping data on movements on and off poultry farms both detailed and up to date.

## Declarations

### Acknowledgements

We would like to thank Defra for funding the project. We also thank Prof George Gettinby and Dr. Louise Kelly for statistical input at the final stages. Further, we gratefully acknowledge the following for helpful comments on the manuscript: Alasdair Cook and Richard Irvine (Animal Health and Veterinary Laboratories Agency) and Victor Del Rio Vilas (Defra). Finally, we would like to thank Lucy Snow (Animal Health and Veterinary Laboratories Agency) for her help with data collection.

## Authors’ Affiliations

## References

- Green DM, Kiss IZ, Kao RR: Modelling the initial spread of foot and mouth disease through livestock movements. Proc R Soc B. 2006, 273: 2729-35. 10.1098/rspb.2006.3648.PubMed CentralPubMedView ArticleGoogle Scholar
- Kao RR, Green DM, Johnson J, Kiss IZ: Disease dynamics over very different time-scales: foot-and-mouth disease and scrapie on the network of livestock movements in the UK. JR Soc Interface. 2007, 4: 907-916. 10.1098/rsif.2007.1129.View ArticleGoogle Scholar
- Kiss I, Green D, Kao RR: The network of sheep movements within Great Britain: Network properties and their implications for infectious disease spread. JR Soc Interface. 2006, 3: 669-677. 10.1098/rsif.2006.0129.View ArticleGoogle Scholar
- Read JM, Keeling MJ: Disease evolution on networks: the role of contact structure. Proc R Soc Lond B. 2003, 270: 699-708. 10.1098/rspb.2002.2305.View ArticleGoogle Scholar
- Defra: Livestock identification and movements. 2011, consulted 06October 2011, [http://www.defra.gov.uk/food-farm/animals/movements]Google Scholar
- Gibbens N: Avian influenza outbreak in Oxfordshire. Vet Rec. 2008, 162 (24): 795.PubMedView ArticleGoogle Scholar
- Manvell RJ, Londt BZ, Ceeraz V, Cox WJ, Essen S, Banks J, Slomka MJ, Pavlidis T, Irvine RM, Wilesmith JW, Sharpe CE, Hurst A, Alexander DJ, Brown IH: Low pathogenic avian influenza in domestic fowl in Norfolk, England, March and April, 2006. Vet Rec. 2008, 162 (9): 278-80. 10.1136/vr.162.9.278.PubMedView ArticleGoogle Scholar
- Alexander DJ: The epidemiology and control of avian influenza and Newcastle disease. Comp Pathol. 1995, 112 (2): 105-126. 10.1016/S0021-9975(05)80054-4.View ArticleGoogle Scholar
- Bahl AK, Langston A, Van Deussen RA, Pomeroy BS, Newman J, Karunakaran D, Halvorson D: Prevention and control of avian influenza in turkeys. Proc Annu Meet US Anim Health Assoc. 1979, 83: 355-363.Google Scholar
- Dent JE, Kao RR, Kiss IZ, Hyder K, Arnold M: Contact structures in the poultry industry in Great Britain: Exploring transmission routes for a potential avian influenza virus epidemic. BMC Veterinary Research. 2008, 4: 27-10.1186/1746-6148-4-27. [http://www.biomedcentral.com/1746-6148/4/27]PubMed CentralPubMedView ArticleGoogle Scholar
- Defra: Welfare of Animals During Transport. consulted 11 March 2010, [http://www.defra.gov.uk/food-farm/animals/welfare/transport/]
- Keeling MJ, Eames KT: Networks and epidemic models. JR Soc Interface. 2005, 2 (4): 295-307. 10.1098/rsif.2005.0051.View ArticleGoogle Scholar
- Gittins J, Canning P: Review of the Poultry Catching Industry in England and Wales. 2006, consulted on 06 October 2011, [http://www.food.gov.uk/multimedia/pdfs/poultycatchreview.pdf]Google Scholar
- Pantin-Jackwood MJ, Kapczynski DR, Wasilenko JL, Sarmento L: Comparison of the pathogenicity of different H5N1 HPAI viruses in chickens and ducks. Proceedings of the American Association of Veterinary Laboratory Diagnosticians 50th Annual Conference. 2007, 125:Google Scholar
- Sharkey KJ, Bowers RG, Morgan KJ, Robinson SE, Christley RM: Epidemiological consequences of an incursion of highly pathogenic avian influenza into the British poultry flock. Proc R Soc B. 2007, 275: 1630-[http://rspb.royalsocietypublishing.org/content/275/1630/19.full]Google Scholar
- Savill NJ, Rose SG, Woolhouse MEJ: Detection of mortality clusters associated with highly pathogenic avian influenza in poultry: a theoretical analysis. JR Soc Interface. 2008, 5: 1049-1419.View ArticleGoogle Scholar
- Stegeman A, Bouma A, Elbers ARW, De Jong MCM, Nodelijk G, De Klerk F, Koch G, van Boven M: Avian influenza A virus (H7N7) epidemic in The Netherlands in 2003: Course of the epidemic and effectiveness of control measures. J Infect Dis. 2004, 190 (12): 20.View ArticleGoogle Scholar
- Yoon H, Park CK, Nam HM, Wee SH: Virus spread pattern within infected chicken farms using regression model: the 2003-2004 HPAI epidemic in the republic of Korea. Journal of Veterinary Medicine Series B. 2005, 52 (10): 428-431. 10.1111/j.1439-0450.2005.00891.x.PubMedView ArticleGoogle Scholar
- Defra: Preliminary epidemiology report: avian influenza outbreak in Suffolk, November 2007 as at 26 November 2007. consulted 06 October 2011, [http://archive.defra.gov.uk/foodfarm/farmanimal/diseases/atoz/ai/documents/ai-prelim-epireport071129.pdf]
- Garske T, Clarke P, Ghani A: The transmissibility of highly pathogenic avian influenza in commercial poultry in industrialised countries. PloS ONE. 2007, 2 (4): e349-10.1371/journal.pone.0000349.PubMed CentralPubMedView ArticleGoogle Scholar
- Evans SJ, Sayers AR: A longitudinal study of Campylobacter infection of broiler flocks in Great Britain. Preventive Veterinary Medicine. 2000, 46: 209-223. 10.1016/S0167-5877(00)00143-4.PubMedView ArticleGoogle Scholar
- Truscott J, Garske T, Chis-Ster I, Guitian J, Pfeiffer D, Snow L, Wilesmith J, Ferguson N, Ghani A: Control of a highly pathogenic H5N1 avian influenza outbreak in the GB poultry flock. Proceedings of the Royal Society B. 2000, 274 (1623): 2287.View ArticleGoogle Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.