BioMed Research International

BioMed Research International / 2015 / Article

Research Article | Open Access

Volume 2015 |Article ID 271728 | https://doi.org/10.1155/2015/271728

Hongfei Duan, Guan Liu, Xiaobo Wang, Yuhong Fu, Qian Liang, Yuanyuan Shang, Naihui Chu, Hairong Huang, "Evaluation of the Ribosomal Protein S1 Gene (rpsA) as a Novel Biomarker for Mycobacterium Species Identification", BioMed Research International, vol. 2015, Article ID 271728, 8 pages, 2015. https://doi.org/10.1155/2015/271728

Evaluation of the Ribosomal Protein S1 Gene (rpsA) as a Novel Biomarker for Mycobacterium Species Identification

Academic Editor: Om Parkash
Received26 Jul 2014
Accepted07 Nov 2014
Published05 Apr 2015

Abstract

Objectives. To evaluate the resolution and reliability of the rpsA gene, encoding ribosomal protein S1, as a novel biomarker for mycobacteria species identification. Methods. A segment of the rpsA gene (565 bp) was amplified by PCR from 42 mycobacterial reference strains, 172 nontuberculosis mycobacteria clinical isolates, and 16 M. tuberculosis complex clinical isolates. The PCR products were sequenced and aligned by using the multiple alignment algorithm in the MegAlign package (DNASTAR) and the MEGA program. A phylogenetic tree was constructed by the neighbor-joining method. Results. Comparative sequence analysis of the rpsA gene provided the basis for species differentiation within the genus Mycobacterium. Slow- and rapid-growing groups of mycobacteria were clearly separated, and each mycobacterial species was differentiated as a distinct entity in the phylogenetic tree. The sequences discrepancy was obvious between M. kansasii and M. gastri, M. chelonae and M. abscessus, M. avium and M. intracellulare, and M. szulgai and M. malmoense, which cannot be achieved by 16S ribosomal DNA (rDNA) homologue genes comparison. 183 of the 188 (97.3%) clinical isolates, consisting of 8 mycobacterial species, were identified correctly by rpsA gene blast. Conclusions. Our study indicates that rpsA sequencing can be used effectively for mycobacteria species identification as a supplement to 16S rDNA sequence analysis.

1. Introduction

Members of the genus Mycobacterium are widespread in nature and range from harmless saprophytic species to strict pathogens that cause serious human and animal diseases. Both slow-growing mycobacteria and rapid-growing mycobacteria can cause human infections. Traditionally, taxonomy based on biochemical characteristics has been used for species determination of mycobacteria, but this approach is limited due to the overlapping biochemical and phenotypic patterns among the different mycobacterial species. Another approach using analysis of cell-wall fatty acid and mycolic acid composition is also limited by profile similarity among some emerging nontuberculosis mycobacteria (NTM) [1]. 16S rDNA homologue gene sequence comparison has been used as an important method for the mycobacterial species identification; however, ambiguous results have been obtained either due to the presence of more than one copy of the 16S rDNA gene within the genome, for example, in M. celatum and M. terrae complex [2, 3], or due to sequence homology between species [4]. Therefore, alternative phylogenetic markers which are capable of complementing 16S rRNA gene would be useful for the phylogenetic study and species identification of the genus Mycobacterium.

Ribosomal protein S1 (rpsA), which is in the 30S ribosome subunit, contains the S1 domain that has been found in a large number of RNA-associated proteins. RpsA is a vital protein involved in protein translation and the ribosome-sparing process of translation. In addition, it has been reported that RpsA is the target of pyrazinoic acid, the active form of the antituberculosis drug pyrazinamide [5]. Mycobacterium tuberculosis has a single copy of the rpsA gene in the genome [6], while the sequence homology among different mycobacteria species in which the gene has already been sequenced is between 86.7 and 100%. This suggests that rpsA may be suitable for phylogenetic study of the genus Mycobacterium. In this paper, we report an evaluation of the rpsA gene as a novel biomarker for mycobacteria species identification.

2. Materials and Methods

2.1. Mycobacterial Reference Strains and Clinical Isolates

42 type and reference strains of the genus Mycobacterium (Table 1), 172 clinical NTM isolates, and 16 M. tuberculosis complex (MTC) isolates were investigated. All the type and reference strains were purchased from the American Type Culture Collection (ATCC) and all the clinical isolates used in this study were obtained from the Clinical Database and Sample Bank of tuberculosis of Beijing, National Clinical Lab on Tuberculosis, Beijing Chest Hospital. All the clinical NTM isolates were identified to the species level by sequence alignment of at least two of the following: 16S rDNA, 16-23S rRNA gene internal transcribed spacer (ITS), and rpoB and hsp65 genes as described before [2, 79]. Following sequencing, the 172 NTM clinical strains were found to include 10 strains of M. avium, 75 strains of M. intracellulare, 23 strains of M. kansasii, 39 strains of M. abscessus, 17 strains of M. fortuitum complex, 7 strains of M. gordonae, and 1 strain of M. neoaurum.


SpeciesType strainFirst match (%)aSecond match (%)a

M.  thermoresistibile ATCC19527M.  senegalense (100%)M.  porcinum (98.7%)
M.  senegalense ATCC35796M.  thermoresistibile (100%) M.  porcinum (98.7%)
M.  porcinum ATCC33776M.  senegalense/M.  thermoresistibile (98.7%) M.  fortuitum (98.3%)
M.  fortuitum ATCC6481M.  porcinum (98.3%) M.  senegalense/M.  thermoresistibile (98.1%)
M.  phlei ATCC11758M.  porcinum (96.5%) M.  fortuitum (96.1%)
M.  gadium ATCC27726M.  porcinum/M.  fortuitum (95.7%) M.  smegmatis (95.3%)
M.  flavescens ATCC14474M.  smegmatis (96.6%) M.  vaccae (96.0%)
M.  smegmatis ATCC19420M.  vaccae/M.  parafortuitum/M.  triviale (96.6%) M.  senegalense/M.  thermoresistibile (96.3%)
M.  vaccae ATCC15483M.  diernhoferi/M.  duvalii (96.8%) M.  parafortuitum/M.  triviale (96.6%)
M.  parafortuitum ATCC19686M.  triviale (100%) M.  terrae/M.  austroafricanum (96.6%)
M.  triviale ATCC23292M.  parafortuitum (100%) M.  terrae/M.  austroafricanum (96.6%)
M.  neoaurum ATCC25795M.  diernhoferi/M.  duvalii (96.4%) M.  porcinum (95.5%)
M.  diernhoferi ATCC19340M.  duvalii (100%) M.  vaccae (96.8%)
M.  duvalii ATCC43910M.  diernhoferi (100%) M.  vaccae (96.8%)
M.  aurum ATCC23366M.  gilvum (93.9%) M.  chubuense/M.  vanbaalenii (93.0%)
M.  chubuense ATCC27278M.  gilvum (96.6%) M.  vanbaalenii (96.0%)
M.  gilvum ATCC43909M.  vanbaalenii (97.2%) M.  terrae/M.  austroafricanum (97.0%)
M.  vanbaalenii PYR-1M.  terrae/M.  austroafricanum (99.8%) M.  gilvum (97.2%)
M.  terrae ATCC15755M.  austroafricanum (100%) M.  vanbaalenii (99.8%)
M.  austroafricanum ATCC33464M.  terrae (100%) M.  vanbaalenii (99.8%)
M.  rhodesiae ATCC27024M.  terrae/M.  austroafricanum (93.9%) M.  obuense (93.7%)
M.  obuense ATCC27023M.  porcinum /M.  fortuitum (96.1%) M.  senegalense/M.  thermoresistibile (96.5%)
M.  aichiense ATCC27280M.  sphagni (95.8%) M.  vaccae (94.7%)
M.  sphagni ATCC33027M.  aichiense (95.8%) M.  porcinum (95.5%)
M.  tokaiense ATCC27282M.  porcinum (95.3%) M.  smegmatis (93.7%)
M.  chelonae ATCC35752M.  obuense (93.9%) M.  phlei (93.5%)
M.  abscessus  subspecies   massiliense str. GO 06M.  abscessus (98.7%) M.  sphagni (91.4%)
M.  abscessus ATCC19977M.  abscessus subspecies massiliense (98.7%) M.  chubuense (92.0%)
M.  tuberculosis  H37RvATCC27294M.  tuberculosis H37Ra/M.  africanum/M.  bovis BCG/M.  bovis (100%) M.  canettii (99.8%)
M.  tuberculosis  H37RaATCC25177M.  tuberculosis H37Rv/M.  africanum/M.  bovis BCG/M.  bovis (100%) M.  canettii (99.8%)
M.  africanum ATCC25420M.  tuberculosis H37Ra /M.  tuberculosis Rv/M.  bovisBCG/M.  bovis (100%) M.  canettii (99.8%)
M.  bovis BCGATCC35735M.  tuberculosis H37Ra /M.  tuberculosis H37Rv /M.  africanum/M.  bovis (100%) M.  canettii (99.8%)
M.  bovis ATCC19210M.  tuberculosis H37Ra/M.  tuberculosis H37Rv/M.  africanum/M.  bovis BCG (100%) M.  canettii (99.8%)
M.  canettii CIPT140010059M.  tuberculosis/M.  africanum/M.  bovis BCG/M.  bovis (99.8%) M.  asiaticum (93.2%)
M.  xenopi ATCC19250M.  asiaticum (93.0%) M.  avium subspecies paratuberculosis (92.4%)
M.  avium ATCC25291M.  avium subspecies paratuberculosis (99.6%) M.  asiaticum (95.8%)
M.  avium subspecies paratuberculosis K-10M.  avium (99.4%) M.  asiaticum (96%)
M.  malmoense ATCC29571M.  intracellulare  (94.3%) M.  szulgai (93.9%)
M.  intracellulare ATCC13950M.  avium (95.6%) M.  avium subspecies paratuberculosis (95.4%)
M.  scrofulaceum ATCC19981M.  asiaticum (95.4%) M.  szulgai (95.1%)
M.  marinum ATCC927M.  ulcerans (99.6%) M.  szulgai (96.0%)
M.   ulcerans Agy99M.  marinum (99.6%) M.  szulgai (96.0%)
M.  szulgai ATCC25799M.  marinum/M.  ulcerans (96.0%) M.  asiaticum (95.6%)
M.  gordonae ATCC14470M.  asiaticum (95.3%) M.  kansasii (95.1%)
M.  asiaticum ATCC25276M.  simiae (97.7%) M.  kansasii (96.4%)
M.  kansasii ATCC12478M.  asiaticum (96.3%) M.  gastri (95.8%)
M.  gastri ATCC15754M.  kansasii (95.8%) M.  asiaticum (94.5%)
Rhodococcus equi 103SM.  smegmatis (90.3%) M.  intracellulare (90.1%)

Beside the queried type strain itself, the first matched or the second matched type strain.

2.2. rpsA Gene Amplification and Sequencing

DNA was released from cultured mycobacteria by boiling the cultured mycobacterial suspension in TE buffer for 10 min. After centrifugation, the supernatant was used for PCR amplification [10]. The primers used were forward primer, 5′-CCCTACATCGGCAAGGAG-3′, position 487–504 in the rpsA gene of Mycobacterium tuberculosis, GenBank accession number NC_000962.2, and reverse primer, 5′-TGTCGATGACCTTGACCATC-3′, position 1032–1051 in the rpsA gene of Mycobacterium tuberculosis, GenBank accession number NC_000962.2. The amplified product was 565 bp. PCR products were purified and sequenced by a commercial company (YINGJUN Biotech Company, Beijing, China) using ABI 3730 DNA Analyzer (Applied Biosystems, California, USA).

2.3. Sequence Analysis and the Phylogenetic Tree Constructions

In addition to the 42 reference strains, the rpsA sequences of 5 other mycobacterial species were obtained directly from GenBank, including M. avium subspecies paratuberculosis (GenBank accession number NC_002944.2), M. ulcerans (GenBank accession number NC_008611.1), M. vanbaalenii (GenBank accession number NC_008726.1), M. abscessus subspecies massiliense (GenBank accession number NC_018150.1), and M. canettii (GenBank accession number NC_015848.1). All sequences were aligned and homology was calculated by using the multiple alignment algorithm in the MegAlign package (Windows version 7.1.0; DNASTAR, Madison, Wis). 527 bp of sequence (excluding 38 nucleotides at each end of the amplicon, corresponding to the primer binding sites) was analyzed for the phylogenetic tree construction by the neighbor-joining method using the MEGA 5.05 package (http://www.megasoftware.net/mega.php). A bootstrap analysis (1000 repeats) using Rhodococcus equi (GenBank accession number NC_006361.1) as the outgroup was performed to evaluate the topology of the phylogenetic tree.

2.4. Species Identification of the Clinical Isolates

The 188 clinical isolates were analyzed blindly. Sequences, minus the known PCR primer sequences, were assembled and edited by using SeqMan software (version 7.1.0; DNASTAR, Madison, WI). Isolates were identified by comparing sequences by using a FASTA BLASTn search with MegAlign (version 7.1.0; DNASTAR, Madison, WI) to an in-house database of sequences consisting of type and reference strains from external culture collections.

3. Results

3.1. rpsA Sequence Alignment of the Reference Strains

Between 85.4% and 100% sequence homology (interspecies divergence, 0% to 14.6%) was observed among the 42 tested reference strains and the 5 additional mycobacterial species whose sequences were obtained from the GenBank database (Table 2). All the rpsA gene sequences of the analyzed mycobacteria strains were distinct from the outgroup strain R. equi. Among the 19 reference strains of the slow-growing Mycobacterium genus, 11 strains were greater than 97% homology, including 5 M. tuberculosis complex strains. Among the 28 reference strains of the rapid-growing Mycobacterium genus, 14 strains were greater than 97% homology (Table 1). The pathogenic M. kansasii was easily differentiated from the nonpathogenic M. gastri (95.8% homology), while those two species were not distinguishable by the 16S rDNA sequence alignment. The sequence homologies between various species were 91.8% between M. chelonae and M. abscessus, 95.6% between M. avium and M. intracellulare, and 93.9% between M. szulgai and M. malmoense. However, the sequence homology between other species was higher; for example, it was 99.6% between M. ulcerans and M. marinum and 98.7% between M. abscessus subspecies massiliense and M. abscessus. All members of the M. tuberculosis complex had identical sequences as did M. senegalense and M. thermoresistibile, M. parafortuitum and M. trivial, M. diernhoferi and M. duvalii, and M. austroafricanum and M. terrae.


Percent identity
123456789101112131415161718192021222324252627282930313233343536373839404142434445464748

1100.098.798.196.095.395.496.496.294.594.595.194.594.592.093.994.194.394.594.592.895.693.794.793.292.488.889.291.191.191.191.191.191.391.593.994.191.393.092.892.692.692.493.493.492.492.889.21M.  thermoresistibile
20.098.798.196.095.395.496.496.294.594.595.194.594.592.093.994.194.394.594.592.895.693.794.793.292.488.889.291.191.191.191.191.191.391.593.994.191.393.092.892.692.692.493.493.492.492.889.22M.  senegalense
31.31.398.396.695.895.896.896.894.794.795.695.495.492.094.195.494.995.195.193.096.294.795.694.593.089.489.891.191.191.191.191.191.392.094.594.792.093.493.593.293.292.893.793.793.092.889.43M.  porcinum
41.91.91.796.295.895.896.296.294.194.195.494.594.592.094.394.594.394.594.593.996.294.994.993.793.790.190.590.990.990.990.990.991.192.094.794.991.393.293.293.093.492.893.593.792.893.089.24M.  fortuitum
54.14.13.53.994.395.496.496.093.993.993.494.394.391.594.394.394.594.794.791.594.394.194.992.891.189.890.390.190.190.190.190.190.391.193.593.791.393.493.093.593.593.093.793.493.092.089.65M.  phlei
64.94.94.34.35.995.195.394.193.793.794.794.194.191.393.793.593.493.593.593.293.793.793.793.092.488.689.090.790.790.790.790.790.991.594.394.190.793.292.091.892.291.593.093.291.891.588.06M.  gadium
74.74.74.34.34.75.196.696.095.495.494.794.794.793.095.195.395.695.895.893.294.194.994.792.492.090.190.592.492.492.492.492.492.692.894.994.793.094.194.993.793.793.593.994.392.693.490.37M.  flavescens
83.73.73.33.93.74.93.596.696.696.695.396.296.292.096.296.296.496.296.293.594.794.594.593.991.590.590.991.191.191.191.191.191.392.293.793.592.094.393.593.293.293.093.493.592.492.290.98M.  smegmatis
93.93.93.33.94.16.14.13.596.696.695.196.896.892.495.895.696.296.496.493.095.694.994.792.492.091.591.890.990.990.990.990.991.192.294.394.592.493.992.893.493.493.293.993.592.892.889.29M.  vaccae
105.75.75.56.16.36.54.73.53.5100.094.996.296.293.094.995.896.496.696.693.094.392.893.092.491.789.890.189.989.989.989.989.990.191.592.692.491.393.292.091.191.191.391.891.590.391.789.010M.  parafortuitum
115.75.75.56.16.36.54.73.53.50.094.996.296.293.094.995.896.496.696.693.094.392.893.092.491.789.890.189.989.989.989.989.990.191.592.692.491.393.292.091.191.191.391.891.590.391.789.011M.  triviale
125.15.14.54.77.05.55.54.95.15.35.396.496.491.894.595.395.195.395.394.195.192.892.292.693.090.390.789.989.989.989.989.990.191.393.292.889.992.092.091.791.791.191.892.490.989.988.212M.  neoaurum
135.75.74.75.75.96.15.53.93.33.93.93.7100.092.695.395.895.695.895.894.194.994.193.792.892.291.391.789.489.489.489.489.489.690.792.091.891.193.091.891.591.591.391.791.891.390.589.813M.  diernhoferi
145.75.74.75.75.96.15.53.93.33.93.93.70.092.695.395.895.695.895.894.194.994.193.792.892.291.391.789.489.489.489.489.489.690.792.091.891.193.091.891.591.591.391.791.891.390.589.814M.  duvalii
158.48.48.48.49.19.37.48.58.07.47.48.67.87.893.093.993.093.293.291.791.891.191.189.491.189.890.189.289.289.289.289.289.488.890.190.591.590.590.589.689.889.889.290.189.689.087.715M.  aurum
166.36.36.15.95.96.55.13.94.35.35.35.74.94.97.496.696.095.895.893.994.794.793.592.892.491.892.490.590.590.590.590.590.791.793.293.091.192.692.492.693.091.392.493.591.590.988.416M.  chubuense
176.16.14.75.75.96.84.93.94.54.34.34.94.34.36.33.597.297.097.094.195.694.594.193.593.491.391.790.390.390.390.390.390.591.393.093.292.092.893.292.693.091.792.892.891.591.190.117M.  gilvum
185.95.95.35.95.77.04.53.73.93.73.75.14.54.57.44.12.999.899.894.194.992.893.492.693.090.791.190.190.190.190.190.190.391.392.892.691.392.492.892.092.091.392.092.491.392.089.218M.  vanbaalenii
195.75.75.15.75.56.84.33.93.73.53.54.94.34.37.24.33.10.2100.093.995.193.093.592.893.290.991.390.390.390.390.390.390.591.593.092.891.592.693.092.292.291.592.292.691.592.289.019M.  terrae
205.75.75.15.75.56.84.33.93.73.53.54.94.34.37.24.33.10.20.093.995.193.093.592.893.290.991.390.390.390.390.390.390.591.593.092.891.592.693.092.292.291.592.292.691.592.289.020M.  austroafricanum
217.67.67.46.39.17.27.26.77.47.47.46.16.16.18.96.36.16.16.36.393.791.891.791.192.688.889.289.489.489.489.489.489.690.391.591.790.190.990.989.990.189.690.791.390.989.989.621M.  rhodesiae
224.54.53.93.95.96.66.15.54.55.95.95.15.35.38.75.54.55.35.15.16.593.794.593.494.190.390.790.390.390.390.390.390.590.593.793.990.392.692.892.292.291.892.693.091.790.989.222M.  obuense
236.56.55.55.36.16.55.35.75.37.67.67.66.16.19.55.55.77.67.47.48.66.595.893.993.091.792.090.990.990.990.990.991.192.293.793.991.793.292.693.493.792.293.793.992.491.789.023M.  aichiense
245.55.54.55.35.36.55.55.75.57.47.48.26.56.59.56.76.17.06.76.78.95.74.393.793.491.391.891.591.591.591.591.591.791.593.794.192.293.093.592.893.093.093.494.193.292.489.224M.  sphagni
257.27.25.76.57.67.48.06.38.08.08.07.87.67.611.57.66.77.87.67.69.57.06.36.592.088.889.489.689.689.689.689.689.890.594.594.590.793.092.891.892.292.092.493.492.291.588.225M.  tokaiense
268.08.07.46.59.58.08.49.18.48.98.97.48.28.29.58.07.07.47.27.27.86.17.47.08.491.591.889.689.689.689.689.689.889.892.292.689.490.591.190.791.189.892.091.790.790.387.126M.  chelonae
2712.112.111.510.611.012.410.610.19.111.011.010.49.39.311.08.69.39.99.79.712.110.48.99.312.19.198.786.786.786.786.786.786.988.488.488.888.488.689.489.689.688.490.189.888.287.187.127M.  abscessus
2811.711.711.010.110.411.910.19.78.610.610.69.98.98.910.68.08.99.59.39.311.79.98.48.611.58.71.387.587.587.587.587.587.788.888.889.288.889.089.889.989.989.090.590.188.687.587.728M.  abscessus
299.59.59.59.710.610.08.09.59.710.810.810.811.511.511.710.110.410.610.410.411.510.49.79.111.311.314.613.7100.0100.0100.0100.099.890.991.792.090.991.792.292.091.892.092.893.492.291.888.029M.  tuberculosis
309.59.59.59.710.610.08.09.59.710.810.810.811.511.511.710.110.410.610.410.411.510.49.79.111.311.314.613.70.0100.0100.0100.099.890.991.792.090.991.792.292.091.892.092.893.492.291.888.030M.  tuberculosis
319.59.59.59.710.610.08.09.59.710.810.810.811.511.511.710.110.410.610.410.411.510.49.79.111.311.314.613.70.00.0100.0100.099.890.991.792.090.991.792.292.091.892.092.893.492.291.888.031M.  bovis
329.59.59.59.710.610.08.09.59.710.810.810.811.511.511.710.110.410.610.410.411.510.49.79.111.311.314.613.70.00.00.0100.099.890.991.792.090.991.792.292.091.892.092.893.492.291.888.032M.  africanum
339.59.59.59.710.610.08.09.59.710.810.810.811.511.511.710.110.410.610.410.411.510.49.79.111.311.314.613.70.00.00.00.099.890.991.792.090.991.792.292.091.892.092.893.492.291.888.033M.  bovis
349.39.39.39.510.49.77.89.39.510.610.610.611.211.211.59.910.210.410.210.211.210.29.58.911.111.014.413.50.20.20.20.20.291.191.892.290.791.892.491.891.791.892.693.292.091.788.034M.  canettii
359.19.18.48.49.59.17.68.28.29.19.19.39.99.912.28.99.39.39.19.110.410.28.29.110.211.012.612.19.89.89.89.89.89.592.292.490.791.792.091.191.190.591.893.091.590.987.735M.  xenopi
366.36.35.75.56.75.95.36.55.97.87.87.28.48.410.67.27.47.67.47.49.16.56.56.55.78.212.612.18.98.98.98.98.98.78.299.4
93.095.694.793.994.394.195.195.894.793.788.636M.  avium
376.16.15.55.36.56.15.56.75.78.08.07.68.68.610.17.47.27.87.67.68.96.36.36.15.77.812.111.78.48.48.48.48.48.28.00.693.595.495.194.594.994.595.496.095.393.988.837M.  avium
389.39.38.49.39.39.97.48.48.09.39.310.89.59.59.19.58.49.39.19.110.610.48.98.29.911.512.612.19.79.79.79.79.710.09.97.46.794.393.793.593.492.493.993.093.793.089.038M.  malmoense
397.47.47.07.27.07.26.15.96.37.27.28.47.47.410.27.87.68.07.87.89.77.87.27.47.410.212.311.98.98.98.98.98.98.78.94.54.75.994.594.194.193.495.194.794.394.390.739M.  intracellulare
407.67.66.77.27.48.45.36.77.68.48.48.48.68.610.18.07.27.67.47.49.77.67.86.77.69.511.511.08.28.28.28.28.28.08.45.55.16.55.794.394.394.995.195.493.993.489.640M.  scrofulaceum
417.87.87.27.46.88.66.57.27.09.59.58.99.19.111.27.87.88.48.28.210.88.27.07.68.79.911.210.88.48.48.48.48.48.79.56.35.76.86.15.999.694.596.095.695.693.789.641M.  marinum
427.87.87.27.06.88.26.57.27.09.59.58.99.19.111.07.47.48.48.28.210.68.26.57.48.29.511.210.88.78.78.78.78.78.99.55.95.37.06.15.90.494.396.095.695.493.989.442M.  ulcerans
438.08.07.67.67.49.16.77.47.29.39.39.59.39.311.09.38.99.39.19.111.28.68.27.48.411.012.611.98.48.48.48.48.48.710.16.15.78.07.05.35.75.994.395.395.193.489.043M.  gordonae
447.07.06.56.76.57.46.37.06.38.68.68.68.98.911.78.07.68.48.28.29.97.86.57.08.08.410.610.17.67.67.67.67.67.88.65.14.76.35.15.14.14.15.995.695.494.590.144M.  szulgai
457.07.06.56.57.07.25.96.76.79.19.18.08.68.610.66.77.68.07.87.89.37.46.36.17.08.911.010.67.07.07.07.07.07.27.44.34.17.45.54.74.54.54.94.596.494.788.845M.  asiaticum
468.08.07.47.67.48.67.88.07.610.410.49.79.39.311.29.19.19.39.19.19.78.98.07.28.29.912.812.38.28.28.28.28.28.59.15.54.96.55.96.34.54.75.14.73.795.888.846M.  kansasii
477.67.67.67.48.49.17.08.27.68.98.910.810.110.111.99.79.58.48.28.210.89.78.98.09.110.414.213.78.78.78.78.78.78.99.76.56.37.45.97.06.66.37.05.75.54.387.947M.  gastri
48   11.711.711.511.711.213.010.49.711.711.911.912.811.011.013.512.610.611.711.911.911.211.711.911.712.814.214.213.513.013.013.013.013.013.013.512.412.111.99.911.211.211.511.910.612.112.113.248Rhodococcus equi

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748

3.2. Phylogenetic Tree Construction

A phylogenetic tree, which provided the basis for species differentiation in the genus Mycobacterium, was constructed (Figure 1). The absolute majority tested species showed good separation. The rapid-growing species were well defined from the slow-growing species in the tree. M. chelonae, M. abscessus subspecies massiliense, and M. abscessus, which are categorized in the pathogenic taxonomic group of rapid-growing mycobacteria, formed a distinctive cluster which was much closer to the slow-growing species compared with the nonpathogenic group of rapid-growing mycobacteria. The reliability of the phylogenetic tree was verified by the bootstrap method, using R. equi as the outgroup.

3.3. Species Identification Outcomes of the Clinical Isolates

rpsA sequence and alignment efficiently identified clinical isolates representing 8 mycobacterial species. We found that a criterion of first distinct rpsA sequence match >97% confirmed the identification of 183 of 188 (97.3%) clinical isolates. 20 M. abscessus strains, 2 M. fortuitum strains, and 10 M. avium strains were identified to the subspecies level. Identification discrepancies between rpsA and other sequences were only encountered with 4 M. gordonae strains and 1 M. kansasii strain (see Table 3). No discrepancies were found in the species submitted as M. intracellulare, M. avium, M. abscessus, and M. tuberculosis complex.


Primary identificationb ()Primary identificationb ()
First choiceSecond choice
Identification ()%matchIdentification ()%match

M.  intracellulare (75)M.  intracellulare (4)100%M.  avium (4)95.6%
M.  intracellulare (1)99.8%M.  avium (1)95.3%
M.  intracellulare (53)98.9%M.  avium (53)95.1%
M.  intracellulare (10)98.7%M.  avium (10)94.9%
M.  intracellulare (1)98.5%M.  avium (1)95.3%
M.  intracellulare (3)98.1%M.  avium (3)94.7%
M.  intracellulare (1)97.5%M.  avium (1)93.2%
M.  intracellulare (2)97.1%M.  avium (2)96.1%

M.  abscessus (39)M.  abscessus (9)100%M.  abscessus  subspecies  massiliense (9)98.7%
M.  abscessus (10)99.8%M.  abscessus  subspecies  massiliense (10)98.5%
M.  abscessus  subspecies  massiliense (20)99.4%M.  abscessus (20)99.1%

M.  kansasii (23)M.  kansasii (22)100%M.  asiaticum (22)96.3%
M.  szulgai (1)96.6%M.  marinum/M.  ulcerans (1)96.0%

M.  avium (10)M.  avium  subspecies paratuberculosis (5)99.6%M.  avium (5)99.4%
M.  avium subspecies paratuberculosis (4)99.8%M.  avium (4)99.2%
M.  avium subspecies paratuberculosis (1)97.5%M.  avium (1)97.0%

M.  fortuitum  complex (17)M.  fortuitum (2)100%M.  senegalense/M.  thermoresistibile (2)98.1%
M.  fortuitum (11)99.8%M.  porcinum (11)98.5%
M.  fortuitum (1)99.6%M.  porcinum (1)98.3%
M.  fortuitum/M.  porcinum (1)98.9%M.  senegalense/M.  thermoresistibile (1)98.3%
M.  thermoresistibile/M.  senegalense (2)99.6%M.  fortuitum (2)98.1%

M.  gordonae (7)M.  gordonae (1)99.2%M.  asiaticum (1)95.5%
M.  gordonae (2)98.5%M.  asiaticum (2)96.3%
M.  gordonae (3)95.9%M.  asiaticum (3)95.5%
M.  szulgai (1)95.9%M.  gordonae (1)94.7%

M.  tuberculosis  complex (16)M.  tuberculosis (2)100%M.  canetti (2)99.8%
M.  tuberculosis (12)99.8%M.  canetti (12)99.6%
M.  tuberculosis (2)99.6%M.  canetti (2)99.4%

M.  neoaurum (1)M.  neoaurum (1)98.7%M.  smegmatis (1)95.4%

Numbers in parentheses () represent the numbers of isolates identified as a particular species.
bIdentification based on sequencing of at least two of the following: 16S rDNA, 16-23S rRNA gene internal transcribed spacer (ITS), and rpoB and hsp65 genes.

The sequence divergence among MTC members was 0.4%. The intraspecies divergence among the 188 clinical strains ranged from 0.6% to 5.3% (Table 3). The sequence diversity among the M. abscessus was 0.6%, while that among the M. gordonae clinical isolates was 5.3%. No M. gordonae clinical strain had an identical sequence with that of the M. gordonae reference strain (ATCC14470). Interestingly, the rpsA sequence of one strain among the 23 M. kansasii isolates exhibited a relatively low level of sequence similarity (95.6%) to that of the reference strain (ATCC 12478), while those of all the other M. kansasii isolates were identical. Since the strain was confirmed as M. kansasii by rpoB and hsp65 gene alignment, it might suggest a distant variant or a new subtype.

4. Discussion

16S rDNA gene sequence alignment has been used as the reference method for mycobacterial species identification. However, it has been reported that, by using the 16S rDNA gene alone for species identification of clinical NTM, 37% of such isolates remained unclassified which illustrates the need for additional molecular tools for proper phylogenetic assignment and accurate NTM identification [11]. The common assumption that bacterial isolates belong to the same species if they have fewer than 5–15 bp differences within the 16S rDNA gene sequence [12] or if they have more than 97% 16S rDNA gene sequence identities [13] may not be applicable to genus of Mycobacterium, whose members are much more closely related to each other. Furthermore, some mycobacterial strains shared uniform sequence of 16S rDNA gene, such as M. intracellulare and M. avium, M. kansasii and M. gastri, M. abscessus and M. chelonae, and M. szulgai and M. malmoense.

In this study, we found that the rpsA gene is useful as a complementary method to the 16S rDNA for mycobacterial species identification. Among the 47 reference strains and totally 1081 paired comparisons {1081 = []; , , , the rpsA gene alone can differentiate 97.2% (1051 out of 1081) of them if using 97% sequence homology as cutoff value. When referring to the phylogenetic tree, the rpsA resolution might increase. For the most commonly seen pathogenic NTM species such as M. intracellulare, M. avium, M. abscessus, M. chelonae, and M. xenopi, rpsA gene sequence could differentiate these species easily. Furthermore, 183 of the 188 (97.3%) clinical isolates, representing 8 mycobacterial species, were identified correctly by rpsA gene blast. Discrepancies of identification were only encountered with 4 M. gordonae strains and 1 M. kansasii strain. Both these species have reportedly more intraspecies sequence divergence [9, 1416]. Three out of the 4 M. gordonae strains with discrepant identification still had the first distinction as M. gordonae by rpsA, but the sequence identity was lower than 97%, which means the identification might be correct when a more complete in-house database is being developed.

rpsA alone cannot differentiate the following: between members of the M. tuberculosis complex, between M. senegalense and M. thermoresistibile, between M. parafortuitum and M. trivial, between M. diernhoferi and M. duvalii, and between M. austroafricanum and M. terrae; however, these species are also difficult to be separated by other markers alone such as 16s rDNA, ITS, rpoB, and hsp65. Several reference strains had similar homology with two or more species which suggests inadequate taxonomy within the currently described species. Based on nucleotide sequences of rpoB, hsp65, and sodA in clinical isolates of the Mycobacterium abscessus group, one-fourth of isolates had discordant identification [16]. Multilocus sequence typing and sequence analysis of several genes facilitate the identification of closely related species or subspecies [17].

In our study, we have demonstrated that the rpsA homology gene is a promising marker for mycobacterial species identification for both fast- and slow-growing mycobacteria. To be a good marker for species differentiation, the target gene should be stable and sequence variations should occur randomly; additionally, an extremely conserved or highly variable gene may not be adequate. As a single-copy house-keeping gene, rpsA gene could work well as a target gene for Mycobacterium species discrimination without ambiguous identification. The 16S rDNA gene has higher homology within the mycobacteria compared to the rpsA gene: 94.3% to 100% compared to 85.4% to 100%. According to our work, when used as sole marker, rpsA had better resolution than 16s rDNA but had similar resolution as ITS, rpoB, and hsp65.

Due to the high resolution power of rpsA for species identification, we presume it has potential of clinical use. From our own experience, we recommend that when performing species identification by homologue DNA sequence comparison, one should start with 16s RNA gene, plus at least one of the other markers such as ITS, rpoB, hsp65, and rpsA. When conflicting outcome or dubious outcomes yield, more markers should be added further. Even though 16s RNA gene is inferior to the other markers considering the capacity, it should be chosen firstly since 16s RNA gene has the most robust sequence database, which can help to avoid big identification errors due to unreliable database of other markers [18].

The main deterrents for the primary use of rpsA sequencing as a routine means of identifying mycobacteria reside in the need for a comprehensive database. Therefore, we constructed our own rpsA sequence database by including as many type strains as possible and integrating rpsA sequences deposited in GenBank such as M. avium subspecies paratuberculosis, M. ulcerans, M. vanbaalenii, M. abscessus subspecies massiliense, and M. canettii. Besides type strains, rpsA sequences of confirmed clinical strains were also included. We will constantly upgrade our database and the capability to identify the most recent described species.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Authors’ Contribution

Hongfei Duan and Guan Liu contributed equally to this study.

Acknowledgments

The work was supported by the research funding from Infectious Diseases Special Project, Ministry of Health of China (2012ZX10003002-09), and Beijing Municipal Administration of Hospitals Clinical Medicine Development of Special Funding (ZYLX201304). All the type strains and clinical isolates used in this study were obtained from the “Beijing Bio-Bank of Clinical Resources on Tuberculosis (D09050704640000)," Beijing Chest Hospital.

References

  1. T. Adékambi, P. Berger, D. Raoult, and M. Drancourt, “rpoB gene sequence-based characterization of emerging non-tuberculous mycobacteria with descriptions of Mycobacterium bolletii sp. nov., Mycobacterium phocaicum sp. nov. and Mycobacterium aubagnense sp. nov,” International Journal of Systematic and Evolutionary Microbiology, vol. 56, no. 1, Article ID 63969, pp. 133–143, 2006. View at: Publisher Site | Google Scholar
  2. B. Ninet, M. Monod, S. Emler et al., “Two different 16S rRNA genes in a mycobacterial strain,” Journal of Clinical Microbiology, vol. 34, no. 10, pp. 2531–2536, 1996. View at: Google Scholar
  3. U. Reischl, K. Feldmann, L. Naumann et al., “16S rRNA sequence diversity in Mycobacterium celatum strains caused by presence of two different copies of 16S rRNA gene,” Journal of Clinical Microbiology, vol. 36, no. 6, pp. 1761–1764, 1998. View at: Google Scholar
  4. D. E. Griffith, T. Aksamit, B. A. Brown-Elliott et al., “An official ATS/IDSA statement: diagnosis, treatment, and prevention of nontuberculous mycobacterial diseases,” American Journal of Respiratory and Critical Care Medicine, vol. 175, no. 4, pp. 367–416, 2007. View at: Publisher Site | Google Scholar
  5. W. Shi, X. Zhang, X. Jiang et al., “Pyrazinamide inhibits trans-translation in Mycobacterium tuberculosis,” Science, vol. 333, no. 6049, pp. 1630–1632, 2011. View at: Publisher Site | Google Scholar
  6. S. T. Cole, R. Brosch, J. Parkhill et al., “Deciphering the biology of Mycobacterium tuberculosis from the complete genome sequence,” Nature, vol. 393, no. 6685, pp. 537–544, 1998. View at: Publisher Site | Google Scholar
  7. A. Roth, M. Fischer, M. E. Hamid, S. Michalke, W. Ludwig, and H. Mauch, “Differentiation of phylogenetically related slowly growing mycobacteria based on 16S-23S rRNA gene internal transcribed spacer sequences,” Journal of Clinical Microbiology, vol. 36, no. 1, pp. 139–147, 1998. View at: Google Scholar
  8. H. Ringuet, A. Honore, A. Varnerot et al., “hsp65 sequencing for identification of rapidly growing mycobacteria,” Journal of Clinical Microbiology, vol. 37, no. 3, pp. 852–857, 1999. View at: Google Scholar
  9. B.-J. Kim, S.-H. Lee, M.-A. Lyu et al., “Identification of mycobacterial species by comparative sequence analysis of the RNA polymerase gene (rpoB),” Journal of Clinical Microbiology, vol. 37, no. 6, pp. 1714–1720, 1999. View at: Google Scholar
  10. J. Amita, T. Vandana, R. S. Guleria, and R. K. Verma, “Qualitative evaluation of mycobacterial DNA extraction protocols for polymerase chain reaction,” Molecular Biology Today, vol. 3, no. 2, pp. 43–49, 2002. View at: Google Scholar
  11. R. J. Pauls, C. Y. Turenne, J. N. Wolfe, and A. Kabani, “A high proportion of novel mycobacteria species identified by 16S rDNA analysis among slowly growing AccuProbe-negative strains in a clinical setting,” American Journal of Clinical Pathology, vol. 120, no. 4, pp. 560–566, 2003. View at: Publisher Site | Google Scholar
  12. G. E. Fox, J. D. Wisotzkey, and P. Jurtshuk Jr., “How close is close: 16S rRNA sequence identity may not be sufficient to guarantee species identity,” International Journal of Systematic Bacteriology, vol. 42, no. 1, pp. 166–170, 1992. View at: Publisher Site | Google Scholar
  13. M. Drancourt, P. Berger, and D. Raoult, “Systematic 16S rRNA gene sequencing of atypical clinical isolates identified 27 new bacterial species associated with humans,” Journal of Clinical Microbiology, vol. 42, no. 5, pp. 2197–2202, 2004. View at: Publisher Site | Google Scholar
  14. P. Kirschner, E. C. Bottger, D. T. Walton, and M. Valesco, “Microheterogeneity within rRNA of Mycobacterium gordonae,” Journal of Clinical Microbiology, vol. 30, no. 4, pp. 1049–1050, 1992. View at: Google Scholar
  15. C. Taillard, G. Greub, R. Weber et al., “Clinical implications of Mycobacterium kansasii species heterogeneity: swiss national survey,” Journal of Clinical Microbiology, vol. 41, no. 3, pp. 1240–1244, 2003. View at: Publisher Site | Google Scholar
  16. E. Macheras, A.-L. Roux, F. Ripoll et al., “Inaccuracy of single-target sequencing for discriminating species of the Mycobacterium abscessus group,” Journal of Clinical Microbiology, vol. 47, no. 8, pp. 2596–2600, 2009. View at: Publisher Site | Google Scholar
  17. L. A. Almeida and R. Araujo, “Highlights on molecular identification of closely related species,” Infection, Genetics and Evolution, vol. 13, no. 1, pp. 67–75, 2013. View at: Publisher Site | Google Scholar
  18. J. Dai, Y. Chen, and M. Lauzardo, “Web-accessible database of hsp65 sequences from Mycobacterium reference strains,” Journal of Clinical Microbiology, vol. 49, no. 6, pp. 2296–2303, 2011. View at: Publisher Site | Google Scholar

Copyright © 2015 Hongfei Duan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


More related articles

1142 Views | 424 Downloads | 7 Citations
 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly as possible. We will be providing unlimited waivers of publication charges for accepted research articles as well as case reports and case series related to COVID-19. Review articles are excluded from this waiver policy. Sign up here as a reviewer to help fast-track new submissions.