BioMed Research International
Volume 2016 (2016), Article ID 6790985, 7 pages
http://dx.doi.org/10.1155/2016/6790985

Research Article

Molecular Typing Characteristic and Drug Susceptibility Analysis of Mycobacterium tuberculosis Isolates from Zigong, China

Hai-Can Liu,¹ Jian-Ping Deng,² Hai-Yan Dong,¹ Ti-Quan Xiao,² Xiu-Qin Zhao,¹ Zheng-Dong Zhang,² Yi Jiang,¹ Zhi-Guang Liu,¹ Qun Li,² and Kang-Lin Wan¹

¹State Key Laboratory for Infectious Diseases Prevention and Control, Collaborative Innovation Center for Diagnosis and Treatment of Infectious Diseases, National Institute for Communicable Disease Control and Prevention, Chinese Center for Disease Control and Prevention, Beijing 102206, China
²Zigong Center for Disease Control and Prevention, Zigong, Sichuan 643000, China

Received 11 November 2015; Revised 5 January 2016; Accepted 14 January 2016

Academic Editor: Peter F. Stadler

Copyright © 2016 Hai-Can Liu 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.

China is one of the 22 countries with high TB burden worldwide, and Sichuan contained the second-largest number of TB cases among all of the Chinese provinces. But the characteristics of Mycobacterium tuberculosis circulated in Zigong, Sichuan, were still unknown. To investigate the character and drug resistance profile, 265 clinical isolates were cultured from tuberculosis patient’s sputum samples in the year of 2010, of which the genetic profile was determined by using Spoligotyping and MIRU-VNTR typing methods, and the drug sensibility testing to the four first-line and four second-line antituberculosis (anti-TB) drugs was performed by using proportion method on Lowenstein-Jensen (L-J) media. The major Spoligotype was Beijing family (143/265, 53.96%), followed by T (80/265, 30.19%) and H (9/265, 3.40%) genotypes; the total Hunter-Gaston discrimination index (HGDI) of the 24 loci MIRU-VNTR was 0.9995. About 27.17% (72/265) of the isolates were resistant to at least one of the eight tested anti-TB drugs, and for Beijing and non-Beijing family isolates the proportion of drug resistance was 28.47% (41/144) and 25.62% (31/121), respectively. That is, the most prevalent genotype here was Beijing family, and the 24 loci VNTR analysis could supply a high resolution for genotyping, and Beijing and non-Beijing isolates had no difference () for drug resistance.

1. Introduction

Tuberculosis (TB) remains an important problem for the public health worldwide, especially in developing countries. Of the five countries with the largest number of incident cases in 2012, China (0.9 million–1.1 million) ranked the second just next to India, which accounted for 12% of the global cases alone [1]. Recently the emerging of drug-resistant TB (DR-TB) is becoming an important threat to TB control and public health worldwide [2, 3], and according to the National Survey of Drug-Resistant Tuberculosis in China, 5.7% of new TB patients and 25.6% of the previously treated cases were multi-drug-resistant TB (MDR-TB) [4]. In the last decade, genotyping of Mycobacterium tuberculosis (M. tuberculosis) has significantly enhanced our understanding of TB epidemiology by detecting suspected outbreaks and by tracing the infection routine. Such molecular methods like spacer oligonucleotide typing (Spoligotyping) and mycobacterial interspersed repetitive unit- (MIRU-) variable number tandem repeats (VNTR) typing have been used to characterize clinical M. tuberculosis isolates in China [5–7]. Obtaining the characteristic of M. tuberculosis isolates circulating in an area is helpful for understanding and controlling the spread of TB infection.

Although molecular epidemiology study of TB has been taken out in some area in China [8–10], few efforts have been made for the collection of DNA fingerprinting data based on the mutation profiles of isolates from the western region of China such as Sichuan province. Sichuan is a province that is located in southwestern China, with a population of approximately 87 million inhabitants. Sichuan contained the second-largest number of TB cases among Chinese provinces, where the prevalence of both TB and drug-resistant TB was much higher than the average level in China [11]. The total cases of TB in Sichuan province are estimated to be about 272,000 in 2006 and the new cases found in 2006 were more than 70,000 [12]. So the genetic and drug-resistant characteristics analysis of clinical M. tuberculosis isolates from this area may give us some benefits on TB control in China.

The aim of this study was to identify and genetically characterize the M. tuberculosis isolates circulating in Zigong, southern Sichuan, and to explore the distribution of drug resistance profiles across the major Spoligotyping-defined M. tuberculosis genotypes.

2. Material and Methods

2.1. Bacterial Strains and DNA Samples

A total of 265 clinical M. tuberculosis isolates were collected from Zigong, Sichuan, in 2010. Chromosomal DNA was extracted from fresh bacteria cultures on the Lowenstein-Jensen (L-J) medium. A loopful of mycobacterial colonies were suspended in 400 μL TE buffer (pH 8.0) and then boiled for 10 min. The suspension was centrifuged at 12,000 rpm for 10 min, and the supernatant was served as PCR template and stored at −20°C.

2.2. Spoligotyping

Spoligotyping was performed as described previously by Kamerbeek et al. [13]. The results were entered into Excel spreadsheets in a binary format and compared with SITVITWEB database [14], which is a publicly available international multimarker database for studying Mycobacterium.

2.3. MIRU-VNTR Typing

VNTR analysis was performed by the amplification of 24 genomic loci, the scheme of which has been described by Supply et al. [15]. Each VNTR locus was amplified individually in 20 μL reaction mixture containing 1x PCR GC Buffer I (Takara Bio Inc.), 5 pmol of each primer set, 200 μM of each dNTP, 0.5 U of DNA Taq Polymerase (Takara Bio Inc.), and 1.5 ng of DNA template. An initial denaturation step of 5 min at 95°C was followed by 35 cycles at 94°C for 1 min, 60°C for 1 min, and 72°C for 1 min, with a final extension at 72°C for 10 min. The PCR products were analyzed by electrophoresis on 2% agarose gel using 100 bp DNA ladder as size markers. The H37Rv strain was included as a positive control for each amplification reaction and run as an additional control for accuracy, and then the copy number of each locus was calculated by using the Image Lab software (Bio-Rad).

2.4. Drug Susceptibility Testing

Drug susceptibility testing of those 265 strains to the four first-line drugs (rifampicin [RFP], isoniazid [INH], streptomycin [SM], and ethambutol [EMB]) and four second-line drugs (amikacin [AMK], kanamycin [KA], ofloxacin [OFLX], and capreomycin [CPM]) was performed by using proportion method recommended by WHO [16]. Multi-drug-resistant TB (MDR-TB) is the isolate that was resistant to at least the two first-line anti-TB drugs, RFP and INH, based on the WHO definition [17].

2.5. Computer Analysis

The Spoligotyping and MIRU-VNTR typing results of all the M. tuberculosis clinical isolates were analyzed by the BioNumerics program (Windows 7, version 5.10; Applied Maths, Kortrijk, Belgium). Clusters were defined as groups that the M. tuberculosis isolates have identical Spoligotyping or MIRU-VNTR patterns. The dendrograms based on the Spoligotyping and 24 VNTR loci data were conducted using UPGMA protocol. The Hunter-Gaston discrimination index (HGDI) for the 24 VNTR loci was calculated using the previously reported formula [18]. The HGDI of each VNTR locus was calculated by using the web based VNTR Diversity and Confidence Extractor (V-DICE http://www.hpa-bioinformatics.org.uk/cgi-bin/DICI/DICI.pl).

2.6. Statistical Analysis

Chi-squared tests were used to compare the proportions of DR-TB and the distribution of genotypes in different groups. All tests of significance were two-sided and the significant threshold was set at 0.05.

3. Results

3.1. Demographic Characters

All of the 265 clinical isolates were cultured from sputum samples of different TB patients in Zigong CDC, Sichuan province, in 2010. And for the 265 TB patients, 44 (16.60%) were previously treated, 139 (52.45%) were in the 30- and 60-year-old age group, 39 (14.72%) female and 209 (78.87%) male patients with other 17 (6.41%) gender uncertain patients, and no difference about these factors between Beijing and non-Beijing family strains was found at level (Table 1).

Table 1: Demographic data of the 265 isolates in Zigong, China.

Factors
Whole
()
Beijing
()
Non-Beijing
()
( value)

Sex

    Female
39
25
14
1.6036 (0.2054)

    Male
209
111
98

    Unknown
17
7
10

Age

    <30
27
19
8
6.9592 (0.0732)

    ⩾30
139
79
60

    ⩾60
69
30
39

    ⩾75
9
6
3

    Unknown
21
9
12

Treated before?

    Yes
44
25
19
0.1732 (0.6773)

    No
221
118
103

: evaluated the difference between Beijing and non-Beijing family strains.

3.2. Spoligotyping

Based on the Spoligotyping results, all of the 265 isolates were subdivided into 64 Spoligotypes, among which 49 patterns were unique whereas other 216 isolates were grouped into 15 clusters containing 2 to 128 isolates. Following comparison with the SITVITWEB database, 239 isolates were clustered into 40 shared international types (SITs), whereas the other 26 isolates that had not been described in the database were referred to as “Orphan.” Family assignment revealed that 143 strains had the classical Beijing family-specific or Beijing family-like Spoligotyping pattern. We treated all these 143 isolates as Beijing family strains (Figure 1). In total, the Beijing family is clearly the most prevalent genetic family in our research setting (143/265, 53.96%), followed by T family (80/265, 30.19%) and H family (9/265, 3.40%).

Figure 1

3.3. MIRU-VNTR Typing

Exploitable VNTR results (defined as less than 4 missing values) were obtained for 262 isolates out of 265 isolates. Using the 24 polymorphic loci, all of the 262 M. tuberculosis isolates were identified to 257 different VNTR types and 3 main clusters. A total of 253 isolates had unique profiles. The remaining 9 isolates formed 4 clusters (2 to 3 isolates) (Figure 2). The total Hunter-Gaston discrimination index (HGDI) for all of the 263 isolates was 0.9995, and that of Beijing and non-Beijing family isolates was 0.9987 and 0.9992, respectively; the HGDI of each of the 24 VNTR loci were shown in Table 2.

Table 2: HGDI of the 24 MIRU-VNTR loci for the whole sample, Beijing and non-Beijing family isolates.

Locus
Whole
()
Beijing
()
Non-Beijing
()

ETRA
0.594
0.262
0.606

ETRB
0.374
0.068
0.504

ETRC
0.110
0.069
0.159

ETRD
0.303
0.028
0.540

ETRE
0.592
0.245
0.202

MIRU02
0.030
0.000
0.065

MIRU10
0.576
0.263
0.242

MIRU16
0.283
0.181
0.392

MIRU20
0.246
0.240
0.254

MIRU23
0.124
0.108
0.143

MIRU24
0.015
0.014
0.017

MIRU26
0.739
0.570
0.598

MIRU27
0.328
0.028
0.523

MIRU39
0.521
0.082
0.097

MIRU40
0.551
0.485
0.619

Mtub04
0.620
0.363
0.622

Mtub21
0.806
0.710
0.571

Mtub29
0.096
0.069
0.127

Mtub30
0.496
0.119
0.182

Mtub34
0.095
0.042
0.157

Mtub39
0.511
0.267
0.647

QUB11b
0.854
0.745
0.796

QUB26
0.723
0.571
0.821

QUB4156
0.052
0.081
0.017

Figure 2

3.4. Drug-Resistant Phenotypic Profiles

Among all the 265 isolates, 72 (27.17%) were resistant to at least one of the eight tested anti-TB drugs, and other 193 (72.83%) were sensitive. Of the 72 drug-resistant isolates, 48 (18.11%) were INH-resistant, 40 (15.09%) were RIF-resistant, 32 (12.08%) were SM-resistant, 17 (6.42%) were EMB-resistant, 32 (12.08%) were OFLX-resistant, and 36 (13.58%) of those drug-resistant strains were MDR-TB; no XDR-TB was found in this study (Table 3). When grouped by gender or age group like previous report [8], no relationship was found with the isolate’s character of drug resistance. But having been treated before may play an important role with the drug resistance, not only for the drug-resistant TB (, ) but also for the MDR-TB (, ) (Table 4).

Table 3: Prevalence of anti-TB drug resistance among 265 isolates in Sichuan province.

Anti-TB drugs
Whole
()
Beijing
()
Non-Beijing
()
( value)

Isoniazid
18.11% ()
19.44% ()
16.53% ()
0.45 (0.5019)

Rifampicin
15.09% ()
18.75% ()
10.74% ()
3.48 (0.0623)

Streptomycin
12.08% ()
13.89% ()
9.92% ()
1.07 (0.3014)

Ethambutol
6.42% ()
7.64% ()
4.96% ()
0.84 (0.3583)

Ofloxacin
12.08% ()
9.03% ()
15.70% ()
2.61 (0.1065)

Resistant to any anti-TB drugs
27.17% ()
28.47% ()
25.62 % ()
0.35 (0.5519)

MDR-TB
13.58% ()
15.97% ()
10.74% ()
1.65 (0.1986)

: all 265 isolates were sensitive to other three kinds of anti-TB drugs, and Pearson and value were evaluated between the Beijing and non-Beijing family strains.

Table 4: Prevalence of anti-TB drug resistance among 265 isolates by factors.

Factors
Sensitive TB ()
Any drug resistance TB ()
( value)
Non-MDR-TB ()
MDR-TB ()
( value)

Sex

    Female
29
10
0.1097 (0.7405)
34
5
0.0095 (0.9225)

    Male
150
59
181
28

    Unknown
14
3

14
3

Age

    <30
22
5
4.1444 (0.2402)
23
4
4.8455 (0.1579)

    ⩾30
93
46
115
24

    ⩾60
52
17
64
5

    ⩾75
8
1
9
0

    Unknown
18
3

18
3

Treated?

    Yes
23
21
11.2680 (0.0008)
29
15
18.8983 (0.0000)

    No
170
51
200
21

: evaluated by the Fisher’s exact test.

3.5. Relationship between Beijing Family and Drug Resistance

According to the demographic data, no matter whether they are grouped by gender, age groups, or even treatment history, there was no difference of the infection rate between Beijing and non-Beijing family isolates observed (Table 1). And for the Spoligotyping patterns and DST results, of the 143 Beijing family isolates, 41 (28.67%) were drug-resistant to at least one of the eight anti-TB drugs of which 23 (16.08%) were MDR-TB. While, in the 122 non-Beijing isolates, 31 (25.41%) were drug-resistant and 13 (10.66%) were MDR-TB, no difference of the drug resistance proportion was found between Beijing and non-Beijing family isolates (Table 3).

4. Discussion

Previous research report has reviewed that Beijing family is prevalent (57.89%) in Sichuan province [6], and according to the Spoligotyping results in this study, the proportion of Beijing family is about 54.34% in southern Sichuan, a little lower than that of the previous report. This may be because our study was just focused on the prevalence isolates in Zigong, and more isolates were included in this study which could afford much more representativeness. But this still conforms to the former reports that the most prevalent isolate here was Beijing family strains and the highest prevalence of Beijing family was found in northern China, followed by central and southern China [6, 7]. Besides, T genotype was the second major epidemic genotype circulating in this area and was also consistent with the previous study [19]. This may indicate that non-Beijing family strains also played an important role in the high tuberculosis burden in Sichuan and some other provinces that are located in southern China [6, 7].

Spoligotyping has been considered as the gold standard for Beijing family identification, as it is simple and efficient. But the discrimination of Spoligotypes for Beijing family is very low and it cannot further differentiate Beijing family strains, as Beijing family strains always show characteristic homological Spoligotyping pattern [20]. For the 24 loci MIRU-VNTR typing method, even the discrimination of each locus is different, but when we combined them together a high resolution was acquired (Table 2). As the MIRU-VNTR genotyping is considered to be faster and easier to perform, the results are in digital format which can be easily compared between laboratories. Thus, if conditions allow, we suggest Spoligotyping and MIRU-VNTR genotyping should be taken out at the same time, to investigate in more detail characters of the selected isolates. But for some loci the HGDI were so low; different countries could select specific loci for genotyping. And for the three main clusters of the UPGMA dendrogram in Figure 2, one formed by Beijing family isolates and other two within non-Beijing family isolates, but no difference was found between them for difference factors, like drug resistance and patient age group.

In our study, 27.17% () of the isolates were resistant to at least one of the eight anti-TB drugs, and 13.58% () were MDR-TB; no XDR-TB was found in this study. These have some difference with the whole nation’s status [4]. Firstly, we focused on the target area here in Sichuan province; the demographic and M. tuberculosis infectious status in this area may have some differences with those of the whole nation. Secondly, the isolates selected into our study may also have some differences with that of the previously report. Also in our study we found that the treatment history may play an important role with the drug resistance, not only for drug-resistant TB but also for MDR-TB (Table 4); thus the management of TB patients should be strengthened. And as the drug-resistant TB may be the most important challenge to achieve the goals of stopping TB, much more attention should be paid to it.

As Beijing family strains were prevalent in this area (143/265, 53.96%) and in some reports, Beijing family has been considered to be related to TB outbreak and drug resistance [21–24], while some are not [25–27]. Referring to our study, the chi-square analysis showed that there was no significant difference between the Beijing and non-Beijing family strains in terms of their drug resistance proportion and patient age group. This may be because Beijing family is the genetic aspect but the drug resistance is a biological aspect of one M. tuberculosis isolate, and there is no constant relationship between them.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors thank the staffs of Zigong Center for Disease Control and Prevention for their excellent contribution to this study. This work was funded by the Chinese National Key Program of Mega Infectious Disease Project 2013ZX10003006-002-001.

References

World Health Organization, Global Tuberculosis Report 2013, World Health Organization, Geneva, Switzerland, 2013.

N. R. Gandhi, P. Nunn, K. Dheda et al., “Multidrug-resistant and extensively drug-resistant tuberculosis: a threat to global control of tuberculosis,” The Lancet, vol. 375, no. 9728, pp. 1830–1843, 2010.

World Health Organization, Multidrug and Extensively Drug-Resistant TB (M/XDR-TB): 2010 Global Report on Surveillance and Response, World Health Organization, 2010.

Y. Zhao, S. Xu, L. Wang et al., “National survey of drug-resistant tuberculosis in China,” The New England Journal of Medicine, vol. 366, no. 23, pp. 2161–2170, 2012.

Y. Jiang, H.-C. Liu, H. Zheng et al., “19-VNTR loci used in genotyping Chinese clinical Mycobacterium tuberculosis complex strains and in association with spoligotyping,” Journal of Basic Microbiology, vol. 53, no. 7, pp. 562–580, 2013.

H. Dong, Z. Liu, B. Lv et al., “Spoligotypes of Mycobacterium tuberculosis from different provinces of China,” Journal of Clinical Microbiology, vol. 48, no. 11, pp. 4102–4106, 2010.

H. Dong, L. Shi, X. Zhao et al., “Genetic diversity of Mycobacterium tuberculosis isolates from Tibetans in Tibet, China,” PLoS ONE, vol. 7, no. 3, Article ID e33904, 2012.

Q. Liu, D. Yang, W. Xu et al., “Molecular typing of Mycobacterium tuberculosis isolates circulating in Jiangsu province, China,” BMC Infectious Diseases, vol. 11, article 288, 2011.

B. Lu, P. Zhao, B. Liu et al., “Genetic diversity of Mycobacterium tuberculosis isolates from Beijing, China assessed by Spoligotyping, LSPs and VNTR profiles,” BMC Infectious Diseases, vol. 12, article 372, 2012.

Q. Yu, Y. Su, B. Lu et al., “Genetic diversity of Mycobacterium tuberculosis isolates from Inner Mongolia, China,” PLoS ONE, vol. 8, no. 5, Article ID e57660, 2013.

L. Wang, S. Cheng, M. Chen et al., “The fifth national tuberculosis epidemiological survey in 2010,” Chinese Journal of Antituberculosis, vol. 34, pp. 485–508, 2012.

J.-H. Guo, W.-L. Xiang, Q.-R. Zhao et al., “Molecular characterization of drug-resistant Mycobacterium tuberculosis isolates from Sichuan Province in China,” Japanese Journal of Infectious Diseases, vol. 61, no. 4, pp. 264–268, 2008.

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.

C. Demay, B. Liens, T. Burguière et al., “SITVITWEB—a publicly available international multimarker database for studying Mycobacterium tuberculosis genetic diversity and molecular epidemiology,” Infection, Genetics and Evolution, vol. 12, no. 4, pp. 755–766, 2012.

P. Supply, E. Mazars, S. Lesjean, V. Vincent, B. Gicquel, and C. Locht, “Variable human minisatellite-like regions in the Mycobacterium tuberculosis genome,” Molecular Microbiology, vol. 36, no. 3, pp. 762–771, 2000.

World Health Organization, Guidelines for Surveillance of Drug Resistance in Tuberculosis, World Health Organization, 2009.

World Health Organization, Guidelines for the Programmatic Management of Drug-Resistant Tuberculosis-2011 Update, World Health Organization, Geneva, Switzerland, 2011.

P. R. Hunter and M. A. Gaston, “Numerical index of the discriminatory ability of typing systems: an application of Simpson's index of diversity,” Journal of Clinical Microbiology, vol. 26, no. 11, pp. 2465–2466, 1988.

K. Wan, J. Liu, Y. Hauck et al., “Investigation on Mycobacterium tuberculosis diversity in china and the origin of the Beijing clade,” PLoS ONE, vol. 6, no. 12, Article ID e29190, 2011.

D. van Soolingen, L. Qian, P. E. W. de Haas et al., “Predominance of a single genotype of Mycobacterium tuberculosis in countries of east Asia,” Journal of Clinical Microbiology, vol. 33, no. 12, pp. 3234–3238, 1995.

K. Kremer, J. R. Glynn, T. Lillebaek et al., “Definition of the Beijing/W lineage of Mycobacterium tuberculosis on the basis of genetic markers,” Journal of Clinical Microbiology, vol. 42, no. 9, pp. 4040–4049, 2004.

Y. Hu, X. Ma, E. A. Graviss, W. Wang, W. Jiang, and B. Xu, “A major subgroup of Beijing family Mycobacterium tuberculosis is associated with multidrug resistance and increased transmissibility,” Epidemiology and Infection, vol. 139, no. 1, pp. 130–138, 2011.

S. Purwar, S. Chaudhari, V. M. Katoch et al., “Determination of drug susceptibility patterns and genotypes of Mycobacterium tuberculosis isolates from Kanpur district, North India,” Infection, Genetics and Evolution, vol. 11, no. 2, pp. 469–475, 2011.

T. Brown, V. Nikolayevskyy, P. Velji, and F. Drobniewski, “Associations between Mycobacterium tuberculosis strains and phenotypes,” Emerging Infectious Diseases, vol. 16, no. 2, pp. 272–280, 2010.

A. Ani, T. Bruvik, Y. Okoh et al., “Genetic diversity of Mycobacterium tuberculosis complex in Jos, Nigeria,” BMC Infectious Diseases, vol. 10, article 189, 2010.

M. K. Kai, W. Y. Chi, W. T. Lai et al., “Utility of mycobacterial interspersed repetitive unit typing for differentiating multidrug-resistant Mycobacterium tuberculosis isolates of the Beijing family,” Journal of Clinical Microbiology, vol. 43, no. 1, pp. 306–313, 2005.

J. Zhang, L. Mi, Y. Wang et al., “Genotypes and drug susceptibility of Mycobacterium tuberculosis Isolates in Shihezi, Xinjiang Province, China,” BMC Research Notes, vol. 5, article 309, 2012.