|Year : 2017 | Volume
| Issue : 2 | Page : 138-141
Molecular typing of Mycobacterium Abscessus isolated from cystic fibrosis patients
Alberto Trovato1, Rossella Baldan1, Danila Costa2, Tullia M Simonetti3, Daniela M Cirillo1, Enrico Tortoli1
1 Emerging Bacterial Pathogens Unit, IRCCS San Raffaele Scientific Institute, Milan, Italy
2 Microbiology Unit, Policlinico University Hospital, Bari, Italy
3 Microbiology and Virology Unit, Careggi University Hospital, Florence, Italy
|Date of Web Publication||19-May-2017|
Emerging Bacterial Pathogens Unit, IRCCS San Raffaele Scientific Institute, Milan
Source of Support: None, Conflict of Interest: None
Background: The possibility of inter-human transmission of Mycobacterium abscessus in cystic fibrosis centres has been recently hypothesized suggesting the need for the molecular characterization of strains isolated from such patients. Materials and Methods: One hundred and forty one isolates of M. abscessus grown from sputum samples of 29 patients with cystic fibrosis were genotyped resorting to variable number of tandem repeats (VNTR) determination and whole genome sequencing (WGS). Results: Out of 29 VNTR profiles, 15 were unique to the same number of patients while seven were shared by multiple patients. WGS showed that only two of the patients sharing common VNTR patterns were indeed infected by the same strain. The shared VNTR patterns were mostly present among the isolates of M. abscessus subsp. abscessus. Conclusion: As expected WGS showed a clearly higher discriminatory power in comparison with VNTR and appeared the only molecular epidemiology tool suitable to effectively discriminate the isolates of M. abscessus subsp. abscessus.
Keywords: Cystic fibrosis, Mycobacterium abscessus, VNRT, whole genome sequencing
|How to cite this article:|
Trovato A, Baldan R, Costa D, Simonetti TM, Cirillo DM, Tortoli E. Molecular typing of Mycobacterium Abscessus isolated from cystic fibrosis patients. Int J Mycobacteriol 2017;6:138-41
|How to cite this URL:|
Trovato A, Baldan R, Costa D, Simonetti TM, Cirillo DM, Tortoli E. Molecular typing of Mycobacterium Abscessus isolated from cystic fibrosis patients. Int J Mycobacteriol [serial online] 2017 [cited 2018 Apr 21];6:138-41. Available from: http://www.ijmyco.org/text.asp?2017/6/2/138/206601
| Introduction|| |
Mycobacterium abscessus is, among the rapidly growing mycobacteria, the species most frequently involved in human infections. In adjunct to the lung, cutis, soft tissues, bone, and joints may be affected. As other nontuberculous mycobacteria, it is present in the environment, is oligotrophic, resistant to chlorination, and is characterized by prolonged survival in biofilms. Its isolation from drinking water has been repeatedly reported., Environmental aerosols are considered the major source of contagion for respiratory infections, but the interhuman transmission has been recently hypothesized as well.M. abscessus is intrinsically resistant to most antimicrobials, the infections are, therefore, hardly treatable, and the eradication cannot be easily achieved.
Cystic fibrosis (CF) patients are particularly susceptible to pulmonary infections caused by mycobacteria and expecially by M. abscessus. In such patients, a prevalence higher than 10% is estimated worldwide. The species is split into three subspecies, i.e., M. abscessus subsp. abscessus (Masa), M. abscessus subsp. bolletii (Masb), and M. abscessus subsp. massiliense (Masm). This taxonomic distinction is also clinically relevant as Masa and Masb have a functional erm (41) gene  which confers inducible macrolide resistance; this gene is instead truncated in Masm which, consequently, lacks inducible resistance.
In the last years, possible outbreaks of M. abscessus in several CF centers have been reported,, suggesting the need for the molecular characterization of strains isolated from such patients.
Different molecular fingerprinting techniques have been used for M. abscessus strains.
Pulsed-field gel electrophoresis of the digested bacterial chromosome yields reproducible patterns of DNA fragments allowing strains differentiation within M. abscessus subspecies., The technique is, however, prone to DNA degradation leading to uninterpretable results.
Multilocus sequence typing is based on the analysis of polymorphisms in the sequences of selected housekeeping genes. The discriminative power is proportional to the number of genes investigated and is further increased when the internal transcribed spacer between 16S and 23S is sequenced as well.
In the last years, the analysis of variable number of tandem repeats (VNTR) present at specific loci has been adapted to M. abscessus, on the wake of the excellent experience accumulated with Mycobacterium tuberculosis. The discriminative potential was at least comparable to that of techniques mentioned above.,,,
Very recently, however, high-throughput sequencing has opened up new horizons, and the comparison of whole genomes, unimaginable before, has become reality. It is, therefore, expected that whole genome sequencing (WGS), thanks to its ultimate potential, will soon become the standard for molecular epidemiology investigations.
The aim of the present study was to reassess with WGS, a cohort of M. abscessus isolated from CF patients previously genotyped with VNTR.
| Materials and Methods|| |
A total of 141 isolates of M. abscessus grown from sputum samples of 29 patients with CF had been identified at subspecies level by sequencing of rpoB gene.
To perform VNTR analysis, we selected 13 M. abscessus loci, characterized by high allelic diversity, stability, and reproducibility, among the ones investigated in previous studies.,,, Loci with repeat size >25 and <50 bp and characterized by high discriminatory power, as determined basing on Hunter–Gaston index, were privileged [Table 1]. PCR amplification and electrophoretic estimation of amplicon size were carried out as reported before.
|Table 1: VNTR loci targeted in previous, and present, studies. In each row are reported the labels used by different authors for each single repeat|
Click here to view
For WGS, we used the Illumina MiSeq sequencer, Nextera XT library preparation kits, and MiSeq reagent kits, following manufacturer's instructions (Illumina, USA). The reads were mapped to the reference genome of M. abscessus ATCC 19977 (GenBank ID: NC_010397.1) with the alignment program BWA, and mappings were refined with the GATK  and SAMtools toolkits. For variant detection in mapped reads were employed custom perl scripts with thresholds of a minimum coverage of four reads in both forward and reverse orientation, four reads calling the allele with at least a Phred score of 20 and 75% allele frequency. Single nucleotide polymorphism (SNP) positions with a reliable base call in at least 95% of the isolates were concatenated to a sequence alignment, excluding SNPs within a window of 12 bp from each other. From the aligned sequences of concatenated SNPs, a maximum parsimony tree was built on 1105001 SNP positions. The isolates fully sequenced were 95; each patient had at least one genome sequenced. For patients with multiple isolates, a proportion of them, including at least the first and the last one, was selected for WGS.
| Results|| |
Seven patients had a single isolate; the remaining 22 had multiple isolates, ranging from 2 to 19. Out of 22 patients with multiple isolates, 19 were persistently infected by a single subspecies: Masa in 13, Masb in 5, and Masm in 1. Three patients yielded isolates of two subspecies: Masm and Masa in two of them, Masb and Masa in the other.
VNTR analysis produced 29 different profiles, 15 patients had the same number of unique VNTR patterns, whereas the remaining five profiles were presented by more than one patient each [Table 2]. The variability of the number of repeats was significantly higher in Masa, with an average of 5.7 repeats/locus, in comparison with Masb (2.4) and Masm (1.6). Masb had steadily one repeat at locus Masb1; all the strains of Masm had one repeat at loci Mab1 and TR155 and two repeats at loci Mab11, Mab14, and Mab21.
In 18 out of 19 patients with multiple isolates belonging to the same subspecies, the VNTR profile remained unchanged over time suggesting persistent infection by a single strain; in the latter patient, a dominant strain, as indicated by five isolates with steady VNTR pattern, was occasionally overgrown by two isolates, each characterized by a unique profile. Interestingly, two patients (BN and NO) shared a VNTR profile presenting the signature of Masm despite the identification as Masa on the basis of the rpoB sequence.
The phylogenetic tree [Figure 1] showed three clearly distinct branches corresponding to the three M. abscessus subspecies; in it the strains of the patients above (BN and NO) belonged to the branch of Masm. The cluster corresponding to the subspecies Masa presented 35 closely related genomes [Figure 1], dotted frame] and 27 dispersed.
|Figure 1: Phylogenetic tree including 95 clinical isolates (plus 3 reference strains) of M. abscessus. White circles indicate M. abscessus subsp. massiliense; gray circles indicate M. abscessus subsp. abscessus; black circles indicates M. abscessus subsp. bolletii; the dotted boxes highlight the dominating clone of M. abscessus subsp. abscessus. Numbers on branches indicate the number of distinct single nucleotide polymorphism positions between isolates. Same-patients isolates are indicated by unique two-letters identifiers followed by a number featuring the isolation order.|
Click here to view
The same patient genomes differed for a number of SNPs ranging from 0 to 21 in 95% of patients; we considered, therefore, thirty SNPs a reasonable cutoff to distinguish isolates of the same strain from those of independent strains. On such basis, two patients (BN and NO) turned out infected by the same strain (difference 24–29 SNPs). WGS analysis showed furthermore that the six patients sharing the same VNTR profile [Table 2] were instead infected by independent, although closely related, strains. They differed in fact for 33–92 SNPs.
| Discussion and Conclusions|| |
As expected, in our comparison, the WGS revealed more discriminative than 13 loci VNTR. The subspecies Masa was the most affected by this deficiency of VNTR. In fact, a large number of Masa strains belonged to one dominating clone circulating worldwide. Out of 11 patients infected by strains belonging to this clone, 7 had identical VNTR genotype [Table 2] and the others (GA, CR, GL, and TR) showed only one locus of difference. It is commonly accepted that VNTR patterns differing for the repeat number at one locus only can be regarded as identical, and we have indeed observed one locus difference in one patient that we considered persistently infected by a single strain. Applying this rule, the number of patients infected by strains indistinguishable by our VNTR analysis, but demonstrated independent by WGS, would increase to 12, corresponding to 75% of the ones infected by Masa. Data excluding any epidemiological link among them further support WGS results. The only true case of transmission concerned a strain of Masm isolated from two patients pertaining to the same CF center and was correctly identified by both VNTR and WGS analysis.
In conclusion, we have observed that SNPs analysis has a higher discriminatory capacity compared to VNTR, in particular for Masa strains belonging to dominant circulating clusters. For other two M. abscessus subspecies, which however included here a limited number of isolates, no significant difference emerged between the two molecular approaches. Both methods revealed more accurate than rpoB gene sequencing in differentiating the strains at subspecies level.
Financial support and sponsorship
This study was supported by research grant FFC #27/2014 from Fondazione Ricerca Fibrosi Cistica.
Conflicts of interest
There are no conflicts of interest.
| References|| |
Tortoli E. Clinical manifestations of nontuberculous mycobacteria infections. Clin Microbiol Infect 2009;15:906-10.
Falkinham JO 3rd
. Ecology of nontuberculous mycobacteria – Where do human infections come from? Semin Respir Crit Care Med 2013;34:95-102.
Thomson R, Tolson C, Sidjabat H, Huygens F, Hargreaves M. Mycobacterium abscessus
isolated from municipal water – A potential source of human infection. BMC Infect Dis 2013;13:241.
Donohue MJ, Mistry JH, Donohue JM, O'Connell K, King D, Byran J, et al.
Increased frequency of nontuberculous mycobacteria detection at potable water taps within the United States. Environ Sci Technol 2015;49:6127-33.
Bryant JM, Grogono DM, Rodriguez-Rincon D, Everall I, Brown KP, Moreno P, et al.
Emergence and spread of a human-transmissible multidrug-resistant nontuberculous mycobacterium. Science 2016;354:751-7.
Benwill JL, Wallace RJ Jr. Mycobacterium abscessus
: Challenges in diagnosis and treatment. Curr Opin Infect Dis 2014;27:506-10.
van Dorn A. Multidrug-resistant Mycobacterium abscessus
threatens patients with cystic fibrosis. Lancet Respir Med 2017;5:15.
Seddon P, Fidler K, Raman S, Wyatt H, Ruiz G, Elston C, et al.
Prevalence of nontuberculous mycobacteria in cystic fibrosis clinics, United Kingdom, 2009. Emerg Infect Dis 2013;19:1128-30.
Tortoli E, Kohl TA, Brown-Elliott BA, Trovato A, Leão SC, Garcia MJ, et al.
Emended description of Mycobacterium abscessus
, Mycobacterium abscessus
subsp. abscessus and Mycobacterium abscessus
and designation of Mycobacterium abscessus
subsp. massiliense comb. nov. Int J Syst Evol Microbiol 2016;66:4471-9.
Nash KA, Brown-Elliott BA, Wallace RJ Jr. 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:1367-76.
Bryant JM, Grogono DM, Greaves D, Foweraker J, Roddick I, Inns T, et al.
Whole-genome sequencing to identify transmission of Mycobacterium abscessus
between patients with cystic fibrosis: A retrospective cohort study. Lancet 2013;381:1551-60.
Jönsson BE, Gilljam M, Lindblad A, Ridell M, Wold AE, Welinder-Olsson C. Molecular epidemiology of Mycobacterium abscessus
, with focus on cystic fibrosis. J Clin Microbiol 2007;45:1497-504.
Matsumoto CK, Chimara E, Bombarda S, Duarte RS, Leão SC. Diversity of pulsed-field gel electrophoresis patterns of Mycobacterium abscessus
type 2 clinical isolates. J Clin Microbiol 2011;49:62-8.
Harris KA, Kenna DT. Mycobacterium abscessus
infection in cystic fibrosis: Molecular typing and clinical outcomes. J Med Microbiol 2014;63(Pt 10):1241-6.
Kim SY, Kang YA, Bae IK, Yim JJ, Park MS, Kim YS, et al.
Standardization of multilocus sequence typing scheme for Mycobacterium abscessus
and Mycobacterium massiliense
. Diagn Microbiol Infect Dis 2013;77:143-9.
Sassi M, Ben Kahla I, Drancourt M. Mycobacterium abscessus
multispacer sequence typing. BMC Microbiol 2013;13:3.
Choi GE, Chang CL, Whang J, Kim HJ, Kwon OJ, Koh WJ, et al.
Efficient differentiation of Mycobacterium abscessus
complex isolates to the species level by a novel PCR-based variable-number tandem-repeat assay. J Clin Microbiol 2011;49:1107-9.
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:1758-61.
Wong YL, Ong CS, Ngeow YF. Molecular typing of Mycobacterium abscessus
based on tandem-repeat polymorphism. J Clin Microbiol 2012;50:3084-8.
Shin SJ, Choi GE, Cho SN, Woo SY, Jeong BH, Jeon K, et al.
Mycobacterial genotypes are associated with clinical manifestation and progression of lung disease caused by Mycobacterium abscessus
and Mycobacterium massiliense
. Clin Infect Dis 2013;57:32-9.
Adékambi T, Colson P, Drancourt M. rpoB
-based identification of nonpigmented and late-pigmenting rapidly growing mycobacteria. J Clin Microbiol 2003;41:5699-708.
Hunter PR, Gaston MA. Numerical index of the discriminatory ability of typing systems: An application of Simpson's index of diversity. J Clin Microbiol 1988;26:2465-6.
Li H, Durbin R. Fast and accurate short read alignment with burrows-wheeler transform. Bioinformatics 2009;25:1754-60.
McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al.
The Genome Analysis Toolkit: A MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res 2010;20:1297-303.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al.
The sequence alignment/map format and SAMtools. Bioinformatics 2009;25:2078-9.
[Table 1], [Table 2]