#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

New Mycobacteroides abscessus subsp. massiliense strains with recombinant hsp65 gene laterally transferred from Mycobacteroides abscessus subsp. abscessus: Potential for misidentification of M. abscessus strains with the hsp65-based method


Authors: Byoung-Jun Kim aff001;  Ga-Na Kim aff001;  Bo-Ram Kim aff001;  Tae-Sun Shim aff002;  Yoon-Hoh Kook aff001;  Bum-Joon Kim aff001
Authors place of work: Department of Microbiology and Immunology, Biomedical Sciences, Liver Research Institute and Cancer Research Institute, College of Medicine, Seoul National University, Seoul, Korea aff001;  Division of Pulmonary and Critical Care Medicine, Department of Internal Medicine, University of Ulsan College of Medicine, Asan Medical Center, Seoul, Korea aff002
Published in the journal: PLoS ONE 14(9)
Category: Research Article
doi: https://doi.org/10.1371/journal.pone.0220312

Summary

It has been reported that lateral gene transfer (LGT) events among Mycobacteroides abscessus strains are prevalent. The hsp65 gene, a chronometer gene for bacterial phylogenetic analysis, is resistant to LGT events, particularly among mycobacterial strains, rendering the hsp65-targeting method the most widely used method for mycobacterial detection. To determine the prevalence of M. abscessus strains that are subject to hsp65 LGT, we applied rpoB typing to 100 clinically isolated Korean strains of M. abscessus that had been identified by hsp65 sequence analysis. The analysis indicated the presence of 2 rough strains, showing a discrepancy between the 2 typing methods. MLST analysis based on the partial sequencing of seven housekeeping genes, erm(41) PCR and further hsp65 PCR-restriction enzyme and polymorphism analysis (PRA) were conducted to identify the two strains. The MLST results showed that the two strains belong to M. abscessus subsp. massiliense and not to M. abscessus subsp. abscessus, as indicated by the rpoB-based analysis, suggesting that their hsp65 genes are subject to LGT from M. abscessus subsp. abscessus. Further analysis of these strains using the hsp65 PRA method indicated that these strains possess a PRA pattern identical to that of M. abscessus subsp. abscessus and distinct from that of M. abscessus subsp. massiliense. In conclusion, we identified two M. abscessus subsp. massiliense rough strains from Korean patients with hsp65 genes that might be laterally transferred from M. abscessus subsp. abscessus. To the best of our knowledge, this is the first demonstration of possible LGT events associated with the hsp65 gene in mycobacteria. Our results also suggest that there is the potential for misidentification when the hsp65-based protocol is used for mycobacterial identification.

Keywords:

Biology and life sciences – Organisms – Research and analysis methods – Molecular biology – Evolutionary biology – Database and informatics methods – Bioinformatics – Sequence analysis – Sequence alignment – Molecular biology techniques – Computer and information sciences – Evolutionary systematics – Phylogenetics – Phylogenetic analysis – Taxonomy – Data management – Medicine and health sciences – Microbiology – Medical microbiology – Microbial pathogens – Bacterial pathogens – Bacteria – Pathology and laboratory medicine – Pathogens – Actinobacteria – Mycobacterium tuberculosis – Artificial gene amplification and extension – Polymerase chain reaction – DNA sequence analysis – Artificial genetic recombination – Gene targeting – Mycobacteria

Introduction

Rapidly growing mycobacteria (RGM) are ubiquitous organisms that have gained increasing attention as important human pathogens [1, 2]. Among RGMs, infections due to the Mycobacteroides abscessus strains have shown increased worldwide clinical importance, and their incidence in cystic fibrosis patients has increased [3]. In South Korea, the incidence of lung diseases caused by M. abscessus has also been increasing, and this organism accounts for 70–80% of the lung disease caused by RGM [47]. M. abscessus can cause lung disease in immunocompetent individuals and shares a number of characteristics with M. tuberculosis, including the ability to induce granulomatous lesions or caseous necrosis [8]. Infections involving M. abscessus are notorious for being difficult to treat due both to the natural broad-spectrum antibiotic resistance of this species and to its acquired resistance, with disparate antibiotic susceptibility patterns being observed among clinical strains [9].

The taxonomy of M. abscessus strains remains problematic. Currently, these strains are divided into two subspecies, M. abscessus subsp. abscessus (the former species Mycobacteroides abscessus) and M. abscessus subsp. bolletii. M. abscessus subsp. bolletii was proposed to combine the two former species, M. massiliense and M. bolletii [10, 11]. M. massiliense can be further subdivided into two genotypes (Type I and Type II) based on hsp65 sequence analysis [1214]. Also, recent phylogenomics and comparative genome analyses on 150 genomes of Mycobacterium species revealed that the genus of Mycobacterium was divided four novel genera. In the case of M. abscessus-chelonae complex, their genus was emended into Mycobacteroides [15].

Lateral gene transfer (LGT) has been proposed as the major driving force for the acquisition of prokaryotic genetic diversity, an attribution that leads to better survival of prokaryotic organisms under harsh environmental conditions [16, 17]. The recent increase in available information on mycobacterial genomes supports the idea that LGT plays an important role in the evolutionary transition of mycobacteria from saprophytic organisms into opportunistic or specialized, highly persisting pathogens [18, 19]. In particular, interspecies or intraspecies LGT events among members of the M. abscessus strains have been reported to be very prevalent [20]. And it has reported that M. abscessus evolution is sporadically punctuated by dramatic genome wide remodelling events [21]. We recently identified six M. abscessus subsp. massiliense strains isolated from Korean patients in which the rpoB gene was laterally transferred from M. abscessus subsp. abscessus [22], suggesting the potential for misidentification when the rpoB tying method is used to the differentiate among M. abscessus strains. In addition, other reports have presented evidence for LGT events involving the rpoB gene, including the observation that the rpoBC operon of the Type I genotype of M. yongoenense has been laterally transferred from a distantly related strain, Mycobacterium parascrofulaceum [23].

The hsp65 (groEL2) gene, another chronometer molecule, has been widely used as a targeting molecule for mycobacterial identification and detection [24, 25] together with the rpoB gene [26, 27]. In this study, we sought to address issues associated with the possibility that LGT events involving the hsp65 gene occur among members of the M. abscessus strains. To this end, we applied rpoB typing (711 bp) to 100 clinically isolated Korean strains of M. abscessus that had already been identified by hsp65 sequence analysis (603 bp). In the strains in which the two methods yielded discordant results, additional phylogenetic analysis was conducted to confirm the authenticity of potential LGT events in their hsp65 genes.

Materials and methods

Mycobacterial strains and culture conditions

Of the total 206 clinical isolated strains, 106 strains were used in the previous paper. These strains were identified as M. massiliense using the hsp65-based method. Among these 106 strains, 6 strains showed different rpoB sequence from M. massiliense, however, these strains showed similar sequence homologies with M. massiliense. So, the possibility of lateral gene transfer in the rpoB gene of these strains were described in the previous study [22].

We used other 100 strains of the total 206 strains in this study which were identified as M. abscessus using the hsp65-based method [12] (Table 1). These M. abscessus strains included both rough and smooth morphotypes (28 and 72 strains, respectively). Among them, two strains, 55262 and 55184 showed different rpoB gene sequence from M. abscessus. Although both studies used samples collected at the same time period and applied similar experimental methods, the two studies were considered to be different, because different groups of samples (which were divided by hsp65-based method) were used and putative lateral gene transfer events were considered to be applied to different genes, respectively.

Tab. 1.

Separation of 100 M. abscessus subsp. abscessus clinical strains into genotype level by sequence analyses based on the partial hsp65 (603 bp) and rpoB (711 bp) gene sequences.

<h2>Separation of 100 <i>M</i>. <i>abscessus subsp</i>. <i>abscessus</i> clinical strains into genotype level by sequence analyses based on the partial <i>hsp65</i> (603 bp) and <i>rpoB</i> (711 bp) gene sequences.</h2>

All clinical strains were collected from the Asan Medical Center (Seoul, Republic of Korea) from January 2004 to June 2011. This work was approved by the institutional review board of the Asan Medical Center (2012–0170) with documentation for waivers of informed consent. Each bacterial isolate was maintained on Middlebrook 7H10 agar plates supplemented with OADC or in Middlebrook 7H9 broth medium supplemented with ADC at 37°C. The bacteria were stored as frozen stocks at -70°C by flash-freezing of intermediate passage samples in 20% glycerol. M. abscessus subsp. abscessus ATCC 19977T (= CIP 104536T), M. abscessus subsp. bolletii CIP 108541T and M. abscessus subsp. massiliense CIP 108297 (= CCUG 48898) were also used as type strains.

DNA extraction, PCR and sequencing

Bacterial DNA was extracted from individual clinical isolates using the bead beater-phenol extraction method [24] and used as templates for PCR amplification. Partial hsp65 (603 bp) and rpoB (711 bp) gene-targeted PCR was conducted in a total of 100 M. abscessus strains as described previously [12, 24]. MLST analyses targeting seven housekeeping genes were applied to investigate the genetic diversity of the two putative recombinant strains. The seven target genes were argH (argininosuccinate lyase), cya (adenylate cyclase), glpK (glycerol kinase), gnd (6-phosphogluconate dehydrogenase), murC (UDP N-acetylmuramate-L-alanine ligase), pta (phosphate acetyltransferase) and purH (phosphoribosylaminoimidazole carboxylase ATPase subunit) [28, 29]. Also, erm(41)-targeted PCR was applied to the reference strain and to two Rec-mas-H strains [30, 31]. Additionally, to confirm the putative recombination site of the Rec-mas-H strain within the hsp65 gene sequence, the entire hsp65 gene sequence of Asan 55262 was sequenced and compared with the hsp65 gene sequences of the M. abscessus subsp. abscessus CIP 104536T, M. abscessus subsp. massiliense CIP 108297T and M. abscessus subsp. massiliense Asan 50594 strains. The complete hsp65 gene sequence of each selected isolate was amplified using 5 primer sets. Detailed information on the primers is provided in S1 Table. All the PCR reactions were conducted as described previously [22].

PCR-based restriction analysis

The partial hsp65 gene sequence was amplified using the primer set Tb11 (5’- ACCAACGATGGTGTGTCCAT-3’) and Tb12 (5’- CTTGTCGAACCGCATACCCT-3’) as described by Telenti et al [25]. Restriction analyses were performed as described by Telenti et al [25]. Briefly, 10 μl of final amplicons were digested with BstEII (NEB, Ipswich, MA, UK; 37°C) and HaeIII (Takara Bio, Shiga, Japan; 60°C) for 2 hours. The two restriction fragments were separated by electrophoresis on 1% agarose gels with a 100-bp DNA ladder as the molecular size standard.

Sequence analyses

The obtained sequences of hsp65 (603 bp), rpoB (711 bp) and of 7 MLST target genes, such as argH (503 bp), cya (541 bp), glpK (563 bp), gnd (494 bp), murC (545 bp), pta (486 bp) and purH (549 bp) from the two putative recombinant strains were aligned with those of M. abscessus complex type strains using the ClustalW algorithm in the MEGA 4.0 program [32]. Phylogenetic trees based on each target gene or concatenated sequence were constructed by the neighbor-joining [33] and maximum parsimony [34] methods with 1,000 replicates [35].

Nucleotide sequence accession numbers

The hsp65, rpoB and 7 MLST gene sequences determined in this study were deposited in GenBank under the accession numbers MH430895—MH430913 and are listed in S2 Table.

Results

Identification of two M. abscessus strains for which discordant results were obtained in hsp65 and rpoB sequence analyses

Among 100 M. abscessus strains that were previously identified as M. abscessus subsp. abscessus by 603-bp hsp65 sequencing analysis for subspecies differentiation of the M. abscessus strains, two strains (2.0%, 2/100) were identified as M. abscessus subsp. massiliense by 711-bp rpoB sequence-based phylogenetic analysis (Table 1). The two Rec-mas-H strains (Asan 55262 and 55184) that yielded results discordant with those obtained by rpoB-based analysis had identical hsp65 sequences, with 1 bp mismatch (99.8% sequence homology) to the hsp65 sequence of M. abscessus subsp. abscessus type strain, CIP 104536T (Fig 1, Table 1). The results showed that no putative recombination site was present in the compared hsp65 gene sequences. Only three base pairs (located at 779, 918, and 1296 nt of the 1,626-bp hsp65 sequence) (S1 Fig) differed between the Asan 55262 and the M. abscessus subsp. abscessus CIP 104536T strains.

<h2>Phylogenetic trees based on the partial <i>hsp65</i> and <i>rpoB</i> gene sequences of two Rec-mas-H strains.</h2>
Fig. 1.

Phylogenetic trees based on the partial hsp65 and rpoB gene sequences of two Rec-mas-H strains.


Phylogenetic trees of 2 Rec-mas-H strains based on (A) the partial sequence of the hsp65 gene (603 bp) and (B) the partial sequence of the rpoB gene (711 bp). The trees were constructed using the neighbor-joining method in the MEGA 4.0 program. The bootstrap values were calculated from 1,000 replications; values <50% are not shown. The black-centered circles indicate that the corresponding clusters were supported with maximum parsimony-based trees. The bar indicates the number of base substitutions per site. The black-centered triangles indicate that the corresponding sequences were sequenced and obtained in this study.

Phylogenetic analysis of two Rec-mas-H strains based on seven different MLST genes

For precise species delineation of the two strains, MLST analyses based on the partial sequencing of seven housekeeping genes, argH, cya, glpK, gnd, murC, pta and purH, that have been previously used for the elucidation of recombination events in M. abscessus strains [29] were also performed in this study. Each single-gene-based tree was built from the sequences of each of the seven genes in the MLST scheme for the separation of the two strains at the subspecies level of the M. abscessus strains (Fig 2). With the exception of the glpK gene, which has a 1-bp difference in the sequence, the sequences of all the MLST genes were identical in the two Rec-mas-H strains. Most of the constructed phylogenetic trees (those for the argH, cya, glpK, gnd, murC, and pta genes) showed a topology similar to that of the rpoB gene sequence-based tree. The two Rec-mas-H strains were clustered together with the M. abscessus subsp. massiliense CIP 108297T strain or the M. abscessus subsp. massiliense Asan 50594 strain, which has been described as the Type II genotype of the M. abscessus subsp. massiliense strain [12, 14]. The argH and cya gene sequences of the two Rec-mas-H strains showed 100% sequence similarity with those of M. abscessus subsp. massiliense CIP 108297T (Table 2, Fig 2A and 2B). However, in the purH gene-based tree, the two Rec-mas-H strains were closely clustered with the M. abscessus subsp. massiliense Asan 50594 strain but not with the M. abscessus subsp. massiliense CIP 108297T strain. In addition, the purH gene sequence similarity between the two Rec-mas-H strains and M. abscessus subsp. massiliense CIP 108297T showed the lowest value of 97.6% among the seven MLST genes.

<h2>Neighbor-joining phylogenetic trees based on the 7 MLST genes of two Rec-mas-H strains.</h2>
Fig. 2.

Neighbor-joining phylogenetic trees based on the 7 MLST genes of two Rec-mas-H strains.


Phylogenetic trees of two Rec-mas-H strains constructed based on the partial sequencing of six MLST genes. (A) argH, (B) cya, (C) glpK, (D) gnd, (E) murC, (F) pta and (G) purH gene sequence based trees were constructed using the neighbor-joining method in the MEGA 4.0 program. The bootstrap values were calculated from 1,000 replications; values <50% are not shown. The black-centered circles indicate that the corresponding clusters were supported by maximum parsimony-based trees. The bar indicates the number of base substitutions per site.

Tab. 2.

Comparison of two Rec-mas-H strains with reference strains of M. abscessus group in sequence similarities of 7 MLST, hsp65 (603 bp) and rpoB (711 bp) gene sequences.

<h2>Comparison of two Rec-mas-H strains with reference strains of <i>M</i>. <i>abscessus</i> group in sequence similarities of 7 MLST, <i>hsp65</i> (603 bp) and <i>rpoB</i> (711 bp) gene sequences.</h2>

With the exception of the purH tree, all of the MLST gene-based trees (argH, cya, glpK, gnd, murC and pta genes) indicated that the two strains were closely related to the M. abscessus subsp. massiliense group, and these topologies were strongly supported by analysis based on the maximum parsimony algorithm (Fig 2). The calculated sequence similarities of all 9 genes (7 MLST genes, the hsp65 gene and the rpoB gene) and their concatenated sequences between the reference strains of the M. abscessus strains and two Rec-mas-H strains are shown in Table 2 and S3 Table.

Phylogenetic analysis of two Rec-mas-H strains using trees based on concatenated sequences

The phylogenetic tree based on the concatenated sequences of the seven MLST genes showed that the two Rec-mas-H strains belong to the M. abscessus subsp. massiliense group and not to M. abscessus subsp. abscessus or M. abscessus subsp. bolletii (Fig 3A) as shown in the rpoB-based analysis (Fig 1B). This suggests that the Rec-mas-H strains may be members of M. abscessus subsp. massiliense rather than M. abscessus subsp. abscessus and that their hsp65 gene may have undergone lateral gene transfer from M. abscessus subsp. abscessus. In the tree based on the concatenated sequences of the seven MLST genes, the two Rec-mas-H strains were closer to the M. abscessus subsp. massiliense Asan 50594 strain than to the M. abscessus subsp. massiliense type strain; this result may be due to the differences in the purH gene sequence (Fig 2G, Table 2 and S3 Table).

<h2>Neighbor-joining phylogenetic trees based on the concatenated sequences of two Rec-mas-H strains.</h2>
Fig. 3.

Neighbor-joining phylogenetic trees based on the concatenated sequences of two Rec-mas-H strains.


Phylogenetic trees based on (A) concatenation of 7 MLST gene sequences, (B) concatenation of 7 MLST gene sequences and the hsp65 gene sequence, and (C) concatenation of 7 MLST gene sequences and the hsp65 and rpoB gene sequences. The trees for all studied strains were generated using the neighbor-joining method. The bootstrap support values (%) from 1,000 replications are indicated for each node; values <50% are not shown. The black-centered circles indicate that the corresponding clusters were supported by maximum parsimony-based trees. The bar indicates the number of base substitutions per site.

Addition of the hsp65 and/or rpoB gene sequences to the 7 MLST concatenated gene sequences did not affect the overall topology of the tree obtained using the 7 MLST concatenated sequences (Fig 3B and 3C). However, the addition of the hsp65 gene sequence slightly affected the sequence similarity between the two Rec-mas-H strains and the M. abscessus subsp. abscessus and M. abscessus subsp. massiliense group strains. The sequence similarity with M. abscessus subsp. abscessus was slightly increased (from 97.6–97.7% to 97.9%), and the similarity with M. abscessus subsp. massiliense group strains was slightly decreased (from 99.5–99.6% to 99.4%) after addition of the hsp65 gene sequence (Table 2).

Separation of 2 Rec-mas-H strains by erm(41) PCR at the subspecies level

The M. abscessus subsp. massiliense erm(41) gene is reported to have a large C-terminal deletion. Therefore, erm(41) PCR can be used as a simple method to differentiate M. abscessus subsp. massiliense from M. abscessus subsp. abscessus and M. abscessus subsp. bolletii species [30, 31]. To further confirm the authenticity of the two Rec-mas-H strains, we applied erm(41) PCR; unlike the M. abscessus subsp. abscessus and M. abscessus subsp. bolletii type strains, which produced a full-sized product (approximately 700 bp), the two Rec-mas-H strains produced a shorter product (approximately 350 bp) that was identical to the product of M. abscessus subsp. massiliense type strain (Fig 4).

<h2>Identification of two Rec-mas-H strains at the subspecies level by PCR targeting the <i>erm</i>(41) gene.</h2>
Fig. 4.

Identification of two Rec-mas-H strains at the subspecies level by PCR targeting the erm(41) gene.


M, 100-bp DNA ladder; Lane 1, M. abscessus subsp. abscessus CIP 104536T; Lane 2, M. abscessus subsp. massiliense CIP 108297T; Lane 3, M. abscessus subsp. bolletii CIP 108541T; Lane 4, Asan 55814; Lane 5, Asan 55262; N, negative control.

Differentiation of 2 Rec-mas-H strains from M. abscessus subsp. massiliense based on restriction patterns of the partial hsp65 gene sequence

PCR-based restriction enzyme and polymorphism analysis (PRA) was performed for further differentiation of the 2 Rec-mas-H strains from M. abscessus subsp. massiliense. Using the two restriction enzymes, BstEII and HaeIII, the partial hsp65 gene sequence was digested [25], and the obtained fragments were compared to those obtained fragments from the partial hsp65 gene sequences of M. abscessus subsp. abscessus and M. abscessus subsp. massiliense by gel electrophoresis. All the BstEII restriction digests showed a similar pattern of fragments 220/245 bp in size. However, HaeIII restriction digests of the 2 Rec-mas-H strains showed patterns identical to those of M. abscessus subsp. abscessus (160/60 bp) but distinct from those of M. abscessus subsp. massiliense (200/60 bp) (Fig 5).

<h2>Differentiation of two Rec-mas-H strains by PCR-restriction enzyme and polymorphism analysis (PRA) of <i>hsp65</i>.</h2>
Fig. 5.

Differentiation of two Rec-mas-H strains by PCR-restriction enzyme and polymorphism analysis (PRA) of hsp65.


Amplified hsp65 gene amplicons were digested with (A) BstEII and (B) HaeIII restriction enzymes. M, 100-bp DNA ladder; Lane 1, M. abscessus subsp. abscessus CIP 104536T; Lane 2, M. abscessus subsp. massiliense CIP 108297T; Lane 3, Asan 55814; Lane 4, Asan 55262.

Discussion

The occurrence of LGT events among species or subspecies in genes encoding chronometer molecules that are used for the diagnosis or identification of pathogenic bacteria may compromise the results obtained when attempting to identify the disease-causing organisms present in infected patients, potentially leading to treatment failure. This problem is especially true in the case of NTM infections, which require a long culture period and often show species- or subspecies-dependent disparities in treatment regimens. Therefore, investigation of LGT events that may affect the diagnosis or identification of target molecules used in the differentiation of mycobacteria is necessary.

The general molecular target for differentiation among species or subspecies in bacterial taxonomy is the 16s rRNA gene [3638]. However, 16s rRNA sequence-based diagnostic and taxonomic methods have some limitations for species differentiation within the genus Mycobacterium, mainly due to the lack of 16S rRNA sequence diversity [3941]. In particular, differentiation of RGM species, including the M. abscessus strains, is almost impossible. Therefore, instead of the 16s rRNA gene, alternative chronometer molecules such as the rpoB and hsp65 genes have been widely used as targets for mycobacterial identification. However, LGT events associated with the rpoB gene have been reported to occur among mycobacteria species or subspecies. For example, the M. yongonense Type I strain carries an rpoBC operon that was laterally transferred from a distantly related species, M. parascrofulaceum [23, 42]. In addition, we recently isolated six M. abscessus subsp. massiliense strains from Korean patients with specific hybrid rpoB genes that were laterally transferred from M. abscessus subsp. abscessus [22]. This suggests that there is a risk of mis-identification when rpoB-based methods are used in mycobacterial diagnosis.

To date, no LGT events associated with another chronometer gene that is used as a mycobacterial target, hsp65, have been reported either in mycobacterial strains or in species within other genera, suggesting that the hsp65 gene is more resistant to LGT than the rpoB gene. However, in this study, we identified for the first time two M. abscessus subsp. massiliense rough strains with hsp65 genes that might have been laterally transferred from M. abscessus subsp. abscessus. These results suggest that the use of hsp65-based diagnosis in mycobacteria also creates a risk of misidentification, at least when attempting to differentiate subspecies within the M. abscessus strains. Indeed, we verified that two hsp65 recombinant M. abscessus subsp. massiliense strains were mis-identified as M. abscessus subsp. abscessus by the hsp65-PRA method targeting the 441-bp Telenti fragment, the most widely used method for NTM differentiation [25]. To the best of our knowledge, this is the first reported case of NTM misidentification by the hsp65-PRA method, and it strongly supports the above notion.

In contrast to the phylogenetic similarity of six recombinant M. abscessus subsp. massiliense strains carrying a hybrid rpoB gene to the smooth M. abscessus subsp. massiliense Type I genotype reported in our previous study [22], our phylogenetic analysis based on MLST sequences indicated that the two M. abscessus subsp. massiliense strains with recombinant hsp65 genes from M. abscessus subsp. abscessus are more closely related to M. abscessus subsp. massiliense Type II than to Type I. The results suggest that these two strains may have descended from the M. abscessus subsp. massiliense Type II genotype. Given that all Type II strains have a rough morphotype due to deletion of the glycopeptidolipid (GPL) gene, the rough colony morphotypes of the 2 hsp65 recombinant strains also support the above notion. However, the two putative recombinant strains showed unique sequences in hsp65, rpoB, glpK, murC, pta and purH genes which were differentiated from M. abscessus subsp. massiliense Type II strain (Figs 1 and 2, S3 Table).

Combined analysis of the results obtained in the present study and in our previous study [22] indicated that LGT events occurred in a total of 8 (3.9%) strains (2 strains with hsp65 recombination and 6 strains with rpoB recombination) in a sample of 206 M. abscessus strains from Korean patients (100 strains in the present study and 106 strains in the recent study). Of note, further MLST analysis showed that all these strains belonged to the M. abscessus subsp. massiliense subspecies and were not M. abscessus subsp. abscessus, suggesting that M. abscessus subsp. massiliense may be more vulnerable to LGT events than M. abscessus subsp. abscessus. This provides a possible explanation of the fact that the genetic and taxonomic diversity among M. abscessus subsp. massiliense strains is higher than that among M. abscessus subsp. abscessus strains.

In conclusion, we identified two M. abscessus subsp. massiliense rough strains from Korean patients with hsp65 genes that might have been laterally transferred from M. abscessus subsp. abscessus. To the best of our knowledge, this is the first report to demonstrate LGT events associated with the hsp65 gene in mycobacteria. This report also suggests that there is potential for misidentification when hsp65-based protocols are used for mycobacterial identification.

Supporting information

S1 Table [docx]
Primer sets used for PCR amplification and sequencing in this study.

S2 Table [docx]
GenBank accession numbers corresponding to the sequences obtained in this study.

S3 Table [docx]
Sequence similarities of the , , and 7 MLST genes and concatenated sequences among . strains.

S1 Fig [docx]
Alignment of the complete gene sequences of . strains and the Asan 55262 strain.


Zdroje

1. Griffith DE, Aksamit T, Brown-Elliott BA, Catanzaro A, Daley C, Gordin F, et al. An official ATS/IDSA statement: diagnosis, treatment, and prevention of nontuberculous mycobacterial diseases. Am J Respir Crit Care Med. 2007;175(4):367–416. doi: 10.1164/rccm.200604-571ST 17277290

2. Daley CL, Griffith DE. Pulmonary disease caused by rapidly growing mycobacteria. Clin Chest Med. 2002;23(3):623–32, vii. 12370998

3. Harris KA, Kenna DT, Blauwendraat C, Hartley JC, Turton JF, Aurora P, et al. Molecular fingerprinting of Mycobacterium abscessus strains in a cohort of pediatric cystic fibrosis patients. J Clin Microbiol. 2012;50(5):1758–61. doi: 10.1128/JCM.00155-12 22403419

4. Maurer FP, Ruegger V, Ritter C, Bloemberg GV, Bottger EC. Acquisition of clarithromycin resistance mutations in the 23S rRNA gene of Mycobacterium abscessus in the presence of inducible erm(41). J Antimicrob Chemother. 2012;67(11):2606–11. doi: 10.1093/jac/dks279 22833642

5. Choi WS, Kim MJ, Park DW, Son SW, Yoon YK, Song T, et al. Clarithromycin and amikacin vs. clarithromycin and moxifloxacin for the treatment of post-acupuncture cutaneous infections due to Mycobacterium abscessus: a prospective observational study. Clin Microbiol Infec. 2011;17(7):1084–90.

6. Nash KA, Brown-Elliott BA, Wallace RJ. A novel gene, erm(41), confers inducible macrolide resistance to clinical isolates of Mycobacterium abscessus but is absent from Mycobacterium chelonae. Antimicrob Agents Chemother. 2009;53(4):1367–76. doi: 10.1128/AAC.01275-08 19171799

7. Kim HY, Kook Y, Yun YJ, Park CG, Lee NY, Shim TS, et al. Proportions of Mycobacterium massiliense and Mycobacterium bolletii strains among Korean Mycobacterium chelonae-Mycobacterium abscessus group isolates. J Clin Microbiol. 2008;46(10):3384–90. doi: 10.1128/JCM.00319-08 18753344

8. Medjahed H, Gaillard JL, Reyrat JM. Mycobacterium abscessus: a new player in the mycobacterial field. Trends Microbiol. 2010;18(3):117–23. doi: 10.1016/j.tim.2009.12.007 20060723

9. Koh WJ, Jeon K, Lee NY, Kim BJ, Kook YH, Lee SH, et al. Clinical Significance of Differentiation of Mycobacterium massiliense from Mycobacterium abscessus. Am J Resp Crit Care. 2011;183(3):405–10.

10. Leao SC, Tortoli E, Euzeby JP, Garcia MJ. Proposal that Mycobacterium massiliense and Mycobacterium bolletii be united and reclassified as Mycobacterium abscessus subsp bolletii comb. nov., designation of Mycobacterium abscessus subsp abscessus subsp nov and emended description of Mycobacterium abscessus. Int J Syst Evol Microbiol. 2011;61:2311–3. doi: 10.1099/ijs.0.023770-0 21037035

11. Tortoli E, Kohl TA, Brown-Elliott BA, Trovato A, Leao SC, Garcia MJ, et al. Emended description of Mycobacterium abscessus, Mycobacterium abscessus subsp. abscessus and Mycobacterium abscessus subsp. bolletii and designation of Mycobacterium abscessus subsp. massiliense comb. nov. Int J Syst Evol Microbiol. 2016;66(11):4471–9. doi: 10.1099/ijsem.0.001376 27499141

12. Kim BJ, Yi SY, Shim TS, Do SY, Yu HK, Park YG, et al. Discovery of a novel hsp65 genotype within Mycobacterium massiliense associated with the rough colony morphology. PLoS One. 2012;7(6):e38420. doi: 10.1371/journal.pone.0038420 22693637

13. Kim BJ, Kim BR, Hong SH, Seok SH, Kook YH. Complete genome sequence of Mycobacterium massiliense clinical strain Asan 50594, belonging to the Type II genotype. Genome Announc. 2013;1(4).

14. Kim BJ, Kim BR, Lee SY, Kook YH, Kim BJ. Rough colony morphology of Mycobacterium massiliense Type II genotype is due to the deletion of glycopeptidolipid locus within its genome. BMC genomics. 2013;14:890. doi: 10.1186/1471-2164-14-890 24341808

15. Gupta RS, Lo B, Son J. Phylogenomics and comparative genomic studies robustly support division of the genus Mycobacterium into an emended genus Mycobacterium and four novel genera. Front Microbiol. 2018;9:67. doi: 10.3389/fmicb.2018.00067 29497402

16. Raz Y, Tannenbaum E. The influence of horizontal gene transfer on the mean fitness of unicellular populations in static environments. Genetics. 2010;185(1):327–37. doi: 10.1534/genetics.109.113613 20194966

17. Ochman H, Lawrence JG, Groisman EA. Lateral gene transfer and the nature of bacterial innovation. Nature. 2000;405(6784):299–304. doi: 10.1038/35012500 10830951

18. Reva O, Korotetskiy I, Ilin A. Role of the horizontal gene exchange in evolution of pathogenic Mycobacteria. BMC Evol Biol. 2015;15. doi: 10.1186/s12862-015-0297-1

19. Krzywinska E, Krzywinski J, Schorey JS. Naturally occurring horizontal gene transfer and homologous recombination in Mycobacterium. Microbiology. 2004;150:1707–12. doi: 10.1099/mic.0.27088-0 15184557

20. Macheras E, Roux AL, Bastian S, Leao SC, Palaci M, Sivadon-Tardy V, et al. Multilocus sequence analysis and rpoB sequencing of Mycobacterium abscessus (Sensu Lato) strains. J Clin Microbiol. 2011;49(2):491–9. doi: 10.1128/JCM.01274-10 21106786

21. Sapriel G, Konjek J, Orgeur M, Bouri L, Frezal L, Roux AL, et al. Genome-wide mosaicism within Mycobacterium abscessus: evolutionary and epidemiological implications. BMC Genomics. 2016;17:118. doi: 10.1186/s12864-016-2448-1 26884275

22. Kim BJ, Kim GN, Kim BR, Shim TS, Kook YH, Kim BJ. Phylogenetic analysis of Mycobacterium massiliense strains having recombinant rpoB gene laterally transferred from Mycobacterium abscessus. PLoS One. 2017;12(6):e0179237. doi: 10.1371/journal.pone.0179237 28604829

23. Kim BJ, Kim BR, Lee SY, Kim GN, Kook YH. Molecular taxonomic evidence for two distinct genotypes of Mycobacterium yongonense via genome-based phylogenetic analysis. PLoS One. 2016;11(3):e0152703. doi: 10.1371/journal.pone.0152703 27031100

24. Kim H, Kim SH, Shim TS, Kim MN, Bai GH, Park YG, et al. Differentiation of Mycobacterium species by analysis of the heat-shock protein 65 gene (hsp65). Int J Syst Evol Microbiol. 2005;55(Pt 4):1649–56. doi: 10.1099/ijs.0.63553-0 16014496

25. Telenti A, Marchesi F, Balz M, Bally F, Bottger EC, Bodmer T. Rapid identification of mycobacteria to the species level by polymerase chain reaction and restriction enzyme analysis. J Clin Microbiol. 1993;31(2):175–8. 8381805

26. Kim BJ, Lee SH, Lyu MA, Kim SJ, Bai GH, Chae GT, et al. Identification of mycobacterial species by comparative sequence analysis of the RNA polymerase gene (rpoB). J Clin Microbiol. 1999;37(6):1714–20. 10325313

27. Adekambi T, Colson P, Drancourt M. rpoB-based identification of nonpigmented and late-pigmenting rapidly growing mycobacteria. J Clin Microbiol. 2003;41(12):5699–708. doi: 10.1128/JCM.41.12.5699-5708.2003 14662964

28. Machado GE, Matsumoto CK, Chimara E, Duarte Rda S, de Freitas D, Palaci M, et al. Multilocus sequence typing scheme versus pulsed-field gel electrophoresis for typing Mycobacterium abscessus isolates. J Clin Microbiol. 2014;52(8):2881–91. doi: 10.1128/JCM.00688-14 24899019

29. Macheras E, Konjek J, Roux AL, Thiberge JM, Bastian S, Leao SC, et al. Multilocus sequence typing scheme for the Mycobacterium abscessus complex. Res Microbiol. 2014;165(2):82–90. doi: 10.1016/j.resmic.2013.12.003 24384536

30. Blauwendraat C, Dixon GLJ, Hartley JC, Foweraker J, Harris KA. The use of a two-gene sequencing approach to accurately distinguish between the species within the Mycobacterium abscessus complex and Mycobacterium chelonae. Eur J Clin Microbiol. 2012;31(8):1847–53.

31. Kim HY, Kim BJ, Kook Y, Yun YJ, Shin JH, Kim BJ, et al. Mycobacterium massiliense is differentiated from Mycobacterium abscessus and Mycobacterium bolletii by erythromycin ribosome methyltransferase gene (erm) and clarithromycin susceptibility patterns. Microbiol Immunol. 2010;54(6):347–53. doi: 10.1111/j.1348-0421.2010.00221.x 20536733

32. Kumar S, Nei M, Dudley J, Tamura K. MEGA: a biologist-centric software for evolutionary analysis of DNA and protein sequences. Brief Bioinform. 2008;9(4):299–306. doi: 10.1093/bib/bbn017 18417537

33. Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987;4(4):406–25. doi: 10.1093/oxfordjournals.molbev.a040454 3447015

34. Fitch WM. Toward defining course of evolution—Minimum change for a specific tree topology. Syst Zool. 1971;20(4):406–16.

35. Felsenstein J. Confidence-limits on phylogenies—an approach using the bootstrap. Evolution. 1985;39(4):783–91. doi: 10.1111/j.1558-5646.1985.tb00420.x 28561359

36. Harmsen D, Karch H. 16S rDNA for diagnosing pathogens: a living tree. Asm News. 2004;70(1):19–24.

37. Kolbert CP, Persing DH. Ribosomal DNA sequencing as a tool for identification of bacterial pathogens. Curr Opin Microbiol. 1999;2(3):299–305. doi: 10.1016/S1369-5274(99)80052-6 10383862

38. Tortoli E. Impact of genotypic studies on mycobacterial taxonomy: the new mycobacteria of the 1990s. Clin Microbiol Rev. 2003;16(2):319–54. doi: 10.1128/CMR.16.2.319-354.2003 12692101

39. Clarridge JE. Impact of 16S rRNA gene sequence analysis for identification of bacteria on clinical microbiology and infectious diseases. Clin Microbiol Rev. 2004;17(4):840–62. doi: 10.1128/CMR.17.4.840-862.2004 15489351

40. Tortoli E. Phylogeny of the genus Mycobacterium: many doubts, few certainties. Infect Genet Evol. 2012;12(4):827–31. doi: 10.1016/j.meegid.2011.05.025 21684354

41. Turenne CY, Tschetter L, Wolfe J, Kabani A. Necessity of quality-controlled 16S rRNA gene sequence databases: identifying nontuberculous Mycobacterium species. J Clin Microbiol. 2001;39(10):3637–48. doi: 10.1128/JCM.39.10.3637-3648.2001 11574585

42. Kim BJ, Hong SH, Kook YH. Molecular evidence of lateral gene transfer in rpoB gene of Mycobacterium yongonense strains via multilocus sequence analysis. PLoS One. 2013;8(1):e51846. doi: 10.1371/journal.pone.0051846 23382812


Článek vyšel v časopise

PLOS One


2019 Číslo 9
Nejčtenější tento týden
Nejčtenější v tomto čísle
Kurzy Podcasty Doporučená témata Časopisy
Přihlášení
Zapomenuté heslo

Zadejte e-mailovou adresu, se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.

Přihlášení

Nemáte účet?  Registrujte se

#ADS_BOTTOM_SCRIPTS#