Tuberculosis is a major health problem in São Paulo, Brazil, which is the most populous and one of the most cosmopolitan cities in South America. To characterize the genetic diversity of Mycobacterium tuberculosis in the population of this city, the genotyping techniques of spoligotyping and MIRU were applied to 93 isolates collected in two consecutive years from 93 different tuberculosis patients residing in São Paulo city and attending the Clemente Ferreira Institute (the reference clinic for the treatment of tuberculosis).
Spoligotyping generated 53 different spoligotype patterns. Fifty-one isolates (54.8%) were grouped into 13 spoligotyping clusters. Seventy- two strains (77.4%) showed spoligotypes described in the international databases (SpolDB4, SITVIT), and 21 (22.6%) showed unidentified patterns. The most frequent spoligotype families were Latin American Mediterranean (LAM) (26 isolates), followed by the T family (24 isolates) and Haarlem (H) (11 isolates), which together accounted for 65.4% of all the isolates. These three families represent the major genotypes found in Africa, Central America, South America and Europe. Six Spoligo-International-types (designated SITs by the database) comprised 51.8% (37/72) of all the identified spoligotypes (SIT53, SIT50, SIT42, SIT60, SIT17 and SIT1). Other SITs found in this study indicated the great genetic diversity of M. tuberculosis, reflecting the remarkable ethnic diversity of São Paulo city inhabitants. The MIRU technique was more discriminatory and did not identify any genetic clusters with 100% similarity among the 93 isolates. The allelic analysis showed that MIRU loci 26, 40, 23 and 10 were the most discriminatory. When MIRU and spoligotyping techniques were combined, all isolates grouped in the 13 spoligotyping clusters were separated.
Our data indicated the genomic stability of over 50% of spoligotypes identified in São Paulo and the great genetic diversity of M. tuberculosis isolates in the remaining SITs, reflecting the large ethnic mix of the São Paulo city inhabitants. The results also indicated that in this city, M. tuberculosis isolates acquired drug resistance independently of genotype and that resistance was more dependent on the selective pressure of treatment failure and the environmental circumstances of patients.
Tuberculosis (TB) is an infectious disease characterized by high morbidity and mortality in developing countries and in urban areas of developed countries . Tuberculosis is a major health problem in São Paulo city with an incidence of 65.2 cases/100,000 inhabitants . The city is considered the most cosmopolitan in South America, owing to its population density of more than 7,000 inhabitants per square kilometre and the very large number of registered immigrants .
To characterize the genetic diversity of the Mycobacterium tuberculosis strains responsible for TB in the population of São Paulo city, the genotyping techniques of spoligotyping and MIRU were applied to 93 strains isolated from clinical samples taken over two consecutive years from 93 different TB patients residing in this city, who were attending the Clemente Ferreira Institute for treatment.
Spoligotyping has been widely used in the molecular epidemiology of TB, to investigate the population structure of M. tuberculosis, focusing on the identification of genotypic lineages and spoligotype families and their geographic distribution . Furthermore, molecular techniques such as MIRU (Mycobacterial Interspersed Repetitive Units) can be used to cluster M. tuberculosis strains from patients with epidemiological links, identify chains of transmission and discriminate isolates at the clonal level . Spoligotyping is based on DNA polymorphisms within the direct repeat (DR) locus of M. tuberculosis, assaying the presence or absence of a set of target spacer sequences between the DRs. The binary pattern of spacers between the conserved DRs in the region is used to differentiate strains of M. tuberculosis . MIRU genotyping is a highly reproducible and fast system, generating genotypes based on a detailed study of 12 loci of the M. tuberculosis genome, with structurally similar or, in some loci, identical Variable Number of Tandem Repeat (VNTR) mini-satellites . MIRU genotyping applied as a first line molecular typing method, in combination with spoligotyping, has proved to be adequate for strain discrimination in most cases, including large-scale studies .
This paper is a report on the genetic diversity of M. tuberculosis isolates collected over 2 years from patients in São Paulo city, in southeast Brazil, who attended the Clemente Ferreira Institute.
From August 2006 to July 2008 a total of 93 M. tuberculosis isolates were obtained from sputum specimens of patients with pulmonary TB admitted to the Clemente Ferreira Institute, a health reference center for TB treatment in São Paulo city. All isolates were cultured on Löwenstein-Jensen (LJ) slants at 37°C, identified as M. tuberculosis by phenotypic  and PCR-IS6110  methods and tested for drug susceptibility by the Becton-Dickinson MGIT 960 system . The study had ethical approval from the Research Ethics Commitee of the School of Pharmaceutical Sciences, São Paulo State University at Araraquara (reference 01/2006).
Briefly, M. tuberculosis cells from LJ slant cultures were lysed by freezing and boiling three times and the DNA (in solution) was purified by addition of phenol-chloroform-isoamyl alcohol (25:24:1, v/v), followed by chloroform-isoamyl alcohol (24:1, v/v) and then precipitated with ice-cold ethanol. The purified DNA pellet was resuspended in Tris-EDTA, pH 8.0 (TE buffer), and stored at -20°C until used .
The extracted DNA was subjected to PCR, to amplify the DR region, using 1 μL of mycobacterial DNA (10 ng) in 24 μL of a reaction mixture containing 0.4 μM of each primer, DRa 5'-GGTTTTGGGTCTGACGAC-3' (biotinylated 5' end) and DRb 5'-CCGAGAGGGGACGGAAAC-3' (both from Bioneer Company, South Korea), and the PCR Master Mix (Fermentas, USA). PCR products were hybridized with a set of 43 spacer oligonucleotides covalently linked to the spoligo-membrane (Pall Biosystems, Portsmouth, UK), according to the manufacturer's instructions. Bound fragments were incubated with streptavidin-peroxidase conjugate (Boehringer, Ingelheim, Germany) and then detected by chemiluminescence and assessed by an enhanced chemiluminescence system (ECL detection kit, Amersham, Buckinghamshire, UK). The autoradiograms were developed after 30 minutes of exposure, using standard photochemical products .
The term Shared Type (ST) designates spoligotypes common to more than one isolate. These have been assigned "Spoligo-International-Type" (SIT) numbers at the international database SITVIT, allowing online query of the current version in the published SpolDB4 database  and at http://www.pasteur-guadeloupe.fr:8081/SITVITDemo/ webcite. The SITs identified in this study were classified into spoligotype families and subfamilies with the help of the SpolDB4 database. For spoligotype patterns not previously reported in SpolDB4, the 'Spotclust' database was used . This model takes into account knowledge of the evolution of the DR region and assigns spoligotype patterns to families and subfamilies, by means of a computer algorithm based on studies of the previous database SpolDB3 http://cgi2.cs.rpi.edu/~bennek/SPOTCLUST.html webcite.
The mycobacterial DNA was extracted by thermolysis . For genotyping, the set of VNTRs, consisting of 12 MIRUs (MIRU loci 2,4,10,16,20,23,24,26,27,31,39,40), was investigated as described in the literature [7,17]. The PCR amplicon was subjected to electrophoresis in 2.0% w/v agarose gel (Invitrogen Life Technologies, São Paulo, Brazil). DNA ladders of 50 and 100 bp (Invitrogen Life Technologies, São Paulo, Brazil) were used as molecular markers. The gels were stained with ethidium bromide and visualized under ultraviolet light and photodocumented with Alpha-imager 2200 (Alpha Innotech Corporation, San Leandro, CA, USA). The size of the PCR fragment for each locus was estimated by visual comparison with the molecular markers and the MIRU allele score was determined from the size as described [7,16]. The results from each of the 12 loci were combined to create a 12-digit allelic profile.
Construction of dendrograms
The spoligotyping data, the MIRU data and a combination of both data sets were analysed with the software program BioNumerics 4.5 (Applied Maths®, Sint Martens Latem, Belgium), using Euclidean Distance to calculate the similarity matrix, and a dendrogram based on this matrix was constructed by the UPGMA method.
Of the 93 patients affected by tuberculosis, 56 (60.2%) were males and 37 (39.8%) females.
The drug resistance profiles of the 93 clinical isolates were as follows: 58.1% showed resistance to isoniazid (INH), 50.5% to rifampicin (RIF), 26.9% to ethambutol (EMB) and 34.4% to streptomycin (ST). Isolates simultaneously resistant to INH and RIF represented 49.5% (46) of the total number of isolates, whereas 25.8% (24) were resistant to INH, RIF and EMB and 17.2% (16) to INH, RIF, EMB and ST.
Colony morphology, biochemical characteristics and the presence of a 245 bp spoligotyping amplification product confirmed (in all 93 clinical isolates) the etiological agent M. tuberculosis.
Spoligotyping revealed that, of the 93 isolates, 72 (77.4%) showed spoligotypes described in the international database (SpolDB4) consisting of 36 different patterns, belonging to six families and 15 subfamilies (T1, T1 (T4-CE1), T2, T2-T3, T3, H1, H3, LAM1, LAM2, LAM3, LAM4, LAM5, LAM6, LAM9, X3, Beijing, S) and undesignated (U and U likely H) (Table 1).
Table 1. Spoligotypes identified in SpolDB4 by a SIT numbera
Among the 72 isolates with spoligotypes described in SpolDB4, the Latin American-Mediterranean (LAM) and T families were the main spoligotype families, with 26 (36.1%) and 24 (33.3%) isolates respectively, followed by H, with 11 (15.3%) isolates. Thus, the LAM and T spoligotype families together totalled 69.4% of the all isolates. The Undesignated group (U) comprised 6.9% of isolates and the S and X families were each represented by a single isolate (Figure 1).
Figure 1. Dendrogram produced with the results obtained by spoligotyping.
Twenty-one isolates (22.6%), representing 17 spoligotype patterns, could not be classified into SITs and families. These unidentified spoligotypes were analyzed by Spotclust, a program enabling the percent similarity of these patterns with existing subfamilies to be determined (Table 2).
Table 2. Spoligotypes analyzed by Spotclust.
The spoligotype dendrogram constructed with BioNumerics 4.5, presented in Figure 1, grouped 51 (54.8%) of the isolates into 13 clonal groups (two of them composed of unidentified clustered strains) with 100% similarity. Six groups consisted of two isolates each, four groups of four isolates each, one of seven isolates, one of eight isolates and one of ten isolates. Regarding the two clusters of unidentified spoligotype patterns, the dendrogram indicated that the spoligotype pattern with four clustered isolates showed 93% similarity with SIT 50 (subfamily H3), differing only by the absence of spacers 20 and 21. The other spoligotype pattern, with two clustered isolates, showed 87% similarity with SIT 130 (subfamily LAM3), differing in just two spacers (absence of spacer 16 and presence of spacer 39).
Genotyping by the MIRU technique generated 93 distinct genotypes and no formation of genetic groups with 100% similarity. The technique also revealed that the MIRU loci 10, 23, 26 and 40 were the most discriminatory. Loci 26 and 40 were the most polymorphic, with eight alleles, followed by loci 23 and 10 with seven alleles each. In all the strains tested, locus 24 showed 0 or 1 allele. Amplification was unsuccessful for locus 23 in one isolate and locus 40 in another, even after five repeats.
The dendrogram resulting from the combined analysis of spoligotyping and MIRU data, performed with the Bionumerics 4.5 program, showed no clonal group, thus generating 93 independent genotypes (data not shown).
Table 3 shows the isolates clustered by spoligotyping and their respective geographic distribution, MIRU and drug resistance profiles against isoniazid (INH), rifampicin (RIF), streptomycin (ST) and ethambutol (EMB).
Table 3. Summarized results of all assays
The evaluation of epidemiological data from 93 patients revealed that 60.2% were male and 39.8% were female. Moreover, 46 (49.5%) of the M. tuberculosis clinical isolates were multi-drug resistant (MDR). In the same health center, it was found that, of 182 patients, 112 (61.5%) were male and that all of them presented MDR tuberculosis .
The main spoligotype families in this study, LAM and T (representing 69.4% of the isolates), were also found in the Hospital Fernando Fonseca in Lisbon, Portugal, where the predominance of the LAM family was noted in 51% of the isolates . The preponderance of these M. tuberculosis families in São Paulo city and Lisbon is probably due to the close relationship between Brazil and Portugal, especially between the two cities. In Brazil, studies carried out in Rio Grande do Sul  and in Maringa (PR)  corroborate the prevalence of these families. The families LAM and T and also Haarlem are most frequently found in Africa, South America and Central Europe .
The other families found in isolates from São Paulo were the Beijing (four isolates), S and X, with one isolate each. The presence of the Beijing family (SIT1), representing 5.6% of all the identified SITs, may be due to the immigration of Asians to São Paulo city. In Portugal, the Beijing spoligotypes are less frequent (3%) . Within the subfamily T1, we can highlight two isolates of a rare genotype SIT1905. Only two isolates of this genotype are cited in SpolDB4 (one in Brazil and the other in Argentina).
Analysing the geographic distribution of the spoligotypes, it appears that SIT1895 is a localized genotype described only in Brazil and that SIT 828 is found only in Brazil and Guinea Bissau. Fourteen spoligotypes (SITs: 33, 93, 157, 1154, 86, 102, 453, 467, 1166, 1214, 1475, 34, 534 and 560) not previously reported in Brazil were also observed. SIT 33 is a spoligotype identified in South Africa, but the presence of this pattern in São Paulo city may not be surprising. The north-eastern and south-eastern regions of Brazil have a strong African influence because of the slave trade in the past and, after its abolition, immigration from that continent. SIT93 is a spoligotype found in Venezuela and SIT 1154 in Haiti and the presence of these SITs in the city of São Paulo may be related to contact with those countries. The spoligotypes with SITs 86, 102, 453, 1166, 1214 and 1475 of the T1 subfamily, SIT 157 of the T3 as well as SIT 34 of the S family, originate in countries that have close immigration links to Brazil. In the mid 19th century, the southeast region (especially São Paulo) was the destination of many European immigrants, including Portuguese, Italians and Spaniards, in addition to immigrants from the Middle East, such as Syrians and Lebanese.
Analyzing the stability of the genotypes, our study indicated that six SITs (17, 60, 42, 53, 50 and 1) comprised 51.8% of all the isolates whose genotypes were identified in SpolDB4, indicating these as the major genotypes responsible for tuberculosis in São Paulo. In southern Brazil, it was observed that about half of the isolates belonged to seven spoligotypes (SITs 17, 20, 33, 42, 50, 53 and 65) . Eight SITs responsible for 50% of the tuberculosis cases were also found in Lisbon, Portugal . In Honduras, the five predominant spoligotypes were SIT33,42,67,53 and 376 . The prevalence of 13.8%, 11.1% and 9.7% respectively of SITs 53, 50 and 42 was to be expected since these SITs are distributed worldwide . SIT 53 is the second most frequently isolated type in SpolDB4, suggesting that this strain is easily transmitted and adaptable . SITs 17 and 60 are also widely distributed around the world, although predominantly in Latin America [6,25]. In Rio Grande do Sul (south Brazil), the two patterns most frequently found were SIT50 and SIT53 . Despite the fact that six SITs represent more than 50% of all spoligotypes identified in São Paulo, the other SITs found show the great genetic diversity of M. tuberculosis isolates, possibly reflecting the large ethnic mix of the São Paulo city inhabitants. Our data indicate that spoligotyping overestimated cluster formation, in agreement with the data in a previously published study  that investigated the genetic diversity of M. tuberculosis in southern Brazil, a region with a high incidence of tuberculosis.
The MIRU technique generated 93 independent genotypes and did not create any groups with 100% similarity, corroborating another study  in which it was shown that spoligotyping had a lower discriminatory power than MIRU. VNTR (minisatellites) are free of selective or functional pressure, so that MIRU profiles can be stable for up to 18 months of consecutive isolation . Despite the two-year follow-up and study of 93 strains, our data can only show a tendency, owing to the small sample size in comparison with the high tuberculosis incidence in São Paulo. However, when applied to M. tuberculosis isolates from a rural population in Kanpur, South Asia, MIRU also failed to reveal genetic group formation . The authors concluded that even in small populations, M. tuberculosis isolates may present high genetic diversity.
There is a hierarchy among the twelve polymorphic MIRU loci and it has been shown that loci 10, 23, 26, 31 and 40 generate the highest allelic diversity . In agreement with that result, in our study loci 10, 23, 26, and 40 were the most discriminatory. However, locus 31 did not show high variation (ranging from alleles 1 to 5). In another study , loci 10, 23, 26 and 40 also showed the greatest discriminatory power.
MIRU locus 24 enables the classification of spoligotyping families into two groups (A and B). Group B has 0 or 1 allele at locus 24 and consists of families of more recent evolution. In our study, all isolates showed 0 or 1 allele at locus 24 and none was found to belong group A, with ancestral families (Bovis, AFRI or EAI) . The isolates found in São Paulo thus belong to families that have evolved more recently and that are apparently more prevalent around the world.
The combination of the two techniques, spoligotyping and MIRU, separated all isolates of the 13 spoligotyping clusters, generating 93 different genotypes. Despite the small number of isolates in each of the spoligotyping clusters, certain MIRU loci were more discriminatory in differentiating SITs and subfamilies. It was observed that MIRU 40 was the most important in the discrimination of LAM, T and X families, and MIRU 26 and MIRU 23 in that of Beijing and Haarlem, respectively. Looking specifically at the 13 SITs clustered by spoligotyping (Table 3), our results agree with the literature with respect to T1 (locus 40) and Beijing (locus 26 for SIT1). However, for SITs 60 and 828 (LAM4), 42 (LAM9) and 17 (LAM2), locus 23 was the most promising, whereas for SIT50 from the Haarlem family, the most important was locus 10 .
From Table 3, we can also see that all isolates with the SIT1 pattern (Beijing family) were MDR strains, except for one isolate resistant only to streptomycin. The Beijing strains may be genetically more invasive and adapted to acquire drug resistance . However, this study found that all the families identified comprised both MDR and non-MDR isolates. In Brazil, the characteristic of MDR in M. tuberculosis is mostly "acquired" as a result of treatment failure, due to irregularity in the taking of medication, neglect and incorrect prescriptions . Our data corroborated these findings, indicating that in São Paulo city M. tuberculosis clinical isolates have acquired resistance to INH and RIF independently of their genotypes. MDR strains may have been generated by the selective pressure due to problems in adherence to treatment and environmental circumstances (high population density and ethnic diversity) of the inhabitants of São Paulo city.
The authors declare that they have no competing interests.
NM carried out molecular genetic studies, comprising spoligotyping, and drafted the manuscript. FM carried out the sampling of Mycobacterium tuberculosis and sent us the isolates. AS carried out the molecular techniques spoligotyping and MIRU-VNTR. JP, LE and FJ carried out the MIRU-VNTR technique. EA carried out the classic epidemiological analysis and SL the resistance profile of the samples. RF carried out the spoligotyping technique and helped in drafting this manuscript. HB helped in a construction of the dendrogram. SD carried out a draft of this manuscript. CL drafted the manuscript and is advisor and head of the laboratory.
All authors read and approved the final manuscript
Acknowledgements and Funding
This study was supported by FAPESP (2009/53292-3) and CNPq (472867/2006-7). We thank careful proofreading of the text by Timothy John C. Roberts (M.Sc., London).
Kamerbeek J, Schouls L, Kolk A, van Agterveld M, van Soolingen D, Kuijper S, Bunschoten A, Molhuizen H, Shaw R, Goyal M, van Embden J: Simultaneous detection and strain differentiation of Mycobacterium tuberculosis for diagnosis and epidemiology.
Cowan LS, Diem L, Monson T, Wand P, Temporado D, Oemig TV, Crawford JT: Evaluation of a two-step approach for large-scale, prospective genotyping of Mycobacterium tuberculosis isolates in the United States.
Sola C, Filliol I, Gutirrez MC, Mokrousov I, Vincent V, Rastogi N: Spoligotype database of Mycobacterium tuberculosis: biogeografic distribution of shared types and epidemiologic and phylogenetic perspectives.
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 2008, 46:1398-1406.
2008PubMed Abstract | Publisher Full Text | PubMed Central Full Text
Leão S, Martin A, Megia GI, Palomino JC, Robledo J, Telles MAS, Portales F: Pratical handbook for the phenotypic and genotypic identification of mycobacteria. [http://www.ESMycobacteriology.eu] webcite
Brugges, Vanden Broelle; :164.
van Embden JD, Cave MD, Crawford JT, Dale JW, Eisenach KD, Gicquel B, Hermans P, Martin C, McAdam R, Shinnick TM, Small PM: Strain identification of Mycobacterium tuberculosis by DNA fingerprinting: recommendations for a standardized methodology.
Brudey K, Driscoll JR, Rigouts L, Prodinger WM, Gori A, Al-Hajoj SA, Allix C, Aristimuno L, Arora J, Baumanis V, Binder L, Cafrune P, Cataldi A, Cheong S, Diel R, Ellermeier C, Evans JT, Fauville-Dufaux M, Ferdinand S, Garcia de Viedma D, Garzelli C, Gazzola L, Gomes HM, Gutierrez MA, Hawkey PM, van Helden PD, Kadival GV, Kreiswirth BN, Kremer K, Kubin M, Kulkarni SP, Liens B, Lillebaek T, Ly HM, Martin C, Moskrousov I, Narvskaïa O, Ngeow YF, Naumann L, Niemann S, Parwati I, Rahim Z, Rasolofo- Razanamparany V, Rasolonavalona T, Rossetti ML, Rüsch-Gerdes S, Sajduda A, Samper S, Shemyakin IG, Singh US, Somoskovi A, Skuce RA, van Soolingen D, Streicher EM, Suffys PN, Tortoli E, Tracevska Tatjana, Vincent V, Victor TC, Warren RM, Yap SF, Zaman K, Portaels F, Rastogi N, Sola C: Mycobacterium tuberculosis complex genetic diversity: mining the fourth international spoligotyping database (SpolDB4) for classification population genetics and epidemiology.
Mazars E, Lesjean S, Banuls AL, Gilbert M, Vincent V, Gicquel B, Tibayrenc M, Locht C, Supply P: High-resolution minisatellite-based typing as a portable approach to global analysis of Mycobacterium tuberculosis molecular epidemiology.
Rev Soc Bras Med Trop 2003, 36:27-34. PubMed Abstract
Silva ABS, Groll AV, Feliz C, Conceição FR, Spies FS, Scaini CJ, ROSSETTI ML, Borsuk S, Dellagostin OA, Silva PEA: Clonal diversity of Mycobacterium tuberculosis isolated in a sea port city in Brazil.
Borsuk S, Dellagostin MM, Madeira Sde G, Lima C, Boffo M, Mattos I, Almeida da Silva PE, Dellagostin AO: Molecular characterization of Mycobacterium tuberculosis isolates in a region of Brazil with a high incidence of tuberculosis.
Microbes Infec 2005, 7:1338-1344. Publisher Full Text
Malaspina AC, Cavalcanti HR, Leite CQ, Machado SM, Viana BH, Silva RM, Hage EF, Figueiredo WM, Marques E, Ferrazoli L, Arbex M, Lessi M, Fonseca LS, Rigouts L, Saad MH: Usefulness of Mycobacterium tuberculosis molecular typing in a tuberculosis low-endemic agro-industrial setting of Brazil.
Filliol I, Driscoll JR, van Soolingen D, Kreiswirth BN, Kremer K, Valétudie G, Anh DD, Barlow R, Banerjee D, Bifani PJ, Brudey k, Cataldi A, Cooksey RC, Cousins DV, Dale JW, Dellagostin OA, Drobniewski f, Engelmann G, Ferdinand S, Binzi DG, Gordon M, Gutierrez C, Haas WH, Heersma H, Kassa-Kelembho E, Ly HO, Makristathis A, Mammina C, Martin G, Mostro¨m P, Mokrousov I, Narbonne V, Narvskaya O, Nastasi A, Niobe-Eyangoh SN, Pape JW, Rasolofo-Razanamparany V, Ridell M, Rossetti ML, Stauffer F, Philip NS, Takiff H, Texier-Maugein J, Vicent V, Waard JH, Sola C, Rastogi N: Global distribution of Mycobacterium tuberculosis spoligotypes.
Sharma P, Chauhan DS, Upadhyay P, Faujdar J, Lavania M, Sachan S, Katoch K, Katoch VM: Molecular typing of Mycobacterium tuberculosis isolates from a rural area of Kanpur by spoligotyping and mycobacterial interspersed repetitive units (MIRUs) typing.
Ferdinand S, Valétudie G, Sola C, Rastogi N: Data mining of Mycobacterium tuberculosis complex genotyping results using mycobacterial interspersed repetitive units validates the clonal structure of spoligotyping-defined families.
Morcillo N, Di Giulio B, Chirico C, Kuriger A, Dolmann A, Alito A, Zumárraga M, van Soolingen D, Kremer K, Cataldi A: Primera descripción del genotipo Beijing de Mycobacterium tuberculosis en la Argentina.
Rev Argent Microbiol 2005, 37:92-95. PubMed Abstract