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.

Evidence for local and international spread of Mycobacterium avium subspecies paratuberculosis through whole genome sequencing of isolates from the island of Ireland.

We describe application of whole genome sequencing (WGS) to a collection of 197 Mycobacterium avium subsp paratuberculosis (MAP) isolates gathered from 122 cattle herds across 27 counties of the island of Ireland. We compare WGS to MAP diversity quantified using mycobacterial interspersed random unit - variable number tandem repeats (MIRU-VNTR). While MIRU-VNTR showed only two major types, WGS could split the 197 isolates into eight major groups. We also found six isolates corresponding to INMV 13, a novel MIRU-VNTR type for Ireland. Evidence for dispersal of MAP across Ireland via cattle movement could be discerned from the data, with mixed infections present in several herds. Furthermore, comparisons of MAP WGS data from Ireland to data from Great Britain and continental Europe revealed many instances of close genetic similarity and hence evidence for international transmission of infection. BEAST MASCOT structured coalescent analyses, with relaxed and strict molecular clocks, estimated the substitution rate to be 0.10-0.13 SNPs/site/year and disclosed greater transitions per lineage per year from Europe to Ireland, indicating transmission into Ireland. Our work therefore reveals new insight into the seeding of MAP infection across Ireland, highlighting how WGS can inform policy formulation to ultimately control MAP transmission at local, national and international scales.

Identifying likely transmissions in Mycobacterium bovis infected populations of cattle and badgers using the Kolmogorov Forward Equations.

Established methods for whole-genome-sequencing (WGS) technology allow for the detection of single-nucleotide polymorphisms (SNPs) in the pathogen genomes sourced from host samples. The information obtained can be used to track the pathogen's evolution in time and potentially identify 'who-infected-whom' with unprecedented accuracy. Successful methods include 'phylodynamic approaches' that integrate evolutionary and epidemiological data. However, they are typically computationally intensive, require extensive data, and are best applied when there is a strong molecular clock signal and substantial pathogen diversity. To determine how much transmission information can be inferred when pathogen genetic diversity is low and metadata limited, we propose an analytical approach that combines pathogen WGS data and sampling times from infected hosts. It accounts for 'between-scale' processes, in particular within-host pathogen evolution and between-host transmission. We applied this to a well-characterised population with an endemic Mycobacterium bovis (the causative agent of bovine/zoonotic tuberculosis, bTB) infection. Our results show that, even with such limited data and low diversity, the computation of the transmission probability between host pairs can help discriminate between likely and unlikely infection pathways and therefore help to identify potential transmission networks. However, the method can be sensitive to assumptions about within-host evolution.

Updated functional annotation of the Mycobacterium bovis AF2122/97 reference genome.

Mycobacterium bovis AF2122/97 is the reference strain for the bovine tuberculosis bacillus. Here we report an update to the M. bovis AF2122/97 genome annotation to reflect 616 new protein identifications that replace many of the old hypothetical coding sequences and proteins of unknown function in the genome. These changes integrate information from functional assignments of orthologous coding sequences in the Mycobacterium tuberculosis H37Rv genome. We have also added 69 additional new gene names.

Mycobacterium bovis genomics reveals transmission of infection between cattle and deer in Ireland.

Control of bovine tuberculosis (bTB), caused by Mycobacterium bovis, in the Republic of Ireland costs €84 million each year. Badgers are recognized as being a wildlife source for M. bovis infection of cattle. Deer are thought to act as spillover hosts for infection; however, population density is recognized as an important driver in shifting their epidemiological role, and deer populations across the country have been increasing in density and range. County Wicklow represents one specific area in the Republic of Ireland with a high density of deer that has had consistently high bTB prevalence for over a decade, despite control operations in both cattle and badgers. Our research used whole-genome sequencing of M. bovis sourced from infected cattle, deer and badgers in County Wicklow to evaluate whether the epidemiological role of deer could have shifted from spillover host to source. Our analyses reveal that cattle and deer share highly similar M. bovis strains, suggesting that transmission between these species is occurring in the area. In addition, the high level of diversity observed in the sampled deer population suggests deer may be acting as a source of infection for local cattle populations. These findings have important implications for the control and ultimate eradication of bTB in Ireland.

Preliminary evaluation of a novel, fully automated, Telenostic device for rapid field-diagnosis of cattle parasites.

New ideas for diagnostics in clinical parasitology are needed to overcome some of the difficulties experienced in the widespread adoption of detection methods for gastrointestinal parasites in livestock. Here we provide an initial evaluation of the performance of a newly developed automated device (Telenostic) to identify and quantify parasitic elements in fecal samples. This study compared the Telenostic device with the McMaster and Mini-FLOTAC for counting of strongyle eggs in a fecal sample. Three bovine fecal samples were examined, in triplicate, on each of the three fecal egg-counting devices. In addition, both manual (laboratory technician) and automated analysis (image analysis algorithm) were performed on the Telenostic device to calculate fecal egg counts (FEC). Overall, there were consistent egg counts reported across the three devices and calculation methods. The Telenostic device compared very favourably to the Mini-FLOTAC and McMaster. Only in sample C, a significant difference (P < 0.05) was observed between the egg counts obtained by Mini-FLOTAC and by the other methods. From this limited dataset it can be concluded that the Telenostic-automated test is comparable to currently used benchmark FEC methods, while improving the workflow, test turn-around time and not requiring trained laboratory personnel to operate or interpret the results.

Combining genomics and epidemiology to analyse bi-directional transmission of Mycobacterium bovis in a multi-host system.

Quantifying pathogen transmission in multi-host systems is difficult, as exemplified in bovine tuberculosis (bTB) systems, but is crucial for control. The agent of bTB, Mycobacterium bovis, persists in cattle populations worldwide, often where potential wildlife reservoirs exist. However, the relative contribution of different host species to bTB persistence is generally unknown. In Britain, the role of badgers in infection persistence in cattle is highly contentious, despite decades of research and control efforts. We applied Bayesian phylogenetic and machine-learning approaches to bacterial genome data to quantify the roles of badgers and cattle in M. bovis infection dynamics in the presence of data biases. Our results suggest that transmission occurs more frequently from badgers to cattle than vice versa (10.4x in the most likely model) and that within-species transmission occurs at higher rates than between-species transmission for both. If representative, our results suggest that control operations should target both cattle and badgers.

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.

HomoplasyFinder: a simple tool to identify homoplasies on a phylogeny.

A homoplasy is a nucleotide identity resulting from a process other than inheritance from a common ancestor. Importantly, by distorting the ancestral relationships between nucleotide sequences, homoplasies can change the structure of the phylogeny. Homoplasies can emerge naturally, especially under high selection pressures and/or high mutation rates, or be created during the generation and processing of sequencing data. Identification of homoplasies is critical, both to understand their influence on the analyses of phylogenetic data and to allow an investigation into how they arose. Here we present HomoplasyFinder, a java application that can be used as a stand-a-lone tool or within the statistical programming environment R. Within R and Java, HomoplasyFinder is shown to be able to automatically, and quickly, identify any homoplasies present in simulated and real phylogenetic data. HomoplasyFinder can easily be incorporated into existing analysis pipelines, either within or outside of R, allowing the user to quickly identify homoplasies to inform downstream analyses and interpretation.

Whole Genome Sequencing for Determining the Source of Mycobacterium bovis Infections in Livestock Herds and Wildlife in New Zealand.

The ability to DNA fingerprint Mycobacterium bovis isolates helped to define the role of wildlife in the persistence of bovine tuberculosis in New Zealand. DNA fingerprinting results currently help to guide wildlife control measures and also aid in tracing the source of infections that result from movement of livestock. During the last 5 years we have developed the ability to distinguish New Zealand (NZ) M. bovis isolates by comparing the sequences of whole genome sequenced (WGS) M. bovis samples. WGS provides much higher resolution than our other established typing methods and greatly improves the definition of the regional localization of NZ M. bovis types. Three outbreak investigations are described and results demonstrate how WGS analysis has led to the confirmation of epidemiological sourcing of infection, to better definition of new sources of infection by ruling out other possible sources, and has revealed probable wildlife infection in an area considered to be free of infected wildlife. The routine use of WGS analyses for sourcing new M. bovis infections will be an important component of the strategy employed to eradicate bovine TB from NZ livestock and wildlife.

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.