Phylogenetic analysis of Mycobacterium massiliense strains having recombinant rpoB gene laterally transferred from Mycobacterium abscessus

Development of multiplex real-time PCR for detection of clarithromycin resistance genes for the Mycobacterium abscessus group

Introduction. The M. abscessus molecular identification and its drug-resistance profile are important to choose the correct therapy.Aim. This work developed a multiplex real-time PCR (mqPCR) for detection of clarithromycin resistance genes for the Mycobacterium abscessus group.Methodology. Isolates received by Adolfo Lutz Institute from 2010 to 2012, identified by PCR restriction enzyme analysis of a fragment of the hsp65 gene (PRA-hsp65) as M. abscessus type 1 (n=135) and 2 (n=71) were used. Drug susceptibility test (DST) for CLA were performed with reading on days 3 and 14. Subespecies identification by hsp65 and rpoB genes sequencing and erm(41) and rrl genes for mutation detection and primer design were performed. erm(41) gene deletion was detected by conventional PCR. Primers and probes were designed for five detections: erm(41) gene full size and with deletion; erm(41) gene T28 and C28; rrl gene A2058.Results. In total, 191/206 (92.7 %) isolates were concordant by all methods and 13/206 (6.3 %) were concordant only between molecular methods. Two isolates (1.0 %) were discordant by mqPCR compared to rrl gene sequencing. The mqPCR obtained 204/206 (99.0 %) isolates in agreement with the gold standard, with sensitivity and specificity of 98 and 100 %, respectively, considering the gold standard method and 92 and 93 % regarding DST.Conclusion. The mqPCR developed by us proved to be an easy-to-apply tool, minimizing time, errors and contamination.

Mycobacterium abscessus: insights from a bioinformatic perspective

Mycobacterium abscessus is a nontuberculous mycobacterium, associated with broncho-pulmonary infections in individuals suffering from cystic fibrosis, bronchiectasis, and pulmonary diseases. The risk factors for transmission include biofilms, contaminated water resources, fomites, and infected individuals. M. abscessus is extensively resistant to antibiotics. To date, there is no vaccine and combination antibiotic therapy is followed. However, drug toxicities, low cure rates, and high cost of treatment make it imperfect. Over the last 20 years, bioinformatic studies on M. abscessus have advanced our understanding of the pathogen. This review integrates knowledge from the analysis of genomes, microbiomes, genomic variations, phylogeny, proteome, transcriptome, secretome, antibiotic resistance, and vaccine design to further our understanding. The utility of genome-based studies in comprehending disease progression, surveillance, tracing transmission routes, and epidemiological outbreaks on a global scale has been highlighted. Furthermore, this review underlined the importance of using computational methodologies for pinpointing factors responsible for pathogen survival and resistance. We reiterate the significance of interdisciplinary research to fight M. abscessus. In a nutshell, the outcome of computational studies can go a long way in creating novel therapeutic avenues to control M. abscessus mediated pulmonary infections.

Clinical Whole Genome Sequencing for Clarithromycin and Amikacin Resistance Prediction and Subspecies Identification of Mycobacterium abscessus

Mycobacterium abscessus infections are an emerging health care concern in patients with chronic pulmonary diseases, leading to high morbidity and mortality. One major challenge is resistance to clarithromycin, a cornerstone antibiotic with high efficacy. Therefore, treatment is primarily guided by phenotypic susceptibility results of clarithromycin, which requires extended incubation to assess for inducible resistance. Resistance mechanisms for clarithromycin include induction of erm(41) and mutations in the 23S rRNA gene (rrl). In addition, mutations in the 16S rRNA encoding gene (rrs) can confer high-level amikacin resistance, another essential drug in the treatment of M. abscessus infections. Herein, we developed a clinical whole genome sequencing (WGS) assay for clarithromycin resistance based on rrl and erm(41) gene sequences and amikacin resistance based on the rrs sequence in M. abscessus, as well as subspecies identification. Genotypic-based predictions were determined for 104 isolates from 68 patients. The overall accuracy of genotypic prediction for clarithromycin compared with phenotypic susceptibility results was 100% (95% CI, 96.45%-100%). For amikacin, we also obtained 100% accuracy (95% CI, 96.52%-100%). The high concordance between the genotypic and phenotypic results demonstrates that a WGS-based assay can be used in a clinical laboratory for determining resistance to clarithromycin and amikacin in M. abscessus isolates. WGS can also provide subspecies identification and high-definition phylogenetic information for more accurate M. abscessus strain typing.

A novel DNA chromatography method to discriminate Mycobacterium abscessus subspecies and macrolide susceptibility

Background

The clinical impact of infection with Mycobacterium (M.) abscessus complex (MABC), a group of emerging non-tuberculosis mycobacteria (NTM), is increasing. M. abscessus subsp. abscessus/bolletii frequently shows natural resistance to macrolide antibiotics, whereas M. abscessus subsp. massiliense is generally susceptible. Therefore, rapid and accurate discrimination of macrolide-susceptible MABC subgroups is required for effective clinical decisions about macrolide treatments for MABC infection. We aimed to develop a simple and rapid diagnostic that can identify MABC isolates showing macrolide susceptibility.

Methods

Whole genome sequencing (WGS) was performed for 148 clinical or environmental MABC isolates from Japan to identify genetic markers that can discriminate three MABC subspecies and the macrolide-susceptible erm(41) T28C sequevar. Using the identified genetic markers, we established PCR based- or DNA chromatography-based assays. Validation testing was performed using MABC isolates from Taiwan.

Finding

We identified unique sequence regions that could be used to differentiate the three subspecies. Our WGS-based phylogenetic analysis indicated that M. abscessus carrying the macrolide-susceptible erm(41) T28C sequevar were tightly clustered, and identified 11 genes that were significantly associated with the lineage for use as genetic markers. To detect these genetic markers and the erm(41) locus, we developed a DNA chromatography method that identified three subspecies, the erm(41) T28C sequevar and intact erm(41) for MABC in a single assay within one hour. The agreement rate between the DNA chromatography-based and WGS-based identification was 99·7%.

Interpretation

We developed a novel, rapid and simple DNA chromatography method for identification of MABC macrolide susceptibility with high accuracy.

Funding

AMED, JSPS KAKENHI

A multilocus sequence typing scheme for Mycobacterium abscessus complex (MAB-multilocus sequence typing) using whole-genome sequencing data

Background

Mycobacterium abscessus is a rapid growing nontuberculous mycobacteria (NTM) and a clinically significant pathogen capable of causing variable infections in humans that are difficult to treat and may require months of therapy/surgical interventions. Like other NTMs, M. abscessus can be associated with outbreaks leading to complex investigations and treatment of affected cases. Typing schemes for bacterial pathogens provide numerous applications; including identifying chain of transmission and tracking genomic evolution, are lacking or limited for many NTMs including M. abscessus.

Methods

We extended the existing scheme from PubMLST using whole-genome data for M. abscessus by extracting data for 15 genetic regions within the M. abscessus genome. A total of 168 whole genomes and 11 gene sequences were used to build this scheme (MAB-multilocus sequence typing [MLST]).

Results

All seven genes from the PubMLST scheme, namely argH, cya, gnd, murC, pta, purH, and rpoB, were expanded by 10, 14, 13, 10, 13, 10, and 9 alleles, respectively. Another eight novel genes were added including hsp 65, erm(41), arr, rrs, rrl, gyrA, gyrB, and recA with 16, 16, 25, 7, 32, 35, 29, and 15 alleles, respectively, with 85 unique sequence types identified among all isolates.

Conclusion

MAB-MLST can provide identification of M. abscessus complex to the subspecies level based on three genes and can provide antimicrobial resistance susceptibility prediction based on results from seven genes. MAB-MLST generated a total of 85 STs, resulting in subtyping of 90 additional isolates that could not be genotyped using PubMLST and yielding results comparable to whole-genome sequencing (WGS). Implementation of a Galaxy-based data analysis tool, MAB-MLST, that simplifies the WGS data and yet maintains a high discriminatory index that can aid in deciphering an outbreak has vast applicability for routine diagnostics.

PCR amplification of the erm(41) gene can be used to predict the sensitivity of Mycobacterium abscessus complex strains to clarithromycin

A worldwide increase in the Mycobacterium abscessus (M. abscessus) complex has been observed. Therefore, the aim of the present study was to investigate the diversity of the rrl and erm(41) genes, both of which are associated with macrolide sensitivity in the M. abscessus complex. The current study also examined the efficacy of mass spectrometry as an alternative to molecular testing to classify subspecies of the M. abscessus complex. A total of 14 strains of the M. abscessus complex were obtained, and based on conventional analyses using housekeeping genes, 57% were determined to be M. abscessus subsp. abscessus, 43% were M. abscessus subsp. massiliense, and none were identified as M. abscessus subsp. bolletii. However, depending on the strain, it was not always possible to distinguish between the subspecies by mass spectrometry. Consequently, PCR products for the rrl and erm(41) genes were directly sequenced. Overall, 7.1% of the strains were identified to have a rrl mutation, and 92.9% carried a T at position 28 of erm(41). Results presented here suggest that the principal cause of treatment failure for M. abscessus complex infections is inducible macrolide resistance encoded by the erm(41) gene. From a strictly pragmatic standpoint, the phenotypic function of a putative erm(41) gene is the most important piece of information required by clinicians in order to prescribe an effective treatment. Although PCR amplification of erm(41) is not sufficient to differentiate between the M. abscessus complex subspecies, PCR can be easily and efficiently used to predict the sensitivity of members of the M. abscessus complex to clarithromycin.

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

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.

Multilocus sequence typing detects new Piscirickettsia salmonis hybrid genogroup in Chilean fish farms: Evidence for genetic diversity and population structure

Piscirickettsia salmonisis the causative bacterial pathogen of piscirickettsiosis, a salmonid disease that causes notable mortalities in the worldwide aquaculture industry. Published research describes the phenotypic traits, virulence factors, pathogenicity and antibiotic-resistance potential for various P. salmonisstrains. However, evolutionary and genetic information is scarce for P. salmonis. The present study used multilocus sequence typing (MLST) to gain insight into the population structure and evolution of P. salmonis. Forty-two Chilean P. salmonisisolates, as well as the type strain LF-89^T , were recovered from diseased Salmo salar, Oncorhynchus kisutchand Oncorhynchus mykissfrom two Chilean Regions. MLST assessed the loci sequences of dnaK, efp, fumC, glyA, murG, rpoD and trpB. Bioinformatics analyses established the genetic diversity among P. salmonis isolates (H = 0.5810). A total of 23 sequence types (ST) were identified, 53.48% of which were represented by ST1, ST5 and ST2. Population structure analysis through polymorphism patterns showed few polymorphic sites (218 nucleotides from 4,010 bp), while dN/dS ratio analysis indicated purifying selection for dnaK, epf, fumC, murG, and rpoD but neutral selection for the trpB loci. The standardized index of association indicated strong linkage disequilibrium, suggesting clonal population structure. However, recombination events were detected in a group of seven isolates. Findings included genogroups homologous to the LF-89^T and EM-90 strains, as well as a seven-isolate hybrid genogroup recovered from both assessed regions (three O. mykiss and four S. salar isolates). The presented MLST scheme has comparative potential, with promising applications in studying distinct P. salmonis isolates (e.g., from different hosts, farms, geographical areas) and in understanding the epidemiology of this pathogen.

Role of the DNA Mismatch Repair Gene MutS4 in Driving the Evolution of Mycobacterium yongonense Type I via Homologous Recombination

We recently showed that Mycobacterium yongonense could be divided into two genotypes: Type I, in which the rpoB gene has been transferred from Mycobacterium parascrofulaceum, and Type II, in which the rpoB gene has not been transferred. Comparative genome analysis of three M. yongonense Type I, two M. yongonense Type II and M. parascrofulaceum type strains were performed in this study to gain insight into gene transfer from M. parascrofulaceum into M. yongonense Type I strains. We found two genome regions transferred from M. parascrofulaceum: one contained 3 consecutive genes, including the rpoBC operon, and the other contained 57 consecutive genes that had been transferred into M. yongonense Type I genomes via homologous recombination. Further comparison between the M. yongonense Type I and II genomes revealed that Type I, but not Type II has a distinct DNA mismatch repair gene (MutS4 subfamily) that was possibly transferred via non-homologous recombination from other actinomycetes. We hypothesized that it could facilitate homologous recombination from the M. parascrofulaceum to the M. yongonense Type I genomes. We therefore generated recombinant Mycobacterium smegmatis containing a MutS4 operon of M. yongonense. We found that the M. tuberculosis rpoB fragment with a rifampin resistance-conferring mutation was more frequently inserted into recombinant M. smegmatis than the wild type, suggesting that MutS4 is a driving force in the gene transfer from M. parascrofulaceum to M. yongonense Type I strains via homologous recombination. In conclusion, our data indicated that MutS4 in M. yongonense Type I genomes may drive gene transfer from M. parascrofulaceum via homologous recombination, resulting in division of M. yongonense into two genotypes, Type I and II.