- Research article
- Open Access
Detection of Beijing strains of MDR M. tuberculosis and their association with drug resistance mutations in katG, rpoB, and embB genes
BMC Infectious Diseases volume 20, Article number: 752 (2020)
Molecular epidemiological studies of Mycobacterium tuberculosis (MTB) are the core of current research to find out the association of the M. tuberculosis genotypes with its outbreak and transmission. The high prevalence of the Beijing genotype strain among multidrug resistance (MDR) TB has already been reported in various studies around India. The overall objective of this study was to detect the prevalence of Beijing genotype strains of MDR M. tuberculosis and their association with the clinical characteristics of TB patients.
In this study 381 M. tuberculosis clinical isolates were obtained from sputum samples from 2008 to 2014. The multiplex-PCR and Spoligotyping (n = 131) methods were used to investigate the prevalence of the Beijing genotype strain by targeting the Rv2820 gene and their association with drug resistance and clinical characteristics of TB patients. The drug susceptibility testing of first-line anti-TB drugs was performed by using the proportion method and MGIT960. A collection of isolates having Beijing and non-Beijing strains were also characterized to see if Beijing genotype strains had a higher rate of mutations at codons 516, 526 and 531 of the 81-bp region of the rpoB gene, codon 315 of the katG gene, and codon 306 of the embB gene.
The sensitivities and specificities of multiplex-PCR assay compared to that of standard Spoligotyping was detected to be 100%. Further, we observe that the multi drug-resistance was significantly associated with Beijing genotype strains (p = 0.03) and a strong correlation between Beijing genotype strains and specific resistance mutations at the katG315, rpoB531, and embB306 codons (p = < 0.0001, < 0.0001 & 0.0014 respectively) was also found.
This rapid, simple, and cost-effective multiplex PCR assay can effectively be used for monitoring the prevalence of Beijing genotype strains in low resource settings. Findings of this study may provide a scientific basis for the development of new diagnostic tools for detection and effective management of DR-TB in countries with a higher incidence rate of Beijing genotype strains.
Beijing genotype of M. tuberculosis (MTB) strain attracted expert’s attention due to its association with multidrug resistance (MDR) and disease outbreaks . Justifications for the rapid dissemination of Beijing genotype strain have included its putative elevated virulence, enhanced transmissibility, greater mutability, ability to escape from BCG vaccine-induced immunity, and its skill to acquire multidrug resistance [2,3,4,5]. Information about the transmission of Beijing-genotype MTB strain circulating in North India is lacking. Therefore, as a part of our overall objective we have attempted to assess the prevalence of Beijing genotype MTB strain and its relation to drug resistance and clinical characteristics of TB patients in North India.
An extensive review of the mechanism behind successful emergence of the Beijing genotype strain has revealed a higher frequency of mutation at the codon 315 of katG gene among MTB Beijing genotype strains . However, the association is less clear for mutations at codons 516, 526, and 531 of the rpoB gene, and at the codon 306 of the embB gene . We characterized a collection of MTB patient isolates to see if the Beijing genotype strain had a higher rate of mutations at the katG315 codon, the rpoB516, rpoB526, and rpoB531 codons of RRDR region which accounts for > 95% of resistance to rifampicin (RIF), and embB306, which is responsible for resistance to ethambutol (EMB) [7, 8]. Further, we also have evaluated a multiplex PCR against standard spoligotyping for the rapid and cost-effective detection of Beijing genotypes of MTB strains.
Banaras Hindu University (BHU) hospital is a tertiary-care hospital with vast catchment area, being the only tertiary care hospital in Eastern Uttar Pradesh (UP), situated in North Central India, providing medical cover to over 200 million population of Eastern UP, Western Bihar and adjoining areas of Madhya Pradesh. It’s a tertiary care center with a referral bias towards non-responding cases.
TB patients, age range 6–90 years, with or without any other additional complications such as HIV-seropositivity, have been included in this study. Other details were received from the medical files and compliance charts. Treatment of TB patients was done according to the standard short-course chemotherapy under DOTS guidelines . The written informed consent was obtained from all participants.
Collection of samples
Four thousand seventy nine sputum specimens were collected during 2008–2014, from the patients of pulmonary tuberculosis, from three TB – centers and analyzed for culture positivity to Mycobacterium tuberculosis.
A total of 381 MTB isolates were isolated from 4079 TB suspected specimens collected over 6 years (2008 to 2014) from TB patients attending different health care centers and BHU hospital of Varanasi, UP, India.
Drug susceptibility testing
A total of 381 MTB isolates were subjected to DST according to the standard proportion method (PM) and commercial MGIT960 system against RIF, INH, streptomycin (STR), and EMB . H37Rv (ATCC 27294) and known RIF, INH, and EMB mono-resistant strains were used as negative and positive controls respectively.
Three rapid HIV test kits based on different antigens/principles were used according to NACO (National AIDS Control Organization) guidelines for the detection of HIV in patients with TB .
Multiplex-allele specific-PCR (MAS-PCR)/DNA sequencing for detection of RIF, INH and EMB resistance determinants
DNA was isolated as described by van Embden et al and a two-step MAS-PCR assay was performed to detect mutations at rpoB codons (516, 526, and 531), katG codon 315 and embB codon 306 [12,13,14]. Automated sequencing was performed using the BigDye Terminator kit v3.1 (Applied Biosystems, USA) as per manufacturer’s protocol. DNA sequence analysis was carried out for 61 strains. Further comparisons of rpoB (n = 128), katG (n = 128) and embB306 (n = 27) were carried out with SeqScape® v2.5 software.
Out of 381 MTB isolates, spoligotyping was performed for 131 isolates according to the standard method of Kamerbeek et al. . The spoligotypes were compared with those contained in the international database SITVIT2, an updated variant of the previously released SpolDB4 database . Isolates of the Beijing genotype were defined by showing hybridization to at least three of the nine spacers within 35–43 and showing the absence of hybridization to spacers 1 to 34. Some Beijing strains lack one or more of the nine signature spacers, due to asymmetrical insertions of IS6110 (Table 3).
Multiplex-PCR for the identification of Beijing genotypes of MTB targeting Rv2820 gene
To differentiate between Beijing and non-Beijing strains of 381 MTB isolate, three sets of PCR primers were used (Supplementary Table 1 [17, 18]). These were based on the information given by Warren et al. that the region covering genes Rv2816 to Rv2819 and portion of Rv2820 is deleted in all Beijing strains of MTB . For multiplex PCR, 10 μl master mix was containing 0.1 U of DyNAzyme EXT (Finnzymes Oy, Espoo, Finland), 250 μM of each dNTP, 2.0 mM MgCl2, 5 μM of each primer, and 1 μl (10 ng) of template DNA in 1× EXT buffer. Reaction conditions included 96 °C for 4 min, followed by 35 cycles of 96 °C for 45 s, 55 °C for 45 s, and 72 °C for 1 min, and a final extension at 72 °C for 5 min. The amplified products were analyzed by electrophoresis on a 3% agarose gel. Beijing strains of MTB will produce one 523-bp internal control amplified product and one 129-bp product while non-Beijing strains will produce the 523-bp and a104-bp product (Fig. 1).
Results were analyzed with the SPSS software (version 16.0.0) (SPSS Inc., Chicago, IL, USA). The features of the two groups were compared using the Z-test and of three groups by the chi-square (χ2) test, for statistical significance assessment.
A total of 381 MTB isolates were isolated. Seventy-six (19.95%) isolates belong to Beijing genotype among which 45 (11.81%) were multidrug-resistant. Beijing genotype strains were significantly more common in women (p = 0.01) than in men (Table 1).
Of the 32 strains that were Beijing genotype as defined by Spoligotyping (total Beijing genotype strains are 76), 30 had all nine characteristic spacers 35–43, corresponding to the shared type SIT1 as defined in SITVIT2. In the remaining two isolates, one isolate had the last six spacers while another had spacers 37–43 and spacer 35 correspondings to SIT250 and SIT621 respectively (Table 2).
Among 76 Beijing strains, 47 were RIF resistant, 50 were resistant to INH, and 45 were found to be MDR. Thirty-seven strains were resistant to STR and 40 were EMB resistant (Supplementary Table 2). While 50 (65.79%) and 45 (59.21%) of the Beijing strains were any drug-resistant and MDR respectively, 144 (47.21%) and 76 (24.92%) of 305 non-Beijing strains were any drug-resistant and MDR during the same period representing a statistically significant difference (p-value 0.004 & 0.032 respectively) (Table 1).
Multiplex-PCR Vs Spoligotyping
Multiplex PCR was performed for all 381 isolates and spoligotyping was performed in 131 (of 381) isolates. When statistically compared, the results of all 131 isolates were the same for the detection of Beijing and non-Beijing isolates with both the methods. The sensitivities and specificities of multiplex-PCR assay compared to that of standard Spoligotyping was observed to be 100%.
Genotyping for RIF, INH and EMB resistance mutations
A mutation in katG gene, showing a nucleotide substitution at codon 315 from AGC to ACC, founded in 43/50 (86%) of the Beijing strains with phenotypic INH resistance. Resistance mutations at the rpoB gene were found in 45/47 (95.74%) Beijing strains with phenotypic RIF resistance. A mutation at rpoB codon 531 was most common (n = 29; 61.70%). Of 44 phenotypically EMB resistant isolates, 21 had mutations at embB codon 306 (Table 3). Compared to the wild-type (wt) isolates, the isolates with the katG315Thr and rpoB531Leu mutations were most frequent (86 and 61.70% respectively) in Beijing genotypes (Table 3). Mutations at katG315, rpoB531, and embB306 were more common among Beijing genotype strains, and these associations were observed to be statistically significant (p = < 0.0001, < 0.0001 and 0.0014 respectively) (Table 4).
Beijing genotype strain was first isolated in Beijing, China in 1995 and now it has been predominantly found in various parts of the world [20, 21]. A systematic review and meta-analysis on the prevalence of the Beijing genotype strain in world population revealed that its prevalence was 44.7% in Asia, 27.9% in Europe, 12.5% in Africa, and 8.9% in America . In India, the Beijing genotype strain is present in a low percentage throughout the country including a previous study in our region [23,24,25] however we observed a comparatively higher incident rate of the Beijing genotype strain (19.95%, 76/381). Possible reason behind this discrepancy could be the low number of MTB strains studied in the previous study in our region.
We found a 100% agreement between multiplex-PCR targeting the Rv2820 gene and standard spoligotyping for the detection of Beijing genotype strains of MTB. Similar results were also observed by Sun et al and Sherafat et al [26, 27].
In an extensive review based on the data available from 49 studies in 35 countries on Beijing genotype MTB strain and drug resistance four different patterns of Beijing genotype strains were suggested worldwide but our study population represents a different group; endemic and associated with drug resistance which is now a common pattern for some other countries too [5, 28, 29].
The worldwide emergence of Beijing genotype strains and their association with multidrug-resistance substantiates the fact that they possess strong capability to acquire drug resistance. This is confirmed by a study, in which 14 drug-resistant clinical isolates of the Beijing genotype were sequenced by whole-genome sequencing to understand the evolution of drug resistance and results showed that the drug resistance among these isolates appeared to be acquired, not clonally derived . Our study observed the similar drift, with 33/50 (66%) drug resistant-Beijing genotype strains and 29/45 (64.44%) MDR-Beijing genotype strains from acquired drug-resistant TB patients.
The higher rate of mutations at rpoB531, katG315, and embB306 codons among MTB, Beijing strains provides further evidence that drug resistance may contribute to the global predominance of Beijing genotype strains. In the present study, mutations at codons katG315, rpoB531 and embB306 were more common among Beijing genotype strains than in non-Beijing genotype strains of MTB (p = < 0.0001; < 0.0001; 0.0014) (Table 4). Similar findings were also observed in other studies [6, 31,32,33,34]. On the basis of our results we consider that recent transmission of Beijing genotype strains could be the driving force of such a high prevalence of the katG315, rpoB531 mutations. Another explanation could be; mutations in putative mutator genes have been found in Beijing genotypes, which results in an altered DNA repair and an increased mutation adaptability rate .
The strong association of Beijing genotype strain and katG, rpoB, and embB mutations among treated cases of TB, reflects the considerable transmission of MDR- Beijing genotype strains in our region which suggests the need of urgent change in the quality of chemotherapy applied in the community.
In this study, we have successfully evaluated a simple, rapid, robust, and cost-effective multiplex PCR assay for the detection of Beijing genotype of M. tuberculosis strains by targeting Rv2820 gene. Study outcomes may provide a scientific basis for the development of new DR-TB diagnostic tools for effective management and control of DR-TB in countries with a higher incidence rate of Beijing genotype strains. Further, a recruitment of large number of representative samples for a systematic study is needed to elucidate the role of Beijing genotype strain in the transmission of TB and their association with DR-TB in the region. Additional studies using MIRU-VNTRs over an extended period of time can be done to fully understand the epidemiology of TB and its transmission dynamics in our region and all over India.
Availability of data and materials
The datasets used and analyzed during the current study available from the corresponding author on reasonable request.
Rifampicin resistance determining region
Banaras Hindu University
Venkataraman P, Paramasivan CN. Drug resistance in tuberculosis and issues related to multidrug resistance in planning for TB control in India. Health Adm Spec Issue Tuberc. 2002;15:127–36.
Glynn JR, Whiteley J, Bifani PJ, Kremer K, van Soolingen D. Worldwide occurrence of Beijing/W strains of mycobacterium tuberculosis: a systematic review. Emerg Infect Dis. 2002;8:843–9.
Tracevska T, Jansone I, Baumanis V, Marga O, Lillebaek T. Prevalence of Beijing genotype in Latvian multidrug-resistant mycobacterium tuberculosis isolates. Int J Tuberc Lung Dis. 2003;7:1097–103.
Pang Y, Song Y, Xia H, Zhou Y, Zhao B, Zhao Y. Risk factors and clinical phenotypes of Beijing genotype strains in tuberculosis patients in China. BMC Infect Dis. 2012;12:354.
Parwati I, van Crevel R, van Soolingen D. Possible underlying mechanisms for successful emergence of the Mycobacterium tuberculosis Beijing genotype strains. Lancet Infect Dis. 2010;10:103–11.
San LL, Aye KS, Oo NAT, Shwe MM, Fukushima Y, Gordon SV, Suzuki Y, Nakajima C. Insight into multidrug-resistant Beijing genotype Mycobacterium tuberculosis isolates in Myanmar. Int J Infect Dis. 2018;76:109–19.
Hazbon MH. Bobadilla del Valle M, Guerrero MI et al. role of embB codon 306 mutations in Mycobacterium tuberculosis revisited: a novel association with broad drug resistance and IS6110 clustering rather than ethambutol resistance. Antimicrob Agents Chemother. 2005;49:3794–802.
Gupta A, Singh SK, Anupurba S. Mutations at embB306 codon and their association with multidrug-resistance. tuberculosis clinical isolates. Ind J Med Microbiol. 2015;33:387–92.
Arora VK, Gupta R. Dots strategy in India - the challenges. Curr Med J North Zone. 2002;VIII:4.
Canetti G, Fox W, Khomenko A, et al. Advances in techniques of testing mycobacterial drug sensitivity, and the use of sensitivity tests in tuberculosis control programs. Bull World Health Organ. 1969;41:21–43.
Guidelines for HIV testing. National AIDS control organisation: Ministry of Health and Family Welfare; 2007. http://www.who.int/hiv/pub/guidelines/india_art.pdf.
Van Embden JDA, Cave MD, Crawford JT, et al. Strain identification of Mycobacterium tuberculosis by DNA fingerprinting: recommendations for a standardized methodology. J Clin Microbiol. 1993;1:406–9.
Gupta A, Prakash P, Singh SK, et al. Rapid genotypic detection of rpoB and katG gene mutations in Mycobacterium tuberculosis clinical isolates from northern India as determined by MAS-PCR. J Clin Lab Anal. 2013;27:31–7.
Mokrousov I, Narvskaya O, Limeschenko E. Detection of ethambutol-resistant Mycobacterium tuberculosis strains by multiplex allele-specific PCR assay targeting embB306 mutations. J Clin Microbiol. 2002;40:1617–20.
Kamerbeek J, Schouls L, Kolk A, et al. Simultaneous detection and strain differentiation of Mycobacterium tuberculosis for diagnosis and epidemiology. J Clin Microbiol. 1997;35:907–14.
Brudey K, Driscoll JR, Rigouts L, 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.
Hillemann D, Warren R, Kubica T, et al. Rapid detection of Mycobacterium tuberculosis Beijing genotype strains by real-time PCR. J Clin Microbiol. 2006;44:302–6.
Plikaytis BB, Marden JL, Crawford JT. Multiplex PCR assay specific for the multidrug-resistant strain W of Mycobacterium tuberculosis. J Clin Microbiol. 1994;32:1542–6.
Warren RM, Victor TC, Streicher EM, et al. Patients with active tuberculosis often have different strains in the same sputum specimen. Am J Respir Crit Care Med. 2004;169:610–4.
Van Soolingen D, Qian L, de Haas PE, Douglas JT, Traore H, Portaels F, et al. Predominance of a single genotype of Mycobacterium tuberculosis in countries of East Asia. J Clin Microbiol. 1995;33:3234–8.
Liu Q, Wang D, Martinez L, Lu P, Zhu L, Lu W, Wang J. Mycobacterium tuberculosis Beijing genotype strains and unfavourable treatment outcomes: a systematic review and meta-analysis. Clin Microbiol Infect. 2020;26(2):180–8.
Ramazanzadeh R, Sayhemiri K. Prevalence of Beijing family in Mycobacterium tuberculosis in world population: systematic review and meta-analysis. Inter J Mycobact. 2014;2:109–13.
Sharma P, Chauhan DS, Upadhyay P, et al. Molecular typing of Mycobacterium tuberculosis isolates from a rural area of Kanpur by spoligotyping and mycobacterial interspersed repetitive units (MIRUs) typing. Infect Genet Evol. 2008;8:621–6.
Stavrum R, Myneedu VP, Arora VK, et al. In-depth molecular characterization of Mycobacterium tuberculosis from New Delhi predominance of drug-resistant isolates of the ‘modern’ (TbD1) type. PLoS One. 2009;4:4540.
Mathuria JP, Srivastava GN, Sharma P, Mathuria BL, Ojha S, Katoch VM, Anupurba S. Prevalence of mycobacterium tuberculosis Beijing genotype and its association with drug resistance in North India. J Infect Pub Heal. 2017;10:409–14.
Sun JR, Lee SY, Dou HY. Using a multiplex polymerase chain reaction for the identification of Beijing strains of Mycobacterium tuberculosis. Eur J Clin Microbiol Infect Dis. 2009;28:105–7.
Sherafat SJ, Goudarzi H, Farnia P. Evaluation of two molecular methods for the detection of Beijing strains of Mycobacterium tuberculosis in pulmonary tuberculosis patients. HealthMED. 2013;7(1):12–6.
European Concerted Action on New Generation Genetic Markers and Techniques for the Epidemiology and Control of Tuberculosis. Beijing/W genotype Mycobacterium tuberculosis and drug resistance. Emerg Infect Dis. 2006;12:736–43.
Jiao W, Liu Z, Han R, et al. A country-wide study of spoligotype and drug resistance characteristics of Mycobacterium tuberculosis isolates from children in China. PLoS One. 2013;30(8):12–e84315.
Ioerger TR, Feng Y, Chen X, et al. The non-clonality of drug resistance in Beijing-genotype isolates of Mycobacterium tuberculosis from the Western cape of South Africa. BMC Genomics. 2010;11:670.
Lipin MY, Stepanshina VN, Shemyakin IG, et al. Association of specific mutations in katG, rpoB, rpsL, and rrs genes with spoligotypes of multidrug-resistant Mycobacterium tuberculosis isolates in Russia. Clin Microbiol Infect. 2007;13:620–6.
Park YK, Shin S, Ryu S, et al. Comparison of drug resistance genotypes between Beijing and non-Beijing family strains of Mycobacterium tuberculosis in Korea. J Microbiol Methods. 2005;63(2):165–72.
Ghebremichael S, Groenheit R, Pennhag A, Koivula T, Andersson E, et al. Drug-resistant mycobacterium tuberculosis of the Beijing genotype does not spread in Sweden. PLoS One. 2010;5(5):e10893. https://doi.org/10.1371/journal.pone.0010893.
Safi H, Sayers B, Hazbón MH. Transfer of embB codon 306 mutations into clinical Mycobacterium tuberculosis strains alters susceptibility to ethambutol, isoniazid, and rifampin. Antimicrob Agents Chemother. 2008;52:2027–34.
Financial assistance to Dr. Anamika Gupta from the Indian Council of Medical Research is also highly acknowledged.
This study was financially supported by the Indian Council of Medical Research, New Delhi. The funders had no role in study design, data collection, and analysis, decision to publish, or preparation of the manuscript.
Ethics approval and consent to participate
This study has been approved by the ethical committee of the Institute of Medical Sciences, Banaras Hindu University, Varanasi, India.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Gupta, A., Sinha, P., Nema, V. et al. Detection of Beijing strains of MDR M. tuberculosis and their association with drug resistance mutations in katG, rpoB, and embB genes. BMC Infect Dis 20, 752 (2020). https://doi.org/10.1186/s12879-020-05479-5
- Beijing genotype
- Drug resistant-TB