European 1: A globally important clonal complex of Mycobacterium bovis

Tuberculosis in Wild Pigs from Argentina

Mycobacterium bovis, a member of the Mycobacterium tuberculosis complex (MTC) and non-tuberculous Mycobacteria (NTM), may infect wild and domestic mammals, including humans. Although cattle are the main hosts and spreaders of M. bovis, many wildlife hosts play an important role worldwide. In Argentina, wild boar and domestic pigs are considered important links in mammalian tuberculosis (mTB) transmission. The aim of this work was to investigate the presence of M. bovis in wild pigs from different regions of Argentina, to characterize isolates of M. bovis obtained, and to compare those with other previously found in vertebrate hosts. A total of 311 samples from wild pigs were obtained, and bacteriological culture, molecular identification and genotyping were performed, obtaining 63 isolates (34 MTC and 29 NTM). Twelve M. bovis spoligotypes were detected. Our findings suggest that wild pigs have a prominent role as reservoirs of mTB in Argentina, based on an estimated prevalence of 11.2 ± 1.8% (95% CI 8.0-14.8) for MTC and the frequency distribution of spoligotypes shared by cattle (75%), domestic pigs (58%) and wildlife (50%). Argentina has a typical scenario where cattle and pigs are farm-raised extensively, sharing the environment with wildlife, creating conditions for effective transmission of mTB in the wildlife-livestock-human interface.

Insight into the genetic diversity of Mycobacterium bovis isolated from cattle in Malawi

We describe the genetic diversity and phylogenetic relationships of Mycobacterium bovis, isolated from cattle in Malawi. Deletion analysis, spoligotyping, and MIRU-VNTR typing were used to genotype the isolates. Combined with a larger dataset from neighboring countries, the overall M. bovis diversity in Southern Africa was contextualized. From the southern and northern regions of Malawi, 24 isolates were confirmed as M. bovis. We pooled data for the central region (60 isolates) from our recent publication to conceptualize the genetic and phylogenetic relationships of M. bovis in Malawi. European 1 was the dominant M. bovis clonal complex, with 10 unique spoligotype patterns, and SB0131 was ubiquitous. High genetic diversity, a low clustering rate, and many singletons, coupled with a low mutation transmission index, infer a low level of recent transmission, and suggest an endemic status of bovine tuberculosis (bTB) in Malawi. M. bovis isolates from Zambia, Mozambique, and South Africa were genetically related to Malawian isolates, whereas Tanzanian isolates were distantly related. The diversity and phylogenetic analysis suggest earlier introductions and maintenance of M. bovis by constant reinfection from reservoir animals. These findings are fundamental to understanding the source and route of infection in order to establish alternative management strategies for bTB.

Deciphering the evolution of the temporal and geographic distribution of French Mycobacterium bovis genotypes using a high throughput SNP-targeted amplicon sequencing method

Mycobacterium bovis, which belongs to the Mycobacterium tuberculosis complex, is a highly clonal pathogen. However, several lineages of M. bovis have been described worldwide and nine different clusters were identified in France. Targeted amplicon sequencing using next-generation sequencing technology of eighty-eight phylogenetically informative single nucleotide polymorphisms (SNPs) were used to infer the phylogenetic relationship of 630 strains of the National Reference Laboratory isolated between 1979 and 2018 from various animal species. This study allowed classifying 618 different genotypic profiles (combination of a spoligotype and 8 loci-MIRU-VNTR profiles) into the nine previously identified clusters. A global analysis of the entire collection of the National Reference Laboratory has made it possible to represent the evolution of clonal complexes and clusters in time and space for better assessing epidemiological changes of bovine tuberculosis in France.

The devil you know and the devil you don't: current status and challenges of bovine tuberculosis eradication in the United States

Having entered into its second century, the eradication program for bovine tuberculosis (bTB, caused by Mycobacterium bovis) in the United States of America occupies a position both enviable and daunting. Excepting four counties in Michigan comprising only 6109 km2 (0.06% of US land area) classified as Modified Accredited, as of April 2022 the entire country was considered Accredited Free of bTB by the US Department of Agriculture for cattle and bison. On the surface, the now well-described circumstances of endemic bTB in Michigan, where white-tailed deer (Odocoileus virginianus) serve as a free-ranging wildlife maintenance host, may appear to be the principal remaining barrier to national eradication. However, the situation there is unique in the U.S., and far-removed from the broader issues of bTB control in the remainder of the country. In Michigan, extensive surveillance for bTB in deer over the last quarter century, and regulatory measures to maximize the harvest of publicly-owned wildlife, have been implemented and sustained. Prevalence of bTB in deer has remained at a low level, although not sufficiently low to eliminate cattle herd infections. Public attitudes towards bTB, cattle and deer, and their relative importance, have been more influential in the management of the disease than any limitations of biological science. However, profound changes in the demographics and social attitudes of Michigan's human population are underway, changes which are likely to force a critical reevaluation of the bTB control strategies thus far considered integral. In the rest of the U.S. where bTB is not self-sustaining in wildlife, changes in the scale of cattle production, coupled with both technical and non-technical issues have created their own substantial challenges. It is against this diverse backdrop that the evolution of whole genome sequencing of M. bovis has revolutionized understanding of the history and ecology of bTB in Michigan, resolved previously undiscernible epidemiological puzzles, provided insights into zoonotic transmission, and unified eradication efforts across species and agencies. We describe the current status of bTB eradication in the U.S., how circumstances and management have changed, what has been learned, and what remains more elusive than ever.

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.

Whole genome sequencing of Ethiopian Brucella abortus isolates expands the known diversity of an early branching sub-Saharan African lineage

Brucellosis remains one of the most significant zoonotic diseases globally, responsible for both considerable human morbidity and economic losses due to its impacts on livestock productivity. Despite this, there remain significant evidence gaps in many low- and middle-income countries, including those of sub-Saharan Africa. Here we report the first molecular characterisation of Brucella sp. from Ethiopia. Fifteen Brucella sp. isolates from an outbreak in cattle from a herd in central Ethiopia were identified as Brucella abortus, using bacterial culture and molecular methods. Sequencing of the Ethiopian B. abortus isolates allowed their phylogenetic comparison with 411 B. abortus strains of diverse geographical origins, using whole genome single nucleotide polymorphisms (wgSNP). The Ethiopian isolates belonged to an early-branching lineage (Lineage A) previously only represented by data from two strains, both of sub-Saharan African origin (Kenya and Mozambique). A second B. abortus lineage (Lineage B), also comprised solely of strains originating from sub-Saharan Africa, was identified. The majority of strains belonged to one of two lineages of strains originating from a much broader geographical range. Further analyses based on multi-locus sequence typing (MLST) and multi-locus variable-number tandem repeat analysis (MLVA) expanded the number of B. abortus strains available for comparison with the Ethiopian isolates and were consistent with the findings from wgSNP analysis. MLST profiles of the Ethiopian isolates expanded the sequence type (ST) diversity of the early branching lineage of B. abortus, equivalent to wgSNP Lineage A. A more diverse cluster of STs, equivalent to wgSNP Lineage B, was comprised solely of strains originating from sub-Saharan Africa. Similarly, analysis of B. abortus MLVA profiles (n = 1891) confirmed that the Ethiopian isolates formed a unique cluster, similar to only two existing strains, and distinct from the majority of other strains of sub-Saharan African origin. These findings expand the known diversity of an under-represented lineage of B. abortus and suggest a potential evolutionary origin for the species in East Africa. In addition to providing information concerning Brucella species extant within Ethiopia this work serves as the basis for further studies on the global population structure and evolutionary history of a major zoonotic pathogen.

Selective sweep sites and SNP dense regions differentiate Mycobacterium bovis isolates across scales

Mycobacterium bovis, a bacterial zoonotic pathogen responsible for the economically and agriculturally important livestock disease bovine tuberculosis (bTB), infects a broad mammalian host range worldwide. This characteristic has led to bidirectional transmission events between livestock and wildlife species as well as the formation of wildlife reservoirs, impacting the success of bTB control measures. Next Generation Sequencing (NGS) has transformed our ability to understand disease transmission events by tracking variant sites, however the genomic signatures related to host adaptation following spillover, alongside the role of other genomic factors in the M. bovis transmission process are understudied problems. We analyzed publicly available M. bovis datasets collected from 700 hosts across three countries with bTB endemic regions (United Kingdom, United States, and New Zealand) to investigate if genomic regions with high SNP density and/or selective sweep sites play a role in Mycobacterium bovis adaptation to new environments (e.g., at the host-species, geographical, and/or sub-population levels). A simulated M. bovis alignment was created to generate null distributions for defining genomic regions with high SNP counts and regions with selective sweeps evidence. Random Forest (RF) models were used to investigate evolutionary metrics within the genomic regions of interest to determine which genomic processes were the best for classifying M. bovis across ecological scales. We identified in the M. bovis genomes 14 and 132 high SNP density and selective sweep regions, respectively. Selective sweep regions were ranked as the most important in classifying M. bovis across the different scales in all RF models. SNP dense regions were found to have high importance in the badger and cattle specific RF models in classifying badger derived isolates from livestock derived ones. Additionally, the genes detected within these genomic regions harbor various pathogenic functions such as virulence and immunogenicity, membrane structure, host survival, and mycobactin production. The results of this study demonstrate how comparative genomics alongside machine learning approaches are useful to investigate further the nature of M. bovis host-pathogen interactions.

Molecular epidemiology of bovine tuberculosis in Northern Ghana identifies several uncharacterized bovine spoligotypes and suggests possible zoonotic transmission

Objective

We conducted an abattoir-based cross-sectional study in the five administrative regions of Northern Ghana to determine the distribution of bovine tuberculosis (BTB) among slaughtered carcasses and identify the possibility of zoonotic transmission.

Methods

Direct smear microscopy was done on 438 tuberculosis-like lesions from selected cattle organs and cultured on Lowenstein-Jensen media. Acid-fast bacilli (AFB) isolates were confirmed as members of the Mycobacterium tuberculosis complex (MTBC) by PCR amplification of IS6110 and rpoß. Characterization and assignment into MTBC lineage and sub-lineage were done by spoligotyping, with the aid of the SITVIT2, miruvntrplus and mbovis.org databases. Spoligotype data was compared to that of clinical M. bovis isolates from the same regions to identify similarities.

Results

A total of 319/438 (72.8%) lesion homogenates were smear positive out of which, 84.6% (270/319) had microscopic grade of at least 1+ for AFB. Two hundred and sixty-five samples (265/438; 60.5%) were culture positive, of which 212 (80.0%) were MTBC. Approximately 16.7% (34/203) of the isolates with correctly defined spoligotypes were negative for IS6110 PCR but were confirmed by rpoß. Spoligotyping characterized 203 isolates as M. bovis (198, 97.5%), M. caprae (3, 1.5%), M. tuberculosis (Mtbss) lineage (L) 4 Cameroon sub-lineage, (1, 0.5%), and M. africanum (Maf) L6 (1, 0.5%). A total of 53 unique spoligotype patterns were identified across the five administrative regions (33 and 28 were identified as orphan respectively by the SITVIT2 and mbovis.org databases), with the most dominant spoligotype being SIT1037/ SB0944 (77/203, 37.93%). Analysis of the bovine and human M. bovis isolates showed 75% (3/4) human M. bovis isolates sharing the same spoligotype pattern with the bovine isolates.

Conclusion

Our study identified that approximately 29% of M. bovis strains causing BTB in Northern Ghana are caused by uncharacterized spoligotypes. Our findings suggest possible zoonotic transmission and highlight the need for BTB disease control in Northern Ghana.

Spatial distribution of Mycobacterium bovis spoligotypes in cattle from the state of Minas Gerais, Brazil

Genotyping methods have led to a better understanding of the epidemiology of Mycobacterium bovis (M. bovis) infection, and its transmission dynamics, as well as the possible phylogenetic relationships between Mycobacterium strains, thus making bovine tuberculosis control programs more efficient. The goal of this study was to characterize the main spoligotypes of M. bovis isolated from cattle in the State of Minas Gerais, Brazil. It was carried out in 28 municipalities of "Triângulo Mineiro" and "Alto Paranaíba" regions of the state. Viscera samples were obtained from 58 bovines positive for tuberculosis according to comparative cervical tests, and from another 100 bovines with lesions suggestive of tuberculosis, which were donated by the National Agricultural Laboratory of Pedro Leopoldo, Minas Gerais. Microbiological isolation was performed in Stonebrink medium, and molecular identification of mycobacteria was performed by PCR. Genotyping was performed using the spoligotyping method at the Agrobiotechnology and Molecular Biology Institute of National Agricultural Technology Institute-National Scientific and Technical Research Council, Buenos Aires, Argentina. Among the 158 viscera samples, we obtained 40 (25%) isolates of M. bovis, and detected 11 spoligotype patterns, with a predominance of SB1142 (37.5%), SB0121 (25.0%), and SB1145 (10.0%). Other standards, SB0295, SB1050, SB0881, SB1144, SB1802, SB0140, SB0120, and SB0849, varied from 2.5 to 7.5%, heterogeneously distributed among the municipalities. The presence of spoligotypes shared with other Brazilian states and different countries indicates their possible exchange through epidemiological relationships, such as the transit of live animals and/or genetic similarity between strains that share a common ancestor.

IS6110 Copy Number in Multi-Host Mycobacterium bovis Strains Circulating in Bovine Tuberculosis Endemic French Regions

IS6110 is an insertion sequence found in the Mycobacterium tuberculosis complex, to which Mycobacterium bovis belongs, which can play a role in genome plasticity and in bacterial evolution. In this study, the abundance and location of IS6110 on M. bovis genomic data of French animal field strains were studied. A first analysis was performed on a panel of 81 strains that reflect the national M. bovis population’s genetic diversity. The results show that more than one-third of them are IS6110 multicopy and that 10% have IS6110 in a high copy number (more than 6 copies). Multicopy strains are those circulating in the regions where prevalence was above the national average. Further study of 93 such strains, with an IS6110 copy number of 10-12, showed stability of IS6110 copy number and genome location over time and between host species. The correlation between M. bovis multicopy strains and high bovine tuberculosis (bTB) prevalence leads us to consider whether their epidemiological success could be partly due to genetic changes originated by IS6110 transposition.

A Bayesian evolutionary model towards understanding wildlife contribution to F4-family Mycobacterium bovis transmission in the South-West of France

In two “départements” in the South-West of France, bovine tuberculosis (bTB) outbreaks due to Mycobacterium bovis spoligotype SB0821 have been identified in cattle since 2002 and in wildlife since 2013. Using whole genome sequencing, the aim of our study was to clarify badger contribution to bTB transmission in this area. We used a Bayesian evolutionary model, to infer phylogenetic trees and migration rates between two pathogen populations defined by their host-species. In order to account for sampling bias, sub-population structure was inferred using the marginal approximation of the structured coalescent (Mascot) implemented in BEAST2. We included 167 SB0821 strains (21 isolated from badgers and 146 from cattle) and identified 171 single nucleotide polymorphisms. We selected a HKY model and a strict molecular clock. We estimated a badger-to-cattle transition rate (median: 2.2 transitions/lineage/year) 52 times superior to the cattle-to-badger rate (median: 0.042 transitions/lineage/year). Using the maximum clade credibility tree, we identified that over 75% of the lineages from 1989 to 2000 were present in badgers. In addition, we calculated a median of 64 transition events from badger-to-cattle (IQR: 10–91) and a median of zero transition event from cattle-to-badger (IQR: 0–3). Our model enabled us to infer inter-species transitions but not intra-population transmission as in previous epidemiological studies, where relevant units were farms and badger social groups. Thus, while we could not confirm badgers as possible intermediaries in farm-to-farm transmission, badger-to-cattle transition rate was high and we confirmed long-term presence of M.bovis in the badger population in the South-West of France.

SITVITBovis—a publicly available database and mapping tool to get an improved overview of animal and human cases caused by Mycobacterium bovis

Limited data are available for bovine tuberculosis and the infections it can cause in humans and other mammals. We therefore constructed a publicly accessible SITVITBovis database that incorporates genotyping and epidemiological data on Mycobacterium bovis. It also includes limited data on Mycobacterium caprae (previously synonymous with the name M. bovis subsp. Caprae) that can infect both animals and humans. SITVITBovis incorporates data on 25,741 isolates corresponding to 60 countries of origin (75 countries of isolation). It reports a total of 1000 spoligotype patterns: 537 spoligotype international types (SITs, containing 25 278 clinical isolates) and 463 orphan patterns, allowing a wide overview of the geographic distribution of various phylogenetical sublineages (BOV_1, BOV_2, BOV_3 and BOV_4-CAPRAE). The SIT identifiers of the SITVITBovis were compared to the SB numbers of the Mbovis.org database to facilitate crosscheck among databases. Note that SITVITBovis also contains limited information on mycobacterial interspersed repetitive units-variable number of tandem repeats when available. Significant differences were observed when comparing age/gender of human isolates as well as various hosts. The database includes information on the regions where a strain was isolated as well as hosts involved, making it possible to see geographic trends. SITVITBovis is publicly accessible at: http://www.pasteur-guadeloupe.fr:8081/SITVIT_Bovis. Finally, a future second version is currently in progress to allow query of associated whole-genome sequencing data.

Database URLhttp://www.pasteur-guadeloupe.fr:8081/SITVIT_Bovis

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.

The open pan-genome architecture and virulence landscape of Mycobacterium bovis

Animal tuberculosis (TB) is an emergent disease caused by Mycobacterium bovis, one of the animal-adapted ecotypes of the Mycobacterium tuberculosis complex (MTC). In this work, whole-genome comparative analyses of 70 M. bovis were performed to gain insights into the pan-genome architecture. The comparison across M. bovis predicted genome composition enabled clustering into the core- and accessory-genome components, with 2736 CDS for the former, while the accessory moiety included 3897 CDS, of which 2656 are restricted to one/two genomes only. These analyses predicted an open pan-genome architecture, with an average of 32 CDS added by each genome and show the diversification of discrete M. bovis subpopulations supported by both core- and accessory-genome components. The functional annotation of the pan-genome classified each CDS into one or several COG (Clusters of Orthologous Groups) categories, revealing ‘transcription’ (total average CDSs, n=258), ‘lipid metabolism and transport’ (n=242), ‘energy production and conversion’ (n=214) and ‘unknown function’ (n=876) as the most represented. The closer analysis of polymorphisms in virulence-related genes in a restrict group of M. bovis from a multi-host system enabled the identification of clade-monomorphic non-synonymous SNPs, illustrating clade-specific virulence landscapes and correlating with disease severity. This first comparative pan-genome study of a diverse collection of M. bovis encompassing all clonal complexes indicates a high percentage of accessory genes and denotes an open, dynamic non-conservative pan-genome structure, with high evolutionary potential, defying the canons of MTC biology. Furthermore, it shows that M. bovis can shape its virulence repertoire, either by acquisition and loss of genes or by SNP-based diversification, likely towards host immune evasion, adaptation and persistence.

Genomic analysis of an outbreak of bovine tuberculosis in a man-made multi-host species system: A call for action on wildlife in Brazil

We report on a 15-year-long outbreak of bovine tuberculosis (bTB) in wildlife from a Brazilian safari park. A timeline of diagnostic events and whole-genome sequencing (WGS) of 21 Mycobacterium bovis isolates from deer and llamas were analyzed. Accordingly, from 2003 to 2018, at least 16 animals, from eight species, died due to TB, which is likely an underestimated number. In three occasions since 2013, the deer presented positive tuberculin tests, leading to the park closure and culling of all deer. WGS indicated that multiple M. bovis strains were circulating, with at least three founding introductions since the park inauguration in 1977. Using a previously sequenced dataset of 71 M. bovis genomes from cattle, we found no recent transmission events between nearby farms and the park based on WGS. Lastly, by discussing socio-economic and environmental factors escaping current regulatory gaps that were determinant of this outbreak, we pledge for the development of a plan to report and control bTB in wildlife in Brazil.

Genome-wide estimation of recombination, mutation and positive selection enlightens diversification drivers of Mycobacterium bovis

Genome sequencing has reinvigorated the infectious disease research field, shedding light on disease epidemiology, pathogenesis, host-pathogen interactions and also evolutionary processes exerted upon pathogens. Mycobacterium tuberculosis complex (MTBC), enclosing M. bovis as one of its animal-adapted members causing tuberculosis (TB) in terrestrial mammals, is a paradigmatic model of bacterial evolution. As other MTBC members, M. bovis is postulated as a strictly clonal, slowly evolving pathogen, with apparently no signs of recombination or horizontal gene transfer. In this work, we applied comparative genomics to a whole genome sequence (WGS) dataset composed by 70 M. bovis from different lineages (European and African) to gain insights into the evolutionary forces that shape genetic diversification in M. bovis. Three distinct approaches were used to estimate signs of recombination. Globally, a small number of recombinant events was identified and confirmed by two independent methods with solid support. Still, recombination reveals a weaker effect on M. bovis diversity compared with mutation (overall r/m = 0.037). The differential r/m average values obtained across the clonal complexes of M. bovis in our dataset are consistent with the general notion that the extent of recombination may vary widely among lineages assigned to the same taxonomical species. Based on this work, recombination in M. bovis cannot be excluded and should thus be a topic of further effort in future comparative genomics studies for which WGS of large datasets from different epidemiological scenarios across the world is crucial. A smaller M. bovis dataset (n = 42) from a multi-host TB endemic scenario was then subjected to additional analyses, with the identification of more than 1,800 sites wherein at least one strain showed a single nucleotide polymorphism (SNP). The majority (87.1%) was located in coding regions, with the global ratio of non-synonymous upon synonymous alterations (dN/dS) exceeding 1.5, suggesting that positive selection is an important evolutionary force exerted upon M. bovis. A higher percentage of SNPs was detected in genes enriched into "lipid metabolism", "cell wall and cell processes" and "intermediary metabolism and respiration" functional categories, revealing their underlying importance in M. bovis biology and evolution. A closer look on genes prone to horizontal gene transfer in the MTBC ancestor and included in the 3R (DNA repair, replication and recombination) system revealed a global average negative value for Taijima's D neutrality test, suggesting that past selective sweeps and population expansion after a recent bottleneck remain as major evolutionary drivers of the obligatory pathogen M. bovis in its struggle with the host.

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.

Mycobacterial Infection of Precision-Cut Lung Slices Reveals Type 1 Interferon Pathway Is Locally Induced by Mycobacterium bovis but Not M. tuberculosis in a Cattle Breed

Tuberculosis exacts a terrible toll on human and animal health. While Mycobacterium tuberculosis (Mtb) is restricted to humans, Mycobacterium bovis (Mb) is present in a large range of mammalian hosts. In cattle, bovine TB (bTB) is a noticeable disease responsible for important economic losses in developed countries and underestimated zoonosis in the developing world. Early interactions that take place between mycobacteria and the lung tissue early after aerosol infection govern the outcome of the disease. In cattle, these early steps remain poorly characterized. The precision-cut lung slice (PCLS) model preserves the structure and cell diversity of the lung. We developed this model in cattle in order to study the early lung response to mycobacterial infection. In situ imaging of PCLS infected with fluorescent Mb revealed bacilli in the alveolar compartment, in adjacent or inside alveolar macrophages, and in close contact with pneumocytes. We analyzed the global transcriptional lung inflammation signature following infection of PCLS with Mb and Mtb in two French beef breeds: Blonde d'Aquitaine and Charolaise. Whereas, lungs from the Blonde d'Aquitaine produced high levels of mediators of neutrophil and monocyte recruitment in response to infection, such signatures were not observed in the Charolaise in our study. In the Blonde d'Aquitaine lung, whereas the inflammatory response was highly induced by two Mb strains, AF2122 isolated from cattle in the UK and Mb3601 circulating in France, the response against two Mtb strains, H37Rv, the reference laboratory strain, and BTB1558, isolated from zebu in Ethiopia, was very low. Strikingly, the type I interferon pathway was only induced by Mb but not Mtb strains, indicating that this pathway may be involved in mycobacterial virulence and host tropism. Hence, the PCLS model in cattle is a valuable tool to deepen our understanding of early interactions between lung host cells and mycobacteria. It revealed striking differences between cattle breeds and mycobacterial strains. This model could help in deciphering biomarkers of resistance vs. susceptibility to bTB in cattle as such information is still critically needed for bovine genetic selection programs and would greatly help the global effort to eradicate bTB.

Genomic and temporal analyses of Mycobacterium bovis in southern Brazil

Mycobacterium bovis is a causal agent of bovine tuberculosis (bTB), one of the most important diseases currently facing the cattle industry worldwide. Tracing the source of M. bovis infections of livestock is an important tool for understanding the epidemiology of bTB and defining control/eradication strategies. In this study, whole genome sequencing (WGS) of 74 M. bovis isolates sourced from naturally infected cattle in the State of Rio Grande do Sul (RS), southern Brazil, was used to evaluate the population structure of M. bovis in the region, identify potential transmission events and date the introduction of clonal complex (CC) European 2 (Eu2). In silico spoligotyping identified 11 distinct patterns including four new profiles and two CCs, European 1 (Eu1) and Eu2. The analyses revealed a high level of genetic diversity in the majority of herds and identified putative transmission clusters that suggested that within- and between-herd transmission is occurring in RS. In addition, a comparison with other published M. bovis isolates from Argentina, Brazil, Paraguay and Uruguay demonstrated some evidence for a possible cross-border transmission of CC Eu1 into RS from Uruguay or Argentina. An estimated date for the introduction of CC Eu2 into RS in the middle of the 19th century correlated with the historical introduction of cattle into RS to improve existing local breeds. These findings contribute to the understanding of the population structure of M. bovis in southern Brazil and highlight the potential of WGS in surveillance and helping to identify bTB transmission.

Population structure and transmission of Mycobacterium bovis in Ethiopia

Bovine tuberculosis (bTB) is endemic in cattle in Ethiopia, a country that hosts the largest national cattle herd in Africa. The intensive dairy sector, most of which is peri-urban, has the highest prevalence of disease. Previous studies in Ethiopia have demonstrated that the main cause is Mycobacterium bovis, which has been investigated using conventional molecular tools including deletion typing, spoligotyping and Mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR). Here we use whole-genome sequencing to examine the population structure of M. bovis in Ethiopia. A total of 134 M. bovis isolates were sequenced including 128 genomes from 85 mainly dairy cattle and six genomes isolated from humans, originating from 12 study sites across Ethiopia. These genomes provided a good representation of the previously described population structure of M. bovis, based on spoligotyping and demonstrated that the population is dominated by the clonal complexes African 2 (Af2) and European 3 (Eu3). A range of within-host diversity was observed amongst the isolates and evidence was found for both short- and long-distance transmission. Detailed analysis of available genomes from the Eu3 clonal complex combined with previously published genomes revealed two distinct introductions of this clonal complex into Ethiopia between 1950 and 1987, likely from Europe. This work is important to help better understand bTB transmission in cattle in Ethiopia and can potentially inform national strategies for bTB control in Ethiopia and beyond.

Molecular epidemiology of Mycobacterium bovis in central parts of Malawi

Bovine tuberculosis (bTB) is a neglected disease that affects cattle and humans. The burden of bTB is higher in developing countries as compared to industrialized countries. The reasons behind this discrepancy include the fact that bTB control measures, such as testing and slaughter of infected cattle and pasteurization of milk, are not usually practised in developing countries largely because of their high cost. To improve our understanding of bTB in developing countries, molecular typing studies are essential, in particular in terms of transmission dynamics, infection sources and knowledge of circulating strains of the principal causative agent, Mycobacterium bovis. In this study, we applied a suite of molecular typing techniques encompassing deletion analysis, spoligotyping and MIRU-VNTR to isolates recovered from samples collected during the routine post-mortem of cattle at the cold storage abattoir in Lilongwe, Malawi. Out of 63 isolates, 51 (81%) belonged to the European 1. M. bovis clonal complex. Spoligotyping identified 8 profiles, with SB0131 being the predominant type (56% of isolates). Spoligotypes SB0273 and SB0425 were identified in 14% and 13%, respectively, of the isolates. MIRU-VNTR showed a high discriminatory power of 0.959 and differentiated the 8 spoligotypes to 31 genotypes. The high diversity of M. bovis within the study area suggests the infection has been circulating in the area for a considerable period of time, likely facilitated by the lack of effective control measures. We also observed genetic similarities between isolates from Malawi (this study) to isolates described in previous studies in Zambia and Mozambique, suggesting transmission links in this region. The information provided by this study provides much needed evidence for the formulation of improved bTB control strategies.

Genetic Diversity and Potential Paths of Transmission of Mycobacterium bovis in the Amazon: The Discovery of M. bovis Lineage Lb1 Circulating in South America

Bovine tuberculosis (bTB) has yet to be eradicated in Brazil. Herds of cattle and buffalo are important sources of revenue to people living in the banks of the Amazon River basin. A better understanding of Mycobacterium bovis (M. bovis) populational structure and transmission dynamics affecting these animals can significantly contribute in efforts to improve their sanitary status. Herein, we sequenced the whole genome of 22 M. bovis isolates (15 from buffalo and 7 from cattle) from 10 municipalities in the region of the Lower Amazon River Basin in Brazil and performed phylogenomic analysis and Single Nucleotide Polymorphism (SNP)-based transmission inference to evaluate population structure and transmission networks. Additionally, we compared these genomes to others obtained in unrelated studies in the Marajó Island (n = 15) and worldwide (n = 128) to understand strain diversity in the Amazon and to infer M. bovis lineages. Our results show a higher genomic diversity of M. bovis genomes obtained in the Lower Amazon River region when compared to the Marajó Island, while no significant difference was observed between M. bovis genomes obtained from cattle and buffalo (p ≥ 0.05). This high genetic diversity is reflected by the weak phylogenetic clustering of M. bovis from the Lower Amazon River region based on geographic proximity and in the detection of only two putative transmission clusters in the region. One of these clusters is the first description of inter-species transmission between cattle and buffalo in the Amazon, bringing implications to the bTB control program. Surprisingly, two M. bovis lineages were detected in our dataset, namely Lb1 and Lb3, constituting the first description of Lb1 in South America. Most of the strains of this study (13/22) and all 15 strains of the Marajó Island carried no clonal complex marker, suggesting that the recent lineage classification better describe the diversity of M. bovis in the Amazon.

Long-term molecular surveillance provides clues on a cattle origin for Mycobacterium bovis in Portugal

Animal tuberculosis (TB), caused by Mycobacterium bovis, is maintained in Portugal in a multi-host system, with cattle, red deer and wild boar, playing a central role. However, the ecological processes driving transmission are not understood. The main aim of this study was thus to contribute to the reconstruction of the spatiotemporal history of animal TB and to refine knowledge on M. bovis population structure in order to inform novel intervention strategies. A collection of 948 M. bovis isolates obtained during long-term surveillance (2002–2016, 15 years) of cattle (n = 384), red deer (n = 303) and wild boar (n = 261), from the main TB hotspot areas, was characterized by spoligotyping and 8 to 12-loci MIRU-VNTR. Spoligotyping identified 64 profiles and MIRU-VNTR distinguished 2 to 36 subtypes within each spoligotype, enabling differentiation of mixed or clonal populations. Common genotypic profiles within and among livestock and wildlife in the same spatiotemporal context highlighted epidemiological links across hosts and regions, as for example the SB0119-M205 genotype shared by cattle in Beja district or SB0121-M34 shared by the three hosts in Castelo Branco and Beja districts. These genomic data, together with metadata, were integrated in a Bayesian inference framework, identifying five ancestral M. bovis populations. The phylogeographic segregation of M. bovis in specific areas of Portugal where the disease persists locally is postulated. Concurrently, robust statistics indicates an association of the most probable ancient population with cattle and Beja, providing a clue on the origin of animal TB epidemics. This relationship was further confirmed through a multinomial probability model that assessed the influence of host species on spatiotemporal clustering. Two significant clusters were identified, one that persisted between 2004 and 2010, in Beja district, with Barrancos county at the centre, overlapping the central TB core area of the Iberian Peninsula, and highlighting a significant higher risk associated to cattle. The second cluster was predominant in the 2012–2016 period, holding the county Rosmaninhal at the centre, in Castelo Branco district, for which wild boar contributed the most in relative risk. These results provide novel quantitative insights beyond empirical perceptions, that may inform adaptive TB control choices in different regions.

Genotype Diversity of Mycobacterium bovis and Pathology of Bovine Tuberculosis in Selected Emerging Dairy Regions of Ethiopia

Bovine tuberculosis (bTB) is endemic in Ethiopia with higher prevalence in cattle, particularly in the central parts. Spread of Mycobacterium bovis (M. bovis) to wider regions is inevitable in uncontrolled conditions. This study was conducted to explore the pathology, characterize M. bovis strains, and describe genotypic diversity to demonstrate possible epidemiological links in emerging dairy areas of Ethiopia, namely, Mekelle and Gondar. Twenty-seven bTB positive cattle identified by the Single Intradermal Comparative Cervical Tuberculin (SICCT) test were subjected to post-mortem inspection to determine lesion distribution and pathological score. Samples from tissues with visible tuberculous or suspected non-visible lesions were processed and cultured following a standard protocol. Isolates identified as M. bovis by Region of Difference (RD)-based Polymerase Chain Reaction (PCR) were also spoligotyped to determine their spoligotype patterns. Post-mortem inspection of visceral organs indicated bTB suggestive lesions in 41% of the animals, with 25% being in the lungs. Lymph nodes from 77% of the animals had lesions. Fifty-five isolates identified from 24 of the slaughtered animals were confirmed as M. bovis. No other mycobacterial species were isolated. Spoligotyping classified strains from 21 of these animals into seven spoligotype patterns: SB0133, SB0134, SB1176, SB2233, SB2290, SB2467, and SB2520. More than one spoligotype were identified from five of these animals, and none of the last four spoligotypes had been reported in Ethiopia before. SB0134 was the most predominant type (47%) followed by SB0133 (25.5%). SB0133, SB2290, SB2467, and SB1176 are spoligotypes lacking spacers 3–7, characteristics of M. bovis strains of the African 2 (Af2) clonal complex, while SB0134, SB2233, and SB2520 do not belong to any of the established clonal complexes and likely to have a different evolutionary history. Despite a small sample size, the present study showed strain diversity with multiple genotypes identified in a single herd and even within a single animal, and the genotypes showed no sign of geographical localization, which could be a consequence of significant movement of bTB diseased cattle around the country, spreading the disease. Therefore, any future control programme of bTB in Ethiopia needs to address the risks of cattle movement.

Mycobacterium tuberculosis complex in wildlife: Review of current applications of antemortem and postmortem diagnosis

Tuberculosis (TB) is a chronic inflammatory and zoonotic disease caused by Mycobacterium tuberculosis complex (MTBC) members, which affects various domestic animals, wildlife, and humans. Some wild animals serve as reservoir hosts in the transmission and epidemiology of the disease. Therefore, the monitoring and surveillance of both wild and domestic hosts are critical for prevention and control strategies. For TB diagnosis, the single intradermal tuberculin test or the single comparative intradermal tuberculin test, and the gamma-interferon test, which is regarded as an ancillary test, are used. Postmortem examination can identify granulomatous lesions compatible with a diagnosis of TB. In contrast, smears of the lesions can be stained for acid-fast bacilli, and samples of the affected organs can be subjected to histopathological analyses. Culture is the gold standard test for isolating mycobacterial bacilli because it has high sensitivity and specificity compared with other methods. Serology for antibody detection allows the testing of many samples simply, rapidly, and inexpensively, and the protocol can be standardized in different laboratories. Molecular biological analyses are also applicable to trace the epidemiology of the disease. In conclusion, reviewing the various techniques used in MTBC diagnosis can help establish guidelines for researchers when choosing a particular diagnostic method depending on the situation at hand, be it disease outbreaks in wildlife or for epidemiological studies. This is because a good understanding of various diagnostic techniques will aid in monitoring and managing emerging pandemic threats of infectious diseases from wildlife and also preventing the potential spread of zoonotic TB to livestock and humans. This review aimed to provide up-to-date information on different techniques used for diagnosing TB at the interfaces between wildlife, livestock, and humans.

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.

Mycobacterium bovis: From Genotyping to Genome Sequencing

Mycobacterium bovis is the main pathogen of bovine, zoonotic, and wildlife tuberculosis. Despite the existence of programs for bovine tuberculosis (bTB) control in many regions, the disease remains a challenge for the veterinary and public health sectors, especially in developing countries and in high-income nations with wildlife reservoirs. Current bTB control programs are mostly based on test-and-slaughter, movement restrictions, and post-mortem inspection measures. In certain settings, contact tracing and surveillance has benefited from M. bovis genotyping techniques. More recently, whole-genome sequencing (WGS) has become the preferential technique to inform outbreak response through contact tracing and source identification for many infectious diseases. As the cost per genome decreases, the application of WGS to bTB control programs is inevitable moving forward. However, there are technical challenges in data analyses and interpretation that hinder the implementation of M. bovis WGS as a molecular epidemiology tool. Therefore, the aim of this review is to describe M. bovis genotyping techniques and discuss current standards and challenges of the use of M. bovis WGS for transmission investigation, surveillance, and global lineages distribution. We compiled a series of associated research gaps to be explored with the ultimate goal of implementing M. bovis WGS in a standardized manner in bTB control programs.

The complete genome sequence of Mycobacterium bovis Mb3601, a SB0120 spoligotype strain representative of a new clonal group

Mycobacterium bovis strain Mb3601 was isolated from the lymph node of an infected bovine in a bovine tuberculosis highly enzoonotic area of Burgundy, France. It was selected to obtain a complete genome for a new clonal complex, mainly constituted by SB0120-spoligotype strains that we propose to name "European 3". It was recently described as "clonal group I" based on whole-genome SNP analysis of 87 French strains. Here we describe the 4,365,068 bp complete genome obtained by the combination of PacBio and Illumina technologies. This genome of 65.64% G + C content includes 4024 predicted protein-coding genes, 52 tRNA, 3 rRNA and 11 copies of IS6110.

Challenges for controlling bovine tuberculosis in South Africa

All effects taken together, bovine tuberculosis (bTB) has a long-term detrimental effect on bovine herds and many wildlife species in South Africa. The disease is not only found in domestic cattle but also in African buffaloes and has to date been diagnosed in 21 wildlife species, including several rare and endangered species, thus having a potentially serious effect on conservation and biodiversity. In cattle, bTB is mostly characterised by sporadic outbreaks, but bovine herds chronically infected with the clinical disease are not uncommon. Presently, the recognised bTB control strategy in South Africa is based on 'test and slaughter', using the intradermal tuberculin test, followed by the slaughter of animals that have tested positive. Affected herds are placed under veterinary quarantine with movement restrictions until the outbreak is eradicated; this can take several years or last indefinitely if the outbreak cannot be eradicated. The same measures apply to infected buffalo populations, often with no prospect of ever being eradicated. This strategy is neither practical nor viable in the context of a communal farming system and becomes unethical when dealing with valuable wildlife reservoir hosts. Transmission of bTB between wildlife and cattle has been demonstrated and emphasises the need for an effective, affordable and culturally acceptable control strategy to curb the spread of bTB in South Africa. In countries with similar challenges, vaccination has been used and found to be promising for treating wild and domestic reservoir species and may hence be of value as a complementary tool for bTB control in South Africa.

Phylogenomic Perspective on a Unique Mycobacterium bovis Clade Dominating Bovine Tuberculosis Infections among Cattle and Buffalos in Northern Brazil

Lack of routine surveillance in countries endemic for bovine tuberculosis (TB) and limited laboratory support contributes to the inability to differentiate the Mycobacterium tuberculosis Complex species, leading to an underestimated burden of the disease. Here, Whole-Genome Sequencing of Mycobacterium bovis isolated from tissues with TB-like lesions obtained from cattle and buffalos at Marajó Island, Brazil, demonstrates that recent transmission of M. bovis is ongoing at distinct sites. Moreover, the M. bovis epidemiology in this setting is herein found to be dominated by an endemic and unique clade composed of strains evolved from a common ancestor that are now genetically differentiated from other M. bovis clades. Additionally, envisioning a rapid strain differentiation and tracing across multiple settings, 28 globally validated strain-specific SNPs were identified, three of which considered as robust markers for the M. bovis Marajó strain. In conclusion, this study contributes with data regarding the identification of a novel M. bovis phylogenetic clade responsible for ongoing transmission events in both cattle and buffalo species in Brazil, provides a framework to investigate the dissemination of this highly prevalent strain and, holds the potential to inform TB control strategies that may help to prevent the spread of bovine and zoonotic TB.

An African origin for Mycobacterium bovis

Background and objectives

Mycobacterium bovis and Mycobacterium caprae are two of the most important agents of tuberculosis in livestock and the most important causes of zoonotic tuberculosis in humans. However, little is known about the global population structure, phylogeography and evolutionary history of these pathogens.

Methodology

We compiled a global collection of 3364 whole-genome sequences from M.bovis and M.caprae originating from 35 countries and inferred their phylogenetic relationships, geographic origins and age.

Results

Our results resolved the phylogenetic relationship among the four previously defined clonal complexes of M.bovis, and another eight newly described here. Our phylogeographic analysis showed that M.bovis likely originated in East Africa. While some groups remained restricted to East and West Africa, others have subsequently dispersed to different parts of the world.

Conclusions and implications

Our results allow a better understanding of the global population structure of M.bovis and its evolutionary history. This knowledge can be used to define better molecular markers for epidemiological investigations of M.bovis in settings where whole-genome sequencing cannot easily be implemented.

Lay summary

During the last few years, analyses of large globally representative collections of whole-genome sequences (WGS) from the human-adapted Mycobacterium tuberculosis complex (MTBC) lineages have enhanced our understanding of the global population structure, phylogeography and evolutionary history of these pathogens. In contrast, little corresponding data exists for M. bovis, the most important agent of tuberculosis in livestock. Using whole-genome sequences of globally distributed M. bovis isolates, we inferred the genetic relationships among different M. bovis genotypes distributed around the world. The most likely origin of M. bovis is East Africa according to our inferences. While some M. bovis groups remained restricted to East and West Africa, others have subsequently dispersed to different parts of the world driven by cattle movements.

Genomic Polymorphism Associated with the Emergence of Virulent Isolates of Mycobacterium bovis in the Nile Delta

Mycobacterium bovis is responsible for bovine tuberculosis in both animals and humans. Despite being one of the most important global zoonotic disease, data related to the ecology and pathogenicity of bovine tuberculosis is scarce, especially in developing countries. In this report, we examined the dynamics of M. bovis transmission among dairy cattle in the Nile Delta of Egypt. Animals belonging to 27 herds from 7 governorates were tested by the Single Intradermal Comparative Skin Tuberculin (SICST), as a preliminary screen for the presence of bovine tuberculosis. Positive SICST reactors were identified in 3% of the animals spread among 40% of the examined herds. Post-mortem examination of slaughtered reactors confirmed the presence of both pulmonary and/or digestive forms of tuberculosis in > 50% of the examined animals. Targeted and whole-genome analysis of M. bovis isolates indicated the emergences of a predominant spoligotype (SB0268) between 2013–2015, suggesting a recent clonal spread of this isolate within the Nile Delta. Surprisingly, 2 isolates belonged to M. bovis BCG group, which are not allowed for animal vaccination in Egypt, while the rest of isolates belonged to the virulent M. bovis clonal complex European 2 present in Latin America and several European countries. Analysis of strain virulence in the murine model of tuberculosis indicated the emergence of a more virulent strain (MBE4) with a specific genotype. More analysis is needed to understand the molecular basis for successful spread of virulent isolates of bovine tuberculosis among animals and to establish genotype/phenotype association.

Identification of clonal complexes of Mycobacterium bovis in Brazil

Bovine tuberculosis is a disease that is widely distributed around the world. Its causative agent, Mycobacterium bovis, has characteristics of a microorganism with clonal multiplication in populations with no evidence of genetic exchange between strains, and, consequently, a group of strains can be identified as descending from a common ancestor. The aim of this study was to investigate the clonal complexes of M. bovis isolated from samples of lesions suggestive of bovine tuberculosis collected from slaughterhouses in various states of Brazil between 2006 and 2012. Ninety samples were analyzed, and it was found that 14.4% belonged to the clonal complex European1 and 81.1% to the clonal complex European2, while 4.65% were not identified as any of the four known complexes.

Molecular epidemiology and whole genome sequencing analysis of clinical Mycobacterium bovis from Ghana

Background

Bovine tuberculosis (bTB) caused by Mycobacterium bovis is a re-emerging problem in both livestock and humans. The association of some M. bovis strains with hyper-virulence, MDR-TB and disseminated disease makes it imperative to understand the biology of the pathogen.

Methods

Mycobacterium bovis (15) among 1755 M. tuberculosis complex (MTBC) isolated between 2012 and 2014 were characterized and analyzed for associated patient demography and other risk factors. Five of the M. bovis isolates were whole-genome sequenced and comparatively analyzed against a global collection of published M. bovis genomes.

Results

Mycobacterium bovis was isolated from 3/560(0.5%) females and 12/1195(1.0%) males with pulmonary TB. The average age of M. bovis infected cases was 46.8 years (7-72years). TB patients from the Northern region of Ghana (1.9%;4/212) had a higher rate of infection with M. bovis (OR = 2.7,p = 0.0968) compared to those from the Greater Accra region (0.7%;11/1543). Among TB patients with available HIV status, the odds of isolating M. bovis from HIV patients (2/119) was 3.3 higher relative to non-HIV patients (4/774). Direct contact with livestock or their unpasteurized products was significantly associated with bTB (p<0.0001, OR = 124.4,95% CI = 30.1–508.3). Two (13.3%) of the M. bovis isolates were INH resistant due to the S315T mutation in katG whereas one (6.7%) was RIF resistant with Q432P and I1491S mutations in rpoB. M. bovis from Ghana resolved as mono-phyletic branch among mostly M. bovis from Africa irrespective of the host and were closest to the root of the global M. bovis phylogeny. M. bovis-specific amino acid mutations were detected among MTBC core genes such as mce1A, mmpL1, pks6, phoT, pstB, glgP and Rv2955c. Additional mutations P6T in chaA, G187E in mgtC, T35A in Rv1979c, S387A in narK1, L400F in fas and A563T in eccA1 were restricted to the 5 clinical M. bovis from Ghana.

Conclusion

Our data indicate potential zoonotic transmission of bTB in Ghana and hence calls for intensified public education on bTB, especially among risk groups.

Wild Animal Tuberculosis: Stakeholder Value Systems and Management of Disease

When human health is put at risk from the transmission of animal diseases, the options for intervention often require input from stakeholders whose differing values systems contribute to decisions on disease management. Animal tuberculosis (TB), caused principally by Mycobacterium bovis is an archetypical zoonotic pathogen in that it can be transmitted from animals to humans and vice versa. Although elimination of zoonotic transmission of TB to humans is frequently promoted as the raison d'être for TB management in livestock, in many countries the control strategies are more likely based on minimizing the impact of sustained infection on the agricultural industry. Where wild animals are implicated in the epidemiology of the disease, the options for control and eradication can require involvement of additional stakeholder groups. Conflict can arise when different monetary and/or societal values are assigned to the affected animals. This may impose practical and ethical dilemmas for decision makers where one or more species of wild animal is seen by some stakeholders to have a greater value than the affected livestock. Here we assess the role of stakeholder values in influencing TB eradication strategies in a number of countries including Ireland, the UK, the USA, Spain, France, Australia, New Zealand and South Africa. What it reveals is that the level of stakeholder involvement increases with the complexity of the epidemiology, and that similar groups of stakeholders may agree to a set of control and eradication measures in one region only to disagree with applying the same measures in another. The level of consensus depends on the considerations of the reservoir status of the infected host, the societal values assigned to each species, the type of interventions proposed, ethical issues raised by culling of sentient wild animals, and the economic cost benefit effectiveness of dealing with the problem in one or more species over a long time frame. While there is a societal benefit from controlling TB, the means to achieve this requires identification and long-term engagement with all key stakeholders in order to reach agreement on ethical frameworks that prioritize and justify control options, particularly where culling of wild animals is concerned.

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.

Genetic profiling of Mycobacterium bovis strains from slaughtered cattle in Eritrea

Mycobacterium bovis (M.bovis) is the main causative agent for bovine tuberculosis (BTB) and can also be the cause of zoonotic tuberculosis in humans. In view of its zoonotic nature, slaughterhouse surveillance, potentially resulting in total or partial condemnation of the carcasses and organs, is conducted routinely. Spoligotyping, VNTR profiling, and whole genome sequencing (WGS) of M. bovis isolated from tissues with tuberculosis-like lesions collected from 14 cattle at Eritrea’s largest slaughterhouse in the capital Asmara, were conducted.The 14 M. bovis isolates were classified into three different spoligotype patterns (SB0120, SB0134 and SB0948) and six VNTR profiles. WGS results matched those of the conventional genotyping methods and further discriminated the six VNTR profiles into 14 strains. Furthermore, phylogenetic analysis of the M. bovis isolates suggests two independent introductions of BTB into Eritrea possibly evolving from a common ancestral strain in Europe.This molecular study revealed the most important strains of M. bovis in Eritrea and their (dis)similarities with the strains generally present in East Africa and Europe, as well as potential routes of introduction of M. bovis. Though the sample size is small, the current study provides important information as well as platform for future in-depth molecular studies on isolates from both the dairy and the traditional livestock sectors in Eritrea and the region. This study provides information onthe origin of some of the M. bovis strains in Eritrea, its genetic diversity, evolution and patterns of spread between dairy herds. Such information is essential in the development and implementation of future BTB control strategy for Eritrea.

The livestock sector plays a major role in poverty and hunger reduction in the vast majority of Africa, as a source of food, cash income, manure, draught power, transportation, savings, insurance and social status. However, for livestock to play this vital role, the impact of diseases of economic and zoonotic importance need to be reduced. Bovine tuberculosis, mainly caused by Mycobacterium bovis, is such an infectious disease. Slaughterhouse (gross pathology) surveillance, followed by bacterial culture and genotyping, are options to identify the disease-causing agents, their distribution, and enabling trace back of the sources of infections, in order to prevent their re-introduction and spread. Unfortunately, genotyping is by far not generally introduced in the continent. In the present study, tissues with tuberculosis-like lesions were collected from the Asmara municipal slaughterhouse, the largest slaughterhouse in Eritrea, and bacterial culture, classical Mycobacterium tuberculosis complex typing (Spoligotyping and VNTR profiling), as well as whole genome sequencing (WGS) were used to gain insight into the spatial and temporal distribution, genetic diversity and evolution of M. bovis strains circulating in Eritrean dairy cattle. The results revealed (dis)similarities of the Eritrean M. bovis strains with the strains generally present in Africa and Europe, potential routes of introduction to Eritrea and genetic diversity of the M. bovis strains. Future in-depth molecular studies including more samples from dairy cattle as well as cattle and goats from the traditional livestock sector are recommended.

Genetic diversity and potential routes of transmission of Mycobacterium bovis in Mozambique

Bovine tuberculosis is a zoonotic disease with largely unknown impact in Africa, with risk factors such as HIV and direct contact with animals or consumption of Mycobacterium bovis infected animal products. In order to understand and quantify this risk and design intervention strategies, good epidemiological studies are needed. Such studies can include molecular typing of M. bovis isolates. The aim of this study was to apply these tools to provide novel information concerning the distribution of bovine tuberculosis in cattle in Mozambique and thereby provide relevant information to guide policy development and strategies to contain the disease in livestock, and reduce the risk associated with transmission to humans. A collection of 178 M. bovis isolates was obtained from cattle in Mozambique. Using spoligotyping and regions of difference analysis, we classified the isolates into clonal complexes, thus reporting the first characterisation of M. bovis strains in this region. Data from MIRU-VNTR typing was used to compare isolates from a number of African countries, revealing a deeply geographically structured diversity of M. bovis. Eastern Africa appears to show high diversity, suggesting deep evolution in that region. The diversity of M. bovis in Africa does not seem to be a function of recent importation of animals, but is probably maintained within each particular region by constant reinfection from reservoir animals. Understanding the transmission routes of M. bovis in Mozambique and elsewhere is essential in order to focus public health and veterinary resources to contain bovine tuberculosis.

Bovine tuberculosis is a rather neglected zoonotic disease caused by Mycobacterium bovis that is of global concern owing to the persistence of the bacillus in reservoirs that can spread bovine tuberculosis between animals and humans. Africa remains understudied regarding this pathogen, and should be an area of concern given that in many regions the consumption of raw milk or meat from infected animals persists and the presence of HIV infection renders the population more susceptible. In order to control the disease, we need to understand M. bovis epidemiology, which includes the sources of infection. The important conclusion drawn from the work presented here is that there is a strong association between M. bovis genetic characteristics and geography. This implies that the diversity of M. bovis isolates in Mozambique does not seem to be caused by recent introductions to the territory, but is probably maintained within reservoirs in each particular region.

Whole genome sequencing of the monomorphic pathogen Mycobacterium bovis reveals local differentiation of cattle clinical isolates

Background

Bovine tuberculosis (bTB) poses serious risks to animal welfare and economy, as well as to public health as a zoonosis. Its etiological agent, Mycobacterium bovis, belongs to the Mycobacterium tuberculosis complex (MTBC), a group of genetically monomorphic organisms featured by a remarkably high overall nucleotide identity (99.9%). Indeed, this characteristic is of major concern for correct typing and determination of strain-specific traits based on sequence diversity. Due to its historical economic dependence on cattle production, Uruguay is deeply affected by the prevailing incidence of Mycobacterium bovis. With the world’s highest number of cattle per human, and its intensive cattle production, Uruguay represents a particularly suited setting to evaluate genomic variability among isolates, and the diversity traits associated to this pathogen.

Results

We compared 186 genomes from MTBC strains isolated worldwide, and found a highly structured population in M. bovis. The analysis of 23 new M. bovis genomes, belonging to strains isolated in Uruguay evidenced three groups present in the country. Despite presenting an expected highly conserved genomic structure and sequence, these strains segregate into a clustered manner within the worldwide phylogeny. Analysis of the non-pe/ppe differential areas against a reference genome defined four main sources of variability, namely: regions of difference (RD), variable genes, duplications and novel genes. RDs and variant analysis segregated the strains into clusters that are concordant with their spoligotype identities. Due to its high homoplasy rate, spoligotyping failed to reflect the true genomic diversity among worldwide representative strains, however, it remains a good indicator for closely related populations.

Conclusions

This study introduces a comprehensive population structure analysis of worldwide M. bovis isolates. The incorporation and analysis of 23 novel Uruguayan M. bovis genomes, sheds light onto the genomic diversity of this pathogen, evidencing the existence of greater genetic variability among strains than previously contemplated.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4249-6) contains supplementary material, which is available to authorized users.

Genotypes of Mycobacterium bovis strains isolated from domestic animals and wildlife in Canada in 1985-2015

Two internationally recognised and standardised genotyping methods, mycobacterial interspersed repetitive unit and variable number tandem repeat analysis (MIRU-VNTR) and spoligotyping, were applied to characterise genetic variations among 137 Mycobacterium bovis isolates recovered from Canadian domestic and wild animals during 1985-2015. Spoligotyping generated seven types that were discriminated further into12 MIRU-VNTR types. The discriminatory power indexes were estimated as 0.71 and 0.77 for spoligotyping and MIRU-VNTR typing approaches, respectively. In total, 6 prominent clusters of isolates were observed by the genotyping schemes. Four genotype clusters were exclusively observed in farmed animals. Three of these four clusters were affiliated with localised tuberculosis outbreaks, and each cluster corresponded to a single specific spoligotype (SB0140, SB0673, and SB1069) and a MIRU-VNTR profile. The fourth genotype cluster, with spoligotype SB0265 which segregated into two MIRU-VNTR types, was associated with bovine tuberculosis outbreaks in several farms across Canada during 1990-2002. Two genotype clusters of M. bovis stains were associated with wildlife reservoirs: a spoligotype SB0130 with 3 unique MIRU-VNTR profiles were observed in wood bison in Wood Buffalo National Park, and unique spoligotypes SB1070 and 1071 represented by four MIRU-VNTR profiles were recovered from cervidae species in and around the Riding Mountain National Park of Manitoba. Genotyping data confirmed M. bovis transmission between wildlife and livestock in Manitoba in 1990-2008. Overall, notwithstanding the low level of genetic diversity of Canadian M. bovis strains, the spoligotyping and MIRU-VNTR typing were useful tools in monitoring transmission of endemic strains and defining new introductions to Canada. The majority of genotypes were most likely introduced into domestic animals through live animal trade, and subsequently eliminated as a result of bovine tuberculosis outbreak investigation and eradication activities.

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.

Molecular epidemiology of Mycobacterium bovis in Cameroon

We describe the largest molecular epidemiological study of Bovine Tuberculosis (bTB) in a sub-Saharan African country with higher spatial resolution providing new insights into bTB. Four hundred and ninety-nine samples were collected for culture from 201 and 179 cattle with and without bTB-like lesions respectively out of 2,346 cattle slaughtered at Bamenda, Ngaoundere, Garoua and Maroua abattoirs between 2012–2013. Two hundred and fifty-five M. bovis were isolated, identified and genotyped using deletion analysis, Hain® Genotype MTBC, spoligotyping and MIRU-VNTR. African 1 was the dominant M. bovis clonal complex, with 97 unique genotypes including 19 novel spoligotypes representing the highest M. bovis genetic diversity observed in Africa to date. SB0944 and SB0953 dominated (63%) the observed spoligotypes. A third of animals with multiple lesions had multiple strain infections. Higher diversity but little evidence of recent transmission of M. bovis was more common in Adamawa compared to the North-West Region. The Adamawa was characterised by a high frequency of singletons possibly due to constant additions from an active livestock movement network compared to the North-West Region where a local expansion was more evident. The latter combined with population-based inferences suggest an unstable and stable bTB-endemic status in the North-West and Adamawa Regions respectively.

Tuberculosis Caused by Mycobacterium orygis in Dairy Cattle and Captured Monkeys in Bangladesh: a New Scenario of Tuberculosis in South Asia

Mycobacterium orygis, commonly known as the oryx bacillus and a newly proposed Mycobacterium tuberculosis complex subspecies, was isolated from 18 cattle in a dairy farm and two captured rhesus monkeys in a zoo in Bangladesh. All the infected animals had tuberculosis lesions in their lungs, suggesting transmission and infection with M. orygis by an airborne route. The 20 isolates were analysed using a range of conventional and molecular typing methods, and RD-deletion typing and sequencing of selected genes confirmed the isolates as M. orygis. Multiple-locus variable-number tandem repeat analysis (MLVA) allowed the isolates to be divided into three clusters based on the relatedness of their MLVA profiles. The two monkey isolates shared the same MLVA pattern with 15 of the cattle isolates, whereas the remaining three cattle isolates had different patterns, even though the latter animals had been kept in the same dairy farm. The diversity observed among isolates may suggest the bacteria have been established in this area for a long period. This study along with other recent findings that report the detection of M. orygis from animals as well as humans originating from South Asia potentially indicate endemic distribution of M. orygis in South Asia.

Molecular Relationship between Strains of M. bovis from Mexico and Those from Countries with Free Trade of Cattle with Mexico

The purpose of this study was to identify relationships between spoligotypes of M. bovis from cattle in Mexico and those reported in countries with free trade of cattle with Mexico: Australia, Canada, New Zealand and the United States of America. Mexican spoligotypes were obtained from isolates collected from cattle in different parts of the country. Spoligotypes from Canada and New Zealand were obtained from different reports in the literature. Those from the United States were obtained from the database of the National Veterinary Services Laboratory in APHIS-USDA. In order to perform the analysis in a single data set, spoligotypes were all converted to binary data and classified according to www.mbovis.org or www.pasteur-guadeloupe.fr:8081. Epidemiologic information included country and species infected. From 3,198 isolates, 174 different spoligotypes were obtained, 95 were orphans. Ninety one percent of the isolates came from the Unites States (n = 1,609) and Mexico (n = 1,323). Spoligotype SB0265 is shared between Canada and the United States in cattle and wildlife. Six spoligotypes, SB0673, SB0121, SB0145, SB0971, SB0140 and SB1165, were frequent in cattle and wildlife in the United States and cattle in Mexico, suggesting wide exchange of strains. Spoligotype SB0669 was found only in Mexico. Spoligotype SB0140 was the most common in Australia and the sixth in the United States and Mexico. In a phylogenetic analysis, spoligotype SB0140 appears as the oldest spoligotype in the data set, suggesting this as the ancestral spoligotype for all spoligotypes in the five countries. Some spoligotypes are shared by animals and humans, corroborating the zoonotic importance of M. bovis.

Molecular typing of Argentinian Mycobacterium avium subsp. paratuberculosis isolates by multiple-locus variable number-tandem repeat analysis

Multiple-locus variable number-tandem repeat analysis (MLVA) of Mycobacterium avium subspecies paratuberculosis (MAP) isolates may contribute to the knowledge of strain diversity in Argentina. Although the diversity of MAP has been previously investigated in Argentina using IS900-RFLP, a small number of isolates were employed, and a low discriminative power was reached. The aim of the present study was to test the genetic diversity among MAP isolates using an MLVA approach based on 8 repetitive loci. We studied 97 isolates from cattle, goat and sheep and could describe 7 different patterns: INMV1, INMV2, INMV11, INMV13, INMV16, INMV33 and one incomplete pattern. INMV1 and INMV2 were the most frequent patterns, grouping 76.3% of the isolates. We were also able to demonstrate the coexistence of genotypes in herds and co-infection at the organism level. This study shows that all the patterns described are common to those described in Europe, suggesting an epidemiological link between the continents.

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.