Population Structure of Mycobacterium bovis in Germany: a Long-Term Study Using Whole-Genome Sequencing Combined with Conventional Molecular Typing Methods

Inside Mycobacterium bovis SB0120 spoligotype circulating in Italy: analysis of the most frequent genotypes by whole genome sequencing

Bovine tuberculosis (bTB) is a chronic inflammatory disease primarily caused by Mycobacterium bovis. The infection affects domestic animals and wildlife, posing a zoonotic risk to humans. To understand the dynamics of transmission and genetic diversity in Italy's M. bovis population, we conducted whole-genome sequencing (WGS) analysis on two prevalent genotypes, belonging to Spoligotype SB0120, identified in different geographical and temporal contexts. By comparing these genomes with international M. bovis isolates, we identified a distinct clade within the lineage La1.2, encompassing the Italian SB0120 isolates, indicating a genomic segregation of Italian M. bovis from other European isolates. Within Italy, a significant level of genetic variability emerged across regions, while isolates within epidemiologically linked outbreaks exhibited minimal genetic diversity. Additionally, isolates derived from cattle and wild boars within a tuberculosis hotspot in Central Italy and from cattle and black pigs in Sicily formed unified clonal clusters. This indicates the presence of persistent strains circulating in the examined regions. The genetic diversity within herds was limited, as specific clones endured over time within certain herds. This research enhances our comprehension of the epidemiology and transmission patterns of bTB in Italy, thereby aiding the development of precise control strategies and disease management. Using WGS and implementing standardized protocols and databases will be pivotal in combating bTB and promoting One-Health approaches to address this noteworthy public health concern.

First Insight into the Whole Genome Sequencing Whole Variations in Mycobacterium bovis from Cattle in Morocco

Six cattle heads which tested positive against bovine tuberculosis (bTB) in Morocco were investigated to confirm the disease and to determine the source(s) of infection. Polymerase Chain Reaction (PCR) was directly performed on tissue samples collected from slaughtered animals. All investigated animals tested positive to PCR for the Mycobacterium bovis sub-type. Bacteriological isolation was conducted according to the technique recommended by WOAH for the cultivation of the Mycobacterium tuberculosis Complex (MBTC). Whole genome sequencing (WGS) was carried out on six mycobacterial isolates and the phylogenic tree was constructed. The six Moroccan isolates fit with clades II, III, IV, V and VII and were confirmed to belong to the clonal complexes Eu2, Unknown 2 and 7 as well as to sublineages La1.7.1, La1.2 and La1.8.2. The significant Single Nucleotide Polymorphism (SNPs) ranged from 84 to 117 between the isolates and the reference M. bovis strain and from 17 to 212 between the six isolates. Considering the high resolution of WGS, these results suggests that the source of infection of the bTB could be linked to imported animals as five of the investigated reactor animals were imported a few months prior. WGS can be a useful component to the Moroccan strategy to control bTB.

Rescue of Mycobacterium bovis DNA Obtained from Cultured Samples during Official Surveillance of Animal TB: Key Steps for Robust Whole Genome Sequence Data Generation

Epidemiological surveillance of animal tuberculosis (TB) based on whole genome sequencing (WGS) of Mycobacterium bovis has recently gained track due to its high resolution to identify infection sources, characterize the pathogen population structure, and facilitate contact tracing. However, the workflow from bacterial isolation to sequence data analysis has several technical challenges that may severely impact the power to understand the epidemiological scenario and inform outbreak response. While trying to use archived DNA from cultured samples obtained during routine official surveillance of animal TB in Portugal, we struggled against three major challenges: the low amount of M. bovis DNA obtained from routinely processed animal samples; the lack of purity of M. bovis DNA, i.e., high levels of contamination with DNA from other organisms; and the co-occurrence of more than one M. bovis strain per sample (within-host mixed infection). The loss of isolated genomes generates missed links in transmission chain reconstruction, hampering the biological and epidemiological interpretation of data as a whole. Upon identification of these challenges, we implemented an integrated solution framework based on whole genome amplification and a dedicated computational pipeline to minimize their effects and recover as many genomes as possible. With the approaches described herein, we were able to recover 62 out of 100 samples that would have otherwise been lost. Based on these results, we discuss adjustments that should be made in official and research laboratories to facilitate the sequential implementation of bacteriological culture, PCR, downstream genomics, and computational-based methods. All of this in a time frame supporting data-driven intervention.

Whole genome sequencing of Mycobacterium bovis directly from clinical tissue samples without culture

Advancement in next generation sequencing offers the possibility of routine use of whole genome sequencing (WGS) for Mycobacterium bovis (M. bovis) genomes in clinical reference laboratories. To date, the M. bovis genome could only be sequenced if the mycobacteria were cultured from tissue. This requirement for culture has been due to the overwhelmingly large amount of host DNA present when DNA is prepared directly from a granuloma. To overcome this formidable hurdle, we evaluated the usefulness of an RNA-based targeted enrichment method to sequence M. bovis DNA directly from tissue samples without culture. Initial spiking experiments for method development were established by spiking DNA extracted from tissue samples with serially diluted M. bovis BCG DNA at the following concentration range: 0.1 ng/μl to 0.1 pg/μl (10^-1 to 10^-4). Library preparation, hybridization and enrichment was performed using SureSelect custom capture library RNA baits and the SureSelect XT HS2 target enrichment system for Illumina paired-end sequencing. The method validation was then assessed using direct WGS of M. bovis DNA extracted from tissue samples from naturally (n = 6) and experimentally (n = 6) infected animals with variable Ct values. Direct WGS of spiked DNA samples achieved 99.1% mean genome coverage (mean depth of coverage: 108×) and 98.8% mean genome coverage (mean depth of coverage: 26.4×) for tissue samples spiked with BCG DNA at 10^-1 (mean Ct value: 20.3) and 10^-2 (mean Ct value: 23.4), respectively. The M. bovis genome from the experimentally and naturally infected tissue samples was successfully sequenced with a mean genome coverage of 99.56% and depth of genome coverage ranging from 9.2× to 72.1×. The spoligoyping and M. bovis group assignment derived from sequencing DNA directly from the infected tissue samples matched that of the cultured isolates from the same sample. Our results show that direct sequencing of M. bovis DNA from tissue samples has the potential to provide accurate sequencing of M. bovis genomes significantly faster than WGS from cultures in research and diagnostic settings.

Features of Mycobacterium bovis Complete Genomes Belonging to 5 Different Lineages

Mammalian tuberculosis (TB) is a zoonotic disease mainly due to Mycobacterium bovis (M. bovis). A current challenge for its eradication is understanding its transmission within multi-host systems. Improvements in long-read sequencing technologies have made it possible to obtain complete bacterial genomes that provide a comprehensive view of species-specific genomic features. In the context of TB, new genomic references based on complete genomes genetically close to field strains are also essential to perform precise field molecular epidemiological studies. A total of 10 M. bovis strains representing each genetic lineage identified in France and in other countries were selected for performing complete assembly of their genomes. Pangenome analysis revealed a "closed" pangenome composed of 3900 core genes and only 96 accessory genes. Whole genomes-based alignment using progressive Mauve showed remarkable conservation of the genomic synteny except that the genomes have a variable number of copies of IS6110. Characteristic genomic traits of each lineage were identified through the discovery of specific indels. Altogether, these results provide new genetic features that improve the description of M. bovis lineages. The availability of new complete representative genomes of M. bovis will be useful to epidemiological studies and better understand the transmission of this clonal-evolving pathogen.

Performance and Agreement Between WGS Variant Calling Pipelines Used for Bovine Tuberculosis Control: Toward International Standardization

Whole genome sequencing (WGS) and allied variant calling pipelines are a valuable tool for the control and eradication of infectious diseases, since they allow the assessment of the genetic relatedness of strains of animal pathogens. In the context of the control of tuberculosis (TB) in livestock, mainly caused by Mycobacterium bovis, these tools offer a high-resolution alternative to traditional molecular methods in the study of herd breakdown events. However, despite the increased use and efforts in the standardization of WGS methods in human tuberculosis around the world, the application of these WGS-enabled approaches to control TB in livestock is still in early development. Our study pursued an initial evaluation of the performance and agreement of four publicly available pipelines for the analysis of M. bovis WGS data (vSNP, SNiPgenie, BovTB, and MTBseq) on a set of simulated Illumina reads generated from a real-world setting with high TB prevalence in cattle and wildlife in the Republic of Ireland. The overall performance of the evaluated pipelines was high, with recall and precision rates above 99% once repeat-rich and problematic regions were removed from the analyses. In addition, when the same filters were applied, distances between inferred phylogenetic trees were similar and pairwise comparison revealed that most of the differences were due to the positioning of polytomies. Hence, under the studied conditions, all pipelines offer similar performance for variant calling to underpin real-world studies of M. bovis transmission dynamics.

A new nomenclature for the livestock-associated Mycobacterium tuberculosis complex based on phylogenomics

Background:  The bacteria that compose the Mycobacterium tuberculosis complex (MTBC) cause tuberculosis (TB) in humans and in different animals, including livestock. Much progress has been made in understanding the population structure of the human-adapted members of the MTBC by combining phylogenetics with genomics. Accompanying the discovery of new genetic diversity, a body of operational nomenclature has evolved to assist comparative and molecular epidemiological studies of human TB. By contrast, for the livestock-associated MTBC members, Mycobacterium bovis, M. caprae and M. orygis, there has been a lack of comprehensive nomenclature to accommodate new genetic diversity uncovered by emerging phylogenomic studies. We propose to fill this gap by putting forward a new nomenclature covering the main phylogenetic groups within M. bovis, M. caprae and M. orygis. Methods:  We gathered a total of 8,736 whole-genome sequences (WGS) from public sources and 39 newly sequenced strains, and selected a subset of 829 WGS, representative of the worldwide diversity of M. bovis, M. caprae and M. orygis. We used phylogenetics and genetic diversity patterns inferred from WGS to define groups. Results:  We propose to divide M. bovis, M. caprae and M. orygis in three main phylogenetic lineages, which we named La1, La2 and La3, respectively. Within La1, we identified several monophyletic groups, which we propose to classify into eight sublineages (La1.1-La1.8). These sublineages differed in geographic distribution, with some being geographically restricted and others globally widespread, suggesting different expansion abilities. To ease molecular characterization of these MTBC groups by the community, we provide phylogenetically informed, single nucleotide polymorphisms that can be used as barcodes for genotyping. These markers were implemented in KvarQ and TB-Profiler, which are platform-independent, open-source tools. Conclusions:  Our results contribute to an improved classification of the genetic diversity within the livestock-associated MTBC, which will benefit future molecular epidemiological and evolutionary studies.

Inferring Mycobacterium bovis transmission between cattle and badgers using isolates from the Randomised Badger Culling Trial

Mycobacterium bovis (M. bovis) is a causative agent of bovine tuberculosis, a significant source of morbidity and mortality in the global cattle industry. The Randomised Badger Culling Trial was a field experiment carried out between 1998 and 2005 in the South West of England. As part of this trial, M. bovis isolates were collected from contemporaneous and overlapping populations of badgers and cattle within ten defined trial areas. We combined whole genome sequences from 1,442 isolates with location and cattle movement data, identifying transmission clusters and inferred rates and routes of transmission of M. bovis. Most trial areas contained a single transmission cluster that had been established shortly before sampling, often contemporaneous with the expansion of bovine tuberculosis in the 1980s. The estimated rate of transmission from badger to cattle was approximately two times higher than from cattle to badger, and the rate of within-species transmission considerably exceeded these for both species. We identified long distance transmission events linked to cattle movement, recurrence of herd breakdown by infection within the same transmission clusters and superspreader events driven by cattle but not badgers. Overall, our data suggests that the transmission clusters in different parts of South West England that are still evident today were established by long-distance seeding events involving cattle movement, not by recrudescence from a long-established wildlife reservoir. Clusters are maintained primarily by within-species transmission, with less frequent spill-over both from badger to cattle and cattle to badger.

Whole-Genome SNP Analysis Identifies Putative Mycobacterium bovis Transmission Clusters in Livestock and Wildlife in Catalonia, Spain

The high-resolution WGS analyses of MTBC strains have provided useful insight for determining sources of infection for animal tuberculosis. In Spain, tuberculosis in livestock is caused by Mycobacterium bovis and Mycobacterium caprae, where wildlife reservoirs play an important role. We analyzed a set of 125 M. bovis isolates obtained from livestock and wildlife from Catalonia to investigate strain diversity and identify possible sources and/or causes of infection. Whole-genome SNP profiles were used for phylogenetic reconstruction and pairwise SNP distance analysis. Additionally, SNPs were investigated to identify virulence and antimicrobial resistance factors to investigate clade-specific associations. Putative transmission clusters (≤12 SNPs) were identified, and associated epidemiological metadata were used to determine possible explanatory factors for transmission. M. bovis distribution was heterogeneous, with 7 major clades and 21 putative transmission clusters. In order of importance, the explanatory factors associated were proximity and neighborhood, residual infection, livestock-wildlife interaction, shared pasture, and movement. Genes related to lipid transport and metabolism showed the highest number of SNPs. All isolates were pyrazinamide resistant, and five were additionally resistant to isoniazid, but no clade-specific associations could be determined. Our findings highlight the importance of high-resolution molecular surveillance to monitor bovine tuberculosis dynamics in a low-prevalence setting.