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

Background Bovine tuberculosis (bTB) is a costly, epidemiologically complex, multi-host, endemic disease. Pathogen whole genome sequencing can improve the resolution of epidemiological tracing. We genome sequenced an exceptional data set of 619 Mycobacterium bovis isolates from badgers and cattle in a 100km2 bTB ‘hotspot’. Historical molecular subtyping data permitted the targeting of an endemic pathogen lineage, whose long-term persistence provided an opportunity to study genome epidemiology in detail. 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. Results Eight lineages of M. bovis were circulating in the study area, seven of which were likely non-endemic, and imported by animal movement. The endemic lineage exhibited low genetic diversity with an average inter-isolate genetic distance of 7.6 SNPs (s.d. ± 4.0), consistent with contemporary transmission. Bayesian phylogenetic methods determined an evolutionary rate of 0.30 substitutions per genome per year for this lineage, estimating its emergence 40-50 years before present, while Bayesian Skyline analysis identified significant population expansion of the endemic lineage in the 1990s and again in 2011-2012. The phylogeny revealed distinct sub-lineages, all of which contained isolates from both cattle and badger hosts, indicative of the sharing of closely related strains and inter-species transmission. However, the presence of significant badger population genetic structure was not associated with the spatial distribution of M. bovis genetic diversity. Conclusions Our data provided unparalleled detail on the evolutionary history of an endemic M. bovis lineage. Findings are consistent with ongoing interspecies transmission in the study area but suggest that badger intra-species transmission may not be a major driver of persistence in this area. In addition, the data collected permitted the tracking of incursions of novel pathogen lineages into the study area and means to determine if they were involved in disease transmission. ### Competing Interest Statement The authors have declared no competing interest.