A phylogenomic approach to bacterial subspecies classification: proof of concept in Mycobacterium abscessus

Identification of Mycobacterium abscessus using the peaks of ribosomal protein L29, L30 and hemophore-related protein by MALDI-MS proteotyping

Mycobacteroides (Mycobacterium) abscessus, which causes a variety of infectious diseases in humans, is becoming detected more frequently in clinical specimens as cases are spreading worldwide. Taxonomically, M. abscessus is composed of three subspecies of M. abscessus subsp. abscessus, M. abscessus subsp. bolletii, and M. abscessus subsp. massiliense, with different susceptibilities to macrolides. In order to identify rapidly these three subspecies, we determined useful biomarker proteins, including ribosomal protein L29, L30, and hemophore-related protein, for distinguishing the subspecies of M. abscessus using the matrix-assisted laser desorption/ionization mass spectrometry (MALDI-MS) profiles. Thirty-three clinical strains of M. abscessus were correctly identified at the subspecies-level by the three biomarker protein peaks. This study ultimately demonstrates the potential of routine MALDI-MS-based laboratory methods for early identification and treatment for M. abscessus infections.

Mycolicibacterium cyprinidarum sp. nov., a rapidly growing species isolated from diseased koi carp, Cyprinus carpio

A rapidly growing nontuberculous mycobacterium was isolated from diseased koi carp in Niigata, Japan, which was identified as representing a novel Mycolicibacterium species through whole genome sequence analysis. The bacterial isolates (NGTWS0302, NGTWS1803T and NGTWSNA01) were found to belong to the genus Mycolicibacterium through phylogenetic analysis using whole genome sequences of mycobacteria species. The bacterial colony was smooth, moist and non-chromogenic on 1% Ogawa medium at 30 °C. In biochemical characteristic tests, the bacterial isolates showed positive reactions for catalase activity, Tween 80 hydrolysis and tellurite reduction. The isolates were sensitive to 2-4 µg ml-1 ampicillin, kanamycin and rifampicin. Based on these results, we propose a novel Mycolicibacterium species, Mycolicibacterium cyprinidarum sp. nov. The type strain is NGTWS1803T (=JCM 35117T=ATCC TSD-289T).

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.

The Use of Comparative Genomic Analysis for the Development of Subspecies-Specific PCR Assays for Mycobacterium abscessus

Mycobacterium abscessus complex (MABC) is an important pathogen of immunocompromised patients. Accurate and rapid determination of MABC at the subspecies level is vital for optimal antibiotic therapy. Here we have used comparative genomics to design MABC subspecies-specific PCR assays. Analysis of single nucleotide polymorphisms and core genome multilocus sequence typing showed clustering of genomes into three distinct clusters representing the MABC subspecies M. abscessus, M. bolletii and M. massiliense. Pangenome analysis of 318 MABC genomes from the three subspecies allowed for the identification of 15 MABC subspecies-specific genes. In silico testing of primer sets against 1,663 publicly available MABC genomes and 66 other closely related Mycobacterium genomes showed that all assays had >97% sensitivity and >98% specificity. Subsequent experimental validation of two subspecies-specific genes each showed the PCR assays worked well in individual and multiplex format with no false-positivity with 5 other mycobacteria of clinical importance. In conclusion, we have developed a rapid, accurate, multiplex PCR-assay for discriminating MABC subspecies that could improve their detection, diagnosis and inform correct treatment choice.

β-Lactam Combinations That Exhibit Synergy against Mycobacteroides abscessus Clinical Isolates

Mycobacteroides abscessus (Mab) is an opportunistic environmental pathogen that can cause chronic pulmonary disease in the setting of structural lung conditions such as bronchiectasis, chronic obstructive pulmonary disease, and cystic fibrosis. These infections are often incurable and associated with rapid lung function decline. Mab is naturally resistant to most of the antibiotics available today, and current treatment guidelines require at least 1 year of daily multidrug therapy, which is often ineffective and is associated with significant toxicities. β-Lactams are the most widely used class of antibiotics and have a demonstrated record of safety and tolerability. Here, using a panel of recent clinical isolates of Mab, we evaluated the in vitro activities of dual-β-lactam combinations to identify new treatments with the potential to treat infections arising from a wide range of Mab strains. The Mab clinical isolates were heterogeneous, as reflected by the diversity of their genomes and differences in their susceptibilities to various drugs. Cefoxitin and imipenem are currently the only two β-lactams included in the guidelines for treating Mab disease, yet they are not used concurrently in clinical practice. However, this dual-β-lactam combination exhibited synergy against 100% of the isolates examined (n = 21). Equally surprising is the finding that the combination of two carbapenems, doripenem and imipenem, exhibited synergy against the majority of Mab isolates. In the setting of multidrug-resistant Mab disease with few therapeutic options, these combinations may offer viable immediate treatment options with efficacy against the broad spectrum of Mab strains infecting patients today.

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.

The recombination-cold region as an epidemiological marker of recombinogenic opportunistic pathogen Mycobacterium avium

Background

The rapid identification of lineage remains a challenge in the genotyping of clinical isolates of recombinogenic pathogens. The chromosome of Mycobacterium avium subsp. hominissuis (MAH), an agent of Mycobacterium avium complex (MAC) lung disease, is often mosaic and is composed of chromosomal segments originating from different lineages. This makes it difficult to infer the MAH lineage in a simple experimental set-up. To overcome this difficulty, we sought to identify chromosomal marker genes containing lineage-specific alleles by genome data mining.

Results

We conducted genetic population structure analysis, phylogenetic analysis, and a survey of historical recombination using data from 125 global MAH isolates. Six MAH lineages (EA1, EA2, SC1, SC2, SC3, and SC4) were identified in the current dataset. One P-450 gene (locus_tag MAH_0788/MAV_0940) in the recombination-cold region was found to have multiple alleles that could discriminate five lineages. By combining the information about allele type from one additional gene, the six MAH lineages as well as other M. avium subspecies were distinguishable. A recombination-cold region of 116 kb contains an insertion hotspot and is flanked by a mammalian cell-entry protein operon where allelic variants have previously been reported to occur. Hence, we speculate that the acquisition of lineage- or strain-specific insertions has introduced homology breaks in the chromosome, thereby reducing the chance of interlineage recombination.

Conclusions

The allele types of the newly identified marker genes can be used to predict major lineages of M. avium. The single nucleotide polymorphism typing approach targeting multiallelic loci in recombination-cold regions will facilitate the epidemiological study of MAC, and may also be useful for equivalent studies of other nontuberculous mycobacteria potentially carrying mosaic genomes.

Improved genome of Agrobacterium radiobacter type strain provides new taxonomic insight into Agrobacterium genomospecies 4

The reported Agrobacterium radiobacter DSM 30174^T genome is highly fragmented, hindering robust comparative genomics and genome-based taxonomic analysis. We re-sequenced the Agrobacterium radiobacter type strain, generating a dramatically improved genome with high contiguity. In addition, we sequenced the genome of Agrobacterium tumefaciens B6^T, enabling for the first time, a proper comparative genomics of these contentious Agrobacterium species. We provide concrete evidence that the previously reported Agrobacterium radiobacter type strain genome (Accession Number: ASXY01) is contaminated which explains its abnormally large genome size and fragmented assembly. We propose that Agrobacterium tumefaciens be reclassified as Agrobacterium radiobacter subsp. tumefaciens and that Agrobacterium radiobacter retains it species status with the proposed name of Agrobacterium radiobacter subsp. radiobacter. This proposal is based, first on the high pairwise genome-scale average nucleotide identity supporting the amalgamation of both Agrobacterium radiobacter and Agrobacterium tumefaciens into a single species. Second, maximum likelihood tree construction based on the concatenated alignment of shared genes (core genes) among related strains indicates that Agrobacterium radiobacter NCPPB3001 is sufficiently divergent from Agrobacterium tumefaciens to propose two independent sub-clades. Third, Agrobacterium tumefaciens demonstrates the genomic potential to synthesize the L configuration of fucose in its lipid polysaccharide, fostering its ability to colonize plant cells more effectively than Agrobacterium radiobacter.

Comparative evaluation of the identification of rapidly growing non-tuberculous mycobacteria by mass spectrometry (MALDI-TOF MS), GenoType Mycobacterium CM/AS assay and partial sequencing of the rpoβ gene with phylogenetic analysis as a reference method

INTRODUCTION

The American Thoracic Society and the Infectious Diseases Society of America recommend that clinically significant non-tuberculous mycobacteria (NTM) should be identified to the species level in order to determine their clinical significance. The aim of this study was to evaluate identification of rapidly growing NTM (RGM) isolated from clinical samples by using MALDI-TOF MS and a commercial molecular system. The results were compared with identification using a reference method.

METHODS

We included 46 clinical isolates of RGM and identified them using the commercial molecular system GenoType^® CM/AS (Hain, Lifescience, Germany), MALDI-TOF MS (Bruker) and, as reference method, partial rpoβ gene sequencing followed by BLAST and phylogenetic analysis with the 1093 sequences available in the GeneBank.

RESULTS

The degree of agreement between GenoType^® and MALDI-TOF MS and the reference method, partial rpoβ sequencing, was 27/43 (62.8%) and 38/43 cases (88.3%) respectively. For all the samples correctly classified by GenoType^®, we obtained the same result with MALDI-TOF MS (27/27). However, MALDI-TOF MS also correctly identified 68.75% (11/16) of the samples that GenoType^® had misclassified (p=0.005).

CONCLUSIONS

MALDI-TOF MS classified significantly better than GenoType^®. When a MALDI-TOF MS score >1.85 was achieved, MALDI-TOF MS and partial rpoβ gene sequencing were equivalent. GenoType^® was not able to distinguish between species belonging to the M. fortuitum complex. MALDI-TOF MS methodology is simple, rapid and associated with lower consumable costs than GenoType^®. The partial rpoβ sequencing methods with BLAST and phylogenetic analysis were not able to identify some RGM unequivocally. Therefore, sequencing of additional regions would be indicated in these cases.

Genomic Comparisons Reveal Microevolutionary Differences in Mycobacterium abscessus Subspecies

Mycobacterium abscessus, a rapid-growing non-tuberculous mycobacterium, has been the cause of sporadic and outbreak infections world-wide. The subspecies in M. abscessus complex (M. abscessus, M. massiliense, and M. bolletii) are associated with different biologic and pathogenic characteristics and are known to be among the most frequently isolated opportunistic pathogens from clinical material. To date, the evolutionary forces that could have contributed to these biological and clinical differences are still unclear. We compared genome data from 243 M. abscessus strains downloaded from the NCBI ftp Refseq database to understand how the microevolutionary processes of homologous recombination and positive selection influenced the diversification of the M. abscessus complex at the subspecies level. The three subspecies are clearly separated in the Minimum Spanning Tree. Their MUMi-based genomic distances support the separation of M. massiliense and M. bolletii into two subspecies. Maximum Likelihood analysis through dN/dS (the ratio of number of non-synonymous substitutions per non-synonymous site, to the number of synonymous substitutions per synonymous site) identified distinct genes in each subspecies that could have been affected by positive selection during evolution. The results of genome-wide alignment based on concatenated locally-collinear blocks suggest that (a) recombination has affected the M. abscessus complex more than mutation and positive selection; (b) recombination occurred more frequently in M. massiliense than in the other two subspecies; and (c) the recombined segments in the three subspecies have come from different intra-species and inter-species origins. The results lead to the identification of possible gene sets that could have been responsible for the subspecies-specific features and suggest independent evolution among the three subspecies, with recombination playing a more significant role than positive selection in the diversification among members in this complex.

Whole Genome Sequencing Allows Better Understanding of the Evolutionary History of Leptospira interrogans Serovar Hardjo

The genome of a laboratory-adapted strain of Leptospira interrogans serovar Hardjo was sequenced and analyzed. Comparison of the sequenced genome with that recently published for a field isolate of the same serovar revealed relatively high sequence conservation at the nucleotide level, despite the different biological background of both samples. Conversely, comparison of both serovar Hardjo genomes with those of L. borgpetersenii serovar Hardjo showed extensive differences between the corresponding chromosomes, except for the region occupied by their rfb loci. Additionally, comparison of the serovar Hardjo genomes with those of different L. interrogans serovars allowed us to detect several genomic features that may confer an adaptive advantage to L. interrogans serovar Hardjo, including a possible integrated plasmid and an additional copy of a cluster encoding a membrane transport system known to be involved in drug resistance. A phylogenomic strategy was used to better understand the evolutionary position of the Hardjo serovar among L. interrogans serovars and other Leptospira species. The proposed phylogeny supports the hypothesis that the presence of similar rfb loci in two different species may be the result of a lateral gene transfer event.

Genome-wide mosaicism within Mycobacterium abscessus: evolutionary and epidemiological implications

Background

In mycobacteria, conjugation differs from the canonical Hfr model, but is still poorly understood. Here, we quantified this evolutionary processe in a natural mycobacterial population, taking advantage of a large clinical strain collection of the emerging pathogen Mycobacterium abscessus (MAB).

Results

Multilocus sequence typing confirmed the existence of three M. abscessus subspecies, and unravelled extensive allelic exchange between them. Furthermore, an asymmetrical gene flow occurring between these main lineages was detected, resulting in highly admixed strains. Intriguingly, these mosaic strains were significantly associated with cystic fibrosis patients with lung infections or chronic colonization. Genome sequencing of those hybrid strains confirmed that half of their genomic content was remodelled in large genomic blocks, leading to original tri-modal ‘patchwork’ architecture. One of these hybrid strains acquired a locus conferring inducible macrolide resistance, and a large genomic insertion from a slowly growing pathogenic mycobacteria, suggesting an adaptive gene transfer. This atypical genomic architecture of the highly recombinogenic strains is consistent with the distributive conjugal transfer (DCT) observed in M. smegmatis. Intriguingly, no known DCT function was found in M. abscessus chromosome, however, a p-RAW-like genetic element was detected in one of the highly admixed strains.

Conclusion

Taken together, our results strongly suggest that MAB evolution is sporadically punctuated by dramatic genome wide remodelling events. These findings might have far reaching epidemiological consequences for emerging mycobacterial pathogens survey in the context of increasing numbers of rapidly growing mycobacteria and M. tuberculosis co-infections.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2448-1) contains supplementary material, which is available to authorized users.

MycoCAP - Mycobacterium Comparative Analysis Platform

Mycobacterium spp. are renowned for being the causative agent of diseases like leprosy, Buruli ulcer and tuberculosis in human beings. With more and more mycobacterial genomes being sequenced, any knowledge generated from comparative genomic analysis would provide better insights into the biology, evolution, phylogeny and pathogenicity of this genus, thus helping in better management of diseases caused by Mycobacterium spp.With this motivation, we constructed MycoCAP, a new comparative analysis platform dedicated to the important genus Mycobacterium. This platform currently provides information of 2108 genome sequences of at least 55 Mycobacterium spp. A number of intuitive web-based tools have been integrated in MycoCAP particularly for comparative analysis including the PGC tool for comparison between two genomes, PathoProT for comparing the virulence genes among the Mycobacterium strains and the SuperClassification tool for the phylogenic classification of the Mycobacterium strains and a specialized classification system for strains of Mycobacterium abscessus. We hope the broad range of functions and easy-to-use tools provided in MycoCAP makes it an invaluable analysis platform to speed up the research discovery on mycobacteria for researchers. Database URL: http://mycobacterium.um.edu.my

Inducible and Acquired Clarithromycin Resistance in the Mycobacterium abscessus Complex

PURPOSE

Clarithromycin was considered the cornerstone for the treatment of Mycobacterium abscessus complex infections. Genetic resistance mechanisms have been described and many experts propose amikacin as an alternative. Nevertheless, clarithromycin has several advantages; therefore, it is necessary to identify the non-functional erm(41) allele to determine the most suitable treatment. The aims of this study were to characterize the molecular mechanisms of clarithromycin resistance in a collection of Mycobacterium abscessus complex isolates and to verify the relationship between these mechanisms and the antibiogram.

MATERIALS AND METHODS

Clinical isolates of M. abscessus complex (n = 22) from 16 patients were identified using four housekeeping genes (rpoB, secA1, sodA and hsp65), and their genetic resistance was characterized by studying erm(41) and rrl genes. Nine strains were recovered from the clinical isolates and subjected to E-test and microdilution clarithromycin susceptibility tests, with readings at 3, 7 and 14 days.

RESULTS

We classified 11/16 (68.8%) M. abscessus subsp. abscessus, 4/16 (25.0%) M. abscessus subsp. bolletii, and 1/16 (6.3%) M. abscessus subsp. massiliense. T28 erm(41) allele was observed in 8 Mycobacterium abscessus subps. abscessus and 3 Mycobacterium abscessus subsp. bolletii. One strain of M. abscessus subsp. bolletii had an erm(41) gene truncated and was susceptible to clarithromycin. No mutations were observed in rrl gene first isolates. In three patients, follow-up of initial rrl wild-type strains showed acquired resistance.

CONCLUSIONS

Most clinical isolates of M. abscessus complex had inducible resistance to clarithromycin and total absence of constitutive resistance. Our findings showed that the acquisition of resistance mutations in rrl gene was associated with functional and non-functional erm(41) gene. Caution is needed when using erm(41) sequencing alone to identify M. abscessus subspecies. This study reports an acquired mutation at position 2057 of rrl gene, conferring medium-low clarithromycin constitutive resistance.

Support from Phylogenomic Networks and Subspecies Signatures for Separation of Mycobacterium massiliense from Mycobacterium bolletii

Mycobacterium abscessus subspecies classification has important clinical implications. We used phylogenomic network and amino acid analyses to provide evidence for the separation of Mycobacterium bolletii and Mycobacterium massiliense into two distinct subspecies which can potentially be differentiated rapidly by their protein signatures.

Isolation and characterization of a native avirulent strain of Streptococcus suis serotype 2: a perspective for vaccine development

Streptococcus suis, an emerging infectious pathogen, is the cause of two large-scale outbreaks of human streptococcal toxic shock syndrome in China, and has attracted much attention from the scientific community. The genetic basis of its pathogenesis remains enigmatic, and no effective prevention measures have been established. To better understand the virulence differentiation of S. suis and develop a promising vaccine, we isolated and sequenced a native avirulent S. suis strain (05HAS68). Animal experiments revealed that 05HAS68 is an avirulent strain and could protect piglets from the attack of virulent strains. Comparative genomics analyses demonstrated the genetic basis for the lack of virulence in 05HAS68, which is characterized by the absence of some important virulence-associated factors and the intact 89K pathogenicity island. Lack of virulence was also illustrated by reduced survival of 05HAS68 compared to a virulent strain in pig whole blood. Further investigations revealed a large-scale genomic rearrangement in 05HAS68, which was proposed to be mediated by transposase genes and/or prophages. This genomic rearrangement may have caused the genomic diversity of S. suis, and resulted in biological discrepancies between 05HAS68 and highly virulent S. suis strains.