Longitudinal dynamics of herd-level Mycobacterium bovis MLVA type surveillance in cattle in Northern Ireland 2003-2016

Whole-Genome sequencing in routine Mycobacterium bovis epidemiology - scoping the potential

Mycobacterium bovis the main agent of bovine tuberculosis (bTB), presents as a series of spatially-localised micro-epidemics across landscapes. Classical molecular typing methods applied to these micro-epidemics, based on genotyping a few variable loci, have significantly improved our understanding of potential epidemiological links between outbreaks. However, they have limited utility owing to low resolution. Conversely, whole-genome sequencing (WGS) provides the highest resolution data available for molecular epidemiology, producing richer outbreak tracing, insights into phylogeography and epidemic evolutionary history. We illustrate these advantages by focusing on a common single lineage of M. bovis (1.140) from Northern Ireland. Specifically, we investigate the spatial sub-structure of 20 years of herd-level multi locus VNTR analysis (MLVA) surveillance data and WGS data from a down sampled subset of isolates of this MLVA type over the same time frame. We mapped 2108 isolate locations of MLVA type 1.140 over the years 2000-2022. We also mapped the locations of 148 contemporary WGS isolates from this lineage, over a similar geographic range, stratifying by single nucleotide polymorphism (SNP) relatedness cut-offs of 15 SNPs. We determined a putative core range for the 1.140 MLVA type and SNP-defined sequence clusters using a 50 % kernel density estimate, using cattle movement data to inform on likely sources of WGS isolates found outside of core ranges. Finally, we applied Bayesian phylogenetic methods to investigate past population history and reproductive number of the 1.140 M. bovis lineage. We demonstrate that WGS SNP-defined clusters exhibit smaller core ranges than the established MLVA type - facilitating superior disease tracing. We also demonstrate the superior functionality of WGS data in determining how this lineage was disseminated across the landscape, likely via cattle movement and to infer how its effective population size and reproductive number has been in flux since its emergence. These initial findings highlight the potential of WGS data for routine monitoring of bTB outbreaks.

Geo-epidemiology of animal tuberculosis and Mycobacterium bovis genotypes in livestock in a small, high-incidence area in Sicily, Italy

Introduction

The persistence of animal tuberculosis (TB) in livestock is a major concern in Sicily, Italy. The objective of this study was to elucidate the transmission dynamics of M. bovis infection in a highly circumscribed, and at the same time geographically diverse, high-risk area of the island through an in-depth geo-epidemiological investigation of TB in cattle and black pigs raised in small-scale extensive farms across the district of Caronia.

Methods

We used genotype analysis coupled with geographic information system (GIS) technology and phylogenetic inference to characterize the spatial distribution of TB and M. bovis genotypes in livestock and the genetic relationships between M. bovis isolates. A total of 589 M. bovis isolates collected from slaughtered cattle (n = 527) and Sicilian black pigs (n = 62) over a 5-year period (2014–2018) were included in the study.

Results

TB was widespread throughout the district and was most frequent in the north-central area of the district, especially along one of the district’s streams. We identified a total of 62 M. bovis genotypes. Identical genetic profiles were isolated from both neighboring and non-neighburing herds. The 10 most frequent genotypes, accounting for 82% of M. bovis isolates, showed geographic specificities in that they tended to cluster in specific spatial niches. The landscape structure of these niches—i.e. steep slopes, rocky ridges, meadows and streams—is likely to have had a significant influence on the distribution of TB among livestock in Caronia. Higher concentrations of TB were observed along streams and in open meadows, while rocky ridges and slopes appeared to have hampered the spread of TB.

Discussion

The geographical distribution of TB cases among livestock in Caronia is consistent with several epidemiological scenarios (e.g., high density of infected herds along the streams or in hilly plateau where livestock share pastures). Landscape structure is likely to play an important role in the transmission and persistence of M. bovis infection across the district. Additional potential risk factors, such as livestock trading and extensive breeding methods, are also discussed. Our results will contribute to the improvement of surveillance, control and eradication activities of TB in Sicily by the implementation of ad hoc TB control measures, especially in farms located along streams, sharing common pastures or with mixed animal species.

Whole Genome Sequencing Refines Knowledge on the Population Structure of Mycobacterium bovis from a Multi-Host Tuberculosis System

Classical molecular analyses of Mycobacterium bovis based on spoligotyping and Variable Number Tandem Repeat (MIRU-VNTR) brought the first insights into the epidemiology of animal tuberculosis (TB) in Portugal, showing high genotypic diversity of circulating strains that mostly cluster within the European 2 clonal complex. Previous surveillance provided valuable information on the prevalence and spatial occurrence of TB and highlighted prevalent genotypes in areas where livestock and wild ungulates are sympatric. However, links at the wildlife-livestock interfaces were established mainly via classical genotype associations. Here, we apply whole genome sequencing (WGS) to cattle, red deer and wild boar isolates to reconstruct the M. bovis population structure in a multi-host, multi-region disease system and to explore links at a fine genomic scale between M. bovis from wildlife hosts and cattle. Whole genome sequences of 44 representative M. bovis isolates, obtained between 2003 and 2015 from three TB hotspots, were compared through single nucleotide polymorphism (SNP) variant calling analyses. Consistent with previous results combining classical genotyping with Bayesian population admixture modelling, SNP-based phylogenies support the branching of this M. bovis population into five genetic clades, three with apparent geographic specificities, as well as the establishment of an SNP catalogue specific to each clade, which may be explored in the future as phylogenetic markers. The core genome alignment of SNPs was integrated within a spatiotemporal metadata framework to further structure this M. bovis population by host species and TB hotspots, providing a baseline for network analyses in different epidemiological and disease control contexts. WGS of M. bovis isolates from Portugal is reported for the first time in this pilot study, refining the spatiotemporal context of TB at the wildlife-livestock interface and providing further support to the key role of red deer and wild boar on disease maintenance. The SNP diversity observed within this dataset supports the natural circulation of M. bovis for a long time period, as well as multiple introduction events of the pathogen in this Iberian multi-host system.

Mycobacterium bovis Population Structure in Cattle and Local Badgers: Co-Localisation and Variation by Farm Type

Bovine tuberculosis surveillance in Northern Ireland includes Multiple-Locus Variable number tandem repeat Analysis (MLVA) to determine the Mycobacterium bovis genetic type present in both cattle and the predominant wildlife host, the European badger (Meles meles). These data are useful for investigating clusters of infection and understanding the scale at which interspecific transmission may occur. We utilised a comprehensive dataset of routinely sampled isolates from infected cattle and from badgers killed in road-traffic accidents to investigate the spatial co-location of MLVA types in, and between, the badger and cattle populations. Furthermore, we investigated the hypothesis that the type of farming enterprise might explain some variation in this relationship. MLVA types were spatially co-localised in cattle and road-traffic accident (RTA) badger hosts, indicative of a shared epidemic. Dairy herds were more likely to have at least one MLVA type in common with nearby RTA badgers, compared to non-dairy herd types. Marginally more MLVA spatial clustering was observed in non-dairy herds, which may be a consequence of relatively more between-herd movements. For the cattle population, local transmission mechanisms such as infection from contiguous herds, infectious wildlife and short-range between-herd cattle movements appear primarily to drive the epidemic: there appears to be a more limited role for long-range movements. Animal management practices are likely to be the driving force behind this observation, as beef rearing is associated with elevated numbers of animal movements compared to dairy herds.