# Analysis of Mycobacterium tuberculosis Genotypic Lineage Distribution in Chile and Neighboring Countries

Jaime Lagos, David Couvin, Loredana Arata, Javier Tognarelli, Carolina Aguayo, Tamara Leiva, Fabiola Arias, Juan Carlos Hormazabal, Nalin Rastogi, Jorge Fernández

► To cite this version:

Jaime Lagos, David Couvin, Loredana Arata, Javier Tognarelli, Carolina Aguayo, et al. Analysis of Mycobacterium tuberculosis Genotypic Lineage Distribution in Chile and Neighboring Countries. PLoS ONE, 2016, 11 (8), pp.160434 - 160434. 10.1371/journal.pone.0160434. pasteur-01652105

HAL Id: pasteur-01652105

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

Submitted on 29 Nov 2017

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
<figure><img src="image_2.png" /></figure>

RESEARCH ARTICLE

# Analysis of *Mycobacterium tuberculosis*
Genotypic Lineage Distribution in Chile and
Neighboring Countries

Jaime Lagos¹, David Couvin², Loredana Arata¹, Javier Tognarelli¹, Carolina Aguayo¹, Tamara Leiva³, Fabiola Arias³, Juan Carlos Hormazabal⁴, Nalin Rastogi²*, Jorge Fernández¹*

1. Subdepartment of Molecular Genetics, Public Health Institute of Chile, Santiago, Chile, 2. WHO Supranational TB Reference Laboratory, TB and Mycobacteria Unit, Institut Pasteur de la Guadeloupe, Guadeloupe, France, 3. Mycobacteria Laboratory, Public Health Institute of Chile, Santiago, Chile, 4. Subdepartment of Infectious Disease, Public Health Institute of Chile, Santiago, Chile

* jfernand@ispch.cl (JF); nrastogi@pasteur-guadeloupe.fr (NR)

## Abstract

Tuberculosis (TB), caused by the pathogen *Mycobacterium tuberculosis* (MTB), remains a disease of high importance to global public health. Studies into the population structure of MTB have become vital to monitoring possible outbreaks and also to develop strategies regarding disease control. Although Chile has a low incidence of MTB, the current rates of migration have the potential to change this scenario. We collected and analyzed a total of 458 *M. tuberculosis* isolates (1 isolate per patient) originating from all 15 regions of Chile. The isolates were genotyped using the spoligotyping method and the data obtained were analyzed and compared with the SITVIT2 database. A total of 169 different patterns were identified, of which, 119 patterns (408 strains) corresponded to Spoligotype International Types (SITs) and 50 patterns corresponded to orphan strains. The most abundantly represented SITs/lineages were: SIT53/T1 (11.57%), SIT33/LAM3 (9.6%), SIT42/LAM9 (9.39%), SIT50/H3 (5.9%), SIT37/T3 (5%); analysis of the spoligotyping minimum spanning tree as well as spoligoforest were suggestive of a recent expansion of SIT42, SIT50 and SIT37; all of which potentially evolved from SIT53. The most abundantly represented lineages were LAM (40.6%), T (34.1%) and Haarlem (13.5%). LAM was more prevalent in the Santiago (43.6%) and Concepción (44.1%) isolates, rather than the Iquique (29.4%) strains. The proportion of X lineage was appreciably higher in Iquique and Concepción (11.7% in both) as compared to Santiago (1.6%). Global analysis of MTB lineage distribution in Chile versus neighboring countries showed that evolutionary recent lineages (LAM, T and Haarlem) accounted together for 88.2% of isolates in Chile, a pattern which mirrored MTB lineage distribution in neighboring countries (n = 7378 isolates recorded in SITVIT2 database for Peru, Brazil, Paraguay, and Argentina; and published studies), highlighting epidemiological advantage of Euro-American lineages in this region. Finally, we also observed exclusive emergence of patterns SIT4014/X1 and SIT4015 (unknown lineage signature) that have hitherto been found exclusively in Chile, indicating that conditions specific to Chile, along

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

OPEN ACCESS

**Citation:** Lagos J, Couvin D, Arata L, Tognarelli J, Aguayo C, Leiva T, et al. (2016) Analysis of *Mycobacterium tuberculosis* Genotypic Lineage Distribution in Chile and Neighboring Countries. PLoS ONE 11(8): e0160434. doi:10.1371/journal.pone.0160434

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

**Received:** April 11, 2016

**Accepted:** July 19, 2016

**Published:** August 12, 2016

**Copyright:** © 2016 Lagos 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:** The author(s) received no specific funding for this work.

**Competing Interests:** The authors have declared that no competing interests exist.
with the unique genetic makeup of the Chilean population, might have allowed for a possible co-evolution leading to the success of these emerging genotypes.

Introduction

*Mycobacterium tuberculosis* (MTB), a Gram-positive bacterium, is the causative agent for tuberculosis (TB); a disease that presents millions of new cases and deaths every year and is considered as a major threat to public health [1]. Although a quick and assertive diagnosis is essential for early treatment and disease management, many people lack proper access to the facilities required for adequate or early diagnosis. In this context, the molecular typing of *MTB* has greatly improved our knowledge and control of tuberculosis by allowing for: (i) the detection of unsuspected transmission, (ii) the identification of false-positive cultures, (iii) discrimination between reinfection and relapse [2], and (iv) the study and analysis of global biodiversity and phylogeographical variations of the tubercle bacilli [3]. In the last 2 decades, the cumbersome “gold standard” IS6110-RFLP methodology necessitating Southern blotting [4,5] was progressively replaced by PCR-based methods, namely spoligotyping and Mycobacterial interspersed repetitive unit-variable number of tandem-repeat (MIRU-VNTR) typing [6–10]. These methods showed satisfactory discriminatory power and reproducibility [9], particularly when used concomitantly that further allows to identify MTB genetic lineages—a prerequisite to study TB population structure at local, regional, and global scale [11].

As far as TB in the Americas is concerned, it accounted for approximately 285,200 new cases of TB in 2013, which equals to 29 new cases per 100,000 inhabitants [12]. Despite two-thirds (69%) of all cases reported being from South American patients, wide epidemiological differences including varying incidence rates were found both country and region wise, highest incidence rates being from Brazil (46/100,000), Bolivia (123/100,000) and Peru (124/100,000) [12]. In Chile different rates of incidence ranging from 7.5/100,000 to 23.7/100,000, were observed, the highest rates being reported from the northern regions of the country [13]. According to the latest available data, the TB incidence in Chile for 2014 was 12.3 per 100,000, which remains significantly higher than the rate of elimination recommended for declassifying TB as a public health problem (<5 per 100,000).

Investigations conducted on MTB genotypic lineage distribution in Southern America have shown that although the Euro-American lineage is the most widely represented, regional differences in the distribution of lineages/sublineages are frequently observed both between as well as within countries [14–22]. According to the SITVIT2, the LAM, Haarlem and T families are the most commonly observed members of the Euro-American lineage in South and Central America and the Caribbean, a distribution profile shared with Europe and Middle Africa [23]. Nonetheless, some regional specificities such as cases involving a LAM family strain designated as RD<sup>Rio</sup> associated with multiple drug resistance (MDR) in Brazil [18], and the almost exclusive “significant” presence of the Beijing family in Peru and Colombia associated both with MDR as well as XDR cases [20,21] may be cited. It is today possible to pinpoint such specificities as well underline the finer differences in MTB population structure at local, regional and macro regional levels, thanks to international genotyping databases and web tools for molecular epidemiology of tuberculosis [23, 24].

Limited MTB genotyping studies in Chile, conducted primarily in metropolitan regions, suggested that Chile is dominated by three major MTB lineages, LAM, T, and Haarlem [25–27]. To provide with a more comprehensive snapshot of MTB lineages circulating in Chile, we
collected and analyzed a total of 458 *M. tuberculosis* isolates (1 isolate per patient) originating from all 15 regions of Chile, including the Metropolitan Region and large cities in north (e.g., Iquique) as well as the south (e.g., Concepción). The data on MTB lineage distribution in Chile versus neighboring countries were compared using the SITVIT2 database [3, 28], so as to highlight exclusive emergence of certain genotypes as well as the potential associated risks in Chile.

# Materials and Methods

## Sample collection

A total of 458 isolates (1 isolate per patient) of *M. tuberculosis* were collected from 15 regions of Chile (S1 Fig). The isolates, obtained from treatment virgin patients (VT), were collected between 2011 (n = 338 or 73.80%) and 2012 (n = 120 or 26.20%), all the isolates were used in the study. Of them, 304 isolates were from male patients (66.38%) and 154 from female patients (33.62%). The isolates were extracted from body fluids (sputum: 95%, bronchoalveolar lavage fluid: 3.7%, tissue: 0.9%, blood 0.2% and pleural fluid: 0.2%). TB was diagnosed in the clinical centers of the respective cities and the diagnosis was confirmed at the Supranational Reference Mycobacteria Laboratory at the Institute of Public Health of Chile which also conducted first and second line drugs susceptibility tests using the Lowenstein-Jensen medium and the proportion method [29].

The strains used in this study were collected routinely during activities of the state TB control program. No patients were contacted to request additional information. The study was reviewed and approved by the review board of Biomedical Department of Public Health Institute who granted permission for use of the MTB isolates and clinical data for the purposes of the study and waived the need for written informed consent from participants. This study is part of the monitoring carried out by the Institute of Public Health of Chile for the diagnosis and characterization of infectious agents.

## DNA extraction and genotyping

MTB isolates obtained from the Lowenstein Jensen medium cultures were resuspended in 500 µL TE (10 mM Tris pH 8.0, 1 mM EDTA pH 8.0). Suspensions were inactivated by heating at 95°C for 15 minutes. Bacterial DNA was isolated by treatment with cetyltrimethylammonium bromide (CTAB) in the presence of 0.7 M NaCl as described previously [30].

Spoligotyping was performed at the Molecular Genetics sub-department of Institute of Public Health of Chile by using the Spoligotyping kit (Ocimum Biosolutions) as described by Kamerbeek *et al.* [8]. Briefly, the 43 spacers between the direct repeats in the target region were amplified using biotinylated primers and the PCR products were then hybridized to a membrane and visualized by chemiluminescence. *M. tuberculosis* H37Rv and *M. bovis* BCG controls, included in the kit, were used as controls for each run. Spoligotyping results were converted into the octal code for comparison with the SITVIT2 proprietary database of Institut Pasteur de la Guadeloupe which is an updated version of the previously released SITVITWEB database (available at: http://www.pasteur-guadeloupe.fr:8081/SITVIT_ONLINE/) [23]. Within this database, a SIT is created when 2 or more isolates share an identical spoligotyping pattern.

## Phylogenetic and statistical analysis

In order to study the phylogenetic relationships of given spoligotypes obtained in this study (n = 458 isolates), a minimum spanning tree (MST) was created using the BioNumerics version 6.6 software (Applied Maths, Sint-Martens-Latem, Belgium; available at: http://www.applied-
maths.com/bionumerics). Furthermore, a spoligoforest tree was constructed using the Fruchterman-Reingold algorithm in the SpolTools software (http://www.emi.unsw.edu.au/spolTools) [31]; the tree was reshaped and colored using the GraphViz software (http://www graphviz.org) [32]. Unlike the MST, the spoligoforest is a directed graph which tentatively highlights parent to descendant relationships between spoligotypes. Lastly, Fisher's Exact and Pearson's Chi-squared tests were calculated using R software. A p-value<0.05 was considered statistically significant.

## Results

The present study identified spoligotypes for clinical isolates of *M. tuberculosis* (n = 458) with a global male/female sex ratio of 304/154 = 1.97. The isolates were collected from different cities of Chile between the years 2011 and 2012. These cities were: Santiago (n = 181), Concepción (n = 34), Iquique (n = 34), Valparaiso (n = 29), Rancagua (n = 20), Temuco (n = 18), Los Angeles (n = 16), Osorno (n = 15), Arica (n = 14), Puerto Montt (n = 12), Talca (n = 12), Talcahuano (n = 12), Antofagasta (n = 11), Coquimbo (n = 11), Castro (n = 9), Punta Arenas (n = 7), Ovalle (n = 6), Valdivia (n = 6), Copiapo (n = 5), Coihaique (n = 2), Viña del Mar (n = 2), Vallenar (n = 1), and Calama (n = 1). Detailed demographic and genotyping information is listed in S1 Table. Analysis of drug resistance profiles showed that most isolates were pan-susceptible (n = 402) while some were resistant to any drug (n = 44) or multidrug resistant (n = 6). For a few cases, data regarding resistance was unavailable (n = 6). Three of the six MDR-TB strains belonged to SIT53/T1 and two were from Santiago. Among strains that were resistant to any antibiotic, our results showed that: (i) of the six belonging to SIT33, five were resistant to streptomycin; (ii) the six strains of SIT42 presented with a heterogeneous resistance, one or a mixture of antibiotics, and (iii) three out of a total of six strains of SIT4015 presented with some form of resistance throughout the study.

Spoligotyping of the 458 isolates revealed a total of 169 different patterns: 119 patterns (408 strains, 89.08% of sample population) corresponded to the shared-types or SITs when compared to the SITVIT2 database, and 50 patterns (50 strains, 10.92% of sample population) corresponded to orphan strains that have not yet been reported in literature. A total of 95/119 SITs (n = 357 isolates) matched a preexisting shared-type recorded in the database, whereas 24/119 SITs (n = 51 isolates) were newly created when they matched orphan strains recorded in the database (SIT4014-SIT4037). The 50 orphan remaining strains were also recorded in the database. A detailed description of the 119 SITs (n = 408 isolates) and their corresponding spoligotyping defined lineages/sublineages are shown in Table 1.

For the 458 isolates analyzed in this study, distribution of lineages was as follow: LAM (n = 186 isolates, or 40.61%); T (n = 156 isolates, or 34.06%); Haarlem (n = 62 isolates, or 13.54%); X (n = 18 isolates, or 3.93%); Beijing (n = 2 isolates, or 0.44%); AFRI (n = 1 isolate, or 0.22%); Cameroon (n = 1 isolate, or 0.22%); and Unknown (n = 32 isolates, or 6.99%) (S1 Table). We observed significant differences in prevalence when the proportion of LAM, T, Haarlem and X lineages were compared within three of Chile's most prominent cities, viz., Santiago, Concepción and Iquique (p-value = 0.0071). The observed proportion of Haarlem was remarkably lower in Concepción city (2.94%) as compared to Santiago (16.57%) and Iquique (11.76%). Regarding other lineages: (i) Beijing lineage was represented by about 3% of strains in Iquique and Santiago city; (ii) LAM lineage represented around 44% of strains in Santiago and Concepción, but only 29% in Iquique city; (iii) X lineage was more noticeable in Iquique and Concepción, representing around 12% of isolates, but in Santiago its incidence was less than 2%; (iv) the unknown strain was conspicuous in Iquique (11.8%) and Concepción (14.7%) but comparatively less prevalent in Santiago (6%) (Fig 1A).
<table><caption>Table 1. Description of 119 shared-types (SITs; n = 408 isolates) and corresponding spoligotyping defined lineages/sublineages starting from a total of 458 *M. tuberculosis* strains isolated in various cities in 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 (0.44)</td><td>0.02</td><td></td><td>Beijing</td><td>Clustered</td></tr><tr><td>2</td><td>000000004020771</td><td>1 (0.22)</td><td>0.21</td><td></td><td>H2</td><td>Unique</td></tr><tr><td>4</td><td>000000007760771</td><td>1 (0.22)</td><td>0.26</td><td></td><td>Unknown</td><td>Unique</td></tr><tr><td>20</td><td>677777607760771</td><td>3 (0.66)</td><td>0.33</td><td></td><td>LAM1</td><td>Clustered</td></tr><tr><td>33</td><td>776177607760771</td><td>44 (9.61)</td><td>3.25</td><td></td><td>LAM3</td><td>Clustered</td></tr><tr><td>37</td><td>77773777760771</td><td>23 (5.02)</td><td>4.23</td><td></td><td>T3</td><td>Clustered</td></tr><tr><td>39</td><td>777777347760471</td><td>1 (0.22)</td><td>0.67</td><td></td><td>T4-CEU1</td><td>Unique</td></tr><tr><td>42</td><td>777777607760771</td><td>43 (9.39)</td><td>1.21</td><td></td><td>LAM9</td><td>Clustered</td></tr><tr><td>44</td><td>777777757760771</td><td>1 (0.22)</td><td>0.46</td><td></td><td>T5</td><td>Unique</td></tr><tr><td>47</td><td>77777774020771</td><td>5 (1.09)</td><td>0.3</td><td></td><td>H1</td><td>Clustered</td></tr><tr><td>50</td><td>77777777720771</td><td>27 (5.9)</td><td>0.68</td><td></td><td>H3</td><td>Clustered</td></tr><tr><td>51</td><td>77777777760700</td><td>3 (0.66)</td><td>0.95</td><td></td><td>T1</td><td>Clustered</td></tr><tr><td>52</td><td>77777777760731</td><td>2 (0.44)</td><td>0.21</td><td></td><td>T2</td><td>Clustered</td></tr><tr><td>53</td><td>77777777760771</td><td>53 (11.57)</td><td>0.8</td><td></td><td>T1</td><td>Clustered</td></tr><tr><td>58</td><td>77777557760771</td><td>4 (0.87)</td><td>2.05</td><td></td><td>T5-Madrid2</td><td>Clustered</td></tr><tr><td>60</td><td>777777607760731</td><td>5 (1.09)</td><td>1.12</td><td></td><td>LAM4</td><td>Clustered</td></tr><tr><td>61</td><td>77777743760771</td><td>1 (0.22)</td><td>0.1</td><td></td><td>Cameroon</td><td>Unique</td></tr><tr><td>62</td><td>77777774020771</td><td>1 (0.22)</td><td>0.17</td><td></td><td>H1</td><td>Unique</td></tr><tr><td>64</td><td>777777607560771</td><td>8 (1.75)</td><td>1.94</td><td></td><td>LAM6</td><td>Clustered</td></tr><tr><td>78</td><td>77777777760771</td><td>1 (0.22)</td><td>1.45</td><td></td><td>T</td><td>Unique</td></tr><tr><td>81</td><td>77737607760771</td><td>1 (0.22)</td><td>2.63</td><td></td><td>LAM9</td><td>Unique</td></tr><tr><td>91</td><td>700036777760771</td><td>5 (1.09)</td><td>1.36</td><td></td><td>X3</td><td>Clustered</td></tr><tr><td>93</td><td>77737607760771</td><td>1 (0.22)</td><td>0.22</td><td></td><td>LAM5</td><td>Unique</td></tr><tr><td>99</td><td>75777777720771</td><td>2 (0.44)</td><td>2.41</td><td></td><td>H3</td><td>Clustered</td></tr><tr><td>111</td><td>776167607760771</td><td>1 (0.22)</td><td>3.13</td><td></td><td>LAM3</td><td>Unique</td></tr><tr><td>118</td><td>77776777760771</td><td>1 (0.22)</td><td>0.6</td><td></td><td>T1</td><td>Unique</td></tr><tr><td>119</td><td>77776777760771</td><td>1 (0.22)</td><td>0.09</td><td></td><td>X1</td><td>Unique</td></tr><tr><td>124</td><td>777777777700771</td><td>1 (0.22)</td><td>1.47</td><td></td><td>Unknown</td><td>Unique</td></tr><tr><td>130</td><td>776177607760731</td><td>2 (0.44)</td><td>1.55</td><td></td><td>LAM3</td><td>Clustered</td></tr><tr><td>137</td><td>777776777760601</td><td>3 (0.66)</td><td>0.29</td><td></td><td>X2</td><td>Clustered</td></tr><tr><td>154</td><td>75777777760771</td><td>1 (0.22)</td><td>0.88</td><td></td><td>T1</td><td>Unique</td></tr><tr><td>162</td><td>777777607760671</td><td>1 (0.22)</td><td>2.86</td><td></td><td>LAM9</td><td>Unique</td></tr><tr><td>163</td><td>777777607760700</td><td>1 (0.22)</td><td>2.5</td><td></td><td>LAM4</td><td>Unique</td></tr><tr><td>180</td><td>67777777720771</td><td>2 (0.44)</td><td>3.23</td><td></td><td>H3</td><td>Clustered</td></tr><tr><td>194</td><td>677737607760731</td><td>2 (0.44)</td><td>6.06</td><td></td><td>LAM2</td><td>Clustered</td></tr><tr><td>195</td><td>677777607660771</td><td>1 (0.22)</td><td>6.25</td><td></td><td>LAM1</td><td>Unique</td></tr><tr><td>211</td><td>776137607760771</td><td>12 (2.62)</td><td>10.17</td><td></td><td>LAM3</td><td>Clustered</td></tr><tr><td>222</td><td>777774077560771</td><td>4 (0.87)</td><td>3.96</td><td></td><td>Unknown</td><td>Clustered</td></tr><tr><td>238</td><td>777777777700171</td><td>1 (0.22)</td><td>5.56</td><td></td><td>Unknown</td><td>Unique</td></tr><tr><td>239</td><td>77777777760031</td><td>1 (0.22)</td><td>1.64</td><td></td><td>T2</td><td>Unique</td></tr><tr><td>241</td><td>77777777760411</td><td>1 (0.22)</td><td>3.57</td><td></td><td>T1</td><td>Unique</td></tr><tr><td>244</td><td>77777777760601</td><td>1 (0.22)</td><td>0.84</td><td></td><td>T1</td><td>Unique</td></tr><tr><td>245</td><td>77777777760671</td><td>1 (0.22)</td><td>4.76</td><td></td><td>T1</td><td>Unique</td></tr><tr><td>283</td><td>777777704020771</td><td>6 (1.31)</td><td>8.11</td><td></td><td>H1</td><td>Clustered</td></tr><tr><td>393</td><td>77757777760771</td><td>1 (0.22)</td><td>2.44</td><td></td><td>T1</td><td>Unique</td></tr><tr><td>430</td><td>704003347760471</td><td>1 (0.22)</td><td>2.7</td><td></td><td>T4-CEU1</td><td>Unique</td></tr><tr><td>434</td><td>76077777740371</td><td>1 (0.22)</td><td>16.67</td><td></td><td>T1</td><td>Unique</td></tr><tr><td>435</td><td>763777607760771</td><td>1 (0.22)</td><td>16.67</td><td></td><td>LAM9</td><td>Unique</td></tr><tr><td>450</td><td>777776770000000</td><td>1 (0.22)</td><td>0.8</td><td></td><td>Unknown</td><td>Unique</td></tr><tr><td>504</td><td>77737737760771</td><td>1 (0.22)</td><td>4.35</td><td></td><td>T3</td><td>Unique</td></tr><tr><td>720</td><td>776177607760371</td><td>6 (1.31)</td><td>54.55</td><td></td><td>LAM3</td><td>Clustered</td></tr><tr><td>753</td><td>477777607760771</td><td>1 (0.22)</td><td>4.55</td><td></td><td>LAM1</td><td>Unique</td></tr></tbody></table>

(Continued)
<table><caption>**Table 1. (Continued)**</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>787</td><td></td><td>77777777760071</td><td>1 (0.22)</td><td>12.5</td><td>T1</td><td>Unique</td></tr><tr><td>798</td><td></td><td>437777777760771</td><td>2 (0.44)</td><td>33.33</td><td>T1</td><td>Clustered</td></tr><tr><td>880</td><td></td><td>777777000760771</td><td>1 (0.22)</td><td>10</td><td>T1</td><td>Unique</td></tr><tr><td>900</td><td></td><td>677777607740011</td><td>1 (0.22)</td><td>25</td><td>LAM</td><td>Unique</td></tr><tr><td>914</td><td></td><td>77637777720771</td><td>1 (0.22)</td><td>1.45</td><td>Unknown</td><td>Unique</td></tr><tr><td>948</td><td></td><td>777777760020611</td><td>2 (0.44)</td><td>14.29</td><td>H3</td><td>Clustered</td></tr><tr><td>1073</td><td></td><td>777637777760771</td><td>1 (0.22)</td><td>12.5</td><td>T1</td><td>Unique</td></tr><tr><td>1105</td><td></td><td>777773777760771</td><td>1 (0.22)</td><td>2.78</td><td>T1</td><td>Unique</td></tr><tr><td>1116</td><td></td><td>607777777720771</td><td>1 (0.22)</td><td>33.33</td><td>H3</td><td>Unique</td></tr><tr><td>1220</td><td></td><td>754377607760700</td><td>1 (0.22)</td><td>16.67</td><td>LAM</td><td>Unique</td></tr><tr><td>1232</td><td></td><td>777774077560731</td><td>1 (0.22)</td><td>10</td><td>T2</td><td>Unique</td></tr><tr><td>1259</td><td></td><td>777740000160771</td><td>1 (0.22)</td><td>25</td><td>T1</td><td>Unique</td></tr><tr><td>1277</td><td></td><td>777777207760771</td><td>6 (1.31)</td><td>16.67</td><td>LAM9</td><td>Clustered</td></tr><tr><td>1476</td><td></td><td>736177607700171</td><td>1 (0.22)</td><td>3.13</td><td>AFRI_2</td><td>Unique</td></tr><tr><td>1560</td><td></td><td>703777777760771</td><td>2 (0.44)</td><td>22.22</td><td>T1</td><td>Clustered</td></tr><tr><td>1626</td><td></td><td>777777776760771</td><td>1 (0.22)</td><td>7.69</td><td>T1</td><td>Unique</td></tr><tr><td>1685</td><td></td><td>776077607760771</td><td>3 (0.66)</td><td>37.5</td><td>LAM3</td><td>Clustered</td></tr><tr><td>1706</td><td></td><td>771777607760771</td><td>1 (0.22)</td><td>10</td><td>LAM9</td><td>Unique</td></tr><tr><td>1710</td><td></td><td>777737607760731</td><td>4 (0.87)</td><td>36.36</td><td>LAM</td><td>Clustered</td></tr><tr><td>1841</td><td></td><td>776177607740771</td><td>1 (0.22)</td><td>16.67</td><td>LAM3</td><td>Unique</td></tr><tr><td>1890</td><td></td><td>777777707760731</td><td>1 (0.22)</td><td>8.33</td><td>T2</td><td>Unique</td></tr><tr><td>2014</td><td></td><td>776177607360771</td><td>1 (0.22)</td><td>20</td><td>LAM3</td><td>Unique</td></tr><tr><td>2015</td><td></td><td>736177607760771</td><td>2 (0.44)</td><td>33.33</td><td>LAM3</td><td>Clustered</td></tr><tr><td>2028</td><td></td><td>777777604020771</td><td>1 (0.22)</td><td>10</td><td>Unknown</td><td>Unique</td></tr><tr><td>2045</td><td></td><td>776037607760771</td><td>1 (0.22)</td><td>6.25</td><td>LAM3</td><td>Unique</td></tr><tr><td>2054</td><td></td><td>677777607760371</td><td>1 (0.22)</td><td>8.33</td><td>LAM</td><td>Unique</td></tr><tr><td>2070</td><td></td><td>777777637760771</td><td>4 (0.87)</td><td>23.53</td><td>T1</td><td>Clustered</td></tr><tr><td>2353</td><td></td><td>776177600000000</td><td>2 (0.44)</td><td>18.18</td><td>Unknown</td><td>Clustered</td></tr><tr><td>2366</td><td></td><td>777777777760171</td><td>1 (0.22)</td><td>16.67</td><td>T1</td><td>Unique</td></tr><tr><td>2402</td><td></td><td>776177607760711</td><td>1 (0.22)</td><td>25</td><td>LAM3</td><td>Unique</td></tr><tr><td>2512</td><td></td><td>657757607760771</td><td>1 (0.22)</td><td>16.67</td><td>LAM1</td><td>Unique</td></tr><tr><td>2745</td><td></td><td>577777607760751</td><td>2 (0.44)</td><td>40</td><td>LAM9</td><td>Clustered</td></tr><tr><td>2885</td><td></td><td>777777477720771</td><td>1 (0.22)</td><td>6.25</td><td>H3</td><td>Unique</td></tr><tr><td>2911</td><td></td><td>777737777700771</td><td>1 (0.22)</td><td>33.33</td><td>Unknown</td><td>Unique</td></tr><tr><td>2941</td><td></td><td>074377607760700</td><td>1 (0.22)</td><td>25</td><td>LAM</td><td>Unique</td></tr><tr><td>3056</td><td></td><td>777777777560731</td><td>1 (0.22)</td><td>16.67</td><td>T2</td><td>Unique</td></tr><tr><td>3127</td><td></td><td>000000017760771</td><td>1 (0.22)</td><td>20</td><td>Unknown</td><td>Unique</td></tr><tr><td>3251</td><td></td><td>776177777720771</td><td>1 (0.22)</td><td>25</td><td>H3</td><td>Unique</td></tr><tr><td>3472</td><td></td><td>776177607760770</td><td>1 (0.22)</td><td>33.33</td><td>LAM3</td><td>Unique</td></tr><tr><td>3491</td><td></td><td>776177607760671</td><td>1 (0.22)</td><td>25</td><td>LAM3</td><td>Unique</td></tr><tr><td>3934</td><td></td><td>437777777760731</td><td>1 (0.22)</td><td>20</td><td>T2</td><td>Unique</td></tr><tr><td>3965</td><td></td><td>707774077560771</td><td>1 (0.22)</td><td>7.69</td><td>T1</td><td>Unique</td></tr><tr><td>3967</td><td></td><td>741774077560771</td><td>1 (0.22)</td><td>33.33</td><td>T1</td><td>Unique</td></tr><tr><td>4014*</td><td></td><td>703776777760771</td><td>7 (1.53)</td><td>100</td><td>X1</td><td>Clustered</td></tr><tr><td>4015*</td><td></td><td>776177703770371</td><td>7 (1.53)</td><td>100</td><td>Unknown</td><td>Clustered</td></tr><tr><td>4016*</td><td></td><td>777777774100771</td><td>2 (0.44)</td><td>100</td><td>Unknown</td><td>Clustered</td></tr><tr><td>4017*</td><td></td><td>777737637760771</td><td>3 (0.66)</td><td>75</td><td>T3</td><td>Clustered</td></tr><tr><td>4018*</td><td></td><td>777377636000771</td><td>2 (0.44)</td><td>100</td><td>Unknown</td><td>Clustered</td></tr><tr><td>4019*</td><td></td><td>776177607720771</td><td>1 (0.22)</td><td>50</td><td>H3</td><td>Unique</td></tr><tr><td>4020*</td><td></td><td>777347777760671</td><td>1 (0.22)</td><td>50</td><td>T1</td><td>Unique</td></tr><tr><td>4021*</td><td></td><td>77737777720671</td><td>1 (0.22)</td><td>50</td><td>H3</td><td>Unique</td></tr><tr><td>4022*</td><td></td><td>777777364020771</td><td>1 (0.22)</td><td>50</td><td>H1</td><td>Unique</td></tr><tr><td>4023*</td><td></td><td>777774077560711</td><td>1 (0.22)</td><td>50</td><td>T1</td><td>Unique</td></tr></tbody></table>

(Continued)
<figure><img src="image_4.png" /><figcaption>**Fig 1. Phylogeographical distribution of *M. tuberculosis* lineages identified in our study and MST illustrating evolutionary relationships.** (A) Phylogeographical distribution of major *M. tuberculosis* lineages in the 3 most important cities in our study (Santiago, Iquique, Concepción) as well as following neighboring countries: Peru, Brazil, and Paraguay (referring to the SITVIT2 database). (B) A minimum spanning tree (MST) illustrating evolutionary relationships between spoligotypes in various cities of Chile (n = 458 isolates). 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 color of the circles indicates the phylogenetic lineage to which the specific pattern belongs. The labels of nodes indicate predominant SITs in study (containing at least 5 or more isolates).</figcaption></figure>

doi:10.1371/journal.pone.0160434.g001

When the proportion of predominant SITs found in Chile (n = 458 isolates) as well as the distribution patterns for major cities (Santiago, Iquique and Concepción) were compared, it was observed that the strains belonging to SIT91/X3 were localized to Iquique, whereas the strains belonging to SIT37/T3, SIT60/LAM4, SIT64/LAM6, SIT283/H1 and SIT720/LAM3 were mostly isolated from the Santiago based population, and strains belonging to newly created SIT4014/X1 were predominantly found in Concepción city (p-value<0.008; Table 3). Significant differences were noted when the distribution patterns of predominant SITs in this study were compared with that of neighboring countries (Peru, Brazil, Paraguay, Argentina) as recorded in the SITVIT2 database (n = 7378 isolates; p-value<0.00001; Table 3 and Fig 1A). In comparison to neighboring countries, SIT33 is well represented in Chile (9.61% of isolates)
<table><caption>Table 2. Description of clusters containing >1% (n = 5 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>44 (9.61)</td><td>3.25</td><td>AMER-S 37.08, AFRI-S 24.0, AMER-N 11.74, EURO-S 10.49, EURO-W 6.06, AMER-C 3.99</td><td>ZAF 24.0, PER 16.4, USA 11.74, BRA 10.56, ESP 6.57, ARG 4.21, CHL 4.06, FXX 3.91, ITA 3.4, HND 3.18</td></tr><tr><td>37 (T3) 777737777760771</td><td>23 (5.02)</td><td>4.23</td><td>AFRI-E 18.38, AMER-S 13.05, EURO-N 11.4, EURO-W 11.03, ASIA-W 10.11, AMER-N 8.27, ASIA-E 5.88, EURO-S 5.52, EURO-E 4.23, ASIA-S 3.49, AFRI-S 3.13</td><td>ETH 15.99, USA 7.35, CHL 5.52, SWE 4.96, SAU 4.78, CHN 4.78, FXX 4.04, ITA 3.86, BRA 3.49, ZAF 3.13, DNK 3.13</td></tr><tr><td>42 (LAM9) 777777607760771</td><td>43 (9.39)</td><td>1.21</td><td>AMER-S 31.99, AMER-N 11.19, EURO-S 10.68, EURO-W 8.94, AFRI-N 8.09, EURO-N 4.62, CARI 4.0, AMER-C 3.35, AFRI-E 3.35</td><td>BRA 12.91, USA 11.19, COL 7.19, MAR 6.65, ITA 6.17, FXX 4.76, PER 3.5, ESP 3.16, VEN 3.13</td></tr><tr><td>47 (H1) 777777774020771</td><td>5 (1.09)</td><td>0.3</td><td>EURO-W 19.03, AMER-S 17.06, AMER-N 14.54, EURO-S 12.63, EURO-N 9.81, EURO-E 6.76, ASIA-W 3.77, AFRI-N 3.41</td><td>USA 14.24, ITA 7.72, PER 7.54, AUT 7.54, BRA 7.12, FXX 6.22, FIN 5.63, CZE 3.53, ESP 3.35, SWE 3.17</td></tr><tr><td>50 (H3) 777777777720771</td><td>27 (5.9)</td><td>0.68</td><td>AMER-S 26.58, EURO-W 14.88, AMER-N 14.88, EURO-S 9.78, CARI 4.93, EURO-E 4.68, EURO-N 4.63, AFRI-N 3.6, AFRI-S 3.43, AFRI-M 3.2</td><td>USA 14.85, PER 13.55, BRA 7.1, FXX 5.83, AUT 5.15, ITA 4.6, ESP 4.6, ZAF 3.43, CMR 3.15, CZE 3.1</td></tr><tr><td>53 (T1) 7777777777760771</td><td>53 (11.57)</td><td>0.8</td><td>AMER-S 15.03, EURO-W 14.88, AMER-N 12.82, EURO-S 8.95, EURO-N 7.12, ASIA-W 6.95, AFRI-S 4.72, AFRI-E 4.42, ASIA-E 4.06, AFRI-N 3.35, EURO-E 3.1, CARI 3.07, AMER-C 3.07</td><td>USA 12.55, FXX 7.49, BRA 5.57, ITA 5.07, ZAF 4.62, PER 3.71, TUR 3.3, AUT 3.26</td></tr><tr><td>60 (LAM4) 777777607760731</td><td>5 (1.09)</td><td>1.12</td><td>AFRI-S 36.63, AMER-S 24.49, EURO-W 6.97, AFRI-N 6.74, AFRI-W 6.07, EURO-S 5.39, AMER-N 4.27</td><td>ZAF 36.63, BRA 14.16, FXX 4.94, USA 4.27, MAR 4.27, GMB 3.82, ITA 3.6, VEN 3.37, PER 3.15</td></tr><tr><td>64 (LAM6) 777777607560771</td><td>8 (1.75)</td><td>1.94</td><td>AMER-S 54.0, AMER-N 22.76, EURO-W 5.81, EURO-S 4.36</td><td>BRA 40.44, USA 22.76, GUF 5.81, PRT 3.39</td></tr><tr><td>91 (X3) 700036777760771</td><td>5 (1.09)</td><td>1.36</td><td>AMER-S 37.88, AMER-N 34.33, CARI 17.98, EURO-S 4.36, EURO-N 3.0</td><td>USA 32.15, PER 27.79, HTI 16.35, ESP 4.36, COL 3.27</td></tr><tr><td>211 (LAM3) 776137607760771</td><td>12 (2.62)</td><td>10.17</td><td>AMER-S 27.12, AMER-N 27.12, AMER-C 18.64, EURO-W 17.8, EURO-S 6.78</td><td>USA 27.12, MEX 18.64, FXX 16.95, CHL 12.71, BRA 10.17, ESP 4.24</td></tr><tr><td>283 (H1) 777777704020771</td><td>6 (1.31)</td><td>8.11</td><td>EURO-N 25.68, EURO-W 20.27, AMER-S 18.92, ASIA-S 10.81, AMER-N 8.11, EURO-S 5.41, EURO-E 4.05, AUST 4.05</td><td>LVA 18.92, IND 9.46, AUT 9.46, USA 8.11, CHL 8.11, PER 6.76, FXX 6.76, SWE 5.41, NLD 4.05, CZE 4.05, AUS 4.05</td></tr><tr><td>720 (LAM3) 776177607760371</td><td>6 (1.31)</td><td>54.55</td><td>AMER-S 100.0</td><td>CHL 54.55, ARG 36.36, PER 9.09</td></tr><tr><td>1277 (LAM9) 777777207760771</td><td>6 (1.31)</td><td>16.67</td><td>AMER-S 47.22, CARI 22.22, AMER-N 8.33, AMER-C 8.33, EURO-S 5.56</td><td>DOM 22.22, CHL 22.22, USA 8.33, MEX 8.33, BRA 8.33, PER 5.56, ESP 5.56, ARG 5.56</td></tr><tr><td>4014 (X1) 703776777760771</td><td>7 (1.53)</td><td>100</td><td>AMER-S 100.0</td><td>CHL 100.0</td></tr><tr><td>4015 (Unknown) 776177703770371</td><td>7 (1.53)</td><td>100</td><td>AMER-S 100.0</td><td>CHL 100.0</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.0160434.t002
<table><caption>Table 3. A comparison of the proportion of predominant SITs found in Chile (n = 458 isolates in total, as well as distribution patterns for major cities, i.e., Santiago, Iquique, and Concepción), as compared to neighboring countries (Peru, Brazil, Paraguay, Argentina), recorded in the SITVIT2 database.</caption><thead><tr><th colspan="8">Major M. tuberculosis lineages in the 3 most important cities in our study (Santiago, Iquique, Concepción)</th><th colspan="7">Comparison with neighboring countries available in SITVIT2 database (n = 7378 isolates)</th></tr><tr><th>SIT (Lineage)</th><th>Iquique (n = 34)</th><th>Santiago (n = 181)</th><th>Concepción (n = 34)</th><th>Other cities (n = 299)</th><th>This study (n = 458)</th><th>%</th><th>PERU (n = 1282)</th><th>%</th><th>BRAZIL (n = 4489)</th><th>%</th><th>PARAGUAY (n = 226)</th><th>%</th><th>ARGENTINA (n = 1381)</th><th>%</th></tr></thead><tbody><tr><td>33 (LAM3)</td><td>2</td><td>23</td><td>3</td><td>16</td><td>44</td><td>9,61</td><td>56</td><td>4,37</td><td>132</td><td>2,94</td><td>2</td><td>0,88</td><td>57</td><td>4,13</td></tr><tr><td>37 (T3)</td><td>0</td><td>14</td><td>0</td><td>9</td><td>23</td><td>5,02</td><td>0</td><td>0,00</td><td>10</td><td>0,22</td><td>0</td><td>0,00</td><td>0</td><td>0,00</td></tr><tr><td>42 (LAM9)</td><td>2</td><td>11</td><td>1</td><td>29</td><td>43</td><td>9,39</td><td>73</td><td>5,69</td><td>401</td><td>8,93</td><td>40</td><td>17,70</td><td>63</td><td>4,56</td></tr><tr><td>47 (H1)</td><td>1</td><td>3</td><td>0</td><td>1</td><td>5</td><td>1,09</td><td>40</td><td>3,12</td><td>110</td><td>2,45</td><td>2</td><td>0,88</td><td>8</td><td>0,58</td></tr><tr><td>50 (H3)</td><td>1</td><td>13</td><td>0</td><td>13</td><td>27</td><td>5,90</td><td>146</td><td>11,39</td><td>239</td><td>5,32</td><td>9</td><td>3,98</td><td>29</td><td>2,10</td></tr><tr><td>53 (T1)</td><td>4</td><td>16</td><td>2</td><td>31</td><td>53</td><td>11,57</td><td>171</td><td>13,34</td><td>322</td><td>7,17</td><td>6</td><td>2,65</td><td>64</td><td>4,63</td></tr><tr><td>60 (LAM4)</td><td>0</td><td>3</td><td>1</td><td>1</td><td>5</td><td>1,09</td><td>2</td><td>0,16</td><td>58</td><td>1,29</td><td>2</td><td>0,88</td><td>2</td><td>0,14</td></tr><tr><td>64 (LAM6)</td><td>1</td><td>5</td><td>0</td><td>2</td><td>8</td><td>1,75</td><td>7</td><td>0,55</td><td>138</td><td>3,07</td><td>0</td><td>0,00</td><td>4</td><td>0,29</td></tr><tr><td>91 (X3)</td><td>3</td><td>0</td><td>0</td><td>2</td><td>5</td><td>1,09</td><td>25</td><td>1,95</td><td>4</td><td>0,09</td><td>0</td><td>0,00</td><td>3</td><td>0,22</td></tr><tr><td>211 (LAM3)</td><td>0</td><td>4</td><td>2</td><td>6</td><td>12</td><td>2,62</td><td>1</td><td>0,08</td><td>10</td><td>0,22</td><td>1</td><td>0,44</td><td>1</td><td>0,07</td></tr><tr><td>283 (H1)</td><td>0</td><td>4</td><td>0</td><td>2</td><td>6</td><td>1,31</td><td>1</td><td>0,08</td><td>1</td><td>0,02</td><td>0</td><td>0,00</td><td>2</td><td>0,14</td></tr><tr><td>720 (LAM3)</td><td>0</td><td>5</td><td>0</td><td>1</td><td>6</td><td>1,31</td><td>1</td><td>0,08</td><td>0</td><td>0,00</td><td>0</td><td>0,00</td><td>4</td><td>0,29</td></tr><tr><td>1277 (LAM9)</td><td>1</td><td>3</td><td>0</td><td>2</td><td>6</td><td>1,31</td><td>0</td><td>0,00</td><td>3</td><td>0,07</td><td>0</td><td>0,00</td><td>2</td><td>0,14</td></tr><tr><td>4014 (X1)</td><td>0</td><td>1</td><td>4</td><td>2</td><td>7</td><td>1,53</td><td>0</td><td>0,00</td><td>0</td><td>0,00</td><td>0</td><td>0,00</td><td>0</td><td>0,00</td></tr><tr><td>4015 (Unknown)</td><td>0</td><td>2</td><td>2</td><td>3</td><td>7</td><td>1,53</td><td>0</td><td>0,00</td><td>0</td><td>0,00</td><td>0</td><td>0,00</td><td>0</td><td>0,00</td></tr></tbody></table>

doi:10.1371/journal.pone.0160434.t003
as compared to others countries. Similar results have also been reported in case of SIT211 and other SITs such as SIT37, SIT64, SIT91, SIT283, SIT720, SIT1277, SIT4014 and SIT4015, which are yet to be reported or present in very low percentages. On the other hand, SIT53, SIT50 and SIT42, which are well represented in our study (11.57%, 5.90%, 9.39% respectively) have also been reported with a relative high percentage in the neighboring countries. In general, it can be concluded that last three SITs are widely distributed. All the analysis described is collated in a phylogeographical distribution map (Fig 1A) of major *M. tuberculosis* lineages in the three of the most important Chilean cities enrolled in this study (Santiago, Iquique and Concepción), as well as the neighboring countries (Brazil, Paraguay and Peru; data extracted from the SITVIT2 database).

The Fig 1B shows the spoligotyping-based MST; it is evident that LAM and T lineages com-
pose two of the most predominant groups in our study and include most of the isolates with
frequently observed patterns SIT33, SIT42, and SIT53, SIT37 respectively. Other isolates were
classified as the Haarlem (SIT50 recorded as most frequent) and X (SIT4014 recorded as most
frequent) lineages. More distance was observed amongst isolates that join in with the Haarlem
than those integrating with LAM or T lineages.

The spoligotype analysis of common SITs (SIT42, SIT53, SIT33, SIT50 and SIT37) clearly demonstrated that they were composed of very closely related isolates and were distinguished from one another by changes in only few alleles/spacers. The spoligoforest tree generated by means of the Fruchterman-Reingold algorithm is illustrated in Fig 2. It allows the visualization of all genetic associations and mutations in our strain sample. In this tree, each node represents a spoligotype colored according to the lineage involved, while any potential mutation event that might have occurred from the parental spoligotype is represented as an edge. Thus spoli-goforest analysis not only allowed to identify genetic associations between different MTB strains but also helped to gain insight into their genetic diversity and evolution as follows. Briefly, it confirmed the dominance of SIT42, SIT33 (LAM) and SIT53 (T), and to a lesser extent that of SIT50 (H) and SIT37 (T). The SIT42/LAM9 cluster was the largest node evolved from SIT53/T1 and multiple spoligotypes were observed to be rising from it, e.g., SIT33/LAM3, while the second largest spoligotype pattern SIT50/H3 appeared to have derived originally from SIT53/T1. Lastly, the newly created SIT4014/X1 (found only in this study) was a relative precursor of SIT91/X3 (predominantly found in USA, Peru, and Haiti).

Discussion

In this study, we analyzed the genetic diversity of MTB isolates obtained from clinical samples encompassing 15 regions of Chile. To the best of our knowledge, this is the first study of its kind aimed at investigating the genetic diversity of MTB at the level of the entire country. Our results indicated that the LAM (40.61%), T (34.06%) and Haarlem (13.54%) genotypes were the most abundant spoligotypes of all of *M. tuberculosis* samples collected and analyzed. Together, these three lineage families accounted for 88.21% of all genotypes while all other families, namely Beijing, X, AFRI and Cameroon, accounted for only 4.81%. Interestingly, it was noted that the X lineage is poorly represented in the Santiago population (1.66%) whereas in Iquique and Concepción it accounts for nearly 11.76% each. This difference can potentially be explained as a result of: (i) immigration patterns, which have greatly increased in recent years from countries like Peru, Colombia or Ecuador, or (ii) a genetic predisposition that makes the population of these cities more susceptible to certain genotypes.

This is consistent with data available from most South American countries which claim that the most represented lineages are LAM, T and Haarlem in proportions that are similar to those reported by us; in some cases the proportions may be subject to change but invariably it is
<figure><img src="image_5.png" /><figcaption>**Fig 2.** Spoligoforest tree based on all spoligotypes (n = 458 isolates). Spoligoforest was drawn using the Fruchterman-Reingold algorithm from the SpolTools software (http://www.emi.unsw.edu.au/spolTools)[29], and reshaped and colored using the GraphViz software (http://www graphviz org)[30]. Each spoligotype pattern from the study is represented by a node with area size being proportional to the total number of isolates with that specific pattern. Changes (loss of spacers) are represented by directed edges between nodes, with the arrowheads pointing to descendant spoligotypes. The heuristic used selects a single inbound edge with a maximum weight using a Zipf model. Solid black lines link patterns that are very similar, i.e., loss of one spacer only (maximum weigh being 1.0), while dashed lines represent links of weight comprised between 0.5 and 1, and dotted lines a weight less than 0.5.</figcaption></figure>

doi:10.1371/journal.pone.0160434.g002

these three families that dominate epidemiologically. In Fig 1A, a representation of data avail-
able from three of our neighboring countries (Peru, Brazil and Paraguay) and three of the
major Chilean cities included in this study are plotted. Several published reports on population
structure of MTB in South America have emphatically proven that each family is adequately represented from an epidemiological perspective [14–22]. In Brazil, Vasconcellos and collaborators demonstrated via spoligotyping that LAM is represented in 43.6%, T in 34.9% and Haarlem 18.3% of case studies. Similar results were also found in a study where isolates from patients belonging to 11 states of Brazil were analyzed. The mentioned study demonstrated that the predominant MTB lineage was LAM (46%) followed by T (18.6%), Haarlem (12.2%), X (4.7%), S (1.9%) and, lastly, the East African Indian (EAI) (0.85%) families [14]. However, in Peru and Bolivia it is reported that the most represented family is Haarlem, followed by the LAM and T families [20,33,34]. For Argentina there is comparatively less information in literature but available data clearly show that Haarlem is predominant except for MDR cases where the LAM family dominates [35,36]. In a recent work that was conducted in the border areas of Brazil, Argentina and Paraguay, a higher presence of the LAM family was reported [37]. All of the above mentioned variance can be potentially accounted for by: i) number of isolates analyzed, ii) region of sampling (restricted to a small region or spanning all country) iii) year of sampling, etc.

Fifteen SITs were noted to be the most represented (n = 5 or more isolates) in this study (Table 2). SIT53/T1 was the most prevalent suggesting that this particular lineage is actively circulating in Chile. This finding is corroborated by findings regarding the worldwide scene which also states that this particular SIT is widely distributed (Table 2) [22]. SIT33/LAM3 (9.61%) and SIT42/LAM9 (9.39%) are present in Brazil in proportions similar to reported by us but in Peru and Paraguay the presence of these two SITs is only about 3% (Table 3). Fig 1B represents the main SITs obtained in this study (labeled circle) as well as the genetic evolution based on distance and allele/spacer change between patterns. In case of the MST, it is possible to ascertain major evolutionary pathways followed by the MTB lineages by observing the similarity between each strain. Further analysis leads us to observe that the Beijing family group is very far from the Euro American lineage present in the central nodes. We also concluded that the strains belonging to the LAM lineage have a high genetic diversity and that the T and Haarlem lineages are genetically closer.

Three additional published studies have analyzed the genetic diversity of MTB in Chile [25–27]. Although all three studies were limited by the fact that they were restricted to a single region of the country; they unanimously concluded that the LAM family was the most predominant lineage in Chile. A recently published work [27] showed that almost two-thirds of the circulating strains in Santiago city could be accounted for by the prevalence of LAM and T genotypes. The most frequent SITs that were seen to cluster spoligotypes were SIT33/LAM3 (10.7%), SIT53/T1 (8.7%), SIT50/H3 (7.8%) and SIT37/T3 (6.8%). The distribution pattern in the above mentioned study goes in parallel with data obtained in our study. Last but not least, our study further corroborated the high presence of SIT37/T3 (n = 23) in the Chilean population, a pattern that was usually found in Eastern Africa in the international database.

Currently, the National Control Program of Human Tuberculosis in Chile [13] has managed to position our country among nations that have managed to control this disease. The incidence of TB in Chile, when compared to other South American countries, is low at only 12.5 per 100,000. During the 2009–2013 period, the incidence value for Brazil was of 46, 44 in Paraguay and 124 in Peru (per 100,000). The situation in Peru, a neighboring country with which we share the border, attracts our attention because the differences in incidence rates are so overwhelming. This could be a possible effect of immigration and public health policies. As with other airborne diseases, the spread of MTB is facilitated by high population densities and crowded indoor environments that optimize the transmission of the pathogen. Other host-related factors that are implicated in increasing the chances of infection are: (i) immune suppression, (ii) smoking, (iii) poor nutrition, (iv) diabetes, and (v) respiratory comorbidities. All
these factors can also result in an increase in the risk of transitioning from latent to active MTB
[38]. This is a very important issue in the developing world where overcrowding, malnutrition
and HIV infection contribute to a high burden of the disease.

Drug resistant tests conducted as part of our study revealed that 87.7% (n = 402) of the isolates were pan-susceptible, 27.8% (n = 44) were resistant to any drugs, and 1.31% (n = 6) were MDR; for 1.31% (n = 6 isolates) the data could not be obtained. Three of the six strains that presented with multi drug resistance belonged to SIT53/T1. Conversely, LAM (which is the most prevalent lineage in our study), was apparently overrepresented among MDR cases in French departments of the Americas-based on data collection over an extended period of time [39], as well as in Peru [40]. Lastly, half of the SIT4015 isolates (of a lineage with unknown signature) showed drug resistance. This pattern should be carefully surveyed in coming years to see if its emergence is not linked to specific host factors since different ethnic groups may have varying susceptibilities to various MTB strains [41]. Indeed, both SIT4014/X1 (n = 7 isolates) and SIT4015 (Unknown signature; n = 7 isolates) have only been found in Chile so far. Considering that most of the samples belonging to these SITs were collected in the city of Concepcion and its vicinity, one should take into account associated factors such as susceptibility genetics of the host when studying emerging MTB clones. A similar case can be made for SIT720/LAM3, which is so far restricted to South America. In summary, our results underline the importance of tracking circulating and emerging MTB clones in conjunction with development of phylogenetical databases in order to maintain proper vigilance. These measures are a powerful way to safeguard against possible future outbreaks.

Although Spoligotyping can contribute to identifying outbreaks and tracking the spread of
disease, it has been known to overestimate clustering of isolates. Therefore, it is our belief that a
more polymorphic technique such as MIRU-VNTR technique would be able to impart greater
in-depth information that will help improve our knowledge of current epidemiological data in
Chile.

To conclude, in this study, spanning 15 regions of Chile, we have shown that the LAM, T, and Haarlem lineages together represent 88.2% of the isolates. This corroborates that the data obtained from the neighboring countries also show that the Euro-American lineages have a clear epidemiological superiority. The exclusive emergence of SIT4014 and SIT4015 indicates that conditions specific to Chile, along with the unique genetic makeup of the Chilean population, might have allowed for a possible co-evolution of certain genotypes leading to their relative success to cause disease.

**Supporting Information**

S1 Fig. Geographic distribution of isolates collected within of 15 regions of Chile (n = 458 strains).

(TIF)

S1 Table. Detailed demographic and genotyping information on *M. tuberculosis* isolates from Chile (n = 458 strains).

(PDF)

**Author Contributions**

Conceived and designed the experiments: JL FA JF.

Performed the experiments: JF TL JT FA JCH CA.

Analyzed the data: JL LA JF DC NR.
**Contributed reagents/materials/analysis tools:** JL JCH JF FA TL DC NR CA.

**Wrote the paper:** JT LA JF DC NR.

## References

1. "Global Tuberculosis Report 2013". WHO report 2013: Available: http://www.who.int/tb/publications/global_report/en/. Accessed 10 March 2015.

2. Barnes PF, Cave MD. Molecular epidemiology of tuberculosis. N Engl J Med. 2003 Sep 18; 349 (12):1149-56. Review. PMID: 13679530.

3. Couvin D, Rastogi N. 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, 2014. 12: p. 36-48. Available: https://pro.anses.fr/euroreference/Documents/ER12-RESEAUX-MycobacteriumEN.pdf

4. Alland D, Kalkut GE, Moss AR, McAdam RA, Hahn JA, Bosworth W, et al. Transmission of tuberculosis in New York City. An analysis by DNA fingerprinting and conventional epidemiologic methods. N Engl J Med. 1994 Jun 16; 330(24):1710-6. PMID: 7993412.

5. van Embden JD, Cave MD, Crawford JT, Dale JW, Eisenach KD, Gicquel B et al. Strain identification of Mycobacterium tuberculosis by DNA fingerprinting: recommendations for a standardized methodology. J Clin Microbiol. 1993 Feb; 31(2):406-9. PMID: 8381814

6. Cronin WA, Golub JE, Magder LS, Baruch NG, Lathan MJ, Mukasa LN, et al. Epidemiologic usefulness of spoligotyping for secondary typing of Mycobacterium tuberculosis isolates with low copy numbers of IS6110. J Clin Microbiol. 2001 Oct; 39(10):3709-11. PMID: 11574598

7. Supply P, Lesjean S, Savine E, Kremer K, van Soolingen D, Locht C. Automated high-throughput genotyping for study of global epidemiology of Mycobacterium tuberculosis based on mycobacterial interspersed repetitive units. J Clin Microbiol. 2001 Oct; 39(10):3563-71. PMID: 11574573

8. Kamerbeek J, Schouls L, Kolk A, van Agterveld M, van Soolingen D, Kuijper S, et al. Simultaneous detection and strain differentiation of Mycobacterium tuberculosis for diagnosis and epidemiology. J Clin Microbiol. 1997 Apr; 35(4):907-14. PMID: 9157152

9. Kremer K, Arnold C, Cataldi A, Gutiérrez MC, Haas WH, Panaiotov S, et al. Discriminatory power and reproducibility of novel DNA typing methods for Mycobacterium tuberculosis complex strains. J Clin Microbiol. 2005 Nov; 43(11):5628-38. PMID: 16272496

10. Supply P, Allix C, Lesjean S, Cardoso-Oelemann M, Rüsch-Gerdes S, Willery E, et al. Proposal for standardization of optimized mycobacterial interspersed repetitive unit-variable-number tandem repeat typing of Mycobacterium tuberculosis. J Clin Microbiol. 2006 Dec; 44(12):4498-510. Epub 2006 Sep 27. PMID: 17005759

11. García de Viedma D, Mokrousov I, Rastogi N. Innovations in the molecular epidemiology of tuberculosis. Enferm Infecc Microbiol Clin. 2011 Mar; 29 Suppl 1:8-13 doi: 10.1016/S0213-005X(11)70012-X PMID: 21420561

12. Tuberculosis in the Americas Regional Report 2014: Epidemiology, Control, and Financing. Available: http://www.paho.org/hq/index.php?option=comTopics&view=readall&cid=5504&Itemid=40242&lang=es. Accessed 14 September 2015.

13. "Tuberculosis Informe de Situación Chile: 2013". September 2014, Ministerio de Salud de Chile: Available: http://web.minsal.cl/sites/default/files/files/TUBERCULOSIS-INFORME%202013.pdf. Accessed 10 September 2015.

14. Gomes HM, Elias AR, Oelemann MA, Pereira MA, Montes FF, Marsico AG, et al. Spoligotypes of Mycobacterium tuberculosis complex isolates from patients residents of 11 states of Brazil. Infect Genet Evol. 2012 Jun; 12(4):649-56. doi: 10.1016/j.meegid.2011.08.027 PMID: 21907830.

15. Abadía E, Sequera M, Ortega D, Méndez MV, Escalona A, Da Mata O, et al. Mycobacterium tuberculosis ecology in Venezuela: epidemiologic correlates of common spoligotypes and a large clonal cluster defined by MIRU-VNTR-24. BMC Infect Dis. 2009 Aug 6; 9:122. doi: 10.1186/1471-2334-9-122 PMID: 19660112

16. Candia N, Lopez B, Zozio T, Carrivale M, Diaz C, Russomando G, et al. First insight into Mycobacterium tuberculosis genetic diversity in Paraguay. BMC Microbiol. 2007 Aug 8; 7:75. PMID: 17686181

17. Cerezo I, Jiménez Y, Hernandez J, Zozio T, Murcia MI, Rastogi N. A first insight on the population structure of Mycobacterium tuberculosis complex as studied by spoligotyping and MIRU-VNTRs in Bogotá, Colombia. Infect Genet Evol. 2012 Jun; 12(4):657-63. doi: 10.1016/j.meegid.2011.07.006 PMID: 21784173.

18. Dalla Costa ER, Lazzarini LC, Perizzolo PF, Díaz CA, Spies FS, Costa LL, et al. Mycobacterium tuberculosis of the RDRio genotype is the predominant cause of tuberculosis and associated with multidrug
resistance in Porto Alegre City, South Brazil. J Clin Microbiol. 2013 Apr; 51(4):1071-7. doi: 10.1128/JCM.01511-12 PMID: 23325819

19. Realpe T, Correa N, Rozo JC, Ferro BE, Gomez V, Zapata E, et al. Population structure among mycobacterium tuberculosis isolates from pulmonary tuberculosis patients in Colombia. PLoS One. 2014 Apr 18; 9(4):e93848. doi: 10.1371/journal.pone.0093848. Erratum in: PLoS One. 2014;9(9):e109456. PMID: 24747767

20. Taype CA, Agapito JC, Accinelli RA, Espinoza JR, Godreuil S, Goodman SJ, et al. Genetic diversity, population structure and drug resistance of Mycobacterium tuberculosis in Peru. Infect Genet Evol. 2012 Apr; 12(3):577-85. doi: 10.1016/j.meegid.2012.02.002 PMID: 22342744

21. Alvarez N, Haft D, Hurtado UA, Robledo J, Rouzaud F. Whole-Genome Sequence of a Beijing Exten-
sively Drug-Resistant Mycobacterium tuberculosis Clinical Isolate from Buenaventura, Colombia.
Genome Announc. 2016 Jan 14; 4(1). pii: e01549-15. doi: 10.1128/genomeA.01549-15 PMID:
26769935; PubMed Central PMCID: PMC4714117.

22. Vasconcellos SE, Acosta CC, Gomes LL, Conceição EC, Lima KV, de Araujo MI, et al. Strain classification of Mycobacterium tuberculosis isolates in Brazil based on genotypes obtained by spoligotyping, mycobacterial interspersed repetitive unit typing and the presence of large sequence and single nucleotide polymorphism. PLoS One. 2014 Oct 14; 9(10):e107747. doi: 10.1371/journal.pone.0107747 PMID: 25314118

23. Demay C, Liens B, Burguière T, Hill V, Couvin D, Millet J, et al. SITVITWEB—a publicly available inter-national multimarker database for studying Mycobacterium tuberculosis genetic diversity and molecular epidemiology. Infect Genet Evol. 2012 Jun; 12(4):755–66. doi: 10.1016/j.meegid.2012.02.004 PMID: 22365971.

24. Shabbeer A, Ozcaglar C, Yener B, Bennett KP (2012) Web tools for molecular epidemiology of tubercu-losis. Infect Genet Evol; 12(4): 767-81. doi: 10.1016/j.meegid.2011.08.019 PMID: 21903179

25. Mancilla EM, Martínez HA, Palavecino BC, Rehren SG, Lucero LP, León RG, Zárraga O AM. [Genetic variants of Mycobacterium tuberculosis isolated from patients of the Xth Region of Chile]. Rev Chilena Infectol. 2006 Sep; 23(3):220-5. Spanish. PMID: 16896494.

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

27. Balcells ME, García P, Meza P, Peña C, Cifuentes M, Couvin D, 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 Feb 11; 10(2):e0118007. doi: 10.1371/journal.pone.0118007 PMID: 25671320

28. Couvin D, Rastogi N. Tuberculosis—A global emergency: Tools and methods to monitor, understand, and control the epidemic with specific example of the Beijing lineage. Tuberculosis (Edinb). 2015 Jun; 95 Suppl 1:S177-89. doi: 10.1016/j.tube.2015.02.023 Epub 2015 Feb 13.

29. Laszlo A, Rahman M, Raviglione M, Bustreo F. 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. 1997 Jun; 1(3):231-8. PMID: 9432369.

30. Araya P, Velasco M, Tognarelli J, Arias F, Leiva T, Sccapatticio A, et al. [Detection of genes associated with drug resistance in Mycobacterium tuberculosis strains isolated in Chile]. Rev Med Chil. 2011 Apr; 139(4):467-73. /S0034-98872011000400008. Epub 2011 Aug 25. Spanish. PMID: 21879185.

31. Reyes JF, Francis AR, Tanaka MM. Models of deletion for visualizing bacterial variation: an application to tuberculosis spoligotypes. BMC Bioinformatics. 2008 Nov 27; 9:496. doi: 10.1186/1471-2105-9-496 PMID: 19036166

32. Ellson J., et al., Graphviz and dynagraph-static and dynamic graph drawing tools. In: Graph Drawing Software. 2004. Jünger M. Mutzel P. Editors. Springer-Verlag. p. 127–148 doi: 10.1007/978-3-642-18638-7

33. Sheen P, Couvin D, Grandjean L, Zimic M, Dominguez M, Luna G, et al. Genetic diversity of Mycobac-
terium tuberculosis in Peru and exploration of phylogenetic associations with drug resistance. PLoS
One. 2013 Jun 24; 8(6):e65873. doi: 10.1371/journal.pone.0065873 PMID: 23826083

34. Monteserin J, Camacho M, Barrera L, Palomino JC, Ritacco V, Martin A. Genotypes of Mycobacterium tuberculosis in patients at risk of drug resistance in Bolivia. Infect Genet Evol. 2013 Jul; 17:195-201. doi: 10.1016/j.meegid.2013.04.010 PMID: 23603419.

35. Gonzalo X, Ambroggi M, Cordova E, Brown T, Poggi S, Drobniewski F. Molecular Epidemiology of *Mycobacterium tuberculosis*, Buenos Aires, Argentina. Emerging Infectious Diseases. 2011; 17 (3):528-531. doi: 10.3201/eid1703100394 PMID: 21392451
36. Ritacco V, López B, Ambroggi M, Palmero D, Salvadores B, Gravina E, et al. HIV infection and geographically bound transmission of drug-resistant tuberculosis, Argentina. Emerg Infect Dis. 2012 Nov; 18(11):1802–10. doi: 10.3201/eid1811.120126 PMID: 23092584

37. Machado LN, Marcondes NR, Leite CQ, Santos AC, Pavan FR, Baldin VP, et al. First baseline of circulating genotypic lineages of Mycobacterium tuberculosis in patients from the Brazilian borders with Argentina and Paraguay. PLoS One. 2014 Sep 9; 9(9):e107106. doi: 10.1371/journal.pone.0107106 PMID: 25202909

38. Zumla A, Nahid P, Cole ST. Advances in the development of new tuberculosis drugs and treatment regimens. Nat Rev Drug Discov. 2013 May; 12(5):388–404. doi: 10.1038/nrd4001 Review. PMID: 23629506

39. Millet J, Streit E, Berchel M, Bomer AG, Schuster F, Paasch D, et al. 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 Res Int. 2014; 2014:689852. doi: 10.1155/2014/689852 Epub 2014 Mar 13. PMID: 24738068

40. Grandjean L, Iwamoto T, Lithgow A, Gilman RH, Arikawa K, Nakanishi N, et al. The Association between Mycobacterium Tuberculosis Genotype and Drug Resistance in Peru. PLoS One. 2015; 10 (5):e0126271. doi: 10.1371/journal.pone.0126271 PMID: 25984723

41. Gagneux S, Small PM. Global phylogeography of *Mycobacterium tuberculosis* and implications for tuberculosis product development. Lancet Infect Dis 2007; 7(5):328–337. PMID: 17448936
<figure><img src="image_6.png" /></figure>

<table><thead><tr><th>Regions of Chile</th><th>City</th><th>N° Samples</th></tr></thead><tbody><tr><td>XV</td><td>Arica</td><td>14</td></tr><tr><td>I</td><td>Iquique</td><td>34</td></tr><tr><td>II</td><td>Antofagasta</td><td>11</td></tr><tr><td>II</td><td>Calama</td><td>1</td></tr><tr><td>III</td><td>Vallenar</td><td>1</td></tr><tr><td>III</td><td>Copiapó</td><td>5</td></tr><tr><td>IV</td><td>Ovalle</td><td>6</td></tr><tr><td>V</td><td>Coquimbo</td><td>11</td></tr><tr><td>V</td><td>Valparaíso</td><td>29</td></tr><tr><td>V</td><td>Viña del Mar</td><td>2</td></tr><tr><td>XIII</td><td>Santiago</td><td>181</td></tr><tr><td>VI</td><td>Rancagua</td><td>20</td></tr><tr><td>VII</td><td>Talca</td><td>12</td></tr><tr><td>VIII</td><td>Concepción</td><td>34</td></tr><tr><td>VIII</td><td>Los Angeles</td><td>16</td></tr><tr><td>VIII</td><td>Talcahuano</td><td>12</td></tr><tr><td>IX</td><td>Osorno</td><td>15</td></tr><tr><td>IX</td><td>Temuco</td><td>18</td></tr><tr><td>XIV</td><td>Valdivia</td><td>6</td></tr><tr><td>X</td><td>Castro</td><td>9</td></tr><tr><td>X</td><td>Pto. Montt</td><td>12</td></tr><tr><td>XI</td><td>Coyhaique</td><td>2</td></tr><tr><td>XII</td><td>Pta. Arenas</td><td>7</td></tr></tbody><tfoot><tr><td colspan="2">TOTAL</td><td>458</td></tr></tfoot></table>
<table><tr><td>1439-2011-S-924</td><td>2011</td><td>Female</td><td>Santiago</td><td>77777720000001</td><td>Or17</td><td>Unknown</td><td>Pansusceptible</td></tr><tr><td>1451-2011-S-934</td><td>2011</td><td>Male</td><td>Iquique</td><td>75777607700771</td><td>Or18</td><td>Unknown</td><td>Pansusceptible</td></tr><tr><td>1468-2011-S-946</td><td>2011</td><td>Female</td><td>Valparaiso</td><td>776177617760771</td><td>Or19</td><td>T1</td><td>Pansusceptible</td></tr><tr><td>1469-2011-S-947</td><td>2011</td><td>Male</td><td>Calama</td><td>777767707760471</td><td>Or20</td><td>T1</td><td>Pansusceptible</td></tr><tr><td>149-2012-S-123</td><td>2012</td><td>Female</td><td>Santiago</td><td>700017377760771</td><td>Or21</td><td>T4</td><td>Pansusceptible</td></tr><tr><td>1492-2011-S-966</td><td>2011</td><td>Male</td><td>Santiago</td><td>71107737760771</td><td>Or22</td><td>T1</td><td>Pansusceptible</td></tr><tr><td>1495-2011-S-969</td><td>2011</td><td>Male</td><td>Valparaiso</td><td>700777777760671</td><td>Or23</td><td>T1</td><td>Pansusceptible</td></tr><tr><td>1829-2011-S-1265</td><td>2011</td><td>Male</td><td>Santiago</td><td>777734077460771</td><td>Or24</td><td>T3</td><td>Pansusceptible</td></tr><tr><td>1830-2011-S-1266</td><td>2011</td><td>Male</td><td>Santiago</td><td>776176207760771</td><td>Or25</td><td>LAM3</td><td>Pansusceptible</td></tr><tr><td>1837-2011-S-1272</td><td>2011</td><td>Female</td><td>Santiago</td><td>757777774120771</td><td>Or26</td><td>H3</td><td>Pansusceptible</td></tr><tr><td>185-2012-S-150</td><td>2012</td><td>Female</td><td>Santiago</td><td>676137607760771</td><td>Or27</td><td>LAM3</td><td>Pansusceptible</td></tr><tr><td>1944-2011-S-1352</td><td>2011</td><td>Male</td><td>Talcahuano</td><td>677777417760771</td><td>Or28</td><td>T1</td><td>Pansusceptible</td></tr><tr><td>1951-2011-S-1355</td><td>2011</td><td>Female</td><td>Osorno</td><td>000000010020771</td><td>Or29</td><td>H3</td><td>Pansusceptible</td></tr><tr><td>1953-2011-S-1357</td><td>2011</td><td>Female</td><td>Osorno</td><td>776177417760731</td><td>Or30</td><td>T2</td><td>Pansusceptible</td></tr><tr><td>1954-2011-S-1358</td><td>2011</td><td>Male</td><td>Osorno</td><td>577777417760751</td><td>Or31</td><td>T1</td><td>Pansusceptible</td></tr><tr><td>1981-2011-S-1377</td><td>2011</td><td>Male</td><td>Valparaiso</td><td>777777767480770</td><td>Or32</td><td>T1</td><td>Pansusceptible</td></tr><tr><td>305-2012-S-232</td><td>2012</td><td>Male</td><td>Concepcion</td><td>777770017560771</td><td>Or33</td><td>T1</td><td>Pansusceptible</td></tr><tr><td>372-2012-S-273</td><td>2012</td><td>Male</td><td>Coquimbo</td><td>723777207760771</td><td>Or34</td><td>LAM9</td><td>R (SM)</td></tr><tr><td>38-2012-S-32</td><td>2012</td><td>Female</td><td>Coquimbo</td><td>557727637740771</td><td>Or35</td><td>T3</td><td>Pansusceptible</td></tr><tr><td>39-2012-S-33</td><td>2012</td><td>Female</td><td>Concepcion</td><td>401767637740771</td><td>Or36</td><td>T1</td><td>Pansusceptible</td></tr><tr><td>50-2012-S-43</td><td>2012</td><td>Male</td><td>Santiago</td><td>777767637700771</td><td>Or37</td><td>Unknown</td><td>Pansusceptible</td></tr><tr><td>61-2012-S-53</td><td>2012</td><td>Female</td><td>Concepcion</td><td>576127607740771</td><td>Or38</td><td>LAM3</td><td>Pansusceptible</td></tr><tr><td>63-2012-S-55</td><td>2012</td><td>Male</td><td>Iquique</td><td>777767637740611</td><td>Or39</td><td>T1</td><td>Pansusceptible</td></tr><tr><td>64-2012-S-56</td><td>2012</td><td>Male</td><td>Puerto Montt</td><td>557767207740771</td><td>Or40</td><td>LAM9</td><td>Pansusceptible</td></tr><tr><td>66-2012-S-58</td><td>2012</td><td>Female</td><td>Valparaiso</td><td>676167607740771</td><td>Or41</td><td>LAM3</td><td>Pansusceptible</td></tr><tr><td>67-2012-S-59</td><td>2012</td><td>Male</td><td>Valparaiso</td><td>577767637740601</td><td>Or42</td><td>T1</td><td>Pansusceptible</td></tr><tr><td>676-2012-S-449</td><td>2012</td><td>Male</td><td>Santiago</td><td>663777777320771</td><td>Or43</td><td>H3</td><td>Pansusceptible</td></tr><tr><td>68-2012-S-60</td><td>2012</td><td>Male</td><td>Valparaiso</td><td>777674037540771</td><td>Or44</td><td>T1</td><td>Pansusceptible</td></tr><tr><td>883-2011-S-536</td><td>2011</td><td>Male</td><td>Arica</td><td>776177637720771</td><td>Or45</td><td>H3</td><td>Pansusceptible</td></tr><tr><td>894-2011-S-547</td><td>2011</td><td>Male</td><td>Antofagasta</td><td>777777637700771</td><td>Or46</td><td>Unknown</td><td>Pansusceptible</td></tr><tr><td>924-2011-S-568</td><td>2011</td><td>Male</td><td>Santiago</td><td>416077703770371</td><td>Or47</td><td>Unknown</td><td>Pansusceptible</td></tr><tr><td>977-2011-S-601</td><td>2011</td><td>Female</td><td>Iquique</td><td>7777376057760771</td><td>Or48</td><td>LAM5</td><td>Pansusceptible</td></tr><tr><td>978-2011-S-602</td><td>2011</td><td>Male</td><td>Iquique</td><td>576177607760731</td><td>Or49</td><td>LAM3</td><td>R (PZA)</td></tr><tr><td>984-2011-S-607</td><td>2011</td><td>Male</td><td>Talca</td><td>776077607720371</td><td>Or50</td><td>H3</td><td>Pansusceptible</td></tr></table>

<footnote>* NEW SITs are followed by an asterisk (*) and highlighted in yellow; Orphan spoligotypes are highlighted in blue.</footnote>

<footnote>** Drug resistance: pansusceptible n=402; NA, not available n=6; MDR, multidrug resistant n=6; R, any resistance, n=44. Abbreviations for first line drugs: INH, Isoniazid; RIF, Rifampin; SM, Streptomycin; PZA, pyrazinamide.</footnote>