Genotyping of Mycobacterium Tuberculosis Isolates From A Low-Endemic Setting in Northwestern State of Paraná in Southern Brazil
Genotyping of Mycobacterium Tuberculosis Isolates From A Low-Endemic Setting in Northwestern State of Paraná in Southern Brazil
Genotyping of Mycobacterium Tuberculosis Isolates From A Low-Endemic Setting in Northwestern State of Paraná in Southern Brazil
The purpose of this study was to provide information about the genetic diversity and prevalent genotype of My-
cobacterium tuberculosis in a low-endemic setting in northwestern state of Paraná in southern Brazil. We employed
spoligotyping and mycobacterial interspersed repetitive units-variable number tandem repeat (MIRU-VNTR) tech-
niques to genotype M. tuberculosis isolates from patients with pulmonary tuberculosis (TB). The 93 isolates ana-
lyzed by spoligotyping were divided into 36 different patterns, 30 of which were described in the SITVIT database.
Latin American and Mediterranean, Haarlem and T families were responsible for 26.9%, 17.2% and 11.8% of TB
cases, respectively. From the 84 isolates analyzed by MIRU-VNTR, 58 shared a unique pattern and the remaining 26
belonged to nine clusters. The MIRU loci 40, 23, 10 and 16 were the most discriminatory. A combination of MIRU-
VNTR and spoligotyping resulted in 85.7% discriminatory power (Hunter-Gaston index = 0.995). Thus, combining
spoligotyping and MIRU-VNTR typing proved to be most useful for epidemiological study in this low-endemic setting
in southern Brazil. The current study demonstrated that there is significant diversity in circulating strains in the city
of Maringá and the surrounding regions, with no single genotype of M. tuberculosis predominating.
Tuberculosis (TB) is the leading cause of death by an isolates of the same strains. Furthermore, the typing
infectious agent and thus an important public health prob- methods employed must be reproducible, discriminatory
lem worldwide, with an estimated 9.27 million new cases and easy to perform.
in 2007. In Brazil, the estimated incidence rate is 48 cases Some polymerase chain reaction (PCR)-based tech-
per 100,000 inhabitants (WHO 2009). Brazil has five dis- niques are now being used to differentiate Mycobacte-
tinct geographic regions, among which the distribution rium tuberculosis isolates. Currently, genotyping ap-
of TB cases varies greatly. In the state of Paraná (PR) in proaches targeting the analysis of the variable number of
southern Brazil, the average rate of TB is 27.52 cases per tandem repeats (VNTR), based on mycobacterial inter-
100,000 inhabitants (Malaghini et al. 2009). Maringá and spersed repetitive units (MIRU), are the most promising.
the small cities surrounding it form a low-endemic area This technique is based on the variability found at 12
of TB in PR, with an incidence rate of 22.62 cases per specific loci interspersed throughout the mycobacterial
100,000 in 2007 (personal communication). genome (Supply et al. 2000). Recently, MIRU-VNTR
Genotyping methods have been applied extensively genotyping strategies using 15 or 24 loci (Supply et al.
in the epidemiological study of TB worldwide (Candia et 2006, Oelemann et al. 2007a, Alonso-Rodríguez et al.
al. 2007, Sharma et al. 2008, Valcheva et al. 2008). These 2008) were evaluated and applied to molecular epide-
studies have been based on the assumption that patients miological typing in mycobacteria.
with genotypically clustered strains are epidemiologi- Spoligotyping is the second most widely used meth-
cally linked and represent recent transmissions. In con- od for M. tuberculosis complex genotype after IS6110-
trast, patients infected with different types of strains are based fingerprinting (Prodinger 2007). It is based on the
not considered indicative of recent transmission. Genetic presence or absence of a set of target sequences in the
markers need to be both sufficiently polymorphic to dis- direct repeat (DR) locus in the M. tuberculosis complex
tinguish unrelated strains and stable enough to identify genome (Kamerbeek et al. 1997) that in combination
with MIRU has been used to replace typing via restric-
tion fragment length polymorphisms (RFLP) based on
the insertion sequence IS6110. The IS6110 method has
been considered the “gold standard” for genotyping M.
tuberculosis since 1993, but it is an expensive, labori-
Financial support: CAPES
ous and lengthy methodology that requires weeks of M.
+ Corresponding author: [email protected] tuberculosis culturing and specific software to analyze
Received 11 January 2010 the RFLP band-patterns, all of which make it difficult
Accepted 13 May 2010 to interpret and exchange data. In addition, the IS6110
online | memorias.ioc.fiocruz.br
780 Genotyping of M. tuberculosis isolates from Brazil • Erika Noda Noguti et al.
method has limited use in genotyping M. tuberculosis were detected by chemiluminescence after incubation
isolates that contain few copies of IS6110. with a streptavidin-peroxidase conjugate (Boehringer
The aim of our study was to use spoligotyping and Ingelheim, Germany) and assessed by an enhanced
MIRU-VNTR typing to provide preliminary informa- chemiluminescence system (GE Healthcare UK Limit-
tion about the genetic diversity and prevalent genotype ed, Buckinghamshire, UK). Spoligotypes were reported
of M. tuberculosis isolates in a low-endemic setting in using an octal code in which the 43-digit binary repre-
northwestern PR. senting the 43 spacers (“1” is hybridization and “0” is
PATIENTS, MATERIALS AND METHODS
no hybridization) was divided into 14 sets of three digits
(spacers 1-42) plus one additional digit (spacer 43). Each
Study isolates - From November 2005-June 2008, a three-digit set was converted to octal code (000 = 0, 001
total of 93 M. tuberculosis isolates were obtained from = 1, 010 = 2, 011 = 3, 100 = 4, 101 = 5, 110 = 6 and 111 =
sputum of patients suspected of TB who were seen at the 7), with the final digit remaining either 1 or 0, yielding a
Clinical Bacteriology Laboratory, Department of Clini- 15-digit octal designation (Dale et al. 2001).
cal Analysis of the State University of Maringá. This lab
is a reference TB laboratory that handles patients from MIRU-VNTR typing - MIRU-VNTR typing was
Maringá and other cities in northwestern PR. All isolates performed in only 84 of the 93 M. tuberculosis isolates
were cultured in DifcoTM Lowenstein Medium Base (Bec- at the Laboratory of Mycobacteria Dr. Hugo David in
ton, Dickinson and Company, Sparks, MD, USA) and the School of Pharmaceutical Science, São Paulo State
identified as M. tuberculosis using the conventional bio- University, Araraquara, São Paulo (SP), Brazil. MIRU-
chemical tests (Kent & Kubica 1985) and molecular biol- VNTR was not performed in nine of the isolates because
ogy (van Embden et al. 1993). The following retrospective there was insufficient DNA for analysis.
demographic and epidemiological data were collected for The isolates were genotyped by PCR amplification
all patients through a review of a national TB notification of the original 12 MIRU-VNTR loci (2, 4, 10, 16, 20, 23,
database (National Diseases Notification System): city 24, 26, 27, 31, 39, 40) as described by Supply et al. (2000)
and zip code of residence at the time of diagnosis, age, and Mazars et al. (2001) in a PTC-100 thermal cycler (MJ
sex, ethnicity, HIV status, sample susceptibility profile, Research, Ramsey, Minnesota, USA). Each locus was
alcoholism and the occurrence of other diseases. amplified individually from 2 µL of mycobacterial DNA
(20 ng) in 23 µL of a reaction mixture containing 0.4 µM
DNA extraction - DNA from M. tuberculosis was of loci-respective primers and PCR Master Mix (Pro-
extracted from a subculture on DifcoTM Lowenstein Me- mega Corporation, Madison, Wisconsin, USA) according
dium Base, as described by González-y-Merchand et al. to the manufacturer’s instructions. The PCR conditions
(1996), with minor modifications. Briefly, a loopful of for each set of primers have been described elsewhere
bacteria was suspended in 6 M guanidine hydrochloride (Pandolfi 2006). PCR products were subjected to electro-
(Sigma Chemical Co, St. Louis, MO, USA) and bacilli phoresis in a 2% weight/volume agarose gel (Invitrogen
were lyzed by freezing at -20°C for 30 min followed by Life Technologies, SP, Brazil). Both 50 and 100-bp DNA
heating at 100ºC for 10 min. This procedure was repeat- Ladders (Invitrogen Life Technologies, SP, Brazil) were
ed twice. DNA was further extracted with two volumes used as molecular markers. The gels were stained with
of phenol-chloroform-isoamyl alcohol (25:24:1, v/v) ethidium bromide and visualized under ultraviolet light,
followed by two steps of extractions with chloroform- then photodocumented with an Alpha-imager 2200 (Al-
isoamyl alcohol (24:1, v/v). DNA was purified by etha- pha Innotech Corporation, San Leandro, CA, USA). PCR
nol precipitation, dissolved in 50 µL of Tris-ethylenedi- fragment size was determined by visual comparison with
amine tetraacetic acid (TE buffer), pH 8.0, and stored at the molecular markers and the MIRU allele scoring was
-20ºC until use. DNA concentration was determined by determined according to Mazars et al. (2001) and Supply
ultraviolet spectrophotometry. et al. (2001). The results from each of the 12 loci were
Spoligotyping - Spoligotyping was performed in all combined to create the 12-digit allelic profiles.
93 M. tuberculosis isolates, using the standard method, Interpretation of genotyping results - The observed
to detect the presence or absence of 43 spacers (Mo- spoligotypes were compared to the international data-
lhuizem et al. 1998). Briefly, the DR region was ampli- base, SITVIT, which is an updated version of the pub-
fied from 1 µL of mycobacterial DNA (10 ng) in 24 µL lished SpolDB4.0 database (Brudey et al. 2006) and it is
of a reaction mixture containing 0.4 µM of each primer available at pasteur-guadeloupe.fr:8081/SITVITDemo/.
and PCR Master Mix (Promega Corporation, Madison, For previously unreported spoligopatterns in the Spol-
WI, USA) according to manufacturer’s instructions. The DB4.0, the “Spotclust” database was used (Vitol et al.
primer sequences were Dra, 5’-GGTTTTGGGTCT- 2006). This model takes into account knowledge of the
GACGAC-3’ (biotinylated 5’ end) and DRb, 5’-CCGA- evolution of the DR region and assigns spoligopatterns
GAGGGGACGGAAAC-3’ (Integrated DNA Technolo- to families and subfamilies, using a computer algorithm
gies, Inc Coralville, USA). The DNA amplification was based on studies of SpolDB3 (cgi2.cs.rpi.edu/~bennek/
carried out in a TC-512 thermal cycler (Techne, UK). SPOTCLUST.html). BioNumerics software (version
PCR products were hybridized with a set of 43 spacer 4.45; Applied Maths, Sint-Martens-Latem, Belgium) was
oligonucleotides covalently linked to the spoligo-mem- used for analysis of spoligotyping and MIRU-VNTR pat-
brane (Isogen Life Sciences, The Netherlands) accord- terns. Dendrograms were constructed for spoligotyping,
ing to the manufacturer’s instructions. Bound fragments MIRU-VNTR and the combination of both methodolo-
Mem Inst Oswaldo Cruz, Rio de Janeiro, Vol. 105(6), September 2010 781
TABLE I
Spoligopatterns and family assignment of 10 isolates of Mycobacterium tuberculosis not described in SITVIT database
ID: identification; LAM: Latin American and Mediterranean; SND: spoligotype not described; T: modern tuberculosis strains.
782 Genotyping of M. tuberculosis isolates from Brazil • Erika Noda Noguti et al.
and 11 (15.3%). The minor families observed in our study Allele polymorphism analysis of 12 MIRU loci re-
were EAI (n = 1, 1.4%), S lineage (n = 1, 1.4%) and H1-S (n vealed that MIRU locus 40 was the most discriminatory
= 1, 1.4%). Eleven M. tuberculosis isolates were not clas- locus with eight alleles, followed by MIRU loci 23, 10
sified into families according to SpolDB4, but did have and 16. In MIRU locus 40, the presence of three alleles
attributed ST numbers according to the SITVIT database was most frequent, followed by a single copy of an al-
(ST2508, ST2512, ST2525, ST2563 and ST2654). lele. The MIRU 20, 26 and 31 loci were moderately dis-
Ten isolates (10.8%), comprising six spoligotyping criminant. Other loci were less polymorphic, with only
patterns not found in the SITVIT database, were as- two alleles for MIRU loci 4, 20 and 39 and one allele for
signed to families and subfamilies by “SpotClust”. The MIRU locus 24, for which a single copy was present in
family assignment of these remaining isolates revealed all 84 M. tuberculosis isolates analyzed (Table II).
that 80% belonged to the LAM family and 20% belonged
Combining spoligotyping and MIRU-VNTR typing -
to the T family (Table I). These six new spoligotypes
Among the 84 M. tuberculosis isolates analyzed by spoligo-
have been submitted to the SITVIT database (pasteur-
typing and MIRU-VNTR typing, 72 (85.7%) distinct geno-
guadeloupe.fr:8081/SITVITDemo/).
types were obtained, resulting in eight clusters with 100%
MIRU-VNTR typing - MIRU-VNTR typing was con- similarity. By considering similarity indices of at least 69%,
ducted for 84 clinical isolates and a total of 67 distinct we observed three distinct clonal groups that represented
MIRU patterns were obtained. Of the 84 isolates, 58 (69%) 92.9% of the isolates analyzed (Figure). Table III summa-
were orphans and the remaining 26 (31%) were grouped rizes the discriminative power of each typing method alone
into nine clusters, each comprising 2-6 isolates (Figure). and of the two methods combined, according to the HGI.
TABLE II
Allelic polymorphism of 12 mycobacterial interspersed repetitive units (MIRU) loci from the 84 Mycobacterium tuberculosis
isolates from patients with tuberculosis in Maringá and region in the northwest of the state of Paraná, Brazil
Allele number
the allelic diversity of the loci was classified as highly discriminant [Hunter-Gaston index (HGI) > 0.6], moderately discriminant
(0.3 ≤ HGI ≤ 0.6) and poorly discriminant (HGI < 0.3), according to Sola et al. (2003).
TABLE III
Discriminatory power of spoligotyping and mycobacterial interspersed repetitive units (MIRU) typing, alone and in association
Supply P, Lesjean S, Savine E, Kremer K, van Soolingen D, Locht van Embden JD, Cave MD, Crawford JT, Dale JW, Eisenach KD,
C 2001. Automated high-throughput genotyping for study of Gicquel B, Hermans P, Martin C, McAdam R, Shinnick TM,
global epidemiology of Mycobacterium tuberculosis based on Small PM 1993. Strain identification of Mycobacterium tuber-
mycobacterial interspersed repetitive units. J Clin Microbiol culosis by DNA fingerprinting: recommendations for a stan-
39: 3563-3571. dardized methodology. J Clin Microbiol 31: 406-409.
Supply P, Mazars E, Lesjean S, Vincent V, Gicquel B, Locht C 2000. Vitol I, Driscoll J, Kreiswirth B, Kurepina N, Bennett KP 2006.
Variable human minisatellite-like regions in the Mycobacteri- Identifying Mycobacterium tuberculosis complex strains fami-
um tuberculosis genome. Mol Microbiol 36: 762-771. lies using spoligotypes. Infect Genet Evol 6: 491-504.
Valcheva V, Mokrousov I, Rastogi N, Narvskaya O, Markova N WHO - World Health Organization 2009. Global tuberculosis con-
2008. Molecular characterization of Mycobacterium tuberculo- trol: epidemiology, strategy, financing. WHO report WHO/
sis isolates from different regions of Bulgaria. J Clin Microbiol HTM/TB/2009.411). Available from: http://www.who.int/tb/
46: 1014-1018. publications/global_report/2009/en/index.html.