Crossroads of highly pathogenic H5N1: overlap between wild and domestic birds in the Black Sea-Mediterranean impacts global transmission.

Understanding transmission dynamics that link wild and domestic animals is a key element of predicting the emergence of infectious disease, an event that has highest likelihood of occurring wherever human livelihoods depend on agriculture and animal trade. Contact between poultry and wild birds is a key driver of the emergence of highly pathogenic avian influenza (HPAI), a process that allows for host switching and accelerated reassortment, diversification, and spread of virus between otherwise unconnected regions. This study addresses questions relevant to the spillover of HPAI at a transmission hotspot: what is the nature of the wild bird-poultry interface in Egypt and adjacent Black Sea-Mediterranean countries and how has this contributed to outbreaks occurring worldwide? Using a spatiotemporal model of infection risk informed by satellite tracking of waterfowl and viral phylogenetics, this study identified ecological conditions that contribute to spillover in this understudied region. Results indicated that multiple ducks (Northern Shoveler and Northern Pintail) hosted segments that shared ancestry with HPAI H5 from both clade 2.2.1 and clade 2.3.4 supporting the role of Anseriformes in linking viral populations in East Asia and Africa over large distances. Quantifying the overlap between wild ducks and H5N1-infected poultry revealed an increasing interface in late winter peaking in early spring when ducks expanded their range before migration, with key differences in the timing of poultry contact risk between local and long-distance migrants.

Do contrasting patterns of migration movements and disease outbreaks between congeneric waterfowl species reflect differing immunity?

Long-distance migrations influence the dynamics of hostpathogen interactions and understanding the role of migratory waterfowl in the spread of the highly pathogenic avian influenza viruses (HPAIV) is important. While wild geese have been associated with outbreak events, disease ecology of closely related species has not been studied to the same extent. The swan goose (Anser cygnoides) and the bar-headed goose (Anser indicus) are congeneric species with distinctly different HPAIV infection records; the former with few and the latter with numerous records. We compared movements of these species, as well as the more distantly related whooper swan (Cygnus cygnus) through their annual migratory cycle to better understand exposure to HPAIV events and how this compares within and between congeneric and noncongeneric species. In spite of their record of fewer infections, swan geese were more likely to come in contact with disease outbreaks than bar-headed geese. We propose two possible explanations: i) frequent prolonged contact with domestic ducks increases innate immunity in swan geese, and/or ii) the stress of high-elevation migration reduces immunity of bar-headed geese. Continued efforts to improve our understanding of species-level pathogen response is critical to assessing disease transmission risk.

Incorporating heterogeneous sampling probabilities in continuous phylogeographic inference - Application to H5N1 spread in the Mekong region.

MOTIVATION: The potentially low precision associated with the geographic origin of sampled sequences represents an important limitation for spatially explicit (i.e. continuous) phylogeographic inference of fast-evolving pathogens such as RNA viruses. A substantial proportion of publicly available sequences is geo-referenced at broad spatial scale such as the administrative unit of origin, rather than more precise locations (e.g. geographic coordinates). Most frequently, such sequences are either discarded prior to continuous phylogeographic inference or arbitrarily assigned to the geographic coordinates of the centroid of their administrative area of origin for lack of a better alternative. RESULTS: We here implement and describe a new approach that allows to incorporate heterogeneous prior sampling probabilities over a geographic area. External data, such as outbreak locations, are used to specify these prior sampling probabilities over a collection of sub-polygons. We apply this new method to the analysis of highly pathogenic avian influenza H5N1 clade data in the Mekong region. Our method allows to properly include, in continuous phylogeographic analyses, H5N1 sequences that are only associated with large administrative areas of origin and assign them with more accurate locations. Finally, we use continuous phylogeographic reconstructions to analyse the dispersal dynamics of different H5N1 clades and investigate the impact of environmental factors on lineage dispersal velocities. AVAILABILITY AND IMPLEMENTATION: Our new method allowing heterogeneous sampling priors for continuous phylogeographic inference is implemented in the open-source multi-platform software package BEAST 1.10. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Optimising the detectability of H5N1 and H5N6 highly pathogenic avian influenza viruses in Vietnamese live-bird markets.

Live bird markets (LBMs) are major targets for avian influenza virus (AIV) surveillance programmes. While sampling the LBM environment has become a widely used alternative to the labour-intensive sampling of live poultry, the design of surveillance programmes and the interpretation of their results are compromised by the lack of knowledge about the effectiveness of these sampling strategies. We used latent class models and a unique empirical dataset collated in Vietnamese LBMs to estimate the sensitivity and specificity of five different sample types for detecting AIVs subtypes H5N1 and H5N6: oropharyngeal duck samples, solid and liquid wastes, poultry drinking water and faeces. Results suggest that the sensitivity of environmental samples for detecting H5N1 viruses is equivalent to that of oropharyngeal duck samples; however, taking oropharyngeal duck samples was estimated to be more effective in detecting H5N6 viruses than taking any of the four environmental samples. This study also stressed that the specificity of the current surveillance strategy in LBMs was not optimal leading to some false positive LBMs. Using simulations, we identified 42 sampling strategies more parsimonious than the current strategy and expected to be highly sensitive for both viruses at the LBM level. All of these strategies involved the collection of both environmental and oropharyngeal duck samples.

Investigating Home Range, Movement Pattern, and Habitat Selection of Bar-headed Geese during Breeding Season at Qinghai Lake, China.

The Bar-headed Goose is the only true goose species or Anserinae to migrate solely within the Central Asian Flyway, and thus, it is an ideal species for observing the effects of both land use and climate change throughout the flyway. In this paper, we investigate the home range, movement pattern, and habitat selection of Bar-headed Geese (Anser indicus) during the breeding season at Qinghai Lake, which is one of their largest breeding areas and a major migration staging area in the flyway. We identified several areas used by the geese during the breeding season along the shoreline of Qinghai Lake and found that most geese had more than one core use area and daily movements that provided insight into their breeding activity. We also observed the intensive use of specific wetlands and habitats near Qinghai Lake. These data provide interesting insights into the movement ecology of this important species and also provide critical information for managers seeking to understand and respond to conservation concerns threatening Bar-headed Geese, such as landscape and habitat changes.

Comparative Epidemiology of Highly Pathogenic Avian Influenza Virus H5N1 and H5N6 in Vietnamese Live Bird Markets: Spatiotemporal Patterns of Distribution and Risk Factors.

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

A Stakeholder Survey on Live Bird Market Closures Policy for Controlling Highly Pathogenic Avian Influenza in Vietnam.

Extensive research in Vietnam and elsewhere has shown that live bird markets (LBMs) play a significant role in the ecology and zoonotic transmission of avian influenzas (AIs) including H5N1 and H7N9. Vietnam has a large number of LBMs reflecting the consumer preferences for live poultry. Under pressure to mitigate risks for H7N9 and other zoonotic AIs, Vietnam is considering, among other mitigation measures, temporary closures of LBMs as a policy to reduce risk of AI outbreaks. However, the efficacy of market closure is debated, particularly because little is known about how poultry traders may react, and whether trading may emerge outside formal marketplaces. Combining efforts of anthropologists, economists, sociologists, and veterinarians can be useful to elucidate the drivers behind poultry traders' reactions and better understanding the barriers to implementing risk mitigation measures. In this paper, we present results from a stakeholder survey of LBM stakeholders in Vietnam. Our qualitative data show that trading outside formal markets is very likely to occur in the event of a temporary LBM market closure. Our data show that the poultry value chain in Vietnam remains highly flexible, with traders willing and able to trade poultry in many possible locations. Our results indicate that simplification of the poultry value chain along with strict enforcement, engagement of stakeholders, and adequate communication would be a necessary prerequisite before market closure could be an effective policy.

Could Changes in the Agricultural Landscape of Northeastern China Have Influenced the Long-Distance Transmission of Highly Pathogenic Avian Influenza H5Nx Viruses?

In the last few years, several reassortant subtypes of highly pathogenic avian influenza viruses (HPAI H5Nx) have emerged in East Asia. These new viruses, mostly of subtype H5N1, H5N2, H5N6, and H5N8 belonging to clade 2.3.4.4, have been found in several Asian countries and have caused outbreaks in poultry in China, South Korea, and Vietnam. HPAI H5Nx also have spread over considerable distances with the introduction of viruses belonging to the same 2.3.4.4 clade in the U.S. (2014-2015) and in Europe (2014-2015 and 2016-2017). In this paper, we examine the emergence and spread of these new viruses in Asia in relation to published datasets on HPAI H5Nx distribution, movement of migratory waterfowl, avian influenza risk models, and land-use change analyses. More specifically, we show that between 2000 and 2015, vast areas of northeast China have been newly planted with rice paddy fields (3.21 million ha in Heilongjiang, Jilin, and Liaoning) in areas connected to other parts of Asia through migratory pathways of wild waterfowl. We hypothesize that recent land use changes in northeast China have affected the spatial distribution of wild waterfowl, their stopover areas, and the wild-domestic interface, thereby altering transmission dynamics of avian influenza viruses across flyways. Detailed studies of the habitat use by wild migratory birds, of the extent of the wild-domestic interface, and of the circulation of avian influenza viruses in those new planted areas may help to shed more light on this hypothesis, and on the possible impact of those changes on the long-distance patterns of avian influenza transmission.

Global mapping of highly pathogenic avian influenza H5N1 and H5Nx clade 2.3.4.4 viruses with spatial cross-validation.

Global disease suitability models are essential tools to inform surveillance systems and enable early detection. We present the first global suitability model of highly pathogenic avian influenza (HPAI) H5N1 and demonstrate that reliable predictions can be obtained at global scale. Best predictions are obtained using spatial predictor variables describing host distributions, rather than land use or eco-climatic spatial predictor variables, with a strong association with domestic duck and extensively raised chicken densities. Our results also support a more systematic use of spatial cross-validation in large-scale disease suitability modelling compared to standard random cross-validation that can lead to unreliable measure of extrapolation accuracy. A global suitability model of the H5 clade 2.3.4.4 viruses, a group of viruses that recently spread extensively in Asia and the US, shows in comparison a lower spatial extrapolation capacity than the HPAI H5N1 models, with a stronger association with intensively raised chicken densities and anthropogenic factors.

Movement analysis of free-grazing domestic ducks in Poyang Lake, China: a disease connection.

Previous work suggests domestic poultry are important contributors to the emergence and transmission of highly pathogenic avian influenza throughout Asia. In Poyang Lake, China, domestic duck production cycles are synchronized with arrival and departure of thousands of migratory wild birds in the area. During these periods, high densities of juvenile domestic ducks are in close proximity to migratory wild ducks, increasing the potential for the virus to be transmitted and subsequently disseminated via migration. In this paper, we use GPS dataloggers and dynamic Brownian bridge models to describe movements and habitat use of free-grazing domestic ducks in the Poyang Lake basin and identify specific areas that may have the highest risk of H5N1 transmission between domestic and wild birds. Specifically, we determine relative use by free-grazing domestic ducks of natural wetlands, which are the most heavily used areas by migratory wild ducks, and of rice paddies, which provide habitat for resident wild ducks and lower densities of migratory wild ducks. To our knowledge, this is the first movement study on domestic ducks, and our data show potential for free-grazing domestic ducks from farms located near natural wetlands to come in contact with wild waterfowl, thereby increasing the risk for disease transmission. This study provides an example of the importance of movement ecology studies in understanding dynamics such as disease transmission on a complicated landscape.

Investigating poultry trade patterns to guide avian influenza surveillance and control: a case study in Vietnam.

Live bird markets are often the focus of surveillance activities monitoring avian influenza viruses (AIV) circulating in poultry. However, in order to ensure a high sensitivity of virus detection and effectiveness of management actions, poultry management practices features influencing AIV dynamics need to be accounted for in the design of surveillance programmes. In order to address this knowledge gap, a cross-sectional survey was conducted through interviews with 791 traders in 18 Vietnamese live bird markets. Markets greatly differed according to the sources from which poultry was obtained, and their connections to other markets through the movements of their traders. These features, which could be informed based on indicators that are easy to measure, suggest that markets could be used as sentinels for monitoring virus strains circulating in specific segments of the poultry production sector. AIV spread within markets was modelled. Due to the high turn-over of poultry, viral amplification was likely to be minimal in most of the largest markets. However, due to the large number of birds being introduced each day, and challenges related to cleaning and disinfection, environmental accumulation of viruses at markets may take place, posing a threat to the poultry production sector and to public health.

Clade-level Spatial Modelling of HPAI H5N1 Dynamics in the Mekong Region Reveals New Patterns and Associations with Agro-Ecological Factors.

The highly pathogenic avian influenza (HPAI) H5N1 virus has been circulating in Asia since 2003 and diversified into several genetic lineages, or clades. Although the spatial distribution of its outbreaks was extensively studied, differences in clades were never previously taken into account. We developed models to quantify associations over time and space between different HPAI H5N1 viruses from clade 1, 2.3.4 and 2.3.2 and agro-ecological factors. We found that the distribution of clades in the Mekong region from 2004 to 2013 was strongly regionalised, defining specific epidemiological zones, or epizones. Clade 1 became entrenched in the Mekong Delta and was not supplanted by newer clades, in association with a relatively higher presence of domestic ducks. In contrast, two new clades were introduced (2.3.4 and 2.3.2) in northern Viet Nam and were associated with higher chicken density and more intensive chicken production systems. We suggest that differences in poultry production systems in these different epizones may explain these associations, along with differences in introduction pressure from neighbouring countries. The different distribution patterns found at the clade level would not be otherwise apparent through analysis treating all outbreaks equally, which requires improved linking of disease outbreak records and genetic sequence data.

Spatial Modeling of Wild Bird Risk Factors for Highly Pathogenic A(H5N1) Avian Influenza Virus Transmission.

One of the longest-persisting avian influenza viruses in history, highly pathogenic avian influenza virus (HPAIV) A(H5N1), continues to evolve after 18 yr, advancing the threat of a global pandemic. Wild waterfowl (family Anatidae) are reported as secondary transmitters of HPAIV and primary reservoirs for low-pathogenic avian influenza viruses, yet spatial inputs for disease risk modeling for this group have been lacking. Using geographic information software and Monte Carlo simulations, we developed geospatial indices of waterfowl abundance at 1 and 30 km resolutions and for the breeding and wintering seasons for China, the epicenter of H5N1. Two spatial layers were developed: cumulative waterfowl abundance (WAB), a measure of predicted abundance across species, and cumulative abundance weighted by H5N1 prevalence (WPR), whereby abundance for each species was adjusted based on prevalence values and then totaled across species. Spatial patterns of the model output differed between seasons, with higher WAB and WPR in the northern and western regions of China for the breeding season and in the southeast for the wintering season. Uncertainty measures indicated highest error in southeastern China for both WAB and WPR. We also explored the effect of resampling waterfowl layers from 1 to 30 km resolution for multiscale risk modeling. Results indicated low average difference (less than 0.16 and 0.01 standard deviations for WAB and WPR, respectively), with greatest differences in the north for the breeding season and southeast for the wintering season. This work provides the first geospatial models of waterfowl abundance available for China. The indices provide important inputs for modeling disease transmission risk at the interface of poultry and wild birds. These models are easily adaptable, have broad utility to both disease and conservation needs, and will be available to the scientific community for advanced modeling applications.

Income Disparities and the Global Distribution of Intensively Farmed Chicken and Pigs.

The rapid transformation of the livestock sector in recent decades brought concerns on its impact on greenhouse gas emissions, disruptions to nitrogen and phosphorous cycles and on land use change, particularly deforestation for production of feed crops. Animal and human health are increasingly interlinked through emerging infectious diseases, zoonoses, and antimicrobial resistance. In many developing countries, the rapidity of change has also had social impacts with increased risk of marginalisation of smallholder farmers. However, both the impacts and benefits of livestock farming often differ between extensive (backyard farming mostly for home-consumption) and intensive, commercial production systems (larger herd or flock size, higher investments in inputs, a tendency towards market-orientation). A density of 10,000 chickens per km2 has different environmental, epidemiological and societal implications if these birds are raised by 1,000 individual households or in a single industrial unit. Here, we introduce a novel relationship that links the national proportion of extensively raised animals to the gross domestic product (GDP) per capita (in purchasing power parity). This relationship is modelled and used together with the global distribution of rural population to disaggregate existing 10 km resolution global maps of chicken and pig distributions into extensive and intensive systems. Our results highlight countries and regions where extensive and intensive chicken and pig production systems are most important. We discuss the sources of uncertainties, the modelling assumptions and ways in which this approach could be developed to forecast future trajectories of intensification.

Molecular evidence of Ebola Reston virus infection in Philippine bats.

BACKGROUND: In 2008-09, evidence of Reston ebolavirus (RESTV) infection was found in domestic pigs and pig workers in the Philippines. With species of bats having been shown to be the cryptic reservoir of filoviruses elsewhere, the Philippine government, in conjunction with the Food and Agriculture Organization of the United Nations, assembled a multi-disciplinary and multi-institutional team to investigate Philippine bats as the possible reservoir of RESTV. METHODS: The team undertook surveillance of bat populations at multiple locations during 2010 using both serology and molecular assays. RESULTS: A total of 464 bats from 21 species were sampled. We found both molecular and serologic evidence of RESTV infection in multiple bat species. RNA was detected with quantitative PCR (qPCR) in oropharyngeal swabs taken from Miniopterus schreibersii, with three samples yielding a product on conventional hemi-nested PCR whose sequences differed from a Philippine pig isolate by a single nucleotide. Uncorroborated qPCR detections may indicate RESTV nucleic acid in several additional bat species (M. australis, C. brachyotis and Ch. plicata). We also detected anti-RESTV antibodies in three bats (Acerodon jubatus) using both Western blot and ELISA. CONCLUSIONS: The findings suggest that ebolavirus infection is taxonomically widespread in Philippine bats, but the evident low prevalence and low viral load warrants expanded surveillance to elaborate the findings, and more broadly, to determine the taxonomic and geographic occurrence of ebolaviruses in bats in the region.

Mapping migratory flyways in Asia using dynamic Brownian bridge movement models.

BACKGROUND: Identifying movement routes and stopover sites is necessary for developing effective management and conservation strategies for migratory animals. In the case of migratory birds, a collection of migration routes, known as a flyway, is often hundreds to thousands of kilometers long and can extend across political boundaries. Flyways encompass the entire geographic range between the breeding and non-breeding areas of a population, species, or a group of species, and they provide spatial frameworks for management and conservation across international borders. Existing flyway maps are largely qualitative accounts based on band returns and survey data rather than observed movement routes. In this study, we use satellite and GPS telemetry data and dynamic Brownian bridge movement models to build upon existing maps and describe waterfowl space use probabilistically in the Central Asian and East Asian-Australasian Flyways. RESULTS: Our approach provided new information on migratory routes that was not easily attainable with existing methods to describe flyways. Utilization distributions from dynamic Brownian bridge movement models identified key staging and stopover sites, migration corridors and general flyway outlines in the Central Asian and East Asian-Australasian Flyways. A map of space use from ruddy shelducks depicted two separate movement corridors within the Central Asian Flyway, likely representing two distinct populations that show relatively strong connectivity between breeding and wintering areas. Bar-headed geese marked at seven locations in the Central Asian Flyway showed heaviest use at several stopover sites in the same general region of high-elevation lakes along the eastern Qinghai-Tibetan Plateau. Our analysis of data from multiple Anatidae species marked at sites throughout Asia highlighted major movement corridors across species and confirmed that the Central Asian and East Asian-Australasian Flyways were spatially distinct. CONCLUSIONS: The dynamic Brownian bridge movement model improves our understanding of flyways by estimating relative use of regions in the flyway while providing detailed, quantitative information on migration timing and population connectivity including uncertainty between locations. This model effectively quantifies the relative importance of different migration corridors and stopover sites and may help prioritize specific areas in flyways for conservation of waterbird populations.

Avian influenza H5N1 viral and bird migration networks in Asia.

The spatial spread of the highly pathogenic avian influenza virus H5N1 and its long-term persistence in Asia have resulted in avian influenza panzootics and enormous economic losses in the poultry sector. However, an understanding of the regional long-distance transmission and seasonal patterns of the virus is still lacking. In this study, we present a phylogeographic approach to reconstruct the viral migration network. We show that within each wild fowl migratory flyway, the timing of H5N1 outbreaks and viral migrations are closely associated, but little viral transmission was observed between the flyways. The bird migration network is shown to better reflect the observed viral gene sequence data than other networks and contributes to seasonal H5N1 epidemics in local regions and its large-scale transmission along flyways. These findings have potentially far-reaching consequences, improving our understanding of how bird migration drives the periodic reemergence of H5N1 in Asia.

Stable Isotopes Suggest Low Site Fidelity in Bar-headed Geese (Anser indicus) in Mongolia: Implications for Disease Transmission.

Population connectivity is an important consideration in studies of disease transmission and biological conservation, especially with regard to migratory species. Determining how and when different subpopulations intermingle during different phases of the annual cycle can help identify important geographical regions or features as targets for conservation efforts and can help inform our understanding of continental-scale disease transmission. In this study, stable isotopes of hydrogen and carbon in contour feathers were used to assess the degree of molt-site fidelity among Bar-headed Geese (Anser indicus) captured in north-central Mongolia. Samples were collected from actively molting Bar-headed Geese (n = 61), and some individual samples included both a newly grown feather (still in sheath) and an old, worn feather from the bird's previous molt (n = 21). Although there was no difference in mean hydrogen isotope ratios for the old and new feathers, the isotopic variance in old feathers was approximately three times higher than that of the new feathers, which suggests that these birds use different and geographically distant molting locations from year to year. To further test this conclusion, online data and modeling tools from the isoMAP website were used to generate probability landscapes for the origin of each feather. Likely molting locations were much more widespread for old feathers than for new feathers, which supports the prospect of low molt-site fidelity. This finding indicates that population connectivity would be greater than expected based on data from a single annual cycle, and that disease spread can be rapid even in areas like Mongolia where Bar-headed Geese generally breed in small isolated groups.

Bird Migration and Avian Influenza: A Comparison of Hydrogen Stable Isotopes and Satellite Tracking Methods.

Satellite-based tracking of migratory waterfowl is an important tool for understanding the potential role of wild birds in the long-distance transmission of highly pathogenic avian influenza. However, employing this technique on a continental scale is prohibitively expensive. This study explores the utility of stable isotope ratios in feathers in examining both the distances traveled by migratory birds and variation in migration behavior. We compared the satellite-derived movement data of 22 ducks from 8 species captured at wintering areas in Bangladesh, Turkey, and Hong Kong with deuterium ratios (δD) of these and other individuals captured at the same locations. We derived likely molting locations from the satellite tracking data and generated expected isotope ratios based on an interpolated map of δD in rainwater. Although δD was correlated with the distance between wintering and molting locations, surprisingly, measured δD values were not correlated with either expected values or latitudes of molting sites. However, population-level parameters derived from the satellite-tracking data, such as mean distance between wintering and molting locations and variation in migration distance, were reflected by means and variation of the stable isotope values. Our findings call into question the relevance of the rainfall isotope map for Asia for linking feather isotopes to molting locations, and underscore the need for extensive ground truthing in the form of feather-based isoscapes. Nevertheless, stable isotopes from feathers could inform disease models by characterizing the degree to which regional breeding populations interact at common wintering locations. Feather isotopes also could aid in surveying wintering locations to determine where high-resolution tracking techniques (e.g. satellite tracking) could most effectively be employed. Moreover, intrinsic markers such as stable isotopes offer the only means of inferring movement information from birds that have died as a result of infection. In the absence of feather based-isoscapes, we recommend a combination of isotope analysis and satellite-tracking as the best means of generating aggregate movement data for informing disease models.