# A first insight on the population structure of Mycobacterium tuberculosis complex as studied by spoligotyping and MIRU-VNTRs in Santiago, Chile.

María Elvira Balcells, Patricia García, Paulina Meza, Carlos Peña, Marcela Cifuentes, David Couvin, Nalin Rastogi

► To cite this version:

María Elvira Balcells, Patricia García, Paulina Meza, Carlos Peña, Marcela Cifuentes, et al.. A first insight on the population structure of Mycobacterium tuberculosis complex as studied by spoligotyping and MIRU-VNTRs in Santiago, Chile.. PLoS ONE, 2015, 10 (2), pp.e0118007. 10.1371/journal.pone.0118007.pone-01153988

HAL Id: pasteur-01153988

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

Submitted on 20 May 2015

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 ARTICLE

# A First Insight on the Population Structure of *Mycobacterium tuberculosis* Complex as Studied by Spoligotyping and MIRU-VNTRs in Santiago, Chile

María Elvira Balcells¹*, Patricia García², Paulina Meza², Carlos Peña³, Marcela Cifuentes³, David Couvin⁴, Nalin Rastogi⁴

1 Infectious Diseases Department, School of Medicine, Pontificia Universidad Católica de Chile, Santiago, Chile, 2 Microbiology Laboratory, Clinical Laboratory Department, School of Medicine, Pontificia Universidad Católica de Chile, Santiago, Chile, 3 Respiratory Division and Microbiology Laboratory, Hospital San Borja Arriarán, Santiago, Chile, 4 WHO Supranational TB Reference Laboratory, Institut Pasteur de la Guadeloupe, Abymes, Guadeloupe, France

* balcells@med.puc.cl

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

CrossMark
click for updates

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

OPEN ACCESS

**Citation:** Balcells ME, García P, Meza P, Peña C, Cifuentes M, Couvin D, et al. (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(2): e0118007. doi:10.1371/journal.pone.0118007

**Academic Editor:** Riccardo Manganelli, University of Padova, Medical School, ITALY

**Received:** October 30, 2014

**Accepted:** January 5, 2015

**Published:** February 11, 2015

**Copyright:** © 2015 Balcells et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files.

**Funding:** This work was supported by funding by Chilean National Funding Agency CONICYT, FONDECYT GRANT N°1130600 (http://www.conicyt.cl/fondecyt/). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** The authors have declared that no competing interests exist.

## Abstract

Tuberculosis (TB) remains a significant public health problem worldwide, but the ecology of the prevalent mycobacterial strains, and their transmission, can vary depending on country and region. Chile is a country with low incidence of TB, that has a geographically isolated location in relation to the rest of South American countries due to the Andes Mountains, but recent migration from neighboring countries has changed this situation. We aimed to assess the genotypic diversity of *Mycobacterium tuberculosis* complex (MTBC) strains in Santiago, Chile, and compare with reports from other Latin-American countries. We analyzed MTBC isolates from pulmonary tuberculosis cases collected between years 2008 and 2013 in Central Santiago, using two genotyping methods: spoligotyping and 12-loci mycobacterial interspersed repetitive unit-variable number of tandem repeats (MIRU-VNTRs). Data obtained were analyzed and compared to the SITVIT2 database. Mean age of the patients was 47.5 years and 61% were male; 11.6% were migrants. Of 103 strains (1 isolate/patient) included, there were 56 distinct spoligotype patterns. Of these, 16 strains (15.5%) corresponded to orphan strains in the SITVIT2 database, not previously reported. Latin American and Mediterranean (LAM) (34%) and T (33%) lineages were the most prevalent strains, followed by Haarlem lineage (16.5%). Beijing family was scarcely represented with only two cases (1.9%), one of them isolated from a Peruvian migrant. The most frequent clustered spoligotypes were SIT33/LAM3 (10.7%), SIT53/T1 (8.7%), SIT50/H3 (7.8%), and SIT37/T3 (6.8%). We conclude that LAM and T genotypes are the most prevalent genotypes of MTBC in Santiago, Chile, and together correspond to almost two thirds of analyzed strains, which is similar to strain distribution reported from other countries of Latin America. Nevertheless, the high proportion of SIT37/T3, which was rarely found in other Latin
American countries, may underline a specific history or demographics of Chile related to
probable human migrations and evolutions.

**Introduction**

Tuberculosis (TB) is an infectious disease that remains present broadly despite intense global efforts for its control and elimination. In 2012, there were 8.6 million incident cases and more than 1.3 million deaths from tuberculosis worldwide [1]. Its causative agent, *Mycobacterium tuberculosis*, has accompanied human beings for thousands of years and, despite this very long period of coexistence, *M. tuberculosis* has remained among the most extreme examples of bacterial genetic homogeneity, with divergence rarely higher than 0.5 single nucleotide polymorphisms per genome per year [2][3]. Accordingly, for decades scientists have believed that different clinical or epidemiological outcomes of this infection in humans were largely related to host susceptibility and environmental factors, and not to the bacteria. However, this view has been challenged from the mid 2000's, when comparative genomic and molecular epidemiological tools identified large sequence polymorphisms (LSPs) that allowed to differentiate *M. tuberculosis* into main six main lineages that defined their origin [4]. These lineages were defined as Indo-Oceanic (lineage I), East Asian (lineage II, including Beijing), East African Indian (lineage III), Euro-American (lineage IV), West-Africa (lineage V, M. *africanum* I) and West Africa (lineage VI; M. *africanum* II). Lineages I, V and VI were considered “ancient” and II to IV “modern” in relation to the presence or absence of the TbD1 genomic region, deleted in modern strains [4].

Independently in an earlier study, a distinction between three genetic groups of *M. tubercu-
losis* was achieved based on two non-synonymous polymorphisms occurring at high
frequency in the genes encoding catalase-peroxidase and the A subunit of gyrase, which led to
a classification in three principal genetic groups (PGG); group 1 bacteria being ancestral to
groups 2 and 3 [5]. Additional genotyping methods such as the presence or absence of spacer
oligonucleotides (spoligotyping), variable number of repeats (MIRU-VNTRs), or single nucle-
otide polymorphisms (SNPs), have not only provided further resolution at the sublineage level
but also corroborated split of *M. tuberculosis* into similar main lineages; interested readers may
refer to a recent review for comparison of spoligotyping-based nomenclature of *M. tuberculosis*
lineages vs. PGG groupings, SNPs and SNP based cluster groups (SCGs), and LSP-based line-
ages [6].

*M. tuberculosis* molecular profiling allows today not only to study the global biodiversity and the phylogeographical variations of tubercle bacilli [6], but also to determine if recurrent disease is caused by endogenous reactivation or exogenous reinfection, to look for associations of particular *M. tuberculosis* strains with their clinical behavior, to determine whether outbreaks in specific geographical areas are caused by the same (cluster) or by multiple strains, if some strains are more transmissible than others, and whether some strains are more prone to develop drug resistance [7][8][9]. It undeniably constitutes an important tool for tuberculosis control and epidemiological studies today.

Chile has reached very low TB rates after decades of sustained decline related to economic
development and a robust and well-coordinated National Tuberculosis Program. However, in
the last 5 years TB incidence rate has become stationary around 13/100,000; with over 2000
new cases diagnosed per year. Newer policies will have to be addressed and implemented in the
country to reverse this situation, such as strategies to detect earlier infections in high-risk
groups, introducing contact tracing and targeted tuberculosis preventive therapy. In this context, the introduction of molecular epidemiology tools is a valuable asset to help understand local TB transmission dynamics. Two previous reports published locally have described genotyping of a small number of MTBC strains collected at the Xth region in year 2006 and Metropolitan Region in 2013, and the most frequently identified spoligotypes were T1 (24%) and LAM (39.5%), respectively [10, 11]. However, none of those reports included phylogenetic analysis or epidemiological data to distinguish local versus imported strains.

We hereby report on a preliminary study to describe genetic diversity and main genotypic lineages of *M. tuberculosis* complex clinical isolates circulating in Santiago, Chile, and compare our findings with MTB genetic epidemiology studies from other South American countries.

# Materials and Methods

## Patients and bacterial isolates

This study describes genotyping results of a total of 103 *M. tuberculosis* isolates collected from August 2008 to September 2013, the majority of strains (84.5%) between 2011 and 2013. The specimens were obtained from the same number of pulmonary TB patients, in 22 outpatient clinics and 3 large hospitals in the central area of Santiago city, with a reference TB incidence of 16.6 per 100,000 population for the year 2012. Our samples represent around 6.44% of all smear-positive TB cases for that period in the whole Metropolitan Region. Basic demographic data (age, gender and country of origin) was obtained from mandatory reports included in the National Tuberculosis Program registry and anonymized. Strains included newly and previously treated cases, and none of the samples were duplicated from a single patient. Culture processing and strain isolating were performed at the central clinical laboratories from the Hospital San Borja Arriarán and Red de Salud de la Pontificia Universidad Católica de Chile. Positive cultures grown over the Lowenstein—Jensen medium or in automatized liquid (MGIT) culture were isolated and DNA extracted by standard methods. Testing for susceptibility to first line drugs was determined over Lowenstein—Jensen medium using the proportion method [12] at the national reference laboratory (ISP). According to programmatic guidelines susceptibility testing was done in patients considered to be at risk for drug resistance, namely patients with TB treatment failure, default or relapse, exposed to a drug-resistant TB case, HIV-associated TB and cases in migrants from countries with high prevalence of drug resistant TB.

## Nucleic acid extraction and genotyping

*M. tuberculosis* complex species were confirmed for all strains by a specific PCR for gyrase B (gyrB) with primers described by Kasai et al. [13]. Briefly each amplification reaction contained: 12,5µl 2X Promega premix, 2µl of each primer 10µM (forward and reverse), 0,2 µl of GoTaq flexi 5U/µl (Promega) DNA polymerase, 5 µl of nuclease-free distilled water and 5 µl of DNA from samples or control, in a final volume of 25µl. Thermal cycling conditions were as follows: 95°C, 2 minutes followed by 10 cycles of 95°C, 30 seconds; annealing temperature step-downs every cycle of 1°C (from 65°C to 55°C); 72°C, 30 seconds. The annealing temperature for the final 35 cycles was 55°C with denaturation and extension phases as above. PCR products were analyzed through electrophoresis in 1.5% agarose gels that were stained with ethidium bromide.

Mycobacterial DNA was prepared by heat inactivation and lysis by sonication, and stored at-20°C until use. Spoligotyping analysis was performed by using a commercial source for membranes and reagents (Ocimum Biosolutions BV) as described by Kamerbeek et al. [14] [15]. A total of 43 spacers between the direct repeats in the target region were amplified by
using DRa biotinylated at 5' end and DRb primers. The PCR products were hybridized to a membrane containing 43 oligonucleotides by reverse line blotting. A positive control with *M. tuberculosis* H37Rv strain and *M. bovis* BCG were used in each run.

Standard 12-loci MIRU-VNTRs typing was manually performed by conventional PCR amplification of 12 of 41 loci: 02, 04, 10, 16, 20, 23, 24, 26, 27, 31, 39 and 40, as previously described [16]. Each locus was amplified separately by simplex PCR with GoTaq flexi DNA polymerase (Promega). Products were analyzed by electrophoresis using 1.5–2% agarose (Promega, Fitchburg, USA) gels. From the gel images, the corresponding MIRU-VNTR bands were interpreted as copy numbers based on a reference table [17].

## Computer assisted genotype analysis and comparison with databases

Spoligotype patterns as octal codes were entered in the SITVIT2 proprietary database of the Institut Pasteur de la Guadeloupe which is an updated version of the previously released SITVITWEB database (available online at http://www.pasteur-guadeloupe.fr:8081/SITVIT_ONLINE/) [18]. At the time of the analysis on August 2014, the database contained genotypes of about 112,000 *M. tuberculosis* clinical isolates from about 170 countries of origin. In this database, Spoligotype International Type (SIT) and MIRU International Type (MIT) designates patterns shared by two or more patient isolates, whereas “orphan” designates patterns reported for a single isolate [18] [19]. Genotypic lineages were assigned according to the rules described in SITVITWEB [18] in which the LSP-based Euro-American lineage (lineage IV) is split in Latin American and Mediterranean (LAM), ill-defined T, Haarlem (H), X, and S lineages. Importantly, the nomenclature “East-African Indian” (EAI or lineage III) denotes two altogether different *M. tuberculosis* clades by spoligotyping vs. LSPs [6]; EAI by spoligotyping was named “Indo-Oceanic” by LSPs, while EAI by LSPs corresponds to Central-Asian (CAS) in the SITVIT2 database [6]. Note that lineages could be subdivided into sublineages.

Data were also analyzed using MIRU-VNTRplus web server (http://www.miru-vntrplus.org/) which is a collection of 186 strains representing the major *M. tuberculosis* complex lineages [19][20]. This database contains information of strain species, lineage, country together with copy numbers of 24 MIRU loci, spoligotyping patterns, regions of difference (RD) profiles, single nucleotide polymorphisms (SNPs), drug susceptibility data, and IS6110-RFLP fingerprint images.

## Phylogenetic analysis

BioNumerics software version 6.6 (Applied Maths, Sint-Martens-Latem, Belgium) was used to compare spoligotypes and MIRU-VNTRs patterns, by drawing Minimum spanning trees (MSTs) in order to visualize evolutionary relationships between the clinical isolates in our study. MSTs are undirected graphs in which all samples are connected together with the fewest possible connections between nearest neighbors. An Unweighted Pair Group Method with Arithmetic mean (UPGMA) tree was also constructed using MLVA Compare software (http://www.ridom.de/mlvacompare/) to visualize deeper associations between spoligotypes and MIRU-VNTRs involved in this study.

## Ethics Statement

Prior to the start of the study, ethical approval was obtained from Ethical Committee of Pontificia Universidad Católica de Chile School of Medicine (CEC-MedUC). The microbiological records and basic demographic information of patients were anonymized and de-identified prior to analysis. Therefore, informed consent was not required and was specifically waived by the Ethical Committee.
## Results

In total, 103 *M. tuberculosis* isolates were collected and analyzed; patients providing samples were 61.2% males and 38.8% females, with a mean age of 45.3 (range 3–93) and 51.05 years old (range 22–93), respectively. Although all patients were current residents in Santiago, Chile, 11.6% were foreign-born migrants with main country of origin being Peru (66.6%) and Bolivia (25%).

Spoligotyping of the 103 isolates revealed a total of 56 different patterns: 16 patterns (15.5%) corresponded to orphan strains in the SITVIT2 database, not yet reported, and 40 patterns (containing 87 strains) corresponded to shared-types or SITs in the SITVIT2 database. Among these 40 shared-types, 2 isolates from the present study in Chile matched an orphan strain recorded in the database from USA, therefore a new shared-type, SIT4013, was created. In addition, a total of ten 12-loci MIRU patterns (from MIT 1670 to MIT 1679) were newly created. Description of the 40 SITs (n = 87 isolates) and corresponding spoligotyping defined lineages/sublineages is shown in Table 1. Among orphan isolates, we noted similar spoligotyping patterns (with absence of spacer 6) among orphans Or08, Or11, and Or12 within our study. (S1 Table)

Distribution of main lineages for the 103 isolates according to SITVIT analysis, showed that 35 (34%) belonged to the LAM family; 34 (33%) belonged to the ill-defined T super-family; 17 (16.5%) belonged to the Haarlem (H) lineage; 3 (2.9%) belonged to the X lineage; 2 belonged to the S lineage and 2 (1.9%) belonged to the Beijing lineage. The rest belonged to “Unknown” or single strain lineages, and included genotypes such as Bovis, CAS, and AFRI (Table 2). Among all LAM strains (n = 35), the main sublineages were LAM3 (42.9%), LAM9 (40%) and LAM6 (8.6%). Among the ill-defined T super-family (n = 34) its prototype T1 (58.8%) was the 1st most frequent pattern in our study, followed by T3 (20.6%). The H lineage represented 16.5% (n = 17) of isolates subdivided in H3 (88.2%), and H1 (11.8%) sublineages. Diversity of strains among migrants (n = 12) were much higher and included lineages such as Beijing and CAS, which were uncommon among non-migrants. Comparison of lineage allocation by SITVIT and MIRU-VNTRplus was highly concordant for the majority of strains with the exception of T lineages: 27 out of 34 T strains on SITVIT analysis were reclassified by MIRU-VNTR, mostly as H (40.7%) and LAM (22.2%) (data shown in S2 Table).

Spoligotype analysis identified clustering (of 3 or more isolates) among 48 of 103 isolates, whereas with MIRU12 identified clustering (of 3 or more isolates) among only 14 of 103 strains. As shown in Table 3, the most frequently clustered spoligotypes were SIT33/LAM3 (n = 11, 10.7%), SIT53/T1 (n = 9, 8.8%), SIT50/H3 (n = 8, 7.8%) and SIT37/T3 (n = 7, 6.8%). The high presence of SIT37/T3 (6.8% of strains) was noteworthy as this SIT is usually more commonly found in Eastern Africa, and migration from Africa is highly rare in Chile. With respect to 12-loci MIRU clustering, the most frequently clustered isolate was 12-MIT585 (n = 6), that in SITVIT2 database corresponds to spoligotypes SIT37/T3 (n = 5) and SIT281/T1 (n = 1). Furthermore, in the SITVIT2 database, 12-MIT585 was present in 6 isolates only: in Spain from two South American patients (from Ecuador and Peru, respectively); and in USA (n = 1), Caribbean (n = 2), and Belgium (n = 1).

Minimum spanning trees (MSTs) illustrating evolutionary relationships between spoligotypes and 12-MIRU-VNTRs were constructed for all strains included in the present study (n = 103). Fig. 1A shows a MST based on spoligotypes in which three major groups belonging to the modern group PGG2/3 (LAM, T, H, X, S) were evident as highly predominant (representing 88.3% of all strains), most notably LAM (n = 35, 34%), T (n = 34, 33%) and Haarlem (n = 17, 16.5%). More distance was evident among isolates belonging to the T sublineages than among those integrating the LAM sublineages (Fig. 1A). In contrast, when MST was
<table><caption>Table 1. Description of 40 shared-types (SITs; n = 87 isolates) and corresponding spoligotyping defined lineages/sublineages starting from a total of 103 *M. tuberculosis* strains isolated in Santiago, Chile.</caption><thead><tr><th>SIT*</th><th>Spoligotype Description</th><th>Octal Number</th><th>Nb. In study (%)</th><th>% in study vs. database</th><th>Lineage**</th><th>Unique vs. Clustered SIT***</th></tr></thead><tbody><tr><td>1</td><td>000000000003771</td><td>2 (1.94)</td><td>0.02</td><td></td><td>Beijing</td><td>Clustered</td></tr><tr><td>20</td><td>67777607760771</td><td>1 (0.97)</td><td>0.11</td><td></td><td>LAM1</td><td>Unique</td></tr><tr><td>33</td><td>776177607760771</td><td>11 (10.68)</td><td>0.84</td><td></td><td>LAM3</td><td>Clustered</td></tr><tr><td>34</td><td>77637777760771</td><td>2 (1.94)</td><td>0.23</td><td></td><td>S</td><td>Clustered</td></tr><tr><td>37</td><td>7773777760771</td><td>7 (6.8)</td><td>1.34</td><td></td><td>T3</td><td>Clustered</td></tr><tr><td>39</td><td>77777347760471</td><td>2 (1.94)</td><td>1.35</td><td></td><td>T4-CEU1</td><td>Clustered</td></tr><tr><td>42</td><td>77777607760771</td><td>7 (6.8)</td><td>0.2</td><td></td><td>LAM9</td><td>Clustered</td></tr><tr><td>49</td><td>77777777720731</td><td>1 (0.97)</td><td>0.39</td><td></td><td>H3</td><td>Unique</td></tr><tr><td>50</td><td>77777777720771</td><td>8 (7.77)</td><td>0.2</td><td></td><td>H3</td><td>Clustered</td></tr><tr><td>52</td><td>77777777760731</td><td>3 (2.91)</td><td>0.31</td><td></td><td>T2</td><td>Clustered</td></tr><tr><td>53</td><td>77777777760771</td><td>9 (8.74)</td><td>0.14</td><td></td><td>T1</td><td>Clustered</td></tr><tr><td>58</td><td>77777557760771</td><td>1 (0.97)</td><td>0.52</td><td></td><td>T5-Madrid2</td><td>Unique</td></tr><tr><td>64</td><td>77777607560771</td><td>1 (0.97)</td><td>0.25</td><td></td><td>LAM6</td><td>Unique</td></tr><tr><td>91</td><td>70003677760771</td><td>1 (0.97)</td><td>0.28</td><td></td><td>X3</td><td>Unique</td></tr><tr><td>161</td><td>77777607740771</td><td>1 (0.97)</td><td>7.69</td><td></td><td>LAM9</td><td>Unique</td></tr><tr><td>180</td><td>67777777720771</td><td>1 (0.97)</td><td>1.67</td><td></td><td>H3</td><td>Unique</td></tr><tr><td>203</td><td>70377740001771</td><td>1 (0.97)</td><td>4.76</td><td></td><td>CAS</td><td>Unique</td></tr><tr><td>211</td><td>776137607760771</td><td>3 (2.91)</td><td>2.83</td><td></td><td>LAM3</td><td>Clustered</td></tr><tr><td>219</td><td>77774077760771</td><td>1 (0.97)</td><td>1.16</td><td></td><td>T1</td><td>Unique</td></tr><tr><td>222</td><td>777774077560771</td><td>2 (1.94)</td><td>2.06</td><td></td><td>Unknown</td><td>Clustered</td></tr><tr><td>281</td><td>77777577760771</td><td>1 (0.97)</td><td>3.23</td><td></td><td>T1</td><td>Unique</td></tr><tr><td>336</td><td>77777677760731</td><td>1 (0.97)</td><td>0.85</td><td></td><td>X1</td><td>Unique</td></tr><tr><td>373</td><td>77777767760771</td><td>1 (0.97)</td><td>1.47</td><td></td><td>T1</td><td>Unique</td></tr><tr><td>390</td><td>77777777620771</td><td>1 (0.97)</td><td>3.03</td><td></td><td>H3</td><td>Unique</td></tr><tr><td>452</td><td>77777606360771</td><td>1 (0.97)</td><td>5.88</td><td></td><td>LAM9</td><td>Unique</td></tr><tr><td>738</td><td>777777607760760</td><td>1 (0.97)</td><td>16.67</td><td></td><td>LAM9</td><td>Unique</td></tr><tr><td>746</td><td>77777777520771</td><td>1 (0.97)</td><td>2.86</td><td></td><td>H3</td><td>Unique</td></tr><tr><td>948</td><td>77777760020611</td><td>1 (0.97)</td><td>8.33</td><td></td><td>H3</td><td>Unique</td></tr><tr><td>1277</td><td>777777207760771</td><td>2 (1.94)</td><td>6.67</td><td></td><td>LAM9</td><td>Clustered</td></tr><tr><td>1293</td><td>676177607760771</td><td>1 (0.97)</td><td>12.5</td><td></td><td>LAM3</td><td>Unique</td></tr><tr><td>1355</td><td>777777407560731</td><td>1 (0.97)</td><td>0.54</td><td></td><td>LAM</td><td>Unique</td></tr><tr><td>1474</td><td>677777627760771</td><td>1 (0.97)</td><td>33.33</td><td></td><td>T1</td><td>Unique</td></tr><tr><td>1624</td><td>776177607560771</td><td>1 (0.97)</td><td>16.67</td><td></td><td>LAM</td><td>Unique</td></tr><tr><td>1877</td><td>77377777760771</td><td>1 (0.97)</td><td>7.14</td><td></td><td>T1</td><td>Unique</td></tr><tr><td>1914</td><td>377777607560771</td><td>1 (0.97)</td><td>16.67</td><td></td><td>LAM6</td><td>Unique</td></tr><tr><td>2030</td><td>777777704020751</td><td>1 (0.97)</td><td>33.33</td><td></td><td>H1</td><td>Unique</td></tr><tr><td>2273</td><td>777777607720731</td><td>1 (0.97)</td><td>25</td><td></td><td>H3</td><td>Unique</td></tr><tr><td>2623</td><td>700000004020771</td><td>1 (0.97)</td><td>33.33</td><td></td><td>H1</td><td>Unique</td></tr><tr><td>2745</td><td>577777607760751</td><td>1 (0.97)</td><td>33.33</td><td></td><td>LAM9</td><td>Unique</td></tr></tbody></table>

(Continued)
<table><caption>Table 3. Description of clusters containing 3 or more isolates in this study, and their worldwide distribution in the SITVIT2 database.</caption><thead><tr><th>SIT (Lineage) Octal Number Spoligotype Description</th><th>Number (%) in study</th><th>% in study vs. database</th><th>Distribution in Regions with ≥3% of a given SIT*</th><th>Distribution in countries with ≥3% of a given SIT**</th></tr></thead><tbody><tr><td>33 (LAM3) 776177607760771</td><td>11 (10.68)</td><td>0,84</td><td>AMER-S 34.96, AFRI-S 24.81, AMER-N 12.14, EURO-S 10.84, EURO-W 6.26, AMER-C 4.12</td><td>ZAF 24.81, PER 16.95, USA 12.14, BRA 10.92, ESP 6.79, ARG 4.35, FXX 4.05, ITA 3.51, HND 3.28</td></tr><tr><td>53 (T1) 777777777760771</td><td>9 (8.74)</td><td>0,14</td><td>EURO-W 15.0, AMER-S 14.35, AMER-N 12.93, EURO-S 9.02, EURO-N 7.17, ASIA-W 7.0, AFRI-S 4.76, AFRI-E 4.46, ASIA-E 4.09, AFRI-N 3.37, EURO-E 3.13, CARI 3.1, AMER-C 3.1</td><td>USA 12.65, FXX 7.55, BRA 5.62, ITA 5.11, ZAF 4.65, PER 3.74, TUR 3.33, AUT 3.28</td></tr><tr><td>50 (H3) 777777777720771</td><td>8 (7.77)</td><td>0,2</td><td>AMER-S 26.08, EURO-W 14.98, AMER-N 14.98, EURO-S 9.84, CARI 4.96, EURO-E 4.71, EURO-N 4.66, AFRI-N 3.63, AFRI-S 3.45, AFRI-M 3.22</td><td>USA 14.96, PER 13.65, BRA 7.15, FXX 5.87, AUT 5.19, ITA 4.63, ESP 4.63, ZAF 3.45, CMR 3.17, CZE 3.12</td></tr><tr><td>37 (T3) 777737777760771</td><td>7 (6.8)</td><td>1,34</td><td>AFRI-E 19.19, EURO-N 11.9, EURO-W 11.52, ASIA-W 10.56, AMER-S 9.21, AMER-N 8.64, ASIA-E 6.14, EURO-S 5.76, EURO-E 4.42, ASIA-S 3.65, AFRI-S 3.26</td><td>ETH 16.7, USA 7.68, SWE 5.18, SAU 4.99, CHN 4.99, FXX 4.22, ITA 4.03, BRA 3.65, ZAF 3.26, DNK 3.26</td></tr><tr><td>42 (LAM9) 777777607760771</td><td>7 (6.8)</td><td>0,2</td><td>AMER-S 31.16, AMER-N 11.33, EURO-S 10.81, EURO-W 9.04, AFRI-N 8.19, EURO-N 4.68, CARI 4.05, AMER-C 3.4, AFRI-E 3.4, AFRI-S 3.02</td><td>BRA 13.07, USA 11.33, COL 7.28, MAR 6.73, ITA 6.25, FXX 4.82, PER 3.54, ESP 3.2, VEN 3.17, ZAF 3.02</td></tr><tr><td>52 (T2) 777777777760731</td><td>3 (2.91)</td><td>0,31</td><td>EURO-W 18.95, ASIA-E 14.97, EURO-N 13.4, AMER-N 11.62, EURO-S 5.13, ASIA-W 5.13, AFRI-M 5.13, AFRI-E 4.82, EURO-E 4.19, AMER-C 3.35, AMER-S 3.04</td><td>CHN 11.62, USA 11.2, SWE 9.42, FXX 8.69, BEL 4.71, CMR 4.19, ITA 3.35, JPN 3.14, ETH 3.14</td></tr><tr><td>211 (LAM3) 776137607760771</td><td>3 (2.91)</td><td>2,83</td><td>AMER-N 30.19, AMER-C 20.76, EURO-W 19.81, AMER-S 18.87, EURO-S 7.55</td><td>USA 30.19, MEX 20.76, FXX 18.87, BRA 11.32, ESP 4.72</td></tr></tbody></table>

* Worldwide distribution is reported for regions with more than 3% of a given SITs as compared to their total number in the SITVIT2 database. The definition of macro-geographical regions and sub-regions (http://unstats.un.org/unsd/methods/m49/m49regin.htm) is according to the United Nations; Regions: AFRI (Africa), AMER (Americas), ASIA (Asia), EURO (Europe), and OCE (Oceania), subdivided in: E (Eastern), M (Middle), C (Central), N (Northern), S (Southern), SE (South-Eastern), and W (Western). Furthermore, CARIB (Caribbean) belongs to Americas, while Oceania is subdivided in 4 sub-regions, AUST (Australasia), MEL (Melanesia), MIC (Micronesia), and POLY (Polynesia). Note that in our classification scheme, Russia has been attributed a new sub-region by itself (Northern Asia) instead of including it among rest of the Eastern Europe. It reflects its geographical localization as well as due to the similarity of specific TB genotypes circulating in Russia (a majority of Beijing genotypes) with those prevalent in Central, Eastern and South-Eastern Asia.

** The 3 letter country codes are according to http://en.wikipedia.org/wiki/ISO_3166-1 alpha-3; countrywide distribution is only shown for SITs with ≥3% of a given SITs as compared to their total number in the SITVIT2 database. Note that FXX code designates Metropolitan France.

doi:10.1371/journal.pone.0118007.t003
<figure><img src="image_4.png" /><figcaption>**Fig 1.** A minimum spanning tree (MST) illustrating evolutionary relationships between spoligotypes and 12-loci MIRU-VNTRs in Santiago, Chile (n = 103 isolates). (A) MST constructed on all spoligotypes; (B) MST constructed on all 12-loci MIRUs alone; and (C) MST constructed on the combination of spoligotypes and 12-loci MIRUs. The phylogenetic tree connects each genotype based on degree of changes required to go from one allele to another. The structure of the tree is represented by branches (continuous vs. dashed and dotted lines) and circles representing each individual pattern. Note that the length of the branches represents the distance between patterns while the complexity of the lines (continuous, gray dashed and gray dotted) denotes the number of allele/spacer changes between two patterns: solid lines, 1 or 2 or 3 changes (thicker ones indicate a single change, while the thinner ones indicate 2 or 3 changes); gray dashed lines represent 4 changes; and gray dotted lines represent 5 or more changes. The size of the circle is proportional to the total number of isolates in our study, illustrating unique isolates (smaller nodes) versus clustered isolates (bigger nodes). The separation inside circle also indicates the number of strains. The color of the circles indicates the phylogenetic lineage to which the specific pattern belongs. The labels of nodes indicate SITs and 12-MITs respectively in Fig. A and B; and main SIT/12-MIT couples (n ≥ 2 isolates) are represented in Fig. C.</figcaption></figure>

doi:10.1371/journal.pone.0118007.g001

Discussion

In the present study, we provided the first insight into the population structure of *M. tuberculosis* isolates in Santiago, Chile, showing the predominance of both the LAM and T lineages. This lineage distribution does not differ significantly as compared to descriptions reported from other Latin American countries where predominant isolates usually belong to three major clades: LAM, T, and Haarlem; in different proportions depending on specific country and year of isolation, but with generally a clear predominance of LAM (Table 4). Thus, in Brazil, the predominant *M. tuberculosis* lineage isolated between 1996 and 2005 was LAM (46%); followed by T (18.6%) and Haarlem (12.2%) families [21]. Similarly, in Argentina between years 2006 and
<table><caption>Table 4. Main published studies (2006–2013) reporting distribution of *Mycobacterium tuberculosis* lineages in Central and South America.</caption><thead><tr><th>Country, Period (years) and Reference</th><th>Numbe of strains</th><th>Harleem (H)</th><th>T</th><th>LAM</th><th>S</th><th>U</th><th>Beijing</th><th>X</th><th>EAI</th><th>Orphans or Unclassified</th></tr></thead><tbody><tr><td rowspan="4">Argentina, 2006–2007 [22]</td><td rowspan="4">157</td><td>37%</td><td>18%</td><td>23%</td><td rowspan="4">NR</td><td rowspan="4">NR</td><td rowspan="4">NR</td><td rowspan="4">NR</td><td rowspan="4">NR</td><td rowspan="4">NA</td></tr><tr><td>H1 10%</td><td>T1 18%</td><td>LAM9 13%</td></tr><tr><td>H2 16%</td><td></td><td>LAM3 10%</td></tr><tr><td>H3 11%</td><td></td><td></td></tr><tr><td rowspan="3">Argentina, 2003–2009 [51]</td><td rowspan="3">787 (MDR only)</td><td>36.3%</td><td>13.9%</td><td>38.8%</td><td>2.8%</td><td>1.7%</td><td>1.5%</td><td>0.9%</td><td>NR</td><td>3.8%</td></tr><tr><td></td><td>T1 8.9%</td><td>LAM3 16.4%</td></tr><tr><td></td><td></td><td></td></tr><tr><td>Argentina, Brasil, Chile, Colombia, Venezuela, 2004–2008 [52]</td><td>951 (MDR only)</td><td>28.9%</td><td>17.4%</td><td>37.2%</td><td>NR</td><td>NR</td><td>1.3%</td><td>NR</td><td>NR</td><td>NR</td></tr><tr><td>French Guyane, 1994–2003 [34]</td><td>744</td><td>21.7%</td><td>30.1%</td><td>19.7%</td><td>1.5%</td><td>NR</td><td>1.0%</td><td>6.5%</td><td>4.8%</td><td>6.5%</td></tr><tr><td rowspan="2">Brasil, 1996–2005 [21]</td><td rowspan="2">1991</td><td>12.2%</td><td>18.6%</td><td>46%</td><td>1.9%</td><td>NR</td><td>NR</td><td>4.7%</td><td>0.85%</td><td>15%</td></tr><tr><td></td><td></td><td>LAM9 10.3%</td></tr><tr><td>Colombia 1995–2007[23]</td><td>152</td><td>25%</td><td>13.8%</td><td>49.3%</td><td>3.3%</td><td>NR</td><td>0.6%</td><td>1.3%</td><td>NR</td><td>6.6%</td></tr><tr><td rowspan="2">Colombia 2005–2008 [33]</td><td rowspan="2">414</td><td>44.3%</td><td></td><td>38.5%</td><td></td><td></td><td>0</td><td></td><td></td><td></td></tr><tr><td>H1 23.4%</td><td></td><td>LAM9:29.9%</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Paraguay, 2003 [39]</td><td>220</td><td>18.2%</td><td>8.6%</td><td>52.3%</td><td>9.5%</td><td>NR</td><td>0.5%</td><td>0.9%</td><td>NR</td><td>NR</td></tr><tr><td rowspan="4">Peru, 1996–2004 [24]</td><td rowspan="4">391 (MDR only)</td><td>11.4%</td><td>27.5%</td><td>33.3%</td><td>NR</td><td>0.3%</td><td>4.9%</td><td>4.4%</td><td>NR</td><td>18%</td></tr><tr><td>H1 5.7%</td><td>T1 22.1%</td><td>LAM5 12.2%</td><td></td><td></td><td></td><td>X3 4.4%</td><td></td><td></td></tr><tr><td>H3 5.7%</td><td>T2 4.9%</td><td>LAM9 9.9%</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>T5 0.5%</td><td>LAM1 6.5%</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Peru, 2004–2006 [25]</td><td>323</td><td>23.8%</td><td>22.3%</td><td>23.8%</td><td>NR</td><td>NR</td><td>9.3%</td><td>NR</td><td>NR</td><td>13.3%</td></tr><tr><td>Peru, 2009 [26]</td><td>199</td><td>34.9%</td><td>15.6%</td><td>12.6%</td><td>1.0%</td><td>8.5%</td><td>14.1%</td><td>3.5%</td><td>NR</td><td>15.1%</td></tr><tr><td rowspan="2">Peru, 1999–2005[27]</td><td rowspan="2">794</td><td>H3 16.4%</td><td>T1 12.3%</td><td>18.7%</td><td>NR</td><td>NR</td><td>5.5%</td><td>2.8%</td><td>NR</td><td>3.4%</td></tr><tr><td>H1 3.9%</td><td></td><td>LAM3 8.3%</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td></td><td></td><td>LAM9 7.4%</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Venezuela, 1998–1999 [38]</td><td>670</td><td>3%</td><td>13%</td><td>74%</td><td>1.0%</td><td>NR</td><td>0%</td><td>NR</td><td>NR</td><td>9%</td></tr><tr><td>Venezuela, 1997–2006 [40]</td><td>1298</td><td>5%</td><td>10.6%</td><td>53%</td><td>1.9%</td><td>9.7%</td><td>0.4%</td><td>1.2%</td><td>0.2%</td><td>17.9%</td></tr><tr><td rowspan="2">Bolivia, 2010 [30]</td><td rowspan="2">99</td><td>39.4%</td><td>22.2%</td><td>26.3%</td><td>2%</td><td>NR</td><td>NR</td><td>1%</td><td>NR</td><td>9.1%</td></tr><tr><td>H3 30.3%</td><td>T1 21.2%</td><td>LAM3 15.2%</td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

doi:10.1371/journal.pone.0118007.t004

2007, the most common lineages were also LAM (50%), followed by T1 (18%) and Haarlem 1 (10%) and in Colombia, for strains isolated between years 1995 and 2007, LAM was present in 49.3% of isolates; Haarlem, in 25% and T group, in 13.8% [22] [23]. On the contrary, in Peru, which is Chile's northern neighbor country, this proportion is lost as LAM strains have been found in lower percentage ranging from only 12.6 to 33.3% [24] [25] [26] [27].

A particular characteristic that has been ascribed to the LAM family is that they are unique to harbor a significant proportion of strains with a specific deletion named RD<sup>Rio</sup> characterized by a single long sequence polymorphism (>26.3 kb) with deletion or modification of 10 genes [28]. LAM1 (SIT20) and LAM2 (SIT17) are exclusively of the RD<sup>Rio</sup> genotype, LAM4, LAM5,
LAM6 and LAM9 are composed of both RD<sup>Rio</sup> and wild type genotypes, and LAM3 composed exclusively of wild type isolates. This RD<sup>Rio</sup> genotype was found in 26.5% (72/270) of all tuberculosis cases in a recent large retrospective study from Brazil, extrapolating that this genotype could be widely spread in that country. RD<sup>Rio</sup> *M. tuberculosis* genotypes were suspected of being more virulent and able to cause disease more efficiently, however, not all researchers have supported this view, and a large survey found disease caused by this strain was not clinically distinctive or more severe than disease caused by non-RD<sup>Rio</sup> [29].

In contrast with other Latin-American countries where LAM is the most prevalent lineage, a shift towards increased proportion of Haarlem lineage was observed in a recent survey of 199 isolates carried out in the North-eastern part of Lima, Peru, where H strains were found in a predominant proportion of 29.6%, followed by T at 15.6%, Beijing at 14.1% and LAM at only 12.6% [26]. Similarly, Bolivia also recently reported a predominance of H genotype in up to 40% of isolates studied in 2010. [30] This Haarlem genotype has the particular significance of having been involved in several outbreaks of MDR-tuberculosis [31] [32]. Furthermore, a recent study conducted in Colombia with strains isolated between 2005 and 2008 has also shown an increase in Haarlem strains actually reaching 44.3% in the city of Medellín, and included a large proportion of MDR among strains belonging to SIT45/H1 [33].

Descriptions for main circulating strains in South-American countries have long described the Euro-American lineage to be widely predominant, and this has been used to support a European dissemination from either early settlement or trade associations [34]. It is of interest to note that this model, however, has been unable to explain the profuse archaeological evidence for the presence of tuberculosis in the Americas before European contact and recent comparative genomics from 1,000-year-old mycobacterial found in Peruvian human skeletons described members of the MTBC most closely related to those adapted to seals and sea lions than those from humans [35]. Worldwide, lineage-specific differences in the virulence and transmissibility of clinical isolates have been reported across independent experimental systems with modern lineages, such as Beijing and Euro-American Haarlem strains believed to exhibit more virulent phenotypes compared with ancient lineages, such as East-African-Indian and *M. africanum* [36]. Strains belonging to the Beijing genotype have been associated with increased transmission risk, large outbreaks and drug resistant TB in many parts of the world. This genotype was only recently described for the first time in Chile and is infrequently observed in TB patients in South America with the exception of Peru [37]. On one hand, in Argentina, Brazil, Colombia, Chile, Ecuador, Paraguay and Venezuela, the prevalence of Beijing lineage strains reported has been generally lesser than 2% [21] [23] [28] [37] [38] [39] [40] [41]. On the other hand, the prevalence of Beijing genotype in Peru is exceptionally high and has been estimated between 4.9% and 14.1%, although not clearly associated with MDR tuberculosis [24] [25] [26] [27]. Ritacco et al. speculated that the Beijing lineage strains were first introduced into Peru, and eventually into other South American countries, when Peru received a significant number of Chinese immigrants in the mid-19th century. The need for vigilance is further warranted by another recently published study that included a total of 200 Beijing strains from Peru where 24 loci VNTR analysis showed a high clustering rate (80.3%) and a high recent transmission index (RTI<sub>n-1</sub> = 0.707), strongly suggesting active and on-going transmission of Beijing lineage strains in the surveyed area [42].

In the present study, susceptibility data was only available for a minority of isolated strains, as susceptibility determination was not required routinely for newly acquired tuberculosis due to the general low MDR rate for the country. With regard to the possibility of strains-other than Beijing- potentially associated with drug resistance, it is of interest to note that in a recent large survey of phylogenetic clades collected over a seventeen-year period in the French Department of the Americas, X and LAM lineages were overrepresented in drug-resistant and MDR-TB cases, respectively, and also, four predominant spoligotypes were significantly
associated with drug resistance corresponding to SIT20/LAM1, SIT64/LAM6, SIT45/H1 and SIT46/undefined lineage [43].

The effects of migration on the epidemiology of tuberculosis has been evaluated by many studies by determining genotypic lineages on the isolates and also providing general information on the phylogeographical origin of the strains circulating in the study area [44] [45]. In the present study, 11.6% of strains belonged to migrant population, which is a slightly higher proportion in relation to the national reports (7%), but consistent with our patients who were recruited in the central area of Santiago where the highest immigration rate is found [46]. Among this group, the most prevalent *M. tuberculosis* lineage found was also LAM, although infrequently isolated strains such as CAS—belonging to a migrant from Nepal—and a Beijing strain belonging to a migrant from Peru were also detected. As there is considerable evidence that genetic susceptibility can influence on tuberculosis acquisition risk, such as the higher susceptibility conferred by individual polymorphisms in genes for natural resistance-associated macrophage protein (NRAMP1), vitamin D receptor, or mannose-binding protein [47] [48], it has been suggested that different ethnic groups may have different susceptibility for particular *M. tuberculosis* strains [49]. This will be an important issue to further explore in Chile and among other populations from Latin America where different ethnic backgrounds coexist and where furthermore, most of ethnic groups are concentrated in areas of poverty and crowding, adding strong confounding factors that can determine a higher prevalence of tuberculosis disease among the most vulnerable.

Finally, it is important to mention that both genotyping methods used in present study are based on analysis of a limited number of loci, and it has been described that VNTR loci as well as particular spoligotyping patterns can exhibit variable discriminatory power in different mycobacterial lineages. Rapid markers evolution with a tendency to converge can lead to significant amounts of homoplasy, and hence the limitations of these tools in defining deep phylogenetic groups in *M. tuberculosis* complex or other bacteria [50]. This can explain discordant findings for T and H lineages classification between both SITVIT2 and MIRU-VNTRplus databases (S2 Table). Future investigations by use of more discriminatory tools such as whole genome sequencing (WGS) might be helpful to identify true phylogenetic relationships and also obtain a more precise evaluation of TB transmission.

Among orphan isolates, we noted similar spoligotyping patterns (with absence of spacer 6) among orphans Or08, Or11, and Or12 within our study (S1 Table). The ancestral genotypes found among orphan patterns in our study are most probably suggestive of ancestral lineages undergoing extinction after an initial distant contact with patient(s) infected with Africanum lineage or the African continent. This finding may suggest a probable ancient link with Africa-num lineage and the African continent. Interestingly thanks to the "SpolSimilaritySearch tool" under development in SITVIT2 database, we detected that a similar spoligotype was found once in Italy: ■■■■□■□□□□■■■■■■■□□□□■■■■■■■■■■■□□■■■■; octal code 76407760777071; AFRI_2. However, one cannot exclude the possibility of a recent introduction of such strains in Chile.

A supplemental figure (S2 Fig.) has been added to support the finding of high presence of SIT37/T3 in Chile. This figure shows geographical maps based on SITVIT2 interrogation prior to entry of present dataset to summarize percentages of SIT37 (A) and strains belonging to T3 sub-lineage (B), by country. Note that this distribution should not be confounded with data in Table 3, which shows proportions of a given spoligotype in the global database. For example, we have 7 strains for SIT37/T3 representing 6.8% of our study sample. However when checked in the database, these 7 strains represent 1.34% of all SIT37/T3 strains present in SITVIT2. Finally, if we look for all identical strains in the database, we see that 19.19% were reported from AFRI-E, 11.9% from EURO-N, 11.52% from EURO-W, etc. Finally, we have also added a new
map (S3 Fig.) to show the global distribution of lineages in USA, Canada, and various Latin American and Caribbean countries as compared to Chile.

In conclusion, we have shown the predominance of LAM and T genotypes corresponding to almost two thirds of circulating stains for the central area of Santiago, Chile, in the study period. We have also underlined the significant presence of SIT37/T3 among Chilean isolates. This high proportion of SIT37/T3 in our study in Chile as well as in other countries of Latin America (Guyana, Cuba, or Dominican Republic) in conjunction with East Africa, and to a lesser extent in other regions and countries—including not negligible proportions in Nepal and China—suggests long-lasting spread of this genotype. Obviously, more data and analysis are needed to clarify probable ancestral human migrations responsible for its actual distribution pattern. Further studies on a larger number of isolates and from patients with larger geographic and ethnic diversity might help to accurately estimate the genetic landscape of prevailing *M. tuberculosis* epidemic in Chile.

**Supporting Information**

**S1 Fig. Unweighted Pair Group Method with Arithmetic Mean (UPGMA) tree based on the combination of spoligotypes and 12-loci MIRUs (using MLVA Compare software).** Labels identifying all strains are as follows: strain number, lineage, SIT, 12-MIT. (TIF)

**S2 Fig. Geographical maps (obtained by interrogating SITVIT2 database prior to entry of this study) representing percentage distributions of strains belonging to SIT37 (A), and strains belonging to T3 sub-lineage (B).** Maps were reproduced from work created and shared by Google (https://google-developers.appspot.com/site-policies) and used according to terms described in the Creative Commons 3.0 Attribution License (http://creativecommons.org/licenses/by/3.0/). (TIF)

**S3 Fig. Geographic distribution of lineages in Latin American and neighboring countries.**
Note that the map file was downloaded under Creative Commons License using the link: http://en.wikipedia.org/wiki/Latin_America and was manually modified for representative purposes only.
(TIF)

**S1 Table. Description of spoligotypes and MIRU12 patterns corresponding to orphans strains (n = 16) in the SITVIT2 database.**
(DOCX)

**S2 Table. Description of spoligotypes, MIRU12 patterns and lineage allocation by SITVIT2 and MIRU-VNTRplus for all the isolates of present study (n = 103).**
(PDF)

**Acknowledgments**

We acknowledge the Hospital San Borja Arriarán and Centro Médico San Joaquín microbiology laboratory personnel who collected and shared the samples used for these analyses.

**Author Contributions**

Conceived and designed the experiments: MEB PG PM CP. Performed the experiments: PG PM MC. Analyzed the data: MEB PG PM DC NR. Contributed reagents/materials/analysis tools: PM MC CP. Wrote the paper: MEB PG DC NR.
References

1. WHO report 2013. "Global Tuberculosis Control". Available: http://www.who.int/tb/publications/global_report/en/http://www.who.int/tb/publications/global_report/en/. Accessed 2014 October 1.

2. Gutierrez MC, Brisse S, Brosch R (2005) Ancient origin and gene mosaicism of the progenitor of Mycobacterium tuberculosis. PLoS pathogens 1 (e5), 55-61 PMID: 16201017

3. Walker TM, Clp C, Harrell R, Evans JT, Kapatai G, et al. (2013) Whole-genome sequencing to delineate Mycobacterium tuberculosis outbreaks: a retrospective observational study. Lancet Infectious Diseases; 13: 137-146

4. Gagneux S, DeRiemer K, Van T, Kato-Maeda M, de Jong BC, et al. (2006) Variable host-pathogen compatibility in Mycobacterium tuberculosis. Proc. Natl Acad. Sci. USA 103, 2869-2873. PMID: 16477032

5. Sreevatsan S, Pan X, Stockbauer KE, Connell ND, Kreiswirth BN, et al. (1997) Restricted structural gene polymorphism in the *Mycobacterium tuberculosis* complex indicates evolutionarily recent global dissemination. Proc Natl Acad Sci USA. 94: 9869-9874. PMID: 9275218

6. Couvin D, Rastogi N (2014) The establishment of databases on circulating genotypes of *Mycobacterium tuberculosis* complex and web tools for an effective response to better monitor, understand and control the tuberculosis epidemic worldwide. EuroReference—Journal of Reference, ISSN 2110-5294. 12:36-48. Available: https://pro.anses.fr/euroreference/Documents/ER12-RESEAUX-MycobacteriumEN.pdf. doi: 10.1186/1471-2407-14-960 PMID: 25523972

7. García de Viedma D, Mokrousov I, Rastogi N (2011) Innovations in the molecular epidemiology of tuberculosis. Enfermedades Infecciosas y Microbiologia Clinica 29 (Suppl.1), 8-13. doi: 10.1016/S0213-005X(11)70012-X PMID: 21420561

8. Kato-Maeda M, Kim EY, Flores L, Jarlsberg LG, Osmond D, et al. (2010) Differences among sublineages of the East-Asian lineage of Mycobacterium tuberculosis in genotypic clustering. Int. J. Tuberc. Lung Dis. 14, 538-544. PMID: 20392345

9. Ford CB, Shah RR, Kato-Maeda M, Gagneux S, Murray MB, et al. (2013) Mycobacterium tuberculosis mutation rate estimates from different lineages predict substantial differences in the emergence of drug-resistant tuberculosis. Nat. Genet. 45, 784-790. doi: 10.1038/ng.2656 PMID: 23749189

10. Mancilla M, Martínez A, Palavecino C, Rehren G, Lucero P, et al. (2006) [Genetic variants of Mycobacterium tuberculosis isolated from patients of the Xth Region of Chile]. Rev Chilena Infectol; 23 (3): 220-225. PMID: 16896494

11. Meza P, Balcells ME, Miranda C, Cifuentes M, Wozniak A, et al. (2014) [Presence of Beijing genotype among Mycobacterium tuberculosis strains in two centres of the Region Metropolitana of Chile]. Rev Chilena Infectol; 31 (1): 21-27. doi: 10.4067/S0716-10182014000100003 PMID: 24740770

12. Laszlo A, Rahman M, Raviglione M, Bustreo F (1997) Quality assurance programme for drug susceptibility testing of Mycobacterium tuberculosis in the WHO/IUATLD Supranational Laboratory Network: first round of proficiency testing. Int J Tuberc Lung Dis. 1, 231-238. PMID: 9432369

13. Kasai H, Ezaki T, Harayama S (2000) Differentiation of phylogenetically related slowly growing mycobacteria by their gyrB sequences. J Clin Microbiol; 38 (1): 301-308. PMID: 10618105

14. Kamerbeek J, Schouls L, Kolk A, van Agterveld M, van Soolingen D, et al. (1997) Simultaneous detection and strain differentiation of Mycobacterium tuberculosis for diagnosis and epidemiology. Journal of Clinical Microbiology; 35 (4): 907-914. PMID: 9157152

15. Abadia E, Zhang J, Ritacco V, Kremer K, Ruimy R, et al. (2011) The use of microbead-based spoligotyping for *Mycobacterium tuberculosis* complex to evaluate the quality of the conventional method: providing guidelines for Quality Assurance when working on membranes. BMC Infect Dis; 11: 110. doi: 10.1186/1471-2334-11-110 PMID: 21527037

16. Cowan LS, Mosher L, Diem L, Massey JP, Crawford JT (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(5): 1592-602 PMID: 11980927

17. Supply P, Allix C, Lesjean S, Cardoso-Oelemann M, Rusch-Gerdes S, et al. (2006) Proposal for standardization of optimized mycobacterial interspersed repetitive unit-variable-number tandem repeat typing of Mycobacterium tuberculosis. J Clin Microbiol 44: 4498-4510. PMID: 17005759

18. Demay C, Liens B, Burguière T, Hill V, Couvin D, et al. (2012) SITVITWEB—A publicly available international multimarker database for studying Mycobacterium tuberculosis genetic diversity and molecular epidemiology. Infect Genet Evol. 12: 755-766. doi: 10.1016/j.meegid.2012.02.004 PMID: 22365971

19. Shabbeer A, Ozcaglar C, Yener B, Bennett KP (2012) Web tools for molecular epidemiology of tuberculosis. Infect Genet Evol; 12(4): 767-81. doi: 10.1016/j.meegid.2011.08.019 PMID: 21903179
20. 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(Web Server issue): W326–W331. doi: 10.1093/nar/gkq351 PMID: 20457747

21. Gomes HM, Elias AR, Oelemann MA, Pereira MA, Montes FF, et al. (2012) Spoligotypes of Mycobacterium tuberculosis complex isolates from patients residents of 11 states of Brazil. Infect Genet Evol 12: 649–656. doi: 10.1016/j.meegid.2011.08.027 PMID: 21907830

22. Gonzalo X, Ambroggi M, Cordova E, Brown T, Poggi S, et al. (2011) Molecular Epidemiology of *Mycobacterium tuberculosis*, Buenos Aires, Argentina. Emerging Infectious Diseases; 17(3): 528–531. doi: 10.3201/eid1703.100394 PMID: 21392451

23. Cerezo I, Jimenez Y, Hernandez J, Zozio T, Murcia MI, et al. (2012) A first insight on the population structure of Mycobacterium tuberculosis complex as studied by spoligotyping and MIRU-VNTRs in Bogota, Colombia. Infect Genet Evol 12: 657–663. doi: 10.1016/j.meegid.2011.07.006 PMID: 21784173

24. Cohen T, Murray M, Abubakar I, Zhang Z, Sloutsky A, et al. (2011) Multiple introductions of multidrug-resistant tuberculosis into households, Lima. Peru. Emerg Infect Dis, 17(6):969–975. doi: 10.3201/eid/1706.101471 PMID: 21749756

25. Taype CA, Agapito JC, Accinelli RA, Espinoza JR, Godreuil S, et al. (2012) Genetic diversity, population structure and drug resistance of Mycobacterium tuberculosis in Peru. Infect Genet Evol 12: 577–585. doi: 10.1016/j.meegid.2012.02.002 PMID: 22342744

26. Barletta F, Otero L, Collantes J, Asto B, de Jong BC, et al. (2013) Genetic variability of Myobacterium tuberculosis complex in patients with no known risk factors for MDR-TB in the North-eastern part of Lima, Peru. BMC Infectious Diseases; 13: 397–343. doi: 10.1186/1471-2334-13-397 PMID: 23984854

27. Sheen P, Couvin D, Grandjean L, Zimic M, Dominguez M, et al. (2013) Genetic Diversity of Mycobacterium tuberculosis in Peru and Exploration of Phylogenetic Associations with Drug Resistance. PLoS ONE 8(6): e65873. doi: 10.1371/journal.pone.0065873 PMID: 23826083

28. Lazzarini LC, Hurad RC, Boechat NI, Gomes HM, Oelemann MC, et al. (2007) Discovery of a novel Mycobacterium tuberculosis lineage that is a major cause of tuberculosis in Rio de Janeiro, Brazil. Journal of Clinical microbiology 45, 3891–3902. PMID: 17898156

29. Barbosa CB, Lazzarini LCO, Elias AR, Leung JAM, Ribeiro SB (2012) Tuberculosis caused by RD<sup>Rio</sup> *Mycobacterium tuberculosis* is not associated with differential clinical features. Int J Tuberc Lung Dis; 16 (10), 1377–1382. doi: 10.5588/ijtld.11.0709 PMID: 22863208

30. Monteserin J, Camacho M, Barrera L, Palomino JC, Ritacco V, et al. (2013) Genotypes of Mycobacterium tuberculosis in patients at risk of drug resistance in Bolivia. Infection, Genetics and Evolution 17; 195–201. doi: 10.1016/j.meegid.2013.04.010 PMID: 23603419

31. Mardassi H, Namouchi A, Haltiti R, Zarrouk M, Mhenni B, et al. (2005) Tuberculosis due to resistant Haarlem strain, Tunisia. Emerg Infect Dis; 11:957–961. PMID: 15963297

32. Jagielski T, Augustynowicz-Kopeć E, Zozio T, Rastogi N, Zwolska Z (2010) Spoligotype-based comparative population structure analysis multidrug-resistant and isoniazid-monoresistant Mycobacterium tuberculosis complex clinical isolates in Poland. J Clin Microbiol; 48(11):3899–3909. doi: 10.1128/JCM.00572-10 PMID: 20810763

33. Realpe T, Correa N, Rozo JC, Ferro BE, Gomez V, et al. (2014) Population Structure among Mycobacterium tuberculosis Isolates from Pulmonary Tuberculosis Patients in Colombia. PLoS ONE 9(4): e93848. doi: 10.1371/journal.pone.0093848 PMID: 24747767

34. Brudey K, Filliol I, Ferdinand S, Guernier V, Duval P, et al. (2006) Long-term population-based genotyping study of *Mycobacterium tuberculosis* complex isolates in the French Departments of the Americas. J Clin Microbiol; 44:183–191. PMID: 16390968

35. Bos KI, Harkins KM, Herbig A, Coscolla M, Weber N, et al. (2014) Pre-Columbian mycobacterial genomes reveal seals as a source of New World human tuberculosis. Nature. doi: 10.1038/nature13591 PMID: 25533962

36. Reiling N, Homolka S, Walter K (2013) Clade-specific virulence patterns of Mycobacterium tuberculosis. MBio 4 e00250–13. doi: 10.1128/mBio.00250-13 PMID: 23900170

37. Ritacco V, Lopez B, Cafrune PI, Ferrazoli L, Suffys PN, et al. (2008) Mycobacterium tuberculosis of the Beijing genotype is rarely observed in tuberculosis patients in South America. Mem Inst Osw Cruz, Rio de Janeiro, Vol 103(5) 489–92. PMID: 18797764

38. Aristimuño L, Armengol R, Cebollada A, España M, Guilarte A, et al. (2006) Molecular characterisation of Mycobacterium tuberculosis isolates in the First National Survey of Anti-tuberculosis Drug Resistance from Venezuela. BMC Microbiol. 6, 90. PMID: 17032442

39. Candia N, Lopez B, Zozio T, Carrivale M, Diaz C, et al. (2007) First insight into Mycobacterium tuberculosis genetic diversity in Paraguay. BMC Microbiol 7: 75. PMID: 17686181
40. Abadia E, Sequera M, Ortega D, Mendez MV, Escalona A, et al. (2009) Mycobacterium tuberculosis ecology in Venezuela: epidemiologic correlates of common spoligotypes and a large clonal cluster defined by MIRU-VNTR-24. BMC Infect Dis 9: 122. doi: 10.1186/1471-2334-9-122 PMID: 19660112

41. Morcillo N, Di Giulio B, Chirico C, Kuriger A, Dolmann A, et al. (2005) First description of Mycobacterium tuberculosis Beijing genotype in Argentina. Rev. Argent. Microbiol. 37, 92–95. PMID: 16178465

42. Iwamoto T, Grandjean L, Arikawa K, Nakanishi N, Caviedes L, et al. (2012) Genetic Diversity and Transmission Characteristics of Beijing Family Strains of *Mycobacterium tuberculosis* in Peru. PLoS One; 7(11): e49651. doi: 10.1371/journal.pone.0049651 PMID: 23185395

43. Millet J, Streit E, Berchel M, Bomer A-G, Schuster F, et al. (2014) A Systematic Follow-Up of *Mycobacterium tuberculosis* Drug-Resistance and Associated Genotypic Lineages in the French Departments of the Americas over a Seventeen-Year Period. BioMed Research International; Volume 2014, Article ID 689852, 14 pages. Available: http://dx.doi.org/10.1155/2014/689852.

44. Filliol L, Driscoll JR, van Soolingen D, Kreiswirth BN, Kremer K, et al. (2003) Snapshot of moving and expanding clones of Mycobacterium tuberculosis and their global distribution assessed by spolygotyping in an international study. J Clin Microbiol; 41, 1963–1971. PMID: 12734235

45. Garzelli C, Rindi L (2012) Molecular epidemiologic approaches to study the epidemiology of tuberculosis in low-incidence settings receiving immigrants. Infection, Genetic and Evolution; 12: 610–18. doi: 10.1016/j.meegid.2011.10.015 PMID: 22041505

46. Tuberculosis Informe De Situación Chile (2012) Programa Nacional de Tuberculosis. MINSAL, Chile. Available: http://web.minsal.cl/sites/default/files/files/tuberculosisInforme2012.pdf. Accessed 2014 October 1.

47. North RJ, Jung Y- J (2004) Immunity to tuberculosis. Annu Rev Immunol; 22:599–623. PMID: 15032590

48. van der Eijk EA, van de Vosse E, Vandenbroucke JP, van Dissel JT (2007) Heredity versus Environment in Tuberculosis in Twins: The 1950s United Kingdom Prophit Survey Simonds and Comstock Revisited. Am J Respir Crit Care Med; 176:1281–1288. PMID: 17823356

49. Gagneux S, Small PM (2007) Global phylogeography of Mycobacterium tuberculosis and implications for tuberculosis product development. Lancet Infect Dis 7(5):328–337. PMID: 17448936

50. 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(11): e7815. doi: 10.1371/journal.pone.0007815 PMID: 19915672

51. Ritacco V, López B, Ambroggi M, Palmero D, Salvadores B, et al. (2012) HIV Infection and Geographically Bound Transmission of Drug-Resistant Tuberculosis, Argentina. Emerg Infect Dis;18 (11):1802–1810. doi: 10.3201/eid1811.120126 PMID: 23092584

52. Ritacco V, Iglesias MJ, Ferrazoli L, Monteserin J, Dalla Costa ER, et al. (2012) Conspicuous multidrug-resistant Mycobacterium tuberculosis cluster strains do not trespass country borders in Latin America and Spain. Infect Genet Evol;12:711–717. doi: 10.1016/j.meegid.2011.06.006 PMID: 21718805
<figure><img src="image_5.png" /></figure>

UPGMA tree based on 55 columns, no missing values, % columns difference

## Supporting Information

S1 Fig. Unweighted Pair Group Method with Arithmetic Mean (UPGMA) tree based on the combination of spoligotypes and 12-loci MIRUs (using MLVA Compare software). Labels identifying all strains are as follows: strain number, lineage, SIT, 12-MIT. Larger TIF image available online at:

http://journals.plos.org/plosone/article/asset?unique&id=info:doi/10.1371/journal.pone.0118007.s001
Percentage of SIT37

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

Percentage of isolates belonging to
T3 sub-lineage

<figure><img src="image_7.png" /><figcaption>**Supplemental Figure S2:** Geographical maps (obtained by interrogating SITVIT2 database prior to entry of this study) representing percentage distributions of strains belonging to SIT37 (A), and strains belonging to T3 sub-lineage (B). Maps were reproduced from work created and shared by Google (https://google-developers.appspot.com/site-policies) and used according to terms described in the Creative Commons 3.0 Attribution License (http://creativecommons.org/licenses/by/3.0/).</figcaption></figure>

**Supporting Information**

S2 Fig. Larger TIF image available online at:

http://journals.plos.org/plosone/article/asset?unique&id=info:doi/10.1371/journal.pone.0118007.s002
<figure><img src="image_8.png" /></figure>

## Supporting Information

S3 Fig. Geographic distribution of lineages in Latin American and neighboring countries. Note that the map file was downloaded under Creative Commons License using the link: http://en.wikipedia.org/wiki/Latin_America and was manually modified for representative purposes only. Larger TIF image available online at: http://journals.plos.org/plosone/article/asset?unique&id=info:doi/10.1371/journal.pone.0118007.s003
<table><caption>Supplemental Table S2. Description of spoligotypes, MIRU12 patterns and lineage allocation by SITVIT2 and MIRU-VNTRplus for all the isolates of present study (n = 103).</caption><thead><tr><th>IsolNumber</th><th>Year</th><th>Spoligotype Description</th><th>Octal code</th><th>Sublineage</th><th>Lineage</th><th>SIT</th><th>MIRU12</th><th>12-MIT</th><th>Lineage (MIRU-VNTR plus)</th><th>Sex / Age</th></tr></thead><tbody><tr><td>CHL022011000041</td><td>2011</td><td></td><td>000000000003771</td><td>Beijing</td><td>Beijing</td><td>1</td><td>224325173533</td><td>229</td><td>Beijing</td><td>M/56</td></tr><tr><td>CHL022011000032</td><td>2011</td><td></td><td>000000000003771</td><td>Beijing</td><td>Beijing</td><td>1</td><td>223325173333</td><td>721</td><td>Beijing</td><td>F/91</td></tr><tr><td>CHL022012000060</td><td>2012</td><td></td><td>677777607760771</td><td>LAM1</td><td>LAM</td><td>20</td><td>224326163321</td><td>375</td><td>LAM</td><td>M</td></tr><tr><td>CHL022011300033</td><td>2011</td><td></td><td>776177607760771</td><td>LAM3</td><td>LAM</td><td>33</td><td>224326133323</td><td>251</td><td>LAM</td><td>M/31</td></tr><tr><td>CHL022012100068</td><td>2012</td><td></td><td>776177607760771</td><td>LAM3</td><td>LAM</td><td>33</td><td>224326133323</td><td>251</td><td>LAM</td><td>M/24</td></tr><tr><td>CHL022011000023</td><td>2011</td><td></td><td>776177607760771</td><td>LAM3</td><td>LAM</td><td>33</td><td>22426133323</td><td>880</td><td>LAM</td><td>F/29</td></tr><tr><td>CHL022010000008</td><td>2010</td><td></td><td>776177607760771</td><td>LAM3</td><td>LAM</td><td>33</td><td>224326133322</td><td>1670</td><td>LAM</td><td>F/63</td></tr><tr><td>CHL022013000091</td><td>2013</td><td></td><td>776177607760771</td><td>LAM3</td><td>LAM</td><td>33</td><td>224326133322</td><td>1670</td><td>LAM</td><td>M/27</td></tr><tr><td>CHL022013000092</td><td>2013</td><td></td><td>776177607760771</td><td>LAM3</td><td>LAM</td><td>33</td><td>224327153324</td><td>1678</td><td>LAM</td><td>F/29</td></tr><tr><td>CHL022011000038</td><td>2011</td><td></td><td>776177607760771</td><td>LAM3</td><td>LAM</td><td>33</td><td>22426133032</td><td>Or01</td><td>LAM</td><td>F/36</td></tr><tr><td>CHL022013000073</td><td>2013</td><td></td><td>776177607760771</td><td>LAM3</td><td>LAM</td><td>33</td><td>224327133324</td><td>Or02</td><td>LAM</td><td>M/56</td></tr><tr><td>CHL022013000074</td><td>2013</td><td></td><td>776177607760771</td><td>LAM3</td><td>LAM</td><td>33</td><td>224126133323</td><td>Or03</td><td>LAM</td><td>M/31</td></tr><tr><td>CHL022013000084</td><td>2013</td><td></td><td>776177607760771</td><td>LAM3</td><td>LAM</td><td>33</td><td>224317133325</td><td>Or04</td><td>LAM</td><td>F/53</td></tr><tr><td>CHL022013000094</td><td>2013</td><td></td><td>776177607760771</td><td>LAM3</td><td>LAM</td><td>33</td><td>224316233325</td><td>Or05</td><td>LAM</td><td>F/27</td></tr><tr><td>CHL022012000050</td><td>2012</td><td></td><td>776377777760771</td><td>S</td><td>S</td><td>34</td><td>233324153224</td><td>Or06</td><td>S</td><td>M/59</td></tr><tr><td>CHL022012100065</td><td>2012</td><td></td><td>776377777760771</td><td>S</td><td>S</td><td>34</td><td>233225143322</td><td>Or07</td><td>S</td><td>M/45</td></tr><tr><td>CHL022013000081</td><td>2013</td><td></td><td>777737777760771</td><td>T3</td><td>T</td><td>37</td><td>224325143325</td><td>220</td><td>X</td><td>F/60</td></tr><tr><td>CHL022013100042</td><td>2013</td><td></td><td>777737777760771</td><td>T3</td><td>T</td><td>37</td><td>223321153324</td><td>585</td><td>H</td><td>M/52</td></tr><tr><td>CHL022013100043</td><td>2013</td><td></td><td>777737777760771</td><td>T3</td><td>T</td><td>37</td><td>223321153324</td><td>585</td><td>H</td><td>F/51</td></tr><tr><td>CHL022012000051</td><td>2012</td><td></td><td>777737777760771</td><td>T3</td><td>T</td><td>37</td><td>223321153324</td><td>585</td><td>H</td><td>M/3</td></tr><tr><td>CHL022012000058</td><td>2012</td><td></td><td>777737777760771</td><td>T3</td><td>T</td><td>37</td><td>223321153324</td><td>585</td><td>H</td><td>F</td></tr><tr><td>CHL022013000095</td><td>2013</td><td></td><td>777737777760771</td><td>T3</td><td>T</td><td>37</td><td>223321153324</td><td>585</td><td>H</td><td>M/50</td></tr><tr><td>CHL022013000078</td><td>2013</td><td></td><td>777737777760771</td><td>T3</td><td>T</td><td>37</td><td>224322153324</td><td>1677</td><td>H</td><td>M</td></tr><tr><td>CHL022011000017</td><td>2011</td><td></td><td>777777347760471</td><td>T4-CEU1</td><td>T</td><td>39</td><td>224325143324</td><td>117</td><td>X</td><td>M/28</td></tr><tr><td>CHL022011000030</td><td>2011</td><td></td><td>777777347760471</td><td>T4-CEU1</td><td>T</td><td>39</td><td>224325143224</td><td>Or08</td><td>LAM</td><td>M/29</td></tr><tr><td>CHL022010000009</td><td>2010</td><td></td><td>777777607760771</td><td>LAM9</td><td>LAM</td><td>42</td><td>124328153326</td><td>180</td><td>LAM</td><td>F/93</td></tr><tr><td>CHL022011000021</td><td>2011</td><td></td><td>777777607760771</td><td>LAM9</td><td>LAM</td><td>42</td><td>124326153220</td><td>190</td><td>LAM</td><td>F/42</td></tr><tr><td>CHL022011000025</td><td>2011</td><td></td><td>777777607760771</td><td>LAM9</td><td>LAM</td><td>42</td><td>124326153220</td><td>190</td><td>LAM</td><td>M/88</td></tr><tr><td>CHL022011000027</td><td>2011</td><td></td><td>777777607760771</td><td>LAM9</td><td>LAM</td><td>42</td><td>124326153220</td><td>190</td><td>LAM</td><td>M/44</td></tr><tr><td>CHL022009100003</td><td>2009</td><td></td><td>777777607760771</td><td>LAM9</td><td>LAM</td><td>42</td><td>124261413220</td><td>Or09</td><td>LAM</td><td>F/74</td></tr><tr><td>CHL022009000010</td><td>2009</td><td></td><td>777777607760771</td><td>LAM9</td><td>LAM</td><td>42</td><td>113326153226</td><td>Or10</td><td>LAM</td><td>F/50</td></tr><tr><td>CHL022011000018</td><td>2011</td><td></td><td>777777607760771</td><td>LAM9</td><td>LAM</td><td>42</td><td>22326163226</td><td>Or11</td><td>LAM</td><td>M/66</td></tr><tr><td>CHL022012000052</td><td>2012</td><td></td><td>77777777720731</td><td>H3</td><td>H</td><td>49</td><td>225325153323</td><td>45</td><td>H</td><td>M/24</td></tr><tr><td>CHL022013100044</td><td>2013</td><td></td><td>77777777720771</td><td>H3</td><td>H</td><td>50</td><td>224325153322</td><td>32</td><td>H</td><td>M/35</td></tr><tr><td>CHL022013100047</td><td>2013</td><td></td><td>77777777720771</td><td>H3</td><td>H</td><td>50</td><td>225325153323</td><td>45</td><td>H</td><td>M/44</td></tr><tr><td>CHL022013100086</td><td>2013</td><td></td><td>77777777720771</td><td>H3</td><td>H</td><td>50</td><td>221326153323</td><td>234</td><td>H</td><td>M/68</td></tr><tr><td>CHL022013000083</td><td>2013</td><td></td><td>77777777720771</td><td>H3</td><td>H</td><td>50</td><td>223326153323</td><td>236</td><td>H</td><td>M/39</td></tr><tr><td>CHL022013000079</td><td>2013</td><td></td><td>77777777720771</td><td>H3</td><td>H</td><td>50</td><td>222326153323</td><td>331</td><td>H</td><td>M/20</td></tr><tr><td>CHL022013000080</td><td>2013</td><td></td><td>77777777720771</td><td>H3</td><td>H</td><td>50</td><td>222326153323</td><td>331</td><td>H</td><td>M/29</td></tr><tr><td>CHL022012000049</td><td>2012</td><td></td><td>77777777720771</td><td>H3</td><td>H</td><td>50</td><td>225326153221</td><td>1673</td><td>H</td><td>M/52</td></tr><tr><td>CHL022013000075</td><td>2013</td><td></td><td>77777777720771</td><td>H3</td><td>H</td><td>50</td><td>221325153323</td><td>Or12</td><td>H</td><td>M/25</td></tr><tr><td>CHL022011000037</td><td>2011</td><td></td><td>77777777760731</td><td>T2</td><td>T</td><td>52</td><td>223125153326</td><td>733</td><td>S</td><td>F/71</td></tr><tr><td>CHL022011000036</td><td>2011</td><td></td><td>77777777760731</td><td>T2</td><td>T</td><td>52</td><td>223325153326</td><td>1041</td><td>Ghana</td><td>F/61</td></tr><tr><td>CHL022013000082</td><td>2013</td><td></td><td>77777777760731</td><td>T2</td><td>T</td><td>52</td><td>143326153321</td><td>Or13</td><td>unknown</td><td>F/67</td></tr><tr><td>CHL022011000028</td><td>2011</td><td></td><td>77777777760771</td><td>T1</td><td>T</td><td>53</td><td>223325153324</td><td>125</td><td>S</td><td>F/30</td></tr><tr><td>CHL022009000004</td><td>2009</td><td></td><td>77777777760771</td><td>T1</td><td>T</td><td>53</td><td>223425153322</td><td>157</td><td>Ugandan</td><td>F/74</td></tr><tr><td>CHL022013100088</td><td>2013</td><td></td><td>77777777760771</td><td>T1</td><td>T</td><td>53</td><td>223327153321</td><td>406</td><td>H</td><td>M/53</td></tr><tr><td>CHL022012300066</td><td>2012</td><td></td><td>77777777760771</td><td>T1</td><td>T</td><td>53</td><td>222325143325</td><td>1185</td><td>Ghana</td><td>M/63</td></tr><tr><td>CHL022010000006</td><td>2010</td><td></td><td>77777777760771</td><td>T1</td><td>T</td><td>53</td><td>223325154324</td><td>1562</td><td>S</td><td>M/66</td></tr><tr><td>CHL022012000057</td><td>2012</td><td></td><td>77777777760771</td><td>T1</td><td>T</td><td>53</td><td>223124153325</td><td>1674</td><td>H</td><td>M/28</td></tr><tr><td>CHL022010000016</td><td>2010</td><td></td><td>77777777760771</td><td>T1</td><td>T</td><td>53</td><td>222321153325</td><td>Or14</td><td>H</td><td>F/76</td></tr><tr><td>CHL022013100045</td><td>2013</td><td></td><td>77777777760771</td><td>T1</td><td>T</td><td>53</td><td>224321153325</td><td>Or15</td><td>X</td><td>M/50</td></tr><tr><td>CHL022013100089</td><td>2013</td><td></td><td>77777777760771</td><td>T1</td><td>T</td><td>53</td><td>223126153322</td><td>Or16</td><td>LAM</td><td>M/5</td></tr><tr><td>CHL022013000093</td><td>2013</td><td></td><td>77777557760771</td><td>TS-Madrid2</td><td>T</td><td>58</td><td>223326153224</td><td>782</td><td>H</td><td>F/22</td></tr><tr><td>CHL022009000005</td><td>2009</td><td></td><td>777777607560771</td><td>LAM6</td><td>LAM</td><td>64</td><td>224326152225</td><td>Or17</td><td>LAM</td><td>F/43</td></tr><tr><td>CHL022012000055</td><td>2012</td><td></td><td>77770036777760771</td><td>X3</td><td>X</td><td>91</td><td>22423143322</td><td>Or18</td><td>unknown</td><td>F/38</td></tr><tr><td>CHL022013000076</td><td>2013</td><td></td><td>777777607740771</td><td>LAM9</td><td>LAM</td><td>161</td><td>124327153220</td><td>Or19</td><td>LAM</td><td>M/37</td></tr><tr><td>CHL022011000039</td><td>2011</td><td></td><td>67777777720771</td><td>H3</td><td>H</td><td>180</td><td>223215143224</td><td>Or20</td><td>Ugandan</td><td>F/29</td></tr><tr><td>CHL022012000056</td><td>2012</td><td></td><td>70377740001771</td><td>CAS</td><td>CAS</td><td>203</td><td>225425163523</td><td>855</td><td>Delhi/CAS</td><td>F/38</td></tr><tr><td>CHL022011000029</td><td>2011</td><td></td><td>776137607760771</td><td>LAM3</td><td>LAM</td><td>211</td><td>224326133323</td><td>251</td><td>LAM</td><td>F/83</td></tr><tr><td>CHL022013000102</td><td>2013</td><td></td><td>776137607760771</td><td>LAM3</td><td>LAM</td><td>211</td><td>224326133323</td><td>251</td><td>LAM</td><td>M/49</td></tr><tr><td>CHL022011000026</td><td>2011</td><td></td><td>776137607760771</td><td>LAM3</td><td>LAM</td><td>211</td><td>22426133323</td><td>880</td><td>LAM</td><td>M/74</td></tr><tr><td>CHL022009000002</td><td>2009</td><td></td><td>777740777760771</td><td>T1</td><td>T</td><td>219</td><td>224315153323</td><td>27</td><td>Cameroon</td><td>F/74</td></tr><tr><td>CHL022010000015</td><td>2010</td><td></td><td>77774077560771</td><td>Unknown</td><td>Unknown</td><td>222</td><td>12426153222</td><td>179</td><td>LAM</td><td>M/36</td></tr><tr><td>CHL022013000098</td><td>2013</td><td></td><td>77774077560771</td><td>Unknown</td><td>Unknown</td><td>222</td><td>125326143221</td><td>Or21</td><td>unknown</td><td>M/56</td></tr></tbody></table>
<table><tr><td>CHL022013000099</td><td>2013</td><td>777775777760771</td><td>T1</td><td>T</td><td>281</td><td>223321153324</td><td>585</td><td>H</td><td>M/46</td></tr><tr><td>CHL02201000007</td><td>2010</td><td>777776777760731</td><td>X1</td><td>X</td><td>336</td><td>224325153325</td><td>35</td><td>X</td><td>F/38</td></tr><tr><td>CHL022010000014</td><td>2010</td><td>7777777767760771</td><td>T1</td><td>T</td><td>373</td><td>223425153322</td><td>157</td><td>Ugandan</td><td>M/61</td></tr><tr><td>CHL022013100100</td><td>2013</td><td>77777777620771</td><td>H3</td><td>H</td><td>390</td><td>223321143324</td><td>1671</td><td>H</td><td>F/25</td></tr><tr><td>CHL022011000020</td><td>2011</td><td>777777606360771</td><td>LAM9</td><td>LAM</td><td>452</td><td>124326154226</td><td>Or22</td><td>LAM</td><td>M/31</td></tr><tr><td>CHL022011000031</td><td>2011</td><td>777777607760760</td><td>LAM9</td><td>LAM</td><td>738</td><td>124326153226</td><td>151</td><td>LAM</td><td>F/73</td></tr><tr><td>CHL022010000012</td><td>2010</td><td>77777777520771</td><td>H3</td><td>H</td><td>746</td><td>226325153323</td><td>110</td><td>H</td><td>M/66</td></tr><tr><td>CHL022011000040</td><td>2011</td><td>777777760020611</td><td>H3</td><td>H</td><td>948</td><td>224325161223</td><td>Or23</td><td>unknown</td><td>F/51</td></tr><tr><td>CHL022012100067</td><td>2012</td><td>7777777207760771</td><td>LAM9</td><td>LAM</td><td>1277</td><td>224327145421</td><td>Or24</td><td>unknown</td><td>M/32</td></tr><tr><td>CHL022013000096</td><td>2013</td><td>777777207760771</td><td>LAM9</td><td>LAM</td><td>1277</td><td>224327155421</td><td>Or25</td><td>LAM</td><td>M/39</td></tr><tr><td>CHL022012100064</td><td>2012</td><td>676177607760771</td><td>LAM3</td><td>LAM</td><td>1293</td><td>224316123325</td><td>Or26</td><td>LAM</td><td>F/25</td></tr><tr><td>CHL022012000053</td><td>2012</td><td>777777407560731</td><td>LAM</td><td>LAM</td><td>1355</td><td>124326153222</td><td>201</td><td>LAM</td><td>M/16</td></tr><tr><td>CHL022011000024</td><td>2011</td><td>677777627760771</td><td>T1</td><td>T</td><td>1474</td><td>224326163321</td><td>375</td><td>LAM</td><td>M/71</td></tr><tr><td>CHL022013100087</td><td>2013</td><td>776177607560771</td><td>LAM</td><td>LAM</td><td>1624</td><td>223326133324</td><td>192</td><td>LAM</td><td>F</td></tr><tr><td>CHL022009000011</td><td>2009</td><td>73737777760771</td><td>T1</td><td>T</td><td>1877</td><td>223321143324</td><td>1671</td><td>Ghana</td><td>F/37</td></tr><tr><td>CHL022013000010</td><td>2013</td><td>377777607560771</td><td>LAM6</td><td>LAM</td><td>1914</td><td>124326163227</td><td>1679</td><td>LAM</td><td>M/48</td></tr><tr><td>CHL022013000077</td><td>2013</td><td>777777704020751</td><td>H1</td><td>H</td><td>2030</td><td>225326153321</td><td>1676</td><td>unknown</td><td>M/34</td></tr><tr><td>CHL022012000061</td><td>2012</td><td>777777607720731</td><td>H3</td><td>H</td><td>2273</td><td>124226153226</td><td>Or27</td><td>LAM</td><td>M/24</td></tr><tr><td>CHL022012000054</td><td>2012</td><td>700000004020771</td><td>H1</td><td>H</td><td>2623</td><td>225313153323</td><td>42</td><td>H</td><td>M/46</td></tr><tr><td>CHL022012000063</td><td>2012</td><td>577777607760751</td><td>LAM9</td><td>LAM</td><td>2745</td><td>124324143320</td><td>Or28</td><td>LAM</td><td>M/31</td></tr><tr><td>CHL022011200034</td><td>2011</td><td>777777777600471</td><td>Unknown</td><td>Unknown</td><td>4013*</td><td>224325143223</td><td>377</td><td>unknown</td><td>M/40</td></tr><tr><td>CHL022012000046</td><td>2012</td><td>777777777600471</td><td>Unknown</td><td>Unknown</td><td>4013*</td><td>222326143321</td><td>1672</td><td>unknown</td><td>M/67</td></tr><tr><td>CHL022008000001</td><td>2008</td><td>7633577775560771</td><td>T1</td><td>T</td><td>Or01</td><td>223425133322</td><td>340</td><td>unknown</td><td>F/74</td></tr><tr><td>CHL022013000013</td><td>2010</td><td>732457667560771</td><td>T1</td><td>T</td><td>Or02</td><td>224326133324</td><td>313</td><td>LAM</td><td>F/37</td></tr><tr><td>CHL022011000019</td><td>2011</td><td>763777777560771</td><td>T1</td><td>T</td><td>Or03</td><td>224225143224</td><td>Or29</td><td>Ugandan</td><td>M/63</td></tr><tr><td>CHL022011000022</td><td>2011</td><td>767777667760771</td><td>T1</td><td>T</td><td>Or04</td><td>124326153220</td><td>190</td><td>LAM</td><td>F/33</td></tr><tr><td>CHL022011000035</td><td>2011</td><td>763377347760471</td><td>T4-CEU1</td><td>T</td><td>Or05</td><td>223325173333</td><td>721</td><td>unknown</td><td>F/36</td></tr><tr><td>CHL022013300048</td><td>2013</td><td>67677377777631</td><td>Unknown</td><td>Unknown</td><td>Or06</td><td>222324253322</td><td>5</td><td>Bovis</td><td>F/84</td></tr><tr><td>CHL022012000059</td><td>2012</td><td>776577607760771</td><td>LAM9</td><td>LAM</td><td>Or07</td><td>224326133324</td><td>313</td><td>LAM</td><td>M/19</td></tr><tr><td>CHL022012100062</td><td>2012</td><td>764077767777771</td><td>Unknown</td><td>Unknown</td><td>Or08</td><td>232224253322</td><td>309</td><td>unknown</td><td>F/23</td></tr><tr><td>CHL022012100069</td><td>2012</td><td>77737775760771</td><td>T1</td><td>T</td><td>Or09</td><td>223425143324</td><td>1675</td><td>Ghana</td><td>M/68</td></tr><tr><td>CHL022012000070</td><td>2012</td><td>777377607560771</td><td>LAM6</td><td>LAM</td><td>Or10</td><td>122326153226</td><td>Or30</td><td>LAM</td><td>M/45</td></tr><tr><td>CHL022012100071</td><td>2012</td><td>664077777777771</td><td>Unknown</td><td>Unknown</td><td>Or11</td><td>232224243322</td><td>990</td><td>unknown</td><td>M/63</td></tr><tr><td>CHL022012100072</td><td>2012</td><td>764077777777671</td><td>AFRI</td><td>AFRI</td><td>Or12</td><td>232224243322</td><td>990</td><td>unknown</td><td>M/93</td></tr><tr><td>CHL022013100085</td><td>2013</td><td>77777775120771</td><td>H3</td><td>H</td><td>Or13</td><td>225326153323</td><td>247</td><td>H</td><td>M/57</td></tr><tr><td>CHL022013000090</td><td>2013</td><td>740136777760771</td><td>X1</td><td>X</td><td>Or14</td><td>224426143322</td><td>Or31</td><td>unknown</td><td>M/79</td></tr><tr><td>CHL022013000097</td><td>2013</td><td>777764003560771</td><td>T1</td><td>T</td><td>Or15</td><td>124226153222</td><td>179</td><td>LAM</td><td>M/22</td></tr><tr><td>CHL022013000103</td><td>2013</td><td>664077777747600</td><td>BOV</td><td>BOV</td><td>Or16</td><td>232224263322</td><td>48</td><td>Bovis</td><td>M/62</td></tr></table>

## Supporting Information

S2 Table. Original PDF version available online at:

http://journals.plos.org/plosone/article/asset?unique&id=info:doi/10.1371/journal.pone.0118007.s005

<footnote>(*) Newly created SIT are highlighted by an asterisk (*)</footnote>