Next Article in Journal
Chronic Manganese Toxicity Associated with Voltage-Gated Potassium Channel Complex Antibodies in a Relapsing Neuropsychiatric Disorder
Next Article in Special Issue
A Simulation-Based Study on the Comparison of Statistical and Time Series Forecasting Methods for Early Detection of Infectious Disease Outbreaks
Previous Article in Journal
Mapping European Welfare Models: State of the Art of Strategies for Professional Integration and Reintegration of Persons with Chronic Diseases
Previous Article in Special Issue
The Role of Culex pipiens L. (Diptera: Culicidae) in Virus Transmission in Europe
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJERPH
Volume 15
Issue 4
10.3390/ijerph15040782
ijerph-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Epidemiological Significance and Temporal Stability of Mycobacterial Interspersed Repetitive Units-Variable Number of Tandem Repeats-Based Method Applied to Mycobacterium tuberculosis in China

by
Yang Li
1,
Yi Hu
1,2,*,
Mikael Mansjö
3,
Qi Zhao
1,
Weili Jiang
1,
Solomon Ghebremichael
3,
Sven Hoffner
2 and
Biao Xu
1,2
1
Department of Epidemiology, China and Key Laboratory of Public Health Safety (Fudan University), School of Public Health, Fudan University, Ministry of Education, Shanghai 200032, China
2
Department of Public Health Sciences, Karolinska Institutet, SE-171 77 Stockholm, Sweden
3
The Public Health Agency of Sweden, SE-171 82 Solna, Sweden
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2018, 15(4), 782; https://doi.org/10.3390/ijerph15040782
Submission received: 8 March 2018 / Revised: 11 April 2018 / Accepted: 11 April 2018 / Published: 17 April 2018
(This article belongs to the Special Issue Emerging Infectious Disease (EID) Research, Management and Response)
Browse Figures
Versions Notes

Abstract

:
This study aimed to validate the epidemiological significance and temporal stability of Mycobacterial Interspersed Repetitive Units-Variable Number of Tandem Repeats (MIRU-VNTR) typing in a genetically and geographically diverse set of clinical isolates from patients diagnosed with pulmonary tuberculosis in China. Between 2010 and 2013, a total of 982 Mycobacterium tuberculosis isolates were collected from four population-based investigations in China. Apart from the currently applied 24-locus MIRU-VNTR, six additional hypervariable loci were analyzed in order to validate the MIRU-VNTR combinations in terms of their epidemiological links, clustering time span, and paired geographic distance. In vitro temporal stability was analyzed for both individual MIRU-VNTR loci, and for several combinations of loci. In the present study, four MIRU-VNTR combinations, including the hypervariable loci 3820, 3232, 2163a, and 4120, were evaluated. All of these combinations obtained a Hunter-Gaston discriminatory index (HGDI) value over 0.9900 with a reduced clustering proportion (from 32.0% to 25.6%). By comparing epidemiological links, clustering time span, and paired geographic distance, we found that the performances of the four MIRU-VNTR combinations were comparable to the insertion sequence 6110 restriction fragment length polymorphism (IS6110-RFLP), and significantly better than that of 24-locus MIRU-VNTR genotyping alone. The proportion of temporally stable loci ranged from 90.5% to 92.5% within the combined MIRU-VNTR genotyping, which is higher than IS6110-RFLP (85.4%). By adding four hypervariable loci to the standard 24-locus MIRU-VNTR genotyping, we obtained a high discriminatory power, stability and epidemiological significance. This algorithm could therefore be used to improve tuberculosis transmission surveillance and outbreak investigation in China.

1. Introduction

Tuberculosis (TB) remains an important threat to public health worldwide. There has been a steady decline in both the incidence and prevalence of TB in China from 2010 to 2016, however TB remains a leading infectious disease [1]. Furthermore, drug-resistant TB, especially multidrug-resistant (MDR)-TB, has made the control of TB in China even more difficult [2]. The world health organization (WHO) estimated that China’s incidence of TB in 2016 was 64 cases per 100,000 people, of which 8.2% were MDR-TB/rifampicin-resistant (RR) TB cases [3].
The molecular epidemiological investigation of Mycobacterium tuberculosis (M. tuberculosis) has significantly increased our understanding of TB epidemiology, and contributed to population-based control of tuberculosis. Molecular genotyping suggests the potential transmission link between TB patients, enabling the identification of outbreaks, and laboratorial cross contamination [4]. Sequencing, e.g., whole genome sequencing, is very promising, but cannot easily be applied in high-burden TB settings, and there is still a demand for alternative molecular typing techniques. The earlier canonical reference method, the insertion sequence 6110 restriction fragment length polymorphism (IS6110-RFLP), cannot be universally applied due to difficulties in comparing results across laboratories, and handling of the strains with either too low [5,6] or too high numbers of IS6110 copies [7]. In addition, the Beijing family strain has been predominant in China, accounting for 63.97% of the strain according to the national TB survey in 2007. Spoligotyping can be used to identify, but not to subdivide, the Beijing family, such as modern Beijing family strains, which limits its usefulness in the Chinese setting. For these reasons, the Polymerase Chain Reaction (PCR) based Mycobacterial Interspersed Repetitive Unit-Variable Number of Tandem Repeats (MIRU-VNTR) method is preferable, given its turnaround time and the convenience of its implementation [8,9,10]. Furthermore, it has a discriminatory power comparable to that of IS6110-RFLP [11].
For epidemiological applications, a typing test must have the ability to quickly and reliably differentiate between related bacterial isolates [12]. Moreover, the genetic stability will determine the applicability of the technology for different epidemiological purposes. In previous studies, several additional MIRU-VNTR hypervariable loci have been applied together with the traditional 24-locus MIRU-VNTR to achieve a higher resolution [13,14,15]. However, the epidemiological significance and temporal stability of this combined MIRU-VNTR genotyping deserves further evaluation due to variations in different populations and settings. For these reasons, we aimed to validate the epidemiological significance and temporal stability of the combined MIRU-VNTR typing in a genetically and geographically diverse set of clinical isolates from TB patients registered between 2010 and 2013 in China.

2. Materials and Methods

2.1. Study Population

Eight counties from four diverse geographical areas of China were chosen as study sites, including two from the Zhejiang province (East), two from the Liaoning Province (North), two from the Sichuan Province (Midwest) and two from the Guangdong Province (South). All of the pulmonary TB patients registered at the TB designated health facilities between 2010 and 2013 were enrolled. Their socio-demographic and clinical information was collected from medical records and the national TB registry system. Meanwhile, sputum samples were collected for culture and microbiological testing. Study subjects were defined as TB patients with positive culture results, whose M. tuberculosis isolates were further grown on Löwenstein–Jensen (LJ) medium for up to 4 weeks. The DNA was extracted according to the CTAB (Cetyltrimethylammonium Ammonium Bromide) method [16,17]. Extracted DNA was then sent to the TB molecular laboratory of Fudan University’s School of Public Health in Shanghai for molecular genotyping.
The study protocol and informed consents were approved by the Institutional Review Board of the School of Public Health, Fudan University (IRB00002408). All patients included in this study provided written informed consent for the use of clinical samples.

2.2. MIRU-VNTR Loci Selection and PCR Amplification

Briefly, a set of 15 discriminatory MIRU-VNTR loci (424, 577, 580, 802, 960, 1644, 1955, 2163b, 2165, 2401, 2996, 3192, 3690, 4052, 4156), 9 ancillary loci (154, 2059, 2347, 2461, 2531, 2687, 3007, 3171, 4348) and 6 hypervariable loci (2163a, 3232, 4120, 1982, 3336, 3820) were analyzed in all studied isolates. PCR was performed as previously reported [18]. The PCR products were identified by electrophoresis in 1.5% agarose gel using a 100-bp and a 50-bp DNA ladder (TAKARA). The size of the PCR fragments was defined using the Quantity One software, version 4.6.2 (BIO-RAD Laboratories, Hercules, CA, USA). DNA from the M. tuberculosis reference strain H37Rv was used as a positive control to assure the quality of amplicon sizing. Meanwhile, sterile water was used as a negative control. The exact number of complete repeats present was assigned using a derived allele-repeat table corresponding to the size of the PCR product. A MIRU-VNTR cluster was defined as two or more samples having the identical number of repeats at each of the MIRU-VNTR loci.

2.3. IS6110-RFLP Genotyping

DNA fingerprinting was conducted as earlier described [19]. In short, 32P-labeled 245 bp IS6110 fragment was applied as a probe. H37Rv strains and sterile water were used as the positive and negative controls, respectively. After hybridization with the DNA probe, IS6110-RFLP was identified by autoradiography. The RFLP type was assigned based on the IS6110 pattern (number and position of bands in the autoradiography). A cluster was defined as M. tuberculosis isolates harboring the identical IS6110RFLP pattern.

2.4. Evaluation of Temporal Stability

A total of 294 clinical isolates were randomly selected in order to evaluate the temporal stability of all 30 of the studied MIRU-VNTR loci. The investigation of clonal stability was performed by isolating all the colonies in a clone 8 times on LJ media at intervals of 2 to 3 months over 2 years. The copy number of the studied MIRU-VNTR loci, as well as their temporal stability, were determined for every subculture from each isolate.

2.5. Data Collection and Measurement

The subjects with clustered isolates were interviewed in terms of socio-demographic characteristics, clinical manifestation, predisposing risk factors, and evidence of contact with patients with active TB. The collected information also included occupational, social, and recreational activities.
The epidemiological links were defined as “Confirmed” in cases of simultaneously staying or working in the same place or by family relationship, while epidemiological links were defined as “Potential” in cases of residing in the same or the neighboring village, albeit in the same township. If the subjects lived in different townships or lacked a common neighborhood, their epidemiological links were regarded as “Not confirmed”. Clustering time span was defined as the period between the diagnosis dates of the initial and final patients in the cluster. Among the patients with an identical genotype, a paired geographic distance was calculated according to the longitude and latitude of their residential addresses.

2.6. Evaluation of Genotyping Methods and Statistical Analysis

The genotyping results were entered and analyzed using the BioNumerics software, version 7.6 (Applied Maths, Kortrijk, Belgium). The clusters were defined using the Jaccard index. The h value was calculated to define the allelic diversity of each locus. The equation was h = 1 − ∑ x j 2 ( n n − 1 ) . The “Hunter-Gaston discrimination index” (HGI) was applied to define the discriminative power of different genotyping strategies for comparison. The equation was: D = 1 − 1 N ( N − 1 ) ∑ j = 1 s X j ( x j − 1 ) [19]. The stability of copies in MIRU-VNTR loci over time was examined using a Generalized Estimating Equation (GEE), which considered the effect caused by relativity between samples, and can be applied for repeated measurements data. The stability of MIRU-VNTR combinations was analyzed using Kaplan-Meier survival analysis. An “event” was defined to occur when the number of repeats at each MIRU-VNTR locus was observed to change. The time at which the “event” occurred was defined to be the midpoint between the two sampling times of the studied isolates. The log rank χ2 test was used to assess the differences in survival curves between methods. The Cochran-Mantel-Haenszel test was applied for categorical variables, while analysis of variance (ANOVA) and a multiplicity test (LSD) were used for continuous variables. Significant level (α) was set at 0.05. Data analysis was run using SAS v9.3 packages (SAS Institute, Cary, NC, USA).

3. Results

3.1. Patient Characteristics

Of the 1065 isolates collected within the study period, 48 non-tuberculosis mycobacterium isolates and three M. bovis isolates were excluded. Of the 1014 M. tuberculosis isolates, 982 had molecular information available for spoligotyping MIRU-VNTR and IS6110-RFLP, and were included in the data analysis.
Patients from the four areas had a similar average age ranging from 45 to 48 years. The male to female ratio was approximately 7:3. A total of 252 (25.7%) patients had cavities on chest radiographs. Nearly sixty percent of them had smear positive TB. In addition, the socio-demographic and clinical characteristics were similar in the four areas (Table 1). By means of spoligotyping, isolates from the Beijing family accounted for 69.0% (n = 678).

3.2. Genetic Diversity of Each MIRU-VNTR Locus

The number of repeats of each MIRU-VNTR locus was shown in Table S1. The genetic heterogeneity of each MIRU-VNTR locus was estimated by calculating h values (Table S2) for the whole sample set, as well as for the for Beijing or non-Beijing families. The h value differed from 0.836 for locus 3820 to 0.020 for locus 3171, with the genetic diversity being similar in the four areas. The loci 2059, 2347, 2687, and 3171 were highly homogeneous within the Beijing family of M. tuberculosis (h < 0.1). Of the six hypervariable loci, a high allelic diversity (h > 0.7) was seen in loci 3820, 3232, 2163a, and 4120.

3.3. Establishment of MIRU-VNTR Combinations

To identify a MIRU-VNTR combination suitable for the Chinese setting, we combined the 24-locus MIRU-VNTR genotyping with the four hypervariable loci (3820, 3232, 2163a, and 4120). As shown in Table 2, the combination of the 24 MIRU-VNTR loci with four additional loci had a higher discriminatory power than the 24-locus MIRU-VNTR genotyping alone (HGDI > 0.9900). Additionally, 123 clusters were identified by 24-locus MIRU-VNTR genotyping including 405 isolates (Table 3). When comparing the four combined MIRU-VNTR genotyping, the number of clusters was reduced from 98 to 82, as was the clustering proportion (from 32.0% to 25.6%).

3.4. Evaluation of Different Genotyping Strategies.

We further assessed the epidemiological links, clustering time span, and paired geographic distance of the MIRU-VNTR combinations compared to 24-locus MIRU-VNTR genotyping and IS6110-RFLP. As shown in Table 4, it was observed that the MIRU-VNTR combination including four hypervariable loci identified more proven epidemiological links (27.9%) and potential epidemiological links (42.2%) compared to the 24-locus MIRU-VNTR genotyping alone (p < 0.001). Moreover, the proportion of clustering time span within 6 months increased significantly from 10.6% to 15.5% (p < 0.001). Among the paired strains with identical 24-locus MIRU-VNTR genotypes, the average geographic distance was 18.9 ± 17.0 km (Figure 1). With the combination including four additional loci, the paired distance was reduced to 13.6 ± 14.6 km (the corresponding value for IS6110-RFLP was 14.7 ± 15.2 km).
Figure legend: In this figure, the distance between parents’ residence among pairs of clustered strains in different genotyping strategies was presented on the box plot, marked with the mean value and standard deviation. C1 refers to 24-locus MIRU VNTR alone; C2 refers to the combination of 24-locus MIRU VNTR plus locus 3820; C3 refers to the combination of 24-locus MIRU VNTR plus loci 3820 and 3232; C4 refers to the combination of 24-locus MIRU VNTR plus loci 3820, 3232, and 2163a; C5 refers to the combination of 24-locus MIRU VNTR plus loci 3820, 3232, 2163a, and 4120; RFLP refers to IS6110-RFLP; Multiplicity test (LSD) of ANOVA showed that the paired geographic distance of C2 to C5 was significantly shorter than that of C1 (p = 0.001~0.011), which was also comparable to the RFLP findings.

3.5. Temporal Stability of MIRU-VNTR Loci and Combinations

In the analysis of clonal isolates, pairs of strained clones of 294 randomly selected strains (including 202 and 92 isolates from the Beijing and non-Beijing families respectively) were analyzed (Table S3). With the exception of locus 2347, the 29 (96.7%) remaining MIRU-VNTR loci showed variant patterns differing by one to two copies of MIRU-VNTR over two years’ observation. Among these loci, 12 out of 29 MIRU-VNTR loci (41.4%) exhibited a variation in the number of repeat copies in one or two isolates within the first 12 months, while 17 out of 29 (58.6%) varied between the 12th and 24th months. Furthermore, the analysis based on the GEE model revealed a significant variation in the locus 3820 (p = 0.015), while the remaining loci were stable over the two years’ observation.
As for the stability of MIRU-VNTR combinations, the Kaplan-Meier survival curves showed no separation of 24-locus MIRU-VNTR from other combinations over the two-year period (Figure 2). The proportion of temporally stable loci ranged from 90.5% to 92.5% within the four MIRU-VNTR combinations, which was comparable to 24 MIRU-VNTR (93.5%). IS6110-RFLP was the least stable, and the probability of no change at 24 months was 85.4%.

4. Discussion

MIRU-VNTR has proven its capacity for typing M. tuberculosis isolates in terms of discriminatory power and convenience. However, the epidemiological implication of MIRU-VNTR on the transmission of TB remains to be empirically evaluated [8,20]. Furthermore, the 24-locus MIRU-VNTR strategy needs to be evaluated in order to determine whether it is sufficiently stable and discriminative to allow meaningful epidemiological study. In this study, we showed that additional hypervariable MIRU-VNTR loci combined with 24-locus MIRU-VNTR not only provided a high discriminatory typing strategy, but also manifested temporal stability, and an improved epidemiological significance.
The heterogeneity of MIRU-VNTR loci has varied in previous studies [8,21,22]. These studies reported that loci 424, 802, 960, 1955, 3007, 3690, and 4348 had higher allelic diversity (h > 0.60). However, in this study, loci 424, 802, 960, 3690, and 4348 were moderately homogeneous (0.2 < h < 0.6), while loci 1955 and 3007 exhibited low allelic diversity (h = 0.144, 0.189). It is possible that the previous studies may have suffered from a problem of representativeness, due to the fact that the sample selection was restricted to only one setting [23]. Moreover, we found that the four hypervariable MIRU-VNTR loci with the highest allelic diversity were loci 3820, 3232, 2163a, and 4120 (h > 0.7), and these were further combined with the traditional 24-locus MIRU-VNTR to improve the typing resolution. However, this combination was inconsistent with Allix’s study of a 4-locus set (loci 1982, 3232, 3820, and 4120) [24], as well as Luo’s study of a 3-locus set (loci 3820, 3232 and 4120) [25]. It indicated the genetic diversity of hypervariable MIRU-VNTR loci varied in a specific setting, and the optimized MIRU-VNTR combination should be systematically evaluated further.
As previously reported, the inclusion of the most hypervariable MIRU-VNTR loci facilitates subdivision of the Beijing family, resulting in a significantly improved discriminatory power [8,25,26]. In this study, we found that the clustering proportion for the combination with all four hypervariable loci was reduced to 25.6%, which was comparable to that of IS6110-RFLP. Therefore, our data reveal that the addition of hypervariable MIRU-VNTR loci to the strategy currently in use could render the discriminative power equal to that of IS6110-RFLP typing.
Most importantly, we analyzed the epidemiological significance of clustering data by correlating them to the contact tracing information. Compared with traditional 24-locus MIRU-VNTR, the proportion of isolates with confirmed transmission links increased from 13.1% to 27.9% when clustered by the modified MIRU-VNTR genotyping. Moreover, the modified assay also showed a shorter clustering time span and geographic distance, which might validate its capability to indicate recent transmission chains in targeted population-based epidemiological studies.
A valid estimation of epidemiologically unrelated TB cases depends on the stability of the biomarker used for the genotyping. An additional aim of this study was to quantify the stability of the MIRU-VNTR combination, in comparison to the currently used genotyping methods, over a two year period, as proposed in previous studies [27,28]. In our observation, 12 loci displayed an altered copy number within the first year, and 17 loci varied between the 12th and 24th months. This result is in accordance with previous findings [29], which indicated that the average time needed for one copy number change in MIRU-VNTR genotyping was 18 months. Furthermore, the four MIRU-VNTR combinations remained stable among more than 90% of the isolates for up to two years, which was comparable to the currently used 24-locus MIRU-VNTR genotyping, and better than IS6110-RFLP typing. This suggests that the 24-locus MIRU-VNTR genotyping plus the additional hypervariable loci will estimate recent transmission in China more accurately and reliably than the 24-MIRU-VNTR alone.
Compared with whole genome sequencing (WGS), MIRU-VNTR typing has a lower discriminatory power [30], however the lower cost of its use facilitates its application and scale-up in most developing countries. Furthermore, previous reviews have also highlighted that the methodological variation in conducting and analyzing WGS data could lead to heterogeneous results that are difficult to compare [31]. Therefore, a genotyping method with a proper MIRU-VNTR combination might be more suitable for population-based molecular epidemiological investigations in developing countries.

5. Conclusions

This study evaluated multiple MIRU-VNTR loci in comparison to IS6110-RFLP in terms of temporal stability and epidemiological significance. The results of the combined MIRU-VNTR genotyping are well in line with those of IS6110 RFLP, and exhibit a relative stability suitable for indicating recent transmission. This study could constitute a basis for further molecular investigation of M. tuberculosis in China.

Supplementary Materials

The following are available online at https://www.mdpi.com/1660-4601/15/4/782/s1, Table S1: The number of repeats at each MIRU-VNTR loci in the studied M. tuberculosis isolates, Table S2: Genetic diversity (h values) of MIRU-VNTR loci in the studied M. tuberculosis isolates, Table S3: Temporal stability of MIRU-VNTR loci within 2 years’ in-vitro culture.

Acknowledgments

We would like to thank the health staff in the municipal and district CDCs for collecting data and clinical samples. The study was financially supported by National Natural Science Foundation of China (Grant No. 81373063), the Sweden-China (VR-NSFC) joint project NSFC Grant No. 81361138019, VR Grant No. 540-2013-8797) and the Fourth Round of Three-Year Public Health Action Plan of Shanghai, China (Grant No. 15GWZK0101).

Author Contributions

Yang Li carried out the data collection, performed the statistical analysis, and drafted the manuscript. Yi Hu developed the design and coordination, performed the molecular genotyping study, and helped to draft and finalize the manuscript. Mikael Mansjö, Solomon Ghebremichael, and Sven Hoffner have revised it critically for important intellectual content. Qi Zhao and Weili Jiang were involved in all the microbiological research. Biao Xu participated in the study design and interpretation. All authors read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cao, K.; Yang, K.; Wang, C.; Guo, J.; Tao, L.; Liu, Q.; Mahara, G.; Zhang, Y.; Guo, X. Spatial-Temporal Epidemiology of Tuberculosis in Mainland China: An Analysis Based on Bayesian Theory. Int. J. Environ. Res. Public Health 2016, 13, 469. [Google Scholar] [CrossRef] [PubMed]
  2. Marais, B.J.; Sintchenko, V. Epidemic spread of multidrug-resistant tuberculosis in China. Lancet Infect. Dis. 2016, 17, 238–239. [Google Scholar] [CrossRef]
  3. World Health Organization. Global Tuberculosis Report 2017. Latent Tuberculosis Infection. Available online: http://www.who.int/tb/publications/global_report/en/ (accessed on 1 November 2017).
  4. Allix-Beguec, C.; Fauville-Dufaux, M.; Supply, P. Three-year population-based evaluation of standardized mycobacterial interspersed repetitive-unit-variable-number tandem-repeat typing of Mycobacterium tuberculosis. J. Clin. Microbiol. 2008, 46, 1398–1406. [Google Scholar] [CrossRef] [PubMed]
  5. Thong-On, A.; Smittipat, N.; Juthayothin, T.; Yanai, H.; Yamada, N.; Yorsangsukkamol, J.; Chaiprasert, A.; Rienthong, D.; Billamas, P.; Palittapongarnpim, P. Variable-number tandem repeats typing of Mycobacterium tuberculosis isolates with low copy numbers of IS6110 in Thailand. Tuberculosis 2010, 90, 9–15. [Google Scholar] [CrossRef] [PubMed]
  6. Warren, R.M.; Victor, T.C.; Streicher, E.M.; Richardson, M.; van der Spuy, G.D.; Johnson, R.; Chihota, V.N.; Locht, C.; Supply, P.; van Helden, P.D. Clonal expansion of a globally disseminated lineage of Mycobacterium tuberculosis with low IS6110 copy numbers. J. Clin. Microbiol. 2004, 42, 5774–5782. [Google Scholar] [CrossRef] [PubMed]
  7. Alonso, H.; Samper, S.; Martín, C.; Otal, I. Mapping IS 6110 in high-copy number Mycobacterium tuberculosis strains shows specific insertion points in the Beijing genotype. BMC Genom. 2013, 14, 422. [Google Scholar] [CrossRef] [PubMed]
  8. Zheng, C.; Zhao, Y.; Zhu, G.; Li, S.; Sun, H.; Feng, Q.; Luo, M.; Wu, F.; Li, X.; Hill, V.; et al. Suitability of IS6110-RFLP and MIRU-VNTR for differentiating spoligotyped drug-resistant Mycobacterium tuberculosis clinical isolates from Sichuan in China. BioMed Res. Int. 2014, 2014, 763204. [Google Scholar] [CrossRef] [PubMed]
  9. Joseph, B.V.; Soman, S.; Radhakrishnan, I.; Hill, V.; Dhanasooraj, D.; Kumar, R.A.; Rastogi, N.; Mundayoor, S. Molecular epidemiology of Mycobacterium tuberculosis isolates from Kerala, India using IS6110-RFLP, spoligotyping and MIRU-VNTRs. Infect. Genet. Evol. 2013, 16, 157–164. [Google Scholar] [CrossRef] [PubMed]
  10. 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, 4498–4510. [Google Scholar] [CrossRef] [PubMed]
  11. Oelemann, M.C.; Diel, R.; Vatin, V.; Haas, W.; Rusch-Gerdes, S.; Locht, C.; Niemann, S.; Supply, P. Assessment of an optimized mycobacterial interspersed repetitive- unit-variable-number tandem-repeat typing system combined with spoligotyping for population-based molecular epidemiology studies of tuberculosis. J. Clin. Microbiol. 2007, 45, 691–697. [Google Scholar] [CrossRef] [PubMed]
  12. Van Belkum, A.; Tassios, P.T.; Dijkshoorn, L.; Haeggman, S.; Cookson, B.; Fry, N.K.; Fussing, V.; Green, J.; Feil, E.; Gerner-Smidt, P.; et al. Guidelines for the validation and application of typing methods for use in bacterial epidemiology. Clin. Microbiol. Infect. 2007, 13, 1–46. [Google Scholar] [CrossRef] [PubMed]
  13. Iwamoto, T.; Yoshida, S.; Suzuki, K.; Tomita, M.; Fujiyama, R.; Tanaka, N.; Kawakami, Y.; Ito, M. Hypervariable loci that enhance the discriminatory ability of newly proposed 15-loci and 24-loci variable-number tandem repeat typing method on Mycobacterium tuberculosis strains predominated by the Beijing family. Fems Microbiol. Lett. 2007, 270, 67–74. [Google Scholar] [CrossRef] [PubMed]
  14. Mokrousov, I.; Narvskaya, O.; Vyazovaya, A.; Millet, J.; Otten, T.; Vishnevsky, B.; Rastogi, N. Mycobacterium tuberculosis Beijing genotype in Russia: In search of informative variable-number tandem-repeat loci. J. Clin. Microbiol. 2008, 46, 3576–3584. [Google Scholar] [CrossRef] [PubMed]
  15. Luo, T.; Yang, C.; Pang, Y.; Zhao, Y.; Mei, J.; Gao, Q. Development of a Hierarchical Variable-Number Tandem Repeat Typing Scheme for Mycobacterium tuberculosis in China. PLoS ONE 2014, 9, e89726. [Google Scholar] [CrossRef] [PubMed]
  16. Somerville, W.; Thibert, L.; Schwartzman, K.; Behr, M.A. Extraction of Mycobacterium tuberculosis DNA: A question of containment. J. Clin. Microbiol. 2005, 43, 2996–2997. [Google Scholar] [CrossRef] [PubMed]
  17. Van Embden, J.D.; Cave, M.D.; Crawford, J.T.; Dale, J.W.; Eisenach, K.D.; Gicquel, B.; Hermans, P.; Martin, C.; McAdam, R.; Shinnick, T.M.; et al. Strain identification of Mycobacterium tuberculosis by DNA fingerprinting: Recommendations for a standardized methodology. J. Clin. Microbiol. 1993, 31, 406–409. [Google Scholar] [PubMed]
  18. Maes, M.; Kremer, K.; van Soolingen, D.; Takiff, H.; de Waard, J.H. 24-locus MIRU-VNTR genotyping is a useful tool to study the molecular epidemiology of tuberculosis among Warao Amerindians in Venezuela. Tuberculosis 2008, 88, 490–494. [Google Scholar] [CrossRef] [PubMed]
  19. Struelens, M.J. Consensus guidelines for appropriate use and evaluation of microbial epidemiologic typing systems. Clin. Microbiol. Infect. 1996, 2, 2–11. [Google Scholar] [CrossRef] [PubMed]
  20. Jonsson, J.; Hoffner, S.; Berggren, I.; Bruchfeld, J.; Ghebremichael, S.; Pennhag, A.; Groenheit, R. Comparison between RFLP and MIRU-VNTR genotyping of Mycobacterium tuberculosis strains isolated in Stockholm 2009 to 2011. PLoS ONE 2014, 9, e95159. [Google Scholar] [CrossRef] [PubMed]
  21. Guo, J.H.; Xiang, W.L.; Zhang, G.; Luo, T.; Xie, N.; Yang, Z.R.; Sun, Q. Mycobacterial interspersed repetitive unit typing in Mycobacterium tuberculosis isolates from Sichuan province in China. Indian J. Med. Res. 2011, 134, 362–368. [Google Scholar] [PubMed]
  22. Liu, Y.; Wang, S.; Lu, H.; Chen, W.; Wang, W. Genetic diversity of the Mycobacterium tuberculosis Beijing family based on multiple genotyping profiles. Epidemiol. Infect. 2016, 144, 1728–1735. [Google Scholar] [CrossRef] [PubMed]
  23. Yun, K.W.; Song, E.J.; Choi, G.E.; Hwang, I.K.; Lee, E.Y.; Chang, C.L. Strain typing of Mycobacterium tuberculosis isolates from Korea by mycobacterial interspersed repetitive units-variable number of tandem repeats. Korean J. Lab. Med. 2009, 29, 314–319. [Google Scholar] [CrossRef] [PubMed]
  24. Allix-Béguec, C.; Wahl, C.; Hanekom, M.; Nikolayevskyy, V.; Drobniewski, F.; Maeda, S.; Campos-Herrero, I.; Mokrousov, I.; Niemann, S.; Kontsevaya, I. Proposal of a Consensus Set of Hypervariable Mycobacterial Interspersed Repetitive-Unit–Variable-Number Tandem-Repeat Loci for Subtyping of Mycobacterium tuberculosis Beijing Isolates. J. Clin. Microbiol. 2014, 52, 164–172. [Google Scholar] [CrossRef] [PubMed]
  25. Lu, W.; Lu, B.; Liu, Q.; Dong, H.; Shao, Y.; Jiang, Y.; Song, H.; Chen, C.; Li, G.; Xu, W.; et al. Genotypes of Mycobacterium tuberculosis isolates in rural China: Using MIRU-VNTR and spoligotyping methods. Scand. J. Infect. Dis. 2014, 46, 98–106. [Google Scholar] [CrossRef] [PubMed]
  26. Sun, Z.; Cao, R.; Tian, M.; Zhang, X.; Zhang, X.; Li, Y.; Xu, Y.; Fan, W.; Huang, B.; Li, C. Evaluation of Spoligotyping and MIRU-VNTR for Mycobacterium bovis in Xinjiang, China. Res. Vet. Sci. 2012, 92, 236–239. [Google Scholar] [CrossRef] [PubMed]
  27. Niemann, S.; Rusch-Gerdes, S.; Richter, E.; Thielen, H.; Heykes-Uden, H.; Diel, R. Stability of IS6110 restriction fragment length polymorphism patterns of Mycobacterium tuberculosis strains in actual chains of transmission. J. Clin. Microbiol. 2000, 38, 2563–2567. [Google Scholar] [PubMed]
  28. Savine, E.; Warren, R.M.; van der Spuy, G.D.; Beyers, N.; van Helden, P.D.; Locht, C.; Supply, P. Stability of variable-number tandem repeats of mycobacterial interspersed repetitive units from 12 loci in serial isolates of Mycobacterium tuberculosis. J. Clin. Microbiol. 2002, 40, 4561–4566. [Google Scholar] [CrossRef] [PubMed]
  29. Mazars, E.; Lesjean, S.; Banuls, A.L.; 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, 1901–1906. [Google Scholar] [CrossRef] [PubMed]
  30. Roetzer, A.; Diel, R.; Kohl, T.A.; Rückert, C.; Nübel, U.; Blom, J.; Wirth, T.; Jaenicke, S.; Schuback, S.; Rüsch-Gerdes, S.; et al. Whole Genome Sequencing versus Traditional Genotyping for Investigation of a Mycobacterium tuberculosis Outbreak: A Longitudinal Molecular Epidemiological Study. PLoS Med. 2013, 10, e1001387. [Google Scholar] [CrossRef] [PubMed]
  31. Hatherell, H.A.; Colijn, C.; Stagg, H.R.; Jackson, C.; Winter, J.R.; Abubakar, I. Interpreting whole genome sequencing for investigating tuberculosis transmission: A systematic review. BMC Med. 2016, 14, 21. [Google Scholar] [CrossRef] [PubMed]
Ijerph 15 00782 g001 550
Figure 1. Paired geographic distance in different genotyping strategies.
Figure 1. Paired geographic distance in different genotyping strategies.
Ijerph 15 00782 g001
Ijerph 15 00782 g002 550
Figure 2. Evaluation of stability in different genotyping strategies applied on 294 randomly selected isolates.
Figure 2. Evaluation of stability in different genotyping strategies applied on 294 randomly selected isolates.
Ijerph 15 00782 g002
Table
Table 1. Clinical characteristics of patients with tuberculosis in different geographic areas of China.
Table 1. Clinical characteristics of patients with tuberculosis in different geographic areas of China.
VariateNumber of Patients (%)
East
n = 385		North
n = 211		South
n = 169		West
n = 217
Age (mean ± SD)48 ± 19.045 ± 18.347 ± 19.348 ± 19.8
Male (Sex)262 (68.1)135 (64.0)128 (75.7)154 (71.0)
Cavity on CXR109 (28.3)50 (23.7)38 (22.5)55 (25.3)
Sputum smear positive233 (60.5)118 (55.9)107 (63.3)137 (63.1)
Table
Table 2. The discriminative power of the different combinations of Mycobacterial Interspersed Repetitive Units-Variable Number of Tandem Repeats (MIRU-VNTR) loci.
Table 2. The discriminative power of the different combinations of Mycobacterial Interspersed Repetitive Units-Variable Number of Tandem Repeats (MIRU-VNTR) loci.
MIRU-VNTR Loci CombinationCumulative Hunter-Gaston Discriminatory Index (HGDI)
TotalStrain Family
BeijingNon-Beijing
24-locus MIRU VNTR0.98870.98280.9922
+38200.99090.98970.9931
+3820 + 32320.99220.99150.9943
+3820 + 3232 + 2163a0.99370.99220.9959
+3820 + 3232 + 2163a + 41200.99540.99380.9961
Table
Table 3. Comparison of technical validation in different genotyping strategies.
Table 3. Comparison of technical validation in different genotyping strategies.
MethodNo. of PatternNo. of Unique PatternNo. of Clustered Pattern (Isolates)Clustering Proportion (%)HGDI
24-locus MIRU VNTR700577123 (405)41.20.9854
+382076666898 (314)32.00.9935
+3820 + 323278669789 (285)29.00.9947
+3820 + 3232 + 2163a79571679 (266)27.10.9952
+3820 + 3232 + 2163a + 412081373182 (251)25.60.9960
IS6110-RFLP81172685 (256)26.10.9958
Table
Table 4. Comparison of epidemiological links and clustering time span in different genotyping strategies.
Table 4. Comparison of epidemiological links and clustering time span in different genotyping strategies.
MethodEpidemiology Link (%)p-Value aTime of Clustering Time Span (%)p-Value a
ConfirmedPotentialNot Confirmed<6 Months6–12 Months12–24 Months24–36 Months≥36 Months
24-locus MIRU VNTR53 (13.1)89 (22.0)263 (64.9)<0.00143 (10.6)40 (9.9)92 (22.7)148 (36.5)82 (20.2)<0.001
+382069 (22.0)110 (35.0)135 (43.0) 43 (13.7)49 (15.6)84 (26.8)100 (31.8)38 (12.1)
+3820 + 323269 (24.2)108 (37.9)108 (37.9) 41 (14.4)48 (16.8)74 (26.0)89 (31.2)33 (11.6)
+3820 + 3232 + 2163a67 (25.2)106 (39.8)93 (35.0) 39 (14.7)47 (17.7)67 (25.2)81 (30.5)32 (12.0)
+3820 + 3232 + 2163a + 412070 (27.9)106 (42.2)75 (29.9) 39 (15.5)43 (17.1)68 (27.1)72 (28.7)29 (11.6)
IS6110-RFLP67 (26.2)105 (41.0)84 (32.8) 4 7 (18.4)45 (17.6)40 (15.6)89 (34.8)35 (13.7)
ap-value was analyzed by Cochran-Mantel-Haenszel test.

Share and Cite

MDPI and ACS Style

Li, Y.; Hu, Y.; Mansjö, M.; Zhao, Q.; Jiang, W.; Ghebremichael, S.; Hoffner, S.; Xu, B. The Epidemiological Significance and Temporal Stability of Mycobacterial Interspersed Repetitive Units-Variable Number of Tandem Repeats-Based Method Applied to Mycobacterium tuberculosis in China. Int. J. Environ. Res. Public Health 2018, 15, 782. https://doi.org/10.3390/ijerph15040782

AMA Style

Li Y, Hu Y, Mansjö M, Zhao Q, Jiang W, Ghebremichael S, Hoffner S, Xu B. The Epidemiological Significance and Temporal Stability of Mycobacterial Interspersed Repetitive Units-Variable Number of Tandem Repeats-Based Method Applied to Mycobacterium tuberculosis in China. International Journal of Environmental Research and Public Health. 2018; 15(4):782. https://doi.org/10.3390/ijerph15040782

Chicago/Turabian Style

Li, Yang, Yi Hu, Mikael Mansjö, Qi Zhao, Weili Jiang, Solomon Ghebremichael, Sven Hoffner, and Biao Xu. 2018. "The Epidemiological Significance and Temporal Stability of Mycobacterial Interspersed Repetitive Units-Variable Number of Tandem Repeats-Based Method Applied to Mycobacterium tuberculosis in China" International Journal of Environmental Research and Public Health 15, no. 4: 782. https://doi.org/10.3390/ijerph15040782

APA Style

Li, Y., Hu, Y., Mansjö, M., Zhao, Q., Jiang, W., Ghebremichael, S., Hoffner, S., & Xu, B. (2018). The Epidemiological Significance and Temporal Stability of Mycobacterial Interspersed Repetitive Units-Variable Number of Tandem Repeats-Based Method Applied to Mycobacterium tuberculosis in China. International Journal of Environmental Research and Public Health, 15(4), 782. https://doi.org/10.3390/ijerph15040782

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop