- Research article
- Open Access
A deletion in the RD105 region confers resistance to multiple drugs in Mycobacterium tuberculosis
BMC Biologyvolume 17, Article number: 7 (2019)
The emergence of drug-resistant strains of Mycobacterium tuberculosis (Mtb), especially those that are multidrug resistant poses a serious threat to global tuberculosis control. However, the mechanism underlying the occurrence of drug resistance against more than one drug is poorly understood. Given that the Beijing/W strains are associated with outbreaks and multidrug resistance, they may harbor a genetic advantage and provide useful insight into the disease. One marker found in all Beijing/W Mtb strains is a deletion of RD105 region that results in a gene fusion, Rv0071/74, with a variable number (3–9 m) of VDP (V: Val, D: Asp; P: Pro) repeats (coded by gtggacccg repeat sequences) at the N-terminal. Here, we report that this variable number of VDP repeats in Rv0071/74 regulates the development of multidrug resistance.
We collected and analyzed 1255 Beijing/W clinical strains. The results showed that the number of VDP repeats in Rv0071/74 was related to the development of multidrug resistance, and the deletion of Rv0071/74-9 m from Beijing/W clinical strain restored drug susceptibility. Rv0071/74-9 m also increased resistance to multiple drugs when transferred to different mycobacterial strains. Cell-free assays indicate that the domain carrying 4–9 VDP repeats (4–9 m) showed a variable binding affinity with peptidoglycan and Rv0071/74 cleaves peptidoglycan. Furthermore, Rv0071/74-9 m increased cell wall thickness and reduced the intracellular concentration of antibiotics.
These findings not only identify Rv0071/74 with VDP repeats as a newly identified multidrug resistance gene but also provide a new model for the development of multiple drug resistance.
Tuberculosis is a global health concern. Almost one third of individuals worldwide are infected with the pathogen Mycobacterium tuberculosis (Mtb) and are at risk of developing tuberculosis disease during their lifetime . The emergence of drug-resistant strains of Mtb, especially those that are multidrug resistant (MDR) and extensively drug resistant (XDR), has posed a serious threat to global tuberculosis control. Given the alarming rise of resistance to tuberculosis drugs worldwide, the identification of resistance gene is critical for the future of tuberculosis control. Several gene mutations in specific loci of the Mtb genome have been reported as the basis for drug resistance and as drug targets for the development of anti-tuberculosis drugs [2,3,4,5]. However, these known resistance genes are mainly for single drug. There is little evidence that single-locus mutations confer resistance to multiple drugs.
A genetically related group of Mtb strains called Beijing/W is widespread in many regions of the world [6,7,8,9]. Strains from this group have been reported to be associated with drug resistance [10, 11]. A deletion of RD105 region is found in all Beijing/W Mtb strains and thus serves as a marker of this group [12, 13]. The RD105 region is 3467 bp in size and includes a truncated C-terminal of Rv0071, the full-length Rv0072, full-length Rv0073, and a truncated N-terminal Rv0074 (Fig. 1a). Deletion of the RD105 in Beijing/W Mtb strain generates a fusion gene Rv0071/74 that contains a 1–84 bp region of Rv0071 and a 288–1236 bp region of Rv0074 (Additional file 1: Figure S1). The 1–84 bp region of Rv0071 includes variable 9-bp sequence (gtggacccg, coding VDP) repeats plus a WxL domain [14,15,16,17] ahead and the 288–1236 bp region of Rv0074, which encodes an amidohydrolase-like family protein . However, the genetic advantage of deletion of the RD105 in Beijing/W strains is unknown. So, we aimed to define the association of RD105 deletion with drug resistance in Mtb.
Fusion gene Rv0071/74 confers resistance to multiple drugs
We analyzed 1508 Mtb clinical strains, from which 1255 of these strains were found to be Beijing/W strains as detected by deletion-targeted multiplex PCR (DTM-PCR) (Additional file 1: Figure S2). In Beijing/W strains, a variable number ranging from 3 to 9 of 9 bp repeats were found at the N-terminal of the Rv0071/74 fusion gene (Fig. 1 and Additional file 1: Figure S1A). Strains carrying 4–9 of 9 bp repeats (4–9 m) showed a different drug resistance, but no detectable drug resistance was observed in the tested Beijing/W clinical strains carrying the 3 m allele (Table 1 and Fig. 2). To understand the correlation between number of tandem repeats and drug-sensitive phenotypes, the phylogeny tree of variable 9-bp sequence (gtggacccg,) repeats was constructed (Fig. 2). Two major groups were clustered according to genotype: 3 m with no drug resistance in the tested strains was divided into a separate group, other 6 alleles (4 m, 5 m, 6 m, 7 m, 8 m, and 9 m) with drug resistance against multiple drugs were grouped together. The closest genetic distance was observed between 8 m (77.63%) and 9 m (100%) with high degree of drug resistance in the tested strains (Fig. 2). Except 4 m allele, the percentage of drug-resistant strains increases as the number of tandem repeats increases. Especially, all the Beijing/W strains carrying the Rv0071/74 fusion gene with 9 of 9 bp repeats (Rv0071/74-9 m) showed a high resistance against multiple drugs, including streptomycin (SM), rifampicin (RFP), isoniazid (isonicotinic acid hydrazide, INH), ethambutol (EMB), ofloxacin (OFX), and amikacin (AMK) (Table 1 and Fig. 2).
To test whether the presence of 9 bp repeats in the 5′ region is linked to the gene expression, we analyzed the mRNA level of Rv0071/74 from Beijing/W clinical strains carrying various 9 bp repeats. Reference strain H37Rv (wild Rv0074) was used as a control. No significant difference was found for the transcriptional level of Rv0071/74 (Additional file 1: Figure S3), suggesting that the number of 9 bp repeats is not associated with the expression of Rv0071/74.
We next examined whether the observed multiple drug resistance of Rv0071/74-9 m is caused by known drug-resistance gene mutations. All Beijing/W strains with Rv0071/74-9 m were selected to be tested. Region of drug-resistant gene for SM (rpsL and rrs1), RFP (rpoB), INH (katG and inhA), OFL (gyrA), and AMK (rrs2) were amplified by conventional PCR, and the point mutations of these genes were analyzed by sequencing assay [19, 20]. We found that some classical drug-resistance-related mutations are absent in the Rv0071/74-9 m strains and yet the strains still showed resistance to multiple drugs (Table 2). These results suggest that the single point mutation of known drug resistance genes may not be responsible for multiple drug resistances of Rv0071/74-9 m.
To further determine whether the Rv0071/74-9 m fusion gene confer the antibiotic resistance in mycobacteria, we cloned the Rv0071/74 fusion gene with 3 m or 9 m into plasmid pVV16 and then transformed these constructs into different mycobacterial strains including M. smegmatis (ATCC 19420), M. marinum (ATCC 927), and M. tuberculosis H37Ra, which lack the Rv0071/74 fusion gene. No significant difference was observed for the growth rate of these recombinant strains (data not shown), but recombinant mycobacterial strains carrying Rv0071/74-9 m showed markedly increased resistance to multiple drugs (Table 3). In contrast, recombinant strains expressing the Rv0071/74 with 3 m or Rv0071/Rv0074 alone showed no significant gain of drug resistance (Additional file 2: Table S3). Lastly, we generated Rv0071/74-9 m-null mutant by the phage transduction method (Additional file 1: Figure S4A-C). Deletion of Rv0071/74-9 m from Beijing/W strains caused no significant difference for growth rate (data not shown), but markedly reduced resistance to multiple anti-tuberculosis drugs (Table 4). The restoration of Rv0071/74-9 m allele in Rv0071/74-9 m-null mutant (△Rv0071/74-9 m:: Rv0071/74-9 m) restored the mutant to resistance to the tested drugs (Table 4). Together, these results suggest that the Rv0071/74 fusion gene with 9 m (Rv0071/74-9 m) may function as a multidrug resistance gene.
Fusion protein Rv0071/74 binds with and cleaves peptidoglycan
Quaternary structure modeling analysis of the Rv0071/74 fusion protein revealed that it has a high homology to a carboxypeptidase from Caulobacter crescentus (Additional file 1: Figure S5A). The carboxypeptidase is associated with cell membranes and cleaves the C-terminal d-alanine residue of UDP-muramyl-pentapeptide of peptidoglycan (PGN) in C. crescentus [21, 22]. Indeed, purified Rv0071/74-9 m protein can directly cleave the PGN isolated from S. aureus or M. smegmatis, as measured by OD450 or liquid chromatography (Fig. 3a–c). Furthermore, carboxypeptidase Cc2672 contains two active-site amino acid residues that are located in β-sheet regions of the protein (Additional file 1: Figure S5B). Structure modeling analysis of the Rv0071/74 protein also revealed two conserved active-site amino acid residues (AS1 and AS2) and analogous active regions similar to those of the carboxypeptidase Cc2672 (Additional file 1: Figure S5B). Point mutations in these two active-site amino acid residues all change the modeled quaternary structure of Rv0071/74 (Additional file 1: Figure S5B). We mutated these two active-site amino acid residues individually in Rv0071/74-9 m. The mutant protein was expressed at the same level as the non-mutant protein in M. smegmatis (data not shown). The purified mutant protein Rv0071/74 lost the carboxypeptidase (Fig. 3a, c). Recombinant M. smegmatis strains carrying these mutated Rv0071/74-9 m alleles lost drug resistance as compared with the strains containing pVV16::9 m (Additional file 2: Table S4).
The Rv0071/74 fusion protein inherits most (including microsatellites) of the WxL domain, which is responsible for the binding to PGN [14,15,16,17] (Fig. 1). To evaluate the association of Rv0071/74-9 m with PGN, we synthesized variable-VDP repeats-associated WxL peptides, MSSITVS-(VDP)n-VDAVVAVGRARR (n = 3–9), and labeled with biotin at the N-terminal. The biotin-labeled peptide was then incubated with PGN from S. aureus and analyzed by a ForteBio’s Octet platform. The WxL domain with 9 VDP repeats showed much higher affinity with PGN as compared to WxL-3 m (Fig. 3d). In addition, WxL domain carrying 4–8 VDP repeats (4-8 m) showed a variable binding affinity with PGN (Additional file 1: Figure S6), which may explain the variable degree of multiple drug resistance observed in strains carrying 4–8 of 9 bp repeats (4–8 m).
Rv0071/74-9 m increases cell wall thickness and reduces intracellular drug concentration
Transmission electron microscopy (TEM) was used to evaluate the differences between M. smegmatis containing pVV16::9 m and control groups containing wild-type M. smegmatis and M. smegmatis with only the pVV16 plasmid or pVV16::3 m. Under TEM, marked differences were observed in the thickness of cell walls: 37.10 ± 1.21 nm for M. smegmatis containing pVV16::9 m, but < 15 nm for the control groups (Fig. 4a, b).
To examine whether the increased thickness of the cell wall in M. smegmatis::Rv0071/74-9 m affects the permeability of drugs, the intracellular concentration of streptomycin or amikacin in recombinant M. smegmatis::Rv0071/74-9 m was analyzed by mass spectrometry. M. smegmatis transformed with pVV16::9 m had a much lower intracellular drug concentration of streptomycin and amikacin than the control M smegmatis carrying the recombinant pVV16::3 m plasmid (Fig. 4c).
To evaluate the differences between Beijing/W Mtb with 9 m and control groups containing H37Rv (wild-type Rv0074) and Beijing/W Mtb with Rv0071/74-3 m or 7 m, cell wall thickness of these different Mtb strains were also analyzed by TEM. The cell walls of Beijing/W Mtb strains with 9 m are much thicker than those of control groups under TEM (p < 0.0001) (Fig. 4d). Deletion of Rv0071/74-9 m from Beijing/W strains significantly decreased the cell wall thickness (p < 0.0001) (Fig. 4d) and then markedly reduced resistance to multiple anti-tuberculosis drugs (Table 4).
Multidrug-resistant Mtb poses the greatest obstacle to the tuberculosis control. A better knowledge of the mechanisms of drug resistance in Mtb and the molecular mechanisms involved will stimulate the exploration of new targets for drug activity and drug development. So far, spontaneous point mutations in the specific gene have been characterized as the main basis for Mtb drug resistance, especially to RIF and INH [2, 5]. More than 95% RIF resistance is mainly caused by mutations of 81 bp RRDR region (resistance-determining region) of rpoB gene [3, 5]. And the S315 T mutation in the katG gene is the predominant mechanism of INH resistance in Mtb [3, 5]. Our data also showed this phenomenon (Table 2). The resistance to EMB, SM, OFL, and AMK are primarily found in mutations in the embB, rpsl, gyrA, and rrs genes, respectively. In summary, mutations in these genes are associated with drug resistance in roughly 70–85% of drug resistance isolates of Mtb [3, 5]. At present, new genetic mutations of drug resistance have been identified by structural analysis of Mtb pan-genome . However, these new resistance genes still require more detailed experimental studies to test their association with resistant phenotypes.
Although point mutations in resistant genes play a major role in Mtb drug resistance, it has been reported that other drug-resistant mechanisms such as drug efflux pumps, changes in cell wall permeability, etc. also play an important role in drug tolerance. The development of drug resistance in mycobacteria is complex and controlled by a variety of mechanisms. Multiple mechanisms of resistance are acting in unison in the emergence of drug resistance. Our study here showed a newly identified multidrug resistance gene Rv0071/74, which is generated from the deletion of the RD105 region in Beijing/W strains.
We have found that the Beijing/W strains carrying 3–9 of VDP repeats (4–9 m) showed a different degree of multiple drug resistance. As the number of repeats increases, the percentage of resistant strains also increases. Based on the quaternary structure modeling of the Rv0071/74 fusion protein, Rv0071/74 (101-214aa) has a high homology to a carboxypeptidase Cc2672 from Caulobacter crescentus Cb15 (PDB ID: 3MTW) Ad3, 152–269 aa). Furthermore, similar to carboxypeptidase Cc2672, Rv0071/74 protein also revealed two conserved active-site amino acid residues (AS1 and AS2) that are located in β-sheet regions of the protein. Purified recombinant Rv0071/74-9 m or wild-type Rv0074 from non-Beijing strains all can directly cleave the peptidoglycan (PGN) isolated from S. aureus or M. smegmati in vitro, suggesting there is no difference in the carboxypeptidase activity of these proteins (P = 0.426). However, compared with wild-type Rv0074 from non-Beijing strain, fusion Rv0071/74 gene generated from the deletion of the RD105 region in Beijing/W strain carries a WxL domain plus a variable number of VDP repeats. WxL domain is known for its ability to interact with PGN [14,15,16,17]. Indeed, WxL domain carrying 4–9 VDP repeats (4–9 m) showed a variable binding affinity with PGN. Thus, the variable degree of multiple drug resistance observed in strains carrying variable number of 9 VDP repeats (3–9 m) could be due to a variable binding affinity WxL domain carrying 3–9 VDP repeats with PGN in vivo.
PGN, an essential component of the bacterial cell envelope, provides structural integrity to mycobacterial cells [24, 25]. Since Rv0071/74-9 m can bind and cleave the PGN, we have evaluated its effect on the morphology characteristics of mycobacteria at the ultrastructural level. Our data showed that the thickness of cell walls of Rv0071/74-9 m strains is markedly increased as compared to those of the control group. Thickening of the bacterial cell wall generally increases drug resistance by reducing drug permeability [26,27,28,29]. Our results also demonstrated that those mycobacterial strains transformed with pVV16::9 m have a much thicker cell wall, and much lower intracellular drug concentration than the control. The previous studies have revealed that the cell walls in resistant (MDR and XDR) strains are much thicker than those of susceptible strains of Mtb under TEM [30, 31]. Moreover, the basal layer of the PG in XDR strains appeared more electron dense and was almost diffused into the intermediate layer, forming the dense thick layer around the plasma membrane . This phenomenon was also observed in M. smegmatis containing pVV16::9 m (Fig. 4a) and Beijing/W Mtb with 9 m under TEM.
Thus, the resistance of Rv0071/74-9 m or to multiple drugs could be due to the reduced permeability of the cell wall to multiple drugs. However, the exact mechanism of how the carboxypeptidase activity of the Rv0071/74 fusion protein change antibiotic permeability is still await for further investigation, one possible explanation is that Rv0071/74-9 m fusion protein may bind to PGN through its WxL domain and cleave PGN through its carboxypeptidase activity. This event may enhance the cross-linking of PGN to increase the thickness of the cell wall in mycobacteria, thus reducing the permeability of the cell wall to multiple drugs.
Tuberculosis remains one of the leading causes of death worldwide. The emergence of drug-resistant strains of Mtb, especially those that are multidrug resistant and extensively drug resistant, has posed a serious threat to global tuberculosis control. Given the alarming rise of resistance to tuberculosis drugs worldwide, the identification of resistance gene is critical for the future of tuberculosis control. However, these known resistance genes are mainly for single drug and there is little evidence that single-locus mutations confer resistance to multiple drugs. Our data suggest that a fused Rv0071/74 gene with 9 of 9-bp repeats (Rv0071/74-9 m) increased resistance to multiple drugs by increasing cell wall thickness and reducing the intracellular concentration of antibiotics. Our findings not only reveal a mechanism for the development of multi-drug resistance but also may provide the basis for the design of more effective anti-tuberculosis drugs.
Strain collection and preparation
A total of 1508 M. tuberculosis strains simultaneously identified by growth characteristics, colony morphology, and growth in the presence of PNB (pnitrobenzoic acid) and TCH (thiophene-2-carboxylic acid hydrazide), and 16SrRNA sequencing  were mainly collected mainly from regions of Eastern China, including Shanghai, Jiangsu, Zhejiang, Shandong, Fujian, Anhui, and Jiangxi, by Laboratory of Tuberculosis, Shanghai Pulmonary Hospital. All strains were grown in Sauton culture medium supplemented with 0.5 g/L KH2PO4, 0.5 g/L MgSO4·7H2O, 2 g/L citric acid, 0.05 g/L ferric ammonium citrate, 4.0 g/L L-asparagine, 6% glycerol and 0.02% Tween 80. Strains were sterilized at 80 °C for 30 min, and collected by centrifugation (12,000g for 5 min). The bacterial pellet was washed three times with sterilized saline and re-collected by centrifugation (12,000g for 10 min each time). The ethics approvals were obtained for this study from Shanghai Pulmonary Hospital Ethics Committee (the permit numbers: 2011-fk-03).
Extraction of bacteria DNA
The total bacterial pellet were resuspended in 50 mL DNA lysis buffer including 10 mmol/L NaCl, 1 mg/mL SDS, 15% Chelex-100, and 1% Tween 20. The mixture was incubated at 50 °C for 1 h, followed by 100 °C for 10 min, then centrifuged (5000g for 10 min) to obtain the aqueous phase containing genomic DNA.
Identification of Beijing strains and microsatellite typing
Beijing strains were identified by deletion-targeted multiplex PCR (DTM-PCR) to detect the genomic deletion RD105, which defines the Beijing family as a separate lineage within Mtb. The DTM-PCR primers BP1 (5′-GGAGTCGTTGAGGGTGTTCATCAGCTCAGTC-3′) and BP2 (5′-CGCCAAGGCCGCATAGTCACGGTCG-3′) were used to amplify a 1466 bp product from the non-W-Beijing strains, while BP1 and BP3 (5′-GGTTGCCCACTGGTCGATATGGTGGACTT-3′) were used to amplify a 761-bp fragment from the Beijing genotype [33, 34]. PCR products were separated on 0.8% agarose gels. To detect polymorphisms of the 9-bp sequence repeats locus within Beijing strains, PCR products from Beijing strains were also sequenced on an ABI 3730xl DNA Analyzer according to the procedures by manufacturers (Sangon, Shanghai, China).
Drug susceptibility test in vitro
The drug susceptibility test of all Beijing strains against the four first-line drugs (SM, INH, RFP, EMB) were carried out by the BACTEC 960 method [35, 36]. The MIC (minimum inhibitory concentration) determination of selected strains was tested in 96-well culture plates containing Middlebrook 7H9 liquid medium supplemented with 0.25% glycerol, 10% oleic acid–albumin-dextrose-catalase, and 0.05% Tween-80, as described by Kumar et al.  .The final concentrations of drugs are provided (Additional file 2: Table S1).
Mutations of drug-resistant genes
Mutation of drug-resistant genes was tested as described previously described [19, 20]. Region of drug-resistant gene for SM (rpsL and rrs1), RFP (rpoB), INH (katG and inhA), OFL (gyrA), and AMK (rrs2) were amplified by conventional PCR. PCR products were sequenced (Sangon, Shanghai, China). Point mutations of drug-resistant genes were analyzed by DNAstar soft. The region and mutations of drug-resistant gene are provided (Additional file 3: Table S2).
Plasmid construction and bacteria transformation
The Rv0071/74 fusion genes with different alleles were cloned into plasmid pVV16  and then transformed into M. smegmatis (ATCC 19420), M. marinum (ATCC927), and H37Ra. The Rv0071/74-9 m were amplified and cloned into the pET28a expressing plasmid and then transformed into E. coli BL21 (DE3) strains.
M. smegmatis and selected recombinant strains were cultured in 96 U well culture plates containing Middlebrook 7H9 liquid medium or on Middlebrook 7H10 agar medium with hygromycin or without hygromycin at 37 °C and were grown until the formation of colonies (about 3–5 days for liquid culture, and 5–7 days for solid culture). Morphological characteristics of tested strains were observed by × 2 magnifier for liquid culture and by an inverted light microscope for solid culture. For ultrastructural characteristics, strains at mid-log phase were collected and analyzed by Tecnai transmission electron microscopy (TEM) with 160 kV according the procedures by manufactures (GOODBIO, Wuhan, China). To analyze ultrastructural morphology of Mtb strains with different alleles, 15 Mtb strains including H37Rv (wild Rv0074) and Beijing/W Mtb with 9 m (n = 5), 7 m (n = 5, MDR strains), 3 m/4 m (n = 4, drug-sensitive strains), and △Rv0071/74-9 m strains were also examined by TEM observation. These tested strains were all grown in Middlebrook 7H9 Broth (Difco) supplemented with 0.2% glycerol and 10% Middlebrook OADC Enrichment (Difco), and then harvested in the exponential growth phase. 107–108 bacterial suspension was used for TEM examination. Two-tailed unpaired t tests were used to analyze the difference of cell wall thickness between recombinant M. smegmatis strains or Mtb strains. Statistical significance was defined as P < 0.05.
Generation of Rv0071/74-9 m mutants
The Rv0071/74-9 m-null mutant was generated by the phage transduction method . To construct a transducing phage for Mtb Rv0071/74-9 m knockout, the left homolog arm was PCR amplification using primers KOP1 and KOP2. The right homolog arm was PCR amplification using primers KOP3 and KOP4. The PCR products were ligated into AflII/XbaI and HindIII/BglII sites of pYUB854. The recombinant transducing phage was used to construct the Rv0071/74-9 m-null MTB as described. The mutant was verified by PCR(P1 and P2) and reverse transcription-PCR (P3 and P4). The Rv0071/74-9 m fusion gene was cloned into plasmid pVV16  and then transformed into the Rv0071/74-9 m-null mutant.
RNA extraction and quantitative real-time PCR
MTB strains (OD600 0.25) at mid-log phase were collected. Total RNA was extracted with TRIzol-Reagent (Invitrogen). Briefly, cell pellet was resuspended in 1 ml Trizol reagent, mixed with 400 ml 0.1 mm Zirconia Beads (Sigma Products) and lysed in a mini-bead beater (Biospec) for three cycles (40 s at maximal speed) with cooling on ice for 1 min between pulses. RNA was extracted according to the protocol of TRIzol-Reagent. The extracted RNA was further followed by DNase I treatment to eliminate DNA contamination. cDNA was synthesized using the PrimerScript II 1st strand cDNA synthesis kit (TaKaRa). Target gene transcript levels were measured by real-time PCR using SYBRH Premix Ex Taq GC (TaKaRa) on applied biosystems 7500 real-time PCR: 95 °C 60 s, 40 cycles of 95 °C 5 s, 62 °C 8 s and 72 °C 20 s, followed by melting curve analysis. The relative transcriptional level was determined by the 2−ΔΔCt method. Reference strain H37Rv (wild-type Rv0074) was as a control. The reference gene was 16S rRNA.
Protein structure analysis
Quaternary structure modeling of Rv0071/74-9 m fusion protein was analyzed online by SWISS-MODEL at http://swissmodel.expasy.org/ . Active regions and active sites were analyzed online in Conserved Domain Database of NCBI (http://www.ncbi.nlm.nih.gov/).
Binding analysis of WxL domain with peptidoglycan
The binding between WxL domain and peptidoglycan (PGN) was measured using interferometry by a ForteBio’s Octet platform. The WxL daomains from Rv0071/74 were synthesized as peptides with the following sequence MSSITVS-(VDP)n-VDAVVAVGRARR (n = 3~9), and biotin-labeled at the N terminal. Streptavidin-labeled biosensor was firstly equilibrated incubated in PBS buffer for 30 s, and then loaded with the biotin-labeled peptide (0.5 g/L) in the PBS buffer for 3 min. After the baseline step in PBS buffer for 30 s, the biosensors with the biotin-labeled peptide was placed in the PGN from S. aureus (Sigma-Aldrich) in PBS buffer (1 g/ L) for 5 min for association, and then placed in the PBS buffer for 5 min for dissociation. The binding of WxL region to PGN was analyzed with subtracted data (obtained by subtracting background data from raw data) using the Octet data analysis software.
Peptidoglycan cleavage assay
Protein was purified to nearly SDS-PAGE homogeneity. Protein concentration was determined by Coomassie brilliant blue G-250 method. The purified Rv0071/74-9 m protein (10 μg/well) was incubated with PGN from S. aureus at 37 °C for 0, 1, 2, 3, 4, 5, and 6 h. OD450 of each time point was tested at a spectrophotometer (BIOTEK uQuant MQX200). The change of OD450 of each time point was analyzed by GraphPad Prism 5.0 statistical software package (http://www.graphpad.com/prism/prism.htm). For liquid chromatography analysis of PGN cleavage, the PGN-S (from S. aureus) or PGN-M (from M. smegmatis) was first incubated with the Rv0071/74 protein or lysozyme, a known glycanohydrolase catalyzing hydrolysis of PGN, as a positive control in our study, at 37 °C for 2 h, the samples were loaded on the filter paper and spread for 1 h, and then stained by 5% Ninhydrin.
Polyclonal antibody preparation and Western blot
Rv0071/74-9 m fusion protein was purified as previously described and used to immunize rabbits for the production of anti-Rv0071/74-9 m fusion protein polyclonal antibody. Standard western blot procedures were used. Cell extracts (50 μg) were separated by 12% SDS–PAGE and then transferred to a polyvinylidene fluoride membrane. The concentration of the primary antibodies used in the corresponding blot were anti-Rv0071/74-9 m (1:1000). Horseradish peroxidase-conjugated goat anti-rabbit polyclonal antibody was used as the secondary antibody at a 1:5000 dilution.
M. smegmatis or recombinant M. smegmatis at mid-log phase were diluted to a calculated starting OD600 of 0.25 and added into 96 well culture plates containing Middlebrook 7H9 liquid medium. OD600 of each time point of each strain was tested in real time by Bioscreen C microplate incubator (FP1100-C, Labsystems, USA) at 37 °C for 7 days. Data shown are representative of three duplicates. Error bars, means ± SD. Two-tailed unpaired t-tests were used to analyze the difference of growth rate between recombinant M. smegmatis strains. Statistical significance was defined as P < 0.05.
Analysis of intracellular concentrations of drugs
M. smegmatis or recombinant M. smegmatis at mid-log phase were diluted to a calculated starting OD600 of 0.25 and grown with 1/4 MIC drugs for 48 h at 37 °C. Strains were collected and washed three times with sterilized ddH2O. Bacterial pellets were resuspended in 200 μL ddH2O. The supernatant containing intracellular drugs was prepared via ultrasonication and loaded on a weak cation-exchange (WCX) solid-phase extraction column. The hydrophilic interaction chromatography (HILIC) was used to retain the drugs on the column for the separation of the analytes. The drugs were detected and quantified with triple quadrupole tandem mass spectrometry with ESI source in the positive ion mode. Quantitation was performed using multiple reaction monitoring (MRM) of the transitions with isepamycin as the internal standard. Data shown are representative of three independent experiments. Error bars, means ± SD. Two-tailed unpaired t-tests were used to analyze the difference of intracellular concentration of drugs between recombinant M. smegmatis strains. Statistical significance was defined as P < 0.05.
Deletion-targeted multiplex PCR
Hydrophilic interaction chromatography
Minimum inhibitory concentration
Multiple reaction monitoring
Sodium dodecyl sulfate polyacrylamide gel electrophoresis
Thiophene-2-carboxylic acid hydrazide
Transmission electron microscopy
Extensively drug resistant
World Health Organization. WHO global tuberculosis report 2017. 2017. Available at: http://www.who.int/tb/publications/global_report/en/.
Walsh C. Molecular mechanisms that confer antibacterial drug resistance. Nature. 2000;406:775–81.
Almeida Da Silva PE, Palomino JC. Molecular basis and mechanisms of drug resistance in Mycobacterium tuberculosis: classical and new drugs. J Antimicrob Chemother. 2011;66:1417–30.
Chang HH, Cohen T, Grad YH, Hanage WP, O'Brien TF, Lipsitch M. Origin and proliferation of multiple-drug resistance in bacterial pathogens. Microbiol Mol Biol Rev. 2015;79:101–16.
Zhang Y, Yew WW. Mechanisms of drug resistance in Mycobacterium tuberculosis: update 2015. Int J Tuberc Lung Dis. 2015;19:1276–89.
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.
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.
Liu Q, Luo T, Dong X, Sun G, Liu Z, Gan M, et al. Genetic features of Mycobacterium tuberculosis modern Beijing sublineage. Emerg Microbes Infect. 2016;5:e14.
Shitikov E, Kolchenko S, Mokrousov I, Bespyatykh J, Ischenko D, Ilina E, et al. Evolutionary pathway analysis and unified classification of East Asian lineage of Mycobacterium tuberculosis. Sci Rep. 2017;7:9227.
de Steenwinkel JE, ten Kate MT, de Knegt GJ, Kremer K, Aarnoutse RE, Boeree MJ, et al. Drug susceptibility of Mycobacterium tuberculosis Beijing genotype and association with MDR TB. Emerg Infect Dis. 2012;18:660–3.
Rodríguez-Castillo JG, Pino C, Niño LF, Rozo JC, Llerena-Polo C, Parra-López CA, et al. Comparative genomic analysis of Mycobacterium tuberculosis Beijing-like strains revealed specific genetic variations associated with virulence and drug resistance. Infect Genet Evol. 2017;54:314–23.
Tsolaki AG, Gagneux S, Pym AS, Goguet de la Salmoniere YO, Kreiswirth BN, Van Soolingen D, et al. Genomic deletions classify the Beijing/W strains as a distinct genetic lineage of Mycobacterium tuberculosis. J Clin Microbiol. 2005;43:3185–91.
Tsolaki AG, Hirsh AE, DeRiemer K, Enciso JA, Wong MZ, Hannan M, et al. Functional and evolutionary genomics of Mycobacterium tuberculosis: insights from genomic deletions in 100 strains. Proc Natl Acad Sci U S A. 2004;101:4865–70.
Brinster S, Furlan S, Serror P. C-terminal WxL domain mediates cell wall binding in Enterococcus faecalis and other gram-positive bacteria. J Bacteriol. 2007;189:1244–53.
Galloway-Peña JR, Liang X, Singh KV, Yadav P, Chang C, La Rosa SL, et al. The identification and functional characterization of WxL proteins from Enterococcus faecium reveal surface proteins involved in extracellular matrix interactions. J Bacteriol. 2015;197:882–92.
Siezen R, Boekhorst J, Muscariello L, Molenaar D, Renckens B, Kleerebezem M. Lactobacillus plantarum gene clusters encoding putative cell-surface protein complexes for carbohydrate utilization are conserved in specific gram-positive bacteria. BMC Genomics. 2006;7:126.
Dramsi S, Dehoux P, Lebrun M, Goossens PL, Cossart P. Identification of four new members of the internalin multigene family of Listeria monocytogenes EGD. Infect Immun. 1997;65:1615–25.
Zhang YY, Chen C, Liu J, Deng HJ, Pan AZ, Zhang LS, et al. Complete genome sequences of Mycobacterium tuberculosis strains CCDC5079 and CCDC5080, which belong to the Beijing family. J Bacteriol. 2011;193:5591–2.
Zheng RJ, Zhu CT, Guo Q, Qin LH, Wang J, Lu JM, et al. Pyrosequencing for rapid detection of tuberculosis resistance in clinical isolates and sputum samples from re-treatment pulmonary tuberculosis patients. BMC Infect Dis. 2014;14:200.
Cui ZL, Li YY, Cheng S, Yang H, Lu JM, Hu ZY, et al. Mutations in the embC-embA intergenic region contribute to Mycobacterium tuberculosis resistance to ethambutol. Antimicrob Agents Chemother. 2014;58:6837–43.
Xiang DF, Patskovsky Y, Xu C, Fedorov AA, Fedorov EV, Sisco AA, et al. Functional identification and structure determination of two novel prolidases from cog1228 in the amidohydrolase superfamily. Biochemistry. 2010;49:6791–803.
Markiewicz Z, Glauner B, Schwarz U. Murein structure and lack of DD- and LD-carboxypeptidase activities in Caulobacter crescentus. J Bacteriol. 1983;156:649–55.
Kavvas E, Catoiu E, Mih N, Yurkovich J, Seif Y, Dillon N, et al. Machine learning and structural analysis of Mycobacterium tuberculosis pan-genome identifies genetic signatures of antibiotic resistance. Nat Commun. 2018;9:4306.
Brennan PJ. Structure, function, and biogenesis of the cell wall of Mycobacterium tuberculosis. Tuberculosis (Edinb). 2003;83:91–7.
Hett EC, Rubin EJ. Bacterial growth and cell division: a mycobacterial perspective. Microbiol Mol Biol Rev. 2008;72:126–56.
Delcour AH. Outer membrane permeability and antibiotic resistance. Biochim Biophys Acta. 2009;1794:808–16.
Hiramatsu K. Vancomycin-resistant Staphylococcus aureus: a new model of antibiotic resistance. Lancet Infect Dis. 2001;1:147–55.
Kawai M, Yamada S, Ishidoshiro A, Oyamada Y, Ito H, Yamagishi J. Cell-wall thickness: possible mechanism of acriflavine resistance in meticillin-resistant Staphylococcus aureus. J Med Microbiol. 2009;58:331–6.
Yuan WC, Hu QW, Cheng H, Shang WL, Liu N, Hua ZY, et al. Cell wall thickening is associated with adaptive resistance to amikacin in methicillin-resistant Staphylococcus aureus clinical isolates. J Antimicrob Chemother. 2013;68:1089–96.
Velayati AA, Farnia P, Masjedi MR, Ibrahim TA, Tabarsi P, Haroun RZ, et al. Totally drug-resistant tuberculosis strains: evidence of adaptation at the cellular level. Eur Respir J. 2009;34:1202–3.
Velayati AA, Farnia P, Ibrahim TA, Haroun RZ, Kuan HO, Ghanavi J, et al. Differences in cell wall thickness between resistant and nonresistant strains of Mycobacterium tuberculosis: using transmission electron microscopy. Chemotherapy. 2009;55:303–7.
Glennon M, Smith T, Cormican M, Noone D, Barry T, Maher M, et al. The ribosomal intergenic spacer region: a target for the PCR based diagnosis of tuberculosis. Tuber Lung Dis. 1994;75:353–60.
Qin LH, Wang J, Zheng RJ, Lu JM, Yang H, Liu ZH, et al. Perspective on sequence evolution of microsatellite locus (CCG)n in Rv0050 gene from Mycobacterium tuberculosis. BMC Evol Biol. 2011;11:247.
Chen J, Tsolaki AG, Shen X, Jiang X, Mei J, Gao Q. Deletion-targeted multiplex PCR (DTM-PCR) for identification of Beijing/W genotypes of Mycobacterium tuberculosis. Tuberculosis (Edinb). 2007;87(5):446–9.
Bemer P, Palicova F, Rüsch-Gerdes S, Drugeon HB, Pfyffer GE. Multicenter evaluation of fully automated BACTEC Mycobacteria Growth Indicator Tube 960 system for susceptibility testing of Mycobacterium tuberculosis. J Clin Microbiol. 2002;40:150–4.
Rüsch-Gerdes S, Pfyffer GE, Casal M, Chadwick M, Siddiqi S. Multicenter laboratory validation of the BACTEC MGIT 960 technique for testing susceptibilities of Mycobacterium tuberculosis to classical second-line drugs and newer antimicrobials. J Clin Microbiol. 2006;44:688–92.
Kumar M, Khan IA, Verma V, Kalyan N, Qazi GN. Rapid, inexpensive MIC determination of Mycobacterium tuberculosis isolates by using microplate nitrate reductase assay. Diagn Microbiol Infect Dis. 2005;53:121–4.
Stover CK, de la Cruz VF, Fuerst TR, Burlein JE, Benson LA, Bennett LT, et al. New use of BCG for recombinant vaccines. Nature. 1991;351:456–60.
Bardarov S, Bardarov S Jr, Pavelka MS Jr, Sambandamurthy V, Larsen M, Tufariello J, et al. Specialized transduction: an efficient method for generating marked and unmarked targeted gene disruptions in Mycobacterium tuberculosis, M bovis BCG and M smegmatis. Microbiology. 2002;148:3007–17.
Arnold K, Bordoli L, Kopp J, Schwede T. The SWISS-MODEL workspace: a web-based environment for protein structure homology modeling. Bioinformatics. 2006;22:195–201.
This work was supported by National Natural Science Foundation of China (81470090 and 81201254) for LHQ, and by National Key Technologies Research and Development Program of China (2012ZX10003002–008) to ZYH.
Availability of data and materials
All data generated or analyzed during this study are included in this published article.
Ethics approval and consent to participate
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.
Figure S1. Sequences of fusion protein Rv0071/74 and fusion protein Rv0071/74-9 m in mycobacterial analyzed by West-blot. Figure S2. Distribution of drug-resistant strains in tested Mtb clinical strains. Figure S3. RT-PCR analysis of Rv0071/74 fusion gene mRNA expression from different Beijing/W clinical strains. Figure S4. Generation of Rv0071/74-9 m mutant. Figure S5. Structure modeling of Rv0071/74. Figure S6. Variable Binding of WxL domain with different alleles with PGN. (PDF 716 kb)
Table S1. Final concentrations of antibiotics in drug susceptibility test. Table S3. Drug susceptibility of recombinant M. smegmatis transformed with Rv0071 or Rv0074. Table S4. Drug susceptibility of recombinant M. smegmatis transformed with Rv0071/74-9 m carrying point mutation of AS1 or AS2. (DOCX 23 kb)
Table S2. The region and mutations of drug-resistant gene. (XLS 32 kb)