Tuberculosis Research and TreatmentVolume 2011 (2011), Article ID 239042, 9 pagesdoi:10.1155/2011/239042

Research Article

Mycobacterium tuberculosis Complex Genotype Diversity and Drug Resistance Profiles in a Pediatric Population in Mexico

Mercedes Macías Parra,¹ Jesús Kumate Rodríguez,² José Luís Arredondo García,³  Yolanda López-Vidal,⁴ Mauricio Castañón-Arreola,⁴ Susana Balandrano,⁵ Nalin Rastogi,⁶ and Pedro Gutiérrez Castrellón1,3

¹Department of Infectious Diseases, Instituto Nacional de Pediatría, Mexico City, Mexico²Instituto Mexicano del Seguro Social Foundation, Mexico City, Mexico³Subdivision of Medical Research, Instituto Nacional de Pediatría, Mexico City, DF, Mexico⁴Programa de Inmunología Molecular Microbiana, Departamento de Microbiología y Parasitología, Facultad de Medicina, Universidad Nacional Autónoma de México, Mexico City, Mexico⁵Department of Micobacteriology, Instituto De Diagnóstico y Referencia de Epidemiológicos, Secretaria de Salud, Mexico City, Mexico⁶WHO Supranational Tuberculosis Reference Laboratory, Institut Pasteur de la Guadeloupe, Abymes, France

Received 17 March 2011; Accepted 22 May 2011

Academic Editor: Prasit Palittapongarnpim

Copyright © 2011 Mercedes Macías Parra et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The aim of this study was to determine the frequency of drug resistance and the clonality of genotype patterns in M. tuberculosis clinical isolates from pediatric patients in Mexico ( patients from 19 states; time period—January 2002 to December 2003). Pulmonary disease was the most frequent clinical manifestation (71%). Children with systemic tuberculosis (TB) were significantly younger compared to patients with localized TB infections (mean years versus years ). Resistance to any anti-TB drug was detected in 24/90 (26.7%) of the isolates; 21/90 (23.3%) and 10/90 (11.1%) were resistant to Isoniazid and Rifampicin, respectively, and 10/90 (11.1%) strains were multidrug-resistant (MDR). Spoligotyping produced a total of 55 different patterns; 12/55 corresponded to clustered isolates (, clustering rate of 52.2%), and 43/55 to unclustered isolates (19 patterns were designated as orphan by the SITVIT2 database). Database comparison led to designation of 36 shared types (SITs); 32 SITs ( isolates) matched a preexisting shared type in SITVIT2, whereas 4 SITs ( isolates) were newly created. Lineage classification based on principal genetic groups (PGG) revealed that 10% of the strains belonged to PGG1 (Bovis and Manu lineages). Among PGG2/3 group, the most predominant clade was the Latin-American and Mediterranean (LAM) in 27.8% of isolates, followed by Haarlem and T lineages. The number of single drug-resistant (DR) and multidrug-resistant (MDR-TB) isolates in this study was similar to previously reported in studies from adult population with risk factors. No association between the spoligotype, age, region, or resistance pattern was observed. However, contrary to a study on M. tuberculosis spoligotyping in Acapulco city that characterized a single cluster of SIT19 corresponding to the EAI2-Manila lineage in 70 (26%) of patients, not a single SIT19 isolate was found in our pediatric patient population. Neither did we find any shared type belonging to the EAI family which represents ancestral PGG1 strains within the M. tuberculosis complex. We conclude that the population structure of pediatric TB in our setting is different from the one prevailing in adult TB patient population of Guerrero.

1. Introduction

Although curable, tuberculosis (TB) remains a global public health concern with a major impact in developing countries [1]. Globally, approximately nine million new TB cases and 1.6 million deaths occur every year. The global epidemiology of TB in children is limited; however, the World Health Organization (WHO) estimates that one million cases occurred in children, comprising up to 10% of the total TB incidence [2–4]. Children rarely transmit the disease, and thus contribute little to the maintenance of the TB epidemic even though they make up a substantial proportion of the global disease burden [5].

The increasing number of drug-resistant, and especially multidrug-resistant (MDR: defined as resistance to at least isoniazid and rifampicin), cases is a serious concern for the Global TB Control Program instituted by the WHO [6–9]. Treatment of drug-resistant TB is expensive and complex because it necessitates the use of second-line drugs, which are associated with a greater incidence of adverse events [10]. The true extent of drug resistance is unknown, but it is estimated that the global proportion of resistance among all cases is 4.8% (95% confidence interval (CI) 4.6 to 6.0%), with an estimated 489,000 cases of MDR-TB occurring worldwide in 2006. It is thus suggested that resistance to anti-TB drugs may be increasing in some geographical areas, with considerable variation between different countries [11].

Mexico has a declining, but still moderate, incidence of TB that is variable among states (WHO 2005). In 2007, the reported incidence was 13.49 cases/100,000 individuals, with a range of 2.35 to 35.34 cases/100,000 inhabitants and an incidence rate of 5.077/100,000 in children under 19 years of age (Indicadores Demográficos 1990–2030 Consejo CONAPO Sistema Nacional de Vigilancia Epidemiológica, Dirección General de Epidemiología, Secretaria de Salud). Unfortunately, Mexico’s surveillance system does not provide an accurate estimation of the drug-resistant strains because the diagnosis of TB is typically performed by direct acid fast bacilli sputum smears. However, limited data from culture-based drug susceptibility testing in 2008 indicate that 2.4% of cases are initially MDR-TB (http://www.cenavece.salud.gob.mx/descargas/pdf/tuberculosis.pdf). There are now up to 479 reported cases of MDR-TB, seven of which were in patients under 19 years of age and 47 that were likely extensively drug-resistant TB (XDR-TB). A strain is considered XDR-TB if it is resistant to isoniazid, rifampicin, any fluoroquinolone and at least one of three injectable second-line drugs (amikacin, kanamycin, or capreomycin); thus, these cases are on the verge of being untreatable [12]. MDR-TB and XDR-TB are among the greatest concerns in the antibiotic resistance pandemic due to the high risk of death [12, 13]. Furthermore, patients can remain infectious for months or even years [14] and may spread MDR- or XDR-TB. Consequently, early detection and an accurate record of individuals the patient has had contact with, are critical in arresting further transmission of the disease and for the proper control of TB.

Study of TB molecular epidemiology through DNA fingerprinting is an important tool contributing to the understanding of the transmission and control of TB [15, 16]. PCR-based spacer oligonucleotide typing (spoligotyping) based on the variability of the direct repeat (DR) locus in the M. tuberculosis complex has emerged as a fast, reliable, and cost-effective alternative to the traditional IS6110 restriction fragment length polymorphism (RFLP) for a first-line genotypic screening of tubercle bacilli [17–19]. More recently, this methodology has been used to provide support for lineage-specific differences and global phylogeography of TB in international databases [20]. In this context, we thought it desirable to focus on pediatric TB since childhood TB is a sentinel event, indicating the ongoing transmission of TB. In addition, these patients are more likely to develop severe cases and extrapulmonary diseases than adults. Unfortunately, there is little information regarding the epidemiology and resistance patterns to anti-TB drugs in children [6]. The aim of this study was to determine the genotype diversity and the extent of drug resistance in M. tuberculosis clinical isolates from pediatric patients in Mexico.

2. Methods

2.1. Patients and Specimens

This study included a total of 90 Mycobacterium tuberculosis complex isolates, recruited from patients who were 18 years of age or younger over a 2-year period (January 2002 to December 2003) from 18 federal States in Mexico (see Table  S1 of Supplementary Material available online at doi: 10.1155/2011/239042). Clinical specimen were processes at the local laboratories for culture, followed by identification and drug susceptibility testing (DST) performed at the mycobacteriology laboratory at Instituto De Diagnóstico y Referencia de Epidemiológicos (InDRE) in Mexico city. Isolates were characterized using several biochemical tests: niacin and nitrate positivity; pyrazinamidase activity; negative for catalase activity at 68°C [18]; in-house PCR-RFLP of hsp65; determination of antimicrobial resistance using the semiautomated radiometric BACTEC 460TB system according to the manufacturer’s instructions (Becton Dickinson, Sparks, MD, USA), with following critical concentrations: streptomycin (STM, 1.0 μg/mL), rifampicin (RMP, 1.0 μg/mL), isoniazid (INH, 0.1 μg/mL), pyrazinamide (PZA, 100 μg/mL), and ethambutol (EMB, 5 μg/mL). Basic demographic data was collected for each patient using a standard questionnaire using the files provided at the local hospitals. The protocol was approved by the Research and Ethics Committee of the Instituto Nacional de Pediatria in Mexico City.

2.2. Spoligotyping and Database Comparison

Spoligotyping was performed at the Immunology and Molecular Microbiology Program at the National Autonomous University of Mexico on DNA extracted from heat inactivated cultures using a previously described protocol [19]. Spoligotypes in octal codes were entered in the SITVIT2 proprietary database of the Pasteur Institute of Guadeloupe, which is an updated version of the previously released SpolDB4 database [20]. A Spoligotype International Type (SIT) number was attributed to each pattern according to the SITVIT2 database. At the time of the present study, SITVIT2 contained more than 3000 SITs with global genotyping information on about 75,000 M. tuberculosis complex (MTC) clinical isolates from 160 countries of origin. Worldwide distribution of spoligotypes for all clustered isolates (SITs representing 2 or more strains) was investigated using the SITVIT2 database and recorded for countries and regions representing ≥5% of a given SIT as compared to their total number in the global database. The various macrogeographical regions and subregions were defined per UN specifications (http://unstats.un.org/unsd/methods/m49/m49regin.htm). In this database, SIT (Spoligotype International Type) designates spoligotyping shared by 2 or more patient isolates, as opposed to “orphan” which designates patterns reported for a single isolate. Major phylogenetic clades were assigned according to signatures provided in SpolDB4, which defined 62 genetic lineages/sublineages [20]. These include specific signatures for various MTC members such as M. bovis, M. caprae, M. microti, M. canetti, M. pinipedi, and M. africanum, as well as rules defining major lineages/sublineages for M. tuberculosis sensu stricto; these include the Beijing clade, the Central Asian (CAS) clade and 2 sublineages, the East African-Indian (EAI) clade and 9 sublineages, the Haarlem (H) clade and 3 sublineages, the Latin American-Mediterranean (LAM) clade and 12 sublineages, the ancestral “Manu” family and 3 sublineages, the S clade, the IS6110-low-banding X clade and 3 sublineages, and an ill-defined T clade with 5 sublineages (as well as further well-characterized phylogeographical specificity for 8 additional spoligotype signatures).

2.3. Distinction of “Ancient” versus “Modern” Lineages

We also compared the overall repartition of isolates according to major M. tuberculosis genotypic families by adding all the shared types for each of the individual family defined (as well as the orphan strains), and further linked the information obtained based on the lineage classification to “ancient” and “modern” lineages of tubercle bacilli as defined by principal genetic groups (PGGs) based on KatG463-gyrA95 polymorphism [21], inferred from the reported linking of specific spoligotype patterns to PGG1, 2 or 3 [22–24].

2.4. Statistical Analyses

SPSS for Windows version 13.0 was used for statistical analyses. Association between categorical variables was assessed by Fisher’s exact test and the Mann-Whitney U test was used to evaluate associations between numerical variables. Differences between groups were detected by univariate analyses and expressed as the odds ratio (OR) with 95% confidence intervals (95% CI). Differences were considered to be significant if values were <0.05.

3. Results

3.1. Study Population

This study included a total of 90 M. tuberculosis complex strains isolated from pediatric patients in 18 federal States over a 2-year period (January 2002 to December 2003), which were subsequently referred to the National Institute for Diagnostic and Epidemiological Reference (InDRE) for identification and drug susceptibility testing. These represented between 0.4 and 15.5% of TB cases in this age group as reported by the General Epidemiology Division by those States during the study period. Detailed information on individual strains obtained during this investigation is summarized in Table  S1. Information regarding each strain includes year of isolation, spoligotype description, shared-type number, genotypic lineage, and drug susceptibility to first-line drugs. Additional demographic and clinical information includes pathological specimen, site of infection, geographical origin, age, and sex of the patients.

3.2. Clinical Manifestations

M. tuberculosis was mostly isolated from sputum samples (77%), mainly in patients with localized TB. In contrast, patients with disseminated disease were diagnosed by spinal fluid (8%), gastric aspirates (4%), urine (2%), pleural fluid (1%), or tissue (8%). Male patients (48) comprised 53% of the population; the mean age of all was 13 ± 5 years, 60% were between 15 and 18 years old and only 8% were under two years of age. Localized TB was diagnosed in 76% of the cases and systemic TB diagnoses corresponded to the remaining 24%. The disease frequently presented as pulmonary TB (65 cases, 71.4%) and was followed by meningeal in 8 cases (8.8%), disseminated in 5 cases (5.5%), miliary in 4 cases (4.4%), abdominal in 3 cases (3.3%), ganglionary in 3 cases (3.3%), urinary in 2 cases (2.2%), and pleural in 1 case (1.1%). Patients with systemic TB were significantly younger ( years) than children with localized TB ( years), . Of this population, in 19 of the patients the presence of prior underlying condition was not known. In 10/71 (14%) of the patients, accompanying pathology was identified: 6 (8.4%) were severely malnourished, 2 (2.8%) had a coinfection with human immunodeficiency virus, 1 (1.4%) patient harboring a SIT1/Beijing genotype had diabetes mellitus, and 1 (1.4%) presented with alcoholism. The patients with diabetes mellitus and alcoholism were 18 years of age. The remaining 61 patients had no associated comorbidity.

3.3. Drug Resistance

As summarized in Table 1, M. tuberculosis isolate showed resistance to at least one anti-TB drug in 24/90 (26.7%) patients; 8 (8.8%) were resistant to a single drug and 4 (4.4%) showed resistance to five first-line TB drugs. The greatest resistance was found against INH (23.3% of cases) followed by RMP (11.1% of cases). Resistance to RMP was always accompanied by INH-resistance leading to 10 (11.1%) cases of MDR-TB. Of 71 patients from whom information was available, 14 patients had previous TB treatment, and secondary resistance was documented in 7/14 (50%) patients in contrast with primary resistance found in 14/57 (24.5%) patients without previous TB treatment ().

Table 1: Characteristics of M. tuberculosis drug resistance to first line antituberculosis drugs.

Number of drugs
Number of isolates
Resistance to antituberculosis drugs (no and %)

INH
STM
RMP
PZA
EMB

1
8
5 (23.8)
3 (21.2)
—
—
—

2
6
6 (28.5)
2 (18.1)
2 (20.0)
1 (12.5)
1 (11.1)

3
5
5 (23.8)
2 (18.1)
3 (30.0)
2 (25.0)
3 (33.3)

4
1
1 (4.7)
—
1 (10.0)
1 (12.5)
1 (11.1)

5
4
4 (19.0)
4 (36.3)
4 (40.0)
4 (50.0)
4 (44.1)

Global resistance
24
21 (23.3)
11 (12.2)
10 (11.1)
8 (8.8)
9 (10.0)

INH: Isoniazid, STM: Streptomycin, RMP: Rifampicin, PZA: Pyrazinamide, EMB: Ethambutol.

3.4. Risk Factors

Bivariate analysis of risks identified previous treatment as a risk factor for developing resistant to M. tuberculosis with an OR = 2.03 (95% CI 1.01 to 4.07, ). Although this result is not statistically significant, it is of clinical relevance. Furthermore, the risk increased significantly in MDR-TB cases with an OR = 5.08 (95% CI 1.5 to 16.5, ). We did not find significant differences between age and the presence or not of drug-resistant M. tuberculosis (resistant years versus non resistant years, ). Nevertheless, age was a significant factor () in patients with MDR-TB ( years) compared to patients with drug-sensitive TB ( years). Drug-resistant isolates were identified in 27.5% of patients with both localized and disseminated disease (). Patients with MDR-TB were less likely to present systemic TB although this was not statistically significant, (localized 13% versus systemic 4.5%, ). The prevalence of drug-resistant strains showed variation between different Mexican states. In seven of the twelve states where drug resistance was documented, MDR-TB strains were identified. Follow-up treatment (6 months for pulmonary TB and 12 months for other clinical forms) was documented in 65 (71%) patients; 6 of these patients (9%) died, 2 (3%) stopped treatment, and 1 (1%) patient reported treatment failure. In this group, the presence of drug resistance or MDR did not contribute to the risk of death.

3.5. Spoligotyping

Spoligotyping produced a total of 55 different patterns for the 90 isolates studied. The succinct analyses obtained of the prevailing genotypes within this subpopulation as compared to the spoligotypes recorded in the SITVIT2 database are summarized in Tables 2, 3, and 4. Nineteen patterns corresponded to orphan strains that were unique among all the strains recorded in the SITVIT2 database (Table 2), as opposed to 36 patterns from 71 patients that corresponded to shared types (Table 3), that is, an identical pattern shared by 2 or more patients worldwide (within this study, or matching another strain in the SITVIT2 database). As shown in Table 3 for the 36 SITs recorded, a total of 32 SITs (containing 65 isolates) matched a preexisting shared type in the SITVIT2 database, whereas 4 SITs (containing 6 isolates) were newly created either within the present study or after a match with an orphan in the database. Irrespective of the database comparison, 12 patterns corresponded to clusters in the present study (Table 4, 12 clusters containing 47 isolates, 2–12 isolates per cluster), accounting for a clustering rate of 52.2% (47/90). However, no statistically significant association between spoligotype patterns, age, region, drug resistance, or HIV serology was observed.

Table 2: Description of orphan strains () and corresponding spoligotyping defined lineages.

ISO number
Year
Strain
Country of origin
Spoligotype description
Octal code
Clade**

MEX0120021C1209
2002
C1209
MEX

640033777777600
BOVIS

MEX0120021C121
2002
C121
MEX
517367636700771
Unk

MEX0120021C1262
2002
C1262
MEX
777767607540131
LAM4

MEX0120021C1365
2002
C1365
MEX
677767606740371
LAM8

MEX0120021C1677
2002
C1677
MEX
640033777777710
Unk

MEX0120021C1842
2002
C1842
MEX
731044606763771
Unk

MEX0120021C521
2002
C521
MEX
777777602003771
Unk

MEX0120021C672
2002
C672
MEX
410222000000031
Unk

MEX0120022C615
2002
C615
MEX
770346627740371
Unk

MEX0120022C890
2002
C890
MEX
777340004020771
H1

MEX0120023C111
2002
C111
MEX
711022103370371
Unk

MEX012002C1519
2002
C1519
MEX
676773777777710
Unk

MEX012002C164
2002
C164
MEX
777367626760771
T1

MEX012002C855
2002
C855
MEX
477777477403701
Unk

MEX0120031C1261
2003
C1261
MEX
700076777740271
X3

MEX0120031C1283
2003
C1283
MEX
777777607560531
LAM6

MEX0120031C54
2003
C54
MEX
674363777776000
BOVIS

MEX0120031C655
2003
C655
MEX
777367623540771
T1

MEX0120031C765
2003
C765
MEX
517367637760771
T1

**Clade designations are according to SITVIT2 database. Unk: unknown patterns within major clades described in SITVIT2.

Table 3: Description of M. tuberculosis complex shared types in pediatric patients ().

SIT
Spoligotype description
Octal number
Total (%) in this study
% in this study versus database
Clade
Clustered versus unique patterns

1

000000000003771
1 (1.1)
0.02
Beijing
Unique

3

000000007720771
1 (1.1)
1.09
H3
Unique

17

677737607760771
1 (1.1)
0.17
LAM2
Unique

20

677777607760771
1 (1.1)
0.14
LAM1
Unique

33

776177607760771
2 (2.2)
0.19
LAM3
Clustered

34

776377777760771
1 (1.1)
0.13
S
Unique

42

777777607760771
12 (13.3)
0.44
LAM9
Clustered

47

777777774020771
4 (4.0)
0.32
H1
Clustered

49

777777777720731
2 (2.2)
1.46
H3
Clustered

50

777777777720771
3 (3.3)
0.11
H3
Clustered

53

777777777760771
9 (10.0)
0.18
T1
Clustered

70

700076777760671
3 (3.3)
2.59
X3
Clustered

81

777377607760771
1 (1.1)
3.23
LAM9
Unique

91

700036777760771
1 (1.1)
0.42
X3
Unique

92

700076777760771
1 (1.1)
0.25
X3
Unique

95

777777607560731
1 (1.1)
2.56
LAM6
Unique

119

777776777760771
4 (4.4)
0.39
X1
Clustered

130

776177607760731
1 (1.1)
1.1
LAM3
Unique

239

777777777760031
1 (1.1)
1.92
T1
Unique

373

777777767760771
1 (1.1)
1.75
T1
Unique

398

777777607760631
1 (1.1)
8.33
LAM4
Unique

450

777776770000000
1 (1.1)
1.19

Unique

804

477777777760771
1 (1.1)
5
T1
Unique

1166

777377777760771
1 (1.1)
2.63
T1
Unique

1193

777777677773771
2 (2.2)
50
Manu1
Clustered

1243

777377777720771
2 (2.2)
14.29
H3
Clustered

1277

777777207760771
1 (1.1)
4.55
LAM9
Unique

1626

777777776760771
1 (1.1)
12.5
T1
Unique

1711

677337607760771
1 (1.1)
25
LAM2
Unique

2188

777377777760731
1 (1.1)
25
T2
Unique

2337

777377774020771
1 (1.1)
16.67
H1
Unique

2496

777367777760771
1 (1.1)
7.14
T1
Unique

3079*

476363636756200
2 (2.2)
100
BOVIS
Clustered

3080*

777377777763771
2 (2.2)
66.67
Manu2
Clustered

3089*

337777777760771
1 (1.1)
50
T1
Unique

3090*

677773777777600
1 (1.1)
50
BOVIS
Unique

Shared types (SIT) followed by an asterisk indicate “newly created” shared types () due to 2 or more strains belonging to an identical new pattern within this study (SIT3079, SIT3080), or a unique strain from this study matching with another orphan in the database (SIT3089 matched an orphan from Peru; SIT3090 matched an orphan from USA). Note that SIT3080 was created by 2 strains belonging to an identical pattern within this study that also matched an orphan from Saudi Arabia.Clade designations are according to SITVIT2 database. Unk: unknown patterns within major clades described in SITVIT2.Clustered versus unique patterns; clustered strains correspond to a similar spoligotype pattern shared by 2 or more strains “within this study”; as opposed to unique strains harboring a spoligotype pattern that does not match with another strain from this study. Unique strains matching a preexisting pattern in the SITVIT2 database are classified as SITs, whereas in case of no match, they are designated as “orphan” (see Table 1).

Table 4: Description of shared types representing clustered strains and their worldwide distribution.

SIT (Clade), octal number spoligotype description
Total (%)in this study
% in this study versus database
Distribution in regions with ≥5% of a given SIT
Distribution in countries with ≥5% of a given SIT

2 (2.2)
0.19
AFRI-S 31.46, AMER-S 21.68, AMER-N 15.78, EURO-S 13.36, EURO-W 5.32
ZAF 31.46, USA 15.39, BRA 8.81, ESP 8.62, ARG 5.52, PER 5.42

12 (13.3)
0.44
AMER-S 28.25, AMER-N 15.85, EURO-S 12.11, AFRI-N 10.44, EURO-W 6.65
USA 14.44, BRA 9.78, MAR 8.58, COL 7.45, ITA 6.36

4 (4.4)
0.32
AMER-N 21.20, EURO-W 20.06, EURO-S 13.73, AMER-S 10.40, EURO-E 8.37
USA 19.33, AUT 10.24, ITA 7.23, BRA 5.85

2 (2.2)
1.46
AMER-N 21.90, EURO-N 21.17, EURO-W 13.14, AMER-S 12.41, EURO-S 10.95, AFRI-M 5.84
USA 18.98, FIN 16.79, PER 5.84, AUT 5.11, PRT 5.11, FXX 5.11

3 (3.3)
0.11
AMER-N 21.94, AMER-S 15.88, EURO-W 15.34, EURO-S 12.62, EURO-E 6.45, AFRI-N 5.16
USA 21.33, AUT 7.38, ESP 6.59, ITA 5.27

9 (10.0)
0.18
AMER-N 19.55, AMER-S 14.24, EURO-W 12.62, EURO-S 9.87, ASIA-W 8.55, AFRI-S 6.40
USA 17.05, ZAF 6.26, ITA 5.05

3 (3.3)
2.59
AMER-N 68.10, CARI 17.24, AMER-S 5.17
USA 60.34, MEX 7.76, BHS 6.90, HTI 6.90, GUF 5.17

4 (4.4)
0.39
AMER-N 70.10, AFRI-S 14.02
USA 61.76, ZAF 14.02, MEX 7.65

2 (2.2)
50.00
ASIA-S 50.00, AMER-N 50.00
MEX 50.00, IND 50.00

2 (2.2)
14.29
EURO-S 35.71, ASIA-W 21.43, AMER-N 14.29, ASIA-SE 7.14, AMER-S 7.14, EURO-W 7.14, ASIA-E 7.14
ESP 28.57, MEX 14.29, TUR 14.29, ITA 7.14, NLD 7.14, SAU 7.14, TWN 7.14, COL 7.14, MYS 7.14

2 (2.2)
100.00
AMER-N 100.00
MEX 100.00

2 (2.2)
66.67
AMER-N 66.67, ASIA-W 33.33
MEX 66.67, SAU 33.33

Worldwide distribution is reported for regions with ≥5% of a given SITs as compared to their total number in the SITVIT2 database. The definition of macrogeographical regions and subregions (http://unstats.un.org/unsd/methods/m49/m49regin.htm) is according to the United Nations; Regions: AFRI (Africa), AMER (Americas), ASIA (Asia), EURO (Europe), and OCE (Oceania), subdivided in: E (Eastern), M (Middle), C (Central), N (Northern), S (Southern), SE (South-Eastern), and W (Western). Furthermore, CARIB (Caribbean) belongs to Americas, while Oceania is subdivided in 4 subregions: AUST (Australasia), MEL (Melanesia), MIC (Micronesia), and POLY (Polynesia). Note that in our classification scheme, Russia has been attributed a new sub-region by itself (Northern Asia) instead of including it among rest of the Eastern Europe. It reflects its geographical localization as well as due to the similarity of specific TB genotypes circulating in Russia (a majority of Beijing genotypes) with those prevalent in Central, Eastern, and South-Eastern Asia.The 3 letter country codes are according to http://en.wikipedia.org/wiki/ISO_3166-1_alpha-3; countrywide distribution is only shown for SITs with ≥5% of a given SITs as compared to their total number in the SITVIT2 database.

3.6. Distinction of “Ancient” versus “Modern” Lineages

Spoligotyping results and clade definitions were linked to the distribution of clinical isolates within the principal genetic group PGG1 versus PGG2/3 (the latter being putatively characterized by the lack of spacers 33–36), it was evident that only 10/90 (11.1%) strains belonged to PGG1 (BOVIS, Manu and Beijing/W lineages), the rest being classified as evolutionary recent PGG2/3 strains. Indeed, LAM was the most predominant clade (25/90 or 27.8% of the isolates, including the major cluster composed of SIT42/LAM9 sublineage, = 12/90 or 13.3% of the isolates), followed by the Haarlem and T families.

4. Discussion

Our study aimed to estimate the genotype diversity and the global prevalence of drug-resistant TB in pediatric patient population. This is the first report of resistance to primary antimicrobial drugs in an exclusively pediatric Mexican population. The extent of childhood TB is unknown and is estimated that has a limited impact on the dynamics of the TB epidemic; however, pediatric TB indicates recent transmission of a disease that has not been previously treated or diagnosed. Therefore, childhood TB can be a reflection on the inefficiency of programs for the control of the disease. In this study, a high frequency of primary antimicrobial resistance was observed (1 out of every 4 isolates). Resistance to INH alone or in combination with other anti-TB drugs was most prevalent, with a global resistance of 23%, a finding which is higher than the 18% prevalence previously reported in population based studies in Mexico [25, 26]. In this respect, our observation is similar to prevalence of 29% to 56% reported in studies in Chiapas and Veracruz (Mexico) for high-risk populations [27–29]. When the prevalence of INH resistance is ≥4%, the use of a fourth drug (myambutol or streptomycin) is recommended to ensure therapeutic success; a standard that was adopted in Mexico in 2000 as part of the TAES strategy—Norma oficial Mexicana NOM-006-SSA2-1993 para la prevención y control de la tuberculosis para la atención primaria de la salud (http://www.salud.gob.mx/unidades/cdi/nom/m006ssa23.html).

Importantly, the magnitude of RMP-resistance in our study (11%) is higher than previously determined level in Mexico (3%). Resistance to RMP was always accompanied by INH-resistance leading to 10 (11.1%) cases of MDR-TB. Although because of the small number of isolates studied, this high rate of MDR-TB observed (11%) might not reflect the exact rates of drug resistance in this age group, this finding is relevant because RMP and INH are among the most important 1st-line drugs involved both in the TAES strategy (see above) as well as in the WHO-recommended Directly Observed Therapy—Short Course (DOTS) strategy (http://www.who.int/tb/dots/en/index.html). MDR-TB is linked to a greater probability of therapeutic failure and a lower survival rate, with a calculated relative risk for mortality of 2.5 (95% CI 1.02 to 6.16, ) [29]. The occurrence of MDR-TB has logistical implications in the necessity of tests for drug resistance and the use of second-line drugs.

Because of the limited number of patients in our study, it is difficult to compare the mortality of MDR-TB, which is lower than reported in other national studies [30, 31]. Distribution clusters of any kind were not observed for resistant M. tuberculosis isolates, however, 20% of patients were retreated due to therapeutic failure in our study, a finding which is in agreement with previous reports showing that previous TB treatment is a predictor for resistance [32, 33]. These observations stress the importance of performing periodical population studies to determine the prevalence and incidence of drug-resistant M. tuberculosis and to establish the most effective strategies for the control TB.

Spoligotyping of the clinical isolates showed the great diversity of circulating genotypes, with 55 different patterns which is common in countries with endemic. The spoligotypes were not associated with the geographical location of the patient or the presence of resistance to one or more antimicrobial drugs. The worldwide distribution of predominant shared types and lineages from our study underlined that many of the genotypes encountered in the pediatric population were the ones most frequently represented in Mexico and neighboring USA, the rest being essentially shared with other Latin American and Caribbean neighbors (Table 4). Thus with the exception of ubiquitous spoligotypes (such as the T clade found throughout the world), the patients in our study mainly harboured M. tuberculosis spoligotypes prevailing in the Americas, and mostly belonging to evolutionary recent PGG2/3 lineages. However, the “T” genotype does not represent a clade in a strict evolutionary sense since it was defined by default to include strains that may not be classified in one of the established genotypic lineages with well-established phylogeographical specificity such as the Haarlem, LAM, CAS, and EAI lineages [20]. Interestingly, out of the 12 sublineages reported worldwide for the LAM clade [20], a total of 6 sublineages were present in our 2-year recruitment. Last but not least, the near absence of Beijing genotype () in this study is noteworthy and none were associated with 2 HIV patients (one of the HIV patients harbored a T1 lineage strain).

A remarkable feature of our study is the presence of a few ancestral Manu lineage strains ( = 4/90 or 4.4%) as clustered strains within our study sample. The “Manu” lineage was initially described as a new family from India in 2004 [34], and later similar strains in small numbers were reported in a study from Madagascar [35]. Soon afterwards, it was tentatively subdivided into Manu-1 (deletion of spacer 34), Manu-2 (deletion of spacers 33-34), and Manu-3 (deletion of spacers 34–36) sublineages, and it was suggested that it could represent an ancestral clone of principal genetic group 1 strains [20]. Manu lineage strains were recently reported from Saudi Arabia [36], Tunisia [37], and more recently in a study from Egypt where it represented as high as 27% of all isolates [38].

Contrary to a recent study on M. tuberculosis spoligotyping in Acapulco city that showed a single cluster of 70 (26%) patients harboring a single spoligotype (SIT19) corresponding to the EAI2-Manila lineage [39], we did not find a single SIT19 isolate in the pediatric patient population of our study. Neither did we find any shared type belonging to the EAI family which represents ancestral PGG1 strains within the M. tuberculosis complex [20]. We therefore conclude that the population structure of pediatric TB in our setting is entirely different from the one prevailing in adult TB patient population of Guerrero with evidence of an ongoing transmission with ancestral EAI2-Manila lineage [39]. In conclusion, our study shows that TB among pediatric patients in Mexico is essentially caused by evolutionary recent genotypes characteristic of the Americas. However, the presence of the ancestral Manu lineage strains, supposed to be a missing link of the split between ancestral and modern tubercle bacilli during M. tuberculosis evolution [38], should be further investigated within larger datasets to know its arrival in Mexico. Lastly, the clustering rate observed in our study (52.2%) is certainly higher than expected since further differentiation of spoligotyping defined clusters was not systematically performed using secondary markers such as MIRU-VNTRs [40], which will be necessary in future studies to draw epidemiological conclusions.

Acknowledgments

The authors are grateful to Dr. Armando Rentería Cárdenas (San Luis Potosí, Mexico) and Dr. Gregorio Cetina (Mérida, Yucatán, Mexico) for providing some of the clinical information. N. Rastogi acknowledges the help of Véronique Hill (Institut Pasteur de Guadeloupe) for helping with SITVIT2 database project which benefited from research Grants from the Regional Council of Guadeloupe (Grant CR/08-1612), and by United States National Institutes of Health (NIH Grant R01LM009731).

References

C. Dye, S. Scheele, P. Dolin, V. Pathania, and M. C. Raviglione, “Consensus statement. Global burden of tuberculosis: estimated incidence, prevalence, and mortality by country. WHO Global Surveillance and Monitoring Project,” The Journal of the American Medical Association, vol. 282, no. 7, pp. 677–686, 1999.

P. Glaziou, K. Floyd, and M. Raviglione, “Global Burden and Epidemiology of Tuberculosis,” Clinics in Chest Medicine, vol. 30, no. 4, pp. 621–636, 2009.

The World Health Organization, “Anti-tuberculosis drug resistance in the world—Fourth Global report,” The WHO/IUATLD Global Project on Anti-tuberculosis Drug Resistance Surveillance 2002–2007, WHO, Geneva, Switzerland, 2008, http://www.who.int/tb/publications/2008/drs_report4_26feb08.pdf.

P. R. Donald and P. D. Van Helden, “The global burden of tuberculosis—combating drug resistance in difficult times,” New England Journal of Medicine, vol. 360, no. 23, pp. 2393–2395, 2009.

T. Walls and D. Shingadia, “Global epidemiology of paediatric tuberculosis,” Journal of Infection, vol. 48, no. 1, pp. 13–22, 2004.

A. Pablos-Méndez, M. C. Raviglione, A. Laszlo et al., “Global surveillance for antituberculosis-drug resistance, 1994-1997,” New England Journal of Medicine, vol. 338, no. 23, pp. 1641–1649, 1998.

M. C. Raviglione, “The global plan to stop TB, 2006-2015,” International Journal of Tuberculosis and Lung Disease, vol. 10, no. 3, pp. 238–239, 2006.

Centers for Disease Control and Prevention, “The global plan to stop TB, 2006-2015: summary,” Morbidity and Mortality Weekly Report, vol. 81, pp. 86–88, 2006.

E. Nathanson, R. Gupta, P. Huamani et al., “Adverse events in the treatment of multidrug-resistant tuberculosis: results from the DOTS-Plus initiative,” International Journal of Tuberculosis and Lung Disease, vol. 8, no. 11, pp. 1382–1384, 2004.

G. Maartens and R. J. Wilkinson, “Tuberculosis,” Lancet, vol. 370, no. 9604, pp. 2030–2043, 2007.

P. Lobue, “Extensively drug-resistant tuberculosis,” Current Opinion in Infectious Diseases, vol. 22, no. 2, pp. 167–173, 2009.

Centers for Disease Control and Prevention, “Extensively drug-resistant tuberculosis—United States, 1993–2006,” Morbidity and Mortality Weekly Report, vol. 56, no. 11, pp. 250–253, 2007.

A. R. De Moura Valim, L. G. Possuelo, P. I. Cafrune et al., “Evaluation and genotyping of multidrug-resistant cases of tuberculosis in southern Brazil,” Microbial Drug Resistance, vol. 12, no. 3, pp. 186–191, 2006.

F. A. Post, P. A. Willcox, B. Mathema et al., “Genetic polymorphism in Mycobacterium tuberculosis isolates from patients with chronic multidrug-resistant tuberculosis,” Journal of Infectious Diseases, vol. 190, no. 1, pp. 99–106, 2004.

A. N. J. Malik and P. Godfrey-Faussett, “Effects of genetic variability of Mycobacterium tuberculosis strains on the presentation of disease,” Lancet Infectious Diseases, vol. 5, no. 3, pp. 174–183, 2005.

C. L. Daley, “Molecular epidemiology: a tool for understanding control of tuberculosis transmission,” Clinics in Chest Medicine, vol. 26, no. 2, pp. 217–231, 2005.

M. Mondragón-Barreto, C. A. Vázquez-Chacón, C. Barrón-Rivero et al., “Comparison among three methods for mycobacteria identification,” Salud Publica de Mexico, vol. 42, no. 6, pp. 484–489, 2000.

S. Kulkarni, C. Sola, I. Filliol, N. Rastogi, and G. Kadival, “Spoligotyping of Mycobacterium tuberculosis isolates from patients with pulmonary tuberculosis in Mumbai, India,” Research in Microbiology, vol. 156, no. 4, pp. 588–596, 2005.

J. Kamerbeek, L. Schouls, A. Kolk et al., “Simultaneous detection and strain differentiation of Mycobacterium tuberculosis for diagnosis and epidemiology,” Journal of Clinical Microbiology, vol. 35, no. 4, pp. 907–914, 1997.

K. Brudey, J. R. Driscoll, L. Rigouts et al., “Mycobacterium tuberculosis complex genetic diversity: mining the fourth international spoligotyping database (SpolDB4) for classification, population genetics and epidemiology,” BMC Microbiology, vol. 6, article no. 23, 2006.

S. Sreevatsan, X. Pan, K. E. Stockbauer et al., “Restricted structural gene polymorphism in the Mycobacterium tuberculosis complex indicates evolutionarily recent global dissemination,” Proceedings of the National Academy of Sciences of the United States of America, vol. 94, no. 18, pp. 9869–9874, 1997.

R. Brosch, S. V. Gordon, M. Marmiesse et al., “A new evolutionary scenario for the Mycobacterium tuberculosis complex,” Proceedings of the National Academy of Sciences of the United States of America, vol. 99, no. 6, pp. 3684–3689, 2002.

H. Soini, X. Pan, A. Amin, E. A. Graviss, A. Siddiqui, and J. M. Musser, “Characterization of Mycobacterium tuberculosis isolates from patients in Houston, Texas, by spoligotyping,” Journal of Clinical Microbiology, vol. 38, no. 2, pp. 669–676, 2000.

N. Rastogi and C. Sola, “Molecular evolution of the Mycobacterium tuberculosis complex,” in Tuberculosis, J. C. Palomino, S. Leao, and V. Ritacco, Eds., pp. 53–91, Amedeo Online Textbooks, 2007, http://www.tuberculosistextbook.com/index.htm.

Centers for Disease Control and Prevention, “Population-based survey for drug resistance of tuberculosis—Mexico, 1997,” Morbidity and Mortality Weekly Report, vol. 47, pp. 371–375, 1998.

J. Sifuentes-Osornio, L. A. Ponce-De-Leon, F. E. Camacho-Mezquita et al., “Drug resistance of M. tuberculosis in Mexican patients. I. Clinical manifestations and risk factors,” Revista de Investigacion Clinica, vol. 47, no. 4, pp. 273–281, 1995.

C. R. Peter, E. Schultz, K. Moser et al., “Drug-resistant pulmonary tuberculosis in the Baja California-San Diego County border population,” Western Journal of Medicine, vol. 169, no. 4, pp. 208–213, 1998.

M. L. García-García, Ma. E. Jiménez-Corona, A. Ponce-De-León et al., “Mycobacterium tuberculosis drug resistance in a suburban community in Southern Mexico,” International Journal of Tuberculosis and Lung Disease, vol. 4, no. 12, supplement 2, pp. S168–S170, 2000.

M. L. García-García, A. Ponce-De-León, M. E. Jiménez-Corona et al., “Clinical consequences and transmissibility of drug-resistant tuberculosis in southern Mexico,” Archives of Internal Medicine, vol. 160, no. 5, pp. 630–636, 2000.

Z. H. Yang, A. Rendon, A. Flores et al., “A clinic-based molecular epidemiologic study of tuberculosis in Monterrey, Mexico,” International Journal of Tuberculosis and Lung Disease, vol. 5, no. 4, pp. 313–320, 2001.

M. L. García-García, J. L. Valdespino-Gómez, M. Palacios-Martínez, M. E. Mayar-Maya, C. García-Sancho, and J. Sepúlveda-Amor, “Tuberculosis and AIDS in Mexico,” Salud Pública de México, vol. 37, no. 6, pp. 539–548, 1995.

M. A. Espinal, A. Laszlo, L. Simonsen et al., “Global trends in resistance to antituberculosis drugs,” New England Journal of Medicine, vol. 344, no. 17, pp. 1294–1303, 2001.

M. C. Becerra, J. Bayona, J. Freeman, P. E. Farmer, and J. Y. Kim, “Redefining MDR-TB transmission ‘hot spots’,” International Journal of Tuberculosis and Lung Disease, vol. 4, no. 5, pp. 387–394, 2000.

U. B. Singh, N. Suresh, N. V. Bhanu et al., “Predominant tuberculosis spoligotypes, Delhi, India,” Emerging Infectious Diseases, vol. 10, no. 6, pp. 1138–1142, 2004.

S. Ferdinand, C. Sola, S. Chanteau et al., “A study of spoligotyping-defined Mycobacterium tuberculosis clades in relation to the origin of peopling and the demographic history in Madagascar,” Infection, Genetics and Evolution, vol. 5, no. 4, pp. 340–348, 2005.

S. A. M. Al-Hajoj, T. Zozio, F. Al-Rabiah et al., “First insight into the population structure of Mycobacterium tuberculosis in Saudi Arabia,” Journal of Clinical Microbiology, vol. 45, no. 8, pp. 2467–2473, 2007.

A. Namouchi, A. Karboul, B. Mhenni et al., “Genetic profiling of Mycobacterium tuberculosis in Tunisia: predominance and evidence for the establishment of a few genotypes,” Journal of Medical Microbiology, vol. 57, no. 7, pp. 864–872, 2008.

Z. H. Helal, M. S. E. D. Ashour, S. A. Eissa et al., “Unexpectedly high proportion of ancestral manu genotype Mycobacterium tuberculosis strains cultured from tuberculosis patients in Egypt,” Journal of Clinical Microbiology, vol. 47, no. 9, pp. 2794–2801, 2009.

E. Nava-Aguilera, Y. López-Vidal, E. Harris et al., “Clustering of Mycobacterium tuberculosis cases in Acapulco: spoligotyping and risk factors,” Clinical and Developmental Immunology, vol. 2011, Article ID 408375, 12 pages, 2011.

P. Supply, C. Allix, S. Lesjean et al., “Proposal for standardization of optimized mycobacterial interspersed repetitive unit-variable-number tandem repeat typing of Mycobacterium tuberculosis,” Journal of Clinical Microbiology, vol. 44, no. 12, pp. 4498–4510, 2006.