# Role of Live-Duck Movement Networks in Transmission of Avian Influenza, France, 2016–2017

Claire Guinat, Benoit Durand, Timothee Vergne, Tifenn Corre, Séverine

Rautureau, Axelle Scoizec, Sophie Lebouquin-Leneveu, Jean-Luc Guérin,

Mathilde Paul

► To cite this version:

Claire Guinat, Benoit Durand, Timothee Vergne, Tifenn Corre, Séverine Rautureau, et al.. Role of Live-Duck Movement Networks in Transmission of Avian Influenza, France, 2016–2017. Emerging Infectious Diseases, 2020, 26 (3), pp.472-480. 10.3201/eid2603.190412 . hal-03207007

HAL Id: hal-03207007

https://hal.inrae.fr/hal-03207007

Submitted on 23 Apr 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Role of Live-Duck Movement Networks in Transmission of Avian Influenza, France, 2016–2017

Claire Guinat, Benoit Durand, Timothee Vergne, Tifenn Corre, Séverine Rautureau, Axelle Scoizec, Sophie Lebouquin-Leneveu, Jean-Luc Guérin, Mathilde C. Paul

The relative roles that movement and proximity networks play in the spread of highly pathogenic avian influenza (HPAI) viruses are often unknown during an epidemic, preventing effective control. We used network analysis to explore the devastating epidemic of HPAI A(H5N8) among poultry, in particular ducks, in France during 2016–2017 and to estimate the likely contribution of live-duck movements. Approximately 0.2% of live-duck movements could have been responsible for between-farm transmission events, mostly early during the epidemic. Results also suggest a transmission risk of 35.5% when an infected holding moves flocks to another holding within 14 days before detection. Finally, we found that densely connected groups of holdings with sparse connections between groups overlapped farmer organizations, which represents important knowledge for surveillance design. This study highlights the importance of movement bans in zones affected by HPAI and of understanding transmission routes to develop appropriate HPAI control strategies.

Infectious diseases commonly spread among animal premises by different transmission pathways, including live-animal movement networks that can cause outbreaks in widespread locations or through proximity networks, leading to spatial clusters of outbreaks (1,2). High-quality data on the spatial distribution of premises have enabled development of transmission models in which the proximity network assumes that any given

infectious premises can infect all susceptible premises within a geographic range (3,4). Increased monitoring of trade-related movement data has enabled the emergence of innovative modeling approaches based on social network analysis (5,6). Such an approach has been widely used to quantify how animal movement networks have contributed to disease transmission between animal premises (7,8) and relies on the assumption that premises intensively connected within the network are more likely to become infected and spread infection. Accordingly, efforts have focused on integrating movement and local spread components into models when the dynamics of past epidemics are explored, the effects of control strategies evaluated, and the pattern of future epidemics predicted (9–11).

However, the relative contribution of movement networks to the overall transmission risks remains poorly understood, compromising assessments of accurate and realistic disease spread modeling and control efforts. First, to assess the likelihood that the infection was acquired from movement networks, tracing of live-animal movement is required but might be challenging, especially in resource-poor settings where movement data are not regularly recorded as part of flock management systems. Then, the order or time at which animal premises become infected must be statistically related to their position in the movement network or in geographic space. However, these dates of infection are often inaccurate because reporting is delayed or completely lacking, particularly when tracking chronic diseases or wildlife populations, or when resources are limited.

In 2016–2017, Europe was hit hard by an unprecedented wave of highly pathogenic avian influenza (HPAI) A(H5N8) outbreaks that had severe socioeconomic effects on poultry production, global trade, and human livelihoods (12). Most outbreaks were reported in France; ducks were the most affected poultry species (13,14). The epidemic was contained by the end of March 2017 in France by timely application

Author affiliations: École Nationale Vétérinaire de Toulouse, Toulouse, France (C. Guinat, T. Vergne, T. Corre, J.-L. Guérin, M.C. Paul); Institut National de Recherche pour l'Agriculture, l'Alimentation et l'Environnement, Toulouse (C. Guinat, T. Vergne, T. Corre, J.-L. Guérin, M.C. Paul); Agence Nationale de Sécurité Sanitaire de l'Alimentation, Université Paris-Est, Maisons-Alfort, France (B. Durand); Direction Générale de l'Alimentation, Paris, France (S. Rautureau); Agence Nationale de Sécurité Sanitaire de l'Alimentation, Ploufragan, France (A. Scoizec, S. Lebouquin-Leneveu)

DOI: https://doi.org/10.3201/eid2603.190412
of measures after detection of the first outbreaks, as
provided for by European Union legislation (15-17).
These measures included culling all birds on the in-
fected holdings, establishing a 3-km protection zone
and 10-km surveillance zone with stringent ingoing
and outgoing movement bans, testing before move-
ments, and increasing biosecurity measures for hold-
ings in these zones. Spatiotemporal analysis of HPAI
outbreaks has shown that the disease spread was
partly driven by transmission events between poul-
try holdings in close proximity in space and time
(14). Although these previous results helped generate
hypotheses about possible routes of infection, they
did not enable weighting their relative contribution.
Duck movement networks were also identified as un-
derlying factors for the spatial distribution of HPAI
outbreaks (18), suggesting that these factors should
be considered to appropriately describe the epidemic
spread. Accurate data on the location and date of sus-
picion (i.e., onset of clinical signs and increased death
rates) of infected holdings and live-duck movements
between holdings were collected in France, provid-
ing a unique opportunity to unravel the spatial and
network dimensions of the epidemic. Our objectives
were to analyze live-duck movement networks dur-
ing the 2016-2017 H5N8 epidemic in France and in-
vestigate their likely contribution to disease spread.

Methods

Data Collection

Outbreak Data

We obtained data on the H5N8 outbreaks in ducks in France during the 2016–2017 epidemic (November 28, 2016–March 23, 2017) from the Direction Générale de l’Alimentation of the French Ministry of Agriculture (Paris, France). An outbreak was defined as detection of ≥1 H5N8-infected animal (confirmed by virus isolation or PCR) in a duck holding. Only outbreaks that occurred in holdings that sent or received duck flocks during the study period were retained for the analysis. Data comprised the list of laboratory-confirmed outbreaks, holding identification number, geographic locations (EPSG:2154/RGF93/Lambert-93 [https://epsg.io/2154]), and date of suspicion available by clinical or active surveillance.

Trade Movement Data

We considered only duck movements because they
represented the most affected poultry species (81.6%)
during the epidemic (14). The French organization of
fattening duck producers (Comité Interprofessionnel

des Palmipèdes à Foie Gras [CIFOG]) requires duck
producers to report movements from and onto their
holdings within 1-2 days of the movement. We thus
obtained data on live-duck movements and hold-
ings' characteristics from the professional database of
the CIFOG, under the appropriate confidential data
transfer agreements. Data included the list of move-
ment records (defined as movement of a flock be-
tween 2 different holdings on the same day), which
consisted of the date of movement, identification
number of the departure and arrival holdings, and
number of ducks moved. The incubation period (i.e.,
time between virus introduction and onset of clini-
cal signs) ranges from ≈1 to ≈5 days at the individual
level and could be longer at the flock level because
of the transmission process (19). Because such dura-
tion is difficult to estimate, a 14-day incubation pe-
riod was assumed at the flock level (we conducted a
sensitivity analysis using a 21-day incubation period
and showed that it did not affect the results). Conse-
quently, movements within 14 days before the detec-
tion of an infected holding might be responsible for
between-holding transmission events (19). Thus, we
retained only movement data during November 1,
2016-March 31, 2017, between holdings for the analy-
sis; movements to slaughterhouses were excluded.
Holdings' characteristics included the geographic
locations (EPSG:2154/RGF93/Lambert-93), group of
farmer organization, and type of production: rearing
(1-day-old ducklings are reared for ≈3 weeks), breed-
ing (1-day to 3-week-old ducks are bred for ≈9-12
weeks), and force-feeding (12-week-old ducks are
force-fed for ≈12 days). For holdings with no avail-
able coordinates (9.5%), we used the coordinates of
the center of the commune (smallest administrative
unit in France, with a median area of 10 km²).

Data Analysis

Spatiotemporal Description of Movements

We first generated descriptive statistics for the num-
ber of active holdings (i.e., holdings that received or
sent ducks during the study period), the number of
flocks moved, and the distances covered by move-
ments (i.e., using Euclidean distance in kilometers
between the departure and arrival holdings) per pair
of holdings. We removed holdings without available
coordinates from the Euclidean distance estimations.
Finally, we mapped the number of movements from/
to holdings between departments by aggregating
movements at the department level (administrative
unit in France corresponding to NUTS [Nomencla-
ture of Territorial Units for Statistics] level 3).
## Network Analysis

We built directed and weighted networks for data from November 1, 2016–March 31, 2017, considering each duck holding as a node and a movement of a flock between 2 holdings as an edge. We assigned directions to each edge according to the date on which ducks were moved between 2 nodes and assigned weights to each edge according to the number of ducks moved between 2 nodes. We identified trade communities (i.e., densely connected groups of nodes, with only sparse connections between groups [20]) over the whole study period using a walktrap algorithm (21) based on random walks through the edges in the network. We selected the 15 largest communities on the basis of their respective numbers of holdings and typed them according to holding production types. We mapped holdings belonging to the 15 largest communities and performed a bootstrapped version of the Fisher exact test with 10,000 replicates (22) to test whether dependence existed between the trade community and the organization to which farmers belong.

Next, we assessed the likely contribution of live-duck movements in the distribution of H5N8 outbreaks in the network using a permutation-based approach (23–25). The rationale behind this approach was that if the outbreaks resulted from disease spread through the movement networks, the mean number of infected holdings in contact with an infected holding in the network would be significantly greater than expected if infected holdings were randomly distributed in the network. Again, duck holdings were assumed to become infected through the movement networks if they had received movements from infected holdings within an at-risk period of 14 days before their date of suspicion. Hence, we assumed the mean number of potential transmission events through the movement networks corresponded with the mean number of at-risk movements defined by movements originating from an infected holding (the sender) within 14 days before its date of suspicion and directed to a distinct infected duck holding (the receiver) within 14 days before the receiver's date of suspicion. We then compared this statistic (i.e., the mean number of transmission events per infected holding) with the distribution of the expected statistic under the null hypothesis according to which the

dates of suspicion were randomly distributed among infected holdings in the network (n = 1,000), with the p value corresponding to the proportion of permutations for which the expected statistic is higher than the observed statistic. Similarly, to assess the role of proximity networks, we also conducted the test by calculating the following statistic: the mean number of infected duck holdings close in time (differences of infection dates within 14 days) and space (both located within a 10-km radius [14]) per infected duck holding. We selected this space-time window on the basis of previous spatiotemporal analysis conducted on the dataset (14). Finally, we identified the likely origins of holding infections by calculating the proportion of infected duck holdings retrieved as receivers in the list of transmission events through the movement network and the proportion of infected duck holdings for which ≥1 infected duck holding close in time and space was retrieved in the proximity network. On the basis of the movement and proximity networks, we thus attributed to each holding a likely origin of infection as follows: ingoing edge in the movement network only, ingoing edge in the proximity network only, ingoing edges in both movement and proximity networks, and no ingoing edge (i.e., other transmission pathways than by movement and proximity; for example, by introduction of infected migratory birds from northern Eurasia [13,26]). Because movement bans were reinforced on February 2, 2017 (27), we retained only movement and outbreak data for November 1, 2016–February 2, 2017, for this analysis. We conducted all analyses in R statistical software version 3.4.2 using the igraph package (28).

## Results

### Spatiotemporal Description

A total of 9,096 movements, involving 10,945,388 ducks moved among 2,098 holdings, occurred during November 1, 2016–March 31, 2017 (Table 1). Most holdings involved in these movements were characterized as force-feeding (48.8%), followed by breeding (35.7%) and breeding plus force-feeding (11.9%). The holdings were located mainly in southwestern and northwestern France (Appendix Figure 1, https://wwwnc.cdc.gov/EID/article/26/3/19-0412-App1.pdf). Overall, most (95.8%) of the flocks were moved

<table><caption>Table 1. Descriptive statistics of duck movements per pair of holdings, France, November 1, 2016–March 31, 2017*</caption><thead><tr><th rowspan="2">Holding type pair</th><th rowspan="2">No. (%) flocks moved</th><th colspan="4">No. ducks moved</th><th colspan="4">Distance moved, km</th></tr><tr><th>Mean</th><th>Median</th><th>IQR</th><th>Max</th><th>Mean</th><th>Median</th><th>IQR</th><th>Max</th></tr></thead><tbody><tr><td>Rearing to breeding</td><td>382 (4.2)</td><td>6,001</td><td>4,773</td><td>3,016–8,991</td><td>15,090</td><td>58</td><td>36</td><td>0.1–101</td><td>213</td></tr><tr><td>Breeding to force-feeding</td><td>8,712 (95.8)</td><td>993</td><td>958</td><td>629–1,188</td><td>8,050</td><td>50</td><td>40</td><td>16–71</td><td>408</td></tr></tbody></table>

<footnote>*IQR, interquartile range; max, maximum.</footnote>
from breeding to force-feeding holdings; only 4.2% of movements occurred from rearing to breeding holdings (Table 1). However, more ducks were moved from rearing to breeding holdings (median 4,773) than from breeding to force-feeding holdings (median 958). Movements clearly clustered in the 2 separate geographic areas, southwestern and northwestern France; a limited amount of movements occurred between these 2 areas (Appendix Figure 2).

## Network Analysis

The network analysis identified 99 trade communities comprising 2,098 holdings during November 1, 2016–March 31, 2017. The 15 largest communities in terms of number of holdings included 91.8% of holdings. These communities showed a relatively distinct spatial distribution in northwestern France but completely overlapped in southwestern France (Figure 1). However, the communities were characterized by similar holding compositions, dominated by breeding and force-feeding holdings. The 15 largest communities overlapped significantly with the 15 largest groups of farmer organizations (p<0.001) (Appendix Figure 3): For example, community 1

included 80.0% of holdings belonging to organization A, community 2 included 51.7% of holdings belonging to organization I and 44.8% to organization J, and community 5 included 63.1% and community 12 36.0% of organization B.

A total of 6,521 movements between 1,988 holdings (involving 104 infected holdings) occurred during November 1, 2016–February 2, 2017. Among the 104 infected holdings, 40 (38.5%) were identified as senders, 36 (34.6%) as receivers, and 28 (26.9%) as senders and receivers during that period. Most (989 [91.8%] of 6,521) movements occurred between non-infected holdings. We identified 16 (0.2%) of 6,521 movements as at risk (i.e., they were compatible with transmission events through the movement networks) (Figure 2). These movements mostly occurred between breeding and force-feeding holdings from the end of November through the beginning of January, before stringent movement bans were implemented, and were directed to areas where most outbreaks were reported during the following weeks (Figure 2). Some of the at-risk movements originated from the first outbreak, reported at the beginning of the epidemic (end of November 2016). A

<figure><img src="image_2.png" /><figcaption>**Figure 1.** Spatial distribution of the 15 largest live-duck trade communities, France, November 1, 2016–March 31, 2017.</figcaption></figure>
**Figure 2.** Spatial distribution of live-duck movements identified as responsible for highly pathogenic avian influenza A(H5N8) transmission events between holdings through the movement networks, France, November 1, 2016–February 2, 2017.

<figure><img src="image_3.png" /></figure>

few movements (0.4%, 29/6,521) occurred between infected holdings and holdings that did not become infected within 14 days after the movements, from the end of November through the beginning of January, before stringent movement bans were implemented. Therefore, transmission risk through live-duck movements was estimated at 35.5% (16/ [16 + 29]), meaning that the likelihood of infection when an infected holding moved flocks to another holding within 14 days before detection was 35.5%. Results from the permutation-based approach indicated the mean number of transmission events per infected holding was significantly greater than under the null hypothesis of an absence of association between movement and infection status (according to which the dates of suspicion should be randomly distributed among network nodes) (p<0.001). Moreover, the mean number of infected holdings close in time and space per infected holding was also significantly greater than expected (p<0.001). By retrieving holding receivers in the list of transmission events through the movement and proximity networks, most sources of holding infection were attributed to proximity networks (66.3%), followed by movement networks (14.4%), and other unknown means

of transmission were possible (23.1%) (Table 2). The 16 at-risk movements could be the likely source of infection for only 15 farms because 1 infected farm received 2 at-risk movements.

## Discussion

Using a detailed analysis of live-duck movements and proximity networks, we unraveled the underlying transmission processes of the H5N8 epidemic in ducks in France during 2016–2017. During November 1, 2016–March 31, 2017, which overlaps the H5N8 epidemic period, we observed the most movements from breeding to force-feeding holdings and the largest duck flocks from the rearing to the breeding stage. These findings are consistent with the production cycle and the high specialization in production within which the number of rearing holdings where ducks are first reared and then sent as large flocks to breeding holdings is limited. Flocks are then divided into small flocks to be moved to force-feeding holdings, resulting in a large number of force-feeding holdings reported in the country. Assuming that movement networks influence disease spread, this structure becomes important in terms of disease prevention and control: the dominant role of such superreceiver and
<table><caption>**Table 2.** Contribution of movement and proximity networks for highly pathogenic avian influenza A(H5N8) transmission events between live-duck holdings, France, November 1, 2016–February 2, 2017</caption><thead><tr><th rowspan="2">Origin of infection</th><th colspan="2">Infected holdings, no. (%), n = 104</th></tr><tr><th>14 d before date of suspicion</th><th>21 d before date of suspicion</th></tr></thead><tbody><tr><td>Movement network</td><td>11 (10.6)</td><td>11 (10.6)</td></tr><tr><td>Proximity network</td><td>65 (62.5)</td><td>72 (69.2)</td></tr><tr><td>Movement and proximity networks</td><td>4 (3.8)</td><td>4 (3.8)</td></tr><tr><td>Other</td><td>24 (23.1)</td><td>17 (16.3)</td></tr></tbody></table>

superspreader holdings indicates that monitoring only a few holdings would be sufficient to reduce disease spread or that targeting sampling in these high-risk holdings would be more effective than random sampling when time and resources are limited (29,30). This structure supports the recent active surveillance campaign of duck flocks before movements between these 2 production stages, implemented as a result of the devastating H5N8 epidemic (31). Overall, most movements were short range (50% cover <40 km and 75% <75 km). This finding is consistent with results from a spatiotemporal analysis (14), which provided evidence that local transmission processes mainly drove the spread. Moreover, our study demonstrated that movements clustered mainly in 2 geographic areas (southwestern and northwestern France) and that a limited number of movements occurred between these 2 areas, potentially explaining why the disease did not spread from south to northwest (14).

The 15 largest trade communities that comprised most (91.8%) holdings clearly overlapped with the 15 largest farmer organizations. Again, this finding is crucial in terms of disease surveillance because it highlights that targeting sampling of holdings belonging to the trade community of infected holdings would be more effective than random sampling to prevent further disease spread. In terms of disease control, these results indicate that trade within a given group of highly connected holdings could be maintained by disrupting epidemiologic links to other groups of holdings at risk, mainly to minimize disruption of global trade during an epidemic (32,33). Moreover, being part of a particular farmer organization implies that holdings are connected by other means than movements of live birds, such as shared transport, equipment, feed, animal staff, or catching teams, that could also facilitate transmission events within the community. Our study also highlighted the important role of the community structure in spreading H5N8: the community to which holdings belong (and thus the farmer organization) was significantly associated with the H5N8 holding infection status (data not shown). Again, trade communities did not overlap between northwestern and southwestern France, which could explain why most of outbreaks

remained clustered in southwestern France during the 2016–2017 epidemic (14).

Results from our permutation-based approach suggested that a limited proportion of holdings (14.4%) became infected through the movement networks before February 2017. We identified some of these transmission events in the movement networks as originating from the first outbreak reported at the beginning of the epidemic (end of November 2016), before stringent movement bans were implemented, and directed to areas where most of the outbreaks were reported during the following weeks (14,34). Therefore, despite their low number, live-duck movements might have played a crucial role in the onset and spatial extent of the 2016–2017 H5N8 epidemic in the country. The limited contribution of movement networks to disease spread is most likely explained by the timely implementation of control strategies and movement bans after the first outbreaks were detected (15,16). This limited contribution also is most likely attributed to the duck production characteristics, highly specialized holdings organized in a small pyramidal structure. Results suggest a transmission risk of 35.5% when an infected holding moves flocks to another holding within 14 days before detection. These findings support efforts by authorities in France in collaboration with the farmer organizations to enhance biosecurity during the transport of ducks (31) after successive waves of HPAI outbreaks within 2 years (14,35). Trucks moving flocks are not allowed to load from several different holdings to minimize the risk for contact infections as trucks travel between holdings. It is likely officials will implement new rules, such as using different sets of trucks and cages to move flocks from breeding to force-feeding and from force-feeding to slaughter. A higher proportion of holdings (66.3%) became infected through proximity networks, consistent with previous work that identified local spread as a predominant transmission pathway in the early stage of the epidemic, that is, before February 2017 (14). As a result, these findings also support the national biosecurity program that was implemented to prevent the introduction and spread of poultry diseases at the holding level (36,37).
The 2 recent devastating epidemics of HPAI in France (2015–2016 and 2016–2017) led to major changes in the collection of movement data. Specifically, the farmer organizations require duck producers to timely and accurately report any details on flock movements, leading to the expectation that underreporting remains limited. Data regarding transport, shared equipment, feed, animal staff, carcass rendering, catching teams, or wild birds were not available (38,39). However, these transmission pathways might be partly reflected by the proximity network (for example, neighboring holdings might share the same equipment or carcass rendering round) or by 23.1% of holdings for which the infection origin was attributed to transmission pathways other than movement or proximity networks (for example, by introduction of infected migratory birds from northern Eurasia [13,26]). Although the epidemiologic mechanisms that could explain some of these transmission events remain to be explored, one could infer that these transmission pathways might have played a larger role in the spread of H5N8 between holdings than movement of live ducks. Recent studies have shown that wild birds are likely to have played a minor role in the spread of H5N8 between holdings (18,40), suggesting that the main driver of the epidemic was holding-to-holding transmission. Further work will compare these results with movement networks during a period with no outbreaks reported as to how outbreaks and intervention strategies have modified the structure of the movement networks.

This study provides insights into the likely contribution of live-duck movement networks into the spread of H5N8 at the beginning of the 2016–2017 epidemic in France. This study also highlights the importance of movement bans in affected zones and that understanding transmission routes is paramount for developing appropriate control strategies for HPAI. A new aspect of this study is the inclusion of a permutation-based approach based on the dates of holding infection to evaluate whether the acquisition of holding infection was consistent with virus transmission through the network. This approach has been limitedly applied in the epidemiology of infectious diseases (23–25), although it outperforms other degree-based statistical methods, such as logistic regression and nonparametric tests. Outcomes about the relative contribution of movement and proximity networks represent a required basis on which predictive models of HPAI spread could be developed. Finally, this study emphasizes the importance of supplementing epidemiologic data with animal movement data and therefore calls for collaborative efforts to report trade

movement data and make them available for appropriately targeting surveillance and interventions during future outbreaks.

## Acknowledgments

We gratefully acknowledge the Direction Générale de l'Alimentation of the French Ministry of Agriculture and the Comité Interprofessionnel des Palmipèdes à Foie Gras (CIFOG, France) for providing the data, particularly Marie Laborde for providing additional detailed information about the data and interpreting the results.

This work was conducted within the framework of the Chaire de Biosécurité at the École Nationale Vétérinaire de Toulouse, which is funded by the French Ministry of Agriculture and supported by the Unité Mixte de Recherche Interactions Hôtes Agents Pathogènes 1225 École Nationale Vétérinaire de Toulouse, Institut National de la Recherche Agronomique. The research leading to the results was funded by the People Programme (Marie Curie Actions) of the European Union's Seventh Framework Programme (FP7/2007-2013) under Research Executive Agency grant agreement no. PCOFUND-GA-2013-609102, through the PRESTIGE program coordinated by Campus France. This work was also financially supported by the FEDER/Région Occitanie Recherche et Sociétés 2018 – AI-TRACK.

## About the Author

Dr. Guinat is a postdoctoral researcher in veterinary epidemiology at the École Nationale Vétérinaire de Toulouse (France). Her major research interest is understanding the spread and maintenance of emerging infectious diseases.

## References

1. Fèvre EM, Bronsvoort BMC, Hamilton KA, Cleaveland S. Animal movements and the spread of infectious diseases. Trends Microbiol. 2006;14:125-31. http://dx.doi.org/10.1016/j.tim.2006.01.004

2. Green DM, Kiss IZ, Mitchell AP, Kao RR. Estimates for local and movement-based transmission of bovine tuberculosis in British cattle. Proc Biol Sci. 2008;275:1001-5. http://dx.doi.org/10.1098/rspb.2007.1601

3. Ferguson NM, Donnelly CA, Anderson RM. The foot-and-mouth epidemic in Great Britain: pattern of spread and impact of interventions. Science. 2001;292:1155-60. http://dx.doi.org/10.1126/science.1061020

4. Keeling MJ, Woolhouse MEJ, Shaw DJ, Matthews L, Chase-Topping M, Haydon DT, et al. Dynamics of the 2001 UK foot and mouth epidemic: stochastic dispersal in a heterogeneous landscape. Science. 2001;294:813-7. http://dx.doi.org/10.1126/science.1065973

5. Eubank S, Guclu H, Kumar VS, Marathe MV, Srinivasan A, Toroczkai Z, et al. Modelling disease outbreaks in realistic urban social networks. Nature. 2004;429:180-4. http://dx.doi.org/10.1038/nature02541

6. Riley S, Ferguson NM. Smallpox transmission and control: spatial dynamics in Great Britain. Proc Natl Acad
Sci U S A. 2006;103:12637-42. http://dx.doi.org/10.1073/pnas.0510873103

7. Kao RR, Danon L, Green DM, Kiss IZ. Demographic structure and pathogen dynamics on the network of livestock movements in Great Britain. Proc Biol Sci. 2006;273:1999-2007. http://dx.doi.org/10.1098/rspb.2006.3505

8. Kiss IZ, Green DM, Kao RR. The network of sheep movements within Great Britain: Network properties and their implications for infectious disease spread. J R Soc Interface. 2006;3:669-77. http://dx.doi.org/10.1098/rsif.2006.0129

9. Fournié G, Guitian J, Desvaux S, Cuong VC, Dung H, Pfeiffer DU, et al. Interventions for avian influenza A (H5N1) risk management in live bird market networks. Proc Natl Acad Sci U S A. 2013;110:9177-82. http://dx.doi.org/10.1073/pnas.1220815110

10. Green DM, Kiss IZ, Kao RR. Modelling the initial spread of foot-and-mouth disease through animal movements. Proc Biol Sci. 2006;273:2729-35. http://dx.doi.org/10.1098/rspb.2006.3648

11. Keeling MJ, Eames KTD. Networks and epidemic models. J R Soc Interface. 2005;2:295-307. http://dx.doi.org/10.1098/rsif.2005.0051

12. Napp S, Majó N, Sánchez-González R, Vergara-Alert J. Emergence and spread of highly pathogenic avian influenza A(H5N8) in Europe in 2016-2017. Transbound Emerg Dis. 2018;65:1217-26. http://dx.doi.org/10.1111/tbed.12861

13. Alarcon P, Brouwer A, Venkatesh D, Duncan D, Dovas CI, Georgiades G, et al. Comparison of 2016-17 and previous epizootics of highly pathogenic avian influenza H5 Guangdong lineage in Europe. Emerg Infect Dis. 2018;24:2270-83. http://dx.doi.org/10.3201/eid2412.171860

14. Guinat C, Nicolas G, Vergne T, Bronner A, Durand B, Courcoul A, et al. Spatio-temporal patterns of highly pathogenic avian influenza virus subtype H5N8 spread, France, 2016 to 2017. Euro Surveill. 2018;23:23. http://dx.doi.org/10.2807/1560-7917.ES.2018.23.26.1700791

15. Direction Générale de l'Alimentation. Arrêté du 4 janvier 2017 définissant les zones géographiques dans lesquelles un abattage préventif est ordonné en application de l'arrêté du 4 janvier 2017 relatif aux mesures complémentaires techniques et financières pour la maîtrise de l'épizootie d'influenza aviaire due au virus H5N8 dans certains départements [cited 2019 May 20]. https://www.legifrance.gouv.fr/affichTexte.do?cidTexte=JORFTEXT000033793904

16. Direction Générale de l'Alimentation. Arrêté du 31 mars 2017 déterminant des dispositions de prévention, de surveillance et lutte complémentaires contre l'influenza aviaire hautement pathogène dans certaines parties du territoire [cited 2019 May 20]. https://www.legifrance.gouv.fr/affichTexte.do?cidTexte=JORFTEXT000034330812

17. Direction Générale de l'Alimentation. Arrêté du 18 janvier 2008 fixant des mesuretechniques et administratives relatives à la lutte contre l'influenza aviaire [cited 2019 Dec 15]. https://www.legifrance.gouv.fr/affichTexte.do?cidTexte=JORFTEXT000017958300

18. Guinat C, Artois J, Bronner A, Guérin JL, Gilbert M, Paul MC. Duck production systems and highly pathogenic avian influenza H5N8 in France, 2016-2017. Sci Rep. 2019;9:6177. http://dx.doi.org/10.1038/s41598-019-42607-x

19. Pantin-Jackwood MJ, Costa-Hurtado M, Shepherd E, Dejesus E, Smith D, Spackman E, et al. Pathogenicity and transmission of H5 and H7 highly pathogenic avian influenza viruses in mallards. J Virol. 2016;90:9967-82. http://dx.doi.org/10.1128/JVI.01165-16

20. Newman MEJ. Modularity and community structure in networks. Proc Natl Acad Sci U S A. 2006;103:8577-82. http://dx.doi.org/10.1073/pnas.0601602103

21. Pons P, Latapy M. Computing communities in large networks using random walks. In: Computer and information sciences. International Symposium on Computer and Information Sciences [cited 2017 Jun 15]. https://link.springer.com/chapter/10.1007/11569596_31

22. Lin J-J, Chang C-H, Pal N. A revisit to contingency table and tests of independence: bootstrap is preferred to Chi-square approximations as well as Fisher's exact test. J Biopharm Stat. 2015;25:438-58. http://dx.doi.org/10.1080/10543406.2014.920851

23. Bouchez-Zacria M, Courcoul A, Durand B. The distribution of bovine tuberculosis in cattle farms is linked to cattle trade and badger-mediated contact networks in south-western France, 2007-2015. Front Vet Sci. 2018;5:173. http://dx.doi.org/10.3389/fvets.2018.00173

24. Soares Magalhães RJ, Ortiz-Pelaez A, Thi KLL, Dinh QH, Otte J, Pfeiffer DU. Associations between attributes of live poultry trade and HPAI H5N1 outbreaks: a descriptive and network analysis study in northern Vietnam. BMC Vet Res. 2010;6:10. http://dx.doi.org/10.1186/1746-6148-6-10

25. VanderWaal K, Enns EA, Picasso C, Packer C, Craft ME. Evaluating empirical contact networks as potential transmission pathways for infectious diseases. J R Soc Interface. 2016;13:20160166. http://dx.doi.org/10.1098/rsif.2016.0166

26. Sims L, Harder TC, Brown I, Gaidet N, Belot G, von Dobschuetz S, et al. Highly pathogenic H5 avian influenza in 2016 and early 2017—observations and future perspectives. Empres Focus On, 2017;11:1-14 [cited 2019 May 20]. http://www.fao.org/3/a-i8068e.pdf

27. Agence nationale de Sécurité Sanitaire de l'Alimentation, de l'Environnement et du Travail. AVIS de l'Agence nationale de sécurité sanitaire de l'alimentation, de l'environnement et du travail relatif au « périmètre optimal de dépeuplement préventif influenza aviaire IA HP H5N8 ≈ [cited 2017 Apr 12]. https://www.anses.fr/fr/system/files/SABA2017SA0011.pdf

28. Csardi G, Nepusz T. The igraph software package for complex network research. InterJournal, Complex Systems. 2006;1695-1-9.

29. Lloyd AL, May RM. Epidemiology. How viruses spread among computers and people. Science. 2001;292:1316-7. http://dx.doi.org/10.1126/science.1061076

30. Pastor-Satorras R, Vespignani A. Epidemic spreading in scale-free networks. Phys Rev Lett. 2001;86:3200-3. http://dx.doi.org/10.1103/PhysRevLett.86.3200

31. Direction Générale de l'Alimentation. Le Pacte de lutte contre l'Influenza Aviaire et de relance de la filière foie gras [cited 2017 May 10]. http://agriculture.gouv.fr/le-pacte-de-lutte-contre-linfluenza-aviaire-et-de-relance-de-la-filiere-foie-gras

32. Ratananakorn L, Wilson D. Zoning and compartmentalisation as risk mitigation measures: an example from poultry production. Rev Sci Tech. 2011;30:297-307. http://dx.doi.org/10.20506/rst.30.1.2029

33. Scott A, Zepeda C, Garber L, Smith J, Swayne D, Rhorer A, et al. The concept of compartmentalisation. Rev Sci Tech. 2006;25:873-9, 881-7, 889-95. http://dx.doi.org/10.20506/rst.25.3.1702

34. Direction Générale de l'Alimentation. Situation de l'Influenza aviaire HP H5N8 en France au 05122016 [cited 2018 Nov 8]. https://www.plateforme-esa.fr/article/situation-de-l-influenza-aviaire-hp-h5n8-en-france-au-05122016
35. Briand F-X, Schmitz A, Ogor K, Le Prioux A, Guillou-Cloarec C, Guillemoto C, et al. Emerging highly pathogenic H5 avian influenza viruses in France during winter 2015/16: phylogenetic analyses and markers for zoonotic potential. Euro Surveill. 2017;22:30473. http://dx.doi.org/10.2807/1560-7917.ES.2017.22.9.30473

36. Delpont M, Blondel V, Robertet L, Duret H, Guerin J-L, Vaillancourt J-P, et al. Biosecurity practices on foie gras duck farms, Southwest France. Prev Vet Med. 2018;158:78-88. http://dx.doi.org/10.1016/j.prevetmed.2018.07.012

37. de l'Alimentation DG. Arrêté du 8 février 2016 relatif aux mesures de biosécurité applicables dans les exploitations de volailles et d'autres oiseaux captifs dans le cadre de la prévention contre l'influenza aviaire [cited 2018 Mar 28]. https://www.legifrance.gouv.fr/affichTexte.do?cidTexte=J0 RFTEXT000032000273

38. 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 Vet Res. 2008;4:27. http://dx.doi.org/10.1186/1746-6148-4-27

39. Nickbakhsh S, Matthews L, Dent JE, Innocent GT, Arnold ME, Reid SWJ, et al. Implications of within-farm transmission for network dynamics: consequences for the spread of avian influenza. Epidemics. 2013;5:67-76. http://dx.doi.org/10.1016/j.epidem.2013.03.001

40. Andronico A, Courcoul A, Bronner A, Scoizec A, Lebouquin-Leneveu S, Guinat C, et al. Highly pathogenic avian influenza H5N8 in south-west France 2016-2017: A modeling study of control strategies. Epidemics. 2019; 28:100340. http://dx.doi.org/10.1016/j.epidem.2019.03.006

Address for correspondence: Claire Guinat, École Nationale Vétérinaire de Toulouse, UMR INRAE-ENVT 1225–IHAP, 23 chemin des Capelles, 31300 Toulouse, France; email: c.guinat@envt.fr

# The Public Health Image Library (PHIL)

<figure><img src="image_4.png" /></figure>

The Public Health Image Library (PHIL), Centers for Disease Control and Prevention, contains thousands of public health–related images, including high-resolution (print quality) photographs, illustrations, and videos.

PHIL collections illustrate current events and articles, supply visual content for health promotion brochures, document the effects of disease, and enhance instructional media.

PHIL images, accessible to PC and Macintosh users, are in the public domain and available without charge.

Visit PHIL at: http://phil.cdc.gov/phil