# Clonal expansion across the seas as seen through CPLP-TB database: A joint effort in cataloguing Mycobacterium tuberculosis genetic diversity in Portuguese-speaking countries

João Perdigão, Carla Silva, Jaciara Diniz, Catarina Pereira, Diana Machado, Jorge Ramos, Hugo Silva, Fernanda Abilleira, Clarice Brum, Ana J Reis, et al.

► To cite this version:

João Perdigão, Carla Silva, Jaciara Diniz, Catarina Pereira, Diana Machado, et al. Clonal expansion across the seas as seen through CPLP-TB database: A joint effort in cataloguing Mycobacterium tuberculosis genetic diversity in Portuguese-speaking countries. Infection, Genetics and Evolution, 2019, Molecular aspects of mycobacteril infections, 72, pp.44-58. 10.1016/j.meegid.2018.03.011. pasteur-01986661

HAL Id: pasteur-01986661

https://riip.hal.science/pasteur-01986661

Submitted on 19 Jan 2019

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
Research paper

# Clonal expansion across the seas as seen through CPLP-TB database: A joint effort in cataloguing *Mycobacterium tuberculosis* genetic diversity in Portuguese-speaking countries

João Perdigão<sup>a,*,1</sup>, Carla Silva<sup>a</sup>, Jaciara Diniz<sup>b</sup>, Catarina Pereira<sup>a</sup>, Diana Machado<sup>c</sup>, Jorge Ramos<sup>c</sup>, Hugo Silva<sup>a</sup>, Fernanda Abilleira<sup>b</sup>, Clarice Brum<sup>b</sup>, Ana J. Reis<sup>b</sup>, Maíra Macedo<sup>b</sup>, João L. Scaini<sup>b</sup>, Ana B. Silva<sup>b</sup>, Leonardo Esteves<sup>d</sup>, Rita Macedo<sup>e</sup>, Fernando Maltez<sup>f</sup>, Sofia Clemente<sup>g</sup>, Elizabeth Coelho<sup>h</sup>, Sofia Viegas<sup>i</sup>, Paulo Rabna<sup>j</sup>, Amabélia Rodrigues<sup>j</sup>, Nuno Taveira<sup>a,k</sup>, Luísa Jordao<sup>e</sup>, Afrânio Kritski<sup>l</sup>, José R. Lapa e Silva<sup>m</sup>, Igor Mokrousov<sup>n</sup>, David Couvin<sup>o</sup>, Nalin Rastogi<sup>o</sup>, Isabel Couto<sup>c</sup>, Arnab Pain<sup>p</sup>, Ruth McNerney<sup>q</sup>, Taane G. Clark<sup>r</sup>, Andrea von Groll<sup>b</sup>, Elis R. Dalla-Costa<sup>d</sup>, Maria Lúcia Rossetti<sup>d,s</sup>, Pedro E.A. Silva<sup>b</sup>, Miguel Viveiros<sup>c</sup>, Isabel Portugal<sup>a,*,1</sup>

## ARTICLE INFO

**Keywords:**
MIRU-VNTR
Spoligotyping
Tuberculosis
Drug resistance
Migration
Mycobacteria
LAM

## ABSTRACT

Tuberculosis (TB) remains a major health problem within the Community of Portuguese Language Speaking Countries (CPLP). Despite the marked variation in TB incidence across its member-states and continued human migratory flux between countries, a considerable gap in the knowledge on the *Mycobacterium tuberculosis* population structure and strain circulation between the countries still exists. To address this, we have assembled and analysed the largest CPLP *M. tuberculosis* molecular and drug susceptibility dataset, comprised by a total of 1447 clinical isolates, including 423 multidrug-resistant isolates, from five CPLP countries. The data herein presented reinforces Latin American and Mediterranean (LAM) strains as the hallmark of *M. tuberculosis* population structure in the CPLP coupled with country-specific differential prevalence of minor clades. Moreover, using high-resolution typing by 24-loci MIRU-VNTR, six cross-border genetic clusters were detected, thus supporting recent clonal expansion across the Lusophone space.

To make this data available to the scientific community and public health authorities we developed CPLP-TB

<footnote><sup>a</sup> iMed.ULisboa – Instituto de Investigação do Medicamento, Faculdade de Farmácia, Universidade de Lisboa, Lisboa, Portugal</footnote>

<footnote><sup>b</sup> Núcleo de Pesquisa em Microbiologia Médica, Faculdade de Medicina, Universidade Federal do Rio Grande, Rio Grande, Rio Grande do Sul, Brazil</footnote>

<footnote><sup>c</sup> Unidade de Microbiologia Médica, Global Health and Tropical Medicine, GHTM, Instituto de Higiene e Medicina Tropical, IHMT, Universidade Nova de Lisboa, UNL, Lisboa, Portugal</footnote>

<footnote><sup>d</sup> Centro de Desenvolvimento Científico e Tecnológico (CDCT), Porto Alegre, Brazil</footnote>

<footnote><sup>e</sup> Departamento de Doenças Infecciosas, Instituto Nacional de Saúde Dr. Ricardo Jorge, Lisboa, Portugal</footnote>

<footnote><sup>f</sup> Serviço de Doenças Infecciosas, Hospital de Curry Cabral, Lisboa, Portugal</footnote>

<footnote><sup>g</sup> Hospital da Divina Providência, Serviço de Doenças Infecciosas, Luanda, Angola</footnote>

<footnote><sup>h</sup> Programa Nacional de Controlo da Tuberculose, Ministério da Saúde de Moçambique, Mozambique</footnote>

<footnote><sup>i</sup> Instituto Nacional de Saúde, Ministério da Saúde de Moçambique, Mozambique</footnote>

<footnote><sup>j</sup> Instituto Nacional de Saúde Pública, Projecto de Saúde de Bandim (INASA/PSB), Bissau, Guinea-Bissau</footnote>

<footnote><sup>k</sup> Centro de Investigação Interdisciplinar Egas Moniz, Instituto Superior de Ciências da Saúde Egas Moniz, Monte de Caparica, Portugal</footnote>

<footnote><sup>l</sup> Academic Tuberculosis Program, School of Medicine, Federal University of Rio de Janeiro, Rio de Janeiro, Brazil</footnote>

<footnote><sup>m</sup> Thoracic Diseases Institute, Federal University of Rio de Janeiro, Rio de Janeiro, Brazil</footnote>

<footnote><sup>n</sup> Laboratory of Molecular Epidemiology and Evolutionary Genetics (former Laboratory of Molecular Microbiology), St. Petersburg Pasteur Institute, St. Petersburg, Russia</footnote>

<footnote><sup>o</sup> WHO Supranational TB Reference Laboratory, Tuberculosis and Mycobacteria Unit, Institut Pasteur de la Guadeloupe, Morne Jolivière Abymes, Guadeloupe, France</footnote>

<footnote><sup>p</sup> Pathogen Genomics Laboratory, BESE Division, King Abdullah University of Science and Technology (KAUST), Thuwal, Saudi Arabia</footnote>

<footnote><sup>q</sup> Lung Infection and Immunity Unit, UCT Lung Institute, University of Cape Town, Groote Schuur Hospital, Observatory, 7925, Cape Town, South Africa</footnote>

<footnote><sup>r</sup> London School of Hygiene & Tropical Medicine, Keppel Street, London WC1E 7HT, United Kingdom</footnote>

<footnote><sup>s</sup> Universidade Luterana do Brasil (ULBRA/RS), Porto Alegre, Brazil</footnote>

<footnote>* Corresponding authors at: Instituto de Investigação do Medicamento, Faculdade de Farmácia, Universidade de Lisboa, Av. Prof. Gama Pinto, 1649-003 Lisboa, Portugal. E-mail addresses: jperdigao@ff.ulisboa.pt (J. Perdigão), isabel.portugal@ff.ulisboa.pt (I. Portugal).</footnote>

<footnote n="1">¹ JP and IP are joint senior authors.</footnote>
(available at http://cplp-tb ff.ulisboa.pt), an online database coupled with web-based tools for exploratory data analysis. As a public health tool, it is expected to contribute to improved knowledge on the *M. tuberculosis* population structure and strain circulation within the CPLP, thus supporting the risk assessment of strain-specific trends.

# 1. Introduction

Tuberculosis (TB) is a devastating disease and a major public health issue, particularly among developing countries. The World Health Organization (WHO) declared TB as a global emergency in 1993 (World Health Organization, 1994). In 1997, new TB cases were estimated to ascend to approximately 8 million new cases per year and TB alone was responsible for 1.9 million deaths (Dye et al., 1999). Almost two decades later, the WHO estimates a total of 10.4 million new incident cases and 1.79 million deaths due to TB in 2015 (World Health Organization, 2016). The increasing trend in TB incidence shown by some countries and the fact that some of these had a larger incidence level than that previously recognized, have contributed to the massive number of TB cases we face today (World Health Organization, 2016).

The Community of Portuguese Speaking Countries (CPLP - https://www.cplp.org) comprises of nine member-states sharing a common trait – Lusophony (i.e., Portuguese-spoke language). Together, the member states account for 10,742,000 Km² (7.2%) of the world's land area and contain around 258,466,109 habitants (approx. 3.5% of the human population). These countries accounted for 348,600 (3.4%) of the new TB cases in 2015 and, together, they show an average incidence rate of 87.5 new cases per 100,000 habitants (World Health Organization, 2016). This high incidence rate is not uniform across all member-states, ranging from 23 (Portugal) to 551 (Mozambique) new cases per 100,000 habitants. Noteworthy, three of these countries (Angola, Brazil and Mozambique) are within the WHO's Top 30 High-Burden countries for TB (World Health Organization, 2016).

*Mycobacterium tuberculosis*, the etiological agent of TB, has shown a highly clonal population structure but even with the present state-of-the-art molecular fingerprinting methods or Whole Genome Sequencing, we still miss some of the links in these transmission chains and occasionally we fail to understand how *M. tuberculosis* is transmitted in the community (van Soolingen, 2014; Walker et al., 2013).

Still, the wide application of CRISPR-based spoligotyping, along with other genetic markers such as IS6110, IS1081 and VNTR loci for isolate fingerprinting has enabled the delineation of broader but distinct genetic lineages and their phylogeographical association and, lineages such as Latin American and Mediterranean (LAM) or the East African-Indian (EAI) were subsequently assigned characteristic spoligotyping profiles (Sola et al., 2001). However, spoligotyping interrogates specific spacer sequences within the Direct Repeat (DR) locus which evolves unidirectionally by independent loss of single spacers or through the deletion of contiguous repeats and spacers (Groenen et al., 1993). This unidirectional and bimodal mode of evolution may be difficult to interpret, model and, foremost, can lead to homoplasy which can encompass a non-negligible effect on strain clustering by spoligotyping alone (Comas et al., 2009; Reyes et al., 2012). In this sense, VNTR-based typing emerged as a genotyping method that is less prone to convergence by relying on the analysis of interspersed multi-allelic loci (Supply et al., 2006). Despite this, some of the spoligotyping main lineages or families were later confirmed and defined using new molecular markers such as the SNP barcode proposed by Coll et al. (2014). The prevalence of these lineages and clades at country level has gradually been made available through databases such as SITVITWEB or MIRU-VNTRplus, enabling the study of the *M. tuberculosis* population structure across different countries (Couvin and Rastogi, 2015; Demay et al., 2012; Weniger et al., 2010).

In the CPLP countries the body of knowledge available in the literature and databases has revealed that the *M. tuberculosis* population

structure contains a high predominance of LAM strains (Demay et al., 2012; Groenheit et al., 2011; Mokrousov et al., 2016; Perdigao et al., 2017; Saifodine et al., 2014; Viegas et al., 2010). However, an intense flow of people has historically existed between these countries and does still presently occur (http://migrationsmap.net). For this reason, additional tools developed for these specific countries are necessary to rapidly disseminate molecular epidemiological data and enable a higher resolution analysis of TB transmission.

The prevalence of drug resistance is highly variable across different CPLP member states. According to the WHO, rifampicin (RIF) resistant and multidrug resistant (MDR) TB, i.e. resistance to at least isoniazid (INH) and RIF, which is associated with poorer treatment outcome and mortality, may have ascended to 14,292 cases across the CPLP in 2015 (World Health Organization, 2016). Moreover, there is mounting evidence that the distribution of drug resistance is associated with specific phylogenetic clades and that specific within-country settings are associated with higher drug resistant TB rates than those currently estimated (Perdigao et al., 2017; Perdigao et al., 2010). The ability to track clones associated with such phylogenetic clades is therefore of the utmost importance in preventing the emergence of drug-resistant strains, including MDR- and XDR-TB, throughout the different countries.

Here, we combined different datasets from several member-states to provide a new framework for tracing specific *M. tuberculosis* strains across the Lusophone space and developed CPLP-TB, an online database aimed at disseminating molecular epidemiological data.

# 2. Methods

## 2.1. Clinical isolates and datasets

Five distinct datasets (Angola, n = 89; Brazil, n = 964; Guinea-Bissau, n = 13; Mozambique, n = 14; and, Portugal, n = 367) totalling 1447 *M. tuberculosis* clinical isolates (all of which recovered from different patients) were studied:

* Angolan dataset (n = 89): comprised by clinical isolates recovered from sputum samples obtained between March to June 2015 from the Hospital of Divina Providência in Luanda which, serves an estimated population of 990,892 inhabitants. This sample comprises approximately 18% of this hospital's yearly diagnosed cases;

* Brazil dataset (n = 964): this dataset is comprised by isolates from Rio Grande recovered at University Federal do Rio Grande (n = 624, 2006-2016) covering 30% of all TB cases and, 100% of the TB/HIV cases in the city of Rio Grande in this period; and, by MDR-TB/INH-resistant isolates from Porto Alegre and surrounding metropolitan area obtained from the LACEN - Rio Grande do Sul (n = 340, 2010-2014), estimated to cover 50% of such cases in the same period;

* Guinea-Bissau (n = 13): includes all MDR-TB strains collected in Bissau, Guinea-Bissau, at Hospital Raoul Follereau, during the period from May to December 2012;

* Mozambique (n = 14): this dataset includes DNA obtained from isolates recovered at Laboratório Nacional de Referência da Tuberculose from patients diagnosed at Hospital Central de Maputo between January to March 2012;

* Portugal (n = 367): includes isolates recovered from nine distinct laboratories and hospitals across Lisbon Health Region between 1995 and 2016. Most cases (n = 349) were obtained between 2007 and 2015 and include 129 MDR-TB cases that are estimated to
correspond to 65% of all the MDR-TB cases reported nationwide for the same period.

The inclusion criterion for this study was the availability of genotypic data: MIRU-VNTR (either 15 or 24 loci) and/or spoligotyping. All methods were performed in accordance with the relevant guidelines and regulations.

## 2.2. Drug susceptibility testing (DST)

DST, when available, was carried out for first-line drugs using the standardized procedure of the BACTEC™ MGIT™ 960 System (Becton Dickinson Diagnostic Systems, Sparks, MD, USA), except for a subset (n = 283/964) of the isolates from Brazil for which phenotypic drug susceptibility was assessed by the resazurin microtiter assay (REMA) (Palomino et al., 2002). Phenotypic data on pyrazinamide (PZA) susceptibility was only available for 90 out of the 713 isolates tested for first-line drug susceptibility in Brazil since PZA susceptibility is not routinely performed in this country. For this same reason, clinical isolates recovered in Portugal before 2000 also lack data on PZA susceptibility.

DST data for second-line drugs (SLDs) was available for 48, 9 and 98 MDR-TB clinical isolates from Brazil, Guinea-Bissau and Portugal, respectively. For this study, only phenotypic susceptibility data for fluoroquinolones (FQ; ofloxacin, ciprofloxacin or moxifloxacin) and second-line injectable drugs (SLIDs; kanamycin, amikacin or capreomycin) were considered.

Clinical isolates susceptible to all first-line drugs tested were classified as pan-susceptible isolates; resistance to at least isoniazid (INH) and rifampicin (RIF) was classified as MDR; any resistance pattern including resistance to one or more antituberculous drugs other than INH and RIF combined were classified as non-MDR; MDR-TB isolates showing cumulative resistance with any FQ and a second-line injectable drug were classified as extensively drug resistant (XDR). MDR-TB isolates lacking data on SLD susceptibility were, nonetheless, classified as MDR-TB.

## 2.3. MIRU-VNTR typing

Fifteen- or 24-loci MIRU-VNTR typing was carried out using the single or multiplex amplification procedure described by Supply et al. (2006). Amplicon size determination was carried out by agarose gel electrophoresis or capillary electrophoresis as described previously (Supply, 2005). Allelic diversity (h) for each locus was calculated using the equation $h = 1 - \sum x_i^2$, where $x_i$ is the relative frequency of the ith allele (Sun et al., 2004). A MIRU-VNTR dendrogram was constructed using R and ape package for phylogenetic analysis (Paradis et al., 2004). A MIRU-VNTR cluster was defined as a group of more than one isolate sharing indistinguishable MIRU-VNTR profiles and a cross-border cluster defined as a MIRU-VNTR cluster where all clinical isolates do not originate from the same country. Isolates showing two or more loci with dual allele profiles (assessed by dual peak or dual banding upon capillary/gel electrophoresis) were removed from the study.

## 2.4. Spoligotyping

Spoligotyping was performed by a single-tube multiplex PCR amplification of 43 spacer regions of the direct repeat (DR) locus using oligonucleotide primers DRa (5'-Biotin-GGTTTTGGGTCTGACGAC-3') and DRb (5'-CCGAGAGGGGACGAAAC-3') and 20 ng of genomic DNA. Amplicons were reverse hybridized on a membrane with amino-linked immobilized probes for each spacer as described previously by Kamerbeek et al. (1997). Detection was performed using the ECL® Chemiluminescence Detection System (GE Healthcare®, Cleveland, OH, USA) as per the manufacturer's instructions.

Spoligotyping profiles were assigned to lineage, clade and shared international type (SIT) using the rules described in SITVITWEB and SITVIT2 international databases (Couvin and Rastogi, 2015; Demay et al., 2012).

SIT assignment was carried out by entering the data in an updated proprietary version of SITVIT (SITVITEXTEND, housed at Institut Pasteur de la Guadeloupe). Minimum Spanning Trees (MST) were drawn using MLVA Compare software (Genoscreen, Lille, France and Ridom Bioinformatics, Münster, Germany).

## 2.5. Database development

CPLP-TB user interface was developed and implemented with Shiny by Rstudio (http://shiny.rstudio.com). Data processing, map and tree plots were implemented using the following R packages: ggplot2, ape, plotGoogleMaps, gtree and Cairo (Paradis et al., 2004). Interactive choropleth mapping of clinical isolates by country of origin was carried out using Google Maps API. Source and data processing code is available upon request from the corresponding author.

The web application contains a left-sided panel that enables the user to apply filters (e.g. year interval, country of origin, drug resistance, SIT/Clade, MIRU-VNTR cluster, etc) to the data being displayed on a main panel. This main panel is automatically updated according to the filtering criteria set on the left panel, and features different page tabs: a data table that lists all isolates' records deposited in this database along with phenotypic and genotypic data, allowing the user to browse through the different clinical isolates and apply additional filters to the data table; different plots that summarize data contained in the database, e.g., barplots showing the distribution per year, resistance type or SITs/clade; an interactive Google Maps API-based choropleth map showing the distribution of isolates per country; MIRU-VNTR and spoligotyping dendrograms; and, a genome-wide SNP-based phylogenetic tree containing the isolates present in the database that have already been subjected to Whole Genome Sequencing, annotated with the different SNP barcode sub-lineages previously identified by Coll et al. (2014). In its first and current release, CPLP-TB accommodates the datasets described in the present study.

## 2.6. Statistical analysis

All statistical analyses were conducted using the IBM® SPSS® Statistics v.21 (IBM Corporation, Armonk, NY, USA) or R (R Foundation for Statistical Computing, Vienna, Austria).

# 3. Results and discussion

## 3.1. Dataset description: country of origin, isolation years and drug resistance

Through a collaborative effort, the largest of this kind aimed at cataloguing the genetic diversity of *M. tuberculosis* across CPLP countries, this study brought together five distinct datasets comprised by *M. tuberculosis* clinical isolates recovered in five distinct CPLP countries: Angola (n = 89); Brazil (n = 964); Guinea-Bissau (n = 13); Mozambique (n = 14); and, Portugal (n = 367). These clinical isolates comprise a convenience sample collected from several laboratories and hospital units (see Methods for further details) between 1995 and 2016 (Table 1). The distribution per year of these clinical isolates and the number of drug resistant isolates to first-line drugs is summarized in Table 1. Briefly, from a total 1447 clinical isolates included in this study, 1182 (81.7%) were subjected to first-line DST, from which 423 (35.8%) were classed as MDR-TB isolates. Despite the availability of FQ and SLID DST data for 48 and 9 MDR-TB isolates from Brazil and Guinea-Bissau, respectively, XDR-TB isolates were only detected in the Portuguese dataset (n = 40) comprising 40.8% of the MDR-TB isolates tested for susceptibility to FQ and SLID in this country. This rate agrees
<table><caption>**Table 1**
Sample characteristics: datasets included in the present study by country of origin and stratification of the number of isolates by isolation year, susceptibility testing and drug resistance.</caption><thead><tr><th rowspan="2">Country of origin</th><th colspan="12">Isolation year</th><th colspan="12">No. of isolates with DST (first/s line)</th></tr><tr><th>1995</th><th>1996</th><th>1997</th><th>1998</th><th>2003</th><th>2005</th><th>2006</th><th>2007</th><th>2008</th><th>2009</th><th>2010</th><th>2011</th><th>2012</th><th>2013</th><th>2014</th><th>2015</th><th>2016</th><th>INH</th><th>RIF</th><th>STR</th><th>EMB</th><th>PZA</th><th>MDR</th><th>XDR</th></tr></thead><tbody><tr><td>Angola (n = 89)</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>89/0</td><td>17 (19.1)</td><td>6 (6.7)</td><td>13 (14.6)</td><td>5 (5.6)</td><td>5 (5.6)</td><td>5 (5.6)</td><td>-</td></tr><tr><td>Brazil (n = 964)</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>50</td><td>65</td><td>37</td><td>38</td><td>71</td><td>123</td><td>202</td><td>161</td><td>149</td><td>46</td><td>22</td><td>713/48</td><td>386 (54.1)</td><td>275 (38.6)</td><td>67 (15.6)</td><td>29 (6.7)</td><td>4 (4.4)</td><td>269 (37.7)</td><td>0 (0)</td></tr><tr><td>Guinea-Bissau (n = 13)</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>13</td><td>0</td><td>0</td><td>0</td><td>0</td><td>13/9</td><td>9 (69.2)</td><td>9 (69.2)</td><td>6 (46.2)</td><td>7 (53.8)</td><td>7 (53.8)</td><td>9 (69.2)</td><td>0 (0)</td></tr><tr><td>Mozambique (n = 14)</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>14</td><td>0</td><td>0</td><td>0</td><td>0/0</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Portugal (n = 367)</td><td>1</td><td>1</td><td>5</td><td>2</td><td>2</td><td>3</td><td>3</td><td>36</td><td>91</td><td>114</td><td>42</td><td>35</td><td>10</td><td>4</td><td>7</td><td>10</td><td>1</td><td>367/98</td><td>174 (47.4)</td><td>67 (15.6)</td><td>150 (40.9)</td><td>97 (26.4)</td><td>112 (31.2)</td><td>140 (38.1)</td><td>40 (40.8)</td></tr><tr><td>Total (n = 1447)</td><td>1</td><td>1</td><td>5</td><td>2</td><td>2</td><td>3</td><td>53</td><td>101</td><td>128</td><td>152</td><td>113</td><td>158</td><td>225</td><td>179</td><td>245</td><td>56</td><td>23</td><td>1182/155</td><td>586 (49.6)</td><td>434 (36.7)</td><td>236 (26.3)</td><td>138 (15.4)</td><td>128 (23.2)</td><td>423 (35.8)</td><td>40 (25.8)</td></tr></tbody></table>

with previous data on this setting and is the cause of major concern regarding ongoing extension of drug resistance to third-line drugs (Perdigao et al., 2013; Perdigao et al., 2016). In the Brazil dataset, no XDR-TB isolate was detected which, is not unexpected given the low number of isolates with data for second-line drug susceptibility testing (n = 48) and is in line with the low rate of XDR-TB estimated by the WHO for Brazil (5.9%) (World Health Organization, 2016).

The convenience-nature of the sampling constitutes a limitation to this study, but this is not an untypical design in a pathogen genetic setting, where clinically and epidemiological relevant isolates are characterized. Isolate inclusion based on availability of genotypic data is a common criterion due to their association with drug resistance. The prevalence of drug resistant strains does not necessarily reflect their local prevalence, except in the Angolan dataset, which was systematically collected over a defined period to assess drug resistance prevalence in a specific study population (Perdigao et al., 2017). Nonetheless, the datasets herein combined comprise the basis for the first attempt to investigate the global molecular epidemiology within the CPLP, enabling the investigation of possible links, and, provide a considerable progress to an expanding body of knowledge concerning *M. tuberculosis* genetic diversity in its member-states. It is, however, necessary to bear in mind potential bias in low sampled settings such as Mozambique and Guinea-Bissau and, furthermore, we stress the importance of future prospective studies that can play a key role in obtaining increasingly robust estimates on drug resistance in countries that do not carry out DST on a routine basis and, evaluate the impact of control programme interventions. Such studies should ideally encompass a uniform DST methodology, both country-wide and CPLP-wide, based on rapid nonradiometric liquid culture (Bemer et al., 2002). Currently, lack of routine culture-based DST poses an enormous limitation to this study and such future studies. Country policy on culture-based DST is not uniform across the CPLP where: Portugal performs it routinely for all culture-positive cases; CPLP African countries due to the high number of cases and lack of adequate laboratory support do not perform it routinely; and, in Brazil, country policy mandates routine culture-based DST to be carried out mainly for HIV co-infection cases, relapse, retreatment or treatment default cases. Herein, a subset of the samples obtained from Brazil were subjected to DST using REMA which should also be considered when analysing the data, albeit the good correlation between both methods on the evaluation of susceptibility to first-line drugs (Coban et al., 2014; Palomino et al., 2002).

## 3.2. Assessing *M. tuberculosis* population structure diversification through spoligotyping

To uncover the *M. tuberculosis* population structure in a country-wise manner, the combined dataset included spoligotyping data for 1239 clinical isolates, which corresponded to: 89 (7.2%) isolates from Angola; 760 (61.3%) from Brazil; 13 (1.0%) from Guinea-Bissau; 14 (1.1%) from Mozambique; and, 363 (29.3%) from Portugal. A total of 248 spoligotyping profiles were obtained, which corresponded to 164 SITs and 62 profiles classified as Orphan (Supplementary Tables S1 and S2). Additionally, 22 new SITs were identified upon data submission to SITVITEXTEND (SIT4139–4160). Overall, the most prevalent SITs found were SITs 20/LAM1 (n = 157; 12.7%), 42/LAM9 (n = 94/7.6%) and 53/T1 (n = 92; 7.4%) while the most prevalent clades encountered were the T1 (n = 211; 17.0%), LAM1 (n = 205; 16.5%), LAM9 (n = 141; 11.4%), H1 (n = 90; 7.3%), H3 (n = 83; 6.7%) and LAM4 (n = 67; 5.4%).

Comparing the distribution of the spoligotyping profiles by country, only the Portuguese dataset included isolates from all major *M. tuberculosis* lineages (Table 2). Remarkably, no Beijing strain was encountered on the Brazilian dataset despite its the large sample size. This is most likely correlated with the low prevalence of Beijing strains in this country as only 10 (0.3%) out of 2992 strains deposited in
<caption>**Table 2**</caption>

– Isolates characterized by spoligotyping, number of obtained SITs and number of isolates classified in the main spoligotyping lineages per dataset.

<table><thead><tr><th rowspan="2">Country</th><th rowspan="2">No. of isolates</th><th rowspan="2">No. of SITs</th><th rowspan="2">No. of orphan profiles</th><th colspan="8">Spoligotyping main lineages - no. of isolates (%)</th></tr><tr><th>Beijing</th><th>CAS</th><th>EAI</th><th>H</th><th>LAM</th><th>S</th><th>T</th><th>X</th></tr></thead><tbody><tr><td>Angola</td><td>89</td><td>25</td><td>8</td><td>0 (0)</td><td>0 (0)</td><td>0 (0)</td><td>1 (1.1)</td><td>55 (61.8)</td><td>0 (0)</td><td>31 (34.8)</td><td>0 (0)</td></tr><tr><td>Brazil</td><td>760</td><td>90</td><td>81</td><td>0 (0)</td><td>0 (0)</td><td>1 (0.1)</td><td>152 (20)</td><td>326 (42.9)</td><td>8 (1.1)</td><td>185 (24.3)</td><td>13 (1.7)</td></tr><tr><td>Guinea-Bissau</td><td>13</td><td>5</td><td>0</td><td>8 (61.5)</td><td>0 (0)</td><td>0 (0)</td><td>0 (0)</td><td>2 (15.4)</td><td>0 (0)</td><td>0 (0)</td><td>0 (0)</td></tr><tr><td>Mozambique</td><td>14</td><td>6</td><td>4</td><td>0 (0)</td><td>0 (0)</td><td>2 (14.3)</td><td>3 (21.4)</td><td>2 (14.3)</td><td>2 (14.3)</td><td>2 (14.3)</td><td>0 (0)</td></tr><tr><td>Portugal</td><td>363</td><td>77</td><td>13</td><td>26 (7.2)</td><td>4 (1.1)</td><td>4 (1.1)</td><td>18 (5)</td><td>226 (62.3)</td><td>3 (0.8)</td><td>42 (11.6)</td><td>17 (4.7)</td></tr></tbody></table>

SITVITWEB belonged to the Beijing clade (Demay et al., 2012). The prevalence in Rio Grande do Sul state, from where the Brazilian dataset originates, is most likely even lower given the fact no Beijing isolate was identified in this study (Table 2). Previously reported data in the literature support this low prevalence of Beijing strains among the *M. tuberculosis* populational structure in Brazil (< 1%) and, similarly, across South American countries such as Argentina, Paraguay and Chile (Balcells et al., 2015; Demay et al., 2012; Gomes et al., 2012). An exception is Peru that has reported a prevalence of Beijing strains ranging

between 5.5 and 9.15% (7.4% on SITVITWEB) (Caceres et al., 2014; Demay et al., 2012; Sheen et al., 2013). In contrast, in Guinea-Bissau, despite the reduced sample size, 8 (61.5%) out of 13 isolates belonged to Beijing clade. Although this fact is not surprising since Beijing isolates were previously detected in this country, the number of isolates retrieved from such a low-sized sample suggests that Beijing isolates may have been introduced more recently and appear to be emerging in this setting, as previously suggested (Rabna et al., 2015). Such emergence of Beijing strains beyond what can be considered its normal

<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Minimum Spanning Tree (MST) based on all available spoligotyping patterns in this study (n = 1239 isolates). Node coloring is shown in function of lineages and thickness and/or shape of branches (continuous, dashed, dotted, black or gray) varied in function of spacer changes between patterns. The number of changes is indicated on branches.</figcaption></figure>
phylogeographical niche may be associated and being driven by (multi) drug resistance and/or HIV co-infection, a phenomenon previously reported in other settings (Lasunskia et al., 2010; Liu et al., 2018; Viegas et al., 2013).

A common denominator appears to be LAM and the “ill-defined” T clade, found to be present across all datasets with diverging prevalence, highlighting specific populational structures across all settings (Table 2). Analysis of a minimum spanning tree based on spoligotypes (Fig. 1) reveals that LAM lineage is particularly diverse in this study. Interestingly, while the prevalence of LAM strains in the Angolan and Portuguese datasets was similar, there was a difference in the distribution of the most prevalent SITs (> 5%) among LAM clades in the Portuguese and Angolan datasets (Table 3). Angola showed a prevalence of SIT42/LAM9 strains that was twice the one in the Portuguese dataset, while the latter had a higher prevalence of SIT20/LAM1 strains. Also, SIT1106/LAM4 strains were only detected in Portugal, stressing the endemicity of this clade and its most likely origin in this country (Supplementary Fig. S1). Although previously reported in Brazil, the number of SIT1106/LAM4 isolates herein included from Portugal exceeds by 1.7-fold the number of these strains present in the SITVIT database (Table 3). Moreover, SIT1106/LAM4 have been associated with M/XDR-TB isolates of the Q1 cluster in Portugal and its global importance has recently been highlighted due to the potential to extend its drug resistance profile towards third-line drugs (Perdigao et al., 2016).

Despite detected in Portugal and Brazil, EAI strains appear to be of particular epidemiological significance in Mozambique where, despite a low sample size, 4 out of 14 (28.6%) strains belonged to this lineage (Table 2). This high prevalence is likely related with the proximity to the Indian subcontinent (Gagneux et al., 2006). In fact, the newly described SIT4140 could be specific to Mozambique and although labelled as an Unknown clade due to the absence of specific spacers it could potentially belong to the EAI lineage. In Brazil, Haarlem strains (SIT45/H1 and SIT50/H3, mostly) are also noteworthy, possibly representing more recently introduced European strains brought by European settlers of diverse origins (e.g., Italy, Poland, Germany and Russia) over two centuries ago (Garzelli et al., 2010; Jagielski et al., 2010; Lari et al., 2007). Still in the Brazilian dataset, the prevalence of T strains appears to be mostly underpinned by two different clades: SIT53/T1 and SIT65/T1. While SIT53/T1 strains are well dispersed through CPLP countries and was detected in all datasets except Guinea-Bissau, SIT65/T1 shows a much more restricted distribution and is in fact more prevalent in Brazil, with our data supporting a much higher prevalence in Rio Grande do Sul state where we obtained a prevalence of 7.9% (Table 3). Finally, 14 out of the 17 main SITs (prevalence ≥ 5%) detected in each country correlate well with the countries associated with at least 2% of the SIT's global distribution, thus strengthening the current understanding of *M. tuberculosis* phylogeography (Table 3).

### 3.3. A MIRU-VNTR framework for the CPLP

Of the 999 (69.0%) isolates genotyped by MIRU-VNTR, 181, 420 and 398 isolates were genotyped by 12, 15 and 24-loci MIRU-VNTR, respectively (Table 4). The MIRU-VNTR loci set used varied according to practice of the laboratories participating in this study. It was not possible to obtain data on MIRU-VNTR typing from the Mozambique dataset.

To understand the local transmission patterns and assess strain identity using a method with a higher discriminative power, we started by considering all clinical isolates genotyped by 15-loci MIRU-VNTR. The latter, can also be obtained from strains genotyped using the 24-loci set and has shown a comparative discriminatory power when compared with the 15-loci set. From the possible 818 MIRU-VNTR profiles available, 15 were excluded due to incomplete profiles, rendering a total of 802 15-loci MIRU-VNTR profiles suitable for analysis. Upon construction of a MIRU-VNTR dendrogram, we have identified

522 unique profiles and 75 genetic clusters encompassing a total of 356 (44.3%) clinical isolates (Supplementary Fig. S2). The remaining 447 (55.7%) isolates were classified as non-clustered (NC) isolates. Combining both country of origin and the genotypic profile generated by 15-loci MIRU-VNTR, 65 clusters were country-specific while the remaining 10 clusters were defined as cross-border clusters (CPLP clusters), containing between 2 and 17 isolates each from distinct countries. The largest CPLP cluster (CPLP-01) included 17 *M. tuberculosis* isolates from three distinct countries: Angola (n = 9), Brazil (n = 6) and Portugal (n = 2) (Supplementary Table S1).

However, despite the high discriminatory power (HGI = 0.993), MIRU-VNTR was not fully congruent with cluster associated SITs as 21 of 51 clusters that had two or more isolates showed mixed spoligotyping profiles (Supplementary Table S3). In nine of these, clustering incongruences cannot be simply explained by subsequent diversification at the DR locus level, suggesting convergence of MIRU-VNTR profiles.

Next, to compare with 24-loci MIRU-VNTR clustering, from the 398 isolates genotyped by this method, 17 failed to amplify at least one locus and a dendrogram was constructed for the remaining 381 isolates: 89 from Angola, 71 from Brazil, 13 from Guinea-Bissau and 208 from Portugal (Fig. 2, Supplementary Fig. S3). Using this method, we have identified 33 clusters encompassing a total of 152 (40.2%) clinical isolates (HGI = 0.986). The lower discriminatory power obtained was mostly due to the high prevalence of MDR- and XDR-TB clones in the Portuguese dataset and because, most of the susceptible isolates from Brazil were genotyped by 15-loci MIRU-VNTR whereas drug resistant isolates were more likely to have been genotyped by 24-loci MIRU-VNTR. Nevertheless, some clusters were further resolved leading to the definition of two additional clusters but supporting at least partial clustering of 32 of the original 75 genetic clusters defined using the 15-loci set (Supplementary Table S4). A subset of 6 clusters defined using the 24-loci set yields further support to clonal expansion at a transnational level. A MIRU-VNTR dendrogram including all isolates genotyped by 15 or 24-loci MIRU-VNTR is also included as Supplemental material (Supplementary Fig. S4).

Beyond the Lusophone space, and concerning the Beijing strains isolated from Guinea-Bissau, analysis by 24-loci allows the assignment of three isolates to type #94-32 (MLVA MtBC15-9 scheme), also termed as Russian/Asian clone CC1 (IS6110-RFLP-defined A0-cluster) (Merker et al., 2015; Mokrousov et al., 2008). These findings place the most likely origin of these strains in the former Soviet Union as type #94-32 is one of the two main Beijing types in Russia and highly prevalent in the former soviet Central Asia (Merker et al., 2015; Mokrousov et al., 2018; Mokrousov et al., 2008; Skiba et al., 2015).

The potential epidemiological links and similar population structures between present Lusophone countries may have its roots back to the 16th century with the historical overseas exploration period that brought together diverse cultures across distinct continents. More specifically, the similar population structure may have its common ground in the transatlantic trade routes and the African diaspora linking Portugal, Brazil and probably having Central Africa, a region that presently includes Angola, as its hub (Patarra and Fernandes, 2011). It is estimated that until 1850, slave trade introduced about 4 million African people in Brazil, most likely shaping the present *M. tuberculosis*

<table><caption>**Table 4**
Number of isolates genotyped by 12, 15 and 24-loci MIRU-VNTR.</caption><thead><tr><th rowspan="2">MIRU-VNTR set</th><th colspan="5">No. of isolates (%)</th><th rowspan="2">Total</th></tr><tr><th>Angola</th><th>Brazil</th><th>Guinea-Bissau</th><th>Mozambique</th><th>Portugal</th></tr></thead><tbody><tr><td>12-loci</td><td>0 (0)</td><td>181 (100)</td><td>0 (0)</td><td>0 (0)</td><td>0 (0)</td><td>181</td></tr><tr><td>15-loci</td><td>0 (0)</td><td>420 (100)</td><td>0 (0)</td><td>0 (0)</td><td>0 (0)</td><td>420</td></tr><tr><td>24-loci</td><td>89 (22.4)</td><td>83 (20.9)</td><td>13 (3.3)</td><td>0 (0)</td><td>213 (53.5)</td><td>398</td></tr></tbody></table>
<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Twenty-four loci MIRU-VNTR tree obtained for 381 isolates annotated with color-coded tip points that represent the country of origin (inner rim) and drug resistance (outer rim). Six cross-border clusters can be observed in the present tree as well as the association of XDR-TB isolates in Portugal with two main MDR-TB genetic clusters (Lisboa3-B and Q1).</figcaption></figure>

population structure of this country since this population may have
reached up to one-third of the total population (Graham, 2002; Prado
Jr., 2012). Furthermore, the importance of this population as an in-
coming TB reservoir is well documented from archaeological specimens
(Jaeger et al., 2013; Jaeger et al., 2012). However, indistinguishable
genetic clusters probably have its origins in more recent migratory
events, difficult to ascertain due to the complexity of the Lusophone
migratory system (Fig. 3), which some consider to be structured around
three main hubs: Angola, Brazil and Portugal (Baganha, 2009; Marques
and Góis, 2011). More recent periods of migration during the third, and
early fourth, quarter of the 20th century, after the independence of
Angola led to the return of approximately 500,000 settlers from former
Portuguese colonies to Portugal, mostly coming from Angola, and a
small minority to Brazil (Baganha, 2009). This mass migration might
have, as well, played a key role in shaping *M. tuberculosis* population
structure in recent times and underpin some of the cross-border clusters

here observed. Nevertheless, these considerations do not exclude the
ongoing migratory flows that take place up to this day and that might
also explain smaller or still undetected clusters. Also, the detection of
other unusual underrepresented clades are likely associated with more
recent events leading to the introduction and emergence of these strains
or, with additional migratory systems to which each country is prone to
be inserted in given the local geographical context (Goldblatt et al.,
2014; Mokrousov, 2016; Tolentino et al., 2011).

3.4. MIRU-VNTR loci and allelic diversities

Concerning the discriminatory capability of the different MIRU-
VNTR loci used, we calculated the allelic diversity (h) for each locus,
overall and across each country set, for the two sets of strains assessed
by 15 and 24-loci MIRU-VNTR (Fig. 4). MIRU-VNTR loci were classified
as poorly (h < 0.3), moderately (h = 0.3–0.6), and highly (h > 0.6)
<figure><img src="image_4.png" /><figcaption>Fig. 3. Lusophone migratory system depicting global migratory flows between the CPLP countries herein studied (orange arrows). Total migratory flow is indicated and proportionally represented by each arrow connecting these countries/regions. The percentage values indicated within parentheses represent the fraction of the total outward migratory flow from the country of origin. It is noteworthy the role played by Portugal as a central hub in this migratory system and the intense migratory events occurring between Portugal, Angola and Brazil. Part of the African diaspora is also represented in the present map as a major mass migratory event, taking place over four centuries, connecting the African continent with Brazil and leading to the introduction of at least 4,000,000 African slaves (yellow arrow). Another more recent and historical mass migratory event, with potential impact on the Portuguese *M. tuberculosis* populational structure in Portugal, occurred during the 70's as the result of the return of Portuguese settlers following the end of the Portuguese Colonial War (green arrow). Source: Trends in International Migrant Stock: Migrants by Destination and Origin. United Nations, Department of Economic and Social Affairs, 2015; figure generated using Microsoft® PowerPoint® 2016 [Version 1707] and Microsoft® Excel® 2016 (Version 1707), incl. Microsoft® Power Map 3D Data Visualization Tool (https://products.office.com/en-us/business/office). (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.)</figcaption></figure>

discriminatory. Globally, *loci* 2163b and 4052 showed the highest dis-
criminatory power (h = 0.77 and 0.81, respectively) while *loci* 580 and
2687 present the lowest allelic diversity (h = 0.07 and 0.05, respec-
tively). The allelic diversity scores found for this sample of isolates were
comparable to other published datasets comprised by clinical isolates of
diverse geographical origin, except for *loci* 580 and 2687 (Mokrousov
et al., 2004; Sola et al., 2003; Supply et al., 2006). These loci exhibited
very low allelic diversities across all datasets that strikingly contrast
with the ones obtained in the aforementioned studies. However, when
comparing with studies of a more geographically localized scope, the
allelic diversities obtained herein tend to be higher (Chen et al., 2012;
Cowan et al., 2002; Mazars et al., 2001; Mokrousov et al., 2004;
Shamputa et al., 2010; Supply et al., 2003; Wada et al., 2009).

3.5. Assigning phylogenetic defined sub-lineages to LAM strains

Given the epidemiological importance of the LAM lineage in the
CPLP, we attempted to assign Spoligotyping and Region of Difference

(RD)-based LAM-sublineages to the isolates genotyped by 24-loci
MIRU-VNTR as the tree obtained for these isolates is more robust from a
phylogenetic standpoint (Supplementary Fig. 2). The resulting tree was
examined for the distribution of specific MIRU-VNTR alleles previously
proposed by Mokrousov et al. (2016) as a proxy for some LAM-sub-
lineages (Mokrousov et al., 2016).

Within the LAM lineage, isolates bearing the 802_1 allele (locus 802,
number of repeats: 1) which was proposed to correlate with RD174
sublineage were found to comprise a monophyletic branch containing
184 isolates in which all, but 1 isolate, exhibit this specific allele. This
sub-lineage harbours the most important strains from an epidemiolo-
gical standpoint, specifically: the Lisboa3-A/B (SIT20/LAM1) or Q1
(SIT1106/LAM4) clusters in Portugal; RS-04 (SIT106/LAM7) cluster in
Brazil; or, the CPLP-1 cluster that is shared at a transnational level.
Supporting this, data on the genome-wide analysis of Lisboa3 and Q1
strains have already shown the association between these clusters and
RD174/RD<sup>RIO</sup> (Perdigao et al., 2014).

Thirty-one isolates were also found to form a monophyletic branch
<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Allelic diversities obtained for the sample subsets genotyped with 15 (A) and 24 (B) loci MIRU-VNTR. Allelic diversities for the combined datasets or for each individual dataset are shown. Horizontal dashed lines are used to highlight the breakpoints between low-to-intermediate and intermediate-to-high allelic diversities.</figcaption></figure>

containing the 154_1 allele as a marker for the RD115 sublineage. This branch included, among others, the “pinnipedii-like” strains SIT863/PINI2 which has been shown by some of us to in fact belong to the LAM lineage (Dalla Costa et al., 2015). Interestingly, the SIT863/PINI2 strains underwent diversification at locus 154 (154_2). The phylogenetic positioning of SIT863/PINI2 in the RD115-sublineage is congruent with previous data on its SNP barcode lineage (4.3.3) as a RD115 sub-lineage branch (Coll et al., 2014; Dalla Costa et al., 2015). Nevertheless, some phylogenetic incongruences were noted as other unrelated strains, e.g., Beijing or T1, were included in this branch, thus exemplifying the problem of MIRU-VNTR convergence.

The “SIT33-sublineage” was also identifiable among LAM strains: characterized by MIRU-VNTR alleles 2401_2 and 3171_1, a branch was identified as associated with these alleles. This specific branch was composed by 18 isolates of which 4 belonged to SIT33/LAM3 while four exhibited putative SIT33-derived spoligotyping profiles. Three isolates were, nevertheless, characterized by spoligotyping profiles apparently unrelated to SIT33. Furthermore, MDR-TB RS-04 cluster (SIT106/LAM7) was found to belong to this sub-lineage highlighting its importance in Brazil.

Additionally, we also examined the tree for allele 2059_1 which is proposed to be associated with the SIT388-sublineage, which was not found for any LAM isolate except for a branch within the RD115 sub-lineage. This marker has been reported as the least robust concerning its Positive Predictive Value but interestingly, this putative sub-branch comprised five isolates, aside an isolate with the 2059_2 allele, that included the SIT863/PINI2 isolates as the only strains bearing a profile compatible with descent from SIT388. However, SIT388-sublineage is proposed by Mokrousov et al. (2016) to be derived directly from SIT42/LAM9/RD115 strains and therefore constitute a Japan-endemic parallel branch to the RD115 sublineage consistent with its association with SNP barcode 4.3.1 (Coll et al., 2014; Mokrousov et al., 2016).

### 3.6. Drug resistance association with MIRU-VNTR clusters and spoligotyping clades

We have further analysed the distribution of drug resistant isolates across the MIRU-VNTR dendrograms and SITs/clades. Intersecting phenotypic drug resistance with the previously defined 24-loci MIRU-VNTR generated clusters, 30 where associated with drug resistant strains. From these, 22 and 4 comprised MDR-TB and XDR-TB strains, respectively. Looking further into the association of these clusters with drug resistance, it is noteworthy the statistically significant higher proportion of XDR-TB within Lisboa3, Q1 and Q1-derived PT-07 (SIT1106/LAM4) clusters in comparison to MDR-TB and/or strains exhibiting any type of drug resistance other than MDR/XDR in the same clusters (z-test, p < 0.05). In Lisboa3 cluster this association was also significant when comparing the higher proportion of MDR-TB strains with the proportion of strains exhibiting any type of drug resistance (z-test, p = 0.02).

Finally, the proportion of non-clustered pan-susceptible isolates was significantly higher than non-clustered isolates exhibiting MDR, XDR or any type of drug resistance (z-test, p < 0.01).

When analysing the dendrogram generated using 15-loci set, which enables the use of a larger subset of the isolates herein analysed, the statistical association is further extended to the clusters PT-03, showing a higher proportion of isolates with any type of drug resistance other than MDR/XDR relative to susceptible isolates (z-test, p = 0.019); and, cluster RS-05, showing a higher proportion of MDR-TB than susceptible isolates (z-test, p = 0.021).

This data provides further support to the epidemiological importance of Lisboa3 and Q1 phylogenetic clades that, despite the inclusion of other datasets from diverse countries, maintain a high degree of endemicity in Portugal and, calls for the adoption of specific control measures (e.g., active case finding) focused on the prevention of transmission of these highly resistant strains. Moreover, these two
<caption>**Fig. 5.** CPLP-TB database, an online tool purposely designed to house data on *M. tuberculosis* clinical isolates' original from the CPLP, therefore enabling the tracking of specific strains across the Lusophone space. The database features a web interface composed of a main panel with six different page tabs and a left-sided panel that enables users to apply filters to data displayed on the main panel. Data displayed on the main panel include a data-table with strain information such as drug resistance, MIRU-VNTR alleles; different data plots showing enable analysis of the selected data: a choropleth map, which enables users to visualize the distribution of clinical isolates; MIRU-VNTR and Spoligotyping dendrograms; and, a genome-wide phylogenetic tree containing the isolates already subjected to WGS and annotated with the SNP barcode nomenclature.</caption>

CPLP-TB

TB Molecular Epidemiological Database for the CPLP

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

CPLP-TB

TB Molecular Epidemiological Database for the CPLP

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

CPLP-TB

TB Molecular Epidemiological Database for the CPLP

<figure><img src="image_8.png" /></figure>
clades (Mtbc15-9 types: 5298-54 [Lisboa3-B] and 10,571-67 [Q1], respectively) have an increasing importance in the European context as these have been recently identified and shortlisted as MDR-TB associated European cross-border clusters, with likely origin in Portugal, through the detection in the United Kingdom (Q1, n = 1) and France (Lisboa3-B, n = 1).(European Centre for Disease Prevention and Control/WHO Regional Office for Europe, 2017).

These cluster level associations were also reflected by specific associations found between the spoligotyping clade and resistance type: a higher proportion of susceptible strains, in comparison to isolates displaying MDR, were found among the H1 (z-test, p < 0.001), LAM3 (z-test, p = 0.002) and LAM9 (z-test, p = 0.024) clades; MDR isolates were found in a higher proportion, in comparison to susceptible isolates, in H3 (z-test, p < 0.001), LAM5 (z-test, p < 0.001), LAM4 (z-test, p < 0.001) and PINI2 (z-test, p = 0.003) clades; and, XDR-TB isolates were found in a higher proportion among LAM1 (z-test, p < 0.001) and LAM4 (z-test, p < 0.001) isolates when compared to MDR or susceptible isolates.

### 3.7. CPLP-TB: a new web-based tool for strain tracking across the CPLP

To facilitate the dissemination and inter-laboratory comparison of the genotypic data generated in the present study we have created CPLP-TB, an online database and web application (available at http://cplp-tb ff.ulisboa pt/) (Fig. 5).

This database will accommodate more data on existing strains and, holds the potential to also accommodate more records on new strains and/or additional genotypic markers (such as SNPs or MIRU-VNTR hypervariable loci). Laboratories across the CPLP can also deposit data from additional settings and/or strains thus improving the comprehensiveness of this new database. With this database and web application it is also possible to track specific *M. tuberculosis* clones or monitor the *M. tuberculosis* population structure of a given country.

The main goal of this new tool is to provide a more comprehensive and global view of TB epidemiology in the CPLP. It holds several advantages in comparison to existing TB databases, which are intimately related with its main mission. CPLP-TB is intended to undergo through more rapid releases and continuous updates, including the addition of strains from underrepresented settings, that will enable to perform or contribute to risk assessment on specific strains at a more early stage; submission of data will be open to other laboratories provided that the strains deposited have been isolated within a CPLP country; it aggregates phenotypic and diverse genotypic data/markers and with the current architecture the database can be easily updated to accommodate additional epidemiological markers; it has been developed specifically for CPLP countries through CPLP-based laboratories, stimulating capacity development by local laboratories; and, CPLP-TB is linked to a biobank distributed across its hosting institution and its associated laboratories, thus enabling fundamental as well as applied research using clinically relevant *M. tuberculosis* strains.

### 3.8. Concluding remarks

The present study constitutes the first and largest collaborative effort to catalogue the *M. tuberculosis* diversity within the Lusophone space. A total of 1447 clinical isolates with origin in five countries, including 423 MDR-TB isolates, were analysed to provide a more global perspective of TB in the CPLP. The population structure covered by this analysis corroborates the high prevalence of LAM strains in most settings and its dispersion argues in favour of common populational structures, with past epidemiological links, driven by the radiation of this specific sub-lineage out of its setting of origin, which is presently unknown (Mokrousov et al., 2016).

Overall, this scenario argues in favour of the importance of global TB surveillance and it is in this context that CPLP-TB database can have a pivotal role in understanding shared transmission networks. In this

regard, 24-loci MIRU-VNTR was found to be particularly useful and superior to the 15-loci set in analysing such a diverse dataset of clinical isolates, as it enabled the discrimination of 15-loci MIRU-VNTR clusters, which was of particular importance among cross-border clusters. Specifically, SIT20/LAM1 strains belonging to a newly identified cluster – CPLP-01 – emerge as the most important strains associated with more recent transnational strain flow. Most of the strains from this cluster were identified in Angola, which, given the present epidemiological scenario, makes it the likely origin of this particular cluster. Nonetheless, despite known limitations of both spoligotyping and MIRU-VNTRs, these typing methods remain the most used, easily-available, robust, and highly affordable methods – which makes them a method of choice for TB molecular epidemiology studies, particularly for developing countries. WGS based genomic epidemiology is not yet widely used in such settings. Future endeavours might consist in building new generation of databases able to establish useful links between data generated by “newer” (WGS) and “older” (spoligotyping/MIRU-VNTRs) molecular typing methods. Such a promising development includes a new tool named “SpoTyping” which allows in silico determination of *M. tuberculosis* spoligotyping patterns with a high accuracy from NGS reads, and subsequent lineage determinations by interrogating the SITVIT resources (Xia et al., 2016). Such developments will allow that the vast data generated using older genotyping methods are not denigrated and are instead useful to highlight the distribution summaries of the meta-data corresponding to *M. tuberculosis* lineages available in CPLP-TB and SITVIT databases. As such, this new tool, along with the several analysis and visualization tools, holds the potential for the analysis of TB transmission in these settings with an unprecedented resolution and can thus comprise a new framework for the analysis of *M. tuberculosis* genetic diversity within the Lusophone context and programmatic TB control.

## Funding

Financial support was provided by the European Society of Clinical Microbiology and Infectious Diseases, for which we would like to would like to acknowledge the Study Group for Mycobacterial Infections; Fundação CAPES [PVE-CAPES. 88881.064961/2014-01- Jose R. Lapa e Silva/UFRJ coordinator]; Genotyping and susceptibility profile of *Mycobacterium tuberculosis* clinical isolates from Rio Grande, Brazil were funded by Apoio a Projetos de Pesquisa em Doenças Negligenciadas, Brazil/MCTI/CNPq/MS-SCTIE – Decit [404081/2012-6] and by Programa Pesquisa para o SUS – PPSUS - FAPERGS/MS/CNPq/SESRS [1193-2551/13-6]; MIRU-VNTR typing and spoligotyping of *Mycobacterium tuberculosis* clinical isolates from Porto Alegre, Brazil were funded by National Council of Research [CNPq/MCTI/Universal - Project number: 441499/2014-7]; Fundação para a Ciência e a Tecnologia (FCT) Portugal [PTDC/SAU-EPI/122400/2010], part of the EDCTP2 program supported by the European Union; Fundação Calouste Gulbenkian, Portugal [Project ref. P-99934]. JP was supported by a post doc fellowship from project [PTDC/SAU-EPI/122400/2010] and by fellowship [SFRH/BPD/95406/2013] from FCT. The phylogenetic analysis work at Nalin Rastogi's lab 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). IM was supported by Russian Science Foundation (grant 14-14-00292). AP was supported by a faculty baseline funding from KAUST [BAS/1/1020-01-01]. DM was supported by FCT fellowship [SFRH/BPD/100688/2014] and DM, IC MV are thankful to [GHTM UID/Multi/04413/20139] from FCT and to projects “ForDILAB-TB” and “A implementação de um novo método de identificação rápida do complexo M. tuberculosis nos Laboratórios de Referência da Tuberculose de Maputo e Beira” from Fundação Calouste Gulbenkian and the Community of the Portuguese Speaking Countries (CPLP). CS was supported by FCT [SFRH/BD/73579/2010]. TC is funded by the Medical Research Council UK (Grant no. MR/K000551/1
and MR/M01360X/1, MR/N010469/1, MC_PC_15103).

## Author contributions

JP, MV and IP designed and directed the study. JP, JD, ED, NT, SC, DM coordinated sample and data collection. JP, CS, JD, CP, DM, JR, HS, FA, CB, AR, MM, JS, AS, LE, RM, FM, SC, EC, SV, PR, AR, LJ, IC, AG, ED, MR, PS, MV and IP undertook sample collection, DNA extraction, genotyping and phenotypic drug susceptibility testing. JP, DC, NR and IM analysed the data and performed bioinformatic and statistical analyses. DC and NR performed comparative analysis against SITVIT. JP, AK, JLS, IM, DC, NR, ED, MR, PS, MV and IP interpreted results. JP wrote the ode of the online database. JP and IP wrote the first draft of the manuscript. TC, AP and RM generated preliminary whole-genome sequence data available online. JP, AK, JLS, NT, IM, DC, NR, ED, MR, PS, MV and IP commented and edited on various versions of the draft manuscript and all authors approved the manuscript. JP, MV and IP compiled the final manuscript.

## Additional information

The authors have no competing financial interests to declare.

## Appendix A. Supplementary data

Supplementary data to this article can be found online at https://doi.org/10.1016/j.meegid.2018.03.011.

## References

Baganha, M.I., 2009. The Lusophone migratory system: patterns and trends. Int. Migr. 47, 5-20.

Balcells, M.E., Garcia, P., Meza, P., Pena, C., Cifuentes, M., Couvin, D., Rastogi, N., 2015. A first insight on the population structure of *Mycobacterium tuberculosis* complex as studied by spoligotyping and MIRU-VNTRs in Santiago, Chile. PLoS One 10, e0118007.

Bemer, P., Palicova, F., Rusch-Gerdes, S., Drugeon, H.B., Pfyffer, G.E., 2002. Multicenter evaluation of fully automated BACTEC mycobacteria growth Indicator tube 960 system for susceptibility testing of *Mycobacterium tuberculosis*. J. Clin. Microbiol. 40, 150-154.

Caceres, O., Rastogi, N., Bartra, C., Couvin, D., Galarza, M., Asencios, L., Mendoza-Ticona, A., 2014. Characterization of the genetic diversity of extensively-drug resistant *Mycobacterium tuberculosis* clinical isolates from pulmonary tuberculosis patients in Peru. PLoS One 9, e112789.

Chen, Y.Y., Chang, J.R., Huang, W.F., Kuo, S.C., Su, I.J., Sun, J.R., Chiueh, T.S., Huang, T.S., Chen, Y.S., Dou, H.Y., 2012. Genetic diversity of the *Mycobacterium tuberculosis* Beijing family based on SNP and VNTR typing profiles in Asian countries. PLoS One 7, e39792.

Coban, A.Y., Deveci, A., Sunter, A.T., Palomino, J.C., Martin, A., 2014. Resazurin microtiter assay for isoniazid, rifampicin, ethambutol and streptomycin resistance detection in *Mycobacterium tuberculosis*: updated meta-analysis. Int. J. Mycobacteriol. 3, 230-241.

Coll, F., McNerney, R., Guerra-Assuncao, J.A., Glynn, J.R., Perdigao, J., Viveiros, M., Portugal, I., Pain, A., Martin, N., Clark, T.G., 2014. A robust SNP barcode for typing *Mycobacterium tuberculosis* complex strains. Nat. Commun. 5, 4812.

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

Couvin, D., Rastogi, N., 2015. Tuberculosis - a global emergency: tools and methods to monitor, understand, and control the epidemic with specific example of the Beijing lineage. Tuberculosis (Edinburgh, Scotland) 95 (Suppl. 1), S177-189.

Cowan, L.S., Mosher, L., Diem, L., Massey, J.P., Crawford, J.T., 2002. Variable-number tandem repeat typing of *Mycobacterium tuberculosis* isolates with low copy numbers of IS6110 by using mycobacterial interspersed repetitive units. J. Clin. Microbiol. 40, 1592-1602.

Dalla Costa, E.R., Vasconcelos, S.E., Esteves, L.S., Gomes, H.M., Gomes, L.L., da Silva, P.A., Perdigao, J., Portugal, I., Viveiros, M., McNerney, R., Pain, A., Clark, T.G., Rastogi, N., Unis, G., Rossetti, M.L., Suffys, P.N., 2015. Multidrug-resistant *Mycobacterium tuberculosis* of the Latin American Mediterranean lineage, wrongly identified as *Mycobacterium pinnipedii* (Spoligotype international type 863 [SIT863]), causing active tuberculosis in South Brazil. J. Clin. Microbiol. 53, 3805-3811.

Demay, C., Liens, B., Burguiere, T., Hill, V., Couvin, D., Millet, J., Mokrousov, I., Sola, C., Zozio, T., Rastogi, N., 2012. SITVITWEB—a publicly available international multi- marker database for studying *Mycobacterium tuberculosis* genetic diversity and mo- lecular epidemiology. Infect. Genet. Evol. 12, 755-766.

Dye, C., Scheele, S., Dolin, P., Pathania, V., Raviglione, M.C., 1999. Consensus statement. Global burden of tuberculosis: estimated incidence, prevalence, and mortality by country. WHO global surveillance and monitoring project. JAMA 282, 677-686.

European Centre for Disease Prevention and Control/WHO Regional Office for Europe, 2017. Molecular Typing for Surveillance of Multidrug-Resistant Tuberculosis in the EU/EEA – March 2017. European Centre for Disease Prevention and Control, Stockholm.

Gagneux, S., DeRiemer, K., Van, T., Kato-Meda, M., de Jong, B.C., Narayanan, S., Nicol, M., Niemann, S., Kremer, K., Gutierrez, M.C., Hilty, M., Hopewell, P.C., Small, P.M., 2006. Variable host-pathogen compatibility in *Mycobacterium tuberculosis*. Proc. Natl. Acad. Sci. U. S. A. 103, 2869-2873.

Garzelli, C., Lari, N., Cuccu, B., Tortoli, E., Rindi, L., 2010. Impact of immigration on tuberculosis in a low-incidence area of Italy: a molecular epidemiological approach. Clin. Microbiol. Infect. 16, 1691-1697.

Goldblatt, D., Rorman, E., Chemtob, D., Freidlin, P.J., Cedar, N., Kaidar-Shwartz, H., Dveyrin, Z., Mor, Z., 2014. Molecular epidemiology and mapping of tuberculosis in Israel: do migrants transmit the disease to locals? Int. J. Tuberc. Lung Dis. 18, 1085-1091.

Gomes, H.M., Elias, A.R., Oelemann, M.A., Pereira, M.A., Montes, F.F., Marsico, A.G., Kritsli, A.L., Filho Ldos, A., Caldas, P.C., Possuelo, L.G., Cafrune, R., Rossetti, M.L., Lucena, N., Saad, M.H., Cavalcanti, H.R., Leite, C.Q., de Brito, R.C., Lopes, M.L., Lima, K., Souza, M., Trindade Rde, C., Zozio, T., Sola, C., Rastogi, N., Suffys, P.N., 2012. Spoligotypes of *Mycobacterium tuberculosis* complex isolates from patients residents of 11 states of Brazil. Infect. Genet. Evol. 12, 649-656.

Graham, R., 2002. Nos tumeiros mais uma vez? o comércio interprovincial de escravos no Brasil. Afro-Asia 121-160.

Groenen, P.M., Bunschoten, A.E., van Soolingen, D., van Embden, J.D., 1993. Nature of DNA polymorphism in the direct repeat cluster of *Mycobacterium tuberculosis*: application for strain differentiation by a novel typing method. Mol. Microbiol. 10, 1057-1065.

Groenheit, R., Ghebremichael, S., Svensson, J., Rabna, P., Colombatti, R., Riccardi, F., Couvin, D., Hill, V., Rastogi, N., Koivula, T., Kallenius, G., 2011. The Guinea-Bissau family of *Mycobacterium tuberculosis* complex revisited. PLoS One 6, e18601.

Jaeger, L.H., Leles, D., Lima Vdos, S., da Silva Lda, P., Dias, O., Iniguez, A.M., 2012. *Mycobacterium tuberculosis* complex detection in human remains: tuberculosis spread since the 17th century in Rio de Janeiro, Brazil. Infect. Genet. Evol. 12, 642-648.

Jaeger, L.H., de Souza, S.M., Dias, O.F., Iniguez, A.M., 2013. *Mycobacterium tuberculosis* complex in remains of 18th-19th century slaves, Brazil. Emerg. Infect. Dis. 19, 837-839.

Jagielski, T., Augustynowicz-Kopec, E., Zozio, T., Rastogi, N., Zwolska, Z., 2010. Spoligotype-based comparative population structure analysis of multidrug-resistant and isoniazid-monoresistant *Mycobacterium tuberculosis* complex clinical isolates in Poland. J. Clin. Microbiol. 48, 3899-3909.

Kamerbeek, J., Schouls, L., Kolk, A., van Agterveld, M., van Soolingen, D., Kuijper, S., Bunschoten, A., Molhuizen, H., Shaw, R., Goyal, M., van Embden, J., 1997). Simultaneous detection and strain differentiation of *Mycobacterium tuberculosis* for diagnosis and epidemiology. J. Clin. Microbiol. 35, 907-914.

Lari, N., Rindi, L., Bonanni, D., Rastogi, N., Sola, C., Tortoli, E., Garzelli, C., 2007. Three-year longitudinal study of genotypes of *Mycobacterium tuberculosis* isolates in Tuscany, Italy. J. Clin. Microbiol. 45, 1851-1857.

Lasunskiaa, E., Ribeiro, S.C., Manicheva, O., Gomes, L.L., Suffys, P.N., Mokrousov, I., Ferrazoli, L., Andrade, M.R., Kritsli, A., Otten, T., Kipnis, T.L., da Silva, W.D., Vishnevsky, B., Oliveira, M.M., Gomes, H.M., Baptista, I.F., Narvskaya, O., 2010. Emerging multidrug resistant *Mycobacterium tuberculosis* strains of the Beijing genotype circulating in Russia express a pattern of biological properties associated with enhanced virulence. Microbes Infect. 12, 467-475.

Liu, H., Zhang, Y., Liu, Z., Liu, J., Hauck, Y., Dong, H., Zhao, X., Lu, B., Jiang, Y., Vergnaud, G., Pourcel, C., Wan, K., 2018. Associations between *Mycobacterium tuberculosis* Beijing genotype and drug resistance to four first-line drugs: a survey in China. Front Med. 12, 92-97.

Marques, J.C., Góis, P., 2011. A evolução do sistema migratório lusófono. Uma análise a partir da imigração e emigração portuguesa. Revista Internacional em Língua Portuguesa 213-231.

Mazars, E., Lesjean, S., Banuls, A.L., Gilbert, M., Vincent, V., Gicquel, B., Tibayrenc, M., Locht, C., Supply, P., 2001. High-resolution minisatellite-based typing as a portable approach to global analysis of *Mycobacterium tuberculosis* molecular epidemiology. Proc. Natl. Acad. Sci. U. S. A. 98, 1901-1906.

Merker, M., Blin, C., Mona, S., Duforet-Frebourg, N., Lecher, S., Willery, E., Blum, M.G., Rusch-Gerdes, S., Mokrousov, I., Aleksic, E., Allix-Beguec, C., Antierens, A., Augustynowicz-Kopec, E., Ballif, M., Barletta, F., Beck, H.P., Barry 3rd, C.E., Bonnet, M., Borroni, E., Campos-Herrero, I., Cirillo, D., Cox, H., Crowe, S., Crudu, V., Diel, R., Drobniewski, F., Fauville-Dufaux, M., Gagneux, S., Ghebremichael, S., Hanekom, M., Hoffner, S., Jiao, W.W., Kalon, S., Kohl, T.A., Kontsevaya, I., Lillebaek, T., Maeda, S., Nikolayevskyy, V., Rasmussen, M., Rastogi, N., Samper, S., Sanchez-Padilla, E., Savic, B., Shamputa, I.C., Shen, A., Sng, L.H., Stakenas, P., Toit, K., Varaine, F., Vukovic, D., Wahl, C., Warren, R., Supply, P., Niemann, S., Wirth, T., 2015. Evolutionary history and global spread of the *Mycobacterium tuberculosis* Beijing lineage. Nat. Genet. 47, 242-249.

Mokrousov, I., 2016. Emerging resistant clone of *Mycobacterium tuberculosis* in West Asia. Lancet Infect. Dis. 16, 1326-1327.

Mokrousov, I., Narvskaya, O., Limeschenko, E., Vyazovaya, A., Otten, T., Vishnevskiy, B., 2004. Analysis of the allelic diversity of the mycobacterial interspersed repetitive units in *Mycobacterium tuberculosis* strains of the Beijing family: practical implications and evolutionary considerations. J. Clin. Microbiol. 42, 2438-2444.

Mokrousov, I., Narvskaya, O., Vyazovaya, A., Millet, J., Otten, T., Vishnevsky, B., Rastogi, N., 2008. *Mycobacterium tuberculosis* Beijing genotype in Russia: in search of
informative variable-number tandem-repeat loci. J. Clin. Microbiol. 46, 3576-3584.

Mokrousov, I., Vyazovaya, A., Iwamoto, T., Skiba, Y., Pole, I., Zhdanova, S., Arikawa, K., Sinkov, V., Umpeleva, T., Valcheva, V., Alvarez Figueroa, M., Ranka, R., Jansone, I., Ogarkov, O., Zhuravlev, V., Narvskaya, O., 2016. Latin-American-Mediterranean lineage of *Mycobacterium tuberculosis*: human traces across pathogen's phylogeny. Mol. Phylogenet. Evol. 99, 133-143.

Mokrousov, I., Chernyaeva, E., Vyazovaya, A., Skiba, Y., Solovieva, N., Valcheva, V., Levina, K., Malakhova, N., Jiao, W.W., Gomes, L.L., Suffys, P.N., Kutt, M., Aitkhozhina, N., Shen, A.D., Narvskaya, O., Zhuravlev, V., 2018. A rapid assay for detection of the epidemiologically important Central Asian/Russian strain of the *Mycobacterium tuberculosis* Beijing genotype. J. Clin. Microbiol. 56 e01551-17.

Palomino, J.C., Martin, A., Camacho, M., Guerra, H., Swings, J., Portaels, F., 2002. Resazurin microtiter assay plate: simple and inexpensive method for detection of drug resistance in *Mycobacterium tuberculosis*. Antimicrob. Agents Chemother. 46, 2720-2722.

Paradis, E., Claude, J., Strimmer, K., 2004. APE: analyses of phylogenetics and evolution in R language. Bioinformatics 20, 289-290.

Patarra, N.L., Fernandes, D., 2011. Brasil: país de imigração? Revista Internacional em Língua Portuguesa 65-96.

Perdigao, J., Macedo, R., Malaquias, A., Ferreira, A., Brum, L., Portugal, I., 2010. Genetic analysis of extensively drug-resistant *Mycobacterium tuberculosis* strains in Lisbon, Portugal. J. Antimicrob. Chemother. 65, 224-227.

Perdigao, J., Macedo, R., Silva, C., Machado, D., Couto, I., Viveiros, M., Jordao, L., Portugal, I., 2013. From multidrug-resistant to extensively drug-resistant tuberculosis in Lisbon, Portugal: the stepwise mode of resistance acquisition. J. Antimicrob. Chemother. 68, 27-33.

Perdigao, J., Silva, H., Machado, D., Macedo, R., Maltez, F., Silva, C., Jordao, L., Couto, I., Mallard, K., Coll, F., Hill-Cawthorne, G.A., McNerney, R., Pain, A., Clark, T.G., Viveiros, M., Portugal, I., 2014. Unraveling *Mycobacterium tuberculosis* genomic diversity and evolution in Lisbon, Portugal, a highly drug resistant setting. BMC Genomics 15, 991.

Perdigao, J., Maltez, F., Machado, D., Silva, H., Pereira, C., Silva, C., Couto, I., Viveiros, M., Portugal, I., 2016. Beyond extensively drug-resistant tuberculosis in Lisbon, Portugal: a case of linezolid resistance acquisition presenting as an iliopsoas abscess. Int. J. Antimicrob. Agents 48, 569-570.

Perdigao, J., Clemente, S., Ramos, J., Masakidi, P., Machado, D., Silva, C., Couto, I., Viveiros, M., Taveira, N., Portugal, I., 2017. Genetic diversity, transmission dynamics and drug resistance of *Mycobacterium tuberculosis* in Angola. Sci. Rep. 7, 42814.

Prado Jr., C., 2012. História Econômica do Brasil. Editora Brasileira, São Paulo.

Rabna, P., Ramos, J., Ponce, G., Sanca, L., Mane, M., Armada, A., Machado, D., Vieira, F., Gomes, V.F., Martins, E., Colombatti, R., Riccardi, F., Perdigao, J., Sotero, J., Portugal, I., Couto, I., Atouguia, J., Rodrigues, A., Viveiros, M., 2015. Direct detection by the Xpert MTB/RIF assay and characterization of multi and poly drug-resistant tuberculosis in Guinea-Bissau, West Africa. PLoS One 10, e0127536.

Reyes, J.F., Chan, C.H., Tanaka, M.M., 2012. Impact of homopylas on variable numbers of tandem repeats and spoligotypes in *Mycobacterium tuberculosis*. Infect. Genet. Evol. 12, 811-818.

Saifodine, A., Fyfe, J., Sievers, A., Coelho, E., Azam, K., Black, J., 2014. Genetic diversity of *Mycobacterium tuberculosis* isolates obtained from patients with pulmonary tuberculosis in Beira city, Mozambique. Int. J. Mycobacteriol. 3, 94-100.

Shamputa, I.C., Lee, J., Allix-Beguec, C., Cho, E.J., Lee, J.I., Rajan, V., Lee, E.G., Min, J.H., Carroll, M.W., Goldfeder, L.C., Kim, J.H., Kang, H.S., Hwang, S., Eum, S.Y., Park, S.K., Lee, H., Supply, P., Cho, S.N., Via, L.E., Barry III, C.E., 2010. Genetic diversity of *Mycobacterium tuberculosis* isolates from a tertiary care tuberculosis hospital in South Korea. J. Clin. Microbiol. 48, 387-394.

Sheen, P., Couvin, D., Grandjean, L., Zimic, M., Dominguez, M., Luna, G., Gilman, R.H., Rastogi, N., Moore, D.A., 2013. Genetic diversity of *Mycobacterium tuberculosis* in Peru and exploration of phylogenetic associations with drug resistance. PLoS One 8, e65873.

Skiba, Y., Mokrousov, I., Ismagulova, G., Maltseva, E., Yurkevich, N., Bismilda, V.,

Chingissova, L., Abildaev, T., Aitkhozhina, N., 2015. Molecular snapshot of *Mycobacterium tuberculosis* population in Kazakhstan: a country-wide study. Tuberculosis (Edinburgh, Scotland) 95, 538-546.

Sola, C., Filliol, I., Legrand, E., Mokrousov, I., Rastogi, N., 2001. *Mycobacterium tuberculosis* phylogeny reconstruction based on combined numerical analysis with IS1081, IS6110, VNTR, and DR-based spoligotyping suggests the existence of two new phylogeographical clades. J. Mol. Evol. 53, 680-689.

Sola, C., Filliol, I., Legrand, E., Lesjean, S., Locht, C., Supply, P., Rastogi, N., 2003. Genotyping of the *Mycobacterium tuberculosis* complex using MIRUs: association with VNTR and spoligotyping for molecular epidemiology and evolutionary genetics. Infect. Genet. Evol. 3, 125-133.

Sun, Y.J., Bellamy, R., Lee, A.S., Ng, S.T., Ravindran, S., Wong, S.Y., Locht, C., Supply, P., Paton, N.I., 2004. Use of mycobacterial interspersed repetitive unit-variable-number tandem repeat typing to examine genetic diversity of *Mycobacterium tuberculosis* in Singapore. J. Clin. Microbiol. 42, 1986-1993.

Supply, P., 2005. Multilocus variable number tandem repeat. In: Genotyping of *Mycobacterium tuberculosis*. Technical Guide. Institut de Biologie/Institut Pasteur de Lille.

Supply, P., Warren, R.M., Banuls, A.L., Lesjean, S., Van Der Spuy, G.D., Lewis, L.A., Tibayrenc, M., Van Helden, P.D., Locht, C., 2003. Linkage disequilibrium between minisatellite loci supports clonal evolution of Mycobacterium tuberculosis in a high tuberculosis incidence area. Mol. Microbiol. 47, 529-538.

Supply, P., Allix, C., Lesjean, S., Cardoso-Oelemann, M., Rusch-Gerdes, S., Willery, E., Savine, E., de Haas, P., van Deutekom, H., Roring, S., Bifani, P., Kurepina, N., Kreiswirth, B., Sola, C., Rastogi, N., Vatin, V., Gutierrez, M.C., Fauville, M., Niemann, S., Skuce, R., Kremer, K., Locht, C., van Soolingen, D., 2006. Proposal for standardization of optimized mycobacterial interspersed repetitive unit-variable-number tandem repeat typing of *Mycobacterium tuberculosis*. J. Clin. Microbiol. 44, 4498-4510.

Tolentino, N., Rocha, C., Tolentino, C., Peixoto, J., 2011. Migração, remessas e desenvolvimento em África: o caso dos países de língua portuguesa. Revista Internacional em Língua Portuguesa 23-45.

van Soolingen, D., 2014. Whole-genome sequencing of *Mycobacterium tuberculosis* as an epidemiological marker. Lancet Respir. Med. 2, 251-252.

Viegas, S.O., Machado, A., Groenheit, R., Ghebremichael, S., Pennhag, A., Gudo, P.S., Cuna, Z., Miotto, P., Hill, V., Marrufo, T., Cirillo, D.M., Rastogi, N., Kallenius, G., Koivula, T., 2010. Molecular diversity of *Mycobacterium tuberculosis* isolates from patients with pulmonary tuberculosis in Mozambique. BMC Microbiol. 10, 195.

Viegas, S.O., Machado, A., Groenheit, R., Ghebremichael, S., Pennhag, A., Gudo, P.S., Cuna, Z., Langa, E., Miotto, P., Cirillo, D.M., Rastogi, N., Warren, R.M., van Helden, P.D., Koivula, T., Kallenius, G., 2013. *Mycobacterium tuberculosis* Beijing genotype is associated with HIV infection in Mozambique. PLoS One 8, e71999.

Wada, T., Iwamoto, T., Maeda, S., 2009. Genetic diversity of the *Mycobacterium tuberculosis* Beijing family in East Asia revealed through refined population structure analysis. FEMS Microbiol. Lett. 291, 35-43.

Walker, T.M., Ip, C.L., Harrell, R.H., Evans, J.T., Kapatai, G., Dedicoat, M.J., Eyre, D.W., Wilson, D.J., Hawkey, P.M., Crook, D.W., Parkhill, J., Harris, D., Walker, A.S., Bowden, R., Monk, P., Smith, E.G., Peto, T.E., 2013. Whole-genome sequencing to delineate *Mycobacterium tuberculosis* outbreaks: a retrospective observational study. Lancet Infect. Dis. 13, 137-146.

Weniger, T., Krawczyk, J., Supply, P., Niemann, S., Harmsen, D., 2010. MIRU-VNTRplus: a web tool for polyphasic genotyping of *Mycobacterium tuberculosis* complex bacteria. Nucleic Acids Res. 38, W326-331.

World Health Organization, 1994. TB: A Global Emergency. World Health Organization, Geneva.

World Health Organization, 2016. Global Tuberculosis Report 2016. World Health Organization, Geneva.

Xia, E., Teo, Y.Y., Ong, R.T., 2016. SpoTyping: fast and accurate in silico Mycobacterium spoligotyping from sequence reads. Genome Med. 8, 19.
<figure><img src="image_9.png" /><figcaption>Supplementary Figure S1 - Minimum Spanning Tree (MST) based on all available spoligotyping patterns in this study (n=1239 isolates). Node coloring is shown in function of country of isolation and thickness and/or shape of branches (continuous, dashed, dotted, black or gray) varied in function of spacer changes between patterns. The number of changes is indicated on branches.</figcaption></figure>
<figure><img src="image_10.png" /><figcaption>Supplementary Figure S2 – MIRU-VNTR dendrogram (15 loci) of 802 M. tuberculosis clinical isolates. Annotation: cluster identification, country of origin (ISO 3166), drug susceptibility data, MIRU-VNTR alleles, spoligotyping hybridization pattern, SIT and clade. First-line drug susceptibility testing: I, INH; R, RIF; S, STP; E, EMB; P, PZA. Second-line drug susceptibility testing: KAN, kanamycin; AMK, amikacin; CAP, capreomycin; OFX, ofloxacin; MOX, moxifloxacin; ETH, ethionamide; PAS, para-amino salicylic acid; CS, cycloserine; na, not available, nd, not determined.</figcaption></figure>
<figure><img src="image_11.png" /></figure>

Supplementary Figure S3 – MIRU-VNTR dendrogram (24 loci) of 381 M. tuberculosis clinical isolates. Annotation: with cluster identification, country of origin (ISO 3166), drug susceptibility data, MIRU-VNTR alleles, spoligotyping hybridization pattern, SIT and clade. First-line drug susceptibility testing: I, INH; R, RIF; S, STP; E, EMB; P, PZA. Second-line drug susceptibility testing: KAN, kanamycin; AMK, amikacin; CAP, capreomycin; OFX, ofloxacin; MOX, moxifloxacin; ETH, ethionamide; PAS, para-amino salicylic acid; CS, cycloserine; na, not available; nd, not determined.
<figure><img src="image_12.png" /><figcaption>Supplementary Figure S4 – MIRU-VNTR dendrogram (combined 24 and 15 loci) of 818 M. tuberculosis clinical isolates. Annotation: with cluster identification, country of origin (ISO 3166), drug susceptibility data, MIRU-VNTR alleles, spoligotyping hybridization pattern, SIT and clade. First-line drug susceptibility testing: I, INH; R, RIF; S, STP; E, EMB; P, PZA. Second-line drug susceptibility testing: KAN, kanamycin; AMK, amikacin; CAP, capreomycin; OFX, ofloxacin; MOX, moxifloxacin; ETH, ethionamide; PAS, para-amino salicylic acid; CS, cycloserine; na, not available; nd, not determined.</figcaption></figure>
<table><caption>Supplementary Table S1: Spoligotype International Types (SIT) recorded in this study.</caption><thead><tr><th>SIT</th><th>Spoligotype Description</th><th>Octal code</th><th>Nb in study</th><th>% in study</th><th>% in study vs. SITVIT</th><th>Lineage</th><th>Country Distribution**</th></tr></thead><tbody><tr><td>1</td><td>0000000000003771</td><td>33</td><td>2.66</td><td>0.31</td><td></td><td>Beijing</td><td></td></tr><tr><td>2</td><td>000000004020771</td><td>1</td><td>0.08</td><td>0.2</td><td></td><td>H2</td><td></td></tr><tr><td>3</td><td>000000007720771</td><td>4</td><td>0.32</td><td>3.64</td><td></td><td>H3</td><td></td></tr><tr><td>4</td><td>000000007760771</td><td>2</td><td>0.16</td><td>0.49</td><td></td><td>Unknown</td><td></td></tr><tr><td>17</td><td>677737607760771</td><td>42</td><td>3.39</td><td>5.62</td><td></td><td>LAM2</td><td></td></tr><tr><td>18</td><td>677776777760601</td><td>1</td><td>0.08</td><td>2.44</td><td></td><td>X2</td><td></td></tr><tr><td>20</td><td>677777607760771</td><td>157</td><td>12.67</td><td>16.86</td><td></td><td>LAM1</td><td></td></tr><tr><td>26</td><td>703777740003771</td><td>2</td><td>0.16</td><td>0.11</td><td></td><td>CAS1-Delhi</td><td></td></tr><tr><td>32</td><td>7760000000000171</td><td>1</td><td>0.08</td><td>0.97</td><td></td><td>Unknown</td><td></td></tr><tr><td>33</td><td>776177607760771</td><td>29</td><td>2.34</td><td>1.99</td><td></td><td>LAM3</td><td></td></tr><tr><td>34</td><td>776377777760771</td><td>5</td><td>0.4</td><td>0.56</td><td></td><td>S</td><td></td></tr><tr><td>36</td><td>777737777720771</td><td>1</td><td>0.08</td><td>0.64</td><td></td><td>H3</td><td></td></tr><tr><td>37</td><td>777737777760771</td><td>1</td><td>0.08</td><td>0.18</td><td></td><td>T3</td><td></td></tr><tr><td>42</td><td>777777607760771</td><td>94</td><td>7.59</td><td>2.51</td><td></td><td>LAM9</td><td></td></tr><tr><td>45</td><td>777777764020771</td><td>58</td><td>4.68</td><td>53.21</td><td></td><td>H1</td><td></td></tr><tr><td>47</td><td>77777774020771</td><td>17</td><td>1.37</td><td>1.0</td><td></td><td>H1</td><td></td></tr><tr><td>49</td><td>77777777720731</td><td>1</td><td>0.08</td><td>0.39</td><td></td><td>H3</td><td></td></tr><tr><td>50</td><td>777777777720771</td><td>59</td><td>4.76</td><td>1.41</td><td></td><td>H3</td><td></td></tr><tr><td>51</td><td>77777777760700</td><td>1</td><td>0.08</td><td>0.31</td><td></td><td>T1</td><td></td></tr><tr><td>52</td><td>77777777760731</td><td>3</td><td>0.24</td><td>0.31</td><td></td><td>T2</td><td></td></tr><tr><td>53</td><td>77777777760771</td><td>92</td><td>7.43</td><td>1.32</td><td></td><td>T1</td><td></td></tr><tr><td>54</td><td>77777777763771</td><td>5</td><td>0.4</td><td>1.79</td><td></td><td>Manu2</td><td></td></tr><tr><td>56</td><td>777777700000000</td><td>1</td><td>0.08</td><td>3.23</td><td></td><td>Unknown</td><td></td></tr><tr><td>58</td><td>777777557760771</td><td>12</td><td>0.97</td><td>5.61</td><td></td><td>T5-Madrid2</td><td></td></tr><tr><td>60</td><td>777777607760731</td><td>16</td><td>1.29</td><td>3.47</td><td></td><td>LAM4</td><td></td></tr><tr><td>61</td><td>777777743760771</td><td>1</td><td>0.08</td><td>0.1</td><td></td><td>Cameroon</td><td></td></tr><tr><td>64</td><td>777777607560771</td><td>15</td><td>1.21</td><td>3.26</td><td></td><td>LAM6</td><td></td></tr><tr><td>65</td><td>77777777760471</td><td>60</td><td>4.84</td><td>42.86</td><td></td><td>T1</td><td></td></tr><tr><td>73</td><td>777737777760731</td><td>13</td><td>1.05</td><td>4.11</td><td></td><td>T</td><td></td></tr><tr><td>74</td><td>77777760016771</td><td>1</td><td>0.08</td><td>6.25</td><td></td><td>Unknown</td><td></td></tr><tr><td>75</td><td>777767777720771</td><td>6</td><td>0.48</td><td>6.98</td><td></td><td>H3</td><td></td></tr><tr><td>81</td><td>777377607760771</td><td>4</td><td>0.32</td><td>10.53</td><td></td><td>LAM9</td><td></td></tr><tr><td>86</td><td>777777737760771</td><td>1</td><td>0.08</td><td>1.03</td><td></td><td>T1</td><td></td></tr><tr><td>90</td><td>67777776000371</td><td>1</td><td>0.08</td><td>2.94</td><td></td><td>Unknown</td><td></td></tr><tr><td>91</td><td>700036777760771</td><td>2</td><td>0.16</td><td>0.54</td><td></td><td>X3</td><td></td></tr><tr><td>92</td><td>700076777760771</td><td>4</td><td>0.32</td><td>0.87</td><td></td><td>X3</td><td></td></tr><tr><td>93</td><td>777737607760771</td><td>28</td><td>2.26</td><td>5.87</td><td></td><td>LAM5</td><td></td></tr><tr><td>95</td><td>777777607560731</td><td>1</td><td>0.08</td><td>1.16</td><td></td><td>LAM6</td><td></td></tr><tr><td>99</td><td>75777777720771</td><td>3</td><td>0.24</td><td>3.66</td><td></td><td>H3</td><td></td></tr><tr><td>105</td><td>776160000000071</td><td>1</td><td>0.08</td><td>1.89</td><td></td><td>Unknown</td><td></td></tr><tr><td>106</td><td>776177400000171</td><td>20</td><td>1.61</td><td>13.7</td><td></td><td>Unknown</td><td></td></tr><tr><td>118</td><td>777777777760771</td><td>4</td><td>0.32</td><td>2.31</td><td></td><td>T1</td><td></td></tr><tr><td>119</td><td>777776777760771</td><td>6</td><td>0.48</td><td>0.54</td><td></td><td>X1</td><td></td></tr><tr><td>129</td><td>700777747413771</td><td>2</td><td>0.16</td><td>3.08</td><td></td><td>EA16-BGD1</td><td></td></tr><tr><td>130</td><td>776177607760731</td><td>1</td><td>0.08</td><td>0.72</td><td></td><td>LAM3</td><td></td></tr><tr><td>137</td><td>777776777760601</td><td>12</td><td>0.97</td><td>1.12</td><td></td><td>X2</td><td></td></tr><tr><td>144</td><td>770000003760771</td><td>5</td><td>0.4</td><td>8.2</td><td></td><td>Unknown</td><td></td></tr><tr><td>150</td><td>777776707760771</td><td>6</td><td>0.48</td><td>5.66</td><td></td><td>LAM9</td><td></td></tr><tr><td>151</td><td>777767774020771</td><td>1</td><td>0.08</td><td>3.7</td><td></td><td>H1</td><td></td></tr><tr><td>153</td><td>75777777760731</td><td>1</td><td>0.08</td><td>1.27</td><td></td><td>T2</td><td></td></tr><tr><td>176</td><td>777737607560771</td><td>2</td><td>0.16</td><td>4.88</td><td></td><td>LAM6</td><td></td></tr><tr><td>177</td><td>377777607760771</td><td>6</td><td>0.48</td><td>4.65</td><td></td><td>LAM9</td><td></td></tr><tr><td>179</td><td>677737607760771</td><td>6</td><td>0.48</td><td>27.27</td><td></td><td>LAM</td><td></td></tr><tr><td>181</td><td>77077777777671</td><td>2</td><td>0.16</td><td>0.69</td><td></td><td>AFRI_1</td><td></td></tr><tr><td>187</td><td>7707777777767671</td><td>1</td><td>0.08</td><td>1.01</td><td></td><td>AFRI_1</td><td></td></tr><tr><td>194</td><td>677737607760731</td><td>1</td><td>0.08</td><td>3.03</td><td></td><td>LAM2</td><td></td></tr><tr><td>195</td><td>677777607660771</td><td>1</td><td>0.08</td><td>6.25</td><td></td><td>LAM1</td><td></td></tr><tr><td>211</td><td>776137607760771</td><td>3</td><td>0.24</td><td>2.52</td><td></td><td>LAM3</td><td></td></tr><tr><td>216</td><td>777717607760771</td><td>2</td><td>0.16</td><td>3.64</td><td></td><td>LAM5</td><td></td></tr><tr><td>217</td><td>777736777760771</td><td>1</td><td>0.08</td><td>2.27</td><td></td><td>X1</td><td></td></tr><tr><td>219</td><td>777747077760771</td><td>1</td><td>0.08</td><td>1.16</td><td></td><td>T1</td><td></td></tr><tr><td>237</td><td>777777777700000</td><td>1</td><td>0.08</td><td>0.76</td><td></td><td>Unknown</td><td></td></tr><tr><td>244</td><td>777777777760601</td><td>15</td><td>1.21</td><td>10.64</td><td></td><td>T1</td><td></td></tr><tr><td>245</td><td>777777777760671</td><td>1</td><td>0.08</td><td>4.35</td><td></td><td>T1</td><td></td></tr><tr><td>250</td><td>0000000000000371</td><td>1</td><td>0.08</td><td>2.56</td><td></td><td>Beijing</td><td></td></tr><tr><td>253</td><td>777777663760771</td><td>1</td><td>0.08</td><td>1.82</td><td></td><td>T1</td><td></td></tr><tr><td>254</td><td>77776007760771</td><td>1</td><td>0.08</td><td>0.53</td><td></td><td>T5-RUS1</td><td></td></tr><tr><td>262</td><td>774777777420771</td><td>3</td><td>0.24</td><td>1.4</td><td></td><td>Ural-1</td><td></td></tr><tr><td>266</td><td>775740003760771</td><td>1</td><td>0.08</td><td>6.67</td><td></td><td>T1</td><td></td></tr><tr><td>290</td><td>777777606760771</td><td>1</td><td>0.08</td><td>3.23</td><td></td><td>LAM8</td><td></td></tr><tr><td>306</td><td>777777601560771</td><td>1</td><td>0.08</td><td>3.33</td><td></td><td>T1</td><td></td></tr><tr><td>333</td><td>777727607760771</td><td>3</td><td>0.24</td><td>25.0</td><td></td><td>LAM5</td><td></td></tr><tr><td>347</td><td>777736777760601</td><td>1</td><td>0.08</td><td>1.92</td><td></td><td>X2</td><td></td></tr><tr><td>373</td><td>777777767760771</td><td>3</td><td>0.24</td><td>4.23</td><td></td><td>T1</td><td></td></tr><tr><td>382</td><td>75777774020771</td><td>1</td><td>0.08</td><td>7.69</td><td></td><td>H1</td><td></td></tr><tr><td>383</td><td>776377774020771</td><td>8</td><td>0.65</td><td>33.33</td><td></td><td>H1</td><td></td></tr></tbody></table>
<table>
  <tbody>
    <tr>
      <td>388</td>
      <td>737777607760771</td>
      <td>2</td>
      <td>0.16</td>
      <td>6.45</td>
      <td>LAM9</td>
      <td></td>
    </tr>
    <tr>
      <td>389</td>
      <td>67777607760771</td>
      <td>3</td>
      <td>0.24</td>
      <td>15.0</td>
      <td>LAM1</td>
      <td></td>
    </tr>
    <tr>
      <td>396</td>
      <td>777777607560371</td>
      <td>3</td>
      <td>0.24</td>
      <td>15.79</td>
      <td>LAM6</td>
      <td></td>
    </tr>
    <tr>
      <td>397</td>
      <td>777777600000771</td>
      <td>1</td>
      <td>0.08</td>
      <td>3.85</td>
      <td>LAM</td>
      <td></td>
    </tr>
    <tr>
      <td>398</td>
      <td>777777607760631</td>
      <td>1</td>
      <td>0.08</td>
      <td>5.0</td>
      <td>LAM4</td>
      <td></td>
    </tr>
    <tr>
      <td>402</td>
      <td>777777600000000</td>
      <td>1</td>
      <td>0.08</td>
      <td>2.22</td>
      <td>Unknown</td>
      <td></td>
    </tr>
    <tr>
      <td>451</td>
      <td>777777477760771</td>
      <td>1</td>
      <td>0.08</td>
      <td>0.68</td>
      <td>T-H37Rv</td>
      <td></td>
    </tr>
    <tr>
      <td>453</td>
      <td>77777777560771</td>
      <td>1</td>
      <td>0.08</td>
      <td>4.76</td>
      <td>T1</td>
      <td></td>
    </tr>
    <tr>
      <td>481</td>
      <td>676773677777600</td>
      <td>1</td>
      <td>0.08</td>
      <td>0.02</td>
      <td>BOV_1</td>
      <td></td>
    </tr>
    <tr>
      <td>482</td>
      <td>67677377777600</td>
      <td>1</td>
      <td>0.08</td>
      <td>0.06</td>
      <td>BOV_1</td>
      <td></td>
    </tr>
    <tr>
      <td>503</td>
      <td>777737377760731</td>
      <td>1</td>
      <td>0.08</td>
      <td>25.0</td>
      <td>T</td>
      <td></td>
    </tr>
    <tr>
      <td>578</td>
      <td>637777607760771</td>
      <td>1</td>
      <td>0.08</td>
      <td>2.78</td>
      <td>LAM1</td>
      <td></td>
    </tr>
    <tr>
      <td>602</td>
      <td>777777770000771</td>
      <td>5</td>
      <td>0.4</td>
      <td>4.81</td>
      <td>Unknown</td>
      <td></td>
    </tr>
    <tr>
      <td>612</td>
      <td>777777777760751</td>
      <td>1</td>
      <td>0.08</td>
      <td>2.94</td>
      <td>T1</td>
      <td></td>
    </tr>
    <tr>
      <td>635</td>
      <td>000000007560771</td>
      <td>2</td>
      <td>0.16</td>
      <td>16.67</td>
      <td>LAM3</td>
      <td></td>
    </tr>
    <tr>
      <td>729</td>
      <td>657777607760771</td>
      <td>3</td>
      <td>0.24</td>
      <td>23.08</td>
      <td>LAM1</td>
      <td></td>
    </tr>
    <tr>
      <td>740</td>
      <td>777777747720771</td>
      <td>1</td>
      <td>0.08</td>
      <td>3.85</td>
      <td>H3</td>
      <td></td>
    </tr>
    <tr>
      <td>806</td>
      <td>75777777413731</td>
      <td>2</td>
      <td>0.16</td>
      <td>3.7</td>
      <td>EAI1-SOM</td>
      <td></td>
    </tr>
    <tr>
      <td>811</td>
      <td>777777604060731</td>
      <td>2</td>
      <td>0.16</td>
      <td>3.51</td>
      <td>LAM11-ZWE</td>
      <td></td>
    </tr>
    <tr>
      <td>815</td>
      <td>777777606060731</td>
      <td>2</td>
      <td>0.16</td>
      <td>1.34</td>
      <td>LAM11-ZWE</td>
      <td></td>
    </tr>
    <tr>
      <td>826</td>
      <td>677737207760771</td>
      <td>1</td>
      <td>0.08</td>
      <td>7.14</td>
      <td>LAM2</td>
      <td></td>
    </tr>
    <tr>
      <td>828</td>
      <td>377777607760731</td>
      <td>1</td>
      <td>0.08</td>
      <td>2.86</td>
      <td>LAM4</td>
      <td></td>
    </tr>
    <tr>
      <td>831</td>
      <td>776367777760771</td>
      <td>1</td>
      <td>0.08</td>
      <td>7.14</td>
      <td>S</td>
      <td></td>
    </tr>
    <tr>
      <td>863</td>
      <td>000000007400000</td>
      <td>28</td>
      <td>2.26</td>
      <td>560.0</td>
      <td>PINI2</td>
      <td></td>
    </tr>
    <tr>
      <td>866</td>
      <td>577777607760771</td>
      <td>1</td>
      <td>0.08</td>
      <td>1.69</td>
      <td>LAM9</td>
      <td></td>
    </tr>
    <tr>
      <td>885</td>
      <td>777773607460771</td>
      <td>1</td>
      <td>0.08</td>
      <td>50.0</td>
      <td>LAM5</td>
      <td></td>
    </tr>
    <tr>
      <td>891</td>
      <td>777777607660771</td>
      <td>1</td>
      <td>0.08</td>
      <td>3.33</td>
      <td>LAM9</td>
      <td></td>
    </tr>
    <tr>
      <td>946</td>
      <td>77777740020771</td>
      <td>1</td>
      <td>0.08</td>
      <td>4.76</td>
      <td>H</td>
      <td></td>
    </tr>
    <tr>
      <td>1074</td>
      <td>777777607760761</td>
      <td>1</td>
      <td>0.08</td>
      <td>3.13</td>
      <td>LAM9</td>
      <td></td>
    </tr>
    <tr>
      <td>1106</td>
      <td>7777777607760400</td>
      <td>38</td>
      <td>3.07</td>
      <td>172.73</td>
      <td>LAM4</td>
      <td></td>
    </tr>
    <tr>
      <td>1154</td>
      <td>777777607760751</td>
      <td>1</td>
      <td>0.08</td>
      <td>5.26</td>
      <td>LAM9</td>
      <td></td>
    </tr>
    <tr>
      <td>1176</td>
      <td>777757607760771</td>
      <td>6</td>
      <td>0.48</td>
      <td>30.0</td>
      <td>LAM9</td>
      <td></td>
    </tr>
    <tr>
      <td>1224</td>
      <td>776177607760700</td>
      <td>1</td>
      <td>0.08</td>
      <td>12.5</td>
      <td>LAM3</td>
      <td></td>
    </tr>
    <tr>
      <td>1241</td>
      <td>777777607700771</td>
      <td>2</td>
      <td>0.16</td>
      <td>2.2</td>
      <td>LAM9</td>
      <td></td>
    </tr>
    <tr>
      <td>1247</td>
      <td>777777607763771</td>
      <td>1</td>
      <td>0.08</td>
      <td>6.67</td>
      <td>Manu2</td>
      <td></td>
    </tr>
    <tr>
      <td>1273</td>
      <td>700076777760731</td>
      <td>1</td>
      <td>0.08</td>
      <td>12.5</td>
      <td>X3</td>
      <td></td>
    </tr>
    <tr>
      <td>1321</td>
      <td>677777607760731</td>
      <td>2</td>
      <td>0.16</td>
      <td>6.9</td>
      <td>LAM</td>
      <td></td>
    </tr>
    <tr>
      <td>1367</td>
      <td>377737607760771</td>
      <td>1</td>
      <td>0.08</td>
      <td>16.67</td>
      <td>LAM5</td>
      <td></td>
    </tr>
    <tr>
      <td>1435</td>
      <td>77477777413731</td>
      <td>1</td>
      <td>0.08</td>
      <td>25.0</td>
      <td>EAI1-SOM</td>
      <td></td>
    </tr>
    <tr>
      <td>1528</td>
      <td>774377607760771</td>
      <td>2</td>
      <td>0.16</td>
      <td>66.67</td>
      <td>LAM9</td>
      <td></td>
    </tr>
    <tr>
      <td>1530</td>
      <td>777777607760711</td>
      <td>1</td>
      <td>0.08</td>
      <td>8.33</td>
      <td>LAM4</td>
      <td></td>
    </tr>
    <tr>
      <td>1535</td>
      <td>777577607760771</td>
      <td>1</td>
      <td>0.08</td>
      <td>7.14</td>
      <td>LAM9</td>
      <td></td>
    </tr>
    <tr>
      <td>1536</td>
      <td>776377607760771</td>
      <td>1</td>
      <td>0.08</td>
      <td>4.76</td>
      <td>LAM9</td>
      <td></td>
    </tr>
    <tr>
      <td>1545</td>
      <td>601775607760771</td>
      <td>1</td>
      <td>0.08</td>
      <td>16.67</td>
      <td>LAM1</td>
      <td></td>
    </tr>
    <tr>
      <td>1548</td>
      <td>777703400000360</td>
      <td>2</td>
      <td>0.16</td>
      <td>25.0</td>
      <td>Unknown</td>
      <td></td>
    </tr>
    <tr>
      <td>1610</td>
      <td>763777607560771</td>
      <td>2</td>
      <td>0.16</td>
      <td>33.33</td>
      <td>LAM6</td>
      <td></td>
    </tr>
    <tr>
      <td>1659</td>
      <td>77770370000000</td>
      <td>5</td>
      <td>0.4</td>
      <td>100.0</td>
      <td>Unknown</td>
      <td></td>
    </tr>
    <tr>
      <td>1662</td>
      <td>600000000060731</td>
      <td>3</td>
      <td>0.24</td>
      <td>150.0</td>
      <td>T2</td>
      <td></td>
    </tr>
    <tr>
      <td>1691</td>
      <td>677727607760771</td>
      <td>2</td>
      <td>0.16</td>
      <td>10.53</td>
      <td>LAM2</td>
      <td></td>
    </tr>
    <tr>
      <td>1750</td>
      <td>777767607760731</td>
      <td>3</td>
      <td>0.24</td>
      <td>33.33</td>
      <td>LAM4</td>
      <td></td>
    </tr>
    <tr>
      <td>1752</td>
      <td>677777607740771</td>
      <td>4</td>
      <td>0.32</td>
      <td>133.33</td>
      <td>LAM1</td>
      <td></td>
    </tr>
    <tr>
      <td>1755</td>
      <td>677777607560771</td>
      <td>1</td>
      <td>0.08</td>
      <td>10.0</td>
      <td>LAM6</td>
      <td></td>
    </tr>
    <tr>
      <td>1757</td>
      <td>757777600000000</td>
      <td>1</td>
      <td>0.08</td>
      <td>33.33</td>
      <td>Unknown</td>
      <td></td>
    </tr>
    <tr>
      <td>1758</td>
      <td>777777607540771</td>
      <td>4</td>
      <td>0.32</td>
      <td>33.33</td>
      <td>LAM9</td>
      <td></td>
    </tr>
    <tr>
      <td>1759</td>
      <td>772177606760771</td>
      <td>2</td>
      <td>0.16</td>
      <td>66.67</td>
      <td>LAM3</td>
      <td></td>
    </tr>
    <tr>
      <td>1779***</td>
      <td>776477777760771</td>
      <td>1</td>
      <td>0.08</td>
      <td>NA</td>
      <td>Unknown</td>
      <td></td>
    </tr>
    <tr>
      <td>1803</td>
      <td>617777607760771</td>
      <td>1</td>
      <td>0.08</td>
      <td>7.14</td>
      <td>LAM1</td>
      <td></td>
    </tr>
    <tr>
      <td>1872</td>
      <td>77777774060771</td>
      <td>2</td>
      <td>0.16</td>
      <td>40.0</td>
      <td>T1</td>
      <td></td>
    </tr>
    <tr>
      <td>1894</td>
      <td>747777607760771</td>
      <td>2</td>
      <td>0.16</td>
      <td>66.67</td>
      <td>LAM9</td>
      <td></td>
    </tr>
    <tr>
      <td>1895</td>
      <td>777777607460731</td>
      <td>1</td>
      <td>0.08</td>
      <td>6.67</td>
      <td>LAM4</td>
      <td></td>
    </tr>
    <tr>
      <td>1908</td>
      <td>777777776020771</td>
      <td>1</td>
      <td>0.08</td>
      <td>16.67</td>
      <td>H3</td>
      <td></td>
    </tr>
    <tr>
      <td>1949</td>
      <td>703777700000771</td>
      <td>1</td>
      <td>0.08</td>
      <td>5.88</td>
      <td>CAS</td>
      <td></td>
    </tr>
    <tr>
      <td>2025</td>
      <td>777777777700371</td>
      <td>2</td>
      <td>0.16</td>
      <td>16.67</td>
      <td>T</td>
      <td></td>
    </tr>
    <tr>
      <td>2073</td>
      <td>776000607760771</td>
      <td>1</td>
      <td>0.08</td>
      <td>25.0</td>
      <td>LAM3</td>
      <td></td>
    </tr>
    <tr>
      <td>2258</td>
      <td>777777777700011</td>
      <td>1</td>
      <td>0.08</td>
      <td>8.33</td>
      <td>Unknown</td>
      <td></td>
    </tr>
    <tr>
      <td>2271</td>
      <td>677717607760771</td>
      <td>3</td>
      <td>0.24</td>
      <td>30.0</td>
      <td>LAM2</td>
      <td></td>
    </tr>
    <tr>
      <td>2375</td>
      <td>76777774020771</td>
      <td>3</td>
      <td>0.24</td>
      <td>21.43</td>
      <td>H1</td>
      <td></td>
    </tr>
    <tr>
      <td>2506</td>
      <td>677757607760771</td>
      <td>2</td>
      <td>0.16</td>
      <td>40.0</td>
      <td>LAM1</td>
      <td></td>
    </tr>
    <tr>
      <td>2512</td>
      <td>657757607760771</td>
      <td>19</td>
      <td>1.53</td>
      <td>316.67</td>
      <td>LAM1</td>
      <td></td>
    </tr>
    <tr>
      <td>2513</td>
      <td>777777777760471</td>
      <td>5</td>
      <td>0.4</td>
      <td>125.0</td>
      <td>T1</td>
      <td></td>
    </tr>
    <tr>
      <td>2514</td>
      <td>774377607760700</td>
      <td>2</td>
      <td>0.16</td>
      <td>28.57</td>
      <td>LAM4</td>
      <td></td>
    </tr>
    <tr>
      <td>2535</td>
      <td>777777207600331</td>
      <td>1</td>
      <td>0.08</td>
      <td>5.88</td>
      <td>Unknown</td>
      <td></td>
    </tr>
    <tr>
      <td>2571</td>
      <td>776340000000171</td>
      <td>2</td>
      <td>0.16</td>
      <td>66.67</td>
      <td>Unknown</td>
      <td></td>
    </tr>
    <tr>
      <td>2572</td>
      <td>667777607760771</td>
      <td>2</td>
      <td>0.16</td>
      <td>50.0</td>
      <td>LAM1</td>
      <td></td>
    </tr>
    <tr>
      <td>2580</td>
      <td>77774001636771</td>
      <td>3</td>
      <td>0.24</td>
      <td>150.0</td>
      <td>T1</td>
      <td></td>
    </tr>
    <tr>
      <td>2640</td>
      <td>677777777720731</td>
      <td>1</td>
      <td>0.08</td>
      <td>50.0</td>
      <td>H3</td>
      <td></td>
    </tr>
    <tr>
      <td>2669</td>
      <td>0000000000000000</td>
      <td>5</td>
      <td>0.4</td>
      <td>4.55</td>
      <td>ATYPIC</td>
      <td></td>
    </tr>
    <tr>
      <td>2815</td>
      <td>770777404177671</td>
      <td>1</td>
      <td>0.08</td>
      <td>50.0</td>
      <td>AFRI_1</td>
      <td></td>
    </tr>
    <tr>
      <td>3040</td>
      <td>773767607760671</td>
      <td>1</td>
      <td>0.08</td>
      <td>50.0</td>
      <td>LAM9</td>
      <td>CO=1, ES=1</td>
    </tr>
    <tr>
      <td>3104</td>
      <td>760377777743671</td>
      <td>1</td>
      <td>0.08</td>
      <td>50.0</td>
      <td>AFRI_1</td>
      <td>GW=2</td>
    </tr>
    <tr>
      <td>3111</td>
      <td>7037777700000000</td>
      <td>1</td>
      <td>0.08</td>
      <td>50.0</td>
      <td>CAS_1-Delhi</td>
      <td>GW=2</td>
    </tr>
  </tbody>
</table>
<table><tr><td>3450</td><td>77777405720771</td><td>1</td><td>0.08</td><td>33.33</td><td>H3</td><td>FI=1, ES=1, FR=1</td></tr><tr><td>3457</td><td>67577760776771</td><td>1</td><td>0.08</td><td>50.0</td><td>LAM1</td><td>FR=1, VE=1</td></tr><tr><td>3585</td><td>476763677777600</td><td>1</td><td>0.08</td><td>12.5</td><td>BOV</td><td>ES=8</td></tr><tr><td>4004</td><td>777777600000331</td><td>3</td><td>0.24</td><td>100.0</td><td>Unknown</td><td>BR=3</td></tr><tr><td>4073</td><td>777777777760661</td><td>1</td><td>0.08</td><td>50.0</td><td>T1</td><td>ZM=1, BR=1</td></tr><tr><td>4078</td><td>777777677560731</td><td>4</td><td>0.32</td><td>200.0</td><td>T2</td><td>BR=2</td></tr><tr><td>4092</td><td>777777717413731</td><td>1</td><td>0.08</td><td>50.0</td><td>EA1-SOM</td><td>MZ=2</td></tr><tr><td>4139*</td><td>777737207400000</td><td>2</td><td>0.16</td><td>100.0</td><td>Unknown</td><td>this study n=2</td></tr><tr><td>4140*</td><td>700777747000071</td><td>2</td><td>0.16</td><td>100.0</td><td>Unknown</td><td>this study n=2</td></tr><tr><td>4141*</td><td>77777777766771</td><td>1</td><td>0.08</td><td>50.0</td><td>Unknown</td><td>this study n=1, NG=1</td></tr><tr><td>4142*</td><td>777777777777671</td><td>1</td><td>0.08</td><td>50.0</td><td>Unknown</td><td>this study n=1, US=1</td></tr><tr><td>4143*</td><td>771777600000000</td><td>1</td><td>0.08</td><td>50.0</td><td>Unknown</td><td>this study n=1, BR=1</td></tr><tr><td>4144*</td><td>557767764000771</td><td>1</td><td>0.08</td><td>50.0</td><td>Unknown</td><td>this study n=1, BR=1</td></tr><tr><td>4145*</td><td>777767764020771</td><td>1</td><td>0.08</td><td>50.0</td><td>H1</td><td>this study n=1, DE=1</td></tr><tr><td>4146*</td><td>577767477760771</td><td>1</td><td>0.08</td><td>50.0</td><td>T-H37Rv</td><td>this study n=1, GB=1</td></tr><tr><td>4147*</td><td>677747607760771</td><td>1</td><td>0.08</td><td>50.0</td><td>LAM1</td><td>this study n=1, AR=1</td></tr><tr><td>4148*</td><td>457747607760771</td><td>2</td><td>0.16</td><td>100.0</td><td>LAM1</td><td>this study n=2</td></tr><tr><td>4149*</td><td>77777740160771</td><td>1</td><td>0.08</td><td>50.0</td><td>T1</td><td>this study n=1, AT=1</td></tr><tr><td>4150*</td><td>776337700760771</td><td>1</td><td>0.08</td><td>50.0</td><td>S</td><td>this study n=1, BR=1</td></tr><tr><td>4151*</td><td>776167400000171</td><td>3</td><td>0.24</td><td>75.0</td><td>Unknown</td><td>this study n=3, ES=1</td></tr><tr><td>4152*</td><td>777037607760771</td><td>1</td><td>0.08</td><td>50.0</td><td>LAM5</td><td>this study n=1, BR=1</td></tr><tr><td>4153*</td><td>747777607560771</td><td>1</td><td>0.08</td><td>50.0</td><td>LAM6</td><td>this study n=1, BR=1</td></tr><tr><td>4154*</td><td>776777607700771</td><td>1</td><td>0.08</td><td>50.0</td><td>Unknown</td><td>this study n=1, BR=1</td></tr><tr><td>4155*</td><td>757757607760771</td><td>2</td><td>0.16</td><td>100.0</td><td>LAM9</td><td>this study n=2</td></tr><tr><td>4156*</td><td>763767607560771</td><td>2</td><td>0.16</td><td>100.0</td><td>LAM6</td><td>this study n=2</td></tr><tr><td>4157*</td><td>776377760760771</td><td>6</td><td>0.48</td><td>100.0</td><td>S</td><td>this study n=6</td></tr><tr><td>4158*</td><td>777777557560771</td><td>1</td><td>0.08</td><td>50.0</td><td>T1</td><td>this study n=1, AR=1</td></tr><tr><td>4159*</td><td>777777670160731</td><td>2</td><td>0.16</td><td>100.0</td><td>T1</td><td>this study n=2</td></tr><tr><td>4160*</td><td>763777677560771</td><td>1</td><td>0.08</td><td>50.0</td><td>T1</td><td>this study n=1, BR=1</td></tr></table>

Newly created SITs are followed by an asterisk (*)

** Country distribution was provided for recently created SITs (the 2 letter country codes are according to http://en.wikipedia.org/wiki/ISO_3166-1 alpha-2)

*** SIT1779 in SITVITWEB was deleted in SITVIT2 database; the pattern is however retained pending its reassignment if new strains are reported in future.
<table><caption>Supplementary Table S2: Description of clusters containing 12 or more isolates in this study, and their worldwide distribution in the SITVIT database.</caption><thead><tr><th>SIT</th><th>Spoligotype Description</th><th>Octal code</th><th>Nb in study</th><th>% in study</th><th>% in study vs. database</th><th>Lineage</th><th>Distribution in Countries with ≥2% of a given SIT*</th><th>Distribution in Regions with ≥2% of a given SIT**</th></tr></thead><tbody><tr><td>1</td><td>0000000000003771</td><td>33</td><td>2.66</td><td>0.31</td><td>Beijing</td><td>CN=19.12, US=18.89, JP=10.86, ZA=7.81, RU=6.6, VN=3.67, IN=3.56, PE=2.96, MY=2.79</td><td>Eastern Asia=32.19, North America=18.93, Southeastern Asia=10.06, Southern Africa=7.81, Northern Asia=6.6, Southern Asia=5.92, Northern Europe=3.45, South America=3.42, Western Asia=2.81, Western Europe=2.42</td></tr><tr><td>17</td><td>677737607760771</td><td>42</td><td>3.39</td><td>5.62</td><td>LAM2</td><td>BR=30.52, VE=23.96, US=15.93, HT=7.36, ES=3.88, GP=2.68, CO=2.41</td><td>South America=59.97, North America=15.93, Caribbean=11.91, Southern Europe=4.95, Southern Africa=2.01</td></tr><tr><td>20</td><td>677777607760771</td><td>157</td><td>12.67</td><td>16.86</td><td>LAM1</td><td>BR=17.62, US=17.08, HT=7.3, NA=6.66, PT=5.26, FR=5.26, VE=4.51, ZA=2.9, ES=2.69, AR=2.36, PE=2.26, CO=2.15, BE=2.04</td><td>South America=31.36, North America=17.08, Western Europe=9.88, Southern Africa=9.56, Caribbean=9.13, Southern Europe=8.49, East Africa=3.97, Northern Europe=2.79, West Africa=2.69</td></tr><tr><td>33</td><td>776177607760771</td><td>29</td><td>2.34</td><td>1.99</td><td>LAM3</td><td>ZA=22.26, PE=14.38, BR=11.1, US=10.89, AR=10.62, ES=6.1, CL=3.77, FR=3.7, IT=3.15, HN=2.95</td><td>South America=41.58, Southern Africa=22.26, North America=10.89, Southern Europe=9.73, Western Europe=5.62, Central Africa=3.7, North Africa=2.19</td></tr><tr><td>42</td><td>777777607760771</td><td>94</td><td>7.59</td><td>2.51</td><td>LAM9</td><td>BR=14.45, US=10.6, CO=7.64, MA=6.3, IT=5.85, FR=4.54, AR=3.53, PE=3.31, ES=2.99, VE=2.96, ZA=2.83, HT=2.16, RU=2.06</td><td>South America=35.02, North America=10.6, Southern Europe=10.12, Western Europe=8.47, North Africa=7.67, Northern Europe=4.38, Caribbean=3.79, East Africa=3.47, Central America=3.18, Southern Africa=2.83, Northern Africa=2.06</td></tr><tr><td>45</td><td>777777764020771</td><td>58</td><td>4.68</td><td>53.21</td><td>H1</td><td>AT=16.51, MQ=11.01, US=9.17, CO=8.26, IT=8.26, BR=8.26, GP=7.34, DE=3.67, NL=2.75, LC=2.75, SE=2.75</td><td>Western Europe=23.85, Caribbean=22.02, South America=19.27, Southern Europe=11.01, North America=9.17, Northern Europe=6.42</td></tr><tr><td>47</td><td>77777774020771</td><td>17</td><td>1.37</td><td>1.0</td><td>H1</td><td>US=13.96, IT=7.57, AT=7.39, BR=7.16, PE=6.92, FR=6.1, FI=5.51, CZ=3.46, ES=3.28, SE=3.11, AR=2.64, MA=2.58, NL=2.29, CF=2.17, TR=2.05</td><td>South America=18.71, Western Europe=18.65, North America=14.25, Southern Europe=12.38, Northern Europe=9.62, Eastern Europe=6.63, Western Asia=3.7, North Africa=3.34, Caribbean=2.52, Central Africa=2.4</td></tr><tr><td>50</td><td>77777777720771</td><td>59</td><td>4.76</td><td>1.41</td><td>H3</td><td>US=14.2, PE=12.6, BR=7.85, FR=5.56, AT=4.92, ES=4.39, IT=4.39, ZA=3.27, AR=3.08, CM=3.01, CZ=2.96, HT=2.77, TR=2.58, SE=2.29</td><td>South America=28.77, Western Europe=14.2, North America=14.2, Southern Europe=9.33, Caribbean=4.7, Eastern Europe=4.46, Northern Europe=4.42, Western Asia=3.63, North Africa=3.44, Southern Africa=3.27, Central Africa=3.06, Eastern Asia=2.03</td></tr><tr><td>53</td><td>777777777760771</td><td>92</td><td>7.43</td><td>1.32</td><td>T1</td><td>US=11.96, FR=7.14, BR=5.87, IT=4.82, TR=4.51, ZA=4.39, PE=4.29, AT=3.09, CN=2.82, FI=2.65, MX=2.58, ET=2.35, SA=2.26, AR=2.19, ES=2.09</td><td>South America=17.09, Western Europe=14.15, North America=12.19, Southern Europe=8.51, Western Asia=7.97, Northern Europe=6.77, East Africa=4.65, Southern Africa=4.49, Eastern Asia=3.86, North Africa=3.18, Eastern Europe=2.95, Central America=2.92, Caribbean=2.88, Southern Africa=2.86</td></tr><tr><td>58</td><td>77777557760771</td><td>12</td><td>0.97</td><td>5.61</td><td>T5-Madrid2</td><td>US=21.96, ES=19.63, PE=11.68, AR=11.68, BR=8.88, FR=6.07, MA=2.8, CL=2.34, CU=2.34, NL=1.87, VE=1.87, TN=1.4, BE=0.93, PT=0.93, IT=0.93, ZA=0.93, AT=0.93, PY=0.47, DZ=0.47, MQ=0.47, GF=0.47, SE=0.47, MX=0.47</td><td>South America=37.38, North America=21.96, Southern Europe=21.5, Western Europe=9.81, North Africa=4.67, Caribbean=2.8</td></tr><tr><td>60</td><td>777777607760771</td><td>16</td><td>1.29</td><td>3.47</td><td>LAM4</td><td>ZA=35.36, BR=15.84, FR=4.77, US=4.12, MA=4.12, GM=3.69, IT=3.47, VE=3.25, PE=3.04</td><td>Southern Africa=35.36, South America=26.9, Western Europe=6.72, North Africa=6.51, West Africa=5.86, Southern Europe=5.21, North America=4.12, East Africa=2.39</td></tr><tr><td>64</td><td>777777607560771</td><td>15</td><td>1.21</td><td>3.26</td><td>LAM6</td><td>BR=44.35, US=20.43, GF=5.0, PT=3.04, PE=2.61, FR=2.61, AR=2.39</td><td>South America=58.26, North America=20.43, Western Europe=5.22, Southern Europe=3.91, Western Asia=2.61</td></tr><tr><td>65</td><td>77777777760471</td><td>60</td><td>4.84</td><td>42.86</td><td>T1</td><td>BR=39.29, HT=21.43, US=9.29, FR=6.43, ES=3.57, AT=2.86, CM=2.14</td><td>South America=40.71, Caribbean=23.57, Western Europe=12.14, North America=9.29, Southern Europe=4.29, Central Africa=2.14, Northern Europe=2.14</td></tr><tr><td>73</td><td>777737777760771</td><td>13</td><td>1.05</td><td>4.11</td><td>T</td><td>BR=18.35, IT=11.71, US=11.39, ZA=8.86, FR=6.33, AR=4.43, CN=3.8, MZ=3.16, MX=2.53, GF=2.22</td><td>South America=26.9, Southern Europe=13.92, Western Europe=11.39, North America=11.39, Southern Africa=8.86, East Africa=6.33, Eastern Asia=3.8, Central America=3.16, West Africa=2.22, Western Asia=2.22</td></tr><tr><td>93</td><td>777737607760771</td><td>28</td><td>2.26</td><td>5.87</td><td>LAM5</td><td>VE=19.71, US=15.72, BR=14.47, PE=13.42, HT=10.06, IT=4.4, GP=3.98, ES=2.31, FR=2.31</td><td>South America=53.88, Caribbean=15.72, North America=15.72, Southern Europe=6.71, Western Europe=3.77</td></tr><tr><td>106</td><td>776177400000171</td><td>20</td><td>1.61</td><td>13.7</td><td>Unknown</td><td>ES=32.88, US=15.75, PA=13.01, BR=9.59, CO=4.79, AR=4.79, PE=3.42, FR=3.42, MA=2.74, DZ=2.05</td><td>Southern Europe=34.25, South America=23.97, North America=15.75, Central America=13.01, North Africa=4.79, Western Europe=4.79</td></tr><tr><td>137</td><td>777776777760601</td><td>12</td><td>0.97</td><td>1.12</td><td>X2</td><td>US=68.28, BR=7.46, ZA=6.06, GB=5.13</td><td>North America=69.31, South America=8.21, Northern Europe=6.81, Southern Africa=6.06, Western Europe=2.61</td></tr><tr><td>244</td><td>777777777760601</td><td>15</td><td>1.21</td><td>10.64</td><td>T1</td><td>BR=23.4, ZA=13.48, PT=11.35, BD=9.93, FR=9.22, ZM=5.67, GW=4.96, TZ=4.26, US=3.55, MZ=2.13</td><td>South America=24.82, East Africa=14.18, Southern Africa=13.48, Southern Europe=12.77, Southern Asia=11.35, Western Europe=10.64, West Africa=5.67, North America=3.55</td></tr><tr><td>863</td><td>000000007400000</td><td>28</td><td>2.26</td><td>560.0</td><td>PINI2</td><td>BR=80.0, US=20.0</td><td>South America=80.0, North America=20.0</td></tr><tr><td>1106</td><td>777777607760400</td><td>38</td><td>3.07</td><td>172.73</td><td>LAM4</td><td>PT=40.91, BR=31.82, ES=18.18, IT=4.55, PE=4.55</td><td>Southern Europe=63.64, South America=36.36</td></tr><tr><td>2512</td><td>657757607760771</td><td>19</td><td>1.53</td><td>316.67</td><td>LAM1</td><td>BR=83.33, CL=16.67</td><td>South America=100.0</td></tr></tbody></table>

<footnote>* The 2 letter country codes are according to http://en.wikipedia.org/wiki/ISO_3166-1 alpha-2; countrywide distribution is only shown for SITs with ≥2% of a given SIT as compared to their total number in the SITVIT database.</footnote>

<footnote>** Worldwide distribution is reported for regions with more than 2% of a given SIT as compared to their total number in the SITVIT database.</footnote>
<table><caption>**Supplementary Table S3** – Fifteen-loci MIRU-VNTR clusters found, main associated SITs and distribution by country of origin.</caption><thead><tr><th rowspan="2">MIRU-VNTR<br>Cluster</th><th rowspan="2">Associated<br>SITs</th><th colspan="6">No. of Isolates</th></tr><tr><th>Angola</th><th>Brazil</th><th>Guinea-<br>Bissau</th><th>Mozambique</th><th>Portugal</th><th>Total</th></tr></thead><tbody><tr><td>AO-01</td><td>20</td><td>2</td><td>0</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>AO-02</td><td>1548</td><td>2</td><td>0</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>AO-03</td><td>244</td><td>2</td><td>0</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>AO-04</td><td>20</td><td>2</td><td>0</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>CPLP-01</td><td>20, 1755, 42</td><td>8</td><td>6</td><td>0</td><td>0</td><td>2</td><td>16</td></tr><tr><td>CPLP-02</td><td>20, 2271, 2572</td><td>1</td><td>0</td><td>0</td><td>0</td><td>2</td><td>3</td></tr><tr><td>CPLP-03</td><td>237</td><td>0</td><td>1</td><td>0</td><td>0</td><td>1</td><td>2</td></tr><tr><td>CPLP-04</td><td>20</td><td>1</td><td>1</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>CPLP-05</td><td>1530</td><td>1</td><td>3</td><td>0</td><td>0</td><td>0</td><td>4</td></tr><tr><td>CPLP-06</td><td>1530, 4152</td><td>0</td><td>1</td><td>0</td><td>0</td><td>1</td><td>2</td></tr><tr><td>CPLP-07</td><td>194, 17</td><td>1</td><td>2</td><td>0</td><td>0</td><td>1</td><td>4</td></tr><tr><td>CPLP-08</td><td>33, 1224</td><td>0</td><td>1</td><td>0</td><td>0</td><td>1</td><td>2</td></tr><tr><td>CPLP-09</td><td>nd</td><td>1</td><td>1</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>CPLP-10</td><td>53</td><td>1</td><td>0</td><td>0</td><td>0</td><td>1</td><td>2</td></tr><tr><td>GW-01</td><td>1</td><td>0</td><td>0</td><td>2</td><td>0</td><td>0</td><td>2</td></tr><tr><td>GW-02</td><td>1</td><td>0</td><td>0</td><td>3</td><td>0</td><td>0</td><td>3</td></tr><tr><td>Lisboa3-A</td><td>20, 42</td><td>0</td><td>0</td><td>0</td><td>0</td><td>11</td><td>11</td></tr><tr><td>Lisboa3-B</td><td>20, 42</td><td>0</td><td>0</td><td>0</td><td>0</td><td>38</td><td>38</td></tr><tr><td>POA-01</td><td>863</td><td>0</td><td>4</td><td>0</td><td>0</td><td>0</td><td>4</td></tr><tr><td>POA-02</td><td>17</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>PT-01</td><td>92</td><td>0</td><td>0</td><td>0</td><td>0</td><td>2</td><td>2</td></tr><tr><td>PT-02</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>3</td><td>3</td></tr><tr><td>PT-03</td><td>1752</td><td>0</td><td>0</td><td>0</td><td>0</td><td>4</td><td>4</td></tr><tr><td>PT-04</td><td>211</td><td>0</td><td>0</td><td>0</td><td>0</td><td>2</td><td>2</td></tr><tr><td>PT-05</td><td>53</td><td>0</td><td>0</td><td>0</td><td>0</td><td>2</td><td>2</td></tr><tr><td>PT-06</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>2</td><td>2</td></tr><tr><td>PT-07</td><td>1106</td><td>0</td><td>0</td><td>0</td><td>0</td><td>4</td><td>4</td></tr><tr><td>PT-08</td><td>20, 2271</td><td>0</td><td>0</td><td>0</td><td>0</td><td>3</td><td>3</td></tr><tr><td>PT-09</td><td>20, 42</td><td>0</td><td>0</td><td>0</td><td>0</td><td>2</td><td>2</td></tr><tr><td>Q1</td><td>290, 42, 1106</td><td>1</td><td>0</td><td>0</td><td>0</td><td>26</td><td>27</td></tr><tr><td>RG-01</td><td>2512, 729</td><td>0</td><td>20</td><td>0</td><td>0</td><td>0</td><td>20</td></tr><tr><td>RG-02</td><td>17</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>RG-03</td><td>65, 45, 2512,<br>4145</td><td>0</td><td>27</td><td>0</td><td>0</td><td>0</td><td>27</td></tr><tr><td>RG-04</td><td>nd</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>RG-05</td><td>45</td><td>0</td><td>13</td><td>0</td><td>0</td><td>0</td><td>13</td></tr><tr><td>RG-06</td><td>nd</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>RG-07</td><td>20</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>RG-08</td><td>45</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>RG-09</td><td>20</td><td>0</td><td>8</td><td>0</td><td>0</td><td>0</td><td>8</td></tr><tr><td>RG-10</td><td>65</td><td>0</td><td>4</td><td>0</td><td>0</td><td>0</td><td>4</td></tr><tr><td>RG-11</td><td>65</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2</td></tr></tbody></table>
<table><tr><td>RG-12</td><td>nd</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>RG-13</td><td>nd</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>RG-14</td><td>602</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>RG-15</td><td>45</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>RG-16</td><td>45</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>RG-17</td><td>137</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>RG-18</td><td>53</td><td>0</td><td>10</td><td>0</td><td>0</td><td>0</td><td>10</td></tr><tr><td>RG-19</td><td>602</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>RG-20</td><td>nd</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>RG-21</td><td>45</td><td>0</td><td>3</td><td>0</td><td>0</td><td>0</td><td>3</td></tr><tr><td>RG-22</td><td>33</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>RG-23</td><td>53, 50</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>RG-24</td><td>60, 177</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>RG-25</td><td>20, 65</td><td>0</td><td>5</td><td>0</td><td>0</td><td>0</td><td>5</td></tr><tr><td>RG-26</td><td>177</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>RG-27</td><td>33</td><td>0</td><td>4</td><td>0</td><td>0</td><td>0</td><td>4</td></tr><tr><td>RG-28</td><td>60</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>RG-29</td><td>65</td><td>0</td><td>4</td><td>0</td><td>0</td><td>0</td><td>4</td></tr><tr><td>RG-30</td><td>65</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>RG-31</td><td>42, 4157</td><td>0</td><td>4</td><td>0</td><td>0</td><td>0</td><td>4</td></tr><tr><td>RG-32</td><td>nd</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>RG-33</td><td>nd</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>RG-34</td><td>4148</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>RG-35</td><td>nd</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>RS-01</td><td>58</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>RS-02</td><td>50</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>RS-03</td><td>65</td><td>0</td><td>24</td><td>0</td><td>0</td><td>0</td><td>24</td></tr><tr><td>RS-04</td><td>106</td><td>0</td><td>4</td><td>0</td><td>0</td><td>0</td><td>4</td></tr><tr><td>RS-05</td><td>863, 4157</td><td>0</td><td>4</td><td>0</td><td>0</td><td>0</td><td>4</td></tr><tr><td>RS-06</td><td>99</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>RS-07</td><td>53</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>RS-08</td><td>53</td><td>0</td><td>3</td><td>0</td><td>0</td><td>0</td><td>3</td></tr><tr><td>RS-09</td><td>179</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>RS-10</td><td>453</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2</td></tr></table>

nd - not determined
<table><caption>**Supplementary Table S4** – Twenty-four-loci MIRU-VNTR clusters found, main associated SITs and distribution by country of origin.</caption><thead><tr><th rowspan="2">MIRU-VNTR Cluster</th><th rowspan="2">Associated SITs</th><th colspan="6">No. of Isolates</th></tr><tr><th>Angola</th><th>Brazil</th><th>Guinea-Bissau</th><th>Mozambique</th><th>Portugal</th><th>Total</th></tr></thead><tbody><tr><td>AO-01</td><td>20</td><td>2</td><td>0</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>AO-02</td><td>1548</td><td>2</td><td>0</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>AO-03</td><td>244</td><td>2</td><td>0</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>AO-04</td><td>20</td><td>2</td><td>0</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>CPLP-01</td><td>20, 42</td><td>6</td><td>1</td><td>0</td><td>0</td><td>1</td><td>8</td></tr><tr><td>CPLP-01B</td><td>20</td><td>1</td><td>0</td><td>0</td><td>0</td><td>1</td><td>2</td></tr><tr><td>CPLP-02</td><td>20, 2271, 2572</td><td>1</td><td>0</td><td>0</td><td>0</td><td>2</td><td>3</td></tr><tr><td>CPLP-03</td><td>237</td><td>0</td><td>1</td><td>0</td><td>0</td><td>1</td><td>2</td></tr><tr><td>CPLP-07</td><td>17, 194</td><td>1</td><td>0</td><td>0</td><td>0</td><td>1</td><td>2</td></tr><tr><td>CPLP_11</td><td>42, 290</td><td>1</td><td>0</td><td>0</td><td>0</td><td>1</td><td>2</td></tr><tr><td>GW-01</td><td>1</td><td>0</td><td>0</td><td>2</td><td>0</td><td>0</td><td>2</td></tr><tr><td>GW-02</td><td>1</td><td>0</td><td>0</td><td>3</td><td>0</td><td>0</td><td>3</td></tr><tr><td>Lisboa3-A</td><td>20, 42</td><td>0</td><td>0</td><td>0</td><td>0</td><td>9</td><td>9</td></tr><tr><td>Lisboa3-B</td><td>20, 42</td><td>0</td><td>0</td><td>0</td><td>0</td><td>37</td><td>37</td></tr><tr><td>POA-01</td><td>863</td><td>0</td><td>3</td><td>0</td><td>0</td><td>0</td><td>3</td></tr><tr><td>POA-02</td><td>17</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>PT-01</td><td>92</td><td>0</td><td>0</td><td>0</td><td>0</td><td>2</td><td>2</td></tr><tr><td>PT-02</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>3</td><td>3</td></tr><tr><td>PT-03</td><td>1752</td><td>0</td><td>0</td><td>0</td><td>0</td><td>4</td><td>4</td></tr><tr><td>PT-04</td><td>211</td><td>0</td><td>0</td><td>0</td><td>0</td><td>2</td><td>2</td></tr><tr><td>PT-05</td><td>53</td><td>0</td><td>0</td><td>0</td><td>0</td><td>2</td><td>2</td></tr><tr><td>PT-06</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>2</td><td>2</td></tr><tr><td>PT-07</td><td>1106</td><td>0</td><td>0</td><td>0</td><td>0</td><td>4</td><td>4</td></tr><tr><td>Q1</td><td>1106</td><td>0</td><td>0</td><td>0</td><td>0</td><td>22</td><td>22</td></tr><tr><td>RS-01</td><td>58</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>RS-02</td><td>50</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>RS-03</td><td>65</td><td>0</td><td>8</td><td>0</td><td>0</td><td>0</td><td>8</td></tr><tr><td>RS-04</td><td>106</td><td>0</td><td>4</td><td>0</td><td>0</td><td>0</td><td>4</td></tr><tr><td>RS-05</td><td>863</td><td>0</td><td>3</td><td>0</td><td>0</td><td>0</td><td>3</td></tr><tr><td>RS-06</td><td>99</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>RS-08</td><td>53</td><td>0</td><td>3</td><td>0</td><td>0</td><td>0</td><td>3</td></tr><tr><td>RS-09</td><td>179</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>RS-10</td><td>453</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2</td></tr></tbody></table>