Field-isolated genotypes of Mycobacterium bovis vary in virulence and influence case pathology but do not affect outbreak size

Serological test performance for bovine tuberculosis in cattle from herds with evidence of on-going infection in Northern Ireland

The ability to accurately identify infected hosts is the cornerstone of effective disease control and eradication programs. In the case of bovine tuberculosis, accurately identifying infected individual animals has been challenging as all available tests exhibit limited discriminatory ability. Here we assess the utility of two serological tests (IDEXX Mycobacterium bovis Ab test and Enfer multiplex antibody assay) and assess their performance relative to skin test (Single Intradermal Comparative Cervical Tuberculin; SICCT), gamma-interferon (IFNγ) and post-mortem results in a Northern Ireland setting. Furthermore, we describe a case-study where one test was used in conjunction with statutory testing. Serological tests using samples taken prior to SICCT disclosed low proportions of animals as test positive (mean 3% positive), despite the cohort having high proportions with positive SICCT test under standard interpretation (121/921; 13%) or IFNγ (365/922; 40%) results. Furthermore, for animals with a post-mortem record (n = 286), there was a high proportion with TB visible lesions (27%) or with laboratory confirmed infection (25%). As a result, apparent sensitivities within this cohort was very low (≤15%), however the tests succeeded in achieving very high specificities (96–100%). During the case-study, 7/670 (1.04%) samples from SICCT negative animals from a large chronically infected herd were serology positive, with a further 17 animals being borderline positive (17/670; 2.54%). Nine of the borderline animals were voluntarily removed, none of which were found to be infected post-mortem (no lesions/bacteriology negative). One serology test negative animal was subsequently found to have lesions at slaughter with M. bovis confirmed in the laboratory.

Global Distribution and Evolution of Mycobacterium bovis Lineages

Mycobacterium bovis is the main causative agent of zoonotic tuberculosis in humans and frequently devastates livestock and wildlife worldwide. Previous studies suggested the existence of genetic groups of M. bovis strains based on limited DNA markers (a.k.a. clonal complexes), and the evolution and ecology of this pathogen has been only marginally explored at the global level. We have screened over 2,600 publicly available M. bovis genomes and newly sequenced four wildlife M. bovis strains, gathering 1,969 genomes from 23 countries and at least 24 host species, including humans, to complete a phylogenomic analyses. We propose the existence of four distinct global lineages of M. bovis (Lb1, Lb2, Lb3, and Lb4) underlying the current disease distribution. These lineages are not fully represented by clonal complexes and are dispersed based on geographic location rather than host species. Our data divergence analysis agreed with previous studies reporting independent archeological data of ancient M. bovis (South Siberian infected skeletons at ∼2,000 years before present) and indicates that extant M. bovis originated between 715 and 3,556 years BP, with later emergence in the New World and Oceania, likely influenced by trades among countries.

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.

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

Investigating genetically-structured diversity in pathogen populations over time is important to better understand disease maintenance and spread. Herd-level surveillance of Mycobacterium bovis genotypes (multi-locus VNTR analysis types, MLVA types) from all culture-confirmed bovine tuberculosis (TB) herd cases was undertaken in Northern Ireland (NI), generating an unparalleled, longitudinal, population-level 14-year survey for this pathogen. Across this population, 295 genetically-distinct M. bovis MLVA types were identified in the 19,717 M. bovis isolates surveyed. Of these, the most frequent was MLVA type 002 (23.0%); 151 MLVA types were represented more than once, in groups ranging from 2 to 4438 isolates. Only 23 MLVA types were isolated in all 14 years. Investigating inter-annual frequency of M. bovis MLVA types, examples of statistically-significant expansions (MLVA types 002, 004, 006, 009 and 027), contractions (MLVA types 001, 007 and 011) and maintenance (MLVA types 003 and 005) were disclosed, during a period of fluctuating bovine TB herd-level incidence at the NI scale. The fixed period frequency distribution of MLVA types remained highly right-skewed. Novel VNTR copy number variant MLVA types (N = 242; an average of 17 per annum) were identified throughout the survey. The MLVA type distribution in the landscape was not random; MLVA types showed statistically-significant geographical localization and strong spatial associations with Divisional Veterinary Office (DVO) regions. There was also evidence of differential risk of particular MLVA types across breeds (Holstein/Friesian vs. other), age-class, and sex and some evidence of an association between the number of animals testing positive for bovine TB during the disclosing test and particular MLVA types, although there was substantial variation.

Genetically Related Mycobacterium bovis Strains Displayed Differential Intracellular Growth in Bovine Macrophages

Molecular typing of bacterial isolates provides a powerful approach for distinguishing Mycobacterium bovis (M. bovis) genotypes. It is known that M. bovis strain virulence plays a role in prevalence and spread of the disease, suggesting that strain virulence and prevailing genotypes are associated. However, it is not well understood whether strain virulence correlates with particular genotypes. In this study, we assessed the in vitro intracellular growth of 18 M. bovis isolates in bovine macrophages as an indicator of bacterial virulence and sought a relationship with the genotype identified by spoligotyping. We found 14 different spoligotypes—11 were already known and three spoligotypes had never been reported before. We identified 2 clusters that were phylogenetically related, containing 10 and 6 strains, respectively, and 2 orphan strains. Intracellular growth and phagocytic rates of 18 M. bovis strains were heterogeneous. Our results suggest that M. bovis intracellular growth and phagocytosis are independent of the bacterial lineage identified by spoligotyping.

Interleukin 10 knock-down in bovine monocyte-derived macrophages has distinct effects during infection with two divergent strains of Mycobacterium bovis

Mycobacterium bovis is the causative agent of bovine tuberculosis (TB), a cattle disease of global importance. M. bovis infects bovine macrophages (Mø) and subverts the host cell response to generate a suitable niche for survival and replication. We investigated the role of the anti-inflammatory cytokine interleukin (IL) 10 during in vitro infection of bovine monocyte-derived Mø (bMDM) with two divergent UK strains of M. bovis, which differentially modulate expression of IL10. The use of IL10-targeting siRNA revealed that IL10 inhibited the production of IL1B, IL6, tumour necrosis factor (TNF) and interferon gamma (IFNG) during infection of bMDM with the M. bovis strain G18. In contrast, IL10 only regulated a subset of these genes; TNF and IFNG, during infection with the M. bovis reference strain AF2122/97. Furthermore, nitric oxide (NO) production was modulated by IL10 during AF2122/97 infection, but not at the nitric oxide synthase 2 (NOS2) mRNA level, as observed during G18 infection. However, IL10 was found to promote survival of both M. bovis strains during early bMDM infection, but this effect disappeared after 24 h. The role of IL10-induced modulation of TNF, IFNG and NO production in M. bovis survival was investigated using siRNA targeting TNF, IFNG receptor 1 (IFNGR1) and NOS2. Knock-down of these genes individually did not promote survival of either M. bovis strain and therefore modulation of these genes does not account for the effect of IL10 on M. bovis survival. However, TNF knock-down was found to be detrimental to the survival of the M. bovis strain G18 during early infection. The results provide further evidence for the importance of IL10 during M. bovis infection of Mø. Furthermore, they highlight M. bovis strain specific differences in the interaction with the infected bMDM, which may influence the course of infection and progression of bovine TB.

Variation in Mycobacterium bovis genetic richness suggests that inwards cattle movements are a more important source of infection in beef herds than in dairy herds

Background

We used genetic Multi-Locus VNTR Analysis (MLVA) data gathered from surveillance efforts to better understand the ongoing bovine tuberculosis (bTB) epidemic in Northern Irish cattle herds. We modelled the factors associated with Mycobacterium bovis MLVA genotype richness at three analytical scales; breakdown level, herd level, and patch level, and compared the results between dairy and non-dairy production types.

Results

In 83% of breakdowns and in 63% of herds, a single MLVA genotype was isolated. Five or more MLVA genotypes were found in less than 3 % of herds. Herd size and the total number of reactors were important explanatory variables, suggesting that increasing MLVA genotype richness was positively related to increases in the number of host animals. Despite their smaller relative size, however, the highest MLVA genotype richness values were observed in non-dairy herds. Increasing inwards cattle movements were important positive predictors of MLVA genotype richness, but mainly in non-dairy settings.

Conclusions

The principal finding is that low MLVA genotype richness indicates that small-scale epidemics, e.g. wildlife, contiguous farms, and within-herd recrudescence, are important routes of M. bovis infection in cattle herds. We hypothesise that these mechanisms will maintain, but may not explicitly increase, MLVA genotype richness. The presence of elevated MLVA richness is relatively rare and likely indicates beef fattening enterprises, which purchase cattle from over long distances. Cattle movements were furthermore an important predictor of MLVA genotype richness in non-dairy herds, but not in dairy herds; this may represent reduced cattle purchasing levels in dairy enterprises, compared to beef. These observations allude to the relative contribution of different routes of bTB infection between production types; we posit that infection associated with local factors may be more evident in dairy herds than beef herds, however in beef herds, inwards movements offer additional opportunities for introducing M. bovis into the herd.

Electronic supplementary material

The online version of this article (10.1186/s12866-019-1530-7) contains supplementary material, which is available to authorized users.

Spatiotemporal analysis of prolonged and recurrent bovine tuberculosis breakdowns in Northern Irish cattle herds reveals a new infection hotspot

Despite a state-led eradication programme, bovine tuberculosis (bTB) remains endemic in Northern Ireland (NI). Of particular concern are "chronic" prolonged and recurrent bTB breakdowns, which represent significant financial and administrative burdens. However, little is known regarding the spatiotemporal distribution of chronic breakdowns in NI. We therefore analysed both the spatial and spatiotemporal distributions of chronic bTB breakdowns between 2004 and 2014. Significantly positive values for Moran's Index of spatial autocorrelation were found, and Local Moran's I clustering was employed to assess for spatial associations in the number and prevalence of chronic bTB breakdowns across NI. Additional spatio-temporal analysis using SaTScan showed that the burden of chronic bTB infection tends to be found where bTB levels are already high. However, a novel hotspot was revealed wherein the prevalence of chronic breakdowns was higher than expected; this should be the subject of follow-up surveillance.

Bovine Tuberculosis in Britain and Ireland - A Perfect Storm? the Confluence of Potential Ecological and Epidemiological Impediments to Controlling a Chronic Infectious Disease

Successful eradication schemes for bovine tuberculosis (bTB) have been implemented in a number of European and other countries over the last 50 years. However, the islands of Britain and Ireland remain a significant aberration to this trend, with the recent exception of Scotland. Why have eradication schemes failed within these countries, while apparently similar programs have been successful elsewhere? While significant socio-economic and political factors have been discussed elsewhere as key determinants of disease eradication, here we review some of the potential ecological and epidemiological constraints that are present in these islands relative to other parts of Europe. We argue that the convergence of these potential factors may interact additively to diminish the potential of the present control programs to achieve eradication. Issues identified include heterogeneity of diagnostic testing approaches, the presence of an abundant wildlife reservoir of infection and the challenge of sustainably managing this risk effectively; the nature, size, density and network structure of cattle farming; potential effects of Mycobacterium bovis strain heterogeneity on disease transmission dynamics; possible impacts of concurrent endemic infections on the disclosure of truly infected animals; climatological differences and change coupled with environmental contamination. We further argue that control and eradication of this complex disease may benefit from an ecosystem level approach to management. We hope that this perspective can stimulate a new conversation about the many factors potentially impacting bTB eradication schemes in Britain and Ireland and possibly stimulate new research in the areas identified.

Variation in the Early Host-Pathogen Interaction of Bovine Macrophages with Divergent Mycobacterium bovis Strains in the United Kingdom

Bovine tuberculosis has been an escalating animal health issue in the United Kingdom since the 1980s, even though control policies have been in place for over 60 years. The importance of the genetics of the etiological agent, Mycobacterium bovis, in the reemergence of the disease has been largely overlooked. We compared the interaction between bovine monocyte-derived macrophages (bMDM) and two M. bovis strains, AF2122/97 and G18, representing distinct genotypes currently circulating in the United Kingdom. These M. bovis strains exhibited differences in survival and growth in bMDM. Although uptake was similar, the number of viable intracellular AF2122/97 organisms increased rapidly, while G18 growth was constrained for the first 24 h. AF2122/97 infection induced a greater transcriptional response by bMDM than G18 infection with respect to the number of differentially expressed genes and the fold changes measured. AF2122/97 infection induced more bMDM cell death, with characteristics of necrosis and apoptosis, more inflammasome activation, and a greater type I interferon response than G18. In conclusion, the two investigated M. bovis strains interact in significantly different ways with the host macrophage. In contrast to the relatively silent infection by G18, AF2122/97 induces greater signaling to attract other immune cells and induces host cell death, which may promote secondary infections of naive macrophages. These differences may affect early events in the host-pathogen interaction, including granuloma development, which could in turn alter the progression of the disease. Therefore, the potential involvement of M. bovis genotypes in the reemergence of bovine tuberculosis in the United Kingdom warrants further investigation.

Modelling the variation in skin-test tuberculin reactions, post-mortem lesion counts and case pathology in tuberculosis-exposed cattle: Effects of animal characteristics, histories and co-infection

Correctly identifying bovine tuberculosis (bTB) in cattle remains a significant problem in endemic countries. We hypothesized that animal characteristics (sex, age, breed), histories (herd effects, testing, movement) and potential exposure to other pathogens (co-infection; BVDV, liver fluke and Mycobacterium avium reactors) could significantly impact the immune responsiveness detected at skin testing and the variation in post-mortem pathology (confirmation) in bTB-exposed cattle. Three model suites were developed using a retrospective observational data set of 5,698 cattle culled during herd breakdowns in Northern Ireland. A linear regression model suggested that antemortem tuberculin reaction size (difference in purified protein derivative avium [PPDa] and bovine [PPDb] reactions) was significantly positively associated with post-mortem maximum lesion size and the number of lesions found. This indicated that reaction size could be considered a predictor of both the extent (number of lesions/tissues) and the pathological progression of infection (maximum lesion size). Tuberculin reaction size was related to age class, and younger animals (<2.85 years) displayed larger reaction sizes than older animals. Tuberculin reaction size was also associated with breed and animal movement and increased with the time between the penultimate and disclosing tests. A negative binomial random-effects model indicated a significant increase in lesion counts for animals with M. avium reactions (PPDb-PPDa < 0) relative to non-reactors (PPDb-PPDa = 0). Lesion counts were significantly increased in animals with previous positive severe interpretation skin-test results. Animals with increased movement histories, young animals and non-dairy breed animals also had significantly increased lesion counts. Animals from herds that had BVDV-positive cattle had significantly lower lesion counts than animals from herds without evidence of BVDV infection. Restricting the data set to only animals with a bTB visible lesion at slaughter (n = 2471), an ordinal regression model indicated that liver fluke-infected animals disclosed smaller lesions, relative to liver fluke-negative animals, and larger lesions were disclosed in animals with increased movement histories.

Bovine tuberculosis visible lesions in cattle culled during herd breakdowns: the effects of individual characteristics, trade movement and co-infection

BACKGROUND

Bovine tuberculosis (bTB), caused by Mycobacterium bovis, remains a significant problem for livestock industries in many countries worldwide including Northern Ireland, where a test and slaughter regime has utilised the Single Intradermal Comparative Cervical Tuberculin (SICCT) test since 1959. We investigated the variation in post-mortem confirmation based on bTB visible lesion (VL) presence during herd breakdowns using two model suites. We investigated animal-level characteristics, while controlling for herd-level factors and clustering. We were interested in potential impacts of concurrent infection, and therefore we assessed whether animals with evidence of liver fluke infection (Fasciola hepatica; post-mortem inspection), M. avium reactors (animals with negative M. bovis-avium (b-a) tuberculin reactions) or Bovine Viral Diarrhoea Virus (BVDV; RT-PCR tested) were associated with bTB confirmation.

RESULTS

The dataset included 6242 animals removed during the 14 month study period (2013-2015). bTB-VL presence was significantly increased in animals with greater b-a reaction size at the disclosing SICCT test (e.g. b-a = 5-9 mm vs. b-a = 0 mm, adjusted Odds ratio (aOR): 14.57; p < 0.001). M. avium reactor animals (b-a < 0) were also significantly more likely to disclose VL than non-reactor animals (b-a = 0; aOR: 2.29; p = 0.023). Animals had a greater probability of exhibiting lesions with the increasing number of herds it had resided within (movement; log-herds: aOR: 2.27-2.42; p < 0.001), if it had an inconclusive penultimate test result (aOR: 2.84-3.89; p < 0.001), and with increasing time between tests (log-time; aOR: 1.23; p = 0.003). Animals were less likely to have VL if they were a dairy breed (aOR: 0.79; p = 0.015) or in an older age-class (e.g. age-quartile 2 vs. 4; aOR: 0.65; p < 0.001). Liver fluke or BVDV variables were not retained in either multivariable model as they were non-significantly associated with bTB-VL status (p > 0.1).

CONCLUSIONS

Our results suggest that neither co-infection of liver fluke nor BVDV had a significant effect on the presence of VLs in this high-risk cohort. M. avium tuberculin reactors had a significantly increased risk of disclosing with a bTB lesion, which could be related to the impact of co-infection with M. avium subsp. paratuberculosis (MAP) affecting the performance of the SICCT however further research in this area is required. Movements, test history, breed and age were important factors influencing confirmation in high-risk animals.

Complete Genome Sequencing of Mycobacterium bovis SP38 and Comparative Genomics of Mycobacterium bovis and M. tuberculosis Strains

Mycobacterium bovis causes bovine tuberculosis and is the main organism responsible for zoonotic tuberculosis in humans. We performed the sequencing, assembly and annotation of a Brazilian strain of M. bovis named SP38, and performed comparative genomics of M. bovis genomes deposited in GenBank. M. bovis SP38 has a traditional tuberculous mycobacterium genome of 4,347,648 bp, with 65.5% GC, and 4,216 genes. The majority of CDSs (2,805, 69.3%) have predictive function, while 1,206 (30.07%) are hypothetical. For comparative analysis, 31 M. bovis, 32 M. bovis BCG, and 23 Mycobacterium tuberculosis genomes available in GenBank were selected. M. bovis RDs (regions of difference) and Clonal Complexes (CC) were identified in silico. Genome dynamics of bacterial groups were analyzed by gene orthology and polymorphic sites identification. M. bovis polymorphic sites were used to construct a phylogenetic tree. Our RD analyses resulted in the exclusion of three genomes, mistakenly annotated as virulent M. bovis. M. bovis SP38 along with strain 35 represent the first report of CC European 2 in Brazil, whereas two other M. bovis strains failed to be classified within current CC. Results of M. bovis orthologous genes analysis suggest a process of genome remodeling through genomic decay and gene duplication. Quantification, pairwise comparisons and distribution analyses of polymorphic sites demonstrate greater genetic variability of M. tuberculosis when compared to M. bovis and M. bovis BCG (p ≤ 0.05), indicating that currently defined M. tuberculosis lineages are more genetically diverse than M. bovis CC and animal-adapted MTC (M. tuberculosis Complex) species. As expected, polymorphic sites annotation shows that M. bovis BCG are subjected to different evolutionary pressures when compared to virulent mycobacteria. Lastly, M. bovis phylogeny indicates that polymorphic sites may be used as markers of M. bovis lineages in association with CC. Our findings highlight the need to better understand host-pathogen co-evolution in genetically homogeneous and/or diverse host populations, considering the fact that M. bovis has a broader host range when compared to M. tuberculosis. Also, the identification of M. bovis genomes not classified within CC indicates that the diversity of M. bovis lineages may be larger than previously thought or that current classification should be reviewed.

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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12917-017-1190-5) contains supplementary material, which is available to authorized users.

Assessment of concurrent infection with bovine viral diarrhoea virus (BVDV) and Mycobacterium bovis: A herd-level risk factor analysis from Northern Ireland

Bovine viral diarrhoea virus (BVDV) is a significant pathogen of cattle, leading to severe economic and animal-welfare impacts. Furthermore, the pathogen has been associated with impacting the progression or spread of other pathogens (e.g. Mycobacterium bovis, the causative agent of bovine tuberculosis (bTB)). During this study we investigated (i) risk factors for BVDV at a herd-level and (ii) whether there was any association between BVDV and herd-level bTB risk. The data for this study were gathered from a voluntary BVDV control programme in Northern Ireland (2013-2015) based on the identification of virus positive animals through tissue tag testing of calves. We assigned a herd-level BVDV status to 2827 participating herds, where a herd was assumed "infected" if one or more animals tested positive for BVDV. Two model suites were developed. Firstly, we assessed risk factors for BVDV herd status using multivariable logit random-effects modelling, aggregating to the calendar year level (2013-2015; n=4828; model 1). Secondly, we aggregated data across the three years of the study to give an overall status for the whole study period (n=2827; logistic model 2). Risk factors included year, herd-type, herd size, number of births, inward trade moves, calf mortality, and region. Furthermore, the herd-level bovine tuberculosis status (based on the single intradermal comparative cervical tuberculin (SICCT) test outcomes, or confirmation at post-mortem), or the size of bTB breakdowns (number of SICCT test positive animals), of herds was also investigated to assess whether there was an association (co-infection) with herd BVDV status. The final models suggested that BVDV herd status was positively associated with increased levels of calf mortality, herd size, number of births, the number of BVDV tests undertaken and the number of animals introduced to the herd. There was a significant univariable positive association between BVDV status, and SICCT breakdown risk, breakdown size and confirmed bTB status in model 2. However, there was no evidence of significant associations between bTB status (using SICTT status, confirmed status or herd breakdown size) and BVDV status in final multivariable models when controlling for other significant confounders. These results provide information for action for the future control and eradication of BVDV in Northern Ireland, though these data provide little support for the hypothesised association between BVDV and bTB status at herd-level. Further animal-level analyses are necessary to investigate whether there is support for a BVD-bTB co-infection association, including the impact of co-infection on the severity of infection.

Use of bacterial whole-genome sequencing to investigate local persistence and spread in bovine tuberculosis

•

We performed whole genome sequencing (WGS) of Mycobacterium bovis for a single molecular (VNTR) type.

•

Under-sampling of one lineage was caused by switching between VNTR-types.

•

Pairwise SNP distances showed a weak genetic isolation by distance pattern.

•

Bayesian phylogeographic inference was feasible despite a low substitution rate.

•

WGS studies of M. bovis need to account for slow evolution and molecular type switching.

Mycobacterium bovis is the causal agent of bovine tuberculosis, one of the most important diseases currently facing the UK cattle industry. Here, we use high-density whole genome sequencing (WGS) in a defined sub-population of M. bovis in 145 cattle across 66 herd breakdowns to gain insights into local spread and persistence. We show that despite low divergence among isolates, WGS can in principle expose contributions of under-sampled host populations to M. bovis transmission. However, we demonstrate that in our data such a signal is due to molecular type switching, which had been previously undocumented for M. bovis. Isolates from farms with a known history of direct cattle movement between them did not show a statistical signal of higher genetic similarity. Despite an overall signal of genetic isolation by distance, genetic distances also showed no apparent relationship with spatial distance among affected farms over distances <5 km. Using simulations, we find that even over the brief evolutionary timescale covered by our data, Bayesian phylogeographic approaches are feasible. Applying such approaches showed that M. bovis dispersal in this system is heterogeneous but slow overall, averaging 2 km/year. These results confirm that widespread application of WGS to M. bovis will bring novel and important insights into the dynamics of M. bovis spread and persistence, but that the current questions most pertinent to control will be best addressed using approaches that more directly integrate WGS with additional epidemiological data.

Genetic evolution of Mycobacterium bovis causing tuberculosis in livestock and wildlife in France since 1978

To study the dynamics of bovine tuberculosis (bTB) in France, 4,654 M. bovis strains isolated mainly from livestock and wildlife since 1978 were characterized by spoligotyping and MLVA based on MIRU-VNTR. In our study spoligotyping allowed the discrimination of 176 types although 3 spoligotypes are predominant and account for more than half of the total strain population: SB0120 (26%), SB0134 (11%) and SB0121 (6%). In addition, 11% of the isolates, principally from Southern France, showing close spoligotypes and MIRU-VNTR types have been gathered in a family designated as the "F4-family". MLVA typing allowed extensive discrimination, particularly for strains with predominant spoligotypes, with a total of 498 genotypes, several of which were highly regionalized. The similarity of the strains' genetic relationships based on spoligotyping and MIRU-VNTR markers supports the co-existence of different clonal populations within the French M. bovis population. A genetic evolution of the strains was observed both geographically and in time. Indeed, as a result of the reduction of bTB due to the national control campaigns, a large reduction of the strains' genetic variability took place in the last ten years. However, in the regions were bTB is highly prevalent at present, cases in both livestock and in wildlife are due to the spread of unique local genotype profiles. Our results show that the highly discriminating genotyping tools used in this study for molecular studies of bTB are useful for addressing pending questions, which would lead to a better insight into the epidemiology of the disease, and for finding proper solutions for its sustainable control in France.

Spatial dynamics of bovine tuberculosis in the Autonomous Community of Madrid, Spain (2010-2012)

Progress in control of bovine tuberculosis (bTB) is often not uniform, usually due to the effect of one or more sometimes unknown epidemiological factors impairing the success of eradication programs. Use of spatial analysis can help to identify clusters of persistence of disease, leading to the identification of these factors thus allowing the implementation of targeted control measures, and may provide some insights of disease transmission, particularly when combined with molecular typing techniques. Here, the spatial dynamics of bTB in a high prevalence region of Spain were assessed during a three year period (2010-2012) using data from the eradication campaigns to detect clusters of positive bTB herds and of those infected with certain Mycobacterium bovis strains (characterized using spoligotyping and VNTR typing). In addition, the within-herd transmission coefficient (β) was estimated in infected herds and its spatial distribution and association with other potential outbreak and herd variables was evaluated. Significant clustering of positive herds was identified in the three years of the study in the same location ("high risk area"). Three spoligotypes (SB0339, SB0121 and SB1142) accounted for >70% of the outbreaks detected in the three years. VNTR subtyping revealed the presence of few but highly prevalent strains within the high risk area, suggesting maintained transmission in the area. The spatial autocorrelation found in the distribution of the estimated within-herd transmission coefficients in herds located within distances <14 km and the results of the spatial regression analysis, support the hypothesis of shared local factors affecting disease transmission in farms located at a close proximity.