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

David Couvin, Iñaki Cervera-Marzal, Audrey David, Yann Reynaud, Nalin Rastogi

► To cite this version:

David Couvin, Iñaki Cervera-Marzal, Audrey David, Yann Reynaud, Nalin Rastogi. SITVITBovis-a publicly available database and mapping tool to get an improved overview of animal and human cases caused by Mycobacterium bovis. Database - The journal of Biological Databases and Curation, 2022, 2022 (2022), pp.baab081. 10.1093/database/baab081. hal-05166350

HAL Id: hal-05166350

https://hal.science/hal-05166350v1

Submitted on 17 Jul 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# SITVITBovis—a publicly available database and mapping tool to get an improved overview of animal and human cases caused by *Mycobacterium bovis*

David Couvin ID*, Iñaki Cervera-Marzal, Audrey David, Yann Reynaud and Nalin Rastogi*

WHO Supranational TB Reference Laboratory-TB and Mycobacteria Unit, Institut Pasteur de Guadeloupe, Abymes 97183, Guadeloupe, France

*Correspondence may also be addressed to David Couvin. Email: dcouvin@pasteur-guadeloupe.fr and Nalin Rastogi. Tel: +590590897661; Email: nrastogi@pasteur-guadeloupe.fr

Citation details: Couvin, D., Cervera-Marzal, I., David, A. et al. SITVITBovis—a publicly available database and mapping tool to get an improved overview of animal and human cases caused by *Mycobacterium bovis*. *Database* (2022) Vol. 2022: article ID baab081; DOI: https://doi.org/10.1093/database/baab081

## Abstract

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 URL: http://www.pasteur-guadeloupe.fr:8081/SITVIT_Bovis

## Introduction

Bovine tuberculosis (bTB) is an infectious disease caused by *Mycobacterium bovis*, a member of the *Mycobacterium tuberculosis* complex (MTBC). *M. bovis* alongside with *M. caprae* infect a wide range of mammalian hosts, including humans. An estimated 140 000 new cases of zoonotic tuberculosis (TB) and 11 400 deaths occurred globally in 2019 due to *M. bovis* (1). Knowing the difficulty to easily distinguish *M. bovis* from *M. tuberculosis*, these numbers possibly represent an underestimation of *M. bovis* cases worldwide.

According to the World Organisation for Animal Health (http://www.oie.int/en/), isolates of *M. bovis* or *M. caprae* were obtained from various animal species such as buffalos, bison, sheep, goats, horses, camels, pigs, wild boars, deer, antelopes, dogs, cats, foxes, minks, badgers, ferrets, rats, primates, llamas, kudu, eland, tapirs, elk, elephants, sitatungas, oryx, addax, rhinoceros, opossums, squirrels, otters, seals, hares, moles, raccoons, coyotes and several predatory felines including lions, tigers, leopards and lynx. Domestication of animals has long been a usual practice for mankind, and this practice seems to have caused irreversible

effects in the evolution of bTB (2). Note that some animals could be infected by their own MTBC species (e.g. *Mycobacterium microti*, *Mycobacterium pinnipedii*, *Mycobacterium mungi*, *Dassie bacillus*, *Mycobacterium oryx* and *Mycobacterium leprae*). Details are provided in Supplementary File 1/Supplementary Table S1. Genotyping methods such as spoligotyping (3) and mycobacterial interspersed repetitive units-variable number of tandem repeats (MIRU-VNTRs) typing (4) have proven to be efficacious to distinguish the strains belonging to MTBC (5-7). Despite the limitations (such as homoplasy) of these polymerase chain reaction-based methods (8), these techniques are well used worldwide for MTBC family identification. Application of fingerprinting tools facilitates analysis of the molecular epidemiology of *M. bovis* in animal-to-human, human-to-human and even animal-to-animal transmission (9). In addition, studies based on whole-genome sequencing (WGS) analyses provide additional details which tend to become the norm in genomic epidemiology and evolutionary studies (10-12).

Development of worldwide or local databases dedicated to animal and human TB improves our global
understanding on evolution and propagation of the disease. Mbovis.org database (www.Mbovis.org) (13) provides authoritative names for spoligotypes (SB numbers) of all members of MTBC isolates of animal origin. MycoDB.es (https://www.visavet.es/mycodb/) is another database focusing on zoonotic TB hosted in Spain (14).

The purpose of this study is to describe the SITVITBovis database containing available genotyping information on 25741 M. bovis or M. caprae isolates. In addition to this information, preliminary data extracted from 188 raw sequence reads allowed to make some links with WGS and classical genotyping data. By publicly releasing this multimarker SITVITBovis database which incorporates user-friendly online tools and interfaces, we hope to serve the global research community with a concerted and coordinated response to monitor and assess global bTB spread.

# Materials and methods

## Ethics statement and collection of data

All the data (human and other mammalian isolates) were obtained from collaborating laboratories as described in SITVIT2 database (15). Data were duly de-identified prior to database entry. SITVITBovis is an excerpt from SITVIT2 database focusing on bTB isolates. It contains additional information on host (e.g. cattle, buffalo, deer, human, etc.) and SB_numbers collected from Mbovis.org database (13). Unlike Mbovis.org database that only provides spoligotypes, SITVITBovis provides available epidemiological data, such as information on host, WGS and MIRU-VNTRs. Many isolates were obtained from MycoDB.es study, accounting for more than 17 000 isolates at the time of the present study (16). Other isolates were obtained from other published studies (14, 17–20). SITVITBovis database aims to grow in the future and further updates will be applied, notably with the release of the upcoming SITVITEXTEND database. A user guide is provided online to facilitate navigation through the internet.

## Genotyping markers and WGS

The genotyping data included in SITVITBovis were similar to those previously described in SITVIT2 database (15). The methods used were spoligotyping (3) and MIRU-VNTRs typing, comprising 5-locus exact tandem repeats (ETR-A to E) (21) and 12- and 15-loci MIRU formats (4). The order of MIRU loci is as follows: 12-loci MIRU patterns—MIRU 2, 4, 10, 16, 20, 23, 24, 26, 27, 31, 39 and 40; 15-loci MIRU patterns—MIRU 4, 10, 16, 26, 31 and 40; ETR-A; ETR-C; QUB-11b; QUB-26; QUB-4156; Mtub04; Mtub21; Mtub30 and Mtub39. Spoligotype international type (SIT) indicates spoligotyping patterns found at least two times in database; VNTR international type (VIT) indicates 5-locus ETR patterns found at least two times in database; and 12- or 15-MIRU international type (12-MIT or 15-MIT) indicates 12- or 15-MIRU-VNTR patterns found at least two times in our database. Mbovis.org website was used to extract the SB numbers matching with some SIT numbers from SITVITBovis database (Supplementary Table S2). bTB sublineages (BOV, BOV_1, BOV_2, BOV_3 and BOV_4 - CAPRAE) have been previously described in SITVIT classification using

revised spoligotyping rules (7, 15). Note that BOV sublineage corresponds to previously labeled *BOV_Like* sublineage in SITVITWEB. To make a link with WGS data, we searched for genomes identified as 'Mycobacterium bovis not BCG' deposited in Sequence Read Archive (SRA) of the National Center for Biotechnology Information (NCBI) as suggested by a recent study (11). Then, an in-house pipeline algorithm (written in Perl language) was used to link WGS data extracted from the European Nucleotide Archive with classical genotyping information and sequencing-based tools. A total of 188 sequence reads were available with easily findable country information. The following tools were used for the analyses: SpoTyping, MIRUReader, SPAdes, SpolLineages, TBProfiler and Fast-lineage-caller (22–27). A preliminary strain list with read accessions is available on the 'query' web page. Formal presentation and query of WGS isolates will be available in an upcoming dedicated database which will contain more isolates.

## Computing approach for database construction

SITVITBovis is a publicly accessible website enabling the viewing of global spread of M. bovis isolates. It has been designed to work optimally with Google Chrome and Firefox. The web interface has been implemented using Java technology (Java Server Pages, Asynchronous JavaScript and Ajax), Google Code API and XML, under the integrated development environment Eclipse (https://www.eclipse.org/). Data were integrated within a MySQL database. Public access to the database is strictly on a read-only basis; therefore, no direct update of the database is allowed from this website. The web application is hosted and deployed on an Apache Tomcat Server (version 6).

## Phylogenetic, statistical and bioinformatics analyses

Existing bioinformatics tools have been used to realize the phylogenetic analysis. Minimum spanning trees (MSTs) based on spoligotypes or MIRU-VNTRs were drawn using BioNumerics software version 6.6 (Applied Maths, Sint-Martens-Latem, Belgium) or MLVA Compare version 1.03 software (GenoScreen; Lille, France). MSTs are undirected connected graphs which link all nodes (representing the isolates) together with the fewest possible linkages between nearest neighbors. Spoligoforests based on spoligotypes were drawn using the SpolTools software (28, 29); they allowed to describe and visualize the potential parent-to-descendant relationships among spoligotypes. In some cases, Spoligoforests have been colored and reshaped using GraphViz software (http://www.graphviz.org/). Contrary to the MST, the spoligoforest trees are directed graphs which only evolve by loss of spacers. In these trees, nodes are not necessarily all connected. In case of too many changes between two strains, there are no edges linking them. STATA software version 12 was used for descriptive and univariate analyses. Pearson's chi-square test and Fisher's exact test were used for comparison of different parameters and P values of <0.05 were considered as statistically significant. V-DICE tool was used to compare discriminatory powers of typing methods, providing Hunter-Gaston Diversity Index (HGDI) or Simpson's diversity index (30).
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Worldwide distribution of bTB isolates and their human or animal hosts by United Nations (UN) sub-region.</figcaption></figure>

# Results

## Structure and main functionalities of SITVITBovis

SITVITBovis was built based on the same architecture as SITVIT2 database. The users can browse through the web pages according to their needs. The following pages are accessible via the tabs ‘Database Description’, ‘Search’, ‘Analysis’, ‘Online Tools’, ‘Statistical Analysis’, ‘WGS data’ and ‘Links and Others’: (i) The ‘Database Description’ page provides general information about genotyping molecular markers contained in the database as well as global distribution of *M. bovis* strains and mammalian hosts represented in the database. (ii) The ‘Search’ page allows users to query SITVITBovis database according to several criteria (such as spoligotype, 12-loci MIRU, SIT, SB number, 12-MIT, lineage, isolation country, drug resistance, host, etc.). Regular expressions can be used. (iii) The ‘Analysis’ page allows users to analyze their own data based on the example file provided on the web page. (iv) The ‘Online Tools’ page allows users to get information on the distribution of isolates in function of several characteristics such as genotyping markers, lineages, hosts, country or city of isolation, etc. (v) The ‘Statistical Analysis’ provides information on the evolution of strains over time and space in function of hosts. (vi) The ‘WGS data’ page contains a preliminary strain list linking read accessions with information on classical genotyping data and drug resistance. (vii) The ‘Links and Others’

page provides supplemental information on interesting links and lab.

## Worldwide diversity of bTB genotypes

The SITVITBovis database provides a broad overview of the geographic distribution of various phylogenetical sublineages belonging to *M. bovis* (*BOV_1*, *BOV_2*, *BOV_3* and *BOV_4* – *CAPRAE*). This web application also helps distinguish between strains isolated from humans and other mammals in order to better understand transmission pathways and other factors that could potentially be responsible for the global spread of bTB. In our study, among the 25 741 isolates, *BOV_1* sublineage was globally the most predominant, representing 67.7% (n = 17 427) of bTB clinical isolates, followed by *BOV_2* (n = 3258 or 12.7%), *BOV* (n = 3157 or 12.3%), *BOV_4* – *CAPRAE* (n = 1735 or 6.7%) and *BOV_3* (n = 164 or 0.6% of isolates). (i) *BOV_1* was notably predominantly found (number of isolates ≥ 25) in Russia, Western and Southern Europe, the whole African continent (North, Middle, West, Austral and East), the Middle-East (Western Asia), Southern Asia and North and South America, with proportions between 43% and 99% (Figure 1). (ii) *BOV_2* sublineage was predominantly found in Northern Europe (84% of clinical isolates), Eastern Asia (68%), Australasia (54%), South America (35%), North America (25%), Central
America (9%) and Austral Africa (8%) (Figure 1). (iii) BOV<sub>3</sub> sublineage was scarcely found in our study, being limited to East Asia (24%), North America (12%) and South America (4%). (iv) BOV_4 – CAPRAE lineage was predominantly found in Eastern Europe (representing 83% of isolates). This lineage represented 8% of clinical isolates in Western Europe and South America and 6% of the isolates in North Africa. (v) BOV sublineage was predominantly found in Central America (59% of isolates), followed by East Africa (46% of strains), and with proportions from 12% to 23% in North, West and Central Africa; North and South America; Western and Southern Europe; and Western and Southern Asia (Figure 1).

The database provided information on strains belonging to human hosts (n = 473 isolates) and other mammalian hosts (mainly cattle; n = 18769 isolates). Most of the bTB isolates from animals were found in Southern Europe (predominantly in Spain), whereas available isolates from humans were mainly found in Western Europe (n = 149), Northern Europe (n = 120) and followed by Northern Asia or Russia (n = 108) (Figure 1). Among human isolates, gender of patients was known for 301 isolates, and the global male/female sex ratio was 173/128 = 1.35. In countries where at least 10 human isolates were recorded, the highest sex ratio was observed in Southern Europe (male/female sex ratio=3.09), as opposed to the lowest ratio in Western

<table><caption>**Table 2.** Diversity of spoligotypes and MIRU-VNTR patterns and loci</caption><thead><tr><th>Spoligotyping and MIRU-VNTRs scheme or loci<sup>*</sup></th><th>Hunter-Gaston Diversity (95% confidence interval)</th><th>Simpson's diversity (95% confidence interval)</th></tr></thead><tbody><tr><td>12-loci MIRU-VNTRs</td><td>0.973 (0.967-0.978)</td><td>0.970 (0.964-0.975)</td></tr><tr><td>MIRU16 (VNTR)<sup>*</sup></td><td>0.623 (0.589-0.657)</td><td>0.621 (0.587-0.655)</td></tr><tr><td>MIRU04 (VNTR)<sup>*</sup></td><td>0.618 (0.575-0.662)</td><td>0.616 (0.573-0.660)</td></tr><tr><td>MIRU31 (VNTR)<sup>*</sup></td><td>0.537 (0.487-0.587)</td><td>0.535 (0.485-0.585)</td></tr><tr><td>MIRU10 (VNTR)<sup>*</sup></td><td>0.534 (0.481-0.587)</td><td>0.532 (0.479-0.585)</td></tr><tr><td>MIRU26 (VNTR)<sup>*</sup></td><td>0.437 (0.376-0.498)</td><td>0.436 (0.375-0.496)</td></tr><tr><td>MIRU23 (VNTR)</td><td>0.265 (0.200-0.329)</td><td>0.264 (0.199-0.328)</td></tr><tr><td>MIRU27 (VNTRc)</td><td>0.259 (0.195-0.323)</td><td>0.258 (0.194-0.322)</td></tr><tr><td>MIRU24 (VNTR)</td><td>0.099 (0.052-0.146)</td><td>0.099 (0.052-0.146)</td></tr><tr><td>MIRU40 (VNTR)<sup>*</sup></td><td>0.072 (0.034-0.109)</td><td>0.071 (0.034-0.109)</td></tr><tr><td>MIRU20 (VNTR)</td><td>0.073 (0.032-0.114)</td><td>0.073 (0.032-0.114)</td></tr><tr><td>MIRU02 (VNTR)</td><td>0.034 (0.005-0.063)</td><td>0.034 (0.005-0.063)</td></tr><tr><td>MIRU39 (VNTR)</td><td>0.000 (0.000-0.025)</td><td>0.000 (0.000-0.025)</td></tr><tr><td>15-loci MIRU-VNTRs</td><td>0.979 (0.962-0.996)</td><td>0.964 (0.947-0.981)</td></tr><tr><td>QUB-11b (VNTR)</td><td>0.668 (0.613-0.723)</td><td>0.658 (0.603-0.713)</td></tr><tr><td>ETR-A (VNTR)<sup>*</sup></td><td>0.653 (0.611-0.695)</td><td>0.651 (0.610-0.693)</td></tr><tr><td>QUB-4156 (VNTR)</td><td>0.544 (0.450-0.639)</td><td>0.536 (0.441-0.630)</td></tr><tr><td>ETR-C (VNTR)<sup>*</sup></td><td>0.481 (0.430-0.532)</td><td>0.480 (0.429-0.531)</td></tr><tr><td>QUB-26 (VNTR)</td><td>0.476 (0.374-0.578)</td><td>0.469 (0.367-0.571)</td></tr><tr><td>Mtub21 (VNTR)</td><td>0.455 (0.339-0.571)</td><td>0.448 (0.332-0.564)</td></tr><tr><td>Mtub04 (VNTR)</td><td>0.441 (0.336-0.546)</td><td>0.434 (0.329-0.539)</td></tr><tr><td>Mtub30 (VNTR)</td><td>0.433 (0.334-0.532)</td><td>0.426 (0.327-0.525)</td></tr><tr><td>Mtub39 (VNTR)</td><td>0.092 (0.000-0.188)</td><td>0.090 (0.000-0.186)</td></tr><tr><td>5-locus ETRs (VNTR)</td><td>0.930 (0.914-0.947)</td><td>0.928 (0.911-0.944)</td></tr><tr><td>ETR-B (VNTR)</td><td>0.559 (0.503-0.615)</td><td>0.558 (0.502-0.614)</td></tr><tr><td>ETR-D</td><td>0.237 (0.180-0.294)</td><td>0.236 (0.179-0.293)</td></tr><tr><td>ETR-E</td><td>0.111 (0.067-0.154)</td><td>0.111 (0.067-0.154)</td></tr><tr><td>Spoligotypes (43 spacers)</td><td>0.930 (0.928-0.932)</td><td>0.930 (0.928-0.932)</td></tr></tbody></table>

Europe (male/female sex ratio=0.94). Regarding the distribution of lineages between human and animal hosts (Supplementary File 1/Supplementary Figure S1), BOV_1 sublineage was globally more common among human hosts (n = 388/473; 82.03%), followed by BOV sublineage (n = 72/473; 15.22%). Considering the diverse sub-regions, BOV sublineage that predominated among animal isolates in Central America (around 88% of isolates) was also largely present among patients from Northern Europe (around 36% of isolates). Lastly, the BOV_3 sublineage was merely visible in patients from Eastern Asia (n = 2) and South America (n = 1). Table 1 shows the distribution of predominant spoligotypes in SITVITBovis.

Globally, the four genotyping methods (spoligotyping, 5-locus ETRs and 12- and 15-loci MIRU-VNTRs) showed good discriminatory power (HGDI in range of 0.930-0.979; Table 2). Regarding the individual VNTR loci, we noted that QUB-11b (VNTR2163), ETR-A (VNTR2165), MIRU16 (VNTR1644) and MIRU04 (VNTR0580) displayed reasonable discriminatory power with HGDI>0.6 (Table 2).

Spoligoforest trees (Supplementary File 2) were drawn to highlight diversity of main spoligotypes (represented by a SB number or a SIT) in Africa (n = 965), The Americas (n = 1662), Asia (n = 289), Australasia (n = 105), and Europe (n = 22701). In these trees, we can visualize the main spoligotypes in each continent. One may notice that, according to our database, SB0121/SIT481 was more predominant in Europe (n = 5083 isolates), followed by SB0134/SIT665 (n = 2021), and SB0140/SIT683 (n = 1696); while SB0140/SIT683 (n = 42) predominated in Australasia, followed by SB0130/SIT691 (n = 31); SB0120/SIT482 was more prevalent in Asia (n = 179), followed by SB0140/SIT683 (n = 85) and SB0130/SIT691 (n = 32); SB0140/SIT683 was more predominant in the Americas (n = 413), followed by SB0121/SIT481 (n = 154), and SB0130/SIT691 (n = 128); finally, SB0944/SIT1037 was more predominant in Africa (n = 136), followed by SB0121/SIT481 (n = 99) and SB1176/SIT3453 (n = 59). Besides, a supplementary file was added in order to allow comparisons between different regions (Supplementary File 3).

## Proposal of an international consensus schema based on MIRU-VNTR loci

Several studies have reported affordable and suitable sets of VNTR loci to discriminate bTB isolates (31–33). However, most of these studies have been based on specific countries. In SITVITBovis database, we have gathered data from various countries, providing a wider view of bTB diversity. According to pooled data in our database, an optimal consensus set could be based on 13 MIRU-VNTR loci showing reasonably high discriminatory power (HGDI>0.4), i.e. QUB-11b (VNTR2163), ETR-A (VNTR2165), MIRU16 (VNTR1644), MIRU04 (VNTR0580), ETR-B (VNTR2461), QUB-4156 (VNTR4156), MIRU31 (VNTR3192), MIRU10 (VNTR0960), ETR-C (VNTR0577), QUB-26 (VNTR4052), Mtub21 (VNTR1955), Mtub04 (VNTR0424) and Mtub30 (VNTR2401) (Table 2).

In addition to the aforementioned VNTRs, we also recommend to include QUB3232 (VNTR3232) as suggested earlier by other studies (31, 32, 34), to further increase the discriminatory power in bTB diversity analyses.

<footnote>*Common loci between 5-locus ETRs and 12- and 15-loci MIRUs.</footnote>
<table><caption>**Table 3.** Comparison of distributions of age groups among isolates of patients infected by *BOV_1* sublineage as compared to other predominant sublineages found in SITVIT2 database</caption><thead><tr><th>Lineage</th><th>Number of strains</th><th>%0-20 yrs</th><th>%21-40 yrs</th><th>%41-60 yrs</th><th>%>60 yrs</th></tr></thead><tbody><tr><td>Beijing</td><td>2640</td><td>9.51</td><td>39.92</td><td>22.20</td><td>28.37</td></tr><tr><td>BOV_1</td><td>162</td><td>22.22</td><td>25.31</td><td>15.43</td><td>37.04</td></tr><tr><td>CAS1-Delhi</td><td>1085</td><td>18.53</td><td>50.05</td><td>22.30</td><td>9.12</td></tr><tr><td>EAI5</td><td>393</td><td>12.21</td><td>51.15</td><td>26.46</td><td>10.18</td></tr><tr><td>H3</td><td>1369</td><td>11.91</td><td>37.84</td><td>22.50</td><td>27.76</td></tr><tr><td>LAM9</td><td>1333</td><td>11.78</td><td>46.59</td><td>25.58</td><td>16.05</td></tr><tr><td>T1</td><td>2641</td><td>11.09</td><td>39.57</td><td>22.95</td><td>26.39</td></tr></tbody></table>

## Distribution of patient age groups

Comparison with SITVIT2 database showed that the distribution of patient age groups of *M. bovis* isolates (belonging to *BOV_1* sublineage) was different from age group distributions in patients infected with other MTBC members (*P*-value < 0.0001). When comparing the distribution of age groups 0–20 years, 21–40 years, 41–60 years and >60 years among patient isolates belonging to H3, LAM9, T1, Beijing, CAS1-Delhi and EAI5 sublineages, we noticed that patient isolates belonging to *BOV_1* sublineage were more common among younger (age group 0–20 years) as well as older patients (>60 years) (Table 3).

Considering only the bTB strains isolated in Europe (north, south and west, with most of the available data), we noted a significant disparity in the distribution of age groups (*P*-value < 0.0001). For all age groups (0–20 years, 21–40 years, 41–60 years and >60 years), we observed that the proportion of patients in Western Europe (EURO-W) was similar in each age group, representing around 20% (from 18 to 24%) of cases, except for the age group >60 years which represented about 35% of patients (Supplementary File 1/Supplementary Figure S2A). On the contrary, a majority of patients in Southern Europe (EURO-S) belonged to the age group 21–40 years (51% of patients), as opposed to the age group >60 years (81% of cases) in Northern Europe (EURO-N). These observations suggest that bovine tuberculosis is a major problem primarily affecting the working population in Southern Europe, while it rather concerns reactivation cases among senior patients in Northern Europe.

## Discussion

Results obtained from this study suggested an important heterogeneity in worldwide distribution of bTB isolates. Both significant cleavages and similarities were observed in the geographical distribution of various bTB sublineages defined in our database as well as for the various bTB genotypes.

Diversity, presence and/or absence of certain spoligotyping spacers have been revealed as being particularly useful for discriminating specific bTB isolates circulating worldwide (16, 13). Our web-based tool allows to have a global overview of bTB genotypes (including spoligotypes, E-locus ETRs and 12- and 15-loci MIRU-VNTRs). By using SITVITBovis tool, users can remarkably visualize the geographical position (at sub-region, country or city level) of specific genotypes in addition to available epidemiological data as well as drug resistance information; the latter being an important health issue in various regions of the world (35).

MIRU-VNTR method obviously increases the discriminatory power of spoligotyping; nonetheless, spoligotyping alone is still a reliable method to differentiate *M. bovis* isolates (Table 2). Several studies have focused on the diversity of MIRU-VNTRs and/or spoligotypes involved in bTB (31, 32). However, the majority of these studies generally focused on limited geographical areas.

In our analysis, spoligotyping was globally considered as a reasonable typing method to discriminate bTB isolates (with a global HGDI of 0.930). However, according to recent research, its discriminatory power could vary in function of countries (33). Despite an international consensus scheme of using MIRU-VNTR loci concomitantly with spoligotyping for improved characterization of MTBC (4–6, 15), only selected groups working on bTB isolates have used such an approach (31–34). Henceforth, the scarcity of MIRU-VNTR data in SITVITBovis should be improved by adding new data in future updates. Nonetheless, with data on 25 741 isolates (75 countries of isolation corresponding to 60 countries of origin), SITVITBovis still incorporates valuable data from other resources. Unfortunately, information on host is not always provided for all collected bTB isolates in published literature, which hampers efforts to precisely calculate the real incidence of bTB among human vs. other animal hosts. Due to the lack of available information, we did not address the issues of clonal complexes of *M. bovis* marked by specific regions of deletions (RDs) or single-nucleotide polymorphisms (SNPs) in the present study. Nevertheless, information on RDs may be obtained by querying the database and exporting results into an Excel file. Last but not least, WGS and SNPs are being frequently used today to decipher bTB transmission and identify lineages involved (36, 37). Studies using WGS already provide significant insights in bTB diversity, distribution and evolution (10–12) and may help identify hosts that could serve as reservoir of bTB infection and spread. As an illustration, a preliminary table linking WGS and classical genotyping data for some bTB isolates is available (Supplementary Table S3).

As the genomic information becomes more easily accessible in conjunction with geographical origin and hosts, we plan to automatically extract from sequencing data a number of variables such as specific phylogenetic markers, drug resistance, genomic determinants of virulence/pathogenicity or other relevant information. This novel generation of database(s) will be able to link the newly generated information with the existing data on classical genotyping information vs. hosts and geographical mapping. The development of all these tools in future bTB genotyping databases would possibly allow anticipating bTB outbreaks among humans and animal hosts alike.
## Conclusion and perspectives

In summary, our study provides a global overview of bTB distribution, highlighting potential relationships between bTB genotypes and affected hosts as well as other epidemiological features. SITVITBovis, with its correlation with other existing resources (such as Mbovis.org), provides a tentative picture of bTB circulation in the world and therefore represents a non-negligible resource for monitoring bTB and make correlations between epidemiological, macro-geographical and other data available.

Future developments aim to enrich the database with WGS data in conjunction with geographical origin and hosts. Development of automated extraction of knowledge related to specific phylogenetic markers, drug resistance, genomic determinants of virulence/pathogenicity or other relevant information will soon make it possible to link the newly generated information with the existing data generated using classical genotyping tools. The development of such a strategy in conjunction with information on hosts and subsequent geographical mapping will boost our efforts to generate a real comprehensive snapshot of global bTB diversity. Being helpful to better identify, treat and control bTB, such a tool will be highly beneficial for both medical and veterinarian specialists as well as public health authorities.

## Supplementary data

Supplementary data are available at Database Online.

## Acknowledgement

We thank more than 300 investigators who provided *M. bovis* data for SpoDB4, SITVITWEB and SITVIT2 databases. We are highly grateful to Igor Mokrousov, Thierry Zozio, Julie Millet, Véronique Hill and Elisabeth Streit for helpful discussions. D.C. was awarded a Ph.D. fellowship by the European Social Funds through the Regional Council of Guadeloupe.

## Competing interests

There is no competing interest.

## Author contributions statement

N.R. designed and supervised the study. D.C., I.C.M. and A.D. developed the database and designed the algorithms allowing to query the web application. Y.R. helped in the analysis of WGS data. D.C. and N.R. wrote and reviewed the manuscript.

## Funding

The work done at Institut Pasteur de la Guadeloupe was supported by a FEDER grant, financed by the European Union and Guadeloupe Region (Programme Opérationnel FEDER-Guadeloupe-Conseil Régional 2014–2020, Grant number 2015-FED-192).

## References

1. World Health Organization (WHO) Global tuberculosis report 2020. (2020), https://www.who.int/publications/i/item/9789240013131, (7 June 2021, accessed).

2. Smith,N.H., Hewinson,R.G., Kremer,K. et al. (2009) Myths and misconceptions: the origin and evolution of Mycobacterium tuberculosis. *Nat. Rev. Microbiol.*, 7, 537–544.

3. Kamerbeek,J., Schouls,L., Kolk,A. et al. (1997) Simultaneous detection and strain differentiation of Mycobacterium tuberculosis for diagnosis and epidemiology. *J. Clin. Microbiol.*, 35, 907–914.

4. Supply,P., Alix,C., Lesjean,S. et al. (2006) Proposal for standardization of optimized mycobacterial interspersed repetitive unit-variable-number tandem repeat typing of Mycobacterium tuberculosis. *J. Clin. Microbiol*, 44, 4498–4510.

5. Rastogi,N. and Sola,C. (2007) Chapter 2 – Molecular evolution of the Mycobacterium tuberculosis complex. Palomino J.C., Leao S. and Ritacco V. (ed), *Tuberculosis 2007: From Basic Science to Patient Care*. Vol. 2007, Amedeo Online Textbooks, pp. 53–91. http://pdf.flyingpublisher.com/tuberculosis2007.pdf.

6. García de Viedma,D., Mokrousov,I. and Rastogi,N. (2011) Innovations in the molecular epidemiology of tuberculosis. *Enferm. Infecc. Microbiol. Clin.*, 29, 8–13.

7. Demay,C., Liens,B., Burguière,T. et al. (2012) SITVITWEB – a publicly available international multimarker database for studying Mycobacterium tuberculosis genetic diversity and molecular epidemiology. *Infect. Genet. Evol.*, 12, 755–766.

8. Comas,I., Homolka,S., Niemann,S. et al. (2009) Genotyping of genetically monomorphic bacteria: DNA sequencing in Mycobacterium tuberculosis highlights the limitations of current methodologies. *PLoS One*, 4, e7815.

9. Pérez-Lago,L., Navarro,Y. and García-de-Viedma,D. (2013) Current knowledge and pending challenges in zoonosis caused by Mycobacterium bovis: a review. *Res. Vet. Sci.*, 97, S94–S100.

10. Lasserre,M., Fresia,P., Greif,G. et al. (2018) Whole genome sequencing of the monomorphic pathogen Mycobacterium bovis reveals local differentiation of cattle clinical isolates. *BMC Genomics*, 19, 2.

11. Zimpel,C.K., Patané,J.S.L., Guedes,A.C.P. et al. (2020) Global distribution and evolution of Mycobacterium bovis lineages. *Front Microbiol.*, 11, 843.

12. Rodrigues,R.A., Ribeiro Araújo,F., Rivera Dávila,A.M. et al. (2021) Genomic and temporal analyses of Mycobacterium bovis in southern Brazil. *Microbial. Genomics*, 7, 5.

13. Smith,N.H. and Upton,P. (2012) Naming spoligotype patterns for the RD9-deleted lineage of the Mycobacterium tuberculosis complex, www.Mbovis.org. *Infect. Genet. Evol.*, 12, 873–876.

14. Rodriguez-Campos,S., Schürch,A.C., Dale,J. et al. (2012) European 2-a clonal complex of Mycobacterium bovis dominant in the Iberian Peninsula. *Infect. Genet. Evol.*, 12, 866–872.

15. Couvin,D., David,A., Zozio,T. et al. (2019) Euro-geographical specificities of the prevailing tuberculosis epidemic as seen through SITVIT2, an updated version of the Mycobacterium tuberculosis genotyping database. *Infect. Genet. Evol.* 72, 31–43.

16. Rodriguez-Campos,S., González,S., de Juan,L. et al. (2012) Spanish network on surveillance monitoring of animal tuberculosis. A database for animal tuberculosis (mycoDB.es) within the context of the Spanish national programme for eradication of bovine tuberculosis. *Infect. Genet. Evol.* 12, 877–882.

17. Smith,N.H., Berg,S., Dale,J. et al. (2011) European 1: a globally important clonal complex of Mycobacterium bovis. *Infect. Genet. Evol.*, 11, 1340–1351.

18. Müller,B., Hilty,M., Berg,S. et al. (1951–1960) African 1, an epidemiologically important clonal complex of Mycobacterium bovis dominant in Mali, Nigeria, Cameroon, and Chad. *J. Bacteriol.*, 191, 2009.

19. Berg,S., Garcia-Pelayo,M.C., Müller,B. et al. (2011) African 2, a Clonal complex of Mycobacterium bovis epidemiologically important in East Africa. *J. Bacteriol.*, 193, 670–678.

20. Aranaz,A., Liébana,E., Mateos,A. et al. (1996) Spacer oligonucleotide typing of Mycobacterium bovis strains from cattle and other animals: a tool for studying epidemiology of tuberculosis. *J. Clin. Microbiol.*, 34, 2734–2740.
21. Frothingham,R. and Meeker-O'Connell,W.A. (1998) Genetic diversity in the Mycobacterium tuberculosis complex based on variable numbers of tandem DNA repeats. *Microbiology*, **144**, 1189-1196.

22. Xia,E., Teo,Y.Y. and Ong,R.T.H. (2016) SpoTyping: fast and accurate in silico Mycobacterium spoligotyping from sequence reads. *Genome Med.*, **8**, 19.

23. Tang,C.Y. and Ong,R.T. (1625-1626) MIRUReader: MIRU-VNTR typing directly from long sequencing reads. *Bioinformatics*, **36**, 2020.

24. Bankevich,A., Nurk,S., Antipov,D. et al. (2012) SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. *J. Comput. Biol.*, **19**, 455-477.

25. Couvin,D., Segretier,W., Stattner,E. et al. (2020) Novel methods included in SpolLineages tool for fast and precise prediction of Mycobacterium tuberculosis complex spoligotype families. *Database (Oxford)*, **2020**, baa108.

26. Phelan,J., O'Sullivan,D.M., Machado,D. et al. (2019) Integrating informatics tools and portable sequencing technology for rapid detection of resistance to anti-tuberculous drugs. *Genome Med.*, **11**, 41.

27. Freschi,L., Vargas Jr,R., Hussain,A. et al. (2020) *Population structure, Biogeography and Transmissibility of Mycobacterium tuberculosis*. *Nat. commun.*, **12**, 1-11.

28. Tang,C., Reyes,J.F., Luciani,F. et al. (2008) SpolTools: online utilities for analyzing spoligotypes of the Mycobacterium tuberculosis complex. *Bioinformatics*, **24**, 2414-2415.

29. Reyes,J.F., Francis,A.R., Tanaka,M.M. (2008) Models of deletion for visualizing bacterial variation: an application to tuberculosis spoligotypes. *BMC Bioinformatics*, **9**, 496.

30. Hunter,P.R. and Gaston,M.A (1988) Numerical index of the discriminatory ability of typing systems: an application of Simpson's index of diversity. *J. Clin. Microbiol.*, **26**, 2465-2466.

31. Duarte,E.L., Domingos,M., Amado,A. et al. (2010) MIRU-VNTR typing adds discriminatory value to groups of Mycobacterium bovis and Mycobacterium caprae strains defined by spoligotyping. *Vet. Microbiol.*, **143**, 299-306.

32. Rodriguez-Campos,S., Navarro,Y., Romero,B. et al. (2013) Splitting of a prevalent Mycobacterium bovis spoligotype by variable-number tandem-repeat typing reveals high heterogeneity in an evolving clonal group. *J. Clin. Microbiol.*, **51**, 3658-3665.

33. Navarro,Y., Herranz,M., Romero,B. et al. (2014) High-throughput multiplex MIRU-VNTR typing of Mycobacterium bovis. *Res. Vet. Sci.*, **96**, 422-425.

34. Lamine-Khemiri,H., Martínez,R., García-Jiménez,W.L. et al. (2014) Genotypic characterization by spoligotyping and VNTR typing of Mycobacterium bovis and Mycobacterium caprae isolates from cattle of Tunisia. *Trop. Anim. Health Prod.*, **46**, 305-311.

35. Vazquez-Chacon,C.A., Martínez-Guarneros,A., Couvin,D. et al. (2015) Human multidrug-resistant Mycobacterium bovis infection in Mexico. *Tuberculosis (Edinb.)*, **95**, 802-809.

36. Kohl,T.A., Utpatel,C., Niemann,S. et al. (2018) Mycobacterium bovis persistence in two different captive wild animal populations in Germany: a longitudinal molecular epidemiological study revealing pathogen transmission by whole-genome sequencing. *J. Clin. Microbiol.*, **56**, e00302-18.

37. Lipworth,S., Jajou,R., de Neeling,A. et al. (2019) SNP-IT tool for identifying subspecies and associated lineages of Mycobacterium tuberculosis complex. *Emerg. Infect. Dis.*, **25**, 482-488.