Molecular typing of mycobacterium tuberculosis isolates circulating in Jiangsu Province, China

Qiao Liu; Dandan Yang; Weiguo Xu; Jianming Wang; Bing LV; Yan Shao; Honghuan Song; Guoli Li; Haiyan Dong; Kanglin Wan; Hua Wang

doi:10.1186/1471-2334-11-288

Research article
Open Access
Open Peer Review

Molecular typing of mycobacterium tuberculosis isolates circulating in Jiangsu Province, China

Qiao Liu†¹,
Dandan Yang†²,
Weiguo Xu²,
Jianming Wang¹,
Bing LV³,
Yan Shao²,
Honghuan Song²,
Guoli Li²,
Haiyan Dong³,
Kanglin Wan³ and
Hua Wang^{1, 2}Email author

†Contributed equally

BMC Infectious Diseases201111:288

https://doi.org/10.1186/1471-2334-11-288

Received: 24 May 2011
Accepted: 26 October 2011
Published: 26 October 2011

Open Peer Review reports

Abstract

Background

Globally, China is the second place with high burden of tuberculosis (TB). To explore the characteristics of the pathogens of Mycobacterium tuberculosis (MTB) circulating in this area is helpful for understanding and controlling the spread of the strains. Recent developments in molecular biology have allowed prompt identification and tracking specific strains of MTB spreading through the population.

Methods

Spacer-oligonucleotide typing (spoligotyping) and mycobacterial interspersed repetitive units variable number tandem repeat (MIRU-VNTR) were performed in combination to yield specific genetic profiles of 260 MTB strains isolated from 30 counties of Jiangsu province in China between June and July 2010. The spoligotyping results were in comparison to the world Spoligotyping Database of Institute Pasteur de Guadeloupe (SpolDB4). Drug susceptibility test (DST) was performed on all strains by proportion method on Lowenstein-Jensen (LJ) culture media.

Results

Based on the spoligotyping method, 246 strains displayed known patterns and 14 were absent in the database. Predominant spoligotypes belonged to the Beijing family (80.4%). By using the 24-loci VNTR typing scheme, 224 different patterns were identified, including 20 clusters and 204 unique patterns. The largest clade comprised 195 strains belonging to the Beijing family. The combination of spoligotyping and 24-loci MIRU-VNTR demonstrated maximal discriminatory power. Furthermore, we observed a significant association between Beijing family strains and drug-resistant phenotypes. The Beijing family strains presented increased risks for developing multi-drug resistant TB, with the OR (95% CI) of 11.07(1.45-84.50).

Conclusions

The present study demonstrated that Beijing family isolates were the most prevalent strains circulating in Jiangsu province of China. The utility of spoligotyping in combination with 24-loci MIRU-VNTR might be a useful tool for epidemiological analysis of MTB transmission.

Keywords

Variable Number Tandem Repeat
Drug Susceptibility Test
Sputum Smear Microscopy
Beijing Genotype
Beijing Family

Background

Tuberculosis (TB) has been a global threat to human health and it remains a major public health burden in developing countries. As a public health dilemma, drug resistance has been an obstacle to achieve the goal of effective TB control. In 2009, it was reported that about 9.4 million incident cases (8.9 million-9.9 million) occurred globally and among them, 250 000 (230 000-270 000) were multi-drug resistant TB (MDR-TB) [1]. China is one of the five countries with the largest number of both incident cases (1.1-1.5 million) and MDR-TB cases (79 000-120 000) in absolute terms in 2008 [1]. The increased rate of drug-resistance and MDR strains of TB remain a serious problem of TB control in China.

MDR-TB may arise by inadequate antituberculosis treatment or by direct transmission of drug resistant strains from one individual to another [2]. The association with MDR-TB has led to concerns that Mycobacterium tuberculosis (MTB) strains may have a predilection for acquiring drug resistance. DNA fingerprinting of MTB isolates is a powerful tool to study the molecular epidemiological characteristics of MDR-TB. The traditional method used to genotype MTB complex (MTBC) is insertion sequence (IS) 6110 restriction fragment-length polymorphism (RFLP), which can be utilized in outbreak investigations, long-term surveillance and the detection of laboratory contamination. Despite the widespread use of IS6110-based RFLP analysis, it remains a labour intensive, technically demanding, and insufficiently discriminating strains with the low copy numbers (fewer than six) of IS6110 [3]. Following IS6110-based RFLP, spoligotyping has been the most widely used method for secondary typing of MTB strains and it has been particularly useful for identifying strains belonging to the Beijing/W family due to the absence of spacers 1-34 in the direct repeat (DR) region of the MTB genome [4]. Variable number tandem repeat (VNTR) is a newly developed typing method, which determines the number of repeated mycobacterial interspersed repetitive units (MIRU) by using small quantities of bacteria [5]. MIRU-VNTR has shown a higher discriminatory power when combining with spoligotyping. The advantages of such PCR-based genotyping techniques include their easy digitalization of the generated profiles and hence easy inter-laboratory comparison, as well as easy creation and maintenance of the database [6].

The Beijing family, first identified in 1995 in Beijing of China, almost omnipresent and significantly prevalent in certain world regions, e.g., East Asia and the former Union of Soviet Socialist Republics(USSR) [7]. Several studies have observed that the Beijing family MTB strains had exhibited important pathogenic features which might be associated with the drug resistant TB in China [6, 8–10]. As the prevalence of drug resistant clones of MTB varies from one area to another, studies in geographical distribution of resistant clones are helpful for understanding the epidemiological characteristics of TB.

Jiangsu is a province locating along the eastern coast of China and it covers an area of 102.6 thousand square kilometers, with a total population of 77 million in 2009. The high prevalence of drug resistance among TB patients in Jiangsu province has been a major challenge for TB control. According to the most recently updated data, the proportion of primary and acquired MDR-TB were 7.63% and 33.7%, respectively [11].

The main goal of this study was to characterize the genotyping of MTB strains circulating in Jiangsu Province and to explore whether there is any evidence that the spreading of Beijing family of MTB strains was associated with the drug resistance.

Methods

Mycobacterial strains

The study was conducted in 30 counties of Jiangsu province, which were sampled according to the "Guidelines for surveillance of drug resistance in tuberculosis" developed by WHO/IUATLD [12]. Two hundred and sixty newly diagnosed sputum smear positive pulmonary TB patients from June 2010 to July 2010 in the study sites were recruited. Totally, there were 216 new cases and 44 previously treated ones. The average age of cases was 48.5 years old. All TB patients were diagnosed by referring to the national guidelines of China and reviewed by Jiangsu provincial Center for Disease Control and Prevention. Each subject present three sputum smear samples with labeled plastic bottles for sputum smear microscopy test. The two sputum samples with the highest bacterial counts were cultured, and one culture was submitted to the provincial reference laboratory for drug susceptibility test (DST). Sputum smear microscopy and culture were performed at the level of county (district) laboratory.

Strains isolation and drug susceptibility test

The sputum samples were cultured and isolated on Lowenstein-Jensen (LJ) culture media. LJ medium was impregnated with isoniazid (INH), rifampicin (RIF), streptomycin (SM), and ethambutol (EMB). The concentrations of anti-tuberculosis drugs were 0.2 μg/ml for INH, 40 μg/ml for RIF, 4 μg/ml for SM, and 2 μg/ml for EMB. The strain was declared resistant to the specific drug; or it was defined as sensitive when the growth rate was < 1% compared to the control. MDR-TB was defined as strains being resistant to at least RIF and INH [11]. Strains isolation, identification and drug susceptibility test (DST) were performed at the provincial reference laboratory.

Genomic DNA extraction

Mycobacterial genomic DNA was extracted from mycobacterial colonies growing on LJ medium by resuspending one loop of mycobacterial colonies in 200 μl TE buffer (10 mM Tris-HCl, 1 mM EDTA) and was incubated at 85°C for 30 minutes to obtain genomic DNA. After centrifugation of the suspension, the supernatant fluid containing DNA was removed and stored at -20°C until further use [13]. Laboratory strain MTB H37Rv was used as a control for all microbiological and genetic procedures.

Genotyping

Spoligotyping of the isolates was performed as described by Kamerbeek et al [4]. In brief, genomic DNA of MTB isolates were amplified using polymerase chain reaction (PCR) with the primers of Dra (biotin labelled) and Drb, and then the PCR products were hybridized to a set of 43 oligonucleotide probes corresponding to each spacer, which were covalently bound to a membrane [4]. Spoligotypes in binary format were compared with the SpolDB4 database. In particular, strains with spoligotype patterns characterized by deletion of spacers 1-34 were defined as "typical" Beijing genotypes, whereas strains with additional deletion of one or more of the last nine spacers were defined as Beijing-like genotypes according to the criteria of the international database SpolDB4 [4, 14–16].

To identify a suitable MIRU-VNTR loci set for genotyping MTB in this area, the number of tandem repeats was determined in 24 MIRU-VNTR genetic loci: ten original MIRU-VNTR loci; six loci of exact tandem repeats (ETRs): ETR-A, -B, -C, -D, -E and -F; five Mtub loci: Mtub4, 21, 30, 38, and 39; and three Queen's University of Belfast (QUBs) loci: QUB-11b, -26, and 4156c. PCR amplification was performed for correspondent loci as described in the previous studies [17–20]. PCR products were analyzed on a 2% agarose gel against a 100-bp DNA ladder, and the copy number at each locus was calculated using the Quantity 1 gel imaging system [10]. Determination of the discriminatory power of the VNTR loci was calculated using the Hunter-Gaston discriminatory Index (HGDI) [21]. The HGDI was calculated using the following formula:

$H G D I = 1 - [\frac{1}{N (N - 1)} \sum_{j = 1}^{s} n j (n j - 1)]$ , where D is the numerical index of discrimination, N is the total number of strains in the typing scheme, s is the total number of different strain types, and nj is the number of strains belonging to the jth type.

Data analysis

Spoligotypes in binary format were entered in an Excel spreadsheet and compared with the spoligotyping database SpolDB4. The patterns were established clusters in the BioNumerics software version 5.0 (Applied Maths, Sint-Martens-Latem, Belgium). The strength of association between genotypes of MTB strains and drug resistance was estimated by odds ratio (OR) and 95% confidence interval (95% CI). All tests of significance were two sided and a significant threshold was set at 0.05. Statistical analyses were carried out by using SPSS software 13.0 (SPSS Inc., USA).

Ethical consideration

This project has been approved by Institutional Review Board of Nanjing Medical University. Written informed consent was obtained from all participants. Ethics has been respected throughout the whole study period.

Results

Spoligotyping analysis

Reproducible results were obtained for 260 strains and 34 distinct spoligotype clusters were obtained as referring to the SpolDB4.0 database. Twenty-seven spoligotypes represented the single isolate and the other 233 isolates were grouped into 7 clusters containing from 2 to 199 isolates (Figure 1). Among them, 14 (5.4%) strains were found to be undefined (absent in the SpolDB4.0 database) while 246(94.6%) were successfully clustered by spoligotyping and divided into 11 types (STs). Out of these, 5 STs had 2 or more strains whereas 6 STs had single isolate each. A determination of predominant STs were Beijing family (n = 209, 80.4%), T1 (n = 23, 8.8%) and T2 (n = 6, 2.3%), representing about 91.5% of total strains (Table 1).

Figure 1
**Samples with MIRU-VNTR and spoligotyping**. Clustering was based upon an average of MIRU-VNTR and spoligotyping, clustered using the dice coefficient and the unweighted pair group method with arithmetic averages (UPGMA) in BioNumerics 5.0.

Table 1

Spoligotypes shared by Mycobacterium tuberculosis strains evaluated in this study

S.NO.	Family	N(%)	Binary
1	Beijing	206(79.2)	□□□□□□□□□□□□□□□□□□□□□□□□□□□□□□□□□□■■■■■■■■■
2	Beijing-like	3(1.2)	□□□□□□□□□□□□□□□□□□□□□□□□□□□□□□□□□□■■■□□□■■■
3	T1	23(8.8)	■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■□□□□■■■■■■■
4	T2	6 (2.3)	■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■□□□□■■■□■■■
5	MANU2	2(0.8)	■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■□□■■■■■■■■■
6	U	1(0.4)	■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■
7	T1-T2	1(0.4)	■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■□□□□■■■□□■■
8	H3	1(0.4)	■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■□■□□□□■■■■■■■
9	H4	1(0.4)	■■■■■■■■■■■■■■■■■■■■■■■■■■■■□□□■□□□□■■■■■■■
10	LAM10_CAM	1(0.4)	■■■■■■■■■■■■■■■■■■■■■■□□□■■■■■■■□□□□■■■□■■■
11	CAS1_DELHI	1(0.4)	■■■□□□□■■■■■□■■■■■■■■■□□□□□□□□□□□□■■■■■■■■■
12	Undefined	1(0.4)	■□■■■■■■■■■■□■■■■■■■■■■■■■■■■■■■□□□□■■■■■■■
13	Undefined	1(0.4)	■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■□□□■■■□■■■
14	Undefined	1(0.4)	■■■■□■■■■■■■□■■■■■■■■■■■■■■■■■■■□□□□■■■■■■■
15	Undefined	1(0.4)	□■■■■■■■■■■■□■■■■■■■■■■■■■■■■■■■□□□□■■■□■■■
16	Undefined	1(0.4)	■□■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■□□□□□□■■■■■
17	Undefined	1(0.4)	■■■■□■■■■■□■■□■■■■■■■■■■■■■■■■■■□□■■■■■■■■■
18	Undefined	1(0.4)	■■■■□■■■■■□■■□■■■■■■■■■■■■■■■■■■□□□□■■■■■■■
19	Undefined	1(0.4)	□□□□□■■■■■■■■■■■■■■■■■■■■■■■■■■■□□□□■■■■■■■
20	Undefined	1(0.4)	□□□■■■■■■■■■□■■■■■■■■■■■■■■■■■■■□□□□■■■■■□□
21	Undefined	2(0.8)	■■□□□■□□■■■■■■■■■■■■■■■■■■■■■■■■□□□□■■■■■■■
22	Undefined	1(0.4)	■■■■■■■■■■■■■■■■■■■■■■■■■■■■□□□□□□□□□□□□□■■
23	Undefined	1(0.4)	■■■■■■■■■■■■■■■■■■■□□□□□■□■□■■■■□□□□■■■■■■■
24	Undefined	1(0.4)	■■■■■■■■■■■■■□□□□■□□□□□□■■■■■■■■□□□□■■■■■■■

(■) Presence of spacer; (□) absence of spacer. N = number of strains.

MIRU-VNTR

Among 260 MTB strains being genotyped, 224 different VNTR genotypes were detected. Two hundred and four (91.1%) were unique (i.e., observed for only one strain) and 56 strains could be grouped into 20 clusters, each including 2 to 8 strains (Additional file 1). Two main clusters which contained 8 (14.3%) and 7 (12.5%) strains showed 4 2 4 3 4 3 2 3 3 2 5 1 5 3 3 3 4 5 4 1 4 5 6 3 and 4 2 4 3 4 3 2 3 3 2 5 1 5 3 3 3 4 5 4 1 4 6 6 3 VNTR profiles, respectively (Figure 1). HGDI scores which were calculated for particular MIRU loci varied significantly from 0.717 of Qub11b to 0 of MIRU24. As suggested by Sola et al [22], the MIRU loci were further classified into highly (> 0.6), moderately (0.3 to 0.6), and poorly (< 0.3) discriminating based on the HGDI scores. Five loci (ETRB, ETRC, MIRU20, MIRU02, MIRU24) showed negligible diversity (HGDI < 0.1) whereas ETRE and Qub11b had highly discriminating ability, especially for Beijing family (Table 2).

Table 2

Variability of 24-VNTR loci among Mycobacterium tuberculosis strains (n = 260)

Alias	HGDI (All strains)	HGDI (Beijing family)	12-loci VNTR	NO. of alleles observed	Size range observed
Qub11b	0.717	0.629	√	8	1-8
ETRE	0.695	0.668	√	5	2-6
MIRU26	0.683	0.560	√	7	1-7
Mtub21	0.660	0.535	√	8	1-9
Qub26	0.646	0.613	√	9	1-9
ETRD	0.592	0.536	√	6	1-6
Mtub04	0.467	0.426	√	4	2-5
MIRU40	0.454	0.276	√	4	1-4
MIRU10	0.446	0.262	√	5	1-5
MIRU39	0.390	0.178	√	4	1-4
Mtub39	0.386	0.213	√	7	2-8
Mtub30	0.365	0.196	√	4	2-5
MIRU16	0.348	0.262		4	1-4
ETRF	0.343	0.243		4	1-4
Mtub38	0.335	0.065		4	1-4
ETRA	0.313	0.201		4	2-5
MIRU23	0.268	0.250		5	3-7
MIRU27	0.243	0.084		4	1-4
Qub4156c	0.179	0.203		5	1-5
ETRB	0.096	0.056		3	1-3
ETRC	0.090	0.066		5	2-6
MIRU20	0.023	0.010		2	1-2
MIRU02	0.008	0.010		2	1-2
MIRU24	0	0		1	1

Drug-susceptibility patterns of the MTB isolates

A total of 260 MTB strains isolated from the sputum samples of TB patients were used for drug susceptibility test. Among them, 13.1% (34/260) were MDR-TB strains, which were resistant to at least INH and RIF, the two most powerful antituberculosis drugs. Furthermore, we compared the distribution of drug resistance between Beijing and non-Beijing genotyping strains. As shown in table 3, among non-Beijing family strains, the proportion of mono-drug resistance was 7.8% (n = 4) to INH, 2.0% (n = 1) to RIF, 9.8% (n = 5) to SM, while 19.6% (n = 10) were resistant to at least one drug and 2.0% (n = 1) were MDR-TB. For Beijing family strains, the proportion of mono-drug resistance was 3.3% (n = 7) to INH, 1.4% (n = 3) to RIF, 5.7% (n = 12) to SM and 0.5% (n = 1) to EMB, respectively, while 11.0% (n = 23) were resistant to at least one drug and 15.8% (n = 33) were MDR-TB. A binary logistic regression model was applied to analyze the factors associated with MDR-TB. The Beijing family strains presented increased risks for developing MDR (15.8% vs.2.0%, OR: 11.07, 95%CI: 1.45-84.50) (Table 4).

Table 3

Differences of characteristics between Beijing and non-Beijing family

Factors	Cases n = 260	Non-Beijing family n = 51, n(%)	Beijing family n = 209, n(%)	χ²	P
Sex
Men	180(69.2)	37(72.5)	143(68.4)
Women	80(30.8)	14(27.5)	66(31.6)	0.328	0.567
Treatment history
No	216(83.1)	41(80.4)	175(83.7)
Yes	44(16.9)	10(19.6)	34(16.3)	0.325	0.568
Pan-susceptible cases	186(71.5)	39(76.5)	147(70.3)	0.758	0.384
Any drug-resistant cases	74(28.5)	12(23.5)	62(29.7)	0.758	0.384
Mono-resistance
INH	11(4.2)	4(7.8)	7(3.3)		0.234*
RIF	4(1.5)	1(2.0)	3(1.4)		0.585*
SM	17(6.5)	5(9.8)	12(5.7)		0.340*
EMB	1(0.4)	0	1(0.5)		1*
Total	33(12.7)	10(19.6)	23(11.0)	2.738	0.098
MDR-TB ^#
INH+RIF	6(2.3)	0	6(2.9)		0.601*
INH+RIF+SM	9(3.5)	0	9(4.3)		0.212*
INH+RIF+EMB	2(0.8)	0	2(1.0)		1*
INH+RIF+SM+EMB	17(6.5)	1(2.0)	16(7.7)		0.208*
MDR	34(13.1)	1(2.0)	33(15.8)	6.897	0.009

INH, Isoniazid; RIF, Rifampicin; EMB, Ethambutol; SM, Streptomycin.

#Defined as resistant to at least INH and RIF.

*Value by Fisher's exact test.

Table 4

Factors associated with multi-drug resistant tuberculosis

Factors	Cases	No-MDR	MDR	cOR(95% CI)	P	aOR(95% CI)	P
	n = 260	n = 226, n(%)	n = 34, n(%)
Gender
Men	180	155(68.6)	25(73.5)	1		1
Women	80	71(31.4)	9(26.5)	0.79(0.35-1.77)	0.561	0.81(0.34-1.96)	0.643
Age(years)
< 30	66	56(24.8)	10(29.4)	1		1
30-	111	93(41.2)	18(52.9)	1.08(0.47-2.51)	0.851	0.68(0.23-1.99)	0.479
60-	57	51(22.6)	6(17.6)	0.66(0.22-1.94)	0.449	0.36(0.08-1.71)	0.199
75+	26	26(11.5)	0	0	0.998	0	0.998
Treatment history ^#
No	216	194(85.8)	22(14.2)	1		1
Yes	44	22(64.7)	12(35.3)	3.31(1.49-7.34)	0.002	3.39(1.37-8.41)	0.008
BCG vaccination
No	142	125(55.3)	17(50.0)	1		1
Yes	118	101(44.7)	17(50.0)	1.24(0.60-2.55)	0.562	0.81(0.26-2.48)	0.711
Migrant population
No	174	151(66.8)	23(67.6)	1		1
Yes	86	75(33.2)	11(32.4)	0.96(0.45-2.08)	0.923	0.94(0.40-2.23)	0.894
Beijing family
No	51	50(22.1)	1(2.9)	1		1
Yes	209	176(77.9)	33(97.1)	9.38(1.25-70.26)	0.009	11.07(1.45-84.50)	0.023

cOR, crude odds ratio; aOR, adjusted odds ratio, adjusted for sex, age, treatment history, BCG vaccination, migrant population, and Beijing family; #:Those having previously received anti-tuberculosis treatment for more than one month was defined as previously treated cases.

Discussion

This study demonstrates that the Beijing genotype is the most predominant lineage of TB strains in Jiangsu province, which is consistent to the findings from other areas of China. As described in table 5, Beijing family accounts for 85.12% to 92.59% of the MTB strains currently epidemic in the Beijing area [23, 24] and it is also prevalent in Heilongjiang (89.50%) [10], Gansu (87.50%) [24], Jilin (89.88%) [24], Ningxia (67.12%) [25], Henan (80.00%) [24], Shandong (85.93%) [9], Shanxi (80.00%) [24], Tianjin(91.07%) [26], Tibet (90.38%) [24], Xinjiang (67.80%) [24], Anhui (85.35%) [24], Fujian (54.50%) [24], Hong Kong (68.45% and 70.00%) [27, 28], Hunan (66.00%) [24], Jiangsu and Zhejiang (69.23%) [29], Shanghai (77.14%) [30], Zhejiang (64.84%) [24], Guangxi (55.29%) [24], Taipei (52.53%) [13], and Sichuan (57.89% to 63.64%) [24, 31]. It was found that amongst all MTB strains studied, the "Beijing" genotype strains were highly prevalent in our geographic area. The proportion of Beijing family circulating in the north area (83.91%) is higher than that in the south area (66.40%) of China (P < 0.001).

Table 5

Distribution of Beijing family strains in China

Area	Strains tested	Beijing family strains	Proportion (%)	Location
Beijing [23]	121	103	85.12	north
Beijing[23, 24]	108	100	92.59	north
Gansu [24]	224	196	87.50	north
Heilongjiang [10]	200	179	89.50	north
Henan [24]	95	76	80.00	north
Jilin [24]	326	293	89.88	north
Ningxia [25]	72	49	67.12	north
Shandong [9]	135	116	85.93	north
Shanxi [24]	115	92	80.00	north
Tianjin [26]	112	102	91.07	north
Tibet [24]	208	188	90.38	north
Xinjiang [24]	205	139	67.80	north
Anhui [24]	157	134	85.35	south
Fujian [24]	433	236	54.50	south
Hong Kong [27]	355	243	68.45	south
Hong Kong [28]	500	350	70.00	south
Hunan [24]	100	66	66.00	south
Jiangsu and Zhejiang [29]	351	243	69.23	south
Zhejiang [24]	91	59	64.84	south
Shanghai [30]	175	135	77.14	south
Guangxi [24]	208	115	55.29	south
Taipei [13]	356	187	52.53	south
Sichuan [24]	76	44	57.89	south
Sichuan [31]	143	91	63.64	south
Jiangsu	260	209	80.38	south

Though world widely the Beijing genotype is apparently the second most prevalent genotype, following after AFR lineage [16], reasons for its successful global spread remain poorly understood. The infective success of this lineage seems to be associated with its effect on the immune response, in that it can control the release of the macrophage-derived cytokines that play a central role in directing the immune response towards a non-protective Th2 phenotype [32, 33]. The widely distributed (but not universal) association of drug resistance and the Beijing genotype suggests that these strains may have a particular propensity for acquiring drug resistance [14]. Besides China, they are increasingly reported in other areas of the world and are frequently associated with outbreaks of TB or MDR-TB [14, 29, 34–37], but not all [9, 10, 27, 38, 39]. In the present study, we observed a significant difference of drug resistance between Beijing and non-Beijing strains, where the Beijing family was associated with an increased risk of MDR-TB. These data are accordant with the results from studies so far performed in China [29, 40]. Definite conclusions on the extent of spread and associations with drug resistance, however, could not yet be drawn. This was due to the limited amount of information available from most areas of the world, the possible biases in many of the published reports, and the absence of standard definitions and study designs.

Both spoligotyping and MIRU-VNTR have different abilities to discriminate MTB strains and can be performed in combination. For example, among 260 strains included in this analysis, 209 (80.4%) were classified as the Beijing genotype. The majority of these Beijing genotype isolates showed the typical spoligotype pattern (hybridization to all of spacers 35-43 and no hybridization to spacers 1-34). Four strains showed the non-Beijing genotype spoligotype pattern, but had specific Beijing genotype signatures of MIRU-VNTR. Furthermore, another four strains showing the typical Beijing genotype spoligotype pattern, but can't be distinguished by MIRU-VNTR (Figure 1). One possible explanation might be the multiple infections of strains in these eight samples [41–43], however, it is rarely found when using the DNA fingerprinting method. Till now, no definition of the Beijing genotype on the basis of genetic markers will be 100% perfect, because there will always be exceptional strains [14]. As shown in our study, when 209 strains in the Beijing clade identified by spoligotyping were further discriminated by MIRU-VNTR, nearly all (205/209, 98.1%) Beijing strains were parts of a clade (Figure 1). This study further demonstrated that MTB strains grouped into the Beijing family by spoligotyping had similar patterns of MIRU-VNTR. This was borne out by the fact that all MIRU-VNTR patterns of Beijing family strains were highly similar.

Spoligotyping has been proved to be the most useful method to recognize Beijing lineage strains as it is the most rapid and easiest method to apply and it correlates well with the other methods [14]. But this method had a poor discriminability when it was used independently, thus should be performed in combination with another independent method (for example MIRU-VNTR typing) to allow effective epidemiological investigation [22]. As shown in the present study, spoligotyping alone yielded 34 genotypes, 7 clusters, with 27 unique genotypes and a HGDI of 0.409. MIRU-VNTR clustered 21.5% of strains with an HGDI of 0.998. The combination of spoligotyping and MIRU-VNTR typing in our population clustered 19 of all 260 strains. Different VNTR typing sets showed various efficiencies in different Beijing genotype strains (Table 6). The discriminatory power of the traditional 12 loci MIRU-VNTR has been found to be insufficient [5], which is often failed to differentiate clonal strains in areas where the Beijing family genotype is predominant [27]. Other sets of MIRU-VNTR loci, such as an optimized set of 24 loci, have also been defined. However, all 24 loci are not required for MTB genotyping in any given situation [5], as the number of loci required depends on the lineage known to be prevalent in the investigated area. We have evaluated MIRUs loci, QUB loci and ETR loci individually for their abilities to differentiate the Beijing and non-Beijing genotype families of MTB. When the top 12 loci (MIRU26,40,10, 39; QUB11b,26; Mtub21,04,39,30; ETR E,D) were selected for typing, the set produced a HGDI value of 0.996 in all stains and 0.994 in Beijing family stains, which was almost equal to that obtained in 24 loci MIRU-VNTR (0.998 and 0.997). This optimized 12 loci VNTR typing set could be an important tool for tracking MTB of the Beijing family in this area.

Table 6

Comparative analysis of genotyping methods applied to Mycobacterium tuberculosis strains (n = 260)

Typing method, strain groups	Genotypes in total	Unique genotypes	clusters	Strains in cluster	HGDI value
Spoligotyping
All strains	34	27	7	2-199	0.409
Beijing family only	7	4	3	3-196	0.093
Non-Beijing strains	27	23	4	2-19	0.856
MIRU-VNTR analysis, All strains
24-VNTR analysis	224	204	20	2-8	0.998
12-VNTR analysis	204	177	27	2-9	0.996
Spoligo+24-VNTR
All strains	227	208	19	2-8	0.998
Beijing family only	178	160	18	2-8	0.997

As observed by 24-MIRU-VNTR analysis, it appeared that individuals infected with strains belonging to the Beijing genotype were more likely to be part of a cluster. Beijing strains were isolated from patients from different counties of Jiangsu Province, suggesting that these virulent strains are spreading throughout the province. In addition to Beijing family strains, we also detected strains belonging to other families such as T1, T2, H3, H4, CAS, LAM, U, and MANU2. The CAS family was previously found primarily in India [44], and also detected in Tibet and Xinjiang [24]. The LAM family was predominantly prevalent in South America and West-Africa, and also found in Taiwan of China [13]. It can be suggested that strains of these families have started spreading across Jiangsu province; of course, this requires further testing using more extended typing of more clinical strains.

Future studies are needed to clarify the characteristics of Beijing family, one important genotype of MTB. Possible variations in host-pathogen interactions among strains of various genotypes needs to be identified and the role of Beijing genotype infection as a possible risk factor for drug resistance and/or treatment failure must urgently be addressed in longitudinal studies in the affected high incident regions. In short, the effect of this genotype on TB control efforts need to be further investigated.

Conclusion

This is the first report on the genotypes of MTB stains using spoligotyping in combination with 24-loci MIRU-VNTR technology in Jiangsu province of China. Based on our preliminary data, Beijing family strains have predominantly prevalent in Jiangsu and shown a high number of clusters in the study population. Furthermore, we also observed that the high prevalence of Beijing genotype may be associated with the risk of MDR-TB. Further genome level studies will be needed to investigate the molecular specificities of the major subgroup.

Notes

Declarations

Acknowledgements and funding

This study is supported by "Transmission Mode of Tuberculosis Project of the National Key Programme of Mega Infectious" (2008ZX10003-010) and "National S&T Major Project Foundation of China" (2009ZX10004-904).

Authors’ Affiliations

(1)

Department of Epidemiology and Biostatistics, School of Public Health, Nanjing Medical University, Nanjing, PR, China

(2)

Department of Chronic Communicable Disease, Center for Disease Control and Prevention of Jiangsu Province, Nanjing, PR, China

(3)

National Institute for Communicable Disease Control and Prevention, Chinese Center for Disease Control and Prevention, State Key Laboratory for Infectious Disease Prevention and Control, Beijing, PR, China

References

Global tuberculosis control - surveillance,planning,financing.WHO Report. 2010Google Scholar
Toungoussova OS, Bjune G, Caugant DA: Epidemic of tuberculosis in the former Soviet Union: social and biological reasons. Tuberculosis (Edinb). 2006, 86 (1): 1-10. 10.1016/j.tube.2005.04.001.View ArticleGoogle Scholar
Christianson S, Wolfe J, Orr P, Karlowsky J, Levett PN, Horsman GB, Thibert L, Tang P, Sharma MK: Evaluation of 24 locus MIRU-VNTR genotyping of Mycobacterium tuberculosis isolates in Canada. Tuberculosis (Edinb). 2010, 90 (1): 31-38. 10.1016/j.tube.2009.12.003.View ArticleGoogle Scholar
Kamerbeek J, Schouls L, Kolk A, van Agterveld M, van Soolingen D, Kuijper S, Bunschoten A, Molhuizen H, Shaw R, Goyal M, et al: Simultaneous detection and strain differentiation of Mycobacterium tuberculosis for diagnosis and epidemiology. Journal of Clinical Microbiology. 1997, 35 (4): 907-914.PubMedPubMed 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. Journal of Clinical Microbiology. 2006, 44 (12): 4498-4510. 10.1128/JCM.01392-06.View ArticlePubMedPubMed CentralGoogle Scholar
Mokrousov I, Narvskaya O, Limeschenko E, Vyazovaya A, Otten T, Vyshnevskiy B: Analysis of the allelic diversity of the mycobacterial interspersed repetitive units in Mycobacterium tuberculosis strains of the Beijing family: practical implications and evolutionary considerations. Journal of Clinical Microbiology. 2004, 42 (6): 2438-2444. 10.1128/JCM.42.6.2438-2444.2004.View ArticlePubMedPubMed CentralGoogle 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. Journal of Clinical Microbiology. 1995, 33 (12): 3234-3238.PubMedPubMed CentralGoogle Scholar
Hu Y, Hoffner S, Jiang W, Wang W, Xu B: Extensive transmission of isoniazid resistant M. tuberculosis and its association with increased multidrug-resistant TB in two rural counties of eastern China: a molecular epidemiological study. BMC Infect Dis. 2010, 10: 43-10.1186/1471-2334-10-43.View ArticlePubMedPubMed CentralGoogle Scholar
Ma X, Wang H, Deng Y, Liu Z, Xu Y, Pan X, Musser JM, Graviss EA: rpoB Gene mutations and molecular characterization of rifampin-resistant Mycobacterium tuberculosis isolates from Shandong Province, China. Journal of Clinical Microbiology. 2006, 44 (9): 3409-3412. 10.1128/JCM.00515-06.View ArticlePubMedPubMed CentralGoogle Scholar
Wang J, Liu Y, Zhang CL, Ji BY, Zhang LZ, Shao YZ, Jiang SL, Suzuki Y, Nakajima C, Fan CL, et al: Genotypes and characteristics of clustering and drug-susceptibility of Mycobacterium tuberculosis isolates in Heilongjiang Province, China. Journal of Clinical Microbiology. 2011Google Scholar
Shao Y, Yang D, Xu W, Lu W, Song H, Dai Y, Shen H, Wang J: Epidemiology of anti-tuberculosis drug resistance in a chinese population: current situation and challenges ahead. BMC Public Health. 2011, 11: 110-10.1186/1471-2458-11-110.View ArticlePubMedPubMed CentralGoogle Scholar
Guidelines for surveillance of drug resistance in tuberculosis (WHO/HTM/TB/2009.422). [http://www.who.int/tb/publications/mdr_surveillance/en/index.html]
Dou HY, Tseng FC, Lin CW, Chang JR, Sun JR, Tsai WS, Lee SY, Su IJ, Lu JJ: Molecular epidemiology and evolutionary genetics of Mycobacterium tuberculosis in Taipei. BMC Infect Dis. 2008, 8: 170-10.1186/1471-2334-8-170.View ArticlePubMedPubMed CentralGoogle Scholar
Kremer K, Glynn JR, Lillebaek T, Niemann S, Kurepina NE, Kreiswirth BN, Bifani PJ, van Soolingen D: Definition of the Beijing/W lineage of Mycobacterium tuberculosis on the basis of genetic markers. Journal of Clinical Microbiology. 2004, 42 (9): 4040-4049. 10.1128/JCM.42.9.4040-4049.2004.View ArticlePubMedPubMed CentralGoogle Scholar
Gori A, Bandera A, Marchetti G, Degli Esposti A, Catozzi L, Nardi GP, Gazzola L, Ferrario G, van Embden JD, van Soolingen D, et al: Spoligotyping and Mycobacterium tuberculosis. Emerging Infectious Diseases. 2005, 11 (8): 1242-1248.View ArticlePubMedPubMed 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
Le Fleche P, Fabre M, Denoeud F, Koeck JL, Vergnaud G: High resolution, on-line identification of strains from the Mycobacterium tuberculosis complex based on tandem repeat typing. BMC Microbiol. 2002, 2: 37-10.1186/1471-2180-2-37.View ArticlePubMedPubMed CentralGoogle Scholar
Meeker-O'Connell RFaWA: Genetic diversity in the Mycobacterium tuberculosis complex based on variable numbers of tandem DNA repeats. Microbiology. 1998, 144: 1189-1196. 10.1099/00221287-144-5-1189.View ArticlePubMedGoogle Scholar
Skuce RA, McCorry TP, McCarroll JF, Roring SM, Scott AN, Brittain D, Hughes SL, Hewinson RG, Neill SD: Discrimination of Mycobacterium tuberculosis complex bacteria using novel VNTR-PCR targets. Microbiology. 2002, 148 (Pt 2): 519-528.View ArticlePubMedGoogle Scholar
Supply P, Lesjean S, Savine E, Kremer K, van Soolingen D, Locht C: Automated high-throughput genotyping for study of global epidemiology of Mycobacterium tuberculosis based on mycobacterial interspersed repetitive units. Journal of Clinical Microbiology. 2001, 39 (10): 3563-3571. 10.1128/JCM.39.10.3563-3571.2001.View ArticlePubMedPubMed CentralGoogle Scholar
Hunter PR: Reproducibility and indices of discriminatory power of microbial typing methods. Journal of Clinical Microbiology. 1990, 28 (9): 1903-1905.PubMedPubMed CentralGoogle Scholar
Sola C, Ferdinand S, Mammina C, Nastasi A, Rastogi N: Genetic diversity of Mycobacterium tuberculosis in Sicily based on spoligotyping and variable number of tandem DNA repeats and comparison with a spoligotyping database for population-based analysis. Journal of Clinical Microbiology. 2001, 39 (4): 1559-1565. 10.1128/JCM.39.4.1559-1565.2001.View ArticlePubMedPubMed CentralGoogle Scholar
Jiao WW, Mokrousov I, Sun GZ, Guo YJ, Vyazovaya A, Narvskaya O, Shen AD: Evaluation of new variable-number tandem-repeat systems for typing Mycobacterium tuberculosis with Beijing genotype isolates from Beijing, China. Journal of Clinical Microbiology. 2008, 46 (3): 1045-1049. 10.1128/JCM.01869-07.View ArticlePubMedPubMed CentralGoogle Scholar
Dong H, Liu Z, Lv B, Zhang Y, Liu J, Zhao X, Wan K: Spoligotypes of Mycobacterium tuberculosis from different Provinces of China. Journal of Clinical Microbiology. 2010, 48 (11): 4102-4106. 10.1128/JCM.00549-10.View ArticlePubMedPubMed CentralGoogle Scholar
Li WM, W S, Pei XY, Liu ZQ, Zhong Q, Qian M, Zhao B, Duan Mu HJ: Mycobacterium tuberculosis strains from Beijing, Guangdong and Ningxia. Chin J Epidemiol. 2003, 24 (5): 381-384.Google Scholar
Xie T, JH , Zhao DF, Mu C, Zhao H, Wu ZP, Wang XX: Polymorphism of Variable-Number Tandem Repeat in Mycobacterium Tuberculosis Strains Isolated from Tianjin. Tianjin Med J. 2010, 38: 940-943.Google Scholar
Kam KM, Yip CW, Tse LW, Wong KL, Lam TK, Kremer K, Au BK, van Soolingen D: Utility of mycobacterial interspersed repetitive unit typing for differentiating multidrug-resistant Mycobacterium tuberculosis isolates of the Beijing family. Journal of Clinical Microbiology. 2005, 43 (1): 306-313. 10.1128/JCM.43.1.306-313.2005.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.View ArticlePubMedPubMed CentralGoogle Scholar
Hu Y, Ma X, Graviss EA, Wang W, Jiang W, Xu B: A major subgroup of Beijing family Mycobacterium tuberculosis is associated with multidrug resistance and increased transmissibility. Epidemiol Infect. 2011, 139 (1): 130-138. 10.1017/S0950268810000890.View ArticlePubMedGoogle Scholar
Qiao KWH, Yang CG, Luo T, Mei J, Gao Q: The application of variable number of tandem repeats in the microevolution study of Mycobacterium tuberculosis strain of Beijing genotype in Chongming Island, Shanghai. Journal of Microbes and Infections. 2010, 5 (4): 208-213.Google Scholar
Peng ZZC, X LL, L RX: Correlation between Beijing genotype strains of Mycobacterium tuberculosis and drug resistant phenotypes in Chongqing children with tuberculosis. Journal of Third Military Medical University. 2010, 32 (23): 2525-2528.Google Scholar
Reed MB, Domenech P, Manca C, Su H, Barczak AK, Kreiswirth BN, Kaplan G, Barry CE: A glycolipid of hypervirulent tuberculosis strains that inhibits the innate immune response. Nature. 2004, 431 (7004): 84-87. 10.1038/nature02837.View ArticlePubMedGoogle Scholar
Manca C, Reed MB, Freeman S, Mathema B, Kreiswirth B, Barry CE, Kaplan G: Differential monocyte activation underlies strain-specific Mycobacterium tuberculosis pathogenesis. Infection and Immunity. 2004, 72 (9): 5511-5514. 10.1128/IAI.72.9.5511-5514.2004.View ArticlePubMedPubMed CentralGoogle Scholar
Purwar S, Chaudhari S, Katoch VM, Sampath A, Sharma P, Upadhyay P, Chauhan DS: Determination of drug susceptibility patterns and genotypes of Mycobacterium tuberculosis isolates from Kanpur district, North India. Infect Genet Evol. 2011, 11 (2): 469-475. 10.1016/j.meegid.2010.12.010.View ArticlePubMedGoogle Scholar
Caws M, Thwaites G, Stepniewska K, Nguyen TN, Nguyen TH, Nguyen TP, Mai NT, Phan MD, Tran HL, Tran TH, et al: Beijing genotype of Mycobacterium tuberculosis is significantly associated with human immunodeficiency virus infection and multidrug resistance in cases of tuberculous meningitis. Journal of Clinical Microbiology. 2006, 44 (11): 3934-3939. 10.1128/JCM.01181-06.View ArticlePubMedPubMed CentralGoogle Scholar
Al Hajoj S, Rastogi N: The emergence of Beijing genotype of Mycobacterium tuberculosis in the Kingdom of Saudi Arabia. Ann Thorac Med. 2010, 5 (3): 149-152. 10.4103/1817-1737.65045.View ArticlePubMedPubMed CentralGoogle Scholar
Brown T, Nikolayevskyy V, Velji P, Drobniewski F: Associations between Mycobacterium tuberculosis strains and phenotypes. Emerging Infectious Diseases. 2010, 16 (2): 272-280.View ArticlePubMedPubMed CentralGoogle Scholar
Ani A, Bruvik T, Okoh Y, Agaba P, Agbaji O, Idoko J, Dahle UR: Genetic diversity of Mycobacterium tuberculosis Complex in Jos, Nigeria. BMC Infect Dis. 2010, 10: 189-10.1186/1471-2334-10-189.View ArticlePubMedPubMed CentralGoogle Scholar
Lai CC, Tan CK, Lin SH, Liao CH, Huang YT, Chou CH, Hsu HL, Wang CY, Lin HI, Hsueh PR: Clinical and genotypic characteristics of extensively drug-resistant and multidrug-resistant tuberculosis. European journal of clinical microbiology & infectious diseases. official publication of the European Society of Clinical Microbiology. 2010, 29 (5): 597-600. 10.1007/s10096-010-0874-6.View ArticleGoogle Scholar
Mokrousov I, Jiao WW, Sun GZ, Liu JW, Valcheva V, Li M, Narvskaya O, Shen AD: Evolution of drug resistance in different sublineages of Mycobacterium tuberculosis Beijing genotype. Antimicrob Agents Chemother. 2006, 50 (8): 2820-2823. 10.1128/AAC.00324-06.View ArticlePubMedPubMed CentralGoogle Scholar
Warren RM, Victor TC, Streicher EM, Richardson M, Beyers N, Gey van Pittius NC, van Helden PD: Patients with active tuberculosis often have different strains in the same sputum specimen. American journal of respiratory and critical care medicine. 2004, 169 (5): 610-614. 10.1164/rccm.200305-714OC.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. Respiratory Research. 2005, 6: 134-10.1186/1465-9921-6-134.View ArticlePubMedPubMed CentralGoogle Scholar
Mallard K, McNerney R, Crampin AC, Houben R, Ndlovu R, Munthali L, Warren RM, French N, Glynn JR: Molecular detection of mixed infections of Mycobacterium tuberculosis strains in sputum samples from patients in Karonga District, Malawi. Journal of Clinical Microbiology. 2010, 48 (12): 4512-4518. 10.1128/JCM.01683-10.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. Emerging Infectious Diseases. 2004, 10 (6): 1138-1142.View ArticlePubMedPubMed CentralGoogle Scholar

Pre-publication history

The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2334/11/288/prepub

Copyright

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.

Molecular typing of mycobacterium tuberculosis isolates circulating in Jiangsu Province, China

Background

Methods

Results

Conclusions

Mycobacterial strains

Strains isolation and drug susceptibility test

Genomic DNA extraction

Genotyping

Data analysis

Ethical consideration

Spoligotyping analysis

MIRU-VNTR

Drug-susceptibility patterns of the MTB isolates

Acknowledgements and funding

Pre-publication history

BMC Infectious Diseases

Contact us