European 2--a clonal complex of Mycobacterium bovis dominant in the Iberian Peninsula

Population structure and history of Mycobacterium bovis European 3 clonal complex reveal transmission across ecological corridors of unrecognized importance in Portugal

Mycobacterium bovis causes animal tuberculosis in livestock and wildlife, with an impact on animal health and production, wildlife management, and public health. In this work, we sampled a multi-host tuberculosis community from the official hotspot risk area of Portugal over 16 years, generating the largest available data set in the country. Using phylogenetic and ecological modeling, we aimed to reconstruct the history of circulating lineages across the livestock-wildlife interface to inform intervention and the implementation of genomic surveillance within the official eradication plan. We find evidence for the co-circulation of M. bovis European 1 (Eu1), Eu2, and Eu3 clonal complexes, with Eu3 providing sufficient temporal signal for further phylogenetic investigation. The Eu3 most recent common ancestor (bovine) was dated in the 1990s, subsequently transitioning to wildlife (red deer and wild boar). Isolate clustering based on sample metadata was used to inform phylogenetic inference, unravelng frequent transmission between two clusters that represent an ecological corridor of previously unrecognized importance in Portugal. The latter was associated with transmission at the livestock-wildlife interface toward locations with higher temperature and precipitation, lower agriculture and road density, and lower host densities. This is the first analysis of M. bovis Eu3 complex in Iberia, shedding light on background ecological factors underlying long-term transmission and informing where efforts could be focused within the larger hotspot risk area of Portugal.

IMPORTANCE

Efforts to strengthen surveillance and control of animal tuberculosis (TB) are ongoing worlwide. Here, we developed an eco-phylodynamic framework based on discrete phylogenetic approaches informed by M. bovis whole-genome sequence data representing a multi-host transmission system at the livestock-wildlife interface, within a rich ecological landscape in Portugal, to understand transmission processes and translate this knowledge into disease management benefits. We find evidence for the co-circulation of several M. bovis clades, with frequent transmission of the Eu3 lineage among cattle and wildlife populations. Most transition events between different ecological settings took place toward host, climate and land use gradients, underscoring animal TB expansion and a potential corridor of unrecognized importance for M. bovis maintenance. Results stress that animal TB is an established wildlife disease without ecological barriers, showing that control measures in place are insufficient to prevent long-distance transmission and spillover across multi-host communities, demanding new interventions targeting livestock-wildlife interactions.

Genomic epidemiology sheds light on the emergence and spread of Mycobacterium bovis Eu2 Clonal Complex in Portugal

ABSTRACTAnimal tuberculosis (TB) remains a serious concern for animal and human health. Mycobacterium bovis circulates in multi-host systems, dominated by the European 2 clonal complex (Eu2) in Iberia. In this work, we use genomic epidemiology to infer the emergence, spread, and spatiotemporal patterns of Eu2 in the official epidemiological risk area of animal TB in Portugal. Phylogenetic analysis of 144 M. bovis whole-genome sequences from cattle, wild boar, and red deer, representing the 2002-2021 period, distinguished three Eu2 clades that evolved independently. The major Eu2 clade underwent phylodynamic inferences to estimate the time and location of outbreaks, host transitions, and spatial diffusion as well. The origin of this Eu2 clade was attributed to the red deer population in the Castelo Branco district, near the border with Spain. Most host transitions were intraspecific (80%), while interspecific transmissions between wildlife species (wild boar-red deer), and between wild boar and cattle, were highly supported. Phylogeographic reconstruction evidenced that most transitions (82%) occur within municipalities, highlighting local transmission corridors.Our study indicates that M. bovis continues to spread at the cattle-wildlife interface within the animal TB hotspot area, possibly driven by the foraging behaviour of wild boar near agricultural lands. Red deer seems to be an important driver of TB within wildlife hosts, while the wild boar links the multi-host wildlife community and livestock. This work highlights the value of combining genomic epidemiology with phylodynamic inference to resolve host jumps and spatial patterns of M. bovis, providing real-time clues about points of intervention.

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.

Unraveling the epidemiology of Mycobacterium bovis using whole-genome sequencing combined with environmental and demographic data

When studying the dynamics of a pathogen in a host population, one crucial question is whether it transitioned from an epidemic (i.e., the pathogen population and the number of infected hosts are increasing) to an endemic stable state (i.e., the pathogen population reached an equilibrium). For slow-growing and slow-evolving clonal pathogens such as Mycobacterium bovis, the causative agent of bovine (or animal) and zoonotic tuberculosis, it can be challenging to discriminate between these two states. This is a result of the combination of suboptimal detection tests so that the actual extent of the pathogen prevalence is often unknown, as well as of the low genetic diversity, which can hide the temporal signal provided by the accumulation of mutations in the bacterial DNA. In recent years, the increased availability, efficiency, and reliability of genomic reading techniques, such as whole-genome sequencing (WGS), have significantly increased the amount of information we can use to study infectious diseases, and therefore, it has improved the precision of epidemiological inferences for pathogens such as M. bovis. In this study, we use WGS to gain insights into the epidemiology of M. bovis in Cameroon, a developing country where the pathogen has been reported for decades. A total of 91 high-quality sequences were obtained from tissue samples collected in four abattoirs, 64 of which were with complete metadata. We combined these with environmental, demographic, ecological, and cattle movement data to generate inferences using phylodynamic models. Our findings suggest M. bovis in Cameroon is slowly expanding its epidemiological range over time; therefore, endemic stability is unlikely. This suggests that animal movement plays an important role in transmission. The simultaneous prevalence of M. bovis in co-located cattle and humans highlights the risk of such transmission being zoonotic. Therefore, using genomic tools as part of surveillance would vastly improve our understanding of disease ecology and control strategies.

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.

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.

Mycobacterium bovis and M. caprae in Bulgaria: insight into transmission and phylogeography gained through whole-genome sequencing

Background

This study aimed to characterize recent Mycobacterium bovis/M. caprae isolates from Bulgaria by whole-genome sequencing (WGS) to gain a first insight into their molecular diversity, transmission, and position within the global phylogeography of this important zoonotic species.

Results

The isolates were obtained from cattle in diverse locations of Bulgaria in 2015-2020 and were identified by microbiological and PCR assays. WGS data were used for phylogenetic analysis that also included M. bovis global dataset. Thirty-seven M. bovis/caprae isolates from Bulgaria were studied and 34 of them were SNP genotyped. The isolates were subdivided into 3 major phylogenetic groups. Type Mbovis-13 (Eu2 complex [western Europe and northern Africa]) included one isolate. Mbovis-37 type included 5 isolates outside of known clonal complexes. The Bulgarian M. caprae isolates formed a sub-group within the Mcaprae-27B cluster which also included 22 M. caprae isolates from Poland, Spain, Germany, and the Republic of Congo. The Bulgarian M. caprae isolates share their latest common ancestors with Spanish isolates. The Mbovis-37 group shares a distant common ancestor (pairwise distance 22-29 SNPs) with an isolate from Poland but was very distant (> 200 SNPs) from the rest of the tree. The Mbovis-13 group shares a common ancestor with two human isolates from Germany. Phylogeographically, both M. bovis clades had limited circulation in northeastern Bulgaria while the majority of the studied isolates (M. caprae) were from central and western provinces. A phylogenetic network-based analysis demonstrated that 11 Bulgarian isolates were separated by 1 to 6 SNPs within four clusters, mostly forming pairs of isolates.

Conclusion

The obtained WGS analysis positioned the Bulgarian isolates within the global phylogeography of M. bovis/M. caprae. Hypothetically, the observed phylogenetic diversity may not have resulted from livestock trade routes, but instead may reflect the deeply rooted M. bovis/M. caprae phylogeography of Europe. A high level of genetic divergence between the majority of the studied isolates suggests limited active transmission of bTB in Bulgaria during the survey period. At the same time, a possibility of the endemic presence of circulating bTB strains in the form of the latent persistent disease cannot be ruled out.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12917-022-03249-w.

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.

Whole Genome Sequencing Refines Knowledge on the Population Structure of Mycobacterium bovis from a Multi-Host Tuberculosis System

Classical molecular analyses of Mycobacterium bovis based on spoligotyping and Variable Number Tandem Repeat (MIRU-VNTR) brought the first insights into the epidemiology of animal tuberculosis (TB) in Portugal, showing high genotypic diversity of circulating strains that mostly cluster within the European 2 clonal complex. Previous surveillance provided valuable information on the prevalence and spatial occurrence of TB and highlighted prevalent genotypes in areas where livestock and wild ungulates are sympatric. However, links at the wildlife-livestock interfaces were established mainly via classical genotype associations. Here, we apply whole genome sequencing (WGS) to cattle, red deer and wild boar isolates to reconstruct the M. bovis population structure in a multi-host, multi-region disease system and to explore links at a fine genomic scale between M. bovis from wildlife hosts and cattle. Whole genome sequences of 44 representative M. bovis isolates, obtained between 2003 and 2015 from three TB hotspots, were compared through single nucleotide polymorphism (SNP) variant calling analyses. Consistent with previous results combining classical genotyping with Bayesian population admixture modelling, SNP-based phylogenies support the branching of this M. bovis population into five genetic clades, three with apparent geographic specificities, as well as the establishment of an SNP catalogue specific to each clade, which may be explored in the future as phylogenetic markers. The core genome alignment of SNPs was integrated within a spatiotemporal metadata framework to further structure this M. bovis population by host species and TB hotspots, providing a baseline for network analyses in different epidemiological and disease control contexts. WGS of M. bovis isolates from Portugal is reported for the first time in this pilot study, refining the spatiotemporal context of TB at the wildlife-livestock interface and providing further support to the key role of red deer and wild boar on disease maintenance. The SNP diversity observed within this dataset supports the natural circulation of M. bovis for a long time period, as well as multiple introduction events of the pathogen in this Iberian multi-host system.

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.

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.

Phylogenetic relationships investigation of Mycobacterium caprae strains from sympatric wild boar and goats based on whole genome sequencing

Tuberculosis (TB) in wildlife challenges epidemiological surveillance and disease control. An outbreak of TB was detected in a free‐ranging wild boar population of a Natural Park in Catalonia (Spain) and the outbreak investigation was conducted in the area. During the study period (2015–2020), 278 wild boars were analysed by gross pathology, histopathology, mycobacterial culture and DVR‐spoligotyping. In addition, all cattle (49) and goat (47) herds of the area were tested with tuberculin skin test. TB compatible lesions were detected in 21 wild boars, and Mycobacterium caprae was isolated in 17 of them with two different spoligotypes: SB0415 (13) and SB1908 (4). Only two goat herds showed TB positive animals that were subsequently slaughtered. M. caprae with the spoligotypes SB0416 and SB0415 were isolated from these animals. To investigate the phylogenetic relationships and the transmission chain of the outbreak, nine strains isolated from six wild boars and three goats of the study area were analysed by whole genome sequencing (WGS) followed by single nucleotide polymorphism (SNP) analysis by maximum likelihood and median‐joining network inference methods. Results indicated that infected wild boars maintained M. caprae strains circulation in their own population and have likely transmitted the infection to goats, thus acting as TB reservoirs, compromising the success of livestock TB eradication campaigns and posing a risk for public health. The results also highlighted the usefulness of WGS followed by SNP analysis in providing relevant epidemiological information when detailed contact data are missing.

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.

Accurate Phylogenetic Relationships Among Mycobacterium bovis Strains Circulating in France Based on Whole Genome Sequencing and Single Nucleotide Polymorphism Analysis

In recent years the diversity of the French Mycobacterium bovis population responsible for bovine tuberculosis (bTB) outbreaks since 1970 has been described in detail. To further understand bTB evolution in France, we used single nucleotide polymorphisms (SNPs) based on whole genome sequence versus classical genotyping methods in order to identify accurate phylogenetic relationships between M. bovis strains. Whole genome sequencing was carried out on a selection of 87 strains which reflect the French M. bovis population’s genetic diversity. Sequences were compared to the M. bovis reference genome AF2122/97. Comparison among the 87 genomes revealed 9,170 sites where at least one strain shows a SNP with respect to the reference genome; 1,172 are intergenic and 7,998 in coding sequences, of which 2,880 are synonymous and 5,118 non-synonymous. SNP-based phylogenetic analysis using these 9,170 SNP is congruent with the cluster defined by spoligotyping and multilocus variable number of tandem repeat analysis typing. In addition, some SNPs were identified as specific to genotypic groups. These findings suggest new SNP targets that can be used for the development of high-resolving methods for genotyping as well as for studying M. bovis evolution and transmission patterns. The detection of non-synonymous SNPs on virulence genes enabled us to distinguish different clusters. Our results seem to indicate that genetically differentiated clusters could also display distinctive phenotypic traits.

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.

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.

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.

Exposure of Threatened Accipitridae to Mycobacterium bovis Calls for Active Surveillance

Anthropogenic activities have cumulatively led to the dramatic decline of world populations of vultures that currently face serious survival challenges in several regions of the world. In Portugal, the three resident species qualify as endangered and are under conservation efforts, mainly in the central east and south-east regions, where habitat protection and artificial feeding stations were implemented. Concurrently, the areas under protection are highly affected by tuberculosis (TB) in cattle and wild ungulates, whose potentially infected carcasses may naturally or artificially be used as feed by local vultures. In this work, we opportunistically surveyed populations of Eurasian griffon (Gyps fulvus) and Eurasian black vulture (Aegypius monachus) for the presence of Mycobacterium bovis. Nine pathogenic mycobacteria, including one M. bovis isolate, were cultured from the oropharynx of nine of the surveyed vultures (n = 55), sampled in recovery centres or in artificial feeding stations. Genotyping of the M. bovis strain indicated spoligotype SB0121, the most frequent type in Portugal, and a unique MIRU-VNTR profile that differed in two loci from the profiles of SB0121 bovine and deer strains from the same geographical area. The M. bovis-positive griffon exhibited poor clinical condition when admitted to the recovery centre; however, clinical evidence of TB was not present. Although the significance of M. bovis isolation in this vulture specimen could not be ascertained and despite the accepted notion that vultures are naturally resistant to microbial pathogens, the sanitary follow-up of Accipitridae vulture populations in TB-hotspot areas is essential to safeguard ongoing conservation efforts and also to evaluate the suitability of standing legislation on deliberate supplementary feeding schemes for menaced birds of prey.

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.

Overview and phylogeny of Mycobacterium tuberculosis complex organisms: implications for diagnostics and legislation of bovine tuberculosis

Members of the Mycobacterium tuberculosis complex (MTBC) cause a serious disease with similar pathology, tuberculosis; in this review, bovine tuberculosis will be considered as disease caused by any member of the MTBC in bovids. Bovine tuberculosis is responsible for significant economic loss due to costly eradication programs and trade limitations and poses a threat to both endangered and protected species as well as to public health. We here give an overview on all members of the MTBC, focusing on their isolation from different animal hosts. We also review the recent advances made in elucidating the evolutionary and phylogenetic relationships of members of the MTBC. Because the nomenclature of the MTBC is controversial, its members have been considered species, subspecies or ecotypes, this review discusses the possible implications for diagnostics and the legal consequences of naming of new species.

Splitting of a prevalent Mycobacterium bovis spoligotype by variable-number tandem-repeat typing reveals high heterogeneity in an evolving clonal group

Mycobacterium bovis populations in countries with persistent bovine tuberculosis usually show a prevalent spoligotype with a wide geographical distribution. This study applied mycobacterial interspersed repetitive-unit-variable-number tandem-repeat (MIRU-VNTR) typing to a random panel of 115 M. bovis isolates that are representative of the most frequent spoligotype in the Iberian Peninsula, SB0121. VNTR typing targeted nine loci: ETR-A (alias VNTR2165), ETR-B (VNTR2461), ETR-D (MIRU4, VNTR580), ETR-E (MIRU31, VNTR3192), MIRU26 (VNTR2996), QUB11a (VNTR2163a), QUB11b (VNTR2163b), QUB26 (VNTR4052), and QUB3232 (VNTR3232). We found a high degree of diversity among the studied isolates (discriminatory index [D] = 0.9856), which were split into 65 different MIRU-VNTR types. An alternative short-format MIRU-VNTR typing targeting only the four loci with the highest variability values was found to offer an equivalent discriminatory index. Minimum spanning trees using the MIRU-VNTR data showed the hypothetical evolution of an apparent clonal group. MIRU-VNTR analysis was also applied to the isolates of 176 animals from 15 farms infected by M. bovis SB0121; in 10 farms, the analysis revealed the coexistence of two to five different MIRU types differing in one to six loci, which highlights the frequency of undetected heterogeneity.

Tuberculosis epidemiology in islands: insularity, hosts and trade

Because of their relative simplicity and the barriers to gene flow, islands are ideal systems to study the distribution of biodiversity. However, the knowledge that can be extracted from this peculiar ecosystem regarding epidemiology of economically relevant diseases has not been widely addressed. We used information available in the scientific literature for 10 old world islands or archipelagos and original data on Sicily to gain new insights into the epidemiology of the Mycobacterium tuberculosis complex (MTC). We explored three nonexclusive working hypotheses on the processes modulating bovine tuberculosis (bTB) herd prevalence in cattle and MTC strain diversity: insularity, hosts and trade. Results suggest that bTB herd prevalence was positively correlated with island size, the presence of wild hosts, and the number of imported cattle, but neither with isolation nor with cattle density. MTC strain diversity was positively related with cattle bTB prevalence, presence of wild hosts and the number of imported cattle, but not with island size, isolation, and cattle density. The three most common spoligotype patterns coincided between Sicily and mainland Italy. However in Sicily, these common patterns showed a clearer dominance than on the Italian mainland, and seven of 19 patterns (37%) found in Sicily had not been reported from continental Italy. Strain patterns were not spatially clustered in Sicily. We were able to infer several aspects of MTC epidemiology and control in islands and thus in fragmented host and pathogen populations. Our results point out the relevance of the intensity of the cattle commercial networks in the epidemiology of MTC, and suggest that eradication will prove more difficult with increasing size of the island and its environmental complexity, mainly in terms of the diversity of suitable domestic and wild MTC hosts.

Molecular discrimination of Mycobacterium bovis in São Paulo, Brazil

Mycobacterium bovis, a member of the Mycobacterium tuberculosis complex, is the most common agent of cattle tuberculosis, a zoonosis that causes losses in meat and milk production in several countries. In order to support epidemiological studies aimed at controlling the disease, several methods for molecular discrimination of M. bovis isolates have recently been developed. The most frequently used are spacer oligonucleotide typing (spoligotyping), mycobacterial interspersed repetitive units (MIRU), and exact tandem repeat (ETR), but they all have different discriminatory power. In the present study, allelic diversity was calculated for each MIRU and ETR locus, and the Hunter-Gaston discriminatory index (HGI) was calculated for spoligotyping, 10 MIRUs, and 3 ETRs, in 116 isolates of M. bovis obtained from cattle. The analysis of allelic diversity indicated that MIRUs 16, 26, and 27, and ETRs A, B, and C, showed the greatest diversity between the assayed loci. The HGIs for each of the techniques were: spoligotyping=0.738381; MIRU=0.829835; and ETR=0.825337. The associations of the methods' improved discriminatory power were: spoligotyping+MIRU=0.930585; spoligotyping+ETR=0.931034; and MIRU+ETR=0.953373. The greatest discriminatory power was obtained when the three techniques were associated (HGI=0.98051). Considering the analyses of the present study, spoligotyping should be the first method to be used because it differentiates M. bovis from the other members of the Mycobacterium tuberculosis complex. As the associations of MIRU and ETR with spoligotyping resulted in nearly identical HGIs, ETR seems to be the best choice after spoligotyping, because it is faster and more economical than MIRU. Finally, MIRU should be the last method used. In spite of this finding, the choice of the method used should be based on the discriminatory power necessary for the objective at hand.

The global distribution and phylogeography of Mycobacterium bovis clonal complexes

The consequences of the clonality of the Mycobacterium tuberculosis complex are described and, in particular, the identification of clonal complexes. Clonal complexes are groups of strains all descended from a single cell that was the most recent common ancestor (MRCA) of the clonal complex and all bearing characteristics derived from the MRCA. Three clonal complexes of Mycobacterium bovis have been identified and called African 1, African 2 and European 1. Members of each clonal complex have a distinct spoligotype signature and are identified by a unique deletion present in each member of the clonal complex. The African 1 and African 2 clonal complexes are geographically localised to Central-West Africa and East Africa, respectively and have not been found in cattle outside of these regions. However, the European 1 clonal complex is globally distributed and has been identified in the British Isles, former British colonies, The Americas as well as Kazakhstan and Korea. It is suggested that modern cattle breeds, such as Herefords, bred in the UK in the 18th Century, would provide a good vehicle for the global distribution of this closely related group of strains. The phylogeography of M. bovis and M. tuberculosis are compared and in particular the diversity of M. tuberculosis in Africa, compared with the localised dominance of M. bovis clonal complexes, is highlighted. Finally, the practical use of M. bovis clonal complexes is reviewed, however, a more significant use of M. bovis clonal complexes is to generate testable hypotheses to understand the pathogenicity and spread of this important veterinary pathogen. This first look at the phylogeography of M. bovis clonal complexes has shown surprising geographical localisation of molecular types but also hints strongly that much of the worldwide distribution of bovine pathogen resulted from human trade in cattle within the last 200 years.