The control of tuberculosis in densely populated cities is complicated by close human-to-human contacts and potential transmission of pathogens from multiple sources. We conducted a molecular epidemiologic analysis of 356 Mycobacterium tuberculosis (MTB) isolates from patients presenting pulmonary tuberculosis in metropolitan Taipei. Classical antibiogram studies and genetic characterization, using mycobacterial interspersed repetitive-unit-variable-number tandem-repeat (MIRU-VNTR) typing and spoligotyping, were applied after culture.
Methods
A total of 356 isolates were genotyped by standard spoligotyping and the strains were compared with in the international spoligotyping database (SpolDB4). All isolates were also categorized using the 15 loci MIRU-VNTR typing method and combin with NTF locus and RD deletion analyses.
Results
Of 356 isolates spoligotyped, 290 (81.4%) displayed known spoligotypes and 66 were not identified in the database. Major spoligotypes found were Beijing lineages (52.5%), followed by Haarlem lineages (13.5%) and EAI plus EAI-like lineages (11%). When MIRU-VNTR was employed, 140 patterns were identified, including 36 clusters by 252 isolates and 104 unique patterns, and the largest cluster comprised 95 isolates from the Beijing family. The combination of spoligotyping and MIRU-VNTR revealed that 236 (67%) of the 356 isolates were clustered in 43 genotypes. Strains of the Beijing family was more likely to be of modern strain and a higher percentage of multiple drug resistance than other families combined (P = 0.08). Patients infected with Beijing strains were younger than those with other strains (mean 58.7 vs. 64.2, p = 0.02). Moreover, 85.3% of infected persons younger than 25 years had Beijing modern strain, suggesting a possible recent spread in the young population by this family of TB strain in Taipei.
Conclusion
Our data on MTB genotype in Taipei suggest that MTB infection has not been optimally controlled. Control efforts should be reinforced in view of the high prevalence of the Beijing strain in young population and association with drug resistance.
Background
Tuberculosis (TB) remains a worldwide healthcare concern and has been characterized as an epidemic by World Health Organization (WHO). It is estimated one third of the world's population has been infected with Mycobacterium tuberculosis (MTB) and that 3 million people will die of the disease per year between now and 2010. The distribution of TB in different geographic regions is characterized by the prevalence of different MTB strains with varied virulence and drug resistance. Both environmental and host factors are responsible for the transmission and prevalence of different MTB strains. Although both the incidence and mortality rates of TB in Taiwan have shown a steady decline since 1950, TB remains a leading notifiable infectious disease on the island. In 2001, 14,486 cases were reported, with a notification rate of 64.9 per 100,000 people.
At the molecular level, the global TB epidemic consists of multiple genotype-specific subepidemics. Different MTB genotypes can be identified by variation in certain well-characterized repetitive sequences, such as the IS6110 transposable element and the direct repeat region [2]. The Beijing genotype family is well recognized as having a distinct genetic signature, and it is genetically highly conserved [3] even though sequence polymorphisms have identified four monophyletic subgroups [4]. It is dispersed worldwide yet predominates in certain geographic areas, particularly in parts of Asia [5, 6] and Russia [7]. Its prevalence in the patient populations of recent studies in Vietnam and Russia suggests the recent spread to those areas [8]. It has been proposed that "Beijing" should be regarded as an emerging genotype family [9].
The association between drug resistance and the Beijing genotype is well documented in recent medical literature [3, 4, 6, 10–13]. The geographic variability observed in this association [4], along with the frequent clustering of resistant genotypes and their successful spread within the Russian prison system [7, 13] suggests recent colonial expansion. This is further supported by the evidence that some strains of the Beijing genotype family retain fitness despite the acquisition of drug resistance [13]. The Haarlem family genotype has a similar relationship with drug resistance and rapid clonal expansion [14]. The association of these genotype families with drug-resistant outbreaks clearly demonstrates their epidemic potential [14, 15]. From a TB-control point of view, it is relevant to understand whether specific genotype families are overrepresented among drug-resistant cases and, in particular, if these resistant strains are successfully transmitted within the community. Taipei is a metropolitan city in northern Taiwan with a population of 2.3 million inhabiting a basin of 272 square km. The population of Taipei includes Han Chinese whose ancestors migrated to this island in the 16th century, the veterans who retreated to the island in the late 1940s during the Chinese civil war, and the Taiwanese Aborigines, who have resided on this island since before the 16th century [16]. The prevalence of TB in large urban areas is complicated by the close human-to-human contacts and potential multiple sources of MTB strains from different ethnic and migratory populations. The goals of this study were therefore to characterize the prevalence of genotypes, cluster pattern, and drug resistance of MTB isolates in Taipei to provide information for potential transmission and formulation of effective infection-control policy.
Methods
Mycobacterial strains and genomic DNA
A total of 356 samples were randomly collected between 2002 and 2004 from 356 patients at the Tri-Service General Hospital, a large medical center that handles a substantial number of TB patients referred from hospitals throughout Taipei. All of the patients were sputum microscopy positive and culture positive. Mycobacterial genomic DNA was extracted from cultured cells as described previously [17, 18]. Resuspending mycobacterial colonies in 100 to 200 μl of distilled H2O and incubating them at 85°C for 30 min obtained genomic DNA. After centrifugation of the suspension, the supernatant containing the DNA was removed and stored at -20°C until further use. The study protocol has been approved by the institutional review board of the National Health Research Institutes, Taiwan.
Spoligotyping and spoligotype analysis
Spoligotyping was carried out according to the manufacturer's instructions (Isogen Bioscience B.V., Maarsen, Netherlands). The resulting spoligotypes were documented under a binary code representing either a positive or negative hybridization result (n and o, respectively) and analyzed using the Excel program for grouping and ordering of the patterns. Spoligotypes common to more than one strain were designated as shared types (ST) and assigned a shared international type number (SIT) according to the updated version of the international spoligotype database SpolDB4 [19].
PCR and MIRU analysis
PCRs were carried out using the PCR reagent system (Gibco-BRL). Sequences of primers used for amplification of 12 MIRU loci and 3 ETR loci (A, B, C) were selected according to descriptions in other studies [20]. Five microliters from fivefold-diluted DNA solutions were added to a final volume of 50 μl containing 0.2 μl of DNA polymerase (1 U); 0.2 mM each of dATP, dCTP, dGTP, and dTTP; 5 μl of PCR buffer; 0.4 μM (2 μM for locus 7) of primers; and 1 to 3.5 mM of MgCl2. The primers and MgCl2 concentrations used were as described by Mazars et al. [21]. The PCR fragments were analyzed by agarose gel electrophoresis with 1.5% agarose. The sizes of the amplicons were estimated by comparison with 50- and 100-bp ladders. The MIRU copy number per locus was calculated by using the conventions described by Supply et al. [22].
TbD1 Analysis
According to Brosch et al. [23], TbD1 is specifically present in the ancestral lineage of MTB. The presence of TbD1 was analyzed by PCR. Briefly, two PCR assays were performed per isolate by using either primers complementary to the sequences flanking the deleted region or primers complementary to the internal sequences. For the isolates that did (TbD1+) or did not (TbD1-) contain the TbD1 region, an amplicon was obtained only with internal primers or only with flanking primers, respectively.
NTFlocus and RD deletion analysis
A multiplex PCR approach was used to determine possible IS6110 insertion(s) in the NTF region of M. tuberculosis strains. The method, including primers within the NTF region as well as the IS6110 sequence and PCR parameters, was adapted from a previously described paper by Plikaytis et al. [24].
A primer set was used to check for the presence or absence of RD105, RD181, RD150, RD142, and RD207. The PCR mixture consisted of 0.2 μg DNA template, 13.9 μl Q buffer, 5 μl 5× buffer, 4 μl 10 mM deoxynucleoside triphosphates, 1 μl of each primer (10 pmol/μl), 1μl DMSO, and 0.6 μl Herculase II Fusion DNA polymerase (STRATAGENE, USA). Sterile water was used to dilute the mixture up to 25 μl. A detailed explanation of this methodology has been described [25–27].
Drug resistance testing
The proportional method for drug susceptibility testing (DST) of MTB was performed as described previously [28]. Briefly, for each drug a 1:10 dilution of standardized suspensions was inoculated onto the control and drug-containing media. The extent of growth in the absence or presence of the drug was compared and expressed as a percentage. If growth at the critical concentration of a drug was >1%, the isolate was considered to be clinically resistant. 7H10 agar with 0.2 or 1 mg/l isoniazid (INH), 1 or 5 mg/l rifampicin, 5 or 10 mg/l ethambutol, and 5 or 10 mg/l streptomycin was used.
Statistical analyses
Frequencies of multiple-drug resistance (MDR) among different genotype families based on spoligotyping were compared with a chi-square test, or a Fisher's exact test when any of the cells had expected counts ≤ 5. The extent of association was expressed as an odds ratio (OR) and 95% confidence interval (95% C.I.). All statistical tests were two-sided; and statistical significance was set at a p-value <0.05.
Patients in this study can be classified into two groups, characterized by clustered and non-clustered MTB isolates. A possible cluster is defined as two or more patients' strains with identical genetic patterns defined by the MIRU-VNTR typing; patients' strains with unmatched genetic profiles were considered non-clustered. Previous literatures have suggested that clusters may be assumed to have arisen from recent transmission; and the clustering rate was used to determine the amount of recent transmission in this population [29, 30]. The patients' strains with the same genetic pattern may represent an epidemiologically linked cluster. Therefore, the minimum estimate of the proportion of M. tuberculosis cases related to recent transmission can be calculated as (number of clustered patients minus the number of clusters)/total number of patients (Additional file 1).
Results
Analysis of the spoligotyping patterns
During the study period, 356 patients were diagnosed with culture-confirmed TB. Molecular analysis showed that all TB cases were caused by M. tuberculosis, except one by Mycobacterium bovis. The median age of these patient was 61.2 years, and 70%(251/356) were male. Spoligotyping and drug susceptibility testing (DST) were performed on all of the specimens. Of the 356 isolates analyzed, spoligotypes from 290 isolates (81.4%) were classified according to SpolDB4 into one of the 47 shared international types (SITs) (Figure 1). Of the remaining 66 isolates, 40 patterns were not identified in the database, 22 were of the East African-India (EAI), 4 were of the Beijing, and the other 40 (11.2%) were orphans (Figure 2). Of the 47 defined spoligotypes, the most frequent strain found was the Beijing spoligotype ST1 (49.43%), followed by ST50 (5.9%) of the Haarlem strains, and ST19 (3.65%) of the EAI_2 Manilla strains (Figure 1). The Beijing family was the most prevalent genotype, identified in 187/356 (52.5%) isolates, followed by the Haarlem family, identified in 48/356 (13.5%) isolates (Table 1). In these novel spoligotypes, 22 strains were found to be TBD1 positive and are further characterized by the absence of DR spacers 29 to 32 and 34 and the presence of spacer 33. Based on this result, these new spoligotypes belong to the East African-India (EAI) family. There are two novel spoligotypes with the RD 105 deletion, indicating their membership in the Beijing family (Figure 2). Of all the isolates studied, the most prevalent subfamilies after the Beijing family (52.5%) were H3 (13.2%) and EAI-like (6.2%) (Figure 3). Among all of the isolates, 63 (18%) displayed unique spoligotypes, and 293 (82%) displayed one of 24 spoligotypes.
Figure 1
Spoligotypes of 290 isolates with a shared international type (SIT) number in SpolDB4. a Shared international type (SIT), international spoligotype database SpolDB4 http://www.pasteur-guadeloupe.fr:8081/SITVITDemo/. b Label representing spoligotype families as assigned in the international spoligotype database SpolDB4. c Number of isolates in this study. d Prevalence, representing the number of isolates with a common SIT relative to the total number of isolates from the same database (356) classified by SIT from Tri-Service General Hospital (expressed as a percentile).
Figure 2
Spoligotypes of 66 orphan strains and clusters of spoligotypes not identified in SpolDB4 by a SIT numbera. a Shared international type (SIT), international spoligotype database SpolDB4 http://www.pasteur-guadeloupe.fr:8081/SITVITDemo/. b TBD1 positive, c RD105 deletion.
Figure 3
Prevalence of the major spoligotype families and subfamilies of 356M. tuberculosisisolates.
Table 1
Mycobacterial genotype and drug resistance in patients with culture-confirmed tuberculosis
No. of isolates with DST result
Genotype family
No. of isolates (%)
MDR(%)
Any one drug(%)
All sensitivity(%)
Beijinga
187(52.5)
8(4.2)
36(19.4)
143(76.4)
Haarlem
48(13.5)
0
9(18.8)
39(81.2)
EAIb
40(11.2)
0
3(7.5)
37(92.5)
T
25(7.1)
0
8(32.0)
17(68.0)
'Others'c(LAM, U, MANU, Bovis1)
16(4.5)
0
1(6.3)
15(93.7)
Unclassifiedd
40(11.2)
2(5)
8(20)
30(75)
Total
356
10(2.8)
65(18.2)
281(79)
a Including Beijing-like strains;
b Including EA-like strains;
c'Others', all genotype families with a frequency of less than 10 cases;
d unclassified, no internationally recognized genotype family assigned, based on the SpolDB4 spoligotype database.
The age distribution of patients in different genotype, which include Beijing, Haarlem, EAT, T and others genotype of average age in year respectively is 58.7, 62.1, 66.9, 61.4, 65.3. Statistical analyses by T-test demonstrated that patient infects with the Beijing family (mean = 58.7) were statistically younger than those infected with other geno-families (mean = 64.2) with a p value of 0.02. Moreover, 68.5% of patients 30 year-old or younger were infected with Beijing strains as compared to 50% of those older than 75 with Beijing strains, giving an odds ratio of 2.18 with 95% C.I. = 1.11–4.28 and p = 0.02 (Table 2). The association between young age and Beijing strains was even stronger in younger ages: 85.3% of people 25 year-old or less had Beijing strain, and when compared with those more than 75 year-old produced an odds ratio of 5.8 with p = 0.0002 (Table 2). These results indicated that young population in Taipei were more likely to be infected with Beijing strains than with other strains and as compared to older age groups, and thus suggested a possible recent spread of the Beijing genotype in the young generation in this area.
Table 2
Beijing family genotype M. tuberculosis isolates from 356 tuberculosis cases
No.(%) isolates No.
(%) of Beijing isolates
Odds Ratio
95% C.I.
p-value
Age group (yr)
356
187 (52.53)
≦ 25
34 (9.55)
29 (85.29)
5.80
2.11–15.98
0.0002
≦ 30
54 (15.17)
37 (68.52)
2.18
1.11–4.28
0.02
31–60
95 (26.69)
50 (52.63)
1.11
0.65–1.90
0.7
61–75
85 (23.88)
39 (45.88)
0.85
0.49–1.48
0.56
≧ 76
122 (34.27)
61 (50.00)
1
reference group
Typing of strains and clustering analysis by MIRU-VNTR
All 356 isolates were also categorized using the MIRU-VNTR typing method, which detected 140 different patterns, comprising 36 clusters formed by 252 isolates and 104 unique patterns formed by 104 isolates (Table 3). The largest cluster comprised 95 from the Beijing family, with an MIRU-VNTR profile of 223325173533. Of 187 Beijing strains as defined by spoligotyping, MIRU-VNTR typing further divided them into 47 different patterns, comprising 15 clusters formed by 155 isolates and 32 unique patterns formed by 32 isolates. Of 48 Harrlem strains, 24 different MIRU-VNTR patterns were found, which consisted of 5 clusters and 19 unique MIRU-VNTR genotypes. The largest cluster comprised 10 in the Haarlem family, with an MIRU profile of 222225153323. The T family was divided into 19 different MIRU genotypes.
Table 3
Degree of discrimination obtained with two typing methods individually and combined
Method
No. of different patterns
No. of unique isolates
No. of clustered isolates
No. of clusters
Spoligotyping alone
87
63(18%)
293(82%)
24
MIRU-VNTR typing alone
140
104(29%)
252(71%)
36
Spoligotyping+ MIRU-VNTR typing
181
120(33%)
236(67%)
43
A test population of 356 MTB isolates classified by spoligotyping analysis into six MTB lineages was typed using 15-locus (including ETR-A, B, C) MIRU-VNTR profiling. The combination of spoligotyping and MIRU-VNTR typing for the 356 isolates revealed that 120 (33%) have unique genotypes and 236 (67%) can be grouped into one of 43 genotypes. The minimum estimate for the proportion of TB in the study population due to recent transmission is estimated around 54% ([236–43]/356) (Table 3).
NTFlocus and RD deletion analysis
The RD105 (regions of difference) LSP appears in all Beijing strains examined to date and can be used to define the set of strains belonging to this lineage. Further subdivisions of the M. tuberculosis Beijing lineage are made on the basis of the variable appearance of the RD207, RD181, RD150, and RD142 deletions. In the present study, we refer to group 1 Beijing strains as those that contain only the RD105 deletion, while group 2 to 7 Beijing strains also contain RD105 and RD207 deletion (group 2); RD105, RD207, RD142, and RD150 deletion (group 3); RD105, RD207, and RD181 deletion (group 4); RD105, RD207, RD181, and RD150 deletion (group 5); RD105, RD207, RD181, and RD142 deletion (group 6); and RD105, RD207, RD181, RD150, and RD142 deletion (group 7). Strains in groups 1–3 have neither deletion of RD181 nor insertion of IS6110 in NTF region, and can be thought of as being "ancient" to the "modern" Beijing lineage (Figure 4). Results of RD analyses suggest that most of Beijing isolates in general population were modern strains (96%, 174/181).
Figure 4
Beijing family subgroup structure. Scheme of the proposed evolutionary pathway of the Beijing lineages, based on the deletion of genomic regions (RD, region of difference) showed inrectangles and on the NTF region with IS6110 insertions. The gray rectangles are lineages that have been reported previously, while the white rectangles represent proposed new lineages that are identified in the current study study.
Drug resistance patterns
Of the 356 strains in this study, 281 isolates (79%) were sensitive to all four of the first-line agents tested, and 75 were resistant to at least one drug, 2.8% are multidrug resistant (MDR) (Table 1). Analysis of the association between MDR and genotypes (as determined by spoligotyping) showed that the Beijing genotype is more likely to be MDR than all other genotypes (Haarlem, T, EAI, others, and orphan combined) [p = 0.08, OR = 3.73, and 95% C.I. = (0.78–17.83)]. The EAI family is significantly more likely to be sensitive to all drugs than are other genotypes [p = 0.02, OR = 3.64, and 95% C.I. = (1.09–12.15)]. Interestingly among the orphan strains, 5% were MDR and 20% were resistant to one drug, showing a distribution similar to that of the Beijing family.
Discussion
Tuberculosis remains a major public health issue in Taiwan and throughout the world. Recent advances in molecular technology such as spoligotyping and MIRU-VNTR typing provide a powerful tool to analyze MTB genotype and transmission patterns, which should be valuable for development of effective infection-control policy.
Among the 356 samples we analyzed in our study, 65% isolates would be considered to be potentially clustered based on 15 loci MIRU-VNTR typing [29]. Such rate is higher than that reported for other major cities such as San Francisco or New York, but similar to that from the Netherlands and Denmark, where the average of recent transmission rate ranged between 35% and 45% [30–33]. In densely populated areas, such as Northern Malawi, South African, clustering could be even higher, approaching 70% in some reports [34–36]. However, the rate of clustering in this study may have been overestimated due to the insufficient number of MIRU-VNTR loci we used, especially in an ares with high prevalence of Beijing strains[37]. Typing by additional loci, such as the 24 loci MIRU, is required to better differentiate the genetic relationship, clustering and possible transmission chain of all Beijing strains. Since typing methods of 24 loci MIRU were published in 2006, we attempted to set up the experiment protocol in our laboratory for the other 8 loci on top of the 16 loci we had used (the 15 shown in this paper plus the "MPTR-A"). Typing for 6 of the 8 loci was successfully established and performed on 52 Beijing strains that were from 15 clusters based on the 15 loci MIRU typing in this study. Our preliminary data using the additional 6 loci revealed that the same clustering patterns as that by the 16 loci were observed in 50% of the 52 Beijing strains analyzed. Therefore, the addition of the new VNTR loci to MIRU analysis is required to clarify the clustering and transmission issue in Taipei city in the future.
The modern MTB strains such as the Beijing, Haarlem, EAI, and T clusters comprise the causal agents of major epidemics. This study revealed the Beijing strain as the dominant pathogen for up to 52.5% of cases in Taipei, similar to the level in a recent report for Taiwan as a whole [38]. The majority of MTB strains of the Beijing family originated from the area in and around Beijing, China, and strains of this family were found to be dominant in neighboring countries such as Indonesia (44%), South Korea (72%), Thailand (44%), and Vietnam (53%) [3, 5, 8]. In contrast to the predominance of the Beijing genotype in many Asian countries, a low frequency (3%) of this genotype was reported among the strains in India [39]. Strains of the Beijing family have also been found in Europe, Africa, and the United States. The W strain, which caused a large outbreak of multidrug-resistant TB in New York and other U.S. cities, also belongs to the Beijing family [40]. In countries neighboring Taiwan, rates of infection with the Beijing family strains are higher than those in more distant countries, suggesting that the Beijing family may have radiated from the Beijing area to other regions. Based on the epidemiologic data, the Beijing strains appear to have a growth advantage over other strains, enabling them to circulate better in the population [3, 5, 8]. Moreover, we also demonstrated that Beijing family strains were associated with MDR phenotypes in this study (p = 0.08), a finding similar to that in the recent report from Taiwan [38]. Association between Beijing strains and MDR varies worldwide. Although such an association was reported in studies in the United States, Estonia, and Vietnam [41], it has not been noted in countries such as China and Indonesia, where representation of Beijing strains in the population is higher [42].
A total of 187 strains in the Beijing cluster identified by spoligotyping were further discriminated by MIRU-VNTR analysis. A total of 155 of the 187 strains were clustered into 15 groups, each consisting of 2 to 95 strains; and remaining 32 strains were found to have unique patterns. This study further showed that MTB isolates grouped into the Beijing family by spoligotyping have a similar grouping pattern when other genetic markers such as MIRU-VNTR typing are used. This was borne out by the fact that all MIRU-VNTR patterns of Beijing family strains were highly similar, differing only in copy numbers for one to three loci. In this study, all isolates that contain ≥ 2 repeats in the MIRU-VNTR locus 24 belong to the ancestral (TbD1+) group; and all but 1 isolate containing two repeat units in locus 24 belong to the modern (TbD1-) groups. In contrast, we found that ST480 of the U lineage has the MIRU-VNTR profile 254326223432 and TBD1+. The genetic characteristics of ST480 are very similar to those of the EAI family. Frequencies of TbD1+/EAI isolates have recently been reported to range from 25% to 50% in Bangladesh [43, 44] and Singapore [45]. A frequency of 8% has been reported in a study that used spoligotyping alone for genetic characterization of 105 isolates from the Delhi area [46]. Our analysis found that 5.5% of the samples were TbD1+/EAI isolates, while another 5.5% were TbD1+/new EAI spoligotype. The Haarlem isolates accounted for 13.5% in this study. Preliminary studies on the MTB strain distribution in eastern Taiwan's Hualien County, where Taiwanese Aborigines comprise a relatively large percentage of the population, showed a predominance of the Haarlem strain of up to 45% [47]. Since Taiwan was colonized by the Dutch in the 17th century, it is conceivable that the Haarlem strain is dominant in Taiwanese Aborigines. In addition to the Beijing and Haarlem strains, we identified the EAI family, the T family, the Latin American-Mediterranean family (LAM), the U family, and the MANU family of MTB in our population in Taipei. In addition to the identified predominant groups, we were able to identify the occurrence of rare clusters or localized STs listed in SpolDB4.0 that had previously been found in America, Australia, and Europe, with more found in countries neighboring Taiwan, such as Vietnam, Malaysia, the Philippines, and China.
TB occurs partly as a primary disease (typically defined as occurring within 5 years of infection) and partly as an endogenous reactivation or exogenous reinfection (occurring >5 years after infection) [48]. With increasing age, a decreasing proportion of cases are due to primary TB. Thus, the association of the Beijing genotype and young age suggests a recent spread of the Beijing genotype in Taipei. Anh et al. [49] reported that M. tuberculosis isolates of the Beijing genotype was less associated with BCG vaccination but was frequently associated with younger age in Vietnam. Lopez et al [50] used the mouse model of pulmonary tuberculosis to investigate the protective efficacy of BCG against these different strains and found that BCG was least protective against the Beijing strain. In contrast, Chan et al. did not find in Hong Kong any association between the Beijing genotype and younger age but did find a weak association with isoniazid (INH) resistance [51]. Although Taiwan executes comprehensively the BCG vaccination for more 40 years, the predominance of Beijing family strain in young population (85%, below 25 years of age) in this study suggest that BCG may fail to protect adequately the young people infected with the Beijing strain MTB.
The Beijing family can be further grouped into ancestral, modern, and recent strains by NTF locus analysis [8] and RD deletion analysis [4, 52], suggesting the strains' temporal evolution or transmission in migratory populations [16]. According to our previous study, the distribution of the Beijing sub-lineage with intact NTF region (ancient) was 19% in the general population, 24% in the veterans, and 50% in Aborigines in Taiwan [47]. We speculate the group 3 needed for 500 years from the evolution by group 1, group 4 to group 7 only need for 50 years, obviously modern strain genome was unstable and perhaps this instability was conducive to its fast spread. According to the view of evolution, existence of RD181 region or not, boundary to become modern and ancient lineage, infer RD181 perhaps (contain Rv2262c, Rv2263) the gene included may relate to the maintenance of genome stability. The hypothetical protein RV2262c may involve in protein modification and repair, and the hypothetical protein RV2263 involve in oxidoreduction http://cmr.jcvi.org/cgi-bin/CMR/shared/GenePage.cgi?locus=NTL02MT02256. The genome instability may be caused by these gene deletions.
Conclusion
This study gives a first overview of the M. tuberculosis strains circulating in metropolitan Taipei. Based on a combination of spoligotyping and MIRU-VNTR, our preliminary data showed that the Beijing strain has a high number of clusters in our sample population and this conclusion should be further clarified in the future using the 24 loci MIRU analysis. The high prevalence of Beijing genotype in young age population warrants a close attention to the control policy and the vaccine strategy. These findings indicate that TB is not optimally controlled in Taipei, and that efforts for control strategies should be enforced. Strain analysis, together with virulence studies, will also helping pinpointing isolates associated with higher morbidity and mortality, with the aim of directing efforts to limit the spread of those strains within the region.
Declarations
Acknowledgements
This project was supported by grants from the National Health Research Institutes and National Science Council (NSC97-3112-B-400-012), and the Department of Health (DOH97-DC-1501-01, from J.J.L.), Taiwan. We thank the mycobacteriology laboratory of Tri-Service General Hospital for providing bacterial isolates. We also thank Dr Daryl Henderson for his kind help in improving the English of this manuscript. All participants of this consortium are acknowledged for valuable discussions.
Authors’ Affiliations
(1)
Division of Clinical Research, National Health Research Institutes, Zhunan, Taiwan
(2)
Division of Clinical Pathology, Department of Pathology, Tri-Service General Hospital and National Defense Medical Center
(3)
Department of Laboratory Medicine, China Medical University Hospital
References
Center for Disease Control, Department of Health: Tuberculosis annual report. Taipei, Taiwan. 2002Google Scholar
Warren RM, Streicher EM, Sampson SL, Spuy van der GD, Richardson M, Nguyen D, Behr MA, Victor TC, van Helden PD: Microevolution of the direct repeat region of Mycobacterium tuberculosis : implications for interpretation of spoligotyping data. J Clin Microbiol. 2002, 40 (12): 4457-4465. 10.1128/JCM.40.12.4457-4465.2002.View ArticlePubMedPubMed CentralGoogle Scholar
Glynn JR, Whiteley J, Bifani PJ, Kremer K, van Soolingen D: Worldwide occurrence of Beijing/W strains of Mycobacterium tuberculosis : a systematic review. Emerg Infect Dis. 2002, 8 (8): 843-849.View ArticlePubMedPubMed CentralGoogle Scholar
Tsolaki AG, Gagneux S, Pym AS, Goguet de la Salmoniere YO, Kreiswirth BN, Van Soolingen D, Small PM: Genomic deletions classify the Beijing/W strains as a distinct genetic lineage of Mycobacterium tuberculosis. J Clin Microbiol. 2005, 43 (7): 3185-3191. 10.1128/JCM.43.7.3185-3191.2005.View ArticlePubMedPubMed CentralGoogle Scholar
Li WM, Wang SM, Li CY, Liu YH, Shen GM, Zhang XX, Niu TG, Gao Q, van Soolingen D, Kremer K, et al: Molecular epidemiology of Mycobacterium tuberculosis in China: a nationwide random survey in 2000. Int J Tuberc Lung Dis. 2005, 9 (12): 1314-1319.PubMedGoogle Scholar
van Soolingen D, Qian L, de Haas PE, Douglas JT, Traore H, Portaels F, Qing HZ, Enkhsaikan D, Nymadawa P, van Embden JD: Predominance of a single genotype of Mycobacterium tuberculosis in countries of east Asia. J Clin Microbiol. 1995, 33 (12): 3234-3238.PubMedPubMed CentralGoogle Scholar
Drobniewski F, Bala banova Y, Nikolayevsky V, Ruddy M, Kuznetzov S, Zakharova S, Melentyev A, Fedorin I: Drug-resistant tuberculosis, clinical virulence, and the dominance of the Beijing strain family in Russia. JAMA. 2005, 293 (22): 2726-2731. 10.1001/jama.293.22.2726.View ArticlePubMedGoogle Scholar
Mokrousov I, Ly HM, Otten T, Lan NN, Vyshnevskyi B, Hoffner S, Narvskaya O: Origin and primary dispersal of the Mycobacterium tuberculosis Beijing genotype: clues from human phylogeography. Genome Res. 2005, 15 (10): 1357-1364. 10.1101/gr.3840605.View ArticlePubMedPubMed CentralGoogle Scholar
Anh DD, Borgdorff MW, Van LN, Lan NT, van Gorkom T, Kremer K, van Soolingen D: Mycobacterium tuberculosis Beijing genotype emerging in Vietnam. Emerg Infect Dis. 2000, 6 (3): 302-305.View ArticlePubMedPubMed CentralGoogle Scholar
Almeida D, Rodrigues C, Ashavaid TF, Lalvani A, Udwadia ZF, Mehta A: High incidence of the Beijing genotype among multidrug-resistant isolates of Mycobacterium tuberculosis in a tertiary care center in Mumbai, India. Clin Infect Dis. 2005, 40 (6): 881-886. 10.1086/427940.View ArticlePubMedGoogle Scholar
Cox HS, Kubica T, Doshetov D, Kebede Y, Rusch-Gerdess S, Niemann S: The Beijing genotype and drug resistant tuberculosis in the Aral Sea region of Central Asia. Respir Res. 2005, 6: 134-10.1186/1465-9921-6-134.View ArticlePubMedPubMed CentralGoogle Scholar
Park YK, Shin S, Ryu S, Cho SN, Koh WJ, Kwon OJ, Shim YS, Lew WJ, Bai GH: Comparison of drug resistance genotypes between Beijing and non-Beijing family strains of Mycobacterium tuberculosis in Korea. J Microbiol Methods. 2005, 63 (2): 165-172. 10.1016/j.mimet.2005.03.002.View ArticlePubMedGoogle Scholar
Toungoussova OS, Caugant DA, Sandven P, Mariandyshev AO, Bjune G: Impact of drug resistance on fitness of Mycobacterium tuberculosis strains of the W-Beijing genotype. FEMS Immunol Med Microbiol. 2004, 42 (3): 281-290. 10.1016/j.femsim.2004.05.012.View ArticlePubMedGoogle Scholar
Mardassi H, Namouchi A, Haltiti R, Zarrouk M, Mhenni B, Karboul A, Khabouchi N, Gey van Pittius NC, Streicher EM, Rauzier J, et al: Tuberculosis due to resistant Haarlem strain, Tunisia. Emerg Infect Dis. 2005, 11 (6): 957-961.View ArticlePubMedPubMed CentralGoogle Scholar
Moss AR, Alland D, Telzak E, Hewlett D, Sharp V, Chiliade P, LaBombardi V, Kabus D, Hanna B, Palumbo L, et al: A city-wide outbreak of a multiple-drug-resistant strain of Mycobacterium tuberculosis in New York. Int J Tuberc Lung Dis. 1997, 1 (2): 115-121.PubMedGoogle Scholar
Kolk AH, Schuitema AR, Kuijper S, van Leeuwen J, Hermans PW, van Embden JD, Hartskeerl RA: Detection of Mycobacterium tuberculosis in clinical samples by using polymerase chain reaction and a nonradioactive detection system. J Clin Microbiol. 1992, 30 (10): 2567-2575.PubMedPubMed CentralGoogle Scholar
Kox LF, Rhienthong D, Miranda AM, Udomsantisuk N, Ellis K, van Leeuwen J, van Heusden S, Kuijper S, Kolk AH: A more reliable PCR for detection of Mycobacterium tuberculosis in clinical samples. J Clin Microbiol. 1994, 32 (3): 672-678.PubMedPubMed CentralGoogle Scholar
Brudey K, Driscoll JR, Rigouts L, Prodinger WM, Gori A, Al-Hajoj SA, Allix C, Aristimuno L, Arora J, Baumanis V, et al: Mycobacterium tuberculosis complex genetic diversity: mining the fourth international spoligotyping database (SpolDB4) for classification, population genetics and epidemiology. BMC Microbiol. 2006, 6: 23-10.1186/1471-2180-6-23.View ArticlePubMedPubMed CentralGoogle Scholar
Supply P, Allix C, Lesjean S, Cardoso-Oelemann M, Rusch-Gerdes S, Willery E, Savine E, de Haas P, van Deutekom H, Roring S, et al: Proposal for standardization of optimized mycobacterial interspersed repetitive unit-variable-number tandem repeat typing of Mycobacterium tuberculosis. J Clin Microbiol. 2006, 44 (12): 4498-4510. 10.1128/JCM.01392-06.View ArticlePubMedPubMed CentralGoogle Scholar
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. Proc Natl Acad Sci USA. 2001, 98 (4): 1901-1906. 10.1073/pnas.98.4.1901.View ArticlePubMedPubMed CentralGoogle Scholar
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 (3): 762-771. 10.1046/j.1365-2958.2000.01905.x.View ArticlePubMedGoogle Scholar
Brosch R, Gordon SV, Marmiesse M, Brodin P, Buchrieser C, Eiglmeier K, Garnier T, Gutierrez C, Hewinson G, Kremer K, et al: A new evolutionary scenario for the Mycobacterium tuberculosis complex. Proc Natl Acad Sci USA. 2002, 99 (6): 3684-3689. 10.1073/pnas.052548299.View ArticlePubMedPubMed CentralGoogle Scholar
Plikaytis BB, Marden JL, Crawford JT, Woodley CL, Butler WR, Shinnick TM: Multiplex PCR assay specific for the multidrug-resistant strain W of Mycobacterium tuberculosis. J Clin Microbiol. 1994, 32 (6): 1542-1546.PubMedPubMed CentralGoogle Scholar
Hirsh AE, Tsolaki AG, DeRiemer K, Feldman MW, Small PM: Stable association between strains of Mycobacterium tuberculosis and their human host populations. Proc Natl Acad Sci USA. 2004, 101 (14): 4871-4876. 10.1073/pnas.0305627101.View ArticlePubMedPubMed CentralGoogle Scholar
Tsolaki AG, Gagneux S, Pym AS, Goguet de la Salmoniere YO, Kreiswirth BN, Van Soolingen D, Small PM: Genomic deletions classify the Beijing/W strains as a distinct genetic lineage of Mycobacterium tuberculosis. J Clin Microbiol. 2005, 43 (7): 3185-3191. 10.1128/JCM.43.7.3185-3191.2005.View ArticlePubMedPubMed CentralGoogle Scholar
Tsolaki AG, Hirsh AE, DeRiemer D, Enciso JA, Wong MZ, Hannan M, Goguet de la Salmoniere YO, Aman K, Kato-Maeda M, Small PM: Functional and evolutionary genomics of Mycobacterium tuberculosis : insights from genomic deletions in 100 strains. Proc Natl Acad Sci USA. 2004, 101 (14): 4865-4870. 10.1073/pnas.0305634101.View ArticlePubMedPubMed CentralGoogle Scholar
National Committee for Clinical Laboratory Standards. Susceptibility testing of mycobacteria, nocardia and other aerobic actinomycetes. Wayne, PA. 2000Google Scholar
Murray M, Alland D: Methodological problems in the molecular epidemiology of tuberculosis. Am J Epidemiol. 2002, 155 (6): 565-571. 10.1093/aje/155.6.565.View ArticlePubMedGoogle Scholar
Small PM, Hopewell PC, Singh SP, Paz A, Parsonnet J, Ruston DC, Schecter GF, Daley CL, Schoolnik GK: The epidemiology of tuberculosis in San Francisco. A population-based study using conventional and molecular methods. N Engl J Med. 1994, 330 (24): 1703-1709. 10.1056/NEJM199406163302402.View ArticlePubMedGoogle Scholar
Alland D, Kalkut GE, Moss AR, McAdam RA, Hahn JA, Bosworth W, Drucker E, Bloom BR: Transmission of tuberculosis in New York City. An analysis by DNA fingerprinting and conventional epidemiologic methods. N Engl J Med. 1994, 330 (24): 1710-1716. 10.1056/NEJM199406163302403.View ArticlePubMedGoogle Scholar
Bauer J, Yang Z, Poulsen S, Andersen AB: Results from 5 years of nationwide DNA fingerprinting of Mycobacterium tuberculosis complex isolates in a country with a low incidence of M. tuberculosis infection. J Clin Microbiol. 1998, 36 (1): 305-308.PubMedPubMed CentralGoogle Scholar
van Soolingen D, Borgdorff MW, de Haas PE, Sebek MM, Veen J, Dessens M, Kremer K, van Embden JD: Molecular epidemiology of tuberculosis in the Netherlands: a nationwide study from 1993 through 1997. J Infect Dis. 1999, 180 (3): 726-736. 10.1086/314930.View ArticlePubMedGoogle Scholar
Glynn JR, Crampin AC, Yates MD, Traore H, Mwaungulu FD, Ngwira BM, Ndlovu R, Drobniewski F, Fine PE: The importance of recent infection with Mycobacterium tuberculosis in an area with high HIV prevalence: a long-term molecular epidemiological study in Northern Malawi. J Infect Dis. 2005, 192 (3): 480-487. 10.1086/431517.View ArticlePubMedGoogle Scholar
Godfrey-Faussett P, Sonnenberg P, Shearer SC, Bruce MC, Mee C, Morris L, Murray J: Tuberculosis control and molecular epidemiology in a South African gold-mining community. Lancet. 2000, 356 (9235): 1066-1071. 10.1016/S0140-6736(00)02730-6.View ArticlePubMedGoogle Scholar
Verver S, Warren RM, Munch Z, Vynnycky E, van Helden PD, Richardson M, Spuy van der GD, Enarson DA, Borgdorff MW, Behr MA, et al: Transmission of tuberculosis in a high incidence urban community in South Africa. Int J Epidemiol. 2004, 33 (2): 351-357. 10.1093/ije/dyh021.View ArticlePubMedGoogle Scholar
Supply P, Allix C, Lesjean S, Cardoso-Oelemann M, Rusch-Gerdes S, Willery E, Savine E, de Haas P, van Deutekom H, Roring S, et al: Proposal for standardization of optimized mycobacterial interspersed repetitive unit-variable-number tandem repeat typing of Mycobacterium tuberculosis. J Clin Microbiol. 2006, 44 (12): 4498-4510. 10.1128/JCM.01392-06.View ArticlePubMedPubMed CentralGoogle Scholar
Jou R, Chiang CY, Huang WL: Distribution of the Beijing family genotypes of Mycobacterium tuberculosis in Taiwan. J Clin Microbiol. 2005, 43 (1): 95-100. 10.1128/JCM.43.1.95-100.2005.View ArticlePubMedPubMed CentralGoogle Scholar
Mistry NF, Iyer AM, D'Souza DT, Taylor GM, Young DB, Antia NH: Spoligotyping of Mycobacterium tuberculosis isolates from multiple-drug-resistant tuberculosis patients from Bombay, India. J Clin Microbiol. 2002, 40 (7): 2677-2680. 10.1128/JCM.40.7.2677-2680.2002.View ArticlePubMedPubMed CentralGoogle Scholar
Kremer K, van Soolingen D, Frothingham R, Haas WH, Hermans PW, Martin C, Palittapongarnpim P, Plikaytis BB, Riley LW, Yakrus MA, et al: Comparison of methods based on different molecular epidemiological markers for typing of Mycobacterium tuberculosis complex strains: interlaboratory study of discriminatory power and reproducibility. J Clin Microbiol. 1999, 37 (8): 2607-2618.PubMedPubMed CentralGoogle Scholar
Caminero JA, Pena MJ, Campos-Herrero MI, Rodriguez JC, Garcia I, Cabrera P, Lafoz C, Samper S, Takiff H, Afonso O, et al: Epidemiological evidence of the spread of a Mycobacterium tuberculosis strain of the Beijing genotype on Gran Canaria Island. Am J Respir Crit Care Med. 2001, 164 (7): 1165-1170.View ArticlePubMedGoogle Scholar
Toungoussova OS, Sandven P, Mariandyshev AO, Nizovtseva NI, Bjune G, Caugant DA: Spread of drug-resistant Mycobacterium tuberculosis strains of the Beijing genotype in the Archangel Oblast, Russia. J Clin Microbiol. 2002, 40 (6): 1930-1937. 10.1128/JCM.40.6.1930-1937.2002.View ArticlePubMedPubMed CentralGoogle Scholar
Banu S, Gordon SV, Palmer S, Islam MR, Ahmed S, Alam KM, Cole ST, Brosch R: Genotypic analysis of Mycobacterium tuberculosis in Bangladesh and prevalence of the Beijing strain. J Clin Microbiol. 2004, 42 (2): 674-682. 10.1128/JCM.42.2.674-682.2004.View ArticlePubMedPubMed CentralGoogle Scholar
Shamputa IC, Rigouts L, Eyongeta LA, El Aila NA, van Deun A, Salim AH, Willery E, Locht C, Supply P, Portaels F: Genotypic and phenotypic heterogeneity among Mycobacterium tuberculosis isolates from pulmonary tuberculosis patients. J Clin Microbiol. 2004, 42 (12): 5528-5536. 10.1128/JCM.42.12.5528-5536.2004.View ArticlePubMedPubMed CentralGoogle Scholar
Sun YJ, Bellamy R, Lee AS, Ng ST, Ravindran S, Wong SY, Locht C, Supply P, Paton NI: 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 (5): 1986-1993. 10.1128/JCM.42.5.1986-1993.2004.View ArticlePubMedPubMed CentralGoogle Scholar
Singh UB, Suresh N, Bhanu NV, Arora J, Pant H, Sinha S, Aggarwal RC, Singh S, Pande JN, Sola C, et al: Predominant tuberculosis spoligotypes, Delhi, India. Emerg Infect Dis. 2004, 10 (6): 1138-1142.View ArticlePubMedPubMed CentralGoogle Scholar
Dou HY, Tseng FC, Lu JJ, Jou R, Tsai SF, Chang JR, Lin CW, Su IJ: Associations of Mycobacterium tuberculosis genotypes with different ethnic and migratory populations in Taiwan. Infection, Genetics and Evolution. 2008, 8: 323-330. 10.1016/j.meegid.2008.02.003.View ArticlePubMedGoogle Scholar
Vynnycky E, Nagelkerke N, Borgdorff MW, van Soolingen D, van Embden JDA, Fine PEM: The effect of age and study duration on the relationship between 'clustering' of DNA fingerprint patterns and the proportion of tuberculosis disease attributable to recent transmission. Epidemiol Infect. 2001, 126: 43-62. 10.1017/S0950268801225424.View ArticlePubMedPubMed CentralGoogle Scholar
Anh DD, Borgdorff MW, Van LN, Lan NTN, van Gorkom T, Kremer K, et al: Mycobacterium tuberculosis genotype Beijing emerging in Vietnam. Emerg Infect Dis. 2000, 6: 302-305.View ArticlePubMedPubMed CentralGoogle Scholar
Lopez B, Aguilar D, Orozco H, Burger M, Espitia C, Ritacco V, Barrera L, Kremer K, Hernandez-Pando R, Huygen K, et al: A marked difference in pathogenesis and immune response induced by different Mycobacterium tuberculosis genotypes. Clin Exp Immunol. 2003, 133 (1): 30-37. 10.1046/j.1365-2249.2003.02171.x.View ArticlePubMedPubMed CentralGoogle Scholar
Chan MY, Borgdorff M, Yip CW, de Haas PE, Wong WS, Kam KM, Van Soolingen D: Seventy percent of the Mycobacterium tuberculosis isolates in Hong Kong represent the Beijing genotype. Epidemiol Infect. 2001, 127 (1): 169-171. 10.1017/S0950268801005659.View ArticlePubMedPubMed CentralGoogle Scholar
Tsolaki AG, Hirsh AE, DeRiemer K, Enciso JA, Wong MZ, Hannan M, Goguet de la Salmoniere YO, Aman K, Kato-Maeda M, Small PM: Functional and evolutionary genomics of Mycobacterium tuberculosis : insights from genomic deletions in 100 strains. Proc Natl Acad Sci USA. 2004, 101 (14): 4865-4870. 10.1073/pnas.0305634101.View ArticlePubMedPubMed CentralGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
