# Investigation of TwoMycobacterium abscessus Outbreaks in Quebec Using Whole Genome Sequencing

Authors: Michelle Wuzinski; Hafid Soualhine; Emilie Valliere; Pierre-Marie Akochy; Nancy Cloutier; Aaron Petkau; Brynn Kaplen; Maryse Duchesne; Jasmin Villeneuve; Meenu K. Sharma
Journal: BioMed Research International
 (2020)

---

## Abstract

In recent decades, nontuberculous mycobacteria (NTM) infections are of emerging public health concern and have contributed towards significant clinical and economic burden globally. One such rapid growing mycobacteria,Mycobacterium abscessus, can cause clonal outbreaks, and these bacteria exhibit a highly resistant antimicrobial susceptibility profile. Here, we present an investigation of two small outbreaks of M. abscessus: first in a pediatric clinic setting and second in a tattoo parlour from Quebec. Two whole genome sequencing approaches were utilized for genotyping: MAB-MLST, a multilocus sequencing typing scheme containing housekeeping, identification, and antimicrobial resistance genes, and SNVPhyl that uses phylogenetics to determine single nucleotide variations between strains. MAB-MLST results showed that the pediatric outbreak strains had two distinct sequence types, demonstrating that one strain did not belong to the outbreak, while all tattoo outbreak isolates belonged to the same sequence type. SNVPhyl results were similar to MAB-MLST results and showed that the pediatric outbreak strains tightly clustered together with 0-1 SNVs between isolates, a sharp contrast between unrelated strains used as controls. Similar results were seen for tattoo outbreak cases with 3-11 SNVs between isolates. NTM infections can be difficult to identify, and outbreak investigations can be complicated. Thus, WGS tools can be used in public health outbreak investigations as they provide high discriminatory power.


---

## Body

## 1. Introduction

Globally, the incidence of nontuberculous mycobacteria (NTMs) includingMycobacterium abscessus infections has been on the rise [1]. NTMs are ubiquitous environmental organisms, and they can be found in air, soil, and water. These mycobacteria can colonize health-care facilities and commercial and public settings and hence can be responsible for numerous small outbreaks [2, 3]. Clinical presentation of one such mycobacteria, M. abscessus, can vary depending on site of infection, but the disease presentation often occurs as soft tissue or pulmonary infections, particularly in those with underlying lung disease such as cystic fibrosis or bronchiectasis. An added disadvantage of M. abscessus infection is its resistance to a large number of common antimycobacterial antibiotics. It is referred to as one of the most resistant bacteria known, and thus, treatment is challenging [4]. M. abscessus is divided into three subspecies, M. abscessus subsp. abscessus, M. abscessus subsp. bolletii, and M. abscessus subsp. massiliense [5, 6] which have different susceptibilities to antibiotics, and thus, differentiation between them is of vital importance for treatment of infection.M. abscessus outbreaks can also result from invasive procedures such as cosmetic surgery, tattooing, and acupuncture. Soft tissue infections and surgeries can be responsible for M. abscessus infection via contamination of disinfectants, saline or other fluids, and medical tools [7–10]. Investigation of these outbreaks often involve cumbersome and time consuming strategies such as specimen tracking, laboratory quality control record checks, or via tests for identification and genotyping such as PCR amplification of 16S, rpoB or hsp65, random amplification of polymorphic DNA (RAPD), Diversilab rep-PCR genotyping, multilocus sequence typing (MLST), or DNA restriction fragment length polymorphism (RFLP) [11–13].Nowadays, outbreak investigations have been revolutionized by the implementation of whole genome sequencing (WGS). WGS provides discrimination that was not possible before its advent and hence makes it a powerful tool for not only identification and outbreak investigation but also for surveillance and diagnostics. As a surveillance tool, it adds to contact tracing information and provides evidence-based link of transmission of mycobacteria. Single nucleotide variant (SNV) phylogenomic analysis, such as SNVPhyl [14], can be utilized to construct a phylogenetic tree of microbial genome collection, in order to cluster and compare genomes and infer relatedness. On the other hand, MAB-MLST (M. abscessus-multilocus sequencing typing developed from our laboratory) is an expansion of a traditional MLST scheme to include antimicrobial resistance, identification, and housekeeping genes to generate allele profiles, called sequence types (ST) [15]. Here, we present the use of two whole genome sequencing-based analysis tools, SNVPhyl and MAB-MLST, to retrospectively investigate two suspected outbreaks in the province of Quebec. In order to delineate the transmission of these outbreaks, the whole genome sequencing data for presumptively linked cases stemming from a tattoo parlour and another separate outbreak from a pediatric clinic was ran against both tools.


## 2. Materials and Methods

### 2.1. Study Strains

M. abscessus isolates were obtained from pediatric patients with otitis of the middle ear and from clients of a tattoo parlour presenting with skin infections. The clinical presentation was skin lesions in the form of multiple 1-2 mm pustules or papules at the tattoo site. Clinical specimens for culture were obtained from patients as described in Table 1. Clinical specimens were decontaminated using standard methodology as described by Clinical and Laboratory Standards Institute [16]. Sediments were cultured to Mycobacterial Growth Indicator Tubes (MGIT) (BD Bactec, UK), Lowenstein Jensen medium, and Middlebrook 7H10 (MK) agar plates and incubated at 37°C. The blood agar plate was used to check for contamination at both 24 hours and 48 hours. Culture was grown on MK plates for 4-7 days, evaluated for purity, before being used for DNA extraction. Culture identification was performed by Laboratoire de Santé Publique du Québec using 16S and hsp65 DNA sequencing.Table 1
Case information of outbreak isolates.
Case #
NRCM #
Age
Gender
Outbreak
Specimen type
Site
A
1801480
3
F
Pediatric
Pus
Ear
B
1801479
6
M
Pediatric
Aspiration
Ear
C
1801477
4
M
Pediatric
Swab
Ear
D
1801478
1
F
Pediatric
Pus
Ear
E
1801481
60
F
NR∗
Sputum
Lung
F
1801485
52
M
NR
Sputum
Lung
G
1801482
20
M
NR
Sputum
Lung
H
1801483
23
F
Tattoo
Biopsy
Skin
I
1801484
30
M
Tattoo
Biopsy
Skin
J
1801486
25
F
Tattoo
Biopsy
Skin
K
2000565
23
F
Tattoo
Biopsy
Skin
∗NR: nonrelated case.Antimicrobial susceptibility testing (AST) was performed using Trek Sensititre broth microdilution method (Sensititre™; Thermo Fisher Scientific, Cleveland, OH, USA) using the RAPMYCO plates, and minimum inhibitory concentrations were documented for recommended drugs as per manufacturer’s instructions. The interpretation criteria were based on CLSI breakpoints [17]. Isolates were shipped to the National Reference Centre for Mycobacteriology for susceptibility testing and whole genome sequencing.


### 2.2. Genomic DNA Extraction

A loopful of solid growth from MK plate was added to a 2 mL tube containing sterile water with glass beads and vortexed. This resulting suspension was adjusted to be 107 CFU/mL by McFarland, and CFUs were confirmed by colony counts on MK agar plates. An aliquot of 1.5 mL of the cell suspension was centrifuged for 30 minutes at 4000 × g. The supernatant was removed, and cell pellets were frozen at -80°C until ready for DNA extraction.M. abscessus DNA extractions were performed using the Qiagen DNeasy UltraClean Microbial Kit (Hilden, Germany). Cell pellets were resuspended in 300 μL of PowerBead Solution, and then, contents were transferred to a PowerBead tube. 50 μL of Solution SL was added and bead beated using the Fast-Prep 24 Tissue and Cell Homogenizer (Hyland Scientific, WA) at 6.5 m/s for 4 cycles of 45 seconds each. Manufacturer’s instructions were followed thereafter except the final elution step was done twice (50μL×2) to bring the final elution volume to 100 μL. Samples were then stored at 4°C until further processing.All genomic DNA extracts were quantified using the Qubit 3 Fluorometer (Thermofisher Scientific, CA, USA) as per manufacturer instructions. Samples were diluted to 0.5 ng/μL, and 4 μL of DNA was used in the sequencing reaction.


### 2.3. Whole Genome Sequencing

Sequencing was performed in house using the NexteraXT DNA Library Preparation Kit (Illumina, CA, USA) and sequenced on the Illumina MiSeq (Illumina, CA, USA) with the v3 600 cycle kit (2×300 bp reads). The sequences have been submitted to GenBank under Bioproject number PRJNA656566. Quality of the sequences was determined using FastQC (http://bioinformatics.babraham.ac.uk/projects/fastqc) and SMALT (http://sanger.ac.uk/tool/smalt-0) to map reads to a reference genome. Sequences that mapped at least 85% of the reference strain (M. abscessus ATCC 19977) and had an average base coverage of at least 50x were included for analysis. Sequences were assembled in the Integrated Rapid Infectious Disease Analysis Project (IRIDA), an open source platform for public health genomics, using the Assembly and Annotation pipeline, which uses Shovill/SPAdes (http://github.com/tseemann/shovill) for genome assembly [18, 19].All sequences were subjected to phylogenetic analysis using the SNVPhyl pipeline in IRIDA withM. abscessus subsp. abscessus ATCC19977 as the reference genome to generate single nucleotide variant (SNV) matrices and phylogenetic tree [14]. Analysis parameters used were minimum SNV coverage of 15, minimum mean mapping quality score of 30, and SNV abundance ratio 0.75. The SNV density filter was disabled to avoid removal of most of the SNVs for the more distantly related genomes. Second, assembled genomes were genotyped as previously described with MAB-MLST [15] which can be run through the Galaxy platform [20]. Identification, antimicrobial susceptibility testing, and MAB-MLST sequence types for the outbreak strains were analyzed (Table S1 and S2). Statistical comparison was conducted between genotyping and phylogenetic analysis, and any discrepancies were identified as previously described [15].


## 2.1. Study Strains

M. abscessus isolates were obtained from pediatric patients with otitis of the middle ear and from clients of a tattoo parlour presenting with skin infections. The clinical presentation was skin lesions in the form of multiple 1-2 mm pustules or papules at the tattoo site. Clinical specimens for culture were obtained from patients as described in Table 1. Clinical specimens were decontaminated using standard methodology as described by Clinical and Laboratory Standards Institute [16]. Sediments were cultured to Mycobacterial Growth Indicator Tubes (MGIT) (BD Bactec, UK), Lowenstein Jensen medium, and Middlebrook 7H10 (MK) agar plates and incubated at 37°C. The blood agar plate was used to check for contamination at both 24 hours and 48 hours. Culture was grown on MK plates for 4-7 days, evaluated for purity, before being used for DNA extraction. Culture identification was performed by Laboratoire de Santé Publique du Québec using 16S and hsp65 DNA sequencing.Table 1
Case information of outbreak isolates.
Case #
NRCM #
Age
Gender
Outbreak
Specimen type
Site
A
1801480
3
F
Pediatric
Pus
Ear
B
1801479
6
M
Pediatric
Aspiration
Ear
C
1801477
4
M
Pediatric
Swab
Ear
D
1801478
1
F
Pediatric
Pus
Ear
E
1801481
60
F
NR∗
Sputum
Lung
F
1801485
52
M
NR
Sputum
Lung
G
1801482
20
M
NR
Sputum
Lung
H
1801483
23
F
Tattoo
Biopsy
Skin
I
1801484
30
M
Tattoo
Biopsy
Skin
J
1801486
25
F
Tattoo
Biopsy
Skin
K
2000565
23
F
Tattoo
Biopsy
Skin
∗NR: nonrelated case.Antimicrobial susceptibility testing (AST) was performed using Trek Sensititre broth microdilution method (Sensititre™; Thermo Fisher Scientific, Cleveland, OH, USA) using the RAPMYCO plates, and minimum inhibitory concentrations were documented for recommended drugs as per manufacturer’s instructions. The interpretation criteria were based on CLSI breakpoints [17]. Isolates were shipped to the National Reference Centre for Mycobacteriology for susceptibility testing and whole genome sequencing.


## 2.2. Genomic DNA Extraction

A loopful of solid growth from MK plate was added to a 2 mL tube containing sterile water with glass beads and vortexed. This resulting suspension was adjusted to be 107 CFU/mL by McFarland, and CFUs were confirmed by colony counts on MK agar plates. An aliquot of 1.5 mL of the cell suspension was centrifuged for 30 minutes at 4000 × g. The supernatant was removed, and cell pellets were frozen at -80°C until ready for DNA extraction.M. abscessus DNA extractions were performed using the Qiagen DNeasy UltraClean Microbial Kit (Hilden, Germany). Cell pellets were resuspended in 300 μL of PowerBead Solution, and then, contents were transferred to a PowerBead tube. 50 μL of Solution SL was added and bead beated using the Fast-Prep 24 Tissue and Cell Homogenizer (Hyland Scientific, WA) at 6.5 m/s for 4 cycles of 45 seconds each. Manufacturer’s instructions were followed thereafter except the final elution step was done twice (50μL×2) to bring the final elution volume to 100 μL. Samples were then stored at 4°C until further processing.All genomic DNA extracts were quantified using the Qubit 3 Fluorometer (Thermofisher Scientific, CA, USA) as per manufacturer instructions. Samples were diluted to 0.5 ng/μL, and 4 μL of DNA was used in the sequencing reaction.


## 2.3. Whole Genome Sequencing

Sequencing was performed in house using the NexteraXT DNA Library Preparation Kit (Illumina, CA, USA) and sequenced on the Illumina MiSeq (Illumina, CA, USA) with the v3 600 cycle kit (2×300 bp reads). The sequences have been submitted to GenBank under Bioproject number PRJNA656566. Quality of the sequences was determined using FastQC (http://bioinformatics.babraham.ac.uk/projects/fastqc) and SMALT (http://sanger.ac.uk/tool/smalt-0) to map reads to a reference genome. Sequences that mapped at least 85% of the reference strain (M. abscessus ATCC 19977) and had an average base coverage of at least 50x were included for analysis. Sequences were assembled in the Integrated Rapid Infectious Disease Analysis Project (IRIDA), an open source platform for public health genomics, using the Assembly and Annotation pipeline, which uses Shovill/SPAdes (http://github.com/tseemann/shovill) for genome assembly [18, 19].All sequences were subjected to phylogenetic analysis using the SNVPhyl pipeline in IRIDA withM. abscessus subsp. abscessus ATCC19977 as the reference genome to generate single nucleotide variant (SNV) matrices and phylogenetic tree [14]. Analysis parameters used were minimum SNV coverage of 15, minimum mean mapping quality score of 30, and SNV abundance ratio 0.75. The SNV density filter was disabled to avoid removal of most of the SNVs for the more distantly related genomes. Second, assembled genomes were genotyped as previously described with MAB-MLST [15] which can be run through the Galaxy platform [20]. Identification, antimicrobial susceptibility testing, and MAB-MLST sequence types for the outbreak strains were analyzed (Table S1 and S2). Statistical comparison was conducted between genotyping and phylogenetic analysis, and any discrepancies were identified as previously described [15].


## 3. Results

Sequences that met quality parameters described in methods were included for bioinformatics analysis. Percent mapping of sequences to the reference (M. abscessus ATCC19977) genome ranged from 86.2% to 95.7%. Average sequence coverage values for all isolates ranged from 90x to 151x.


### 3.1. Pediatric Outbreak

Outbreak isolates B, C, and D were identified asM. abscessus subsp. massiliense both by MAB-MLST (using hsp65, erm(41), and rpoB alleles) and traditional PCR sequencing of 16S and hsp65. Isolate A and a nonrelated sample E were identified as M. abscessus subsp. abscessus. A total of 126299 valid positions (of 138733 total identified positions) were used to generate the SNVPhyl phylogeny shown in Figure 1. The number of single nucleotide differences between pediatric outbreak isolates is shown in Table 2. Isolates B, C, and D clustered closely in the phylogenetic tree, belonged to sequence type 100, and had 0-1 SNVs. A and E did not cluster with the outbreak isolates and belonged to sequence types 102 and 101, respectively. Four percent of the total identified positions were filtered out due to coverage.Figure 1
Phylogenetic analysis of isolates of suspected outbreak originating from a pediatric clinic. Tree was generated as output from the SNVPhyl pipeline run through IRIDA and visualized in MicroReact. The reference genome refers to ATCC19977. Default SNVPhyl parameters were used with a minimum coverage depth of 15 and mean mapping quality score of 30. Timeline shows date of sample isolation.Table 2
Matrix of single nucleotide differences between pediatric outbreak isolates.
Strain
A
B
C
D
E
Reference
A
0
108610
108610
108611
28650
29500
B
108610
0
0
1
107181
107907
C
108610
0
0
1
107181
107907
D
108611
1
1
0
107182
107908
E
28650
107181
107181
107182
0
19259
Reference
29500
107907
107907
107908
19259
0
Single nucleotide variations between study isolates as determined by SNVPhyl. Reference isM. abscessus ATCC19977.


### 3.2. Tattoo Outbreak

All tattoo isolates were identified asM. abscessus subsp. abscessus by MAB-MLST and PCR sequencing of 16S and hsp65. In addition to the four suspected outbreak cases (isolates H, I, J, K) thought to have originated from the tattoo parlour, unrelated outliers F and G were included in the analysis. Isolates H, I, J, and K clustered closely and belonged to sequence type 103. Isolates F and G did not cluster with the outbreak isolates and belonged to sequence type 1 and 83, respectively. The phylogenetic tree presented in Figure 2 was built from 46464 valid positions identified by SNVPhyl (of 50803 total identified positions). Eight percent of total identified positions were filtered out due to coverage or otherwise invalid. As seen in Table 3, the outbreak isolates have between 3 and 11 SNVs between each other in contrast to the outlier which has between 45933 and 45934 SNVs from the outbreak isolates (Table 3). Further, all suspected outbreak cases were assigned the sequence type (ST) as seen in Table S1 and S2. Drug susceptibility profile for isolate K was susceptible to clarithromycin with MIC of 2 μg/mL, while the related isolates were intermediary to clarithromycin with MIC of 4 μg/mL (Table S1).Figure 2
Phylogenetic tree of isolates belonging to the tattoo parlour outbreak. Tree generated as output from the SNVPhyl pipeline run through IRIDA and visualized in MicroReact. The reference genome refers to ATCC19977. Default parameters used with minimum coverage depth of 15 and minimum mean mapping quality score of 30. Timeline shows month sample isolation.Table 3
SNV matrix of tattoo isolates representing nucleotide differences between isolates.
Strain
H
I
J
K
G
F
Reference
H
0
3
4
7
45715
45909
45701
I
3
0
3
10
45714
45908
45700
J
4
3
0
11
45715
45909
45701
K
7
10
11
0
45714
45908
45700
G
45715
45714
45715
45714
0
1270
36
F
45909
45908
45909
45908
1270
0
1256
Reference
45701
45700
45701
45700
36
1256
0
Single nucleotide variations between study isolates as determined by SNVPhyl. Reference isM. abscessus ATCC19977.


## 3.1. Pediatric Outbreak

Outbreak isolates B, C, and D were identified asM. abscessus subsp. massiliense both by MAB-MLST (using hsp65, erm(41), and rpoB alleles) and traditional PCR sequencing of 16S and hsp65. Isolate A and a nonrelated sample E were identified as M. abscessus subsp. abscessus. A total of 126299 valid positions (of 138733 total identified positions) were used to generate the SNVPhyl phylogeny shown in Figure 1. The number of single nucleotide differences between pediatric outbreak isolates is shown in Table 2. Isolates B, C, and D clustered closely in the phylogenetic tree, belonged to sequence type 100, and had 0-1 SNVs. A and E did not cluster with the outbreak isolates and belonged to sequence types 102 and 101, respectively. Four percent of the total identified positions were filtered out due to coverage.Figure 1
Phylogenetic analysis of isolates of suspected outbreak originating from a pediatric clinic. Tree was generated as output from the SNVPhyl pipeline run through IRIDA and visualized in MicroReact. The reference genome refers to ATCC19977. Default SNVPhyl parameters were used with a minimum coverage depth of 15 and mean mapping quality score of 30. Timeline shows date of sample isolation.Table 2
Matrix of single nucleotide differences between pediatric outbreak isolates.
Strain
A
B
C
D
E
Reference
A
0
108610
108610
108611
28650
29500
B
108610
0
0
1
107181
107907
C
108610
0
0
1
107181
107907
D
108611
1
1
0
107182
107908
E
28650
107181
107181
107182
0
19259
Reference
29500
107907
107907
107908
19259
0
Single nucleotide variations between study isolates as determined by SNVPhyl. Reference isM. abscessus ATCC19977.


## 3.2. Tattoo Outbreak

All tattoo isolates were identified asM. abscessus subsp. abscessus by MAB-MLST and PCR sequencing of 16S and hsp65. In addition to the four suspected outbreak cases (isolates H, I, J, K) thought to have originated from the tattoo parlour, unrelated outliers F and G were included in the analysis. Isolates H, I, J, and K clustered closely and belonged to sequence type 103. Isolates F and G did not cluster with the outbreak isolates and belonged to sequence type 1 and 83, respectively. The phylogenetic tree presented in Figure 2 was built from 46464 valid positions identified by SNVPhyl (of 50803 total identified positions). Eight percent of total identified positions were filtered out due to coverage or otherwise invalid. As seen in Table 3, the outbreak isolates have between 3 and 11 SNVs between each other in contrast to the outlier which has between 45933 and 45934 SNVs from the outbreak isolates (Table 3). Further, all suspected outbreak cases were assigned the sequence type (ST) as seen in Table S1 and S2. Drug susceptibility profile for isolate K was susceptible to clarithromycin with MIC of 2 μg/mL, while the related isolates were intermediary to clarithromycin with MIC of 4 μg/mL (Table S1).Figure 2
Phylogenetic tree of isolates belonging to the tattoo parlour outbreak. Tree generated as output from the SNVPhyl pipeline run through IRIDA and visualized in MicroReact. The reference genome refers to ATCC19977. Default parameters used with minimum coverage depth of 15 and minimum mean mapping quality score of 30. Timeline shows month sample isolation.Table 3
SNV matrix of tattoo isolates representing nucleotide differences between isolates.
Strain
H
I
J
K
G
F
Reference
H
0
3
4
7
45715
45909
45701
I
3
0
3
10
45714
45908
45700
J
4
3
0
11
45715
45909
45701
K
7
10
11
0
45714
45908
45700
G
45715
45714
45715
45714
0
1270
36
F
45909
45908
45909
45908
1270
0
1256
Reference
45701
45700
45701
45700
36
1256
0
Single nucleotide variations between study isolates as determined by SNVPhyl. Reference isM. abscessus ATCC19977.


## 4. Discussion

Laboratory methods used for infectious disease outbreak investigations have changed drastically over the years as WGS becomes more widespread and integrated as an epidemiological tool [21]. Methods of utilizing WGS data with a high discriminatory index that are fast, reliable, and easy to interpret are invaluable in public health systems. This study describes two separate small outbreaks caused by nontuberculous mycobacteria, Mycobacterium abscessus. One outbreak was detected at a pediatric clinic (n=4) and a second originated from a tattoo parlour (n=4). Both were identified by Public Health Laboratory in Quebec and sent to the National Microbiology Laboratory to be investigated through WGS. Two different methods, MAB-MLST and SNVPhyl, were used to investigate these outbreaks. MAB-MLST is a novel MLST-like scheme that authors previously developed, containing 14 genes involved in both housekeeping and antimicrobial resistance, and generates an output of a sequence type ranging from 1 to 116 [15]. SNVPhyl is an in-house developed IRIDA pipeline at the National Microbiology Laboratory that scans genomes for single nucleotide variants (SNVs), validates them, and uses these variants to build a phylogenetic tree [14].The first outbreak was suspected in a pediatric clinic in Montreal, Quebec, where four children aged 1-6 presented with symptoms of chronic ear infections between October 2017 and January 2018. Children are known to be more susceptible to infections by NTMs. Fomites, biofilms, water sources, and human to human transmission have all been suggested as source of infection [22]. M. abscessus has been found to inhabit municipal water supplies which can be a major issue in health-care facilities [2]. Due to their inherent resistance to many disinfectants, these organisms can survive on medical devices and spread infection and hence result in a health-care setting outbreak. For this outbreak, it was hypothesized that a medical device inappropriately decontaminated at a private clinic was the source of infection. An outlier acquired from another hospital (E) was included as a negative control. Of the isolates from the 4 submitted patients, phylogenetic tree and SNV matrix suggest that three isolates belong to the outbreak, and isolate A was acquired from an alternative source. This particular isolate was not only a different sequence type (ST) with different alleles at all 14 loci but belonged to a different subspecies. In contrast, the 3 isolates belonging to the outbreak clustered together and had 0-1 SNVs between them (Table 2; Figure 1). One T to G mutation at position 3378819 was identified between D and the other outbreak strains B and C; however, the comparison with M. abscessus ATCC19977 strain suggests that this is not an annotated gene. Due to the near 100% similarity of the three isolates as well as a single source linked to a medical device from one private clinic, it is concluded that all patients acquired the same strain from the same source at a private clinic.The second outbreak investigation stemmed from a tattoo parlour. Four patients presented with skin infections withM. abscessus. The time frame of this outbreak was 6 weeks. M. abscessus is a common cause of infections following cosmetic procedures [8]. The SNVs between outbreak strains H, I, J, and K ranged from 3 to 11. As per Bryant et al., this high degree of genetic similarity suggests that these infections originated from the same source [23]. Strains H, I, J, and K have same ST by MAB-MLST, but susceptibility testing data shows one MIC difference for clarithromycin, making K strain sensitive to the antibiotic while other 3 strains remain intermediate. This one MIC dilution difference is considered acceptable in a laboratory [16, 17]. As the clinical data highly suggests the same source for all strains and the number of SNVs between these strains are very low and meet cluster definition by Bryant et al., it can be safely concluded that these strains are related. Bryant et al. used a threshold of less than 25 SNVs to define a cluster; all outbreak strains fall well within this threshold [23]. Of the SNVs unique to the outbreak isolates in Table 3, they were located in 12 chromosomal positions, with the majority found in conserved hypothetical proteins. One of the SNVs found in H was located in a putative aminoglycoside phosphotransferase (MAB_0951). While M. abscessus is inherently resistant to aminoglycosides through the mechanism of aminoglycoside phosphotransferases, the presence of this SNV showed no effect on amikacin resistance [24]. Multiple SNVs were located different families of transcriptional regulators including the AraC (MAB_1337c), TetR (MAB_1337c and MAB_1873c), and the LysR family (MAB_0998). LysR is a diverse family of transcriptional regulators involved in many functions including virulence and metabolism [25]. In mycobacteria, a LysR type regulator, OxyS, was found to regulate the catalase-peroxidase gene and thus is a mechanism for adapting to oxidative stress by regulating katG [26]. This gene is involved in susceptibility to reactive oxygen species and isoniazid in Mycobacterium tuberculosis. However, there was no observable effect on isoniazid antimicrobial susceptibility data in the mutated M. abscessus strain I, as its wild type is intrinsically resistant to isoniazid as well (Table S1).A recent study utilized WGS to investigate an outbreak of tattoo-associated outbreak ofM. abscessus and was able to isolate the source as a bottle of grey wash ink [27]. In our study, the source of infection could not be identified in tattoo outbreak; however, it was hypothesized that tap water bottle was likely a cause since it was routinely used for tattoo ink dilution at the parlour. Drastic infection control measure was put into place including increased hygiene, otoscope disinfection and replacements, medical device disinfection protocol revision at the private clinic, disposal material, and reagent changes and disinfection. After a year, no further cases have occurred in either location. Laboratory analysis, time period, SNV, and outbreak origin suggest that these infections were acquired from the same source.Both WGS-based methods, MAB-MLST and SNVPhyl, showed concordance between the two outbreaks even with a limited number of study isolates. SNVPhyl data provided more discrimination in the tattoo outbreak. The isolate K had the same sequence type as the other tattoo outbreak isolates (H, I, and J) but had more SNVs identified by SNVPhyl. The number of SNVs between H, I, and J was 3-4 but between H, I, J, and K was 7-11 SNVs. As per Bryant et al., this difference is less than 25 SNVs, hence does not impact clinical outcome, and all 4 isolates belong to the same cluster [23].The major difference between MAB-MLST and SNVPhyl lies between the output of each respective pipeline. SNVPhyl provides more information on mutated genes, a phylogenetic tree, SNP matrix, or number of mutations present in the entire genome. The significance of each mutation in a specific gene is not easily determined. This data is valuable in research, not easily readable and requires bioinformatics knowledge to decipher this data. MAB-MLST provides direct input on STs as well as which gene profile is present. This gene profile is linked to antibiotic resistance data or to other identification markers since MAB-MLST has housekeeping genes, antimicrobial resistance markers, and markers of identification. This tool is easy to use; output is easy to understand and does not require bioinformatics expertise and hence can be utilized by general technicians in a clinical laboratory. MAB-MLST outputs tabular data with an allele number of each gene and associated ST. This provides an added advantage of directly outputting which genes are mutated and can have implications in easily predicting antimicrobial resistance. In addition, MAB-MLST also has the capability to differentiate the three subspecies by analyzing thehsp65, rpoB, and erm(41) genes of which are all utilized in subspecies identification [6, 28]. SNVPhyl has higher discriminatory power than MAB-MLST as it uses more genes. It generates multiple file types; the most prominent of which is a phylogenetic tree. While this data is useful in research, knowledge of phylogenetics for data interpretation can pose technical challenges for a diagnostic laboratory worker. SNVPhyl holds an advantage as it provides additional information on relatedness of outbreak strains. It produces a SNV matrix with the number of variants between each strain while MAB-MLST provides sequence type number and lacks SNV information. The applicability of either method is highly dependent on research or diagnostic question.


## 5. Conclusion

Infections caused by nontuberculous mycobacteria, particularlyM. abscessus, continue to be a growing health concern. As whole genome sequencing (WGS) becomes commonplace in diagnostic laboratories, its applicability in public health will continue to expand. Both tools utilized WGS data and were used to investigate two suspected outbreaks. MAB-MLST was able to identify outbreak isolates by assigning the same species and sequence type. However, SNVPhyl was able to go a step further and show that the pediatric outbreak isolates were more closely related than the tattoo isolates by showing single nucleotide differences not identified by MAB-MLST. Utilization of high throughput sequencing for outbreak investigation serves as an invaluable tool for routine mycobacteriology as it provides identification and antimicrobial susceptibility prediction as well as source case identification or genotyping by use of a single laboratory test. Both of these WGS-based methods are useful and can be easily adopted by different laboratories for NTM outbreak investigations.



---

*Source: 7092053-2020-11-09.xml*