# 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

Sidra E. G. Vasconcellos, Chyntia Carolina Acosta, Lia Lima Gomes, Emily Costa Conceicao, Karla Valéria Lima, Marcelo Ivens de Araujo, Maria de Lourdes Leite, Flavio Tannure, Paulo Cesar de Souza Caldas, Harrison M Gomes, et al.

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

► To cite this version:

Sidra E. G. Vasconcellos, Chyntia Carolina Acosta, Lia Lima Gomes, Emily Costa Conceicao, Karla Valéria Lima, 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, 9 (10), pp.e107747. 10.1371/journal.pone.0107747. pasteur-01153012

HAL Id: pasteur-01153012

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

Submitted on 19 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_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# 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

Sidra E. G. Vasconcellos¹,⁸, Chyntia Carolina Acosta², Lia Lima Gomes¹, Emily Costa Conceição³, Karla Valéria Lima³, Marcelo Ivens de Araujo¹, Maria de Lourdes Leite⁴, Flávio Tannure⁴, Paulo Cesar de Souza Caldas⁵, Harrison M. Gomes¹, Adalberto Rezende Santos¹, Michel K. Gomgnimbou⁶, Christophe Sola⁶, David Couvin⁷, Nalin Rastogi⁷, Neio Boechat⁸,⁹, Philip Noel Suffys¹*

¹ Laboratory of Molecular Biology Applied to Mycobacteria, Oswaldo Cruz Institute, FIOCRUZ, Rio de Janeiro, Rio de Janeiro, Brazil, ² Laboratory of Cellular Microbiology, Oswaldo Cruz Institute, FIOCRUZ, Rio de Janeiro, Rio de Janeiro, Brazil, ³ Instituto Evandro Chagas, Section of Bacteriology and Mycology, Belém, Pará, Brazil, ⁴ Hospital Municipal Rafael de Paula Souza, Municipal Secretary of Health, Rio de Janeiro, Rio de Janeiro, Brazil, ⁵ Centro de Referência Professor Hélio Fraga, Escola Nacional de Saúde Pública Sergio Arouca, FIOCRUZ, Rio de Janeiro, Rio de Janeiro, Brazil, ⁶ CNRS–Université Paris–Sud, Institut de Génétique et Microbiologie–Infection Genetics Emerging Pathogens Evolution Team, Orsay, France, ⁷ Supranational TB Reference Laboratory, Unité de la Tuberculose et des Mycobactéries, Institut Pasteur de Guadeloupe, Abymes, Guadeloupe, France, ⁸ Multidisciplinary Research Laboratory, University Hospital Clementino Fraga Filho – HUCEE, Federal University of Rio de Janeiro, Rio de Janeiro, Rio de Janeiro, Brazil, ⁹ Graduate Program in Clinical Medicine, Faculty of Medicine, University Hospital Clementino Fraga Filho, Rio de Janeiro, Rio de Janeiro, Brazil

## Abstract

Rio de Janeiro is endemic for tuberculosis (TB) and presents the second largest prevalence of the disease in Brazil. Here, we present the bacterial population structure of 218 isolates of *Mycobacterium tuberculosis*, derived from 186 patients that were diagnosed between January 2008 and December 2009. Genotypes were generated by means of spoligotyping, 24 MIRU-VNTR typing and presence of fbpc¹⁰³, RDR² and RD174. The results confirmed earlier data that predominant genotypes in Rio de Janeiro are those of the Euro American Lineages (99%). However, we observed differences between the classification by spoligotyping when comparing to that of 24 MIRU-VNTR typing, being respectively 43.6% vs. 62.4% of LAM, 34.9% vs. 9.6% of T and 18.3% vs. 21.5% of Haarlem. Among isolates classified as LAM by MIRU typing, 28.0% did not present the characteristic spoligotype profile with absence of spacers 21 to 24 and 32 to 36 and we designated these conveniently as "LAM-like", 79.3% of these presenting the LAM-specific SNP fbpc¹⁰³. The frequency of RDR² and RD174 in the LAM strains, as defined both by spoligotyping and 24 MIRU-VNTR loci, were respectively 11% and 15.4%, demonstrating that RD174 is not always a marker for LAM/RDR² strains. We conclude that, although spoligotyping alone is a tool for classification of strains of the Euro-American lineage, when combined with MIRU-VNTRs, SNPs and RD typing, it leads to a much better understanding of the bacterial population structure and phylogenetic relationships among strains of *M. tuberculosis* in regions with high incidence of TB.

**Citation:** Vasconcellos SEG, Acosta CC, Gomes LL, Conceição EC, Lima KV, et al. (2014) 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 9(10): e107747. doi:10.1371/journal.pone.0107747

**Editor:** Igor Mokrousov, St. Petersburg Pasteur Institute, Russian Federation

**Received:** January 1, 2014; **Accepted:** August 21, 2014; **Published:** October 14, 2014

**Copyright:** © 2014 Vasconcellos 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.

**Funding:** This work was supported by funding by Brazilian funding agencies PAPES/CNPq (407624/2012-0), PROEP/CNPq of the Oswaldo Cruz Institute (IOC), and CNPq PhD grant (142958/2008-5). 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.

* Email: psuffys@gmail.com

## Introduction

Tuberculosis (TB) is an infectious disease with an effective treatment but remains an important cause of morbidity and mortality in many countries. In Brazil, 73.000 new TB cases and 4.600 cases that deceased were registered in 2011. The southeast region has 44.8% of all cases reported in the country and the state of Rio de Janeiro has the highest disease incidence (72.3/100,000) and mortality rate (5.7/100,000) in the country [1].

A couple of years ago, *Mycobacterium tuberculosis* (Mtb) strains designated as RDR² were reported as being a predominant genotype in Rio de Janeiro [2]. These strains have a deletion of 26.3 kb and seem to be restricted to the Latin American-Mediterranean (LAM) family. The RDR² strains have been associated with higher levels of recent transmission and of Multi-Drug Resistance (MDR) but there are contradictory data about their relation with disease severity [3], [4], [5], [6]. Another Region of Difference called RD174 has been described as a co-
marker of RD<sup>Rio</sup> and as a marker for the LAM type [7], [8], [9], [10] and strains with the RD174 deletion were found to have an increased secondary case rate ratio in San Francisco [10].

The most commonly used genotyping methods to characterize *M. tuberculosis* Complex (MTBC) isolates are IS6110-RFLP, 24-locus MIRU-VNTR, spoligotyping, detection of Single Nucleotide Polymorphisms (SNPs) and Large Sequence Polymorphisms (LSP). Typing of MTBC by IS6110-RFLP was considered the gold-standard during more than a decade [11], but is labor intensive. Currently, the two most commonly used methods are spoligotyping and Mycobacterium Interspersed Repetitive Units – Variable Number of Tandem Repeats (MIRU-VNTR) typing, based on the variability of the direct repeat (DR) locus [12] and of minisatellites [13], [14], respectively. However, for epidemiologic studies, spoligotyping overestimates the proportion of clustered strains and should be used in combination with MIRU-VNTR typing [15]. Spoligotyping alone is also unable to reveal exact phylogenetic relationships between MTBC strains, particularly for the classification of Euro-American strains, part of the PGG2/PGG3 group [9], [16], [17], [18]. These phylogenetic subtypes include the genotype families LAM, H, T, S and X, frequently observed in Brazil, in South America and elsewhere [2], [4]. [8], [18], [19], [20], [21], [22], [23].

The generalized use of spoligotyping and MIRU-VNTR resulted in the construction of large international databases of genotypes, allowing the study of global distribution and phylogenetic analysis of the distribution of *M. tuberculosis* worldwide [23], [24], [25], [26], [27], [28], [29]. Combined MIRU-VNTR typing and spoligotyping also helps in revealing epidemiologically meaningful clonal diversity of *M. tuberculosis* strain lineages and is useful to explore internal phylogenetic ramifications [8]. In addition, SNPs and LSPs represent robust markers for inferring phylogenies and for strain classification [30], [31], [32], [33], [34], [35]. Besides the presence of RDs, that of determinate SNPs is also a characteristic of LAM strains, [8], [9], [10] such as *fbpC*<sup>103</sup> in codon 103 (G to A) of the gene encoding antigen 85 Complex Ag85<sup>c</sup> (Rv0129c) resulting in a Glu103Asp amino acid replacement in a protein that is involved in biosynthesis of cell wall components of *M. tuberculosis* [7], [36]. Other LAM-specific SNPs are (i) the silent C to G *ligB* mutation at genome position 3426795 [37], validated by the Abadia *et al.* [17] and (ii) the SNP at 240 codon of *mgtC* (C to T) that was discovered by Homolka *et al.* [38].

The objective of this study was to evaluate 24 MIRU-VNTR typing for alternative classification of strain families/lineages by means of phylogenetic tree building and comparison with MIRU databases (MIRU-VNTRPlus), as compared with spoligotypes-based classification (SITVITWEB and SITVIT2), and increase the consistency of the analysis by including a SNP and two LSPs. In addition, the differentiating power of both techniques in a population that is almost exclusively of the Euro-American lineages.

# Materials and Methods

## Study setting

The study was based on a convenience sampling of TB patients diagnosed between January 2008 and December 2009 at the “Hospital Municipal Raphael de Paula e Sousa”, Curicica, Rio de Janeiro, Brazil. Two hundred eighteen clinical isolates were obtained from 188 patients and 28 cases had more than one isolate, including 25 patients with two and three patients with three isolates. Multiple isolates from the same patient were

included to evaluate the consistency of the results of genotyping and to verify eventual multiple infections.

## Ethics statement

This study was approved by the Research Ethics Committee of the Municipal Health and Civil Defense of the city of Rio de Janeiro Number 160/09 CAAE: 0182.0.314.000-9. Isolates of this study were obtained through bacteriological culture from clinical specimens of patients diagnosed at the Hospital Raphael and Paula Souza as part of routine diagnosis and drug susceptibility testing. The ethical committee stated that, as genotyping of *M. tuberculosis* isolates is complementary to routine diagnosis, there was no need for a written or verbal consent. No information other than that provided for the diagnosis was used.

## Culture, DNA extraction and identification of Mycobacterium tuberculosis complex

Sputum samples were cultured on Löwenstein-Jensen (LJ) medium following standard microbiological laboratory procedures as a routine of the Clinical Analysis Laboratory of “Hospital Raphael de Paula e Souza”. After cultivation, bacterial mass was re-suspended in 400 µl of sterile distilled water and heat inactivated at 90°C for one hour, followed by DNA extraction and purification using the CTAB method [39].

## Spoligotyping

Spoligotyping was performed either as described by Kamerbeek *et al.* (1997) [12], using commercially available kit from Ocimum Biosolutions (Hyderabad, India) according to the manufacturer's instructions or using microbead-based hybridization assay as described by Zhang *et al.* (2009) [40].

## MIRU-VNTR typing

Amplification of 24 MIRU-VNTR *loci* was performed by using a commercial typing kit (Genoscreen, Lille, France) and automated MIRU-VNTR analysis performed as previously described [14]. Fragment size of the amplicons was analyzed on a ABI 3730 DNA sequence analyzer (Applied Biosystems, California, USA) and number of copies of each locus was determined by automated assignment using the Genemapper 4.0 software (Applied Biosystems, California, USA). In the case of doubtful results, the size of the repeats was double checked by size estimation as compared to a DNA ladder (50 and 100 bp) and the positive control (H37Rv) on agarose gels and by comparing to a reference table as described [13].

## PCR-RFLP of *fbpC*<sup>103</sup>

For characterization of SNP *fbpC*<sup>103</sup>, we adapted the procedure described by Gibson *et al.* in 2008 [7]. Amplifications were performed in 50 µl reactions containing 40 pmol each of the primers Ag85C103F (5-CTG GCT GTT GCC CTG ATA CTG CGA GGG CCA-3') and Ag85C103R (5-CGA GCA GCT TCT GCG GCC ACA ACG TT-3'), 2 mM MgCl2, 0.2 mM dNTPs, 1 U Taq DNA polymerase (Invitrogen, Brazil), buffer (10 mM Tris-HCl, 1, 5 mM MgCl2, 50 mM KCl, pH 8.3), 5% DMSO (v/v) and 10 ng of target DNA. Amplification was performed in a Vieriti Thermal Cycler (Applied Biosystems, Foster City, CA), starting with denaturation at 95°C for 5 min, followed by 45 cycles of 1 min at 95°C, 1 min at 60°C, 4 min at 72°C and final extension for 10 min at 72°C. The amplified products of 519 bp were analyzed on 2% agarose gels after staining with ethidium bromide. Fifteen microliters of the amplified products were subjected to enzymatic digestion with 1 U of *MnI* (New England
observed in Brazil. Presently, we also observed five SITs, 136 (LAM5), 189 (T1), 268 (H3), 2107 (unknown) and 2322 (T2), that had not been reported before in Brazil. When comparing our results with the previous and the present version of the SITVIT database, we noticed some differences in classification (Table S1). Thirty-six isolates could be classified using the new version and among these nine new SITs were observed: 3125 (LAM6), 3504 (LAM6), 3904 (LAM6) 3905 (T1) 3906 (T3), 3907 (T1), 3908 (LAM2; originally designated as SIT2560 Unknown in SITVIT-WEB), 3909 (LAM8) and 3910 (T1). In addition, nineteen orphan profiles were classified and SIT268, SIT1241 and SIT2539, without classification in SITVITWEB, were now classified respectively as H3, LAM9 and LAM1. In addition, some SIT showed changes in clade and/or subclade classification (SITs 73, 102, 159, 451, 742 and 777). Further information is detailed in Table S1.

## Classification based on MIRU-VNTR analysis

The 218 isolates were also submitted to 24 loci MIRU-VNTR typing with the objective to confirm or not spoligotype-based classification of existing and new sublineages and to correct eventual miss classification that could have occurred using the latter technique. We observed 208 different patterns and seven clusters with two to four strains (clustering rate = 1.37%) (Figure S1). The allelic diversity of each MIRU-VNTR locus for 218 isolates was evaluated and classified as either highly (HGDI>0.6), moderately (0.6≤0.3) and poorly discriminative (HGDI<0.3), according to Sola *et al.* [15], as summarized in Table 2. Eight loci were highly discriminatory (QUB4156, QUB11, MTUB04, MIRU10, MIRU40, MIRU 23, MIRU26 and MIRU21) ten were moderate (MTUB30, MTUB39, MIRU16, ETR-C, ETR-A, QUB26, MTUB34, MIRU31, MIRU29 AND ETR-B) and six were poorly discriminative (MIRU24, MIRU29, MIRU04, MIRU27, MIRU20 and MIRU2). In the Table 2, we can observe a difference in discriminatory power of some loci according to lineage by Spoligotyping.

The patterns were classified based on the database tool that allows construction of a Neighbor-Joining based phylogenetic tree, visualizing proximity of a particular genotype with that of a set of reference strains to the genotype family level. One hundred thirty-six isolates (62.4%) were classified as LAM (127 patterns), 47 (21.5%) as Haarlem (45 patterns), 21 (9.6%) as T (21 patterns), 11 (5%) as S (11 patterns), and one isolate (0.45%) each as Beijing and EAI. Considerable differences were observed between Spoligotyping and 24 MIRU-VNTR loci-based classifications, even after excluding eventual small typing errors by repeating the assays. The differences observed between the two classifications are presented in the table 3 and in the Confusion Matrix (Figure 1). The precision, accuracy, sensitivity and error rate were respectively 0.74, 0.64, 0.82 and 0.36. The T lineage showed the highest incongruence rate related to classification (sensitivity = 0.26).

## Combined analysis of 24 MIRU-VNTR and spoligotyping

**The LAM family.** For better understanding of the population structure of the Mtb strains, dendrograms of spoligotyping and 24 MIRU-VNTRs, either separately or as combined patterns were constructed using BioNumerics. Among the 136 (62.4%) isolates that were classified as LAM based on MIRU-VNTR analysis, 99 (72.8%) exhibited the characteristic LAM spoligotype profile with absence of spacers 21 to 24 and 33–36. Among these however, four spoligotype isolates were designated to the T clade by SITVIT2 (SITs: 306/n = 3, 1688/n = 1, 3905/n = 1) and one as T2/orphan. Another thirty-seven strains (27.2%) defined as LAM by MIRU-VNTRs did not present the lack of spacers 21–24 and

were therefore designated as “LAM-like”. Among isolates initially classified as LAM by spoligotyping (n = 95), 91 were confirmed as LAM (95.8%/TPR = 0.96) by 24 loci MIRU-VNTR (Figure 1). Six spoligotypes referred to as “Unknown” in SITVIT2 (SITs 132/n = 1; 1952/n = 1; 2107/n = 2; 2511/n = 2; 2535/n = 1 and 2548/n = 1) were also classified as LAM. Using SpotClust, SIT132 and SIT2511 were classified as EAI and SIT1952 as Haarlem. All Orphans patterns classified as LAM by SITVIT2 (n = 10) were also classified as LAM according to the MIRU-VNTR based phylogenetic tree.

Signatures based on MIRU copy number for “real” LAM (n = 99) types showed that the majority had two copies of MIRU04 (94; 95%) and two MIRU20 (n = 96; (97%), one copy of MIRU24 (n = 91; 92%), two copies of MIRU39 and ETR-A (n = 93; 94%) and 86 (87%) one copy of MTUB30. This indicates that the 221221 combination of these six alleles is representative for LAM strains, except for LAM3, presenting two copies of MTUB30 (221222).

**The “LAM-like” isolates.** The thirty-seven isolates that did not present the typical LAM spoligotype signature (absence of spacers 21 to 24 and 33–36) were positioned within the LAM branch in the neighbor-joining tree using MIRU-VNTRPlus and were conveniently designated as “LAM-like”. These had been classified by spoligotyping earlier as T (n = 24), H (n = 10) and three “Unknown” patterns (SIT2511/n = 2 and SIT1952/n = 1; classified by SpotClust as EAI and H1 respectively). Ninety-one percent of these strains showed one copy of MIRU24 and two copies of MIRU04. Intriguingly, when including both 24 MIRU-VNTRs and spoligotypes for tree building, these isolates are localized among those classified as T and H.

**The T family.** Among the isolates initially classified as T by spoligotyping, only 26.3% (20/76) was confirmed by MIRU-VNTR 24 loci (TPR: 0.2631) (Table 3 and figure 1). By a Neighbor-Joining based phylogenetic tree, 31 of these isolates (40.9%) were positioned in LAM branch, 16 (21%) in the Haarlem branch and 9 (11.8%) in the S branch. The MIRU-VNTR characteristic signature for the T lineage confirmed by both techniques, was one copy of MIRU24 (95.8%), two copies each of MIRU02 (95.8%), MIRU04 (87.5%), MIRU20 (95.8%), MIRU39 (100%) and three copies each of MIRU27 (87.5%) and MTUB34 (87.5%). Unlike LAM isolates that present a single copy of MTUB30, isolates of the T family presented two or four copies of this allele.

**The Haarlem family.** Forty-seven isolates were classified as belonging to the H family representing 21.5% of all isolate, 85% of the these isolates initially classified as Haarlem by Spoligotyping were confirmed by MIRU-VNTRs 24 *loci* (TPR = 0.65, TNR = 0,88) (Table 3 and figure 1). Among these isolated 16 (34%) presented spoligotype profile compatible with T family, 3 (6.4%) with LAM family and 2 (2% each) isolates, respectively compatible with X and EAI. The MIRU-VNTRs characteristic for Haarlem are two copies of MIRU02 (100%) and MIRU39 (89.1%) and three copies of MTUB34 (100%), MIRU16 (91.3%), MIRU27 (97.8%), ETR-A (95.6%) and ETR-C (91.3%). Different from LAM isolates having a single copy of MTUB30, H isolates present either three or four copies of this MIRU. Combined analysis of the MIRU-24 and spoligotyping showed that isolates of genotype H1 (absences of spacers 26–31) have two copies of MIRU04, ETR-B and MTUB29 and three copies of MIRU16, MIRU31, ETR-A and ETR-C. Isolates of genotype H3 (absences of spacer 31) share two copies of MIRU02, ETR-C, MIRU4 and MTUB29 and three copies of MIRU27, ETR-A, ETR-B and MTUB34.
<table><caption>**Table 1.** Description of clusters containing 3 or more isolates in this study and their worldwide distribution in the SITVIT2 database (interrogation made on September 25th 2013).</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 SITs *</th><th>Distribution in countries with ≥3% of a given SITs **</th></tr></thead><tbody><tr><td>53 (T1) 77777777777760771</td><td>26 (11.93)</td><td>0.41</td><td>EURO-W 15.65, AMER-N 13.48, AMER-S 12.45, EURO-S 9.41, EURO-N 7.48, ASIA-W 7.31, AFRI-S 4.97, AFRI-E 4.65, ASIA-E 4.23, AFRI-N 3.52, EURO-E 3.26, AMER-C 3.23</td><td>USA 13.2, FXX 7.88, BRA 5.54, ITA 5.33, ZAF 4.86, TUR 3.47, AUT 3.42, CHN 3.09</td></tr><tr><td>42 (LAM9) 77777777777760771</td><td>18 (8.26)</td><td>0.54</td><td>AMER-S 29.7, AMER-N 11.85, EURO-S 11.31, EURO-W 9.46, AFRI-N 8.57, EURO-N 4.9, AMER-C 3.55, AFRI-E 3.55, AFRI-S 3.16</td><td>BRA 12.51, USA 11.85, COL 7.61, MAR 7.05, ITA 6.54, FXX 5.05, ESP 3.34, VEN 3.31, ZAF 3.16</td></tr><tr><td>50 (H3) 777777777777720771</td><td>16 (7.34)</td><td>0.47</td><td>AMER-S 18.16, EURO-W 17.54, AMER-N 17.54, EURO-S 11.53, EURO-E 5.51, EURO-N 5.45, AFRI-N 4.25, AFRI-S 4.04, CARI 3.48, ASIA-W 3.21</td><td>USA 17.51, BRA 7.52, FXX 6.87, AUT 6.07, ITA 5.43, ESP 5.43, PER 4.72, ZAF 4.04, CZE 3.66</td></tr><tr><td>SIT64 (LAM6) 77777777607560771</td><td>6 (2.75)</td><td>1.6</td><td>AMER-S 49.33, AMER-N 25.07, EURO-W 6.4, EURO-S 4.8, ASIA-W 3.2</td><td>BRA 38.4, USA 25.07, GUF 6.13, PRT 3.73</td></tr><tr><td>SIT47 (H1) 7777777777774020771</td><td>5 (2.29)</td><td>0.32</td><td>EURO-W 20.29, AMER-N 15.51, EURO-S 13.47, AMER-S 12.64, EURO-N 10.47, EURO-E 7.21, ASIA-W 4.02, AFRI-N 3.64</td><td>USA 15.19, ITA 8.23, AUT 8.04, BRA 7.34, FXX 6.64, FIN 6.0, CZE 3.77, ESP 3.57, SWE 3.38, PER 3.06</td></tr><tr><td>SIT60 (LAM4) 7777777760777031</td><td>5 (2.29)</td><td>1.17</td><td>AFRI-S 38.08, AMER-S 21.5, EURO-W 7.24, AFRI-N 7.01, AFRI-W 6.31, EURO-S 5.61, AMER-N 4.4</td><td>ZAF 38.08, BRA 14.72, FXX 5.14, USA 4.44, MAR 4.44, GMB 3.97, ITA 3.74, VEN 3.51</td></tr><tr><td>SIT73 (T) 777773777777607731</td><td>5 (2.29)</td><td>1.94</td><td>EURO-S 17.05, EURO-W 13.95, AMER-N 13.95, AMER-S 11.24, AFRI-S 10.85, AFRI-E 7.36, ASIA-E 4.65, AMER-C 3.88</td><td>ITA 14.34, USA 13.95, ZAF 10.85, FXX 7.75, BRA 6.59, CHN 4.65, MOZ 3.88, MEX 3.1</td></tr><tr><td>4 (Unknown) 000000007760771</td><td>4 (1.83)</td><td>1.08</td><td>EURO-S 16.49, AMER-S 14.32, AMER-N 12.43, EURO-W 10.81, ASIA-W 10.54, AFRI-S 8.11, AFRI-E 8.11, EURO-N 6.22, EURO-E 5.68</td><td>USA 9.73, ITA 8.92, BRA 8.92, ZAF 8.11, TUR 7.84, BGR 5.41, ETH 4.32, FXX 4.05, ALB 4.05</td></tr><tr><td>SIT17 (LAM2) 677737777777607771</td><td>4 (1.83)</td><td>0.6</td><td>AMER-S 59.58, AMER-N 17.95, CARI 8.9, EURO-S 5.58</td><td>BRA 28.36, VEN 27.0, USA 17.95, ESP 4.37, HTI 3.32, GLP 3.32</td></tr><tr><td>SIT33 (LAM3) 776177777777607771</td><td>4 (1.83)</td><td>0.35</td><td>AFRI-S 28.76, AMER-S 25.31, AMER-N 14.07, EURO-S 12.57, EURO-W 7.26, AMER-C 4.78</td><td>ZAF 28.76, USA 14.07, BRA 12.04, ESP 7.88, PER 6.02, ARG 5.04, FXX 4.69, ITA 4.07, HND 3.81</td></tr></tbody></table>
<table><caption>Table 1. Cont.</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 SITs *</th><th>Distribution in countries with ≥3% of a given SITs **</th></tr></thead><tbody><tr><td>SIT3907 (T1) 777703777760631</td><td>4 (1.83)</td><td>100</td><td>AMER-S 100</td><td>BRA 100</td></tr><tr><td>SIT306 (T1) 777777601560771</td><td>3 (1.38)</td><td>11.54</td><td>EURO-W 38.46, AMER-N 26.92, AMER-S 19.23, EURO-N 7.69, AFRI-M 7.69</td><td>USA 26.92, BEL 26.92, BRA 15.39, FXX 11.54, GBR 7.69, CAF 7.69, GUF 3.85</td></tr><tr><td>SIT828 (LAM4) 3777776077607731</td><td>3 (1.38)</td><td>9.38</td><td>AMER-S 90.63, EURO-S 3.13, AFRI-W 3.13, AFRI-S 3.13</td><td>BRA 84.38, PRY 6.25, ZAF 3.13, ITA 3.13, GNB 3.13</td></tr><tr><td>SIT1711 (LAM2) 677337607760771</td><td>3 (1.38)</td><td>37.5</td><td>AMER-S 75.0, EURO-S 12.5, AMER-C 12.5</td><td>BRA 37.5, VEN 25.0, MEX 12.5, ESP 12.5, COL 12.5</td></tr><tr><td>SIT2498 (H3) 777760077620771</td><td>3 (1.38)</td><td>27.27</td><td>AMER-S 100</td><td>BRA 100</td></tr><tr><td>SIT3908 (LAM2) 677730207760771</td><td>3 (1.38)</td><td>75</td><td>AMER-S 100</td><td>BRA 100</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 nose 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.

doi:10.1371/journal.pone.0107747.t001
<table><caption>**Table 2.** Allelic diversity of 24 MIRU-VNTRs loci on 218 *Mycobacterium tuberculosis* strains isolated from pulmonary tuberculosis patients in Rio de Janeiro, Brazil.</caption><thead><tr><th rowspan="2">MIRU-VNTR <i>loci</i></th><th colspan="4">Genotype families by Spoligotyping</th></tr><tr><th>LAM</th><th>H</th><th>T</th><th>All families</th></tr></thead><tbody><tr><td>MIRU02</td><td>0.362</td><td>0.148</td><td>0.353</td><td>0.325</td></tr><tr><td>MIRU04</td><td>0.061</td><td>0.153</td><td>0.225</td><td>0.139</td></tr><tr><td>MIRU10</td><td>0.56</td><td>0.687</td><td>0.692</td><td>0.692</td></tr><tr><td>MIRU16</td><td>0.636</td><td>0.243</td><td>0.439</td><td>0.544</td></tr><tr><td>MIRU20</td><td>0.081</td><td>0.499</td><td>0.24</td><td>0.27</td></tr><tr><td>MIRU23</td><td>0.541</td><td>0.529</td><td>0.584</td><td>0.688</td></tr><tr><td>MIRU24</td><td>0.021</td><td>0.101</td><td>0.081</td><td>0.071</td></tr><tr><td>MIRU26</td><td>0.723</td><td>0.559</td><td>0.687</td><td>0.639</td></tr><tr><td>MIRU27</td><td>0.314</td><td>0.145</td><td>0.105</td><td>0.227</td></tr><tr><td>MIRU31</td><td>0.429</td><td>0.357</td><td>0.51</td><td>0.445</td></tr><tr><td>MIRU39</td><td>0.138</td><td>0.251</td><td>0.027</td><td>0.124</td></tr><tr><td>MIRU40</td><td>0.517</td><td>0.623</td><td>0.76</td><td>0.689</td></tr><tr><td>ETR-A</td><td>0.14</td><td>0.417</td><td>0.625</td><td>0.537</td></tr><tr><td>ETR-B</td><td>0.407</td><td>0.237</td><td>0.396</td><td>0.389</td></tr><tr><td>ETR-C</td><td>0.357</td><td>0.447</td><td>0.573</td><td>0.541</td></tr><tr><td>MTUB04</td><td>0.612</td><td>0.53</td><td>0.68</td><td>0.716</td></tr><tr><td>MTUB21</td><td>0.436</td><td>0.587</td><td>0.673</td><td>0.624</td></tr><tr><td>MTUB29</td><td>0.253</td><td>0.432</td><td>0.266</td><td>0.397</td></tr><tr><td>MTUB30</td><td>0.269</td><td>0.494</td><td>0.652</td><td>0.591</td></tr><tr><td>MTUB34</td><td>0.617</td><td>0.305</td><td>0.617</td><td>0.452</td></tr><tr><td>MTUB39</td><td>0.391</td><td>0.573</td><td>0.391</td><td>0.577</td></tr><tr><td>QUB11</td><td>0.697</td><td>0.702</td><td>0.697</td><td>0.737</td></tr><tr><td>QUB26</td><td>0.384</td><td>0.63</td><td>0.384</td><td>0.522</td></tr><tr><td>QUB4156</td><td>0.842</td><td>0.842</td><td>0.842</td><td>0.84</td></tr></tbody></table>

doi:10.1371/journal.pone.0107747.t002

**The S family.** Eleven isolates were classified as belonging to the S family and four of these exhibit SIT4, that in SpolDB4 was classified as LAM3-S Convergent, having absence of spacers 1–24 and 33–34 (The characteristic pattern of S family is absence of spacers 9–10, and 33–34). Among the others isolates, we observed SITs 53/n = 1, 102/n = 1, 378/n = 1, 2500/n = 1, 3907/n = 2 and 3909/n = 1; only two of these had the characteristic of S family (SIT3909/LAM8 and SIT2500/Unknown by SITVIT2 and H1 by Spotclust). Isolates of the S family shared their copy number in six loci, being one copy of MIRU24 and MTUB 21, two copies

each of MIRU20 and MIRU39 and three copies MIRU10 and MIRU27.

### RD<sup>Rio</sup>, RD174 and *fbpC*<sup>103</sup> analysis

The genotypes defined by the presence of RD<sup>Rio</sup>, RD174 and SNP *fbpC*<sup>103</sup> were added to the classification based on 24-MIRU-VNTR typing (Table S2). Thirty-six isolates (16.5%) were excluded from the final analysis either because of showing genotypes suggestive for multiple infection of because of failure in at least one of the three genotype procedures. The results of 182

<table><caption>**Table 3.** Frequencies of strains according to classification by Spoligotyping (SITVIT2) and 24 loci MIRU-VNTR (MIRUVNTRPlus).</caption><thead><tr><th>Lineage</th><th colspan="2">Spoligotyping (%)</th><th colspan="2">24 loci MIRU-VNTR (%)</th></tr></thead><tbody><tr><td>LAM</td><td>95</td><td>(43.58)</td><td>136</td><td>(62.39)</td></tr><tr><td>T</td><td>76</td><td>(34.86)</td><td>21</td><td>(9.63)</td></tr><tr><td>H</td><td>40</td><td>(18.35)</td><td>47</td><td>(21.56)</td></tr><tr><td>S</td><td>0</td><td>0</td><td>11</td><td>(5.05)</td></tr><tr><td>X</td><td>1</td><td>(0.46)</td><td>1</td><td>(0.46)</td></tr><tr><td>EAI</td><td>5</td><td>(2.29)</td><td>1</td><td>(0.46)</td></tr><tr><td>BEIJING</td><td>1</td><td>(0.46)</td><td>1</td><td>(0.46)</td></tr></tbody></table>

doi:10.1371/journal.pone.0107747.t003
<figure><img src="image_3.png" /></figure>

<table><thead><tr><th colspan="2">MIRU-VNTR 24 loci vs Spoligotyping</th></tr></thead><tbody><tr><td>Positive Predictive Value (PPV) or Precision</td><td>0.7433</td></tr><tr><td>Accuracy (AAC)</td><td>0.6376</td></tr><tr><td>True positive rate (TPR) or sencivity</td><td>0.8176</td></tr><tr><td>Erro Rate</td><td>0.3623</td></tr></tbody></table>

<table><caption>**Figure 1.** Confusion Matrix comparing the classifications obtained by Spoligotyping and MIRU-VNTR. <fnref n="1" /> Classification according to SITIVIT2. For unknown Spoligotypes, we used SpotClust. <fnref n="2" /> Patterns were classified based on a VNTRplus database that allows construction of a Neighbor-Joining based phylogenetic tree, visualizing proximity of a particular genotype with that of a set of reference strains to the genotype family level.
doi:10.1371/journal.pone.0107747.g001</caption><thead><tr><th></th><th>LAM</th><th>T</th><th>Haarlen</th></tr></thead><tbody><tr><td>Positive Predictive Value (PPV) or Precision</td><td>0.6691</td><td>0.9523</td><td>0.5531</td></tr><tr><td>Accuracy (AAC)</td><td>0.7752</td><td>0.7385</td><td>0.8394</td></tr><tr><td>True positive rate (TPR) or sencivity</td><td>0.9578</td><td>0.2631</td><td>0.6500</td></tr><tr><td>True Negative Rate (TNR) or Specificity</td><td>0.4193</td><td>0.9929</td><td>0.8820</td></tr><tr><td>Negative Predictive Value (NPV)</td><td>0.9512</td><td>0.7157</td><td>0.9181</td></tr><tr><td>Positive Predictive Value (PPV)</td><td>0.669</td><td>0.9523</td><td>0.5531</td></tr><tr><td>Erro Rate</td><td>0.2247</td><td>0.2614</td><td>0.3500</td></tr></tbody></table>

strains are summarized in Table 3 and for simplification of interpretation, we defined three groups, being LAM (n = 77), LAM-like (n = 33) and non-LAM (n = 72) as determined by spoligotyping and MIRU-VNTR typing.

The SNP $fbpC^{103}$ was observed in 107 (58.8%) of the isolates, including 97.3% (74/77) of LAM, 69.6% (23/29) of “LAM-like” and 13.9% (10/72) of non-LAM. The total frequency of RD$^{Rio}$ was 11% (20/182), with 20.7% (16/77) among LAM, 6.0% (2/29) among “LAM-like” and 2.7% (2/72) among non-LAM isolates. The overall frequency of RD174 was 15.3% (28/182), being 33.7% (25/77) in the LAM genotype and 6.8% (2/29) in “LAM-like”. All LAM/RD$^{Rio}$ isolates presented the SNP $fbpC^{103}$ (Table 4 and Table 5), but not all were RD174; two isolates classified as LAM2 SIT3908 were not had RD174 (isolated from different patients). The frequency of RD174 (n = 25) was therefore higher than that of RD$^{Rio}$ (32.5% vs 20.7%; p = 0.053 and $X^2 = 2.69$) in LAM. Eleven LAM isolates (14.3%) presented RD174 but not RD$^{Rio}$, all were positive for the $fbpC^{103}$.

The allelic diversity of the 24 loci MIRU-VNTR loci in LAM/RDRio isolates (Figure S2) showed that the copy number of

combining of MIRU04, MIRU20, MIRU24, MIRU31, ETR-A, MTUB21 and MTUB30 loci was a signature of this genotype (2213231) (Table 4). In general, these loci present low variability in LAM, except for MTUB21, is being moderately variable in such isolates (table 2). Upon comparing 24 MIRU-VNTR signatures of LAM, LAM-like and non-LAM strains (Table 4), we observed two isolates that presented the hypothetical ancestral MIRU-VNTR signature (224226153321) for RD$^{Rio}$ that was suggested by Lazzarini et al. (2007) [2]. One isolate (C2009) was classified as LAM1 SIT2539 and the other (C1966) as H3 SIT50 by spoligotyping and “LAM-like” by 24 MIRU-VNTR; both presented SNP $fbpC^{103}$ and were deleted for RD$^{Rio}$ and RD174. Among the LAM/RD$^{Rio}$ strains, frequency of LAM subtypes as defined by spoligotyping using SITVIT2 was 31.5% of LAM2, 25% of LAM9, 12.5% each of LAM6 and LAM5 and 6.3% each of LAM4 and LAM1. Figure 2 is a graphical representation of these three markers in the LAM strains as defined by 24 loci MIRU-VNTRs.

<footnote n="1">¹ Classification according to the SITIVIT2. For unknown Spoligotypes, we used the SpotClust.</footnote>

<footnote n="2">² The patterns were classified based on MIRU-VNTRPlus tool that allows construction of a Neighbor-Joining based phylogenetic tree, visualizing proximity of a particular genotype with that of a set of reference strains to the genotype family level.</footnote>
<table><caption>Table 4. Cont.</caption><thead><tr><th colspan="11">24 MIRU-VNTRs</th></tr><tr><th colspan="2">MIRU</th><th colspan="2">ETR</th><th colspan="2">MTUB</th><th colspan="2">QUB</th></tr><tr><th colspan="2">12 MIRU-VNTRs**</th><th colspan="2">Additional 12 MIRU-VNTRs</th><th colspan="2"></th><th colspan="2"></th></tr><tr><th>Strain</th><th>Year</th><th>SITVIT2</th><th>SIT2</th><th>Lineage*</th><th>RD<sup>Rio</sup></th><th>RD174</th><th>fbpC<sup>103</sup></th><th></th><th></th><th></th><th></th></tr></thead><tbody><tr><td>C2142A2</td><td>2009</td><td>H3</td><td>50</td><td>H</td><td>+</td><td>-</td><td>-</td><td>3</td><td>2</td><td>4</td><td>3</td><td>1</td></tr><tr><td>C5817</td><td>2008</td><td></td><td></td><td>Unknown</td><td>4</td><td>5</td><td></td><td>3</td><td>3</td><td>4</td><td>2</td><td>2</td></tr></tbody></table>

*Classification based on construction of phylogenetic tree using the N-J Algorithm and evaluating proximity of particular isolate with a set of 182 reference strains in the MIRU-VNTRPlus. Note: Strains with RD<sup>Rio</sup> (LAM, Like-LAM and Non LAM).

**In bold the same number of copy of hypothetic Ancestral RD<sup>Rio</sup> as suggested by Iazzarini et al, 2007 [2] and n the probable number of copy of hypothetic Ancestral RD<sup>Rio</sup> in the additional 12 loci.
doi:10.1371/journal.pone.0107747.t004

## Discussion

Rio de Janeiro is the capital of the state of Rio de Janeiro, located in the southeast of Brazil, has a population of 16 million habitants, six and a half million of these living in the capital, being the second largest city and a major touristic attraction in Brazil (Census 2010 http://cidades.ibge.gov.br/xtras/perfil.php?lang=&codmun=330455 accessed 12/18/2013). According to the Ministry of Health, 11.639 new TB cases were recorded in 2011 in the state with an incidence of 72.3/100.000 habitants and the highest mortality rate (5.1/100.000 habitants) at the national level [1]. Rio de Janeiro city has strong economic and social contrasts, with a large portion of the population living in numerous suburbs, consisting of “communities”, urban areas where housing conditions, health, education and security are extremely precarious. These factors are directly related to the number of TB cases detected.

In the present study, we performed spoligotyping and 24 MIRU-VNTR typing and characterized a SNP in *fbpC*<sup>103</sup> and the status of RD174 and RD<sup>Rio</sup>, to decipher the population structure and the phylogenetic relationships of the MTBC strains of TB cases in this particular population that attended a single reference center in the city of Rio de Janeiro. *Mycobacterium tuberculosis* is classified into six phylogenetic lineages, each of which can be divided into sublineages [32]. Members of the Euro-American lineages represented by the LAM, H, T, S and X Spoligotyping families, which are genetically closely related are the most common in South America [19], [20], [21] and this was confirmed recently also in Brazil [22] and in different states of the country [8], [43], [44], [45], [46], [47] [48], but differences in the frequency of these families are observed, these difference could also be due to differences in population and immigration history in each region or in period associated genotype frequencies or in differences in sample size. We here confirm the predominance of the Euro-American lineage (also known as lineage 4) with only two isolates classified as Indo-Oceanic lineage (EAI5) and East Asian lineage (Beijing), the first being more prevalent in east Africa, southeast Asia and in south India and the second in east Asia, Russia and South Africa. [32]. Recently, the frequency of such strains in Brazil was described as being less than 1% except for the higher frequency in the state of Pará, north Brazil [22].

The most prevalent families as defined by spoligotyping were LAM (43.6%), T (34.9%) and H (18.3%); however, classification based on 24-MIRU-VNTR and Neighbor-Joining based phylogenetic tree building using the database tool (www.miruvntrplus.org) showed considerable difference with that obtained by spoligotyping, being respectively 62.4%, 9.6% and 21.6% (error rate of 0.36). Among the isolates with discordant results, mostly isolates initially classified as T (n = 31), H (n = 11) and EAI (n = 3) by spoligotyping were classified as LAM by 24 MIRU-VNTRs typing, for isolates that did not show LAM prototype (absence of spacers 21–24 and 33–36), conveniently named LAM-like. Subsequently, 78.3% these isolates confirmed to be LAM by the presence of the LAM-specific SNP *fbpC*<sup>103</sup> [7], [38] and when considering the presence of this SNP as an absolute marker for LAM family the frequency of this lineage is 58.9%. We also have preliminary data showing that the supposed LAM-specific marker *ligB*<sup>404</sup> was observed in isolates that had been defined by spoligotyping as being T or H (data not shown). Very recently, Mokrousov et al. reported difference in classification of LAM strains as defined by spoligotyping and SNP analysis (*fbpC*<sup>103</sup> and *ligB*<sup>404</sup>) [49] although both SNPs were previously validated in different collection of clinical isolates and reference strains as to be specific for the LAM family [7], [38].
<table><caption>**Table 5.** The frequency of *fbpC*<sup>103</sup>, RD<sup>Rio</sup> and RD174 in LAM, LAM-like and Non LAM isolates, as designated by spoligotyping and 24 MIRU-VNTR typing.</caption><thead><tr><th>SNP and RD</th><th>LAM (%)<br>(n=77)</th><th>LAM-Like (%)<br>(n=29)</th><th>Non LAM (%)<br>(n=75)</th><th>Total Frequency (%)<br>(n=182)</th></tr></thead><tbody><tr><td>fbpC<sup>103</sup></td><td>74 (97.3)</td><td>23 (79.3)</td><td>10 (13.3%)</td><td>107 (58.8)</td></tr><tr><td>RD<sup>Rio</sup></td><td>16 (20.7)</td><td>2 (6.9%)</td><td>2 (2.6%)</td><td>20 (11)</td></tr><tr><td>RD174</td><td>25 (32.5)</td><td>2 (6.9%)</td><td>1 (1.3%)</td><td>28 (15.4)</td></tr></tbody></table>

doi:10.1371/journal.pone.0107747.t005

It is now common knowledge that spoligotyping has limitations as a tool prediction of the exact phylogenetic relationships between strains of the MTBC, particularly in modern strains (Euro American, East Asian and Indian e East African) [16], [18], [24], [48], [49], [50], [51] [52], mainly due to homoplasy [16], [24]. The accuracy of the phylogenetic grouping by MIRU-VNTR is more exact than that of spoligotyping but depends on the number of loci included in the analysis [49] and classification errors are reduced when analyzing 24 loci [16]. Indeed, several authors have suggested that SNPs are more suitable than spoligotyping and MIRU-VNTR settings for phylogenetic classification [16], [38], [53], [54], [55], [56]. The *fbpC*<sup>103</sup> is described as a good marker for the LAM family [7], [38], [49], [57], [58], [59], and our original intention was to evaluate if strains RD<sup>Rio</sup> was still prevalent in Rio de Janeiro as seen previously by Lazzarini in clinical isolates collected between 2002 and 2003 [2], for this

purpose we believed only a marker that could differentiate LAM and non-LAM associated with the 24 MIRU-VNTR and Spoligotyping could be sufficient, however the scenario was observed more complex.

We here present the first data comparing classification by spoligotyping and MIRU-VNTRs and *fbpC*<sup>103</sup> in Euro American lineage prevalent in Rio de Janeiro and Brazil. Among these different types of markers, we observed four major groups: (i) strains classified as LAM by spoligotyping and MIRU-VNTRs 24 loci without the LAM-characteristic SNP *fbpC*<sup>103</sup>, (ii) strains not classified as LAM by spoligotyping and MIRU-VNTRs 24 loci but carrying the *fbpC*<sup>103</sup> SNP, (iii) strains classified by spoligotyping as non-LAM but as LAM by MIRU-VNTRs 24 loci (LAM-like) and with SNP *fbpC*<sup>103</sup> and (iv), strains classified by spoligotyping as non-LAM and LAM by MIRU-VNTRs 24 loci (LAM-like) but not presenting the SNP *fbpC*<sup>103</sup>. These different scenarios could be explained by convergent evolution of spoligotypes and of MIRU-VNTRs loci (even including 24 alleles) because a limited number of loci were evaluated that might evolve rapidly and therefore susceptible to pronounced convergence [53] and/or because the existence of ancestral progenitor of Euro-American lineages containing SNP *fbpC*<sup>103</sup>. In addition, a possible important limitation of the current classification based on MIRU-VNTRs and their similarity with genotypes present in MIRU-VNTRplus is that, despite including well characterized strains, the database contains a limited number of strains that does not reflect the real genetic diversity of isolates belonging to the MTBC; another limitation of this study is that no additional specific SNPs for H, T, S and X and T were investigated. The Whole Genome Sequencing (WGS) is superior to conventional genotyping for MTBC [60] and has been used in areas not yet studied, from global (phylogeography) for site (transmission chains and diversity of circulating strains), for single patient (clonal diversity) and the bacterium itself (evolutionary studies) [61]. We intend to compare through WGS (developing) these isolates to a better comprehension of the evolution of lineages Euro American, such as the development of new methodologies that allow a more rapid and accurate typing.

Analyzing the data obtained in this study, the spoligotypes of the T family showed the largest number of divergent results when compared to classification by 24-MIRU (TPR = 0,26/ TNR = 0,99). The prototype of the T family is characterized by the absence of spacers 33–36 only [48] and have been observed in almost every country, representing 20% of all isolates deposited in the database SITVITWEB. Despite the high frequency, this genotype family is still considered as “ill-defined” and includes non-monophyletic groups [54] [55]. In South America, the frequency of this family is 26.7% and in Brazil 18.6% (370/ 1991), with the T1 SIT53 subfamily being the most prevalent SIT [22], as observed also in this study. Here, among the 76 isolates classified as T by spoligotyping, only 21 (27.6%) were confirmed

<figure><img src="image_4.png" /><figcaption>**Figure 2.** Venn diagram illustrating the different markers in isolated classified as LAM by 24 loci MIRU-VNTR. Notes: The Venn diagram was constructed based on LAM isolates defined by 24 loci MIRU-VNTR. The sizes of the circles is not proportional to the real frequency of these markers. Twenty isolates classified with LAM by MIRU-VNTRs were removed from the final analysis because of showing a mixed genotype or failure in at least one of the three genotype procedures. * absence of spacers 33–36 in the spoligotyping profile, ** absence of spacers 21–24 in the spoligotyping profile, ***Two isolates RD<sup>Rio</sup> but not RD174 (LAM2 SIT3908). doi:10.1371/journal.pone.0107747.g002</figcaption></figure>
by 24 MIRU-VNTRs typing, the rest was reclassified as LAM (n = 31), S (n = 9) and H (n = 16). Interestingly, SIT53 was associated with mixed infection in a study conducted in South Africa, a country characterized by a high prevalence of TB [57] and Lazzarini *et al.* [62], using a computational approach, showed that this is the most frequent false spoligotype derived from mixed infections. However, among isolates with this SIT, we did not observe double signals during 24 MIRU typing indicative for mixed infections. We also observed that spoligotypes, indicative for the T family, were sometimes grouped with spoligotypes of the H family by MIRU-VNTR typing. This could be related with the fact that the prototype spoligotype defining T1 SIT53 and H3 SIT 50 differ only in the absence of the spacer 31. We also observed that in our population, the absence of spacer 31 is not crucial for classification as being either H or T; what differentiates between the two is the number of copies of ETR-A, ETR-B, ETR-C. The H3 subfamily is characterized by the combination 323 of these alleles while T strains present considerable diversity of these loci (one to three copies of ETR-A, two or three copies of ETR-B and three to five copies of ETR-C).

The RD nominated RD<sup>Rio</sup> was first reported as new *M. tuberculosis* lineage in Rio de Janeiro in 2007, Lazzarini *et al.* [2] and is a deletion of 26.3 kb restricted to the LAM family and in particular, in subfamilies LAM1, LAM2, LAM4, LAM5, LAM6 and LAM9. This deletion affects two genes, including two genes encoding Proline-proline-glutamic Acid Proteins (PPE) [2] and association between RD<sup>Rio</sup> strains and high prevalence may be related to virulence and/or adaptation specifies the Latin American and European population-based epidemiological and clinical characteristics; however, studies have proven to be inconclusive or contradictory [4], [63], [64]. The lineage RD<sup>Rio</sup>, was identified in different countries [5], [6], [7], [49] and in Brazil, has been described in different states, including Rio de Janeiro, Rio Grande do Sul, Minas Gerais, Espírito Santo and Rio Grande do Sul [3], [8], [63], [64]. The frequency ranges from 30 to 38% of isolates tested and was associated with MRD-TB [4], [63], [64] and with genotype clustering [5], [63], indicating a higher rate of recent infection and transmissibility. Another deletion, RD174, initially described as a marker for the LAM family and as a co- marker for RD<sup>Rio</sup> [7], [8] was associated with high transmissibility [10].

The frequency of strains RD<sup>Rio</sup> in the present study, with an overall frequency of 11% and 20.7% in LAM, is lower than that observed in other studies and even when including the eight isolates that showed mixed RD<sup>Rio</sup>/WT signals, resulting frequencies of 14.7% and 28.2% still lower than previously observed. This difference could be related with the relative low frequency of LAM1, LAM2, LAM5, LAM6 and LAM9 related subfamilies in our study sample. Earlier studies on classification of RD<sup>Rio</sup> strains was based on genotypes defined by spoligotyping and/or 12 MIRUs only [2], [5], [6], [7], [8], [63], [64] and again, this is the first study that used 24 MIRU-VNTR typing to conduct a more detailed phylogenetic analysis. We verified the signature of MIRU-VNTR loci for RD<sup>Rio</sup> and RD174, as compared with that of the hypothetical ancestral of RD<sup>Rio</sup> as defined by 12 MIRUs [2] and that, besides sharing two copies of MIRU04 and MIRU20, they also share one copy of MIRU24, three copies of MIRU31, two copies of ETR-A, three copies of MTUB21 and one copy of MTUB30, yielding 2213231 as fingerprint for RD<sup>Rio</sup>. All LAM/RD174 isolates, with or without RD<sup>Rio</sup>, carried two copies of MIRU20 and ETRA and one copy of MTUB31. We observed that LAM3 has two copies of MTUB30 (2213232) and we propose

that this subfamily that do not carry RD<sup>Rio</sup>, has evolved independently. We also observed RD<sup>Rio</sup> in two isolates (2.6%) with a spoligotype not indicative for being LAM; a small number of such cases have also been observed by other research groups [5], [7], [64]. One of these isolates that had been classified by spoligotyping as H3 SIT53 was reclassified by 24 MIRU as being "LAM-Like" and called attention because it had the MIRU signature of the hypothetical ancestor of RD<sup>Rio</sup> and carried the SNP *fbpC*<sup>103</sup> and RD174. Overall, we observed four scenarios: (i) isolates with the RD<sup>Rio</sup> and RD174, (ii) isolates showing only RD<sup>Rio</sup> (iii) isolates showing only RD174 and (iv) isolates that did not carry any of the deletions, suggesting that both markers evolved on different time points. This is different from earlier data [7] that claim that RD174 is an absolute co-marker for RD<sup>Rio</sup> and therefore, studies that use the presence of RD174 to infer the presence of RD<sup>Rio</sup> [8] may overestimate the frequency of the latter. In 2007, Lazzarini *et al.* [2] suggested that RD<sup>Rio</sup> arose by homologous recombination between genes and although all neighboring sequences were identical, such event could have happened more than once. This suggests that the RD174 deletion occurred before RD<sup>Rio</sup> but this needs to be confirmed as we also observed RD<sup>Rio</sup> strains that had intact RD174. In Figure S3, we propose a possible evolution of members of the Euro-American lineages but with the limitation that the spoligotype defined lineages S, X and other are not represented and this concerns a sampling only from Rio de Janeiro.

In a recent study, Hill *et al.* (2012) [53], mentions the difficulty in studying the evolution of the Euro-American lineage (LAM, Haarlem, T, X and S) using spoligotyping due to the large number of IS*6110* copies in such strains that may result in IS*6110* mediated deletions in the DR locus. This might be the reason why bacterial evolution exclusively based on spoligotyping is not robust and the wide range of profiles reported as unclassified in SITVITWEB. Our approach, combining spoligotyping, MIRU-VNTRs, SNPs and RDs allowed the reclassification of 13 SITs that did not rank in SITVITWEB, allowing definition of 29 new spoligotypes and refine classification of isolates belonging to Euro-American lineage. Our data are also support the idea that absence of spacers 21–24 is not sufficient for classification as LAM and of spacer 31 to differentiate T and H, the latter indicative for subfamilies H3 SIT50 and T1 SIT53. Possible explanations are that ancestral lineages are currently circulating (plesiomorphic state) or that the isolates are suffering homoplasic evolution and reversion into the plesiomorphic state.

## Supporting Information

**Figure S1** Dendrogram constructed with BioNumerics software version 6.6 with the results of MIRU-VNTRs 24 loci and spoligotyping by similarity coefficient for categorical data and the neighbor-joining algorithm. 1st column (after spoligotypes): number of isolated label; 2nd column: patients who have more than one isolate in the study (n = 27) received numbering 1–27, and the different strains present the same numbering; 3rd column: International Spolotype Types (SIT); 4th column: classification obtained through SITVITWEB (family and subfamily).
(DOC)

**Figure S2** Allelic diversity of the 24 loci MIRU-VNTR loci in LAM/RDRio isolates.
(DOC)
**Figure S3** Possible evolution of *M. tuberculosis* lineage Euro-American (LAM, T and H) according to the markers analyzed in this study. The Euro-American Ancestral evolved into two distinct groups: LAM ancestral and T/H ancestral, characterized by absences of spacers 33–36 and one copy of MIRU24 that is common to all Euro-American lineages. The ancestral LAM strains have *fbpC*¹⁰³ and this was the basis of LAM A (LAM9), with additional absence of spacers 21–24, one copy of MTUB30 and two copies of MIRU04 and ETR-A. LAM A on its turn is the basis for two other groups: LAM B (LAM9-LAM4, LAM1-LAM2-LAM5 and LAM6) and LAM C (LAM3). The LAM B evolved from LAM B1, characterized by a deleted RD174 and on its turn to LAM B2, with both deleted RD174 and RD*Rio*. The H/T Ancestral lineage is the origin of both groups H and T (difference only in MIRU-VNTR copies), showing absence of spacers 33–36; additional loss of spacer 31 led to subtype H A, observed in high frequency in this study. (DOC)

<caption>**Table S1** Detailed genotyping results and associated demographic, epidemiologic and nomenclature information on 218 *Mycobacterium tuberculosis* strains</caption>

**References**

1. Secretaria de Vigilância em Saúde – Ministério da Saúde, Especial Tuberculose (2012) Boletim Epidemiológico, vol. 43.

2. Lazzarini LC, Huard RC, Boechat NL, 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. J Clin Microbiol 45: 3891-3902.

3. Lazzarini LC, Spindola SM, Bang H, Gibson AL, Weisenberg S, et al. (2008) RDRio Mycobacterium tuberculosis infection is associated with a higher frequency of cavitary pulmonary disease. J Clin Microbiol 46: 2175-2183.

4. Dalla Costa ER, Lazzarini LC, Perizzolo PF, Diaz CA, Spies FS, et al. (2013) 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 51: 1071-1077.

5. Weisenberg SA, Gibson AL, Huard RC, Kurepina N, Bang H, et al. (2012) Distinct clinical and epidemiological features of tuberculosis in New York City caused by the RD(Rio) Mycobacterium tuberculosis sublineage. Infect Genet Evol 12: 664-670.

6. David S, Duarte EL, Leite CQ, Ribeiro JN, Maio JN, et al. (2012) Implication of the RD(Rio) Mycobacterium tuberculosis sublineage in multidrug resistant tuberculosis in Portugal. Infect Genet Evol 12: 1362-1367.

7. Gibson AL, Huard RC, Gey van Pittius NC, Lazzarini LC, Driscoll J, et al. (2008) Application of sensitive and specific molecular methods to uncover global dissemination of the major RDRio Sublineage of the Latin American-Mediterranean Mycobacterium tuberculosis spoligotype family. J Clin Microbiol 46: 1259-1267.

8. Cardoso Oelemann M, Gomes HM, Willery E, Possuelo L, Batista Lima KV, et al. (2011) The forest behind the tree: phylogenetic exploration of a dominant Mycobacterium tuberculosis strain lineage from a high tuberculosis burden country. PLoS One 6: e18256.

9. Rindi L, Lari N, Garzelli C (2012) Large Sequence Polymorphisms of the Euro-American lineage of Mycobacterium tuberculosis: a phylogenetic reconstruction and evidence for convergent evolution in the DR locus. Infect Genet Evol 12: 1551-1557.

10. de Jong BC, Antonio M, Awine T, Ogungbemi K, de Jong YP, et al. (2009) Use of spoligotyping and large sequence polymorphisms to study the population structure of the Mycobacterium tuberculosis complex in a cohort study of consecutive smear-positive tuberculosis cases in The Gambia. J Clin Microbiol 47: 994-1001.

11. 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.

12. 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. J Clin Microbiol 35: 907-914.

13. Supply P, Mazars E, Lesjean S, Vincent V, Gicquel B, et al. (2000) Variable human minisatellite-like regions in the Mycobacterium tuberculosis genome. Mol Microbiol 36: 762-771.

14. Chacón-Salinas R, Serafin-López J, Ramos-Payán R, Méndez-Aragón P, Hernández-Pando R, et al. (2005) Differential pattern of cytokine expression by macrophages infected in vitro with different Mycobacterium tuberculosis genotypes. Clin Exp Immunol 140: 443-449.

isolated from pulmonary tuberculosis patients in Rio de Janeiro, Brazil.
(XLS)

<caption>**Table S2** Detailed genotyping results of *fbpC*¹⁰³, RD174 and RD*Rio* on 182 *M. tuberculosis* strains isolated from pulmonary tuberculosis patients in Rio de Janeiro, Brazil.
(XLS)</caption>

**Acknowledgments**

We thank the Genomic Platform of Network Technology Platforms (PDTIS) Institute Oswaldo Cruz/FIOCRUZ.RPT01D - Analysis of Fragments Oswaldo Cruz Institute in Rio de Janeiro.

**Author Contributions**

Conceived and designed the experiments: SEGV HMG ARS NB PNS. Performed the experiments: SEGV HMG ARS NB PNS. Analyzed the data: SEGV CCA HMG LLG ECC MIA DC MKG. Contributed reagents/materials/analysis tools: PNS KVL MLL FT PCSC NR CS. Wrote the paper: SEGV CCA PNS NR CS.

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

16. Cowan LS, Diem L, Monson T, Wand P, Temporado D, et al. (2005) Evaluation of a two-step approach for large-scale, prospective genotyping of *Mycobacterium tuberculosis* isolates in the United States. J Clin Microbiol 43: 688-695.

17. Dale JW, Nor RM, Ramayah S, Tang TH, Zainuddin ZF (1999) Molecular epidemiology of tuberculosis in Malaysia. J Clin Microbiol 37: 1265-1268.

18. Kato-Meda M, Gagneux S, Flores LL, Kim EY, Small PM, et al. (2011) Strain classification of *Mycobacterium tuberculosis*: congruence between large sequence polymorphisms and spoligotypes. Int J Tuberc Lung Dis 15: 131-133.

19. Abadia E, Sequera M, Ortega D, Méndez 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.

20. Cerezo I, Jiménez 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 Bogotá, Colombia. Infect Genet Evol 12: 657-663.

21. Gonzalo X, Ambroggi M, Cordova E, Brown T, Poggi S, et al. (2011) Molecular epidemiology of *Mycobacterium tuberculosis*, Buenos Aires, Argentina. Emerg Infect Dis 17: 528-531.

22. 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.

23. 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.

24. Brudey K, Driscoll JR, Rigouts L, Prodinger WM, Gori A, et al. (2006) Mycobacterium tuberculosis complex genetic diversity: mining the fourth international spoligotyping database (SpolDB4) for classification, population genetics and epidemiology. BMC Microbiol 6: 23.

25. Filliol I, Driscoll JR, Van Soolingen D, Kreiswirth BN, Kremer K, et al. (2002) Global distribution of *Mycobacterium tuberculosis* spoligotypes. Emerg Infect Dis 8: 1347-1349.

26. Sola C, Filliol I, Gutierrez MC, Mokrousov I, Vincent V, et al. (2001) Spoligotype database of *Mycobacterium tuberculosis*: biogeographic distribution of shared types and epidemiologic and phylogenetic perspectives. Emerg Infect Dis 7: 390-396

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

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

29. Allix-Béguec C, Fauville-Dufaux M, Supply P (2008) Three-year population-based evaluation of standardized mycobacterial interspersed repetitive-unit-variable-number tandem-repeat typing of *Mycobacterium tuberculosis*. J Clin Microbiol 46: 1398-1406.
30. García de Viedma D, Mokrousov I, Rastogi N (2011) Innovations in the molecular epidemiology of tuberculosis. Enferm Infecc Microbiol Clin 29 Suppl 1: 8-13.

31. Goh KS, Rastogi N, Berchel M, Huard RC, Sola C (2005) Molecular evolutionary history of tubercle bacilli assessed by study of the polymorphic nucleotide within the nitrate reductase (narGHJI) operon promoter. J Clin Microbiol 43: 4010-4014.

32. Gagneux S, Small PM (2007) Global phylogeography of Mycobacterium tuberculosis and implications for tuberculosis product development. Lancet Infect Dis 7: 328-337.

33. Huard RC, Fabre M, de Haas P, Lazzarini LC, van Soolingen D, et al. (2006) Novel genetic polymorphisms that further delineate the phylogeny of the Mycobacterium tuberculosis complex. J Bacteriol 188: 4271-4287.

34. Vasconcellos SE, Huard RC, Niemann S, Kremer K, Santos AR, et al. (2010) Distinct genotypic profiles of the two major clades of Mycobacterium africanum. BMC Infect Dis 10: 80.

35. Stucki D, Gagneux S (2013) Single nucleotide polymorphisms in Mycobacterium tuberculosis and the need for a curated database. Tuberculosis (Edinb) 93: 30-39.

36. Musser JM, Amin A, Ramaswamy S (2000) Negligible genetic diversity of mycobacterium tuberculosis host immune system protein targets: evidence of limited selective pressure. Genetics 155: 7-16.

37. Dos Vultos T, Mestre O, Rauzier J, Golec M, Rastogi N, et al. (2008) Evolution and diversity of clonal bacteria: the paradigm of Mycobacterium tuberculosis. PLoS One 3: e1538.

38. Homolka S, Projahn M, Feuerriegel S, Ubben T, Diel R, et al. (2012) High resolution discrimination of clinical Mycobacterium tuberculosis complex strains based on single nucleotide polymorphisms. PLoS One 7: e39855.

39. van Embeden JD, Cave MD, Crawford JT, Dale JW, Eisenach KD, et al. (1993) Strain identification of Mycobacterium tuberculosis by DNA fingerprinting: recommendations for a standardized methodology. J Clin Microbiol 31: 406-409.

40. Zhang J, Abadia E, Refregier G, Tafaj S, Boschiroli ML, et al. (2010) Mycobacterium tuberculosis complex CRISPR genotyping: improving efficiency, throughput and discriminative power of 'spoligotyping' with new spacers and a microbead-based hybridization assay. J Med Microbiol 59: 285-294.

41. Graur D, W-H Li (2000) Dynamics of genes in populations, p. 58. In D. Graur and W.-H. Li (ed.), Fundamentals of molecular evolution. Sinauer Associates, Sunderland, Mass.

42. Hunter PR, Gaston MA (1988) Numerical index of the discriminatory ability of typing systems: an application of Simpson's index of diversity. J Clin Microbiol 26: 2465-2466.

43. Noguti EN, Leite CQ, Malaspina AC, Santos AC, Hirata RD, et al. (2010) Genotyping of Mycobacterium tuberculosis isolates from a low-endemic setting in northwestern state of Paraná in Southern Brazil. Mem Inst Oswaldo Cruz 105: 779-785.

44. Mendes NH, Melo FA, Santos AC, Pandolfi JR, Almeida EA, et al. (2011) Characterization of the genetic diversity of Mycobacterium tuberculosis in São Paulo city, Brazil. BMC Res Notes 4: 269.

45. Miranda SS, Carvalho WaS, Suffys PN, Kritski AL, Oliveira M, et al. (2011) Spoligotyping of clinical Mycobacterium tuberculosis isolates from the state of Minas Gerais, Brazil. Mem Inst Oswaldo Cruz 106: 267-273.

46. Luiz Ro S, Suffys P, Barroso EC, Kerr LR, Duarte CR, et al. (2013) Genotyping and drug resistance patterns of Mycobacterium tuberculosis strains observed in a tuberculosis high-burden municipality in Northeast, Brazil. Braz J Infect Dis 17: 338-345.

47. Perizzolo PF, Dalla Costa ER, Ribeiro AW, Spies FS, Ribeiro MO, et al. (2012) Characteristics of multidrug-resistant Mycobacterium tuberculosis in southern Brazil. Tuberculosis (Edinb) 92: 56-59.

48. 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 U S A 94: 9869-9874.

49. Mokrousov I, Vyazovaya A, Narvskaya O (2014) Mycobacterium tuberculosis Latin American-Mediterranean family and its sublineages in the light of robust evolutionary markers. J Bacteriol 196: 1833-1841.

50. Warren RM, Streicher EM, Sampson SL, van der Spuy GD, Richardson M, et al. (2002) Microevolution of the direct repeat region of Mycobacterium tuberculosis: implications for interpretation of spoligotyping data. J Clin Microbiol 40: 4457-4465.

51. Kato-Meda M, Gagneux S, Flores LL, Kim EY, Small PM, et al. (2011) Strain classification of Mycobacterium tuberculosis: congruence between large sequence polymorphisms and spoligotypes. Int J Tuberc Lung Dis 15: 131-133.

52. Hill V, Zozio T, Sadikalay S, Viegas S, Streit E, et al. (2012) MLVA based classification of Mycobacterium tuberculosis complex lineages for a robust phylogeographic snapshot of its worldwide molecular diversity. PLoS One 7: e41991.

53. Filliol I, Motiwala AS, Cavatore M, Qj W, Hazbón MH, et al. (2006) Global phylogeny of Mycobacterium tuberculosis based on single nucleotide polymorphism (SNP) analysis: insights into tuberculosis evolution, phylogenetic accuracy of other DNA fingerprinting systems, and recommendations for a minimal standard SNP set. J Bacteriol 188: 759-772.

54. Borile C, Labarre M, Franz S, Sola C, Refrégier G (2011) Using affinity propagation for identifying subspecies among clonal organisms: lessons from M. tuberculosis. BMC Bioinformatics 12: 224.

55. Borile C, Labarre M, Franz S, Sola C, Refrégier G (2011) Using affinity propagation for identifying subspecies among clonal organisms: lessons from M. tuberculosis. BMC Bioinformatics 12: 224.

56. Nakanishi N, Wada T, Arikawa K, Millet J, Rastogi N, et al. (2013) Evolutionary robust SNPs reveal the misclassification of Mycobacterium tuberculosis Beijing family strains into sublineages. Infect Genet Evol 16: 174-177.

57. Stavrum R, Mphahlele M, Ovreås K, Muthivhi T, Fourie PB, et al. (2009) High diversity of Mycobacterium tuberculosis genotypes in South Africa and preponderance of mixed infections among ST53 isolates. J Clin Microbiol 47: 1848-1856.

58. Lopes JS, Marques I, Soares P, Nebenzahl-Guimaraes H, Costa J, et al. (2013) SNP typing reveals similarity in Mycobacterium tuberculosis genetic diversity between Portugal and Northeast Brazil. Infect Genet Evol 18: 238-246.

59. Chuang PC, Chen YM, Chen HY, Jou R (2010) Single nucleotide polymorphisms in cell wall biosynthesis-associated genes and phylogeny of Mycobacterium tuberculosis lineages. Infect Genet Evol 10: 459-466.

60. Roetzer A, Diel R, Kohl TA, Rückert C, Nübel U, et al. (2013) Whole genome sequencing versus traditional genotyping for investigation of a Mycobacterium tuberculosis outbreak: a longitudinal molecular epidemiological study. PLoS Med 10: e1001387.

61. Ford C, Yusim K, Ioerger T, Feng S, Chase M, et al. (2012) Mycobacterium tuberculosis–heterogeneity revealed through whole genome sequencing. Tuberculosis (Edinb) 92: 194-201.

62. Lazzarini LC, Rosenfeld J, Huard RC, Hill V, Lapa e Silva JR, et al. (2012) Mycobacterium tuberculosis spoligotypes that may derive from mixed strain infections are revealed by a novel computational approach. Infect Genet Evol 12: 798-806.

63. Von Groll A, Martin A, Felix C, Prata PF, Honscha G, et al. (2010) Fitness study of the RDRio lineage and Latin American-Mediterranean family of Mycobacterium tuberculosis in the city of Rio Grande, Brazil. FEMS Immunol Med Microbiol 58: 119-127.

64. Vinhas SA, Palaci M, Marques HS, Lobo de Aguiar PP, Ribeiro FK, et al. (2013) Mycobacterium tuberculosis DNA fingerprint clusters and its relationship with RD(Rio) genotype in Brazil. Tuberculosis (Edinb) 93: 207-212.
<figure><img src="image_5.png" /><figcaption>**Figure S1:** Dendrogram constructed with BioNumerics software version 6.6 with the results of MIRU-VNTRs 24 loci and spoligotyping by similarity coefficient for categorical data and the neighbor-joining algorithm. 1st column (after spoligotypes): number of isolated label; 2nd column: patients who have more than one isolate in the study (n = 27) received numbering 1–27, and the different strains present the same numbering; 3rd column: International Spolotype Types (SIT); 4th column: classification obtained through SITVITWEB (family and subfamily).</figcaption></figure>
2 2 2 1 2 4 1 4 2 3 2 1 2 1 4 4 3 4 1 5 2 4 7 2
2 2 2 1 2 4 1 4 2 3 2 1 2 1 4 4 3 4 1 5 2 4 4 2
2 2 3 1 2 6 1 4 2 3 2 1 2 1 4 3 3 4 1 4 2 2 8 2
2 2 2 1 2 4 1 4 2 3 2 1 2 1 4 4 3 4 1 5 2 4 7 2
2 2 4 2 2 6 1 5 3 3 2 1 2 1 4 4 3 4 1 4 2 3 1 2
2 2 4 2 2 6 1 5 3 3 2 1 2 1 4 4 3 4 1 2 3 2 7 2
2 2 4 2 2 6 1 5 3 3 2 1 2 4 4 3 3 4 1 5 2 3 7 2
2 4 2 2 6 1 5 3 3 2 1 2 1 4 4 3 4 1 5 2 3 7 2
2 2 4 2 2 6 1 5 3 3 2 5 2 1 4 5 3 4 1 5 3 1 2
2 2 5 2 2 6 1 5 3 3 2 3 2 1 4 4 3 4 1 5 2 2 1 2
1 2 2 5 1 4 3 2 1 2 1 4 4 3 4 1 2 3 1 2
2 2 3 2 2 6 1 4 3 3 2 1 2 1 4 4 3 5 1 5 2 4 5 2
2 2 2 2 2 6 1 4 3 3 2 1 2 1 4 4 3 4 1 6 2 4 1
2 2 2 2 2 6 1 4 3 3 2 1 2 1 4 4 3 4 1 6 2 4 2
2 2 4 1 1 5 1 7 3 3 2 1 2 2 3 4 3 4 1 3 2 6 3
2 2 4 2 2 5 1 6 3 3 2 1 2 2 4 4 3 4 1 3 2 2 6 2
2 2 4 2 2 5 1 3 3 3 1 2 2 4 4 3 4 1 3 2 2 2 2
2 2 2 2 2 5 1 5 3 3 2 1 2 2 4 4 3 4 1 2 2 3 4 2
2 2 7 2 2 5 1 6 3 3 2 1 2 2 4 4 3 4 1 3 2 2 6 2
2 2 5 3 1 3 1 4 3 3 2 1 2 2 2 4 2 2 1 3 2 1 5 2
2 2 3 2 3 1 1 3 3 2 3 3 2 3 2 2 2 4 3 3 2 5 2
2 2 5 3 1 3 1 5 3 3 2 3 3 2 3 1 3 2 4 3 2 2 3 3
2 2 2 3 2 5 1 4 3 3 2 1 3 2 3 2 4 4 4 3 3 3 1 3
2 2 3 1 2 5 1 5 3 3 2 1 14 2 4 1 4 2 4 5 2
2 2 3 3 2 5 1 3 3 2 1 2 4 2 4 3 4 1 2 2 3 4 2
2 2 3 3 2 5 1 8 3 5 3 3 2 0 4 4 3 2 6 1
2 2 5 3 2 3 1 5 3 3 2 3 3 2 3 2 3 2 4 3 3 3 1 3
2 2 5 3 2 3 1 5 3 3 2 3 2 2 3 1 3 2 4 3 3 3 7 3
3 3 3 3 2 5 1 6 3 2 2 2 4 2 4 3 1 4 1 3 3 3 1 2
3 3 3 2 2 5 1 6 3 3 2 2 4 3 4 3 1 2 2 2 3 5
3 3 3 2 2 5 1 6 3 2 2 2 4 2 4 3 1 4 2 3 0 3 0 3
3 3 3 2 2 5 1 6 3 2 2 2 4 2 4 3 1 4 2 3 0 3 7 2
2 2 3 3 2 5 1 2 3 3 2 1 2 2 4 2 2 3 2 3 7 2 6 2
2 2 2 2 3 1 5 3 2 1 2 2 4 4 2 4 1 3 2 2 3 3
2 2 5 2 3 2 2 3 1 2 1 2 2 3 4 2 4 1 1 4 2 2 4
2 2 4 2 5 1 5 3 3 2 1 3 1 4 1 4 4 1 2 2 8 2
2 2 3 1 2 6 1 4 2 3 2 1 2 1 3 5 2 4 1 4 2 2 4 2
2 2 3 2 2 6 1 5 3 3 2 1 1 1 4 3 2 4 1 4 2 3 6 3
2 2 3 3 2 5 1 4 3 3 2 5 2 1 4 2 1 4 2 2 3 3 2
2 2 4 2 2 6 1 3 3 3 2 4 2 2 4 4 2 4 2 1 4 3 7 0
2 2 4 3 2 6 1 3 3 3 2 3 2 4 4 2 4 2 1 1 3 3 2
2 2 4 3 2 6 1 3 3 3 2 3 2 4 4 2 4 2 1 1 3 3 2
2 2 4 3 2 6 1 3 3 3 2 3 2 4 4 2 4 2 1 1 3 3 2
2 2 4 3 2 3 1 5 3 3 2 1 2 2 3 3 3 2 1 3 2 4 5 3
2 2 5 3 1 3 3 5 3 3 2 3 3 2 3 2 2 2 4 3 3 3 5 2
2 2 5 3 1 3 1 5 3 3 2 1 3 2 3 2 2 2 4 3 3 1 5 2
2 2 2 2 2 5 1 3 3 3 2 3 3 2 3 5 2 4 4 3 0 2 3 3
2 2 5 3 2 5 1 3 3 3 2 3 3 2 3 4 2 4 2 3 3 2 4 2
3 2 4 3 2 2 1 6 3 3 3 3 2 3 2 2 2 4 3 7 5 1 2
2 2 3 3 2 4 1 5 3 3 2 2 3 2 3 2 2 2 4 3 3 2 1 3
2 2 5 3 1 3 1 5 2 3 3 2 3 3 2 3 2 2 2 4 3 3 5 2
2 2 5 3 1 3 1 5 3 3 2 1 3 2 3 2 2 2 4 3 3 1 5 2
2 2 3 3 1 3 1 5 3 3 2 2 3 2 3 1 2 2 4 3 3 4 7 3
2 2 5 3 2 3 2 5 3 3 2 3 3 2 3 2 3 2 4 3 3 6 7 3
2 2 5 3 1 3 1 6 3 3 2 3 3 2 3 2 2 2 4 3 0 5 7 3
2 2 5 3 2 3 1 6 3 3 2 3 3 2 3 2 2 2 4 3 3 5 7 2
2 2 1 1 4 1 2 2 5 2 2 4 5 4 4 1 3 2 2 2
2 2 4 0 1 1 3 1 5 3 3 2 1 2 2 3 2 2 4 1 3 2 2 4
2 2 5 3 2 3 1 5 3 2 2 3 2 2 4 4 2 2 1 3 2 2 5 2
2 2 4 2 3 1 5 3 2 2 1 2 2 4 2 2 4 1 3 2 5
1 2 4 3 1 6 5 3 2 2 1 2 2 3 2 3 2 1 3 2 3 3 1
1 2 4 3 1 6 1 3 3 2 2 2 1 3 4 3 3 2 1 3 2 2 3 2
2 2 4 3 1 6 1 5 3 2 2 1 2 2 4 3 4 1 3 2 2 1 2
2 2 4 3 1 6 1 5 3 2 2 6 2 2 4 4 3 4 2 3 2 2 4 2
2 2 4 3 1 6 3 5 4 2 2 7 2 2 4 4 5 4 4 3 3 3 2 2
2 2 3 3 1 6 1 5 3 4 2 3 2 2 2 4 4 3 4 2 1 1 3 4 2
2 2 4 2 2 5 1 5 3 3 6 1 2 2 2 4 1 4 1 3 3 2 7
1 3 3 2 5 1 5 3 3 2 2 2 2 4 4 3 4 1 2 2 3 2
2 2 3 3 2 6 1 5 3 2 1 2 1 4 4 3 4 1 3 2 3
2 2 3 3 2 5 1 5 3 3 2 1 2 1 4 4 2 4 1 5 3 6 3
2 2 4 2 2 6 1 5 3 3 2 1 1 1 4 3 3 4 1 4 2 4 1 3
2 2 4 2 2 6 1 5 3 3 2 1 2 1 4 4 3 4 1 5 2 3 6 2
2 2 5 1 2 5 1 1 3 0 2 2 2 1 3 2 3 4 2 4 2 2 4 3
1 2 2 1 2 3 3 4 2 1 3 4 2 4 2 1 3 3 5 2
2 3 3 2 5 1 2 3 3 2 4 2 4 2 2 2 2
2 2 4 1 2 5 1 5 2 2 2 1 3 3 4 3 2 4 1 3 2 4 7 3
2 2 4 3 2 5 1 6 3 4 2 1 3 3 4 3 2 4 4 3 2 2 8 3
2 2 4 3 2 5 1 5 3 2 2 1 3 3 3 3 2 4 4 3 2 2 5 3
2 2 4 3 2 5 1 6 3 2 2 1 3 3 4 3 4 4 5 2 1 6 3
2 2 4 3 2 5 1 6 3 2 2 1 1 3 4 3 2 4 4 3 2 3 5 3
2 2 4 3 2 5 1 6 3 2 2 1 3 3 4 3 2 4 4 3 2 3 1 3
2 2 3 4 5 2 6 1 3 2 1 3 4 3 2 4 2 3 3 7 4
2 2 4 3 2 7 4 1 6 3 2 9 1 3 3 4 3 3 4 3 2 3 1

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

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

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

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

<figure><img src="image_10.png" /></figure>
<table><tr><td>2</td><td>2</td><td>2</td><td>3</td><td>2</td><td>5</td><td>1</td><td>5</td><td>3</td><td>3</td><td>2</td><td>1</td><td>2</td><td>2</td><td>3</td><td>2</td><td>2</td><td>4</td><td>1</td><td>3</td><td>3</td><td>2</td><td>3</td></tr><tr><td>5</td><td>2</td><td>5</td><td>3</td><td>2</td><td>2</td><td>1</td><td>5</td><td>3</td><td>3</td><td>2</td><td>4</td><td>3</td><td>2</td><td>3</td><td>2</td><td>2</td><td>4</td><td>4</td><td>3</td><td>3</td><td>4</td><td>6</td><td>3</td></tr><tr><td>2</td><td>2</td><td>5</td><td>3</td><td>2</td><td>5</td><td>1</td><td>5</td><td>3</td><td>3</td><td>2</td><td>4</td><td>3</td><td>2</td><td>3</td><td>2</td><td>2</td><td>4</td><td>1</td><td>3</td><td>2</td><td>4</td><td>6</td><td>2</td></tr><tr><td>2</td><td>2</td><td>4</td><td>3</td><td>2</td><td>5</td><td>1</td><td>5</td><td>3</td><td>3</td><td>2</td><td>3</td><td>2</td><td>4</td><td>2</td><td>2</td><td>2</td><td>6</td><td>4</td><td>3</td><td>3</td><td>3</td><td>4</td><td>4</td></tr><tr><td>2</td><td>2</td><td>5</td><td>3</td><td>2</td><td>5</td><td>1</td><td>3</td><td>2</td><td>2</td><td>3</td><td>3</td><td>2</td><td>3</td><td>4</td><td>2</td><td>2</td><td>4</td><td>3</td><td>3</td><td>2</td><td>4</td><td>3</td><td>4</td></tr><tr><td>2</td><td>2</td><td>4</td><td>3</td><td>2</td><td>5</td><td>1</td><td>3</td><td>2</td><td>2</td><td>3</td><td>3</td><td>3</td><td>2</td><td>3</td><td>4</td><td>4</td><td>3</td><td>2</td><td>3</td><td>7</td><td>3</td><td>3</td><td>7</td><td>1</td></tr><tr><td>2</td><td>2</td><td>3</td><td>2</td><td>5</td><td>1</td><td>5</td><td>3</td><td>3</td><td>2</td><td>3</td><td>2</td><td>3</td><td>2</td><td>3</td><td>4</td><td>4</td><td>3</td><td>3</td><td>3</td><td>7</td><td>1</td><td>3</td><td>1</td><td>3</td></tr><tr><td>2</td><td>1</td><td>6</td><td>3</td><td>2</td><td>5</td><td>1</td><td>5</td><td>3</td><td>3</td><td>2</td><td>3</td><td>2</td><td>3</td><td>2</td><td>4</td><td>4</td><td>3</td><td>3</td><td>2</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td></tr><tr><td>2</td><td>2</td><td>4</td><td>3</td><td>2</td><td>5</td><td>1</td><td>4</td><td>3</td><td>2</td><td>3</td><td>3</td><td>2</td><td>3</td><td>2</td><td>4</td><td>4</td><td>3</td><td>3</td><td>8</td><td>3</td><td>3</td><td>3</td><td>5</td><td>7</td><td>3</td></tr><tr><td>2</td><td>2</td><td>4</td><td>3</td><td>2</td><td>5</td><td>1</td><td>4</td><td>3</td><td>3</td><td>2</td><td>3</td><td>3</td><td>2</td><td>3</td><td>4</td><td>4</td><td>4</td><td>3</td><td>3</td><td>5</td><td>7</td><td>3</td><td>3</td><td>5</td><td>7</td><td>2</td></tr><tr><td>2</td><td>2</td><td>4</td><td>3</td><td>2</td><td>5</td><td>1</td><td>4</td><td>3</td><td>3</td><td>2</td><td>3</td><td>3</td><td>2</td><td>3</td><td>2</td><td>4</td><td>2</td><td>3</td><td>3</td><td>6</td><td>8</td><td>3</td><td>6</td><td>8</td><td>3</td><td>6</td><td>2</td></tr><tr><td>2</td><td>2</td><td>5</td><td>3</td><td>2</td><td>3</td><td>1</td><td>5</td><td>3</td><td>3</td><td>2</td><td>3</td><td>3</td><td>2</td><td>3</td><td>2</td><td>4</td><td>4</td><td>2</td><td>3</td><td>3</td><td>6</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td></tr><tr><td>2</td><td>2</td><td>5</td><td>3</td><td>2</td><td>3</td><td>1</td><td>5</td><td>3</td><td>3</td><td>2</td><td>3</td><td>3</td><td>2</td><td>3</td><td>2</td><td>2</td><td>2</td><td>1</td><td>3</td><td>1</td><td>2</td><td>4</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td></tr><tr><td>2</td><td>2</td><td>5</td><td>4</td><td>1</td><td>3</td><td>1</td><td>5</td><td>3</td><td>3</td><td>2</td><td>3</td><td>3</td><td>2</td><td>3</td><td>1</td><td>2</td><td>2</td><td>4</td><td>3</td><td>2</td><td>7</td><td>3</td><td>3</td><td>7</td><td>3</td><td>7</td><td>2</td></tr><tr><td>2</td><td>2</td><td>5</td><td>3</td><td>2</td><td>3</td><td>3</td><td>5</td><td>2</td><td>1</td><td>1</td><td>3</td><td>3</td><td>3</td><td>2</td><td>3</td><td>2</td><td>4</td><td>3</td><td>3</td><td>6</td><td>7</td><td>3</td><td>3</td><td>6</td><td>7</td><td>3</td><td>7</td><td>2</td></tr><tr><td>2</td><td>2</td><td>3</td><td>3</td><td>1</td><td>3</td><td>1</td><td>5</td><td>3</td><td>3</td><td>1</td><td>3</td><td>2</td><td>3</td><td>2</td><td>3</td><td>2</td><td>4</td><td>3</td><td>2</td><td>6</td><td>7</td><td>3</td><td>3</td><td>6</td><td>7</td><td>3</td><td>7</td><td>2</td></tr><tr><td>2</td><td>2</td><td>5</td><td>3</td><td>1</td><td>3</td><td>1</td><td>5</td><td>3</td><td>3</td><td>1</td><td>1</td><td>3</td><td>2</td><td>3</td><td>2</td><td>3</td><td>2</td><td>4</td><td>3</td><td>3</td><td>6</td><td>7</td><td>3</td><td>3</td><td>6</td><td>7</td><td>3</td><td>7</td><td>2</td></tr><tr><td>2</td><td>2</td><td>5</td><td>3</td><td>1</td><td>3</td><td>1</td><td>5</td><td>3</td><td>3</td><td>1</td><td>1</td><td>3</td><td>2</td><td>3</td><td>2</td><td>3</td><td>2</td><td>4</td><td>3</td><td>3</td><td>3</td><td>7</td><td>3</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td></tr><tr><td>2</td><td>2</td><td>5</td><td>3</td><td>1</td><td>3</td><td>1</td><td>5</td><td>3</td><td>2</td><td>2</td><td>3</td><td>2</td><td>2</td><td>4</td><td>4</td><td>2</td><td>4</td><td>2</td><td>3</td><td>3</td><td>5</td><td>7</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td></tr><tr><td>2</td><td>2</td><td>5</td><td>2</td><td>2</td><td>3</td><td>1</td><td>5</td><td>3</td><td>3</td><td>2</td><td>3</td><td>3</td><td>2</td><td>3</td><td>2</td><td>2</td><td>2</td><td>4</td><td>3</td><td>4</td><td>7</td><td>2</td><td>4</td><td>7</td><td>2</td><td>2</td><td>2</td><td>2</td></tr><tr><td>3</td><td>2</td><td>5</td><td>3</td><td>1</td><td>2</td><td>1</td><td>2</td><td>3</td><td>3</td><td>2</td><td>2</td><td>3</td><td>2</td><td>3</td><td>2</td><td>2</td><td>2</td><td>4</td><td>3</td><td>3</td><td>4</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td></tr><tr><td>2</td><td>2</td><td>5</td><td>3</td><td>1</td><td>3</td><td>1</td><td>2</td><td>3</td><td>3</td><td>2</td><td>2</td><td>3</td><td>2</td><td>3</td><td>2</td><td>2</td><td>2</td><td>4</td><td>3</td><td>3</td><td>2</td><td>7</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td></tr><tr><td>2</td><td>2</td><td>5</td><td>3</td><td>1</td><td>3</td><td>1</td><td>5</td><td>3</td><td>3</td><td>2</td><td>3</td><td>3</td><td>2</td><td>3</td><td>2</td><td>2</td><td>2</td><td>4</td><td>3</td><td>4</td><td>7</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td></tr><tr><td>2</td><td>2</td><td>5</td><td>3</td><td>1</td><td>3</td><td>1</td><td>5</td><td>3</td><td>3</td><td>2</td><td>3</td><td>3</td><td>2</td><td>3</td><td>2</td><td>2</td><td>2</td><td>4</td><td>3</td><td>3</td><td>5</td><td>7</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td></tr><tr><td>2</td><td>2</td><td>5</td><td>3</td><td>1</td><td>3</td><td>1</td><td>5</td><td>3</td><td>2</td><td>2</td><td>3</td><td>2</td><td>3</td><td>2</td><td>2</td><td>2</td><td>4</td><td>3</td><td>3</td><td>5</td><td>7</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td></tr><tr><td>2</td><td>2</td><td>5</td><td>3</td><td>1</td><td>3</td><td>1</td><td>5</td><td>3</td><td>2</td><td>2</td><td>3</td><td>2</td><td>3</td><td>2</td><td>2</td><td>2</td><td>4</td><td>3</td><td>3</td><td>5</td><td>7</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td></tr><tr><td>2</td><td>2</td><td>2</td><td>3</td><td>1</td><td>3</td><td>1</td><td>5</td><td>3</td><td>3</td><td>2</td><td>3</td><td>3</td><td>2</td><td>3</td><td>2</td><td>2</td><td>2</td><td>4</td><td>3</td><td>3</td><td>5</td><td>7</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td></tr><tr><td>2</td><td>2</td><td>2</td><td>3</td><td>2</td><td>6</td><td>1</td><td>6</td><td>3</td><td>3</td><td>2</td><td>3</td><td>3</td><td>2</td><td>4</td><td>2</td><td>2</td><td>4</td><td>2</td><td>3</td><td>2</td><td>5</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td></tr><tr><td>2</td><td>2</td><td>2</td><td>1</td><td>2</td><td>5</td><td>1</td><td>5</td><td>3</td><td>2</td><td>2</td><td>3</td><td>3</td><td>2</td><td>4</td><td>2</td><td>1</td><td>4</td><td>3</td><td>4</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td></tr><tr><td>2</td><td>2</td><td>3</td><td>3</td><td>2</td><td>5</td><td>1</td><td>4</td><td>3</td><td>3</td><td>2</td><td>1</td><td>3</td><td>2</td><td>4</td><td>2</td><td>5</td><td>4</td><td>2</td><td>3</td><td>2</td><td>2</td><td>2</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td></tr><tr><td>2</td><td>2</td><td>3</td><td>3</td><td>2</td><td>5</td><td>1</td><td>1</td><td>3</td><td>3</td><td>2</td><td>1</td><td>3</td><td>2</td><td>4</td><td>2</td><td>3</td><td>4</td><td>2</td><td>3</td><td>0</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td></tr><tr><td>2</td><td>2</td><td>3</td><td>3</td><td>2</td><td>6</td><td>1</td><td>5</td><td>3</td><td>3</td><td>2</td><td>3</td><td>2</td><td>2</td><td>3</td><td>4</td><td>4</td><td>2</td><td>1</td><td>...</td><td>3</td><td>5</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td></tr><tr><td>2</td><td>2</td><td>3</td><td>3</td><td>2</td><td>5</td><td>1</td><td>5</td><td>3</td><td>3</td><td>2</td><td>2</td><td>3</td><td>2</td><td>4</td><td>3</td><td>8</td><td>4</td><td>2</td><td>3</td><td>0</td><td>4</td><td>5</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td></tr><tr><td>2</td><td>2</td><td>3</td><td>3</td><td>2</td><td>5</td><td>1</td><td>5</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td><td>2</td><td>4</td><td>2</td><td>4</td><td>2</td><td>3</td><td>4</td><td>4</td><td>2</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td></tr><tr><td>2</td><td>2</td><td>4</td><td>2</td><td>2</td><td>6</td><td>1</td><td>1</td><td>3</td><td>3</td><td>2</td><td>3</td><td>3</td><td>2</td><td>4</td><td>2</td><td>2</td><td>4</td><td>2</td><td>1</td><td>1</td><td>4</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td></tr><tr><td>2</td><td>2</td><td>3</td><td>3</td><td>2</td><td>5</td><td>1</td><td>5</td><td>3</td><td>3</td><td>2</td><td>4</td><td>3</td><td>2</td><td>5</td><td>2</td><td>2</td><td>4</td><td>2</td><td>1</td><td>2</td><td>4</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td></tr><tr><td>2</td><td>2</td><td>3</td><td>5</td><td>2</td><td>5</td><td>1</td><td>5</td><td>3</td><td>3</td><td>2</td><td>4</td><td>2</td><td>2</td><td>4</td><td>2</td><td>2</td><td>4</td><td>2</td><td>3</td><td>3</td><td>5</td><td>5</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td></tr><tr><td>2</td><td>2</td><td>4</td><td>4</td><td>2</td><td>5</td><td>1</td><td>6</td><td>3</td><td>3</td><td>2</td><td>4</td><td>3</td><td>2</td><td>4</td><td>2</td><td>2</td><td>4</td><td>2</td><td>3</td><td>4</td><td>4</td><td>5</td><td>4</td><td>4</td><td>5</td><td>4</td><td>4</td><td>5</td></tr><tr><td>2</td><td>2</td><td>3</td><td>1</td><td>2</td><td>5</td><td>1</td><td>2</td><td>3</td><td>3</td><td>2</td><td>4</td><td>3</td><td>2</td><td>5</td><td>2</td><td>2</td><td>4</td><td>2</td><td>3</td><td>3</td><td>3</td><td>5</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td></tr><tr><td>2</td><td>2</td><td>3</td><td>1</td><td>2</td><td>5</td><td>1</td><td>2</td><td>3</td><td>3</td><td>2</td><td>4</td><td>3</td><td>2</td><td>5</td><td>2</td><td>2</td><td>4</td><td>2</td><td>3</td><td>3</td><td>3</td><td>5</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td></tr><tr><td>2</td><td>2</td><td>3</td><td>2</td><td>6</td><td>1</td><td>5</td><td>3</td><td>3</td><td>2</td><td>4</td><td>2</td><td>2</td><td>2</td><td>2</td><td>3</td><td>4</td><td>2</td><td>3</td><td>3</td><td>6</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td></tr><tr><td>2</td><td>2</td><td>3</td><td>3</td><td>2</td><td>5</td><td>1</td><td>4</td><td>3</td><td>3</td><td>2</td><td>3</td><td>2</td><td>2</td><td>4</td><td>3</td><td>2</td><td>4</td><td>2</td><td>3</td><td>2</td><td>5</td><td>7</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td></tr><tr><td>5</td><td>3</td><td>4</td><td>2</td><td>2</td><td>1</td><td>1</td><td>5</td><td>3</td><td>3</td><td>4</td><td>4</td><td>2</td><td>4</td><td>4</td><td>1</td><td>4</td><td>2</td><td>3</td><td>4</td><td>5</td><td>3</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td></tr><tr><td>2</td><td>3</td><td>4</td><td>4</td><td>2</td><td>5</td><td>1</td><td>5</td><td>3</td><td>4</td><td>2</td><td>5</td><td>4</td><td>2</td><td>4</td><td>4</td><td>1</td><td>4</td><td>2</td><td>3</td><td>5</td><td>5</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td></tr><tr><td>2</td><td>1</td><td>4</td><td>2</td><td>5</td><td>1</td><td>4</td><td>3</td><td>3</td><td>2</td><td>3</td><td>4</td><td>2</td><td>3</td><td>4</td><td>1</td><td>4</td><td>2</td><td>3</td><td>3</td><td>2</td><td>4</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td></tr><tr><td>2</td><td>3</td><td>3</td><td>3</td><td>2</td><td>5</td><td>1</td><td>4</td><td>3</td><td>3</td><td>2</td><td>3</td><td>4</td><td>2</td><td>4</td><td>4</td><td>1</td><td>4</td><td>2</td><td>3</td><td>1</td><td>3</td><td>8</td><td>3</td><td>8</td><td>3</td><td>8</td><td>3</td><td>8</td></tr><tr><td>2</td><td>2</td><td>3</td><td>3</td><td>2</td><td>5</td><td>1</td><td>2</td><td>3</td><td>3</td><td>2</td><td>4</td><td>2</td><td>2</td><td>4</td><td>2</td><td>2</td><td>2</td><td>3</td><td>2</td><td>3</td><td>7</td><td>2</td><td>6</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td></tr><tr><td>2</td><td>2</td><td>4</td><td>3</td><td>2</td><td>6</td><td>1</td><td>5</td><td>3</td><td>2</td><td>6</td><td>1</td><td>2</td><td>4</td><td>0</td><td>2</td><td>4</td><td>1</td><td>3</td><td>2</td><td>2</td><td>7</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td></tr><tr><td>2</td><td>5</td><td>4</td><td>3</td><td>2</td><td>6</td><td>2</td><td>2</td><td>3</td><td>5</td><td>1</td><td>3</td><td>2</td><td>4</td><td>2</td><td>3</td><td>8</td><td>3</td><td>2</td><td>3</td><td>2</td><td>3</td><td>3</td><td>1</td><td>3</td><td>2</td><td>3</td><td>2</td><td>3</td></tr><tr><td>1</td><td>4</td><td>2</td><td>7</td><td>1</td><td>5</td><td>3</td><td>2</td><td>2</td><td>2</td><td>2</td><td>4</td><td>3</td><td>4</td><td>1</td><td>3</td><td>2</td><td>2</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td></tr><tr><td>1</td><td>2</td><td>4</td><td>3</td><td>2</td><td>6</td><td>1</td><td>5</td><td>3</td><td>3</td><td>2</td><td>1</td><td>2</td><td>2</td><td>4</td><td>1</td><td>4</td><td>4</td><td>4</td><td>3</td><td>0</td><td>2</td><td>1</td><td>4</td><td>2</td><td>4</td><td>2</td><td>4</td><td>2</td></tr><tr><td>2</td><td>2</td><td>4</td><td>3</td><td>2</td><td>5</td><td>1</td><td>5</td><td>3</td><td>3</td><td>2</td><td>3</td><td>3</td><td>2</td><td>4</td><td>1</td><td>4</td><td>4</td><td>1</td><td>3</td><td>3</td><td>2</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td></tr><tr><td>2</td><td>2</td><td>4</td><td>3</td><td>2</td><td>5</td><td>1</td><td>5</td><td>3</td><td>3</td><td>2</td><td>3</td><td>3</td><td>2</td><td>4</td><td>2</td><td>2</td><td>2</td><td>4</td><td>3</td><td>4</td><td>3</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td></tr><tr><td>2</td><td>2</td><td>2</td><td>3</td><td>2</td><td>5</td><td>1</td><td>5</td><td>3</td><td>2</td><td>2</td><td>3</td><td>3</td><td>2</td><td>4</td><td>2</td><td>4</td><td>4</td><td>3</td><td>2</td><td>3</td><td>4</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td></tr></table>

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

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

<table><tr><td>C7765</td><td>1926</td><td>T1</td></tr><tr><td>c0635</td><td>521</td><td>T1</td></tr><tr><td>C0665</td><td>53</td><td>T1</td></tr><tr><td>C6697</td><td>25</td><td>Orphan T3</td></tr><tr><td>C1579</td><td>1163</td><td>T3</td></tr><tr><td>C7321</td><td>8</td><td>T1</td></tr><tr><td>C2104</td><td>53</td><td>T1</td></tr><tr><td>C1926</td><td>746</td><td>H3</td></tr><tr><td>C0395</td><td>53</td><td>T1</td></tr><tr><td>C0837</td><td>65</td><td>T1</td></tr><tr><td>C1171</td><td>27</td><td>T3</td></tr><tr><td>C0982</td><td>157</td><td>T3</td></tr><tr><td>C1842</td><td>50</td><td>H3</td></tr><tr><td>C1120</td><td>53</td><td>T1</td></tr><tr><td>c0821</td><td>50</td><td>H3</td></tr><tr><td>C0825</td><td>50</td><td>H3</td></tr><tr><td>C0859</td><td>50</td><td>H3</td></tr><tr><td>C0888</td><td>50</td><td>H3</td></tr><tr><td>c0317</td><td>50</td><td>H3</td></tr><tr><td>C2213A2</td><td>741</td><td>H3</td></tr><tr><td>C2142A211</td><td>50</td><td>H3</td></tr><tr><td>C2142</td><td>11</td><td>H3</td></tr><tr><td>C0766</td><td>777</td><td>H4/Ural-1</td></tr><tr><td>C1433</td><td>27</td><td>H3</td></tr><tr><td>C1427</td><td>5</td><td>T1</td></tr><tr><td>C1999</td><td>50</td><td>H3</td></tr><tr><td>C2454</td><td>50</td><td>H3</td></tr><tr><td>C2144</td><td>1227</td><td>T5-Madrid2</td></tr><tr><td>C2159</td><td>1227</td><td>T5-Madrid2</td></tr><tr><td>C4509</td><td>53</td><td>T1</td></tr><tr><td>C7448</td><td>53</td><td>T1</td></tr><tr><td>C8190</td><td>53</td><td>T1</td></tr><tr><td>C0874</td><td>53</td><td>T1</td></tr><tr><td>C0459</td><td>53</td><td>T1</td></tr><tr><td>c0627A2</td><td>17</td><td>H3</td></tr><tr><td>C7829</td><td>53</td><td>T1</td></tr><tr><td>C8123</td><td>53</td><td>T1</td></tr><tr><td>C1130</td><td>6</td><td>T1</td></tr><tr><td>C0746</td><td>53</td><td>T1</td></tr><tr><td>C1696</td><td>53</td><td>T1</td></tr><tr><td>C7157</td><td>291</td><td>T1</td></tr><tr><td>C0807</td><td>52</td><td>T2</td></tr><tr><td>C0696</td><td>52</td><td>T2</td></tr><tr><td>C1210</td><td>153</td><td>T2</td></tr><tr><td>C2463</td><td>12</td><td>T1</td></tr><tr><td>C2484</td><td>12</td><td>T1</td></tr><tr><td>C1839</td><td>189</td><td>T1</td></tr><tr><td>C1474</td><td>102</td><td>T</td></tr><tr><td>C3501</td><td>3907</td><td>T1</td></tr><tr><td>C1378</td><td>10</td><td>T1</td></tr><tr><td>C7674</td><td>Orphan EAI5</td><td></td></tr><tr><td>C0901</td><td>Orphan T4</td><td></td></tr><tr><td>C0517</td><td>Orphan T1</td><td></td></tr><tr><td>C0627</td><td>17</td><td>T1</td></tr><tr><td>C7195</td><td>25</td><td>X2</td></tr><tr><td>C1754</td><td>Orphan T1</td><td></td></tr></table>
<figure><img src="image_13.png" /><figcaption>**Figure S2:** Allelic diversity of the 24 loci MIRU-VNTR loci in LAM/RDRio isolates.</figcaption></figure>
**Figure S3:** Possible evolution of *M. tuberculosis* lineage Euro-American (LAM, T and H) according to the markers analyzed in this study. The Euro-American Ancestral evolved into two distinct groups: LAM ancestral and T/H ancestral, characterized by absences of spacers 33–36 and one copy of MIRU24 that is common to all Euro-American lineages. The ancestral LAM strains have fbpC103 and this was the basis of LAM A (LAM9), with additional absence of spacers 21–24, one copy of MTUB30 and two copies of MIRU04 and ETR-A. LAM A on its turn is the basis for two other groups: LAM B (LAM9-LAM4, LAM1-LAM2-LAM5 and LAM6) and LAM C (LAM3). The LAM B evolved from LAM B1, characterized by a deleted RD174 and on its turn to LAM B2, with both deleted RD174 and RDRio. The H/T Ancestral lineage is the origin of both groups H and T (difference only in MIRU-VNTR copies), showing absence of spacers 33–36; additional loss of spacer 31 led to subtype H A, observed in high frequency in this study.

<figure><img src="image_14.png" /></figure>
<table><caption>Supplemental table S2: Detailed genotyping results of *fbpC*<sup>103</sup>, RD174 and RD<sup>Rio</sup> on 182 *M. tuberculosis* strains isolated from pulmonary tuberculosis patients in Rio de Janeiro, Brazil.</caption><thead><tr><th>Strain</th><th colspan="12">Spoligotype Description</th><th>MIRU-VNTR 24 loci</th><th>Clade SITVIT2</th><th>SIT SITVIT2</th><th>designation</th><th>fbpC<sup>103</sup></th><th>DR174</th><th>RD<sup>Rio</sup></th></tr></thead><tbody><tr><td>C2170A2</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C2013</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C1559</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C1599</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C2009</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C2349</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C2559</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C3554</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C7516</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C1772A2</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C2024</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C2046</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C2017</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C1602</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C0879</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C1318</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C1358</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C1106</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C1020</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C1256</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C4991</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C1231A2</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C3198</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C6949</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C2519</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C0470A2</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>c0817</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C0963</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C1598</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C1473</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C6921</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>c7823</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C2314</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C0438A3</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C6582</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C2629</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C2537</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C1377A2</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C2028</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C4641</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C0669A2</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C0143</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C2768</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C7891</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C1559A2</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C1377</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C1938</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C7670A2</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C0438</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C0537A2</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C0996</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C2105A2</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C2105</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C0818</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr><tr><td>C3165</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td><td>■</td></tr></tbody></table>