Abstract
Background:
Rifampicin resistant tuberculosis is a serious problem faced by tuberculosis control in China, and rapid detection of rifampicin resistance is urgently needed.Objectives:
This study aimed to describe the molecular characteristics and frequency of RNA polymerase β subunit (rpoB) gene mutations in rifampicin-resistant tuberculosis (RR-TB) in the Anqing area.Methods:
The rpoB gene fragment was amplified by polymerase chain reaction (PCR), and all isolates were sequenced for mutations in the rpoB gene. The mutations were obtained by comparing the sequencing results with the MUBII database. In addition, logistic regression was used to analyze the relationship between rpoB mutations and rifampicin (RIF) resistance.Results:
There were 152 males and 42 females in this study, and the mean age was 56.60 ± 17.91 years. Mutations in the rpoB gene were a risk factor for rifampicin resistance (β = 5.271, P < 0.001 OR = 195.192). Among the 19 RR-TB strains, 16 (84.21%) had mutations in the ropB gene, and three (1.71%) of 175 rifampicin-sensitive strains were mutated. The mutation sites of five strains (31.58%) were at the codon 526 and five strains (31.58%) at the codon 531. However, there were two strains at the codon 513 and two strains at the codon 533 (15.79%), and two strains (10.53%) were double mutations.Conclusions:
The mutation characteristics of the rpoB gene in the Anqing area are complex, and rpoB mutation detection can be used as an indicator to screen drug resistance of RIF.Keywords
1. Background
In 2020, the World Health Organization (WHO) reported that 78% of rifampicin-resistant tuberculosis (RR-TB) strains were multidrug-resistant TB (MDR-TB) and that the prevalence of MDR-TB was 7.1% and 21% among new and previously treated TB cases in China in 2018 (1). WHO announced that MDR-TB was resistant to at least isoniazid (INH) and rifampicin (RIF) of anti-TB drugs. According to previous studies, most RIF-resistant isolates are MDR and also resistant to INH. Thus, RIF resistance is thought to be a reliable biomarker for MDR-TB detection (2-4).
Previous studies had confirmed that RR-TB was mainly associated with RNA polymerase β subunit (rpoB) gene mutation. Nearly 96% of RR-TB cases can be attributed to a point mutation in a specific region of this gene, known as the RIF resistance determining region (RRDR) between the rpoB codons 507 and 533 of Mycobacterium tuberculosis (MTB) (5, 6). The rpoB gene mutations are usually single, or multiple base mutations, and about 80%-86% of base mutations occur in three amino acid sites: 531 Serine to Leucine, 526 Histidine to Tyrosine, and 526 Histidine to other amino acids. Besides, different mutations in this region could lead to different degrees of drug resistance. For example, the mutation at 518, 516, 511, and 522 only produces low-grade drug resistance.
The Xpert MTB/RIF system based on the mutation of the rpoB gene to detect RR-TB has been accurately and rapidly approved by WHO (7). Thus, few sites have been used to study the relation between rpoB mutations and rifampicin resistance, except for 531 and 526 codon mutations. It is worth exploring whether there are other sites in the rpoB gene that can be used for improving the sensitivity and specificity of the prediction of rifampicin resistance. Compared with other areas in China, the investigation of mutation characteristics of drug-resistant genes was inadequate in Anhui Province (8). Besides, mutations in rpoB gene studies of Anhui Province only involved several places, such as Maan Shan and Hefei (9), and the results of the actual rpoB gene mutation level need to be clarified. While drug resistance often shows regional differences, such as biological characteristics and functions of bacteria (10), investigating the rpoB gene mutation characteristics of clinical strains of MTB in Anqing of Anhui Province is still of important practical significance for the prevention and treatment of tuberculosis (TB) in Anhui Province.
2. Objectives
This study aimed to describe the molecular characteristics and frequency of RNA polymerase β subunit (rpoB) gene mutations in rifampicin-resistant tuberculosis (RR-TB) in the Anqing area.
3. Methods
3.1. Study Setting and Samples Collection
A total of 194 DNA samples from sputum were collected from tuberculosis patients who went to Anqing Center for Diseases Prevention and Control (CDC) for drug sensitivity tests from September 2018 to December 2019. Verbal consent was obtained from all participants and legal guardians of participants under 18. The H37Rv standard strain was produced by the Center of Biological Products of the Department of Health (Beijing, China). We completed all procedures in the P2 Laboratory of Anqing CDC, which met the second biosafety level.
3.2. Drug Sensitivity Test
A total of 194 MTB isolates successfully cultured from sputum samples were used in this study. According to the "diagnostic laboratory test procedures for tuberculosis," the drug sensitivity test of MTB clinical isolates was carried out to determine the primary susceptibility to RIF by proportional method on Lowenstein-Jensen (L-J) medium containing the corresponding anti-TB drugs (11).
3.3. DNA Extraction and PCR Amplification
Genomic DNA of MTB isolates was extracted from L-J cultures using a mericon TM DNA Bacteria Kit. The designed primers were synthesized using SanGonBiotech. The sequences of primers were as follows (5’→3’): Forward: GTTGTAAAACGACGGCCAGGGATGACCACCCAGGA Reverse: CAGGAAACAGCTATGACGGTTTAGATCGGCACAT. The PCR conditions were: 95˚C for 5 min of denaturation, then 98°C for 10 seconds, 50 - 60°C for 30 seconds, and 72°C for 1 min, for 40 cycles in total. After the last 5 minutes at 72°C, the PCR products were taken out and stored in the refrigerator at 4°C (Figure 1).
Detection of the amplification products of rpoB by a polymerase chain reaction in clinical sputum specimens. *Lane Marker, DL500 marker (500, 400, 300, 250, 200, 150, and 100 bp); lane 1 - 4, clinical rifampicin-resistant strains with rpoB mutation.
3.4. DNA Sequencing
The sequence determination of the DNA sequences of the PCR product was detected by Illumina paired-end sequencing. The accession number of the mutated gene in NCBI was OP784366-OP784384. DNAStar software was used to splice forward and reverse sequences, and the MUBII-TB-D (12) (https://umr5558-bibiserv.univ-lyon1.fr/mubii/mubii-in.cgi) database was used to analyze the mutation sites of rpoB genes.
3.5. Statistical Analysis
We collected and analyzed the data in Statistical Package for the Social Sciences (SPSS) version 18.0. The demographics of the participants were determined as mean (S.D). Logistic regression was used to analyze the relationship between rpoB mutations and rifampicin resistance. A P-value less than 0.05 was considered to be significant for all analyses.
4. Results
In the total 194 cases of MTB isolated from the 200 sputum samples, 19 RIF-resistant strains and 175 RIF-sensitive strains were detected. Males and females accounted for 78.35% and 21.65% of these patients, respectively. The mean age of the patients was 56.60±17.91 years, with extreme ages at 4 and 90 years. Among the 19 RIF-resistant isolates identified, we detected 14 (73.68%) with a single mutation, two (10.53%) with double mutations, and three (15.79%) without mutation in the ropB gene (Table 1). Nucleotide substitutions in codons 531 and 526 were recorded as single mutation (31.25%). The other hot spot codons, namely 513 and 533, were mutated in 12.5% of the cases, respectively. The double mutations recorded were founded in codons 516 and 518 and codons 505 and 511 (Figure 2). Of note, the mutations of codons 505 and 518 were novel observations. Among 305 RIF-sensitive cases, codon changes in 526, 531, and 533 were recorded. Previous studies have shown that rpoB mutation is associated with RIF resistance, and we have verified this conclusion based on our experiments. Logistic regression was used to analyze the results. It was shown that mutations in the rpoB gene were a risk factor for RIF resistance (β = 5.274, P < 0.001, OR = 195.192) (Table 2).
The rpoB DNA codon mutation information of clinical isolates of Mycobacterium tuberculosis. *Due to the influence of the naming method, for example, the mutation of 1295 in this picture is located at codon 432, which needs to add 81, to become codon 513 in Table 1.
Distribution of rpoB Gene Mutations Among 27 RIF-resistant Strains of Mycobacterium tuberculosis
Gene | Position | Amino Acid Changes | Codon Changes | Strains | Frequency (%) |
---|---|---|---|---|---|
rpoB | 513 | Gln-Pro | CAA-CCA | 1 | 6.25 |
513 | Gln-Glu | CAA-GAA | 1 | 6.25 | |
526 | His-Gly | CAC-GGC | 1 | 6.25 | |
526 | His-Asn | CAC-AAC | 1 | 6.25 | |
526 | His-Asp | CAC-GAC | 1 | 6.25 | |
526 | His-Tyr | CAC-TAC | 2 | 12.50 | |
531 | Ser-Leu | TCG-TTG | 5 | 31.25 | |
533 | Leu-Pro | CTG-CCG | 2 | 12.50 | |
516,518 | Asp-Tyr, Asn-His | GAC-TAC, AAC-CAC | 1 | 6.25 | |
505,511 | Phe-Leu, Leu-Pro | TTC-TTG, CTG-CCG | 1 | 6.25 |
The Logistic Regression Analysis Between rpoB Mutations and Rifampicin-drug Resistance
Variables | β | Wald | P-Value | OR | 95% CI |
---|---|---|---|---|---|
rpoB mutation | 5.274 | 51.557 | < 0.001 | 195.192 | 46.265 - 823.522 |
Age | -0.020 | 1.721 | 0.190 | 0.980 | 0.950 - 1.010 |
Gender | -0.411 | 0.382 | 0.536 | 0.663 | 0.180 - 2.440 |
Constant | -2.004 | 5.074 | 0.024 | 0.135 |
5. Discussion
Drug-resistant pulmonary tuberculosis requires extended treatment and more expensive drugs, which has become a severe public health issue in many developing countries (13). TB patients who do not take medicine in time, change the medicine or stop the medicine without permission may lead to drug resistance of MTB (14). There were fewer rifampicin-resistant strains in this study because the patients from which the strains originated were mainly newly treated cases. There are noticeable regional differences in the characteristics of rpoB gene mutation. The variation spectrum and variation rate of the rpoB gene of RR-TB in different countries and regions show different characteristics.
In China, most rpoB gene mutations were concentrated in RRDR, including 533, 526, 531, 530, 526, 522, 516, 513, 511, and 508, of which 531 and 526 were the most common codon mutations (15). Up to now, more than 80 mutations of unit point codon, 70 combined mutations of two or three codons, more than 20 deletion mutations, and four insertion mutation types have been found (16). The present studies found that the mutations in this area were mainly located in codons 531, 526, and other codons of RRDR. The primary mutation type was single base substitution (17).
Our results were similar to previous findings. The most common codon mutations were 531 (Ser→Leu) and 526 (His→Gly/Asn/Asp/Tyr), occurring by base substitutions. The same situation was discovered in Hebei, Jiangxi, and Huainan of Anhui province (15, 18, 19). The former site included mainly five strains of Ser→Leu (TCG→TTG), while the latter contained two strains of His→Tyr (CAC→TAC), one of His→Gly (CAC→GGC), one strain of His→Asp (CAC→AAC) and one of His→Asp (CAC→GAC) substitutions. In the current study, 89.5% of the isolates had mutations in the RRDR, which was lower than the results reported in Shanghai and Jiangxi provinces.
A study in Hunan province reported a lower mutation frequency in the 81-bp core region of rpoB than our result. However, the results remained consistent with most places worldwide (19). In addition, we also identified some new mutations, such as AAC518CAC (Asn→His) and TTC505TTG (Phe→Leu). The types of mutations at codons 518 (AAC→GAC/ATC) had been reported in other areas of China (20), and it is worth noting that the mutation type (AAC→CAC) of codons 518 was first found as well as the 505 codons in China.
In this study, we also found that 40.74% (11/27) of the rifampicin-resistant strains did not have RRDR mutations. Present studies suggested a hypothesis that the mutations in these resistant strains may be caused by mutations of genes outside the rpoB gene or due to other drug-resistant mechanisms (21, 22). There are two main theories about rifampicin resistance: one is the mutation of RNA polymerase β subunit caused by mutations of the rpoB gene; the other is that the permeability of the cell wall changes (23). Our experiment may be due to the latter, and further research is needed. Besides, mutations of the rpoB gene were found in three rifampicin-sensitive strains, which were noticed at codons 533, 526, and 531, respectively. These mutations may cause by low-concentration rifampicin resistance. Some strains sensitive to proportioning or absolute concentration methods may have a low resistance to RIF. In this situation, the drug susceptibility testing result could be sensitive, but the rpoB gene occurred mutations (24). In this study, rifampicin-resistant strains only accounted for a small part of all MTB strains, and we did not involve relevant research on the minimum inhibitory concentration. We will complement relevant research in a future study.
5.1. Conclusions
The use of rpoB gene sequencing revealed the molecular characteristics of rifampicin resistance of MTB in Anqing of China. Our data indicated that rpoB gene mutation could be a risk factor for rifampicin resistance in Anhui Province. In this study, we found standard types of mutations and two new types of mutations. It could provide data support for developing a new molecular detection method for MDR-TB and be of great significance to better adjuvant clinical medication.
References
-
1.
World Health Organisation. Global tuberculosis report 2020. World Health Organisation; 2020.
-
2.
Ahmad S, Mokaddas E, Fares E. Characterization of rpoB mutations in rifampin-resistant clinical Mycobacterium tuberculosis isolates from Kuwait and Dubai. Diagn Microbiol Infect Dis. 2002;44(3):245-52. [PubMed ID: 12493171]. https://doi.org/10.1016/s0732-8893(02)00457-1.
-
3.
Karimi H, Oudghiri A, En-Nanei L, Mzibri ME, Laglaoui A, Chaoui I, et al. Frequency of genomic mutations mediating resistance of Mycobacterium tuberculosis isolates to rifampicin in Northern Morocco. Rev Inst Med Trop Sao Paulo. 2020;62. e37. [PubMed ID: 32520212]. [PubMed Central ID: PMC7274765]. https://doi.org/10.1590/S1678-9946202062037.
-
4.
Valim AR, Rossetti ML, Ribeiro MO, Zaha A. Mutations in the rpoB gene of multidrug-resistant Mycobacterium tuberculosis isolates from Brazil. J Clin Microbiol. 2000;38(8):3119-22. [PubMed ID: 10921994]. [PubMed Central ID: PMC87207]. https://doi.org/10.1128/JCM.38.8.3119-3122.2000.
-
5.
Zhang Y, Yew WW. Mechanisms of drug resistance in Mycobacterium tuberculosis: update 2015. Int J Tuberc Lung Dis. 2015;19(11):1276-89. [PubMed ID: 26467578]. https://doi.org/10.5588/ijtld.15.0389.
-
6.
Ramaswamy S, Musser JM. Molecular genetic basis of antimicrobial agent resistance in Mycobacterium tuberculosis: 1998 update. Tuber Lung Dis. 1998;79(1):3-29. [PubMed ID: 10645439]. https://doi.org/10.1054/tuld.1998.0002.
-
7.
Steingart KR, Schiller I, Horne DJ, Pai M, Boehme CC, Dendukuri N. Xpert(R) MTB/RIF assay for pulmonary tuberculosis and rifampicin resistance in adults. Cochrane Database Syst Rev. 2014;2014(1). CD009593. [PubMed ID: 24448973]. [PubMed Central ID: PMC4470349]. https://doi.org/10.1002/14651858.CD009593.pub3.
-
8.
Yao S, Yan J, Li L, Ma D, Liu J, Wang Q, et al. Determining Mycobacterium tuberculosis drug resistance and risk factors for multidrug-resistant tuberculosis in sputum smear-positive tuberculosis outpatients in Anhui Province, China, 2015-2016. Infect Drug Resist. 2020;13:1023-32. [PubMed ID: 32308444]. [PubMed Central ID: PMC7156340]. https://doi.org/10.2147/IDR.S244482.
-
9.
Zhao Y, Xu S, Wang L, Chin DP, Wang S, Jiang G, et al. National survey of drug-resistant tuberculosis in China. N Engl J Med. 2012;366(23):2161-70. [PubMed ID: 22670902]. https://doi.org/10.1056/NEJMoa1108789.
-
10.
Jing W, Pang Y, Zong Z, Wang J, Guo R, Huo F, et al. Rifabutin resistance associated with double mutations in rpoB gene in Mycobacterium tuberculosis isolates. Front Microbiol. 2017;8:1768. [PubMed ID: 28959248]. [PubMed Central ID: PMC5603767]. https://doi.org/10.3389/fmicb.2017.01768.
-
11.
Neshani A, Kamali Kakhki R, Sankian M, Zare H, Hooshyar Chichaklu A, Sayyadi M, et al. Modified genome comparison method: a new approach for identification of specific targets in molecular diagnostic tests using Mycobacterium tuberculosis complex as an example. BMC Infect Dis. 2018;18(1):517. [PubMed ID: 30314460]. [PubMed Central ID: PMC6186143]. https://doi.org/10.1186/s12879-018-3417-x.
-
12.
Flandrois JP, Lina G, Dumitrescu O. MUBII-TB-DB: a database of mutations associated with antibiotic resistance in Mycobacterium tuberculosis. BMC Bioinformatics. 2014;15:107. [PubMed ID: 24731071]. [PubMed Central ID: PMC4021062]. https://doi.org/10.1186/1471-2105-15-107.
-
13.
Biru D, Woldesemayat EM. Determinants of drug-resistant tuberculosis in Southern Ethiopia: A case-control study. Infect Drug Resist. 2020;13:1823-9. [PubMed ID: 32606831]. [PubMed Central ID: PMC7305943]. https://doi.org/10.2147/IDR.S256536.
-
14.
Gunther G, Ruswa N, Keller PM. Drug-resistant tuberculosis: advances in diagnosis and management. Curr Opin Pulm Med. 2022;28(3):211-7. [PubMed ID: 35220372]. [PubMed Central ID: PMC9415219]. https://doi.org/10.1097/MCP.0000000000000866.
-
15.
Li Q, Wang Y, Li Y, Gao H, Zhang Z, Feng F, et al. Characterisation of drug resistance-associated mutations among clinical multidrug-resistant Mycobacterium tuberculosis isolates from Hebei Province, China. J Glob Antimicrob Resist. 2019;18:168-76. [PubMed ID: 30926465]. https://doi.org/10.1016/j.jgar.2019.03.012.
-
16.
Lu J, Ye S, Li CP, Qin R, Li WP. [Sequence analysis on drug-resistant gene of rpoB in MDR-TB among pneumoconiosis patients complicated with tuberculosis in Huainan mining district]. Zhonghua Lao Dong Wei Sheng Zhi Ye Bing Za Zhi. 2012;30(8):579-81. Chinese. [PubMed ID: 22931875].
-
17.
Lu J, Jiang S, Ye S, Deng Y, Ma S, Li CP. Sequence analysis of the drug‑resistant rpoB gene in the Mycobacterium tuberculosis L‑form among patients with pneumoconiosis complicated by tuberculosis. Mol Med Rep. 2014;9(4):1325-30. [PubMed ID: 24535137]. https://doi.org/10.3892/mmr.2014.1948.
-
18.
Luo D, Chen Q, Xiong G, Peng Y, Liu T, Chen X, et al. Prevalence and molecular characterization of multidrug-resistant M. tuberculosis in Jiangxi province, China. Sci Rep. 2019;9(1):7315. [PubMed ID: 31086215]. [PubMed Central ID: PMC6513856]. https://doi.org/10.1038/s41598-019-43547-2.
-
19.
Naghdi S, Moradi T, Tavangar F, Bahrami G, Shahboulaghi M, Ghiasvand H. The barriers to achieve financial protection in Iranian health system: A qualitative study in a developing country. Ethiop J Health Sci. 2017;27(5):491-500. [PubMed ID: 29217954]. [PubMed Central ID: PMC5615010]. https://doi.org/10.4314/ejhs.v27i5.7.
-
20.
Sheng J, Li J, Sheng G, Yu H, Huang H, Cao H, et al. Characterization of rpoB mutations associated with rifampin resistance in Mycobacterium tuberculosis from eastern China. J Appl Microbiol. 2008;105(3):904-11. [PubMed ID: 18422555]. https://doi.org/10.1111/j.1365-2672.2008.03815.x.
-
21.
Wang M, Xu G, Lu L, Xu K, Chen Y, Pan H, et al. Genetic polymorphisms of IL-17A, IL-17F, TLR4 and miR-146a in association with the risk of pulmonary tuberculosis. Sci Rep. 2016;6:28586. [PubMed ID: 27339100]. [PubMed Central ID: PMC4919632]. https://doi.org/10.1038/srep28586.
-
22.
Postma DS, Bush A, van den Berge M. Risk factors and early origins of chronic obstructive pulmonary disease. Lancet. 2015;385(9971):899-909. [PubMed ID: 25123778]. https://doi.org/10.1016/S0140-6736(14)60446-3.
-
23.
Uddin MKM, Rahman A, Ather MF, Ahmed T, Rahman SMM, Ahmed S, et al. Distribution and frequency of rpoB mutations detected by Xpert MTB/RIF assay among Beijing and Non-Beijing rifampicin resistant Mycobacterium tuberculosis isolates in Bangladesh. Infect Drug Resist. 2020;13:789-97. [PubMed ID: 32210593]. [PubMed Central ID: PMC7073589]. https://doi.org/10.2147/IDR.S240408.
-
24.
Bemer-Melchior P, Bryskier A, Drugeon HB. Comparison of the in vitro activities of rifapentine and rifampicin against Mycobacterium tuberculosis complex. J Antimicrob Chemother. 2000;46(4):571-6. [PubMed ID: 11020254]. https://doi.org/10.1093/jac/46.4.571.