https://orcid.org/0000-0002-1961-4454
ABSTRACT
Reports on antimicrobial resistance (AMR) of Mycobacterium leprae, relationship with bacteriological index (BI), and transmission in China are limited. We investigated the emergence of AMR mutations, the relationship between BI and AMR in complete, moderate and lack of BI decline cases, and molecular epidemiological features of AMR cases by enrolling 290 leprosy cases from four endemic provinces. Seven (2.41%), one (0.34%), five (1.72%), one (0.34%), and one (0.34%) strains had single mutations in folP1, rpoC, gyrA, gyrB, and 23S rRNA, respectively. Double mutations in folP1 and gyrA, rpoB and gyrA, and gyrA and 23S rRNA were observed in one (0.34%) strain each. Mutated strains occurred in three out of 81 (95% CI−0.005-0.079, p = 0.083) cases with complete BI decline, in seven out of 103 (95% CI 0.018-0.117, p = 0.008) cases with moderate BI decline, and in four out of 34 (95% CI 0.003-0.231, p = 0.044) cases with lack of BI decline. Most of these mutated strains were geographically separated and diverged genotypically. AMR mutations may not be the main cause of the lack of BI decline. The low transmission of AMR strains at the county level indicates an ongoing transmission at close contact levels.
Introduction
Leprosy is a chronic granulomatous disease that mainly affects the skin, peripheral nerves, and mucous membrane [Citation1]. The introduction of multidrug therapy (MDT) with dapsone, rifampicin, and clofazimine has decreased the prevalence of leprosy worldwide. However, more than 216,108 new leprosy cases were reported in 2016, indicating the continuous transmission of the disease [Citation2]. According to National reports, China estimated 634 new cases in 2017 [Citation2]. The reduction in leprosy prevalence has been observed in recent years in China; however, a considerable number of new, relapse, and drug-resistant cases occurred representing a major public health concern [Citation3].
Drug resistance detection of WHO-recommended DRDRs containing folP1, rpoB, and gyrA associated with dapsone, rifampicin, and ofloxacin resistance respectively is well established and has been incorporated in the WHO guidelines for AMR surveillance. There are several AMR reports available regarding dapsone [Citation4], rifampicin [Citation5], and ofloxacin [Citation6] resistance in leprosy. In addition to mutations in WHO-recommended DRDRs, compensatory AMR-associated mutations, including nth (DNA repair), rpoA (rifampicin), rpoB, rpoC (rifampicin), gyrA, gyrB (ofloxacin), and 23S rRNA (clarithromycin) have been reported in Mycobacterium tuberculosis and Mycobacterium leprae. However, the roles of these mutations in M. leprae are yet to be confirmed [Citation7–9].
The BI measurement in the slit-skin smear of patients with potential leprosy is used for diagnosis, classification, and clinical recovery and relapse detection. BI at diagnosis and during treatment, release from treatment (RFT), and follow-up have previously been studied to justify the efficacy of MDT. Generally, 6-months and 1-year MDT courses are recommended for paucibacillary (PB) and multibacillary (MB) cases respectively. Nevertheless, no effective BI change was observed in some leprosy cases, even if MDT was prescribed for more than 1-year (Supplementary Table 1) [Citation10–19]. The factors influencing the lack of BI decline remain unknown, but high initial BI and no decline in BI at RFT are possible signs of AMR [Citation12].
Molecular epidemiological studies for strain typing of M. leprae have been well developed to track strain transmission and geographical distribution in endemic and non-endemic regions [Citation20]. M. leprae isolates having identical genotypes are expected to be associated with recent transmission. To our knowledge, no studies on the transmission of AMR strains of M. leprae at the county level in China have been performed.
In the current study, we analysed the AMR of M. leprae by direct PCR sequencing of the WHO-recommended DRDRs containing folP1, rpoB, and gyrA. We also investigated the sequence of regions and genes of extended DRDRs including nth, rpoA, rpoB, rpoC, gyrA, gyrB, and 23S rRNA to evaluate the necessity to complement the current surveillance system for AMR mutations. AMR mutations were further studied for their relationship with BI and to identify mutation trends in new and relapse cases. Finally, we performed a population-based study of leprosy by using genomic VNTR and epidemiological features of wild-type and mutated strains of M. leprae at the county level to examine the transmission of mutated strains of leprosy in China.
Materials and methods
Study sites
The detailed study profile is depicted in . The study was conducted under the national AMR surveillance programme from 2013 to 2017 in four endemic provinces namely Guizhou, Hunan, Sichuan, and Yunnan. These provinces account for > 80% of all leprosy cases in China. Some of the districts of these provinces have a high burden of leprosy with a prevalence rate is about > 1/100,000. The study area is located between the 20° N and 30° N latitudes and has tropical and subtropical climate conditions ().
Figure 1. Schematic procedure and structure of the study. BI status at diagnosis was compared to BI at 1-year post-RFT, which identified three categorized patients: complete decline, moderate decline, and lack of BI decline. Patients under treatment were excluded from this grouping. M. leprae VNTR typing and molecular drug resistance analysis were performed for all 290 patients.
Figure 2. Study sites and geographical distribution of AMR strains of M. leprae. Four provinces shown in colour were included in the surveillance for AMR of leprosy. AMR mutations were not detected in Sichuan province. All AMR mutation types detected in Yunnan, Guizhou, and Hunan provinces are indicated.
Patient enrolment, BI monitoring, and follow-up of cases
During the study period, 270 newly diagnosed and 20 relapse cases of leprosy in four provinces were enrolled. For relapse screening, retrospective data of RFT patients from 2005 to 2010 were collected. Diagnosis and prescription of PB and MB-MDT regimen were performed as per WHO guidelines. Initial BI was recorded at the time of MDT regimen initiation, and cases were followed up by professional health workers periodically until the completion of 6-months or 1-year MDT. The BI was assessed at the time of RFT and 1-year post-RFT. Three categories namely complete BI decline, moderate BI decline, and lack of BI decline were categorized according to BI comparison at initial and 1-year post-RFT.
PCR and sequencing of DRDRs
The presence of M. leprae genomic DNA was confirmed by PCR detection of the 16S rRNA gene [Citation21]. PCR and nested PCR were performed to amplify the WHO-recommended DRDRs (folP1, rpoB, and gyrA) and regions and genes of extended DRDRs including nth, rpoA, rpoB, rpoC, gyrA, gyrB, and 23S rRNA by using primers and conditions listed in Supplementary Tables 2–6. The quality of amplified products was assessed using 1–2% agarose gel. PCR amplicons were purified and sequenced at the local commercial sequencing company (Tsingke Biological Technology Co., Ltd., Nanjing, China). Alignments to NCBI nucleotide databases were performed using the basic local alignment sequence tool (NCBI-BLAST) to detect the mutations.
Molecular epidemiological analysis
The strain types and geographical distances within and between counties and provinces were analysed. Counties with drug-resistant strains were further studied and correlated with other counties, forming clear clusters with drug-resistant strains. Seventeen VNTR loci, AC8b, (GTA)9, (GGT)5, (AT)17, rpoT, 21-3, (AC)9, (AT)15, (AC)8a, 27-5, 6-7, (TA)18, (TTC)21,18-8, 12-5, 23-3, and (TA)10 were used to perform genotyping of M. leprae strains [Citation22]. A phylogenetic tree was constructed using the ggtree package of R using the neighbour joining method.
Statistical analysis
Data were statistically analysed using IBM SPSS software version 22. The Chi-square test was used to compare the parametric continuous data amongst different study groups. P < 0.05 was considered statistically significant.
Results
Clinical and BI details of MDT completed cases and cases under treatment
The demographical and clinical details of the patients from the four provinces are shown in . Amongst the 290 cases, 218 had completed MDT, and the rest were under treatment. Comparison of initial BI to BI at RFT for each of the 218 patients indicated three trends: complete decline, moderate decline, and lack of decline (including increase). Initial BIs (measured on a scale of 0–6) were stratified by one unit intervals, and their levels at RFT and follow-up are plotted in Supplementary Figure 1. Most of the BI distributions showed a gradual BI decline pattern from initial to 1-year post-RFT, except 5 to < 6 BI distribution, exhibiting a rapid BI decline pattern from initial to 1-year post-RFT. On the other hand, we observed a gradual BI decline pattern from RFT to 1-year post-RFT in all BI distributions, including 5 to < 6. The mean initial and RFT BIs of 81 cases (including 38 cases with initial BI zero), whose complete BI declines were 1.08 (range 0–4.8) and 0.21 (range 0–3.6) respectively were found to be BI negative and non-symptomatic at 1-year post-RFT. In 103 cases with moderate BI decline, the mean initial BI was 3.82 (range 1.2–5.7), which declined to 2.5 (range 0.4–4.8) at RFT and further to 2.04 (range 0.2–4.3) at 1-year post-RFT. A total of 34 cases (including 23 cases with increased BI) had BIs that did not decline, and they did not exhibit any clinical improvement. The mean BI at initial, RFT, and 1-year post-RFT were 2.65 (range 0–5.8), 3.45 (range 1–5.8), and 3.45 (range 1–5.8) respectively (). A different plot shows BI distribution at three-time points for complete, moderate, and lack of BI decline (Supplementary Figure 2). In the complete BI decline group, BI rapidly declined from the initial to the RFT and remained stable towards the 1-year post-RFT. The stable BI decline pattern at three time periods, including initial, RFT, and 1-year post-RFT were observed in the moderately BI decline group. Instead of stable BI pattern from initial to RFT in the lack of BI decline group, we observed a rise in BI from initial to RFT due to increased BI load during the treatment in some cases of lack of BI decline (Supplementary Figure 2). Other cases are currently under treatment, and their relevant details are summarized in Supplementary Table 6.
Table 1. Clinical and other relative information of all cases from four provinces.
Table 2. Clinical characteristics and AMR results of BI completely decline, BI moderately decline, and lack of BI decline cases of leprosy.
Molecular AMR analysis
WHO-recommended DRDRs: AMR strains of M. leprae were identified only in Guizhou, Hunan, and Yunnan, and none of them was identified in Sichuan ( and ). Amongst the 290 cases, 11 (3.7%) strains showed drug resistance, and the remaining (96.3%) were sensitive to all the drugs. Drug resistance analysis of the 290 cases revealed eight single drug-resistant (seven, 2.4% dapsone and one, 0.34% ofloxacin) and two double drug-resistant (one, 0.34% each to dapsone and ofloxacin, rifampicin and ofloxacin) strains of M. leprae. Another double drug resistance to ofloxacin and clarithromycin was observed in one strain (0.34%). No correlation of drug resistance proportion was identified in male vs. female (9/200 vs. 2/90; p = 0.512), nerve vs. no nerve involvement (10/220 vs. 1/70; p = 0.470), deformity vs. no deformity (3/103 vs. 8/187; p = 0.751), and reaction vs. non-reaction (1/40 vs. 10 /250; p = 1.000) cases. None of the mutated strains was harboured in PB. The AMR strains were detected in only LL and BL cases (5/72 vs. 6/136; p = 0.518). The frequency of drug resistance in relapse cases was higher compared to that of new cases (3/20 vs. 8/270; p = 0.032). Interestingly, strains from new cases showed single drug resistance to dapsone or ofloxacin, while relapse cases exhibited double drug resistance to dapsone and ofloxacin, rifampicin and ofloxacin, and ofloxacin and clarithromycin ().
Table 3. Clinical characteristics and other relevant information of all AMR cases of leprosy.
Extended DRDRs: Among 290 cases, no mutations were identified in nth, rpoA, and rpoB. The four single mutations in gyrA were detected in different strains (1.37%), whereas each of the rpoC, gyrB, and 23S rRNA mutations were detected in one strain each (0.34%) ().
The screening of WHO-recommended and extended DRDRs revealed 18 mutated strains among 290 cases of leprosy. Similar to WHO-recommended DRDRs, no correlation of drug resistance proportion was observed in male vs. female (15/200 vs. 3/90; p = 0.200), nerve vs. no nerve involvement (16/220 vs. 2/70; p = 0.258), deformity vs. no deformity (4/103 vs. 14/187; p = 0.310), and reaction vs. non-reaction (4/40 vs. 14/250; p = 0.288) cases of leprosy. No mutations existed in PB, whereas all mutations were assigned to MB (18/270) cases. The TT and BB cases did not have mutations, whereas LL (6/72) followed by BL (10/136) and BT (2/53) harboured mutations (p = 0.583). The mutation frequency was significantly higher in relapse cases than that in new cases (5/20 vs. 13/270; p = 0.004) ().
Relationship between BI and AMR
Analysis of the relationship between the BI and AMR mutations in WHO-recommended DRDRs revealed one folP1 mutated strain amongst complete BI decline cases, five strains of folP1 mutations and one strain with double mutation in rpoB and gyrA amongst moderate BI decline cases, and another gyrA mutated strain amongst cases of lack of BI decline (p = 0.251) ( and ). The mutated strain analysis of cases under treatment is shown in Supplementary Table 6.
The relationship between the BI and AMR mutations in DRDRs (WHO and extended) identified mutations in three out of 81 cases (folP1, rpoC, and gyrA mutated strain in each) with complete BI decline, seven out of 103 cases (five folP1, one 23S rRNA, and one both rpoB and gyrA) with moderate BI decline, and four out of 34 cases (three gyrA and one gyrB) with lack of BI decline (p = 0.268) ( and ).
The emergence of novel mutations
In addition to previously described mutations, one strain presented mutation Asp698Thr in rpoC, conferring rifampicin resistance. Gly362Asp in gyrA and Val214Gly in gyrB were novel quinolone resistance mutations identified in four and one strains of M. leprae respectively. Other mutations at A2142C and A2143C in 23S rRNA conferring clarithromycin resistance were identified in one strain each. We further retrieved the mutations in previously treated cases and identified the emergence of Thr53Ile mutation in folP1 in one patient previously treated with DDS monotherapy. Mutation analysis of other previous cases treated with MDT drugs and newly diagnosed cases with no previous history of MDT treatment revealed acquisition of different mutations, confirming the emergence of drug resistance of leprosy ().
Geographical and genetic clustering analysis
Guizhou, Sichuan, and Yunnan have common geographical borders, whereas Hunan shares border only with Guizhou (). A total of 18 strains of M. leprae with mutations were distributed across three provinces in 15 counties geographically well separated. The Qixingguan, Huayuan, and Guangnan counties had two mutated strains each identified in different villages and families. Further investigations revealed no evidence of the migration of any of the cases infected with these 18 strains at the village, county, and province levels.
The VNTR genetic distances of M. leprae revealed that most of the mutated strains had unique genotypes and were also genetically well separated from wild-type strains (). Genotype cluster analysis of wild-type and mutated strains of M. leprae at county level was performed using strains from patients residing in the counties where mutated strains were detected. We observed that the mutated strains had no clear cluster pattern, indicating genetic diversity even in the same county (). Interestingly, one clear cluster was observed with wild-type (n = 5) and mutated strains (n = 2) in the Huayuan county of Hunan; however, mutated strains of this cluster originated from a different village, and no migration information of the patients was observed ().
Figure 3. Phylogeny of wild-type and mutated strains from all four provinces. Genetic clustering pattern based on variable number tandem repeats for all wild-type and mutated strains (290 strains) of the counties of studied provinces. All branches are highlighted with four colours indicate four provinces. Mutated strains are represented with coloured rectangle columns. Each sample information is represented with three letters, where the 1st, 2nd, and 3rd letter specifies province, city, and county respectively.
Figure 4. Phylogeny of wild-type and mutated strains from counties where drug resistance strains were detected. In counties where DR strains were detected, a phylogenetic relationship to local wild-type strains (94 strains) was examined based on their variable number tandem repeat genotypes. All branches are highlighted with four colours designate four provinces. Mutated strains are represented with coloured rectangle columns. Each sample information is represented with three letters, where the 1st, 2nd, and 3rd letter specifies province, city, and county respectively.
Discussion
WHO declares that MDT treatment is effectively controlling leprosy incidence and contributing to the elimination of the leprosy burden in several countries, including China. The emergence of drug resistance in leprosy is a major concern for the implementation of disease intervention programmes. As very limited studies described the prevalence of gene mutations associated with AMR, the present study highlighted the characterization of the emergence of drug resistance by analysing DRDRs and some extended potential DRDR genes of M. leprae in China.
In newly diagnosed cases of leprosy, the high frequency of AMR to dapsone (six, 2.2%) in WHO-recommended DRDRs was identified in the present study, which was less than that in India [Citation23], Korea [Citation24], Myanmar, Indonesia, Philippines [Citation25], Vietnam [Citation26], and Brazil [Citation3] but higher than that in Colombia [Citation12] and China [Citation27]. In this study, no mutation associated with rifampicin resistance was confirmed, whereas rifampicin resistance cases were confirmed in India [Citation23], Brazil [Citation28], China [Citation27], Myanmar, Indonesia, and Philippines [Citation25]. We reported that the AMR frequency for ofloxacin (two, 0.7%) was less in the current study leprosy population than that of Colombia [Citation12], India [Citation3], and China [Citation27] leprosy population.
In relapse cases, the AMR to dapsone (one, 5%) and double mutation in both dapsone and ofloxacin (one, 5%) and rifampicin and ofloxacin (one, 5%) were less frequent than those found in Brazil [Citation29] and Korea [Citation24] but higher than those in India [Citation23] and Colombia [Citation12]. Other studies from Myanmar, Indonesia, Philippines [Citation25], Vietnam [Citation26], and South America [Citation30] reported higher dapsone AMR cases than those observed in the present study. A previous study from China reported higher ofloxacin resistance than that described in the present study [Citation27].
There was no significant correlation in AMR mutations in M. leprae in male vs. female, nerve vs. no nerve involvement, deformity vs. no deformity, and reaction vs. non-reaction. However, higher AMR mutations were recorded in male, nerve involvement, no deformity, and reaction cases of leprosy, in agreement with the findings of the other studies [Citation12,Citation16,Citation28,Citation29]. Similar to studies reported previously, no mutations in M. leprae were identified in PB, whereas all mutations were assigned to MB cases of leprosy [Citation12,Citation16]. In accordance with other studies, the TT and BB cases did not exhibit any mutations, whereas higher mutation frequency was identified in LL cases of leprosy [Citation16,Citation28,Citation29]. The reported AMR mutations in M. leprae in different clinical forms of leprosy suggest that the situation needs to be carefully monitored with an effective surveillance system.
Mutations in extended DRDRs genes of M. leprae, including rpoA, rpoB, and rpoC were observed in M. tuberculosis [Citation7], whereas only rpoB and rpoC mutations were found in M. leprae [Citation8]. We identified only one rpoC mutation in M. leprae. Mutations at 362 of gyrA and 214 of gyrB conferring fluoroquinolone AMR were the first identified mutations in the Chinese leprosy population. We also reported the 23S rRNA mutation conferring clarithromycin resistance in Chinese leprosy patients for the first time. The role of these mutations still requires further study to determine their phenotypic outcomes in M. leprae.
Considering well established national policies and socioeconomic factors, reduced transmission of both AMR and non-AMR cases of leprosy have been observed in China. Though dapsone monotherapy has been used in the pre-MDT era across China [Citation31], the irregular use of this drug resulted in high dapsone resistance. It was believed that these resistant strains are continuously transmitted in the community [Citation32]. Fortunately, rifampicin drug resistance is very low in the current study population. Patients administered with rifampicin for the treatment of other diseases and irregular MDT may explain the rifampicin resistance in China [Citation33]. Fluoroquinolones are promising and widely used antibiotics introduced into routine clinical practice to treat infectious diseases in China. Given the easy availability and inappropriate use of these drugs, AMR to fluoroquinolones has emerged in China [Citation34]. Clarithromycin is another antibiotic used to treat various bacterial infections [Citation35,Citation36], and it is used as an alternative or combinational drug that can result in clarithromycin resistance of leprosy [Citation37].
Generally, patients who remain BI zero or become negative do not develop any positive indication for AMR. However, in the current study, we observed folP1, rpoC, and gyrA mutated strains in three different cases with BI zero. These findings revealed that BI negative cases might also have active bacteria showing drug resistance. Hence, careful monitoring of BI negative cases of leprosy is necessary for drug resistance surveillance. Generally, a lack of BI decline is an indicator of drug resistance, but the present study demonstrated that drug resistance might not be the main factor for lack of BI decline. BI decline is gradual and can take several years to become negative. Therefore, examining bacterial viability by mouse footpad culture, morphological index, and molecular tests may be necessary for establishing the lack of BI decline in leprosy.
Genotype analysis can be used to trace the transmission of M. leprae strains [Citation38]. Isolates forming clear clusters of a certain strain type commonly demonstrate the recent and ongoing transmission of leprosy [Citation38]. Most of the strains of M. leprae from Guizhou, Yunnan, and Hunan exhibited more clear clusters than those from Sichuan at the province level (), implying a recent transmission and being the major drive for maintaining leprosy at the province level in China. The genotyping cluster analysis of both wild-type and mutated strains of M. leprae in leprosy cases obtained from the counties where mutated strains were detected revealed that the mutated strains had less cluster pattern with wild-type strains at the county level, indicating the existence of genotypically diverged AMR strains of M. leprae even in the same county-level (). Hence, genotypically diverged strains with less clustering patterns denote no active transmission of AMR strains of leprosy at the county level. Interestingly, one clear cluster was observed with wild-type (n = 5) and mutated strains (n = 2) from Huayuan county of Hunan. These cases were analysed for geographical landscape and migration history at the village and county levels. The patients carrying the mutated strains of this cluster were from a different village, and no migration information of the patients amongst these villages was observed. Remaining patients carrying wild-type strains were also evaluated to the mutated cases for their contact and migration information, and there was no knowledge of contact and migration information. Thus, these cluster cases might not have the same infectious source and are supposed to have acquired the infection at an earlier stage or by a divergent transmission line at the village and county levels where many generations of leprosy cases resided. China is operating very effective leprosy control programmes at the county level to eliminate leprosy, led to decreased transmission and thus further lined up the less cluster formation of strains and AMR strains of leprosy at county levels. Thus, molecular epidemiological analysis revealed low ongoing transmission of AMR strains of leprosy at the village and county levels in China.
Our AMR surveillance study provides additional information to the leprosy control programmes in China. The emergence of drug resistance in newly diagnosed cases with no previous treatment history of any drugs confirms the ongoing transmission of primary AMR strains of leprosy in China. Given the high mutation frequency in relapse cases, the surveillance of relapse cases should be highly prioritized. Based on our results, we recommend systemic screening of AMR associated mutations when prescribing either WHO-recommended MDT or alternative drug treatment for leprosy. While considering alternative drugs to dapsone and rifampicin in MDT, clarithromycin may be a promising choice as it is associated with low AMR mutation frequency of leprosy in China. We also reported that AMR mutation might not be the main factor for the lack of BI decline in leprosy. Thus, the mechanism underlying the phenomenon of lack of BI decline in leprosy should be investigated. Furthermore, molecular epidemiological analysis revealed low ongoing transmission of mutated strains at the county and village levels in China. However, interventions need to be implemented for tracing close contact of leprosy to prevent further transmission of mutated strains.
Ethics approval and consent to participate
This study was approved by the institutional ethical committee and informed consent was obtained from all participants.
Supplemental Material
Download MS Word (2.5 MB)Acknowledgment
We thank Jianping Shen, Guocheng Zhang, Ning Yong, Liu Yangying and all other technical and non-technical staff of Yunnan, Guizhou, Sichuan and Hunan provinces for their support in sample and data collection.
Disclosure statement
No potential conflict of interest was reported by the authors.
Additional information
Funding
Related Research Data
References
- Scollard DM, Adams LB, Gillis TP, et al. The continuing challenges of leprosy. Clin Microbiol Rev. 2006;19(2):338–3381.
- Sun PW, Wang L, Yu MW, et al. Leprosy statistics in China, 2017. Int J Dermatol Venereol. 2019;2(1):1–5.
- Cambau E, Saunderson P, Matsuoka M, et al. Antimicrobial resistance in leprosy: results of the first prospective open survey conducted by a WHO surveillance network for the period 2009–15. Clin Microbiol Infect. 2018;24(12):1305–1310.
- Dallas WS, Gowen JE, Ray PH, et al. Cloning, sequencing, and enhanced expression of the dihydropteroate synthase gene of Escherichia coli MC4100. J Bacteriol. 1992;174(18):5961–5970.
- Ebenezer GJ, Norman G, Joseph GA, et al. Drug resistant-Mycobacterium leprae–results of mouse footpad studies from a laboratory in south India. Indian J Lepr. 2002;74(4):301–312.
- Cambau E, Perani E, Guillemin I, et al. Multidrug-resistance to dapsone, rifampicin, and ofloxacin in Mycobacterium leprae. Lancet. 1997;349(9045):103–104.
- Comas I, Borrell S, Roetzer A, et al. Whole-genome sequencing of rifampicin-resistant Mycobacterium tuberculosis strains identifies compensatory mutations in RNA polymerase genes. Nat Genet. 2011;44(1):106–110.
- Benjak A, Avanzi C, Singh P, et al. Phylogenomics and antimicrobial resistance of the leprosy bacillus Mycobacterium leprae. Nat Commun. 2018;9(1):352.
- Malik S, Willby M, Sikes D, et al. New insights into fluoroquinolone resistance in Mycobacterium tuberculosis: functional genetic analysis of gyrA and gyrB mutations. PLoS One. 2012;7:e39754.
- Shen J, Liu M, Zhang J, et al. Relapse in MB leprosy patients treated with 24 months of MDT in south west China: a short report. Lepr Rev. 2006;77(3):219–224.
- Guerrero MI, Muvdi-Arenas S, León-Franco CI. Relapses in multibacillary leprosy patients: a retrospective cohort of 11 years in Colombia. Lepr Rev. 2012;83(3):247–260.
- Beltrán-Alzate C, López Díaz F, Romero-Montoya M, et al. Leprosy drug resistance surveillance in Colombia: The Experience of a Sentinel Country. PLoS Negl Trop Dis. 2016;10(10):e0005041.
- Job CK, Jayakumar J, Aschhoff M, et al. Viability of Mycobacterium leprae in skin and peripheral nerves and persistence of nerve destruction in multibacillary patients after 2 years of multidrug therapy. Int J Lepr Other Mycobact Dis. 1996;64(1):44–50.
- Hasanoor Reja AH, Biswas N, Biswas S, et al. Report of rpoB mutation in clinically suspected cases of drug resistant leprosy: a study from Eastern India. Indian J Dermatol Venereol Leprol. 2015;81(2):155–261.
- Sekar B, Arunagiri K, Kumar BN, et al. Detection of mutations in folp1, rpoB and gyrA genes of M. leprae by PCR- direct sequencing–a rapid tool for screening drug resistance in leprosy. Lepr Rev. 2011;82(1):36–45.
- Lavania M, Jadhav RS, Chaitanya VS, et al. Drug resistance patterns in Mycobacterium leprae isolates from relapsed leprosy patients attending The leprosy Mission (TLM) Hospitals in India. Lepr Rev. 2014;85(3):177–185.
- Matsuoka M, Suzuki Y, Garcia IE, et al. Possible mode of emergence for drug-resistant leprosy is revealed by an analysis of samples from Mexico. Jpn J Infect Dis. 2010;63(6):412–416.
- Sapkota BR, Ranjit C, Neupane KD, et al. Development and evaluation of a novel multiple-primer PCR amplification refractory mutation system for the rapid detection of mutations conferring rifampicin resistance in codon 425 of the rpoB gene of Mycobacterium leprae. J Med Microbiol. 2008;57(2):179–184.
- Balagon MF, Cellona RV, Cruz E, et al. Long-term relapse risk of multibacillary leprosy after completion of 2 years of multiple drug therapy (WHO-MDT) in Cebu, Philippines. Am J Trop Med Hyg. 2009;81(5):895–899.
- Weng X, Wang Z, Liu J, et al. Identification and distribution of Mycobacterium leprae genotypes in a region of high leprosy prevalence in China: a 3-year molecular epidemiological study. J Clin Microbiol. 2007;45(6):1728–1734.
- Jadhav RS, Kamble RR, Shinde VS, et al. Use of reverse transcription polymerase chain reaction for the detection of Mycobacterium leprae in the slit-skin smears of leprosy patients. Indian J Lepr. 2005;77(2):116–127.
- Kimura M, Sakamuri RM, Groathouse NA, et al. Rapid variable-number tandem-repeat genotyping for Mycobacterium leprae clinical specimens. J Clin Microbiol. 2009;47(6):1757–1766.
- Lavania M, Singh I, Turankar RP, et al. Molecular detection of multidrug-resistant Mycobacterium leprae from Indian leprosy patients. J Glob Antimicrob Resist. 2018;12:214–219.
- You EY, Kang TJ, Kim SK, et al. Mutations in genes related to drug resistance in Mycobacterium leprae isolates from leprosy patients in Korea. J Infect. 2005;50(1):6–11.
- Matsuoka M, Budiawan T, Aye KS, et al. The frequency of drug resistance mutations in Mycobacterium leprae isolates in untreated and relapsed leprosy patients from Myanmar, Indonesia and the Philippines. Lepr Rev. 2007;78(4):343–352.
- Kai M, Nguyen Phuc NH, Nguyen HA, et al. Analysis of drug-resistant strains of Mycobacterium leprae in an endemic area of Vietnam. Clin Infect Dis. 2011;52(5):127–132.
- Liu D, Zhang Q, Sun Y, et al. Drug resistance in Mycobacterium leprae from patients with leprosy in China. Clin Exp Dermatol. 2015;40(8):908–911.
- Contreras Mejía Mdel C, Porto Dos Santos M, Villarouco da Silva GA, et al. Identification of primary drug resistance to rifampin in Mycobacterium leprae strains from leprosy patients in Amazonas State, Brazil. J Clin Microbiol. 2014;52(12):4359–4360.
- Da Silva Rocha A, Cunha M, Diniz LM, et al. Drug and multidrug resistance among Mycobacterium leprae isolates from Brazilian relapsed leprosy patients. J Clin Microbiol. 2012;50(6):1912–1917.
- Singh P, Busso P, Paniz-Mondolfi A, et al. Molecular drug susceptibility testing and genotyping of Mycobacterium leprae strains from South America. Antimicrob Agents Chemother. 2011;56(6):2971–2983.
- Shen JP, Gupte MD, Jiang C, et al. Trends of case detection and other indicators of leprosy in China during 1985–2002. Chin Med Sci J. 2005;20(2):77–82.
- Jing Z, Zhang R, Zhou D, et al. Twenty-five years follow up of MB leprosy patients retreated with a modified MDT regimen after a full course of dapsone mono-therapy. Lepr Rev. 2009;80(2):170–186.
- Shen J, Yan L, Sun P. Clinical features of relapse after multidrug therapy for leprosy in China. Lepr Rev. 2015;86(2):165–179.
- Xiao YH, Wang J, Li Y. Bacterial resistance surveillance in China: a report from Mohnarin 2004–2005. Eur J Clin Microbiol Infect Dis. 2008;27:697–708.
- Fontana C1, Favaro M, Minelli S, et al. New site of modification of 23S rRNA associated with clarithromycin resistance of Helicobacter pylori clinical isolates. Antimicrob Agents Chemother. 2002;46(2):3765–3769.
- De Francesco V1, Zullo A, Giorgio F, et al. Change of point mutations in Helicobacter pylori rRNA associated with clarithromycin resistance in Italy. J Med Microbiol. 2014;63(3):453–467.
- Gunawan H, Sasmojo M, Putri HE, et al. Clinical pilot study: clarithromycin efficacy in multibacillary leprosy therapy. Int J Mycobacteriol. 2018;7(2):152–155.
- Smith CM, Maguire H, Anderson C, et al. Multiple large clusters of tuberculosis in London: a cross-sectional analysis of molecular and spatial data. ERJ Open Res. 2017;3(1):00098–02016.