Skip to main content
  • Research article
  • Open access
  • Published:

Genetic diversity and molecular epidemiology of multidrug-resistant Mycobacterium tuberculosis in Minas Gerais State, Brazil

Abstract

Background

We aimed to characterize the genetic diversity of drug-resistant Mycobacterium tuberculosis (MTb) clinical isolates and investigate the molecular epidemiology of multidrug-resistant (MDR) tuberculosis from Minas Gerais State, Brazil.

Methods

One hundred and four MTb clinical isolates were assessed by IS6110-RFLP, 24-locus mycobacterial interspersed repetitive units variable-number tandem repeats (MIRU-VNTR), TB-SPRINT (simultaneous spoligotyping and rifampicin-isoniazid drug-resistance mutation analysis) and 3R-SNP-typing (analysis of single-nucleotide polymorphisms in the genes involved in replication, recombination and repair functions).

Results

Fifty-seven different IS6110-RFLP patterns were found, among which 50 had unique patterns and 17 were grouped into seven clusters. The discriminatory index (Hunter and Gaston, HGDI) for RFLP was 0.9937. Ninety-nine different MIRU-VNTR patterns were found, 95 of which had unique patterns and nine isolates were grouped into four clusters. The major allelic diversity index in the MIRU-VNTR loci ranged from 0.6568 to 0.7789. The global HGDI for MIRU-VNTR was 0.9991. Thirty-two different spoligotyping profiles were found: 16 unique patterns (n = 16) and 16 clustered profiles (n = 88). The HGDI for spoligotyping was 0.9009. The spoligotyped clinical isolates were phylogenetically classified into Latin-American Mediterranean (66.34 %), T (14.42 %), Haarlem (5.76 %), X (1.92 %), S (1.92 %) and U (unknown profile; 8.65 %). Among the U isolates, 77.8 % were classified further by 3R-SNP-typing as 44.5 % Haarlem and 33.3 % LAM, while the 22.2 % remaining were not classified. Among the 104 clinical isolates, 86 were identified by TB-SPRINT as MDR, 12 were resistant to rifampicin only, one was resistant to isoniazid only, three were susceptible to both drugs, and two were not successfully amplified by PCR. A total of 42, 28 and eight isolates had mutations in rpoB positions 531, 526 and 516, respectively. Correlating the cluster analysis with the patient data did not suggest recent transmission of MDR-TB.

Conclusions

Although our results do not suggest strong transmission of MDR-TB in Minas Gerais (using a classical 100 % MDR-TB identical isolates cluster definition), use of a smoother cluster definition (>85 % similarity) does not allow us to fully eliminate this possibility; hence, around 20–30 % of the isolates we analyzed might be MDR-TB transmission cases.

Peer Review reports

Background

Multidrug-resistant (MDR) tuberculosis (TB) is an increasingly serious global public health threat that requires robust, efficient and quick actions to improve the control and spread of drug-resistant clinical isolates. TB isolates with resistance to isoniazid (INH) and rifampin (RIF), defined as MDR-TB, are prone to sequential accumulation of mutations in the target genes that confer resistance to them [1].

Brazil ranks sixteenth among the world’s 22 countries with high TB burdens; here, the TB prevalence is 92,000 cases, the incidence rate is 35.4 per 100,000 per year, and the mortality rate is 4.9 per 100,000 of the population according to World Health Organization estimates [2]. Minas Gerais State has the fourth lowest TB incidence (17.9/100 000/year) in Brazil [3]. Although the current prevalence of primary MDR-TB is relatively low in Brazil, it has potential to become a major public health issue, because resistance to more than one drug has been shown to be strongly associated with household contact, as suggested in a recent study conducted in Amazonia [4].

Mycobacterium tuberculosis complex (MTBC) genotyping methods have been widely used for investigating epidemics involving MDR-TB [5]. These methods help to define the recent transmission factors for MDR-TB isolates and enable better control programs to be initiated to avoid MDR-TB expansion at local or global population levels.

Among the various genetic markers available for studying genetic polymorphism in drug-resistant M. tuberculosis (MTb), restriction fragment length polymorphism (RFLP) analysis of the IS6110 insertion sequence is a “gold standard” for MTBC typing [57]; this method has been used widely to identify and investigate TB transmission and re-infection rates, and was used at the end of the 1990s to identify cross-contamination in laboratories [7].

However, the IS6110-RFLP technique has many disadvantages in routine practice: it is laborious, requires trained staff, replies on microgram quantities of purified DNA, and has poor discriminatory power when applied to isolates with low IS6110 copy numbers [5, 6, 8]. The method based on mycobacterial interspersed repetitive units variable-number tandem repeats (MIRU-VNTR) has progressively replaced IS6110-RFLP. MIRU-VNTR, a very powerful technique, provides adequate discrimination between MTb clinical isolates and is comparable to IS6110-RFLP in terms of its accuracy for estimating TB outbreaks and for use in phylogenetic investigations [5, 813]. MIRU-VNTR analysis also generates readily comparable numerical values, a useful feature for interlaboratory studies; hence, over the last 10 years it has come into standard use in TB research [9], but may progressively be replaced by whole genome sequencing (WGS) [14].

Spacer oligonucleotide typing (spoligotyping), another well-established (1996) technique, is based on the polymorphisms found in the clustered regularly interspersed short palindromic repeats (CRISPR) of MTb. Spoligotyping detects the presence or absence of 43 spacer sequences in the CRISPR region of MTb [5, 15, 16]. Spoligotyping data can be expressed in an octal or binary format. Additionally, spoligotyping requires small amounts of crude or purified DNA and its results are as portable as MIRU-VNTR, making data from it easily shared between laboratories. Spoligotyping has provided a lot of highly informative results on the phylogeographic distribution of genotypic diversity in MTb. Furthermore, spoligotyping in combination with MIRU-VNTR has excellent discriminatory power for cluster analysis of tubercle bacilli genomics, making it a valuable tool for the epidemiology and evolutionary biology of MTb [1719]. However, because the discriminatory power of spoligotyping is generally inferior to that of IS6110-based RFLP, it cannot be used alone for molecular epidemiology studies [2022].

The use of single-nucleotide polymorphisms (SNPs) as markers of genetic variation for phylogenetic analysis has been described in many studies [2328]. Because SNPs offer important advantages for high-throughput analyses, they are the markers of choice in genetics research, likewise are regions where deletions in the genome occur [29]. Despite MTBC structural genes exhibiting very low levels of polymorphism among strains [3033], higher polymorphism levels were found recently in several genes among which were those involved in replication, recombination and repair functions (3R genes) [34]. Indeed, the high-throughput 3R-SNP-typing method is able to classify undefined spoligotype signatures making it an efficient, easy to use tool for evolutionary studies on MTBC clinical isolates [23].

For genetic characterization of TB drug resistance, molecular detection tests currently search for known mutations in different TB-specific target genes [33]. For MDR-TB, the most frequent mutations associated with RIF and INH resistance can be assessed by sequencing, line-probe assays or other methods [28, 29]. RIF resistance is mainly (95 %) caused by the 81-bp rifampin resistance-determining region (RRDR) of the rpoB gene [35, 36]. INH resistance is often caused by mutations in katG (codon 315), inhA (positions −15 and −8 in the inhA promoter sequence), and in other genes [35, 37]. Phenotypic TB culture-based drug susceptibility testing (DST), however, remains the gold standard for diagnosis of MDR-TB [38]. Tuberculosis-spoligo-rifampin-isoniazid typing (TB-SPRINT), a 59-plex multiplexed microbead-based, high-throughput DNA array method, provides simultaneous spoligotyping and mutation analysis of the most common resistance-associated SNPs for RIF (rpoB RRDR direct and indirect coverage) and INH resistance (katG, inhA) [36, 3840]. This method will soon be replaced by an improved 77-plex version (TB-SPRINT-plus) capable of identifying mutations in target genes conferring resistance to some second-line drugs (Molina et al., unpublished observations).

Here, we aimed to investigate the genetic profile of MDR-TB clinical isolates from Minas Gerais, a Brazilian state, using the following four molecular techniques: IS6110-RFLP, MIRU-VNTR, TB-SPRINT and 3R-SNPs-typing.

Methods

Clinical isolates and drug susceptibility testing

One hundred and four MDR-TB clinical isolates were collected from 2008 to 2013 in Minas Gerais, Brazil, each one corresponding to a unique TB patient. All the isolates were obtained by culture of respiratory samples and represent distinct local laboratories in Minas Gerais. The isolates were referred to the Ezequiel Dias Foundation (FUNED) for culture, identification and drug susceptibility testing. Samples were transferred to the Research Laboratory of Mycobacteria of the Faculty of Medicine of the Federal University of Minas Gerais (UFMG) where they now belong to the MTb clinical isolate collection. The samples are representative of MDR-TB in Minas Gerais. All the clinical isolates were from patients diagnosed with MDR-TB by the BACTEC™ MGIT™ 960 System [41]. Demographic data were obtained from the Information System on Diseases of Compulsory Declaration of Brazil (otherwise known as SINAN).

Genomic DNA extraction

MTb genomic DNA was extracted from mycobacterial colonies subcultured on Löwenstein-Jensen (LJ) medium. One loopful of mycobacterial colonies was collected in a tube containing 500 μL of TE buffer (10 mM Tris-Cl, 1 mM EDTA) and then incubated at 80 °C for 60 min. Lysozyme 10 mg/mL (70 μL) was added to each tube, followed by incubation at 65 °C for 15 min with occasional mixing, after which 5 M NaCl (100 μL) and 10 % cetyltrimethylammonium bromide (CTAB) (100 μL) were added to each sample. After adding 70 μL of 10 % SDS and 6 μL of proteinase K (10 mg/mL) the samples were vortexed briefly and then incubated at 65 °C for 15 min. Chloroform/isoamyl alcohol (24:1 v/v) (700 μL) was added to each tube, and the solution was centrifuged for 20 min at 12,000 rpm at 4 °C in a microcentrifuge. The supernatant was transferred to a new 1.5 mL microcentrifuge tube, 450 μL of ice-cold isopropanol was added, and the tube was inverted 20 times to precipitate the nucleic acids. Samples were incubated overnight at −20 °C and then centrifuged at 12,000 rpm for 30 min at 4 °C in a microcentrifuge, after which the supernatant was discarded. The pellet was air-dried for 2 h and then resuspended in 60 μL of TE buffer (10 mM Tris-Cl, 1 mM EDTA).

Genotyping methods

IS6110-RFLP

DNA samples were typed by IS6110-RFLP analysis in accordance with the standardized protocol described by van Embden et al. [7] and van Soolingen et al. [42]. The reference strain used was Mt 14323.

MIRU-VNTR

The standard 24 MIRU-VNTR loci method [18] was performed based on agarose gel electrophoresis. The simplex PCR product size was determined as previously reported [43].

TB-SPRINT

High-throughput TB-SPRINT was performed at the Institute of Genetics and Microbiology at the University Paris-Sud, France, on a Luminex 200™ flow cytometry device (Luminex Corp, Austin, TX) as previously described, using a microbead-based DNA array method [4446]. The TB-SPRINT analysis was performed according to the standardized protocol recommended by Gomgnimbou et al. [37].

3R-SNP typing

The 3R-SNP typing was performed as described by Abadia et al. [23]. This seven gene multiplex-PCR method uses primers designed on the dual-priming oligonucleotide principle, which has been shown to strongly increase the mutated to wild-type signal ratio [47].

Bioinformatic cluster analysis

All results (except IS6110-RFLP) were entered into Excel® spreadsheets, and then transferred to BioNumerics™ software version 6.6 (Applied Maths, Sint-Martens-Latem, Belgium). IS6110 RFLP fingerprints were digitalized and compared using the Dice coefficient and the unweighted-pair group method using average linkage (UPGMA) according to the manufacturer’s instructions [44]. MIRU-VNTR data were analyzed using the categorical coefficient and UPGMA [45]. TB-SPRINT and 3R-SNP-typing data were analyzed in BioNumerics™ using the Jaccard index and UPGMA [37]. Spoligotyping data were also analyzed using the minimum spanning tree (MST) method, as shown in Fig. 1. A composite data set for the four methods mentioned above and a composite dendrogram were also built (Fig. 2) [46]. Cluster definition was based on identical patterns using the above four methods (tighter definition) or by setting the percentage similarity at >85 % (smoother definition) [48, 49]. The recent transmission index was determined by computing the n and (n minus 1) index [50, 51].

Fig. 1
figure 1

Minimum Spanning Tree (MST) obtained from the spoligotyping dataset (n = 104 MDR-TB isolates) identified from 2008 to 2013 in Minas Gerais, Brazil. Color code indicates major subclades found within the Lineage 4 (Euro-American), LAM, T, H, X, S and U

Fig. 2
figure 2

Results Matrix of the 104 studied clinical isolates. a: final composite dendrogram built using UPGMA using five results data set (b: IS6110-RFLP; c: MIRU-VNTR; d: spoligotyping; e: RIF-INH SNPs typing; f: 3R-SNPs-typing; g: subclade color code of Fig. 1; h: (Key) Unique Clinical isolate number; i: city of isolation; j: (SIT) Spoligo-international-type label

Spoligotyping patterns were assigned a SIT (spoligotyping-international-type) using the SITVITWEB website (http://www.pasteur-guadeloupe.fr:8081/SITVIT_ONLINE/) [19]. MIRU-VNTR data were analyzed using the MIRU-VNTRplus website (http://www.miru-vntrplus.org/).

The VNTR allelic diversity index (h) was used to evaluate the allelic diversity of the various VNTR loci. The value of h was calculated using the formula described by Selander et al. [52]. The discriminatory power of each typing method was also computed and compared by the Hunter-Gaston discriminatory index (HGDI) [53, 54]. The HGDI was calculated using the discriminatory power calculator available at http://insilico.ehu.es/mini_tools/discriminatory_power/index.php.

To correlate the molecular and patient data, an analysis of basic demographic patient data (city and address) and familiar data (mother’s name) was performed.

Ethics statement

The study was approved by the Ethics Committee of the Federal University of Minas Gerais (number 122.941; CAAE 06611912.8.0000.5149). Isolates from this study were obtained by culturing stock clinical isolates.

Results and discussion

This is the first molecular characterization of MDR-TB clinical isolates from Minas Gerais State. The study is based on 104 clinical isolates obtained from patients in Minas Gerais State with MDR-TB-positive cultures between 2008 and 2013. The isolates, from 46 cities in the State of Minas Gerais, are considered to be representative of the MDR-TB strains in this region of Brazil. The mean age of the patients was 43.24 years and the male-to-female ratio was 1.66:1. Drug susceptibility testing (BACTEC™ MGIT™ 960 System method) showed all of 104 isolates were resistant to rifampin and isoniazid at least.

IS6110-RFLP typing

The IS6110 copy number from each isolate was assessed from the number of bands hybridizing with the probe. The 104 clinical isolates were typed and a total of 67 fingerprint patterns were obtained (65 %). The majority of them (94.03 %) had multiple IS6110 copies (7–15). This high degree of IS6110 polymorphism is in accordance with the results observed in drug-susceptible clinical isolates and suggests a recent low rate of MDR-TB transmission [5560]. That we observed a low frequency of isolates with low IS6110 copy numbers (5.97 %) demonstrates the excellent discriminatory power of IS6110-RFLP in our setting. Thirty-seven isolates (35 %) lacked IS6110-RFLP profiles, probably resulting from poor DNA quantity. Fifty-seven different IS6110 RFLP patterns were identified, 50 of which were unique, and 17 isolates were found in seven clusters (HGDI = 0.9937).

MIRU-VNTR

The 104 isolates were successfully typed and 99 different MIRU-VNTR patterns were found. Among these patterns, 95 were unique, and nine isolates belonged to four clusters (HGDI = 0.9991). Lineage signature was performed by MIRU-VNTRplus best-match labeling using 24 global MIRU-VNTR loci and the following six main lineages/sublineages were observed: Cameroon (3.85 %), CAS/Delhi (0.96 %), Haarlem (17.3 %), LAM (73.08 %), S (1.92 %) and T2-Uganda (2.88 %). The high HGDI of the 24-MIRU-VNTR confirms the importance of using this technique, either in association with RFLP or with spoligotyping [17, 18, 51, 6062].

The allelic diversity of each MIRU-VNTR locus in our setting was evaluated and classified into highly (HGDI >0.6), moderately (0.6> HGDI <0.3) or poorly discriminative (HGDI <0.3) [52], as summarized in Table 1. The highest allelic diversity indexes were for Qub 26, Mtub 04, MIRU 26, MIRU 16, Qub 11, and MIRU 10. The allelic diversity index was low (h ≤ 0.3) for five of the 24 loci. As supported by this study and others, partial MIRU-VNTR genotyping could be sufficient to define epi-linked clusters after first-line and high-throughput spoligotyping [63]. Ali et al. described seven loci with the highest discriminatory level that could be used preferentially to investigate possible transmission events [61].

Table 1 Allelic diversity of each MIRU-VNTR locus and their discriminatory power

TB-SPRINT typing

All the clinical isolates were successfully typed by spoligotyping and were classified phylogenetically into five lineages and 11 sublineages, as shown in Table 2. Thirty-two different spoligotyping patterns were found; 16 of them had unique patterns, and 88 isolates were grouped into 16 clusters (HGDI = 0.9009). A minimum spanning tree (MST) was built (Fig. 1). The single lineage is lineage four (Euro-American) with a majority of LAM, T, H and a minority of S and X2. The two major sublineages (LAM: n = 69 and T: n = 15) are in central positions of the MST (Fig. 1). They show tight links between patterns and represent the major proportion of the clinical isolates in the Minas Gerais State (n = 84/104 or 80.76 % of the clinical isolates). The continuous transmission of MTB in certain settings is strongly affected by the prevailing population structure of tubercle bacilli [63], which induces the predominance of a homogeneous group, such as the LAM family in South America and similarly for the Beijing family in Asia [19, 57, 6466]. LAM, T and H clinical isolates are more likely to become MDR in Brazil [67, 68].

Table 2 Phylogenetical classification of multidrug-resistant Mycobacterium tuberculosis (n = 104) by TB-SPRINT method

A further characterization of the RIF-INH typing scheme (16-plex) was performed for 102 of 104 (98.7 %) isolates, either on the 81 base-pair RRDR for RIF or in inhA and katG for INH. For two isolates, none of the positions tested were amplifiable. For the isolates typed at the RRDR locus, 98 (96.07 %) had mutated genotypes, among which 78 (76.4 % of total) carried defined targeted mutations conferring RIF resistance. However, 20 isolates with undefined RIF mutations were detected and four samples had no mutations. A total of 42, 28 and eight isolates had mutations at rpoB positions 531 (38 samples rpoB 531TTG and four samples rpoB 531TGG), 526 (19 samples rpoB 526GAC and nine samples rpoB 526TAC) and 516 (eight samples rpoB 516GTC), respectively, and only one of these samples carried mutations in the three rpoB positions. Out of 102 isolates typed at loci known to be involved in INH resistance, 87 (85.29 %) had a mutated genotype at either katG315 or within the inhA promoter. Among them, 82 (94.2 %) carried the following defined mutations conferring INH resistance: 54 samples harbored katG 315ACC (65.85 %), four samples katG 315AAC (4.87 %), nine samples inhA-15 (10.97 %), 14 samples harbored a double katG 315ACC + inhA-15 mutation and one sample harbored a double katG 315AAC + inhA-15 mutation. These results are in accordance with other first-line TB drug resistance studies in Brazil [67, 69, 70]. Altogether, out of 102 isolates typed for both loci, 86 were resistant to both drugs (INH and RIF). Forty-five different RIF-INH profiles were obtained for the isolates and 25 of them had a unique profile, while 17 profiles clustered into 72 isolates (HGDI = 0.9561).

3R-SNP-typing

For some spoligotyping patterns (n = 9), it was not possible to assign their lineages/sublineages. To assign such isolates (called “U” isolates in SpolDB4 and/or SITVITWEB), a 3R-SNP scheme was used successfully to solve 77.8 % of the cases [23]. The patterns found were Haarlem (n = 4) and LAM (n = 3), and these signatures were confirmed by MIRU-VNTR typing. Only two samples remained unclassified.

Out of 104 isolates, the 3R-SNP-typing method allowed us to find mutations in 94 (90.38 %) isolates associated with four specific genotype families: LAM (n = 67), Haarlem (n = 12), X (n = 2), T2 (n = 1) and an unknown lineage (n = 12). The ten isolates remaining (9.62 %) did not amplify successfully. Use of the 3R-SNP-based method helped to clarify the infra-specific taxonomy of our sampling, thus improving our confidence in the evolutionary analysis of our data [23, 71].

Comparison of the discriminatory powers of the genotyping methods

The discriminatory powers of the IS6110-RFLP, MIRU-VNTR, TB-SPRINT (spoligotyping and RIF-INH-typing) and 3R-SNP methodologies are shown in Table 3. Use of a combination of different techniques is important for improved epidemiological and phylogeographical interpretation of molecular results [17, 18, 51, 61, 62]. MIRU-VNTR has the highest discriminatory power followed by IS6110-RFLP and TB-SPRINT. 3R-SNP-typing has a lower discriminatory power because only seven SNPs were used in the current format.

Table 3 Discriminatory index of IS6110-RFLP, MIRU-VNTR, TB-SPRINT (Spoligotyping and RIF-INH typing), 3R-SNPs

Molecular epidemiology in the Minas Gerais State

Among the 104 clinical isolates, 71 displayed a low similarity index (<85 %) and 33 a high similarity index (>85 %). Twelve clusters without any obvious epidemiological link were observed. According to the (n minus 1) or n Recent Transmission Index [50, 51] definition (where there is the choice to diminish or not to diminish all the clusters by one index case) and using the smooth cluster definition (>85 %), the maximum transmission rate of MDR-TB in Minas Gerais State would be (33 minus 12/104), or 20 % using the (n minus 1) method, and 33/104 or 31 % using the n method. However, if we assume a 100 % identity cluster definition, no 100 % identity cluster was found, which suggests that no cases of MDR-TB transmission occurred in Minas Gerais State. Our results clearly point to an extended classical “stone in the pound” epidemiological analysis of the 12 suspected clusters, for which MDR-TB transmission remains likely [72].

Our results are possibly explained by the fact that MDR-TB transmission in Minas Gerais is individually acquired (i.e., there are no primary MDR-TB cases). An alternative explanation is that many cases of MDR-TB were missed, but this seems unlikely because the sampling is representative of MDR-TB cases in Minas Gerais. The situation for Minas Gerais differs from that of other Brazilian studies [67, 73] in that the other studies did not use all the techniques used herein, which may have increased the discriminatory power of our analysis. This discordance could also be explained by the global differences in TB prevalence between different regions of Brazil.

When looking more closely at the geographical origin of the samples from Minas Gerais, we could identify only one factor that made us suspicious of an epidemiological link (in eight patients from the same city). Our results suggest that clustered genotypes indicative of recent MDR-TB transmission should be interpreted with caution, unless direct evidence of epidemiological links between clustered cases can be demonstrated [74].

Despite the low number of samples, this collection of MTBC isolates is likely to be representative of the confirmed MDR-TB cases in Minas Gerais State. One limitation of the present study, however, is the absence of clinical epidemiological data. Also, contact tracing for individual patients could not be performed.

Future long-term studies are necessary to identify the possible risk factors for the emergence of drug resistance and/or treatment failure. Additionally, longitudinal studies in regions of Brazil with a high incidence of MDR-TB are now urgently needed.

Conclusions

To sum up, use of four different discriminant genotyping techniques (IS6110-RFLP, MIRU-VNTR, TB-SPRINT and 3R-SNP-typing) provided useful data for phylogenetic evaluation and fine taxonomic characterization of MDR-TB clinical isolates from Minas Gerais State, Brazil. The most common MDR-TB isolates belonged to the LAM lineage and approximately two thirds of them did not provide evidence for recent transmission of MDR-TB. Our data indicate that MDR-TB in Minas Gerais State is caused by clinical isolates that were not transmitted in recent years or that the outbreak is driven by individually acquired resistance and endogenous reactivation. This situation contrasts with the findings from other Brazilian studies, which all reported a high transmission rate for MDR-TB. Such an important issue requires locally-adapted solutions and state-specific control measures in Brazil. Continuous surveillance of MDR-TB transmission could be improved by introduction of new diagnostic tools and epidemiological research using WGS methods.

Abbreviations

MTb :

Mycobacterium tuberculosis

MDR:

Multidrug-resistant

TB:

Tuberculosis

HGDI:

Hunter-Gaston discriminatory index

LAM:

Latin-American Mediterranean

U:

Unknown profile

INH:

Isoniazid

RIF:

Rifampin

MTBC:

Mycobacterium tuberculosis complex

RFLP:

Restriction fragment length polymorphism

IS:

Insertion sequence

MIRU-VNTR:

Mycobacterial interspersed repetitive units variable-number tandem repeats

WGS:

Whole genome sequencing

Spoligotyping:

Spacer oligonucleotide typing

CRISPR:

Clustered regularly interspersed short palindromic repeats

PCR:

Polymerase chain reaction

RRDR:

Rifampin resistance-determining region

DST:

Drug susceptibility testing

SNPs:

Single-nucleotide polymorphism

TB-SPRINT:

Tuberculosis-spoligo-rifampin-isoniazid typing

3R:

Replication, recombination and repair functions

FUNED:

Ezequiel Dias Foundation

UFMG:

Federal University of Minas Gerais

SINAN:

Information system on diseases of compulsory declaration

LJ:

Löwenstein-Jensen

DPO:

Dual-priming oligonucleotide

UPGMA:

Unweighted-pair group method using average linkage

MST:

Minimum spanning trees

SIT:

Spoligotyping-international-type

h:

Allelic diversity index

RTI:

Recent transmission index

References

  1. Brossier F, Sola C, Millota G, Jarliera V, Veziris N, Sougakoff W. Comparison of a Semiautomated Commercial Repetitive-Sequence-Based PCR Method with Spoligotyping, 24-Locus Mycobacterial Interspersed Repetitive-Unit–Variable-Number Tandem-Repeat Typing, and Restriction Fragment Length Polymorphism-Based Analysis of IS6110 for Mycobacterium tuberculosis Typing. J Clin Microbiol. 2014;52:4082–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. World Health Organization. WHO report 2013. Global tuberculosis control: surveillance, planning, financing. WHO/HTM/ TB/2013.306. Geneva, Switzerland: WHO; 2013.

    Google Scholar 

  3. Brazil. Tuberculose: casos confirmados notificados no Sistema de Informação de Agravos de Notificação (Sinan). Brasilia, Brazil: Ministério da Saúde, 2004. http://dtr2004.saude.gov.br/sinanweb/tabnet/dh?sinan/tuberculose/bases/tubercbr.def Accessed January 2015. [Portuguese]

  4. da Silva GM, Ramasawmy R, Perez-Porcuna TM, Zaranza E, Chrusciak Talhari A, Martinez-Espinosa FE, et al. Primary drug resistance among pulmonary treatment-naive tuberculosis patients in Amazonas State. Brazil Int J Tuberc Lung Dis. 2014;18:559–63.

    Article  Google Scholar 

  5. Sougakoff W. Molecular epidemiology of multidrug-resistant clinical isolates of Mycobacterium tuberculosis. Clin Microbiol Infect. 2011;17:800–5. doi:10.1111/j.1469-0691.2011.03577.

    Article  CAS  PubMed  Google Scholar 

  6. McEvoy CR, Falmer AA, van Pittius NC G, Victor TC, van Helden PD, Warren RM. The role of IS6110 in the evolution of Mycobacterium tuberculosis. Tuberculosis. 2007;87:393–404. doi:10.1016/j.tube.2007.05.010.

    Article  CAS  PubMed  Google Scholar 

  7. van Embden JD, Cave MD, Crawford JT, Dale JW, Eisenach KD, Gicquel B, et al. Strain identification of Mycobacterium tuberculosis by DNA fingerprinting: recommendations for a standardized methodology. J Clin Microbiol. 1993;31:406–9.

    PubMed  PubMed Central  Google Scholar 

  8. Sun YJ, Bellamy R, Lee ASG, Ng ST, Ravindran S, Wong SY, et al. Use of Mycobacterial Interspersed Repetitive Unit-Variable-Number Tandem Repeat Typing to examine genetic diversity of Mycobacterium tuberculosis in Singapore. J Clin Microbiol. 2004;42:1986–93. doi:10.1128/JCM.42.5.1986-1993.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Supply P, Mazars E, Lesjean S, Vincent V, Gicquel B, Locht C. Variable human minisatellite‐like regions in the Mycobacterium tuberculosis genome. Mol Microbiol. 2000;36:762–71.

    Article  CAS  PubMed  Google Scholar 

  10. Cowan LS, Mosher L, Diem L, Massey JP, Crawford JT. Variable-number tandem repeat typing of Mycobacterium tuberculosis isolates with low copy numbers of IS6110 by using mycobacterial interspersed repetitive units. J Clin Microbiol. 2002;40:1592–602.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Mazars E, Lesjean S, Banuls AL, Gilbert M, Vincent V, Gicquel B, et al. High-resolution minisatellite-based typing as a portable approach to global analysis of Mycobacterium tuberculosis molecular epidemiology. Proc Natl Acad Sci U S A. 2001;98:1901–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Sola C, Filliol I, Legrand E, Lesjean S, Locht C, Supply P, et al. Genotyping of the Mycobacterium tuberculosis complex using MIRUs: association with VNTR and spoligotyping for molecular epidemiology and evolutionary genetics. Infect Genet Evol. 2003;3:125–33.

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Roetzer A, Diel R, Kohl TA, Rückert C, Nübel U, Blom J, et al. Whole Genome Sequencing versus Traditional Genotyping for Investigation of a Mycobacterium tuberculosis Outbreak: A Longitudinal Molecular Epidemiological Study. PLoSMed. 2013;10(2), e1001387.

    Google Scholar 

  15. Zhang J, Abadia E, Refregier G, Tafaj S, Boschiroli ML, Guillard B, et al. 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. 2010;59:285–94.

    Article  PubMed  Google Scholar 

  16. Kato‐Maeda M, Metcalfe JZ, Flores L. Genotyping of Mycobacterium tuberculosis: application in epidemiologic studies. Future Microbiol. 2011;6:203–16.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Sola C. Clustered regularly interspersed short palindromic repeats (CRISPR) genetic diversity studies as a mean to reconstruct the evolution of the Mycobacterium tuberculosis complex. Tuberculosis (Edinb). 2015;95:S159–66.

    Article  CAS  Google Scholar 

  18. Supply P, Allix C, Lesjean S, Cardoso-Oelemann M, Rusch-Gerdes S, Willery E, et al. Proposal for Standardization of Optimized Mycobacterial Interspersed Repetitive Unit-Variable-Number Tandem Repeat Typing of Mycobacterium tuberculosis. J Clin Microbiol. 2006;44:4498–510.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Demay C, Liens B, Burguière T, Hill V, Couvin D, Millet J, et al. SITVITWEB – A publicly available international multimarker database for studying Mycobacterium tuberculosis genetic diversity and molecular epidemiology. Infect Genet Evol. 2012;12(4):755–66.

    Article  CAS  PubMed  Google Scholar 

  20. Kurepina N, Likhoshvay E, Shashkina E, Mathema B, Kremer K, van Soolingen D, et al. Targeted hybridization of IS6110 fingerprints identifies the W-Beijing Mycobacterium tuberculosis clinical isolates among clinical isolates. J Clin Microbiol. 2005;43:2148–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Driscoll JR. Spoligotyping for Molecular Epidemiology of the Mycobacterium tuberculosis Complex Molecular Epidemiology of Microorganisms. Methods in Molecular Biology™. 2009;551:117–28.

    Article  CAS  Google Scholar 

  22. Kamerbeek J, Schouls L, Kolk A, van Agterveld M, van Soolingen D, Kuijper S, et al. Simultaneous detection and strain differentiation of Mycobacterium tuberculosis for diagnosis and epidemiology. J Clin Microbiol. 1997;35:907–14.

    CAS  PubMed  PubMed Central  Google Scholar 

  23. Abadia E, Zhang J, dos Vultos T, Ritacco V, Kremer K, Aktas E, et al. Resolving lineage assignation on Mycobacterium tuberculosis clinical isolates classified by spoligotyping with a new high-throughput 3R SNPs based method. Infect Genet Evol. 2010;10:1066–74.

    Article  CAS  PubMed  Google Scholar 

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

    Article  PubMed  PubMed Central  Google Scholar 

  25. Coll F, McNerney R, Guerra-Assuncao JA, Glynn JR, Perdigao J, Viveiros M, et al. A robust SNP barcode for typing Mycobacterium tuberculosis complex clinical isolates. Nat Commun. 2014;5:4812.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Homolka S, Projahn M, Feuerriegel S, Ubben T, Diel R, Nubel U, et al. High resolution discrimination of clinical Mycobacterium tuberculosis complex clinical isolates based on single nucleotide polymorphisms. PLoS One. 2012;7:e39855.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Feuerriegel S, Koser CU, Niemann S. Phylogenetic polymorphisms in antibiotic resistance genes of the Mycobacterium tuberculosis complex. J Antimicrob Chemother. 2014;69:1205–10.

    Article  CAS  PubMed  Google Scholar 

  28. Stucki D, Malla B, Hostettler S, Huna T, Feldmann J, Yeboah-Manu D, et al. Two new rapid SNP-typing methods for classifying Mycobacterium tuberculosis complex into the main phylogenetic lineages. PLoS One. 2012;7:e41253.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Schork NJ, Fallin D, Lanchbury JS. Single nucleotide polymorphisms and the future of genetic epidemiology. Clin Genet. 2000;58:250–64.

    Article  CAS  PubMed  Google Scholar 

  30. Musser JM, Amin A, Ramaswamy S. Negligible genetic diversity of Mycobacterium tuberculosis host immune system protein targets: evidence of limited selective pressure. Genetics. 2000;155:7–16.

    CAS  PubMed  PubMed Central  Google Scholar 

  31. Sreevatsan S, Pan X, Stockbauer KE, Connell ND, Kreiswirth BN, Whittam TS, et al. Restricted structural gene polymorphism in the Mycobacterium tuberculosis complex indicates evolutionarily recent global dissemination. Proc Natl Acad Sci U S A. 1997;94:9869–74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Kapur V, Whittam TS, Musser JM. Is Mycobacterium tuberculosis 15,000 years old? J Infect Dis. 1994;170:1348–9.

    Article  CAS  PubMed  Google Scholar 

  33. Achtman M. Evolution, population structure, and phylogeography of genetically monomorphic bacterial pathogens. Annu Rev Microbiol. 2008;62:53–70.

    Article  CAS  PubMed  Google Scholar 

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

    Article  PubMed  PubMed Central  Google Scholar 

  35. Gomgnimbou MK, Abadia E, Zhang J, Refregier G, Panaiotov S, Bachiyska E, et al. “Spoligoriftyping”, a dual-priming-oligonucleotide-based direct-hybridization assay for tuberculosis control with a multianalyte microbead-based hybridization system. J Clin Microbiol. 2012;50:3172–9.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Heep M, Brandstätter B, Rieger U, Lehn N, Richter E, Rüsch-Gerdes S, et al. Frequency of rpoB mutations inside and outside the cluster I region in rifampin-resistant clinical Mycobacterium tuberculosis isolates. J Clin Microbiol. 2001;39:107–10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Gomgnimbou M, Hernandez-Neuta I, Panaiotov S, Bachyiska E, Palomino JC, Martin A, et al. TB-SPRINT: TuBerculosis-SPoligo-Rifampin-IsoNiazid Typing”, a one-in all-assay technique for surveillance and control of multi-drug resistant tuberculosis on Luminex® devices J. Clin Microbiol. 2013;51:3527–34.

    Article  CAS  Google Scholar 

  38. Ocheretina O, Escuyer VE, Mabou MM, Royal-Mardi G, Collins S, Vilbrun SC, et al. Correlation between Genotypic and Phenotypic Testing for Resistance to Rifampin in Mycobacterium tuberculosis Clinical Isolates in Haiti: Investigation of Cases with Discrepant Susceptibility Results. PLoS One. 2014;9:e90569.

    Article  PubMed  PubMed Central  Google Scholar 

  39. WANG X, Jiaoa J, Xua W, Chaib X, Lia Z, Wanga Q. A simple, rapid and economic method for detecting multidrug-resistant tuberculosis. Braz J Infect Dis. 2013;17:667–71. doi:10.1016/j.bjid.2013.04.008.

    Article  PubMed  Google Scholar 

  40. Nikolayevskyy V, Balabanova Y, Simak T, Malomanova N, Fedorin I, Drobniewski F. Performance of the Genotype MTBDRPlus assay in the diagnosis of tuberculosis and drug resistance in Samara. Russian Federation BMC Clin Pathol. 2009;9:2.

    Article  PubMed  Google Scholar 

  41. Tortoli E, Benedetti M, Fontanelli A, Simonetti MT. Evaluation of automated BACTEC MGIT 960 system for testing susceptibility of Mycobacterium tuberculosis to four major antituberculous drugs: comparison with the radiometric BACTEC 460 TB method and the agar plate method of proportion. J Clin Microbiol. 2002;40:607–10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. van Soolingen D, de Haas PE, Hermans PW, van Embden JD. DNA fingerprinting of Mycobacterium tuberculosis. Methods Enzymol. 1994;235:196–205.

    Article  CAS  PubMed  Google Scholar 

  43. Bonura C, Gomgnimbou MK, Refrégier G, Aleo A, Fasciana T, Giammanco A, et al. Molecular epidemiology of tuberculosis in Sicily, Italy: what has changed after a decade? BMC Infect Dis. 2014;14:602.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Hwang HY, Chang CY, Chang LL, Chang SF, Chang YH, Chen YJ. Characterization of rifampicin-resistant Mycobacterium tuberculosis in Taiwan. J Med Microbiol. 2003;52:239–45.

    Article  CAS  PubMed  Google Scholar 

  45. Roetzer A, Schuback S, Diel R, Gasau F, Ubben T, di Nauta A, et al. Evaluation of Mycobacterium tuberculosis typing methods in a 4-year study in Schleswig-Holstein. Northern Germany J Clin Microbiol. 2011;49:4173–8.

    Article  PubMed  Google Scholar 

  46. Varma-Basil M, Kumar S, Arora J, Angrup A, Zozio T, Banavaliker JN, et al. Comparison of spoligotyping, mycobacterial interspersed repetitive units typing and IS6110-RFLP in a study of genotypic diversity of Mycobacterium tuberculosis in Delhi. North India Mem Inst Oswaldo Cruz. 2011;106:524–35.

    Article  PubMed  Google Scholar 

  47. Chun JY, Kim KJ, Hwang IT, Kim YJ, Lee DH, Lee IK, et al. Dual priming oligonucleotide system for the multiplex detection of respiratory viruses and SNP genotyping of CYP2C19 gene. Nucleic Acids Res. 2007;35(6):e40.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Borgdorff MW, van de Hof S, Kalisvaart N, Kremer K, van Soolingen D. Influence of Sampling on clustering and Associations With Risk Factors in the Molecular Epidemiology of Tuberculosis. Am J Epidemiol. 2011;174:243–51.

    Article  PubMed  Google Scholar 

  49. Fok A, Numata Y, Schulzer M, Fitzgerald MJ. Risk factors for clustering of tuberculosis cases: a systematic review of population-based molecular epidemiology studies. Int J Tuberc Lung Dis. 2008;2:480–92.

    Google Scholar 

  50. Alland D, Kalkut GE, Moss AR, et al. Transmission of tuberculosis in New York City. An analysis by DNA fingerprinting and conventional epidemiologic methods. N Engl J Med. 1994;330:1710–16.

    Article  CAS  PubMed  Google Scholar 

  51. Murray M, Alland D. Methodological Problems in the Molecular Epidemiology of Tuberculosis. Am J Epidemiol. 2002;155:565–71.

    Article  PubMed  Google Scholar 

  52. Selander RK, Caugant DA, Ochman H, Musser JM, Gilmour MN, Whittam TS. Methods of multilocus enzyme electrophoresis for bacterial population genetics and systematics. Appl Environ Microbiol. 1986;51:873–84.

    CAS  PubMed  PubMed Central  Google Scholar 

  53. Hunter PR, Gaston MA. Numerical index of the discriminatory ability of typing systems: an application of Simpson’s index of diversity. J Clin Microbiol. 1988;26:2465–6.

    CAS  PubMed  PubMed Central  Google Scholar 

  54. Nakagawa T, Shibayama T, Uchiya K, Nikai T, Ogawa K. Comparison of a Variable-Number Tandem-Repeat (VNTR) Method for Typing Mycobacterium avium with Mycobacterial Interspersed Repetitive-Unit–VNTR and IS1245 Restriction Fragment Length Polymorphism Typing. J Clin Microbiol. 2009;47:2156–64.

    Article  PubMed  PubMed Central  Google Scholar 

  55. Zheng C, Zhao E, Zhu G, Li S, Sun H, Feng Q, et al. Suitability of IS6110-RFLP and MIRU-VNTR for Differentiating Spoligotyped Drug-Resistant Mycobacterium tuberculosis Clinical Isolates from Sichuan in China. Biomed Res Int. 2014;2014:763204.

    PubMed  PubMed Central  Google Scholar 

  56. Chauhan A, Chauhan DS, Parashar D, Gupta P, Sharma VD, Sachan AS, et al. DNA Fingerprinting of Mycobacterium tuberculosis isolates from AGRA region by IS6110 probe. Indian J Med Microbiol. 2004;22:238–40.

    CAS  PubMed  Google Scholar 

  57. Niemann S, Sch-Gerdes S, Virarichter A. IS6110 Fingerprinting of drug-resistant Mycobacterium tuberculosis clinical isolates isolated in Germany during 1995. J Clin Microbiol. 1997;35:3015–20.

    CAS  PubMed  PubMed Central  Google Scholar 

  58. Hermans PWM, Messadi F, Guebrexabher H, van Soolingen D, de Haas PEW, Heersma H, et al. Analysis of the population structure of Mycobacterium tuberculosis in Ethiopia, Tunisia, and the Netherlands: usefulness of DNA typing for global tuberculosis epidemiology. J Infect Dis. 1995;171:1504–13.

    Article  CAS  PubMed  Google Scholar 

  59. Strassle A, Putnik J, Weber R, Fehr-Merhof A, Wust J, Pfyffer GE. Molecular epidemiology of Mycobacterium tuberculosis clinical isolates isolated from patients in a human immunodeficiency virus cohort in Switzerland. J Clin Microbiol. 1997;35:374–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  60. Van Soolingen D, Hermans PWM, de Haas PEW, Soll DR, van Embden JDA. Occurrence and stability of insertion sequences in Mycobacterium tuberculosis complex clinical isolates: evaluation of an insertion sequence-dependent DNA polymorphism as a tool in the epidemiology of tuberculosis. J Clin Microbiol. 1991;29:2578–86.

    PubMed  PubMed Central  Google Scholar 

  61. Ali A, Hasan Z, Tanveer M, Siddiqui AR, Ghebremichael S, Kallenius G, et al. Characterization of Mycobacterium tuberculosis Central Asian Strain 1 using mycobacterial interspersed repetitive unit genotyping. BMC Microbiol. 2007;7:76.

    Article  PubMed  PubMed Central  Google Scholar 

  62. Jonsson J, Hoffner S, Berggren I, Bruchfeld J, Ghebremichael S, Pennhag A, et al. Comparison between RFLP and MIRU-VNTR Genotyping of Mycobacterium tuberculosis Clinical isolates Isolated in Stockholm 2009 to 2011. Plos one. 2014;9:e95159.

    Article  PubMed  PubMed Central  Google Scholar 

  63. Asante-Poku A, Nyaho MS, Borrell S, Comas I, Gagneux S, Yeboah-Manu D. Evaluation of Customised Lineage-Specific Set of MIRU-VNTR Loci for Genotyping Mycobacterium tuberculosis Complex Isolates in Ghana. PlosOne. 2014;9:e92675.

    Article  Google Scholar 

  64. Mokrousov I. The quiet and controversial: ural family of Mycobacterium tuberculosis. Infect Genet Evol. 2012;12:619–29.

    Article  PubMed  Google Scholar 

  65. Millet J, Baboolal S, Streit E, Akpaka PE, Rastogi N. A First Assessment of Mycobacterium tuberculosis Genetic Diversity and Drug-Resistance Patterns in Twelve Caribbean Territories. BioMed Res Int. 2014;2014:718496.

    PubMed  PubMed Central  Google Scholar 

  66. Miranda SS, Carvalho WS, Suffys PN, Kritski AL, Oliveira M, Zarade N, et al. Spoligotyping of clinical Mycobacterium tuberculosis isolates from the State of Minas Gerais. Brazil Mem Inst Oswaldo Cruz. 2011;106:267–73.

    Article  PubMed  Google Scholar 

  67. Perizzolo PF, Dalla Costa ER, Ribeiro AW, Spies FS, Ribeiro MO, Dias CF, et al. Characteristics of multidrug-resistant Mycobacterium tuberculosis in southern Brazil. Tuberculosis. 2012;92:56–9.

    Article  PubMed  Google Scholar 

  68. Mendes NH, Melo FAF, Santos ACB, Pandolfi JRC, Almeida EA, Cardoso RF, et al. Characterization of the genetic diversity of Mycobacterium tuberculosis in São Paulo city, Brazil. BMC Research Notes. 2011;4:269.

    Article  PubMed  PubMed Central  Google Scholar 

  69. de Oliveira MM, da Silva AR, Cardoso MO, Gomes HM, Fonseca L, Werneck-Barreto AM, et al. Rapid detection of resistance against rifampicin in isolates of Mycobacterium tuberculosis from Brazilian patients using a reverse-phase hybridization assay. J Microbiol Methods. 2003;53:335–42.

    Article  PubMed  Google Scholar 

  70. Jardim PCR, Zamarioli LA, Coelho AGV, Figueiredo TR, Rozman MA. Resistência do Mycobacterium tuberculosis às drogas no município de São Vicente. Rev Inst Adolfo Lutz. 2001;60:119–23.

    Google Scholar 

  71. Filliol I, Motiwala AS, Cavatore M, Qi W, Hernando MH, Del Valle MB, et al. 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. 2006;188:759–72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  72. Veen J. Microepidemics of tuberculosis: the stone-in-the-pond principle. Tubercle Lung Dis. 1992;73:73–6.

    Article  CAS  Google Scholar 

  73. Ferrazoli L, Palaci M, Marques LR, Jamal LF, Afiune JB, Chimara E, et al. Transmission of tuberculosis in an endemic urban setting in Brazil. Int J Tuberc Lung Dis. 2000;4:18–25.

    CAS  PubMed  Google Scholar 

  74. Narvskaya N, Mokrousov I, Otten T, Vyshnevskiy B. Molecular Markers: Application for Studies of Mycobacterium tuberculosis Population in Russia. In: Read MM, editor. Trends in DNA Fingerprinting Research. New York: Nova; 2005. p. 111–25.

    Google Scholar 

Download references

Acknowledgments

We are grateful to the Research Group in Mycobacteriosis, in Belo Horizonte, Brazil who contributed with intellectual discussion and data acquisition; to Molecular Biology Applied to Mycobacteria team, in Rio de Janeiro for the expert support; to Infection Genetics Emerging Pathogens Evolution team in Orsay, who has made substantial contributions to acquisition, analysis and interpretation of data and provided expert support; to Coordination for the Improvement of Higher Level Personnel (CAPES) who provided funding for research fellowship; to National Research Council (CNPq) and Foundation for Research Support of Minas Gerais (FAPEMIG) (4467962014/APQ 03266–13 and APQ00094-12) for the research financial support; and to Post-Graduate Program in Infectious Diseases and Tropical Medicine of UFMG for research support.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Silvana Spíndola de Miranda.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

Conceived and designed the experiments: NGTD CJA SSM WSC CS. Performed the experiments: NGTD INA LJA CJA HMG. Analyzed the data: NGTD SSM MKG GR CS. Wrote the paper: NGTD SSM HMG PNS CS. All authors read and approved the final manuscript.

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.

The Creative Commons Public Domain Dedication waiver (https://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Dantas, N.G.T., Suffys, P.N., Carvalho, W.d.S. et al. Genetic diversity and molecular epidemiology of multidrug-resistant Mycobacterium tuberculosis in Minas Gerais State, Brazil. BMC Infect Dis 15, 306 (2015). https://doi.org/10.1186/s12879-015-1057-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12879-015-1057-y

Keywords