Assessing the importance of demographic risk factors across two waves of SARS-CoV-2 using fine-scale case data.

For the long term control of an infectious disease such as COVID-19, it is crucial to identify the most likely individuals to become infected and the role that differences in demographic characteristics play in the observed patterns of infection. As high-volume surveillance winds down, testing data from earlier periods are invaluable for studying risk factors for infection in detail. Observed changes in time during these periods may then inform how stable the pattern will be in the long term. To this end we analyse the distribution of cases of COVID-19 across Scotland in 2021, where the location (census areas of order 500-1,000 residents) and reporting date of cases are known. We consider over 450,000 individually recorded cases, in two infection waves triggered by different lineages: B.1.1.529 ("Omicron") and B.1.617.2 ("Delta"). We use random forests, informed by measures of geography, demography, testing and vaccination. We show that the distributions are only adequately explained when considering multiple explanatory variables, implying that case heterogeneity arose from a combination of individual behaviour, immunity, and testing frequency. Despite differences in virus lineage, time of year, and interventions in place, we find the risk factors remained broadly consistent between the two waves. Many of the observed smaller differences could be reasonably explained by changes in control measures.

Ascertainment rate of SARS-CoV-2 infections from healthcare and community testing in the UK.

The proportion of SARS-CoV-2 infections ascertained through healthcare and community testing is generally unknown and expected to vary depending on natural factors and changes in test-seeking behaviour. Here we use population surveillance data and reported daily case numbers in the United Kingdom to estimate the rate of case ascertainment. We mathematically describe the relationship between the ascertainment rate, the daily number of reported cases, population prevalence, and the sensitivity of PCR and Lateral Flow tests as a function time since exposure. Applying this model to the data, we estimate that 20%-40% of SARS-CoV-2 infections in the UK were ascertained with a positive test with results varying by time and region. Cases of the Alpha variant were ascertained at a higher rate than the wild type variants circulating in the early pandemic, and higher again for the Delta variant and Omicron BA.1 sub-lineage, but lower for the BA.2 sub-lineage. Case ascertainment was higher in adults than in children. We further estimate the daily number of infections and compare this to mortality data to estimate that the infection fatality rate increased by a factor of 3 during the period dominated by the Alpha variant, and declined in line with the distribution of vaccines. This manuscript was submitted as part of a theme issue on "Modelling COVID-19 and Preparedness for Future Pandemics".

A new phylodynamic model of Mycobacterium bovis transmission in a multi-host system uncovers the role of the unobserved reservoir.

Multi-host pathogens are particularly difficult to control, especially when at least one of the hosts acts as a hidden reservoir. Deep sequencing of densely sampled pathogens has the potential to transform this understanding, but requires analytical approaches that jointly consider epidemiological and genetic data to best address this problem. While there has been considerable success in analyses of single species systems, the hidden reservoir problem is relatively under-studied. A well-known exemplar of this problem is bovine Tuberculosis, a disease found in British and Irish cattle caused by Mycobacterium bovis, where the Eurasian badger has long been believed to act as a reservoir but remains of poorly quantified importance except in very specific locations. As a result, the effort that should be directed at controlling disease in badgers is unclear. Here, we analyse densely collected epidemiological and genetic data from a cattle population but do not explicitly consider any data from badgers. We use a simulation modelling approach to show that, in our system, a model that exploits available cattle demographic and herd-to-herd movement data, but only considers the ability of a hidden reservoir to generate pathogen diversity, can be used to choose between different epidemiological scenarios. In our analysis, a model where the reservoir does not generate any diversity but contributes to new infections at a local farm scale are significantly preferred over models which generate diversity and/or spread disease at broader spatial scales. While we cannot directly attribute the role of the reservoir to badgers based on this analysis alone, the result supports the hypothesis that under current cattle control regimes, infected cattle alone cannot sustain M. bovis circulation. Given the observed close phylogenetic relationship for the bacteria taken from cattle and badgers sampled near to each other, the most parsimonious hypothesis is that the reservoir is the infected badger population. More broadly, our approach demonstrates that carefully constructed bespoke models can exploit the combination of genetic and epidemiological data to overcome issues of extreme data bias, and uncover important general characteristics of transmission in multi-host pathogen systems.

Site Specific Relationships between COVID-19 Cases and SARS-CoV-2 Viral Load in Wastewater Treatment Plant Influent.

Wastewater based epidemiology (WBE) has become an important tool during the COVID-19 pandemic, however the relationship between SARS-CoV-2 RNA in wastewater treatment plant influent (WWTP) and cases in the community is not well-defined. We report here the development of a national WBE program across 28 WWTPs serving 50% of the population of Scotland, including large conurbations, as well as low-density rural and remote island communities. For each WWTP catchment area, we quantified spatial and temporal relationships between SARS-CoV-2 RNA in wastewater and COVID-19 cases. Daily WWTP SARS-CoV-2 influent viral RNA load, calculated using daily influent flow rates, had the strongest correlation (ρ > 0.9) with COVID-19 cases within a catchment. As the incidence of COVID-19 cases within a community increased, a linear relationship emerged between cases and influent viral RNA load. There were significant differences between WWTPs in their capacity to predict case numbers based on influent viral RNA load, with the limit of detection ranging from 25 cases for larger plants to a single case in smaller plants. SARS-CoV-2 viral RNA load can be used to predict the number of cases detected in the WWTP catchment area, with a clear statistically significant relationship observed above site-specific case thresholds.

Predictive analysis across spatial scales links zoonotic malaria to deforestation.

The complex transmission ecologies of vector-borne and zoonotic diseases pose challenges to their control, especially in changing landscapes. Human incidence of zoonotic malaria ( Plasmodium knowlesi) is associated with deforestation although mechanisms are unknown. Here, a novel application of a method for predicting disease occurrence that combines machine learning and statistics is used to identify the key spatial scales that define the relationship between zoonotic malaria cases and environmental change. Using data from satellite imagery, a case-control study, and a cross-sectional survey, predictive models of household-level occurrence of P. knowlesi were fitted with 16 variables summarized at 11 spatial scales simultaneously. The method identified a strong and well-defined peak of predictive influence of the proportion of cleared land within 1 km of households on P. knowlesi occurrence. Aspect (1 and 2 km), slope (0.5 km) and canopy regrowth (0.5 km) were important at small scales. By contrast, fragmentation of deforested areas influenced P. knowlesi occurrence probability most strongly at large scales (4 and 5 km). The identification of these spatial scales narrows the field of plausible mechanisms that connect land use change and P. knowlesi, allowing for the refinement of disease occurrence predictions and the design of spatially-targeted interventions.

Disease management at the wildlife-livestock interface: Using whole-genome sequencing to study the role of elk in Mycobacterium bovis transmission in Michigan, USA.

The role of wildlife in the persistence and spread of livestock diseases is difficult to quantify and control. These difficulties are exacerbated when several wildlife species are potentially involved. Bovine tuberculosis (bTB), caused by Mycobacterium bovis, has experienced an ecological shift in Michigan, with spillover from cattle leading to an endemically infected white-tailed deer (deer) population. It has potentially substantial implications for the health and well-being of both wildlife and livestock and incurs a significant economic cost to industry and government. Deer are known to act as a reservoir of infection, with evidence of M. bovis transmission to sympatric elk and cattle populations. However, the role of elk in the circulation of M. bovis is uncertain; they are few in number, but range further than deer, so may enable long distance spread. Combining Whole Genome Sequences (WGS) for M. bovis isolates from exceptionally well-observed populations of elk, deer and cattle with spatiotemporal locations, we use spatial and Bayesian phylogenetic analyses to show strong spatiotemporal admixture of M. bovis isolates. Clustering of bTB in elk and cattle suggests either intraspecies transmission within the two populations, or exposure to a common source. However, there is no support for significant pathogen transfer amongst elk and cattle, and our data are in accordance with existing evidence that interspecies transmission in Michigan is likely only maintained by deer. This study demonstrates the value of whole genome population studies of M. bovis transmission at the wildlife-livestock interface, providing insights into bTB management in an endemic system.

Quantifying the roles of host movement and vector dispersal in the transmission of vector-borne diseases of livestock.

The role of host movement in the spread of vector-borne diseases of livestock has been little studied. Here we develop a mathematical framework that allows us to disentangle and quantify the roles of vector dispersal and livestock movement in transmission between farms. We apply this framework to outbreaks of bluetongue virus (BTV) and Schmallenberg virus (SBV) in Great Britain, both of which are spread by Culicoides biting midges and have recently emerged in northern Europe. For BTV we estimate parameters by fitting the model to outbreak data using approximate Bayesian computation, while for SBV we use previously derived estimates. We find that around 90% of transmission of BTV between farms is a result of vector dispersal, while for SBV this proportion is 98%. This difference is a consequence of higher vector competence and shorter duration of viraemia for SBV compared with BTV. For both viruses we estimate that the mean number of secondary infections per infected farm is greater than one for vector dispersal, but below one for livestock movements. Although livestock movements account for a small proportion of transmission and cannot sustain an outbreak on their own, they play an important role in establishing new foci of infection. However, the impact of restricting livestock movements on the spread of both viruses depends critically on assumptions made about the distances over which vector dispersal occurs. If vector dispersal occurs primarily at a local scale (99% of transmission occurs <25 km), movement restrictions are predicted to be effective at reducing spread, but if dispersal occurs frequently over longer distances (99% of transmission occurs <50 km) they are not.

Plasmodium knowlesi invasion following spread by infected mosquitoes, macaques and humans.

Plasmodium knowlesi is increasingly recognized as a major cause of malaria in Southeast Asia. Anopheles leucosphyrous group mosquitoes transmit the parasite and natural hosts include long-tailed and pig-tailed macaques. Despite early laboratory experiments demonstrating successful passage of infection between humans, the true role that humans play in P. knowlesi epidemiology remains unclear. The threat posed by its introduction into immunologically naïve populations is unknown despite being a public health priority for this region. A two-host species mathematical model was constructed to analyse this threat. Global sensitivity analysis using Monte Carlo methods highlighted the biological processes of greatest influence to transmission. These included parameters known to be influential in classic mosquito-borne disease models (e.g. vector longevity); however, interesting ecological components that are specific to this system were also highlighted: while local vectors likely have intrinsic preferences for certain host species, how plastic these preferences are, and how this is shaped by local conditions, are key determinants of parasite transmission potential. Invasion analysis demonstrates that this behavioural plasticity can qualitatively impact the probability of an epidemic sparked by imported infection. Identifying key vector sub/species and studying their biting behaviours constitute important next steps before models can better assist in strategizing disease control.

Using whole genome sequencing to investigate transmission in a multi-host system: bovine tuberculosis in New Zealand.

BACKGROUND: Bovine tuberculosis (bTB), caused by Mycobacterium bovis, is an important livestock disease raising public health and economic concerns around the world. In New Zealand, a number of wildlife species are implicated in the spread and persistence of bTB in cattle populations, most notably the brushtail possum (Trichosurus vulpecula). Whole Genome Sequenced (WGS) M. bovis isolates sourced from infected cattle and wildlife across New Zealand were analysed. Bayesian phylogenetic analyses were conducted to estimate the substitution rate of the sampled population and investigate the role of wildlife. In addition, the utility of WGS was examined with a view to these methods being incorporated into routine bTB surveillance. RESULTS: A high rate of exchange was evident between the sampled wildlife and cattle populations but directional estimates of inter-species transmission were sensitive to the sampling strategy employed. A relatively high substitution rate was estimated, this, in combination with a strong spatial signature and a good agreement to previous typing methods, acts to endorse WGS as a typing tool. CONCLUSIONS: In agreement with the current knowledge of bTB in New Zealand, transmission of M. bovis between cattle and wildlife was evident. Without direction, these estimates are less informative but taken in conjunction with the low prevalence of bTB in New Zealand's cattle population it is likely that, currently, wildlife populations are acting as the main bTB reservoir. Wildlife should therefore continue to be targeted if bTB is to be eradicated from New Zealand. WGS will be a considerable aid to bTB eradication by greatly improving the discriminatory power of molecular typing data. The substitution rates estimated here will be an important part of epidemiological investigations using WGS data.

Relative abundance of Mycobacterium bovis molecular types in cattle: a simulation study of potential epidemiological drivers.

BACKGROUND: The patterns of relative species abundance are commonly studied in ecology and epidemiology to provide insights into underlying dynamical processes. Molecular types (MVLA-types) of Mycobacterium bovis, the causal agent of bovine tuberculosis, are now routinely recorded in culture-confirmed bovine tuberculosis cases in Northern Ireland. In this study, we use ecological approaches and simulation modelling to investigate the distribution of relative abundances of MVLA-types and its potential drivers. We explore four biologically plausible hypotheses regarding the processes driving molecular type relative abundances: sampling and speciation; structuring of the pathogen population; historical changes in population size; and transmission heterogeneity (superspreading). RESULTS: Northern Irish herd-level MVLA-type surveillance shows a right-skewed distribution of MVLA-types, with a small number of types present at very high frequencies and the majority of types very rare. We demonstrate that this skew is too extreme to be accounted for by simple neutral ecological processes. Simulation results indicate that the process of MVLA-type speciation and the manner in which the MVLA-typing loci were chosen in Northern Ireland cannot account for the observed skew. Similarly, we find that pathogen population structure, assuming for example a reservoir of infection in a separate host, would drive the relative abundance distribution in the opposite direction to that observed, generating more even abundances of molecular types. However, we find that historical increases in bovine tuberculosis prevalence and/or transmission heterogeneity (superspreading) are both capable of generating the skewed MVLA-type distribution, consistent with findings of previous work examining the distribution of molecular types in human tuberculosis. CONCLUSION: Although the distribution of MVLA-type abundances does not fit classical neutral predictions, our simulations show that increases in pathogen population size and/or superspreading are consistent with the pattern observed, even in the absence of selective pressures acting on the system.

Broadwick: a framework for computational epidemiology.

BACKGROUND: Modelling disease outbreaks often involves integrating the wealth of data that are gathered during modern outbreaks into complex mathematical or computational models of transmission. Incorporating these data into simple compartmental epidemiological models is often challenging, requiring the use of more complex but also more efficient computational models. In this paper we introduce a new framework that allows for a more systematic and user-friendly way of building and running epidemiological models that efficiently handles disease data and reduces much of the boilerplate code that usually associated to these models. We introduce the framework by developing an SIR model on a simple network as an example. RESULTS: We develop Broadwick, a modular, object-oriented epidemiological framework that efficiently handles large epidemiological datasets and provides packages for stochastic simulations, parameter inference using Approximate Bayesian Computation (ABC) and Markov Chain Monte Carlo (MCMC) methods. Each algorithm used is fully customisable with sensible defaults that are easily overridden by custom algorithms as required. CONCLUSION: Broadwick is an epidemiological modelling framework developed to increase the productivity of researchers by providing a common framework with which to develop and share complex models. It will appeal to research team leaders as it allows for models to be created prior to a disease outbreak and has the ability to handle large datasets commonly found in epidemiological modelling.

A descriptive analysis of the growth of unrecorded interactions amongst cattle-raising premises in Scotland and their implications for disease spread.

BACKGROUND: Individual animal-level reporting of cattle movements between agricultural holdings is in place in Scotland, and the resulting detailed movement data are used to inform epidemiological models and intervention. However, recent years have seen a rapid increase in the use of registered links that allow Scottish farmers to move cattle between linked holdings without reporting. RESULTS: By analyzing these registered trade links as a number of different networks, we find that the geographical reach of these registered links has increased over time, with many holdings linked indirectly to a large number of holdings, some potentially geographically distant. This increase was not linked to decreases in recorded movements at the holding level. When combining registered links with reported movements, we find that registered links increase the size of a possible outward chain of infection from a Scottish holding. The impact on the maximum size is considerably greater than the impact on the mean. CONCLUSIONS: We outline the magnitude and geographic extent of that increase, and show that this growth both has the potential to substantially increase the size of epidemics driven by livestock movements, and undermines the extensive, invaluable recording within the cattle tracing system in Scotland and, by extension, the rest of Great Britain.

Limitations of variable number of tandem repeat typing identified through whole genome sequencing of Mycobacterium avium subsp. paratuberculosis on a national and herd level.

BACKGROUND: Mycobacterium avium subsp. paratuberculosis (MAP), the causative bacterium of Johne's disease in dairy cattle, is widespread in the Canadian dairy industry and has significant economic and animal welfare implications. An understanding of the population dynamics of MAP can be used to identify introduction events, improve control efforts and target transmission pathways, although this requires an adequate understanding of MAP diversity and distribution between herds and across the country. Whole genome sequencing (WGS) offers a detailed assessment of the SNP-level diversity and genetic relationship of isolates, whereas several molecular typing techniques used to investigate the molecular epidemiology of MAP, such as variable number of tandem repeat (VNTR) typing, target relatively unstable repetitive elements in the genome that may be too unpredictable to draw accurate conclusions. The objective of this study was to evaluate the diversity of bovine MAP isolates in Canadian dairy herds using WGS and then determine if VNTR typing can distinguish truly related and unrelated isolates. RESULTS: Phylogenetic analysis based on 3,039 SNPs identified through WGS of 124 MAP isolates identified eight genetically distinct subtypes in dairy herds from seven Canadian provinces, with the dominant type including over 80% of MAP isolates. VNTR typing of 527 MAP isolates identified 12 types, including "bison type" isolates, from seven different herds. At a national level, MAP isolates differed from each other by 1-2 to 239-240 SNPs, regardless of whether they belonged to the same or different VNTR types. A herd-level analysis of MAP isolates demonstrated that VNTR typing may both over-estimate and under-estimate the relatedness of MAP isolates found within a single herd. CONCLUSIONS: The presence of multiple MAP subtypes in Canada suggests multiple introductions into the country including what has now become one dominant type, an important finding for Johne's disease control. VNTR typing often failed to identify closely and distantly related isolates, limiting the applicability of using this typing scheme to study the molecular epidemiology of MAP at a national and herd-level.

Whole genome sequencing reveals local transmission patterns of Mycobacterium bovis in sympatric cattle and badger populations.

Whole genome sequencing (WGS) technology holds great promise as a tool for the forensic epidemiology of bacterial pathogens. It is likely to be particularly useful for studying the transmission dynamics of an observed epidemic involving a largely unsampled 'reservoir' host, as for bovine tuberculosis (bTB) in British and Irish cattle and badgers. BTB is caused by Mycobacterium bovis, a member of the M. tuberculosis complex that also includes the aetiological agent for human TB. In this study, we identified a spatio-temporally linked group of 26 cattle and 4 badgers infected with the same Variable Number Tandem Repeat (VNTR) type of M. bovis. Single-nucleotide polymorphisms (SNPs) between sequences identified differences that were consistent with bacterial lineages being persistent on or near farms for several years, despite multiple clear whole herd tests in the interim. Comparing WGS data to mathematical models showed good correlations between genetic divergence and spatial distance, but poor correspondence to the network of cattle movements or within-herd contacts. Badger isolates showed between zero and four SNP differences from the nearest cattle isolate, providing evidence for recent transmissions between the two hosts. This is the first direct genetic evidence of M. bovis persistence on farms over multiple outbreaks with a continued, ongoing interaction with local badgers. However, despite unprecedented resolution, directionality of transmission cannot be inferred at this stage. Despite the often notoriously long timescales between time of infection and time of sampling for TB, our results suggest that WGS data alone can provide insights into TB epidemiology even where detailed contact data are not available, and that more extensive sampling and analysis will allow for quantification of the extent and direction of transmission between cattle and badgers.

A restatement of the natural science evidence base relevant to the control of bovine tuberculosis in Great Britain.

Bovine tuberculosis (bTB) is a very important disease of cattle in Great Britain, where it has been increasing in incidence and geographical distribution. In addition to cattle, it infects other species of domestic and wild animals, in particular the European badger (Meles meles). Policy to control bTB is vigorously debated and contentious because of its implications for the livestock industry and because some policy options involve culling badgers, the most important wildlife reservoir. This paper describes a project to provide a succinct summary of the natural science evidence base relevant to the control of bTB, couched in terms that are as policy-neutral as possible. Each evidence statement is placed into one of four categories describing the nature of the underlying information. The evidence summary forms the appendix to this paper and an annotated bibliography is provided in the electronic supplementary material.

An online survey of horse-owners in Great Britain.

BACKGROUND: Contingency planning for potential equine infectious disease outbreaks relies on accurate information on horse location and movements to estimate the risk of dissemination of disease(s). An online questionnaire was used to obtain unique information linking owner and horse location to characteristics of horse movements within and outwith Great Britain (GB). RESULTS: This online survey yielded a strong response, providing more than four times the target number of respondents (1000 target respondents) living in all parts of GB. Key demographic findings of this study indicated that horses which were kept on livery yards and riding schools were likely to be found in urban environments, some distance away from the owner's home and vaccinated against influenza and herpes virus. Survey respondents were likely to travel greater than 10 miles to attend activities such as eventing or endurance but were also likely to travel and return home within a single day (58.6%, 2063/3522). This may affect the geographical extent and speed of disease spread, if large numbers of people from disparate parts of the country are attending the same event and the disease agent is highly infectious or virulent. The greatest risk for disease introduction and spread may be represented by a small proportion of people who import or travel internationally with their horses. These respondents were likely to have foreign horse passports, which were not necessarily recorded in the National Equine Database (NED), making the location of these horses untraceable. CONCLUSIONS: These results illustrate the difficulties which exist with national GB horse traceability despite the existence of the NED and the horse passport system. This study also demonstrates that an online approach could be adopted to obtain important demographic data on GB horse owners on a more routine and frequent basis to inform decisions or policy pertaining to equine disease control. This represents a reasonable alternative to collection of GB horse location and movement data given that the NED no longer exists and there is no immediate plan to replace it.

Empirical distributions of time intervals between COVID-19 cases and more severe outcomes in Scotland.

A critical factor in infectious disease control is the risk of an outbreak overwhelming local healthcare capacity. The overall demand on healthcare services will depend on disease severity, but the precise timing and size of peak demand also depends on the time interval (or clinical time delay) between initial infection, and development of severe disease. A broader distribution of intervals may draw that demand out over a longer period, but have a lower peak demand. These interval distributions are therefore important in modelling trajectories of e.g. hospital admissions, given a trajectory of incidence. Conversely, as testing rates decline, an incidence trajectory may need to be inferred through the delayed, but relatively unbiased signal of hospital admissions. Healthcare demand has been extensively modelled during the COVID-19 pandemic, where localised waves of infection have imposed severe stresses on healthcare services. While the initial acute threat posed by this disease has since subsided with immunity buildup from vaccination and prior infection, prevalence remains high and waning immunity may lead to substantial pressures for years to come. In this work, then, we present a set of interval distributions, for COVID-19 cases and subsequent severe outcomes; hospital admission, ICU admission, and death. These may be used to model more realistic scenarios of hospital admissions and occupancy, given a trajectory of infections or cases. We present a method for obtaining empirical distributions using COVID-19 outcomes data from Scotland between September 2020 and January 2022 (N = 31724 hospital admissions, N = 3514 ICU admissions, N = 8306 mortalities). We present separate distributions for individual age, sex, and deprivation of residing community. While the risk of severe disease following COVID-19 infection is substantially higher for the elderly and those residing in areas of high deprivation, the length of stay shows no strong dependence, suggesting that severe outcomes are equally severe across risk groups. As Scotland and other countries move into a phase where testing is no longer abundant, these intervals may be of use for retrospective modelling of patterns of infection, given data on severe outcomes.

Modeling the effectiveness of targeting Rift Valley fever virus vaccination using imperfect network information.

Livestock movements contribute to the spread of several infectious diseases. Data on livestock movements can therefore be harnessed to guide policy on targeted interventions for controlling infectious livestock diseases, including Rift Valley fever (RVF)-a vaccine-preventable arboviral fever. Detailed livestock movement data are known to be useful for targeting control efforts including vaccination. These data are available in many countries, however, such data are generally lacking in others, including many in East Africa, where multiple RVF outbreaks have been reported in recent years. Available movement data are imperfect, and the impact of this uncertainty in the utility of movement data on informing targeting of vaccination is not fully understood. Here, we used a network simulation model to describe the spread of RVF within and between 398 wards in northern Tanzania connected by cattle movements, on which we evaluated the impact of targeting vaccination using imperfect movement data. We show that pre-emptive vaccination guided by only market movement permit data could prevent large outbreaks. Targeted control (either by the risk of RVF introduction or onward transmission) at any level of imperfect movement information is preferred over random vaccination, and any improvement in information reliability is advantageous to their effectiveness. Our modeling approach demonstrates how targeted interventions can be effectively used to inform animal and public health policies for disease control planning. This is particularly valuable in settings where detailed data on livestock movements are either unavailable or imperfect due to resource limitations in data collection, as well as challenges associated with poor compliance.

Genomic epidemiology of Mycobacterium bovis infection in sympatric badger and cattle populations in Northern Ireland.

Bovine tuberculosis (bTB) is a costly, epidemiologically complex, multi-host, endemic disease. Lack of understanding of transmission dynamics may undermine eradication efforts. Pathogen whole-genome sequencing improves epidemiological inferences, providing a means to determine the relative importance of inter- and intra-species host transmission for disease persistence. We sequenced an exceptional data set of 619 Mycobacterium bovis isolates from badgers and cattle in a 100 km2 bTB 'hotspot' in Northern Ireland. Historical molecular subtyping data permitted the targeting of an endemic pathogen lineage, whose long-term persistence provided a unique opportunity to study disease transmission dynamics in unparalleled detail. Additionally, to assess whether badger population genetic structure was associated with the spatial distribution of pathogen genetic diversity, we microsatellite genotyped hair samples from 769 badgers trapped in this area. Birth death models and TransPhylo analyses indicated that cattle were likely driving the local epidemic, with transmission from cattle to badgers being more common than badger to cattle. Furthermore, the presence of significant badger population genetic structure in the landscape was not associated with the spatial distribution of M. bovis genetic diversity, suggesting that badger-to-badger transmission is not playing a major role in transmission dynamics. Our data were consistent with badgers playing a smaller role in transmission of M. bovis infection in this study site, compared to cattle. We hypothesize, however, that this minor role may still be important for persistence. Comparison to other areas suggests that M. bovis transmission dynamics are likely to be context dependent, with the role of wildlife being difficult to generalize.

Summary of current knowledge of the size and spatial distribution of the horse population within Great Britain.

BACKGROUND: Robust demographic information is important to understanding the risk of introduction and spread of exotic diseases as well as the development of effective disease control strategies, but is often based on datasets collected for other purposes. Thus, it is important to validate, or at least cross-reference these datasets to other sources to assess whether they are being used appropriately. The aim of this study was to use horse location data collected from different contributing industry sectors ("Stakeholder horse data") to calibrate the spatial distribution of horses as indicated by owner locations registered in the National Equine Database (the NED). RESULTS: A conservative estimate for the accurately geo-located NED horse population within GB is approximately 840,000 horses. This is likely to be an underestimate because of the exclusion of horses due to age or location criteria. In both datasets, horse density was higher in England and Wales than in Scotland. The high density of horses located in urban areas as indicated in the NED is consistent with previous reports indicating that owner location cannot always be viewed as a direct substitute for horse location. Otherwise, at a regional resolution, there are few differences between the datasets. There are inevitable biases in the stakeholder data, and leisure horses that are unaffiliated to major stakeholders are not included in these data. Despite this, the similarity in distributions of these datasets is re-assuring, suggesting that there are few regional biases in the NED. CONCLUSIONS: Our analyses suggest that stakeholder data could be used to monitor possible changes in horse demographics. Given such changes in horse demographics and the advantages of stakeholder data (which include annual updates and accurate horse location), it may be appropriate to use these data for future disease modelling in conjunction with, if not in place of the NED.