1887

Abstract

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 (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 . 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.

Funding
This study was supported by the:
  • Department of Agriculture, Environment and Rural Affairs, UK Government (Award 21/3/05)
    • Principle Award Recipient: AdrianAllen
  • This is an open-access article distributed under the terms of the Creative Commons Attribution License.
Loading

Article metrics loading...

/content/journal/mgen/10.1099/mgen.0.001185
2024-02-14
2024-05-04
Loading full text...

Full text loading...

/deliver/fulltext/mgen/10/2/mgen001185.html?itemId=/content/journal/mgen/10.1099/mgen.0.001185&mimeType=html&fmt=ahah

References

  1. Skuce RA, Byrne AW, Lahuerta A, Allen AR. Mycobacterium bovis molecular typing and surveillance. In Chambers M, Gordon SV, Olea-Popelka F, Barrow P. eds Bovine Tuberculosis CABI; 2018 [View Article]
    [Google Scholar]
  2. Kamerbeek J, Schouls L, Kolk A, van Agterveld M, van Soolingen D et al. Simultaneous detection and strain differentiation of Mycobacterium tuberculosis for diagnosis and epidemiology. J Clin Microbiol 1997; 35:907–914 [View Article] [PubMed]
    [Google Scholar]
  3. Smith NH, Upton P. Naming spoligotype patterns for the RD9-deleted lineage of the Mycobacterium tuberculosis complex; www.Mbovis.org. Infect Genet Evol 2012; 12:873–876 [View Article] [PubMed]
    [Google Scholar]
  4. Skuce RA, McDowell SW, Mallon TR, Luke B, Breadon EL et al. Discrimination of isolates of Mycobacterium bovis in Northern Ireland on the basis of variable numbers of tandem repeats (VNTRs). Vet Rec 2005; 157:501–504 [View Article] [PubMed]
    [Google Scholar]
  5. Olea-Popelka FJ, Flynn O, Costello E, McGrath G, Collins JD et al. Spatial relationship between Mycobacterium bovis strains in cattle and badgers in four areas in Ireland. Prev Vet Med 2005; 71:57–70 [View Article] [PubMed]
    [Google Scholar]
  6. Skuce RA, Mallon TR, McCormick CM, McBride SH, Clarke G et al. Mycobacterium bovis genotypes in Northern Ireland: herd-level surveillance (2003 to 2008). Vet Rec 2010; 167:684–689 [View Article] [PubMed]
    [Google Scholar]
  7. Skuce R, Breadon E, Allen A, Milne G, McCormick C et al. Longitudinal dynamics of herd-level Mycobacterium bovis MLVA type surveillance in cattle in Northern Ireland 2003-2016. Infect Genet Evol 2020; 79:104131 [View Article] [PubMed]
    [Google Scholar]
  8. Meehan CJ, Moris P, Kohl TA, Pečerska J, Akter S et al. The relationship between transmission time and clustering methods in Mycobacterium tuberculosis epidemiology. EBioMedicine 2018; 37:410–416 [View Article] [PubMed]
    [Google Scholar]
  9. Pearson T, Okinaka RT, Foster JT, Keim P. Phylogenetic understanding of clonal populations in an era of whole genome sequencing. Infect Genet Evol 2009; 9:1010–1019 [View Article] [PubMed]
    [Google Scholar]
  10. Gardy JL, Loman NJ. Towards a genomics-informed, real-time, global pathogen surveillance system. Nat Rev Genet 2018; 19:9–20 [View Article] [PubMed]
    [Google Scholar]
  11. Kao RR, Price-Carter M, Robbe-Austerman S. Use of genomics to track bovine tuberculosis transmission. Rev Sci Tech 2016; 35:241–258 [View Article] [PubMed]
    [Google Scholar]
  12. Biek R, O’Hare A, Wright D, Mallon T, McCormick C et al. Whole genome sequencing reveals local transmission patterns of Mycobacterium bovis in sympatric cattle and badger populations. PLoS Pathog 2012; 8:e1003008 [View Article] [PubMed]
    [Google Scholar]
  13. Crispell J, Zadoks RN, Harris SR, Paterson B, Collins DM et al. Using whole genome sequencing to investigate transmission in a multi-host system: bovine tuberculosis in New Zealand. BMC Genomics 2017; 18:180 [View Article] [PubMed]
    [Google Scholar]
  14. Salvador LCM, O’Brien DJ, Cosgrove MK, Stuber TP, Schooley AM et al. Disease management at the wildlife-livestock interface: using whole-genome sequencing to study the role of elk in Mycobacterium bovis transmission in Michigan, USA. Mol Ecol 2019; 28:2192–2205 [View Article] [PubMed]
    [Google Scholar]
  15. Crispell J, Benton CH, Balaz D, De Maio N, Ahkmetova A et al. Combining genomics and epidemiology to analyse bi-directional transmission of Mycobacterium bovis in a multi-host system. Elife 2019; 8:e45833 [View Article] [PubMed]
    [Google Scholar]
  16. Crispell J, Cassidy S, Kenny K, McGrath G, Warde S et al. Mycobacterium bovis genomics reveals transmission of infection between cattle and deer in Ireland. Microb Genom 2020; 6:mgen000388 [View Article] [PubMed]
    [Google Scholar]
  17. van Tonder AJ, Thornton MJ, Conlan AJK, Jolley KA, Goolding L et al. Inferring Mycobacterium bovis transmission between cattle and badgers using isolates from the Randomised Badger Culling Trial. PLoS Pathog 2021; 17:e1010075 [View Article] [PubMed]
    [Google Scholar]
  18. Rossi G, Crispell J, Brough T, Lycett SJ, White PCL et al. Phylodynamic analysis of an emergent Mycobacterium bovis outbreak in an area with no previously known wildlife infections. J Appl Ecol 2022; 59:210–222 [View Article]
    [Google Scholar]
  19. Canini L, Modenesi G, Courcoul A, Boschiroli M, Durand B et al. Deciphering the role of host species for two Mycobacterium bovis genotypes from the European 3 clonal complex circulation within a cattle‐badger‐wild boar multihost system. MicrobiologyOpen 2023; 12: [View Article]
    [Google Scholar]
  20. Akhmetova A, Guerrero J, McAdam P, Salvador LCM, Crispell J et al. Genomic epidemiology of Mycobacterium bovis infection in sympatric badger and cattle populations in Northern Ireland. Microb Genom 2023; 9:mgen001023 [View Article] [PubMed]
    [Google Scholar]
  21. Biek R, Pybus OG, Lloyd-Smith JO, Didelot X. Measurably evolving pathogens in the genomic era. Trends Ecol Evol 2015; 30:306–313 [View Article] [PubMed]
    [Google Scholar]
  22. Lewis PO, Holder MT, Holsinger KE. Polytomies and Bayesian phylogenetic inference. Syst Biol 2005; 54:241–253 [View Article]
    [Google Scholar]
  23. Didelot X, Fraser C, Gardy J, Colijn C. Genomic infectious disease epidemiology in partially sampled and ongoing outbreaks. Mol Biol Evol 2017; 34:msw075 [View Article]
    [Google Scholar]
  24. Nelson KN, Talarico S, Poonja S, McDaniel CJ, Cilnis M et al. Mutation of Mycobacterium tuberculosis and implications for using whole-genome sequencing for investigating recent tuberculosis transmission. Front Public Health 2021; 9:790544 [View Article] [PubMed]
    [Google Scholar]
  25. Firth C, Kitchen A, Shapiro B, Suchard MA, Holmes EC et al. Using time-structured data to estimate evolutionary rates of double-stranded DNA viruses. Mol Biol Evol 2010; 27:2038–2051 [View Article] [PubMed]
    [Google Scholar]
  26. Perets V, Allen A, Crispell J, Cassidy S, O’Connor A et al. Evidence for local and international spread of Mycobacterium avium subspecies paratuberculosis through whole genome sequencing of isolates from the island of Ireland. Vet Microbiol 2022; 268:109416 [View Article] [PubMed]
    [Google Scholar]
  27. Allen AR, Dale J, McCormick C, Mallon TR, Costello E et al. The phylogeny and population structure of Mycobacterium bovis in the British Isles. Infect Genet Evol 2013; 20:8–15 [View Article] [PubMed]
    [Google Scholar]
  28. Courcier EA, Menzies FD, Strain SAJ, Skuce RA, Robinson PA et al. Mycobacterium bovis surveillance in Eurasian badgers (Meles Meles) killed by vehicles in Northern Ireland between 1998 and 2011. Vet Rec 2018259–265 [View Article] [PubMed]
    [Google Scholar]
  29. Parish T, Stoker NG. Mycobacterium Tuberculosis Protocols New Jersey: Humana Press; 2001 [View Article]
    [Google Scholar]
  30. van Soolingen D, de Haas PE, Kremer K. Restriction fragment length polymorphism typing of mycobacteria. Methods Mol Med 2001; 54:165–203 [View Article] [PubMed]
    [Google Scholar]
  31. Malone KM, Farrell D, Stuber TP, Schubert OT, Aebersold R et al. Updated reference genome sequence and annotation of Mycobacterium bovis AF2122/97. Genome Announc 2017; 5:14 [View Article] [PubMed]
    [Google Scholar]
  32. Edwards DB, Pope B, Holt K. RedDog mapping based phylogenomics pipeline; 2016 https://github.com/katholt/RedDog accessed 1 September 2023
  33. Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods 2012; 9:357–359 [View Article] [PubMed]
    [Google Scholar]
  34. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J et al. The sequence alignment/Map format and SAMtools. Bioinformatics 2009; 25:2078–2079 [View Article] [PubMed]
    [Google Scholar]
  35. Schliep KP. phangorn: phylogenetic analysis in R. Bioinformatics 2011; 27:592–593 [View Article] [PubMed]
    [Google Scholar]
  36. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria; 2008 http://www.R-project.org
  37. Kozlov AM, Darriba D, Flouri T, Morel B, Stamatakis A. RAxML-NG: a fast, scalable and user-friendly tool for maximum likelihood phylogenetic inference. Bioinformatics 2019; 35:4453–4455 [View Article] [PubMed]
    [Google Scholar]
  38. Yu G. ggtree: Elegant Graphics for Phylogenetic Tree Visualization and Annotation 2020. https://guangchuangyu.github.io/ggtree-book/short-introduction-to-r.html
  39. Walker TM, Ip CLC, Harrell RH, Evans JT, Kapatai G et al. Whole-genome sequencing to delineate Mycobacterium tuberculosis outbreaks: a retrospective observational study. Lancet Infect Dis 2013; 13:137–146 [View Article] [PubMed]
    [Google Scholar]
  40. Bryant JM, Schürch AC, van Deutekom H, Harris SR, de Beer JL et al. Inferring patient to patient transmission of Mycobacterium tuberculosis from whole genome sequencing data. BMC Infect Dis 2013; 13:110 [View Article] [PubMed]
    [Google Scholar]
  41. Roetzer A, Diel R, Kohl TA, Rückert C, Nübel U et al. Whole genome sequencing versus traditional genotyping for investigation of a Mycobacterium tuberculosis outbreak: a longitudinal molecular epidemiological study. PLoS Med 2013; 10:e1001387 [View Article] [PubMed]
    [Google Scholar]
  42. Yang C, Luo T, Shen X, Wu J, Gan M et al. Transmission of multidrug-resistant Mycobacterium tuberculosis in Shanghai, China: a retrospective observational study using whole-genome sequencing and epidemiological investigation. Lancet Infect Dis 2017; 17:275–284 [View Article] [PubMed]
    [Google Scholar]
  43. Jajou R, de Neeling A, van Hunen R, de Vries G, Schimmel H et al. Epidemiological links between tuberculosis cases identified twice as efficiently by whole genome sequencing than conventional molecular typing: a population-based study. PLoS One 2018; 13:e0195413 [View Article] [PubMed]
    [Google Scholar]
  44. Balaban M, Moshiri N, Mai U, Jia X, Mirarab S. TreeCluster: clustering biological sequences using phylogenetic trees. PLoS One 2019; 14:e0221068 [View Article] [PubMed]
    [Google Scholar]
  45. Péron G. Modified home range kernel density estimators that take environmental interactions into account. Mov Ecol 2019; 7:16 [View Article]
    [Google Scholar]
  46. O’Hagan MJH, Gordon AW, McCormick CM, Collins SF, Trimble NA et al. Effect of selective removal of badgers (Meles meles) on ranging behaviour during a “Test and Vaccinate or Remove” intervention in Northern Ireland. Epidemiol Infect 2021; 149:e125 [View Article] [PubMed]
    [Google Scholar]
  47. Wickham H. ggplot2. In Ggplot2: Elegant Graphics for Data Analysis Cham: Springer-Verlag; 2016 [View Article]
    [Google Scholar]
  48. Otto J, Kahle D. ggdensity: improved bivariate density visualization in R. The R Journal 2023; 15:220–236 [View Article]
    [Google Scholar]
  49. Sheather SJ, Jones MC. A reliable data‐based bandwidth selection method for kernel density estimation. J Roy Stat Soc B 1991; 53:683–690 [View Article]
    [Google Scholar]
  50. Harpole JK, Woods CM, Rodebaugh TL, Levinson CA, Lenze EJ. How bandwidth selection algorithms impact exploratory data analysis using kernel density estimation. Psychol Methods 2014; 19:428–443 [View Article] [PubMed]
    [Google Scholar]
  51. Pebesma E. Simple features for R: standardized support for Spatial vector data. The R Journal 2018; 10:439 [View Article]
    [Google Scholar]
  52. Bivand RS, Pebesma E, Gómez-Rubio V. Applied Spatial Data Analysis with R. In Applied Spatial Data Analysis with R, Second edn New York, NY: Springer; 2013 [View Article]
    [Google Scholar]
  53. Houston R. A computerised database system for bovine traceability. Rev Sci Tech 2001; 20:652–661 [View Article] [PubMed]
    [Google Scholar]
  54. Rambaut A, Lam TT, Max Carvalho L, Pybus OG. Exploring the temporal structure of heterochronous sequences using TempEst (formerly Path-O-Gen). Virus Evol 2016; 2:vew007 [View Article] [PubMed]
    [Google Scholar]
  55. Rieux A, Khatchikian CE. tipdatingbeast: an r package to assist the implementation of phylogenetic tip-dating tests using beast. Mol Ecol Resour 2017; 17:608–613 [View Article] [PubMed]
    [Google Scholar]
  56. Bouckaert R, Heled J, Kühnert D, Vaughan T, Wu C-H et al. BEAST 2: a Software platform for Bayesian evolutionary analysis. PLoS Comput Biol 2014; 10:e1003537 [View Article]
    [Google Scholar]
  57. Rambaut A, Drummond AJ, Xie D, Baele G, Suchard MA. Posterior summarization in Bayesian phylogenetics using tracer 1.7. Syst Biol 2018; 67:901–904 [View Article] [PubMed]
    [Google Scholar]
  58. Drummond AJ, Rambaut A, Shapiro B, Pybus OG. Bayesian coalescent inference of past population dynamics from molecular sequences. Mol Biol Evol 2005; 22:1185–1192 [View Article] [PubMed]
    [Google Scholar]
  59. Boskova V, Bonhoeffer S, Stadler T. Inference of epidemiological dynamics based on simulated phylogenies using birth-death and coalescent models. PLoS Comput Biol 2014; 10:e1003913 [View Article]
    [Google Scholar]
  60. Smith NH, Dale J, Inwald J, Palmer S, Gordon SV et al. The population structure of Mycobacterium bovis in Great Britain: clonal expansion. Proc Natl Acad Sci U S A 2003; 100:15271–15275 [View Article] [PubMed]
    [Google Scholar]
  61. Roper T. Badger London: Harper Collins; 2010
    [Google Scholar]
  62. Allen A, Milne G, McCormick C, Collins S, O’Hagan M et al. European badger (Meles meles) responses to low‐intensity, selective culling: using mark–recapture and relatedness data to assess social perturbation. Ecol Sol and Evidence 2022; 3: [View Article]
    [Google Scholar]
  63. Milne G, Graham J, McGrath J, Kirke R, McMaster W et al. Investigating farm fragmentation as a risk factor for bovine tuberculosis in cattle herds: a matched case-control study from Northern Ireland. Pathogens 2022; 11:299 [View Article] [PubMed]
    [Google Scholar]
  64. Madden JM, McGrath G, Sweeney J, Murray G, Tratalos JA et al. Spatio-temporal models of bovine tuberculosis in the Irish cattle population, 2012-2019. Spat Spatiotemporal Epidemiol 2021; 39:100441 [View Article] [PubMed]
    [Google Scholar]
  65. Brown E, Marshall AH, Mitchell HJ, Byrne AW. Cattle movements in Northern Ireland form a robust network: implications for disease management. Prev Vet Med 2019; 170:104740 [View Article] [PubMed]
    [Google Scholar]
  66. Malone FE, Wilson EC, Pollock JM, Skuce RA. Investigations into an outbreak of tuberculosis in a flock of sheep in contact with tuberculous cattle. J Vet Med Series B 2003; 50:500–504 [View Article]
    [Google Scholar]
  67. Muñoz-Mendoza M, Romero B, Del Cerro A, Gortázar C, García-Marín JF et al. Sheep as a potential source of bovine TB: epidemiology, pathology and evaluation of diagnostic techniques. Transbound Emerg Dis 2016; 63:635–646 [View Article] [PubMed]
    [Google Scholar]
  68. Lee J, Hanna R, Hill R, McCormick CM, Skuce RA. Bovine tuberculosis in an Eurasian otter. Vet Rec 2009; 164:727–728 [View Article] [PubMed]
    [Google Scholar]
  69. Allen AR. One bacillus to rule them all? – investigating broad range host adaptation in Mycobacterium bovis. Infection Genet Evol 2017; 53:68–76 [View Article]
    [Google Scholar]
  70. Zimpel CK, Patané JSL, Guedes ACP, de Souza RF, Silva-Pereira TT et al. Global distribution and evolution of Mycobacterium bovis lineages. Front Microbiol 2020; 11:843 [View Article]
    [Google Scholar]
  71. de la Rua-Domenech R, Goodchild AT, Vordermeier HM, Hewinson RG, Christiansen KH et al. Ante mortem diagnosis of tuberculosis in cattle: a review of the tuberculin tests, gamma-interferon assay and other ancillary diagnostic techniques. Res Vet Sci 2006; 81:190–210 [View Article] [PubMed]
    [Google Scholar]
  72. Bermingham ML, Handel IG, Glass EJ, Woolliams JA, de Clare Bronsvoort BM et al. Hui and Walter’s latent-class model extended to estimate diagnostic test properties from surveillance data: a latent model for latent data. Sci Rep 2015; 5:11861 [View Article] [PubMed]
    [Google Scholar]
  73. Zeineldin M, Camp P, Farrell D, Lehman K, Thacker T. Whole genome sequencing of Mycobacterium bovis directly from clinical tissue samples without culture. Front Microbiol 2023; 14:1141651 [View Article] [PubMed]
    [Google Scholar]
  74. Abernethy DA, Denny GO, Menzies FD, McGuckian P, Honhold N et al. The Northern Ireland programme for the control and eradication of Mycobacterium bovis. Vet Microbiol 2006; 112:231–237 [View Article] [PubMed]
    [Google Scholar]
  75. Rambaut A, Holmes EC, O’Toole Á, Hill V, McCrone JT et al. A dynamic nomenclature proposal for SARS-CoV-2 lineages to assist genomic epidemiology. Nat Microbiol 2020; 5:1403–1407 [View Article] [PubMed]
    [Google Scholar]
  76. O’Toole Á, Scher E, Underwood A, Jackson B, Hill V et al. Assignment of epidemiological lineages in an emerging pandemic using the pangolin tool. Virus Evol 2021; 7:veab064 [View Article] [PubMed]
    [Google Scholar]
  77. Zwyer M, Çavusoglu C, Ghielmetti G, Pacciarini ML, Scaltriti E et al. A new nomenclature for the livestock-associated Mycobacterium tuberculosis complex based on phylogenomics. Open Res Eur 2021; 1:100 [View Article] [PubMed]
    [Google Scholar]
http://instance.metastore.ingenta.com/content/journal/mgen/10.1099/mgen.0.001185
Loading
/content/journal/mgen/10.1099/mgen.0.001185
Loading

Data & Media loading...

Supplements

Supplementary material 1

PDF

Supplementary material 2

EXCEL
This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error