The resistomes of Mycobacteroides abscessus complex and their possible acquisition from horizontal gene transfer.

BACKGROUND: Mycobacteroides abscessus complex (MABC), an emerging pathogen, causes human infections resistant to multiple antibiotics. In this study, the genome data of 1,581 MABC strains were downloaded from NCBI database for phylogenetic relatedness inference, resistance profile identification and the estimation of evolutionary pressure on resistance genes in silico. RESULTS: From genes associated with resistance to 28 antibiotic classes, 395 putative proteins (ARPs) were identified, based on the information in two antibiotic resistance databases (CARD and ARG-ANNOT). The ARPs most frequently identified in MABC were those associated with resistance to multiple antibiotic classes, beta-lactams and aminoglycosides. After excluding ARPs that had undergone recombination, two ARPs were predicted to be under diversifying selection and 202 under purifying selection. This wide occurrence of purifying selection suggested that the diversity of commonly shared ARPs in MABC have been reduced to achieve stability. The unequal distribution of ARPs in members of the MABC could be due to horizontal gene transfer or ARPs pseudogenization events. Most (81.5%) of the ARPs were observed in the accessory genome and 72.2% ARPs were highly homologous to proteins associated with mobile genetic elements such as plasmids, prophages and viruses. On the other hand, with TBLASTN search, only 18 of the ARPs were identified as pseudogenes. CONCLUSION: Altogether, our results suggested an important role of horizontal gene transfer in shaping the resistome of MABC.

Genomic diversity of SARS-CoV-2 in Malaysia.

BACKGROUND: More than a year after its first appearance in December 2019, the COVID-19 pandemic is still on a rampage in many parts of the world. Although several vaccines have been approved for emergency use, the emergence and rapid spread of new SARS-CoV-2 variants have sparked fears of vaccine failure due to immune evasion. Massive viral genome sequencing has been recommended to track the genetic changes that could lead to adverse consequences. METHODS: We sequenced SARS-CoV-2 respiratory isolates from the National Public Health Laboratory, Malaysia and examined them together with viral genomes deposited in GISAID by other Malaysian researchers, to understand the evolutionary trend of the virus circulating in the country. We studied the distribution of virus lineages and site-wise mutations, analysed genetic clustering with the goeBURST full Minimum Spanning Tree algorithm, examined the trend of viral nucleotide diversity over time and performed nucleotide substitution association analyses. RESULTS: We identified 22 sub-lineages, 13 clonal complexes, 178 sequence types and seven sites of linkage disequilibrium in 277 SARS-CoV-2 genomes sequenced between January and December 2020. B.1.524 was the largest lineage group. The number of mutations per genome ranged from 0 to 19. The mean genomic diversity value over 12 months was 3.26 × 10-4. Of 359 mutations detected, 60.5% of which were non-synonymous, the most frequent were in the ORF1ab (P4715L), S (D614G and A701V) and N (S194L) genes. CONCLUSION: The SARS-CoV-2 virus accumulated an abundance of mutations in the first year of the COVID-19 pandemic in Malaysia. Its overall genetic diversity, however, is relatively low compared to other Asian countries with larger populations. Continuous genomic and epidemiological surveillance will help to clarify the evolutionary processes determining viral diversity and impacting on human health.

Genome sequence analysis of multidrug-resistant Mycobacterium tuberculosis from Malaysia.

Mycobacterium tuberculosis (MTB) is commonly used as a model to study pathogenicity and multiple drug resistance in bacteria. These MTB characteristics are highly dependent on the evolution and phylogeography of the bacterium. In this paper, we describe 15 new genomes of multidrug-resistant MTB (MDRTB) from Malaysia. The assessments and annotations on the genome assemblies suggest that strain differences are due to lineages and horizontal gene transfer during the course of evolution. The genomes show mutations listed in current drug resistance databases and global MTB collections. This genome data will augment existing information available for comparative genomic studies to understand MTB drug resistance mechanisms and evolution.

Tigecycline resistance may be associated with dysregulated response to stress in Mycobacterium abscessus.

Previously, we characterized 7C, a laboratory-derived tigecycline-resistant mutant of Mycobacterium abscessus ATCC 19977, and found that the resistance was conferred by a mutation in MAB_3542c, which encodes an RshA-like protein. In M. tuberculosis, RshA is an anti-sigma factor that negatively regulates the SigH-dependent heat/oxidative stress response. We hypothesized that this mutation in 7C might dysregulate the stress response which has been generally linked to antibiotic resistance. In this study, we tested this hypothesis by subjecting 7C to transcriptomic dissection using RNA sequencing. We found an over-expression of genes encoding the SigH ortholog, chaperones and oxidoreductases. In line with these findings, 7C demonstrated better survival against heat shock when compared to the wild-type ATCC 19977. Another interesting observation from the RNA-Seq analysis was the down-regulation of ribosomal protein-encoding genes. This highlights the possibility of ribosomal conformation changes which could negatively affect the binding of tigecycline to its target, leading to phenotypic resistance. We also demonstrated that transient resistance to tigecycline could be induced in the ATCC 19977 by elevated temperature. Taken together, these findings suggest that dysregulated stress response may be associated with tigecycline resistance in M. abscessus.

A mutation in anti-sigma factor MAB_3542c may be responsible for tigecycline resistance in Mycobacterium abscessus.

In this study, we characterized 7C, a spontaneous mutant selected from tigecycline-susceptible Mycobacterium abscessus ATCC 19977. Whole-genome sequencing (WGS) was used to identify possible resistance determinants in this mutant. Compared to the wild-type, 7C demonstrated resistance to tigecycline as well as cross-resistance to imipenem, and had a slightly retarded growth rate. WGS and subsequent biological verifications showed that these phenotypes were caused by a point mutation in MAB_3542c, which encodes an RshA-like protein. In Mycobacterium tuberculosis, RshA is an anti-sigma factor that negatively regulates the heat/oxidative stress response mechanisms. The MAB_3542c mutation may represent a novel determinant of tigecycline resistance. We hypothesize that this mutation may dysregulate the stress-response pathways which have been shown to be linked to antibiotic resistance in previous studies.

Chromosomal rearrangements and protein globularity changes in Mycobacterium tuberculosis isolates from cerebrospinal fluid.

BACKGROUND: Meningitis is a major cause of mortality in tuberculosis (TB). It is not clear what factors promote central nervous system invasion and pathology but it has been reported that certain strains of Mycobacterium tuberculosis (Mtb) might have genetic traits associated with neurotropism. METHODS: In this study, we generated whole genome sequences of eight clinical strains of Mtb that were isolated from the cerebrospinal fluid (CSF) of patients presenting with tuberculous meningitis (TBM) in Malaysia, and compared them to the genomes of H37Rv and other respiratory Mtb genomes either downloaded from public databases or extracted from local sputum isolates. We aimed to find genomic features that might be distinctly different between CSF-derived and respiratory Mtb. RESULTS: Genome-wide comparisons revealed rearrangements (translocations, inversions, insertions and deletions) and non-synonymous SNPs in our CSF-derived strains that were not observed in the respiratory Mtb genomes used for comparison. These rearranged segments were rich in genes for PE (proline-glutamate)/PPE (proline-proline-glutamate), transcriptional and membrane proteins. Similarly, most of the ns SNPs common in CSF strains were noted in genes encoding PE/PPE proteins. Protein globularity differences were observed among mycobacteria from CSF and respiratory sources and in proteins previously reported to be associated with TB meningitis. Transcription factors and other transcription regulators featured prominently in these proteins. Homologs of proteins associated with Streptococcus pneumoniae meningitis and Neisseria meningitidis virulence were identified in neuropathogenic as well as respiratory mycobacterial spp. examined in this study. DISCUSSION: The occurrence of in silico genetic differences in CSF-derived but not respiratory Mtb suggests their possible involvement in the pathogenesis of TBM. However, overall findings in this comparative analysis support the postulation that TB meningeal infection is more likely to be related to the expression of multiple virulence factors on interaction with host defences than to CNS tropism associated with specific genetic traits.

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.

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.

Identification of new genomospecies in the Mycobacterium terrae complex.

Members of the Mycobacterium terrae complex are slow-growing, non-chromogenic acid-fast bacilli found in the natural environment and occasionally in clinical material. These genetically closely-related members are difficult to differentiate by conventional phenotypic and molecular tests. In this paper we describe the use of whole genome data for the identification of four strains genetically similar to Mycobacterium sp. JDM601, a newly identified member of the M. terrae complex. Phylogenetic information from the alignment of genome-wide orthologous genes and single nucleotide polymorphisms show consistent clustering of the four strains together with M. sp. JDM601 into a distinct clade separate from other rapid and slow growing mycobacterial species. More detailed inter-strain comparisons using average nucleotide identity, tetra-nucleotide frequencies and analysis of synteny indicate that our strains are closely related to but not of the same species as M. sp. JDM601. Besides the 16S rRNA signature described previously for the M. terrae complex, five more hypothetical proteins were found that are potentially useful for the rapid identification of mycobacterial species belonging to the M. terrae complex. This paper illustrates the versatile utilization of whole genome data for the delineation of new bacterial species and introduces four new genomospecies to add to current members in the M. terrae complex.

Comparative genomic analysis of Mycobacterium iranicum UM_TJL against representative mycobacterial species suggests its environmental origin.

Mycobacterium iranicum is a newly reported mycobacterial species. We present the first comparative study of M. iranicum UM_TJL and other mycobacteria. We found M. iranicum to have a close genetic association with environmental mycobacteria infrequently associated with human infections. Nonetheless, UM_TJL is also equipped with many virulence genes (some of which appear to be the consequence of transduction-related gene transfer) that have been identified in established human pathogens. Taken all together, our data suggest that M. iranicum is an environmental bacterium adapted for pathogenicity in the human host. This comparative study provides important clues and forms the basis for future functional studies on this mycobacterium.

Genome Sequence of Fusobacterium nucleatum Strain W1481, a Possible New Subspecies Isolated from a Periodontal Pocket.

Fusobacterium nucleatum is a bacterial species commonly detected in dental plaque within the human oral cavity, with some strains associated with periodontal disease, one of the most common clinical bacterial infections in the human body. The exact mechanisms of its pathogenesis are still not completely understood. In this study, we present the genome sequence and annotation of F. nucleatum strain W1481, isolated from a periodontal pocket of a dental patient at the University of Bristol, United Kingdom, the 16S rRNA gene sequencing of which showed it to be markedly different from the five previously named subspecies.

Genomic reconnaissance of clinical isolates of emerging human pathogen Mycobacterium abscessus reveals high evolutionary potential.

Mycobacterium abscessus (Ma) is an emerging human pathogen that causes both soft tissue infections and systemic disease. We present the first comparative whole-genome study of Ma strains isolated from patients of wide geographical origin. We found a high proportion of accessory strain-specific genes indicating an open, non-conservative pan-genome structure, and clear evidence of rapid phage-mediated evolution. Although we found fewer virulence factors in Ma compared to M. tuberculosis, our data indicated that Ma evolves rapidly and therefore should be monitored closely for the acquisition of more pathogenic traits. This comparative study provides a better understanding of Ma and forms the basis for future functional work on this important pathogen.

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

BACKGROUND: Mycobacterium abscessus is a rapidly growing mycobacterium that is often associated with human infections. The taxonomy of this species has undergone several revisions and is still being debated. In this study, we sequenced the genomes of 12 M. abscessus strains and used phylogenomic analysis to perform subspecies classification. RESULTS: A data mining approach was used to rank and select informative genes based on the relative entropy metric for the construction of a phylogenetic tree. The resulting tree topology was similar to that generated using the concatenation of five classical housekeeping genes: rpoB, hsp65, secA, recA and sodA. Additional support for the reliability of the subspecies classification came from the analysis of erm41 and ITS gene sequences, single nucleotide polymorphisms (SNPs)-based classification and strain clustering demonstrated by a variable number tandem repeat (VNTR) assay and a multilocus sequence analysis (MLSA). We subsequently found that the concatenation of a minimal set of three median-ranked genes: DNA polymerase III subunit alpha (polC), 4-hydroxy-2-ketovalerate aldolase (Hoa) and cell division protein FtsZ (ftsZ), is sufficient to recover the same tree topology. PCR assays designed specifically for these genes showed that all three genes could be amplified in the reference strain of M. abscessus ATCC 19977T. CONCLUSION: This study provides proof of concept that whole-genome sequence-based data mining approach can provide confirmatory evidence of the phylogenetic informativeness of existing markers, as well as lead to the discovery of a more economical and informative set of markers that produces similar subspecies classification in M. abscessus. The systematic procedure used in this study to choose the informative minimal set of gene markers can potentially be applied to species or subspecies classification of other bacteria.

Genome Sequence of Parvimonas micra Strain A293, Isolated from an Abdominal Abscess from a Patient in the United Kingdom.

Parvimonas micra is an important oral microbe that has the ability to grow and proliferate within oral biofilms and is involved in periodontal disease, leading to gingival bleeding, gingival recession, alveolar bone loss, and tooth mobility. However, occasionally these normally oral pathogens can cause infections at other sites in the body. We present the genome sequence of Parvimonas micra strain A293, a smooth Parvimonas micra strain isolated from an abdominal abscess from a patient at Barts Hospital, London, United Kingdom.

First Whole-Genome Sequence of Mycobacterium iranicum, a Newly Reported Mycobacterial Species.

Mycobacterium iranicum is a new species of nontuberculous mycobacterium reported in 2013. Here, we describe the first whole-genome sequence of this species, that of M. iranicum strain UM_TJL, isolated from a patient in Malaysia.

Draft Genome Sequence of Mycobacterium massiliense Strain M159, Showing Phenotypic Resistance to ß-Lactam and Tetracycline Antibiotics.

Mycobacterium massiliense is a nontuberculous mycobacterium associated with human infections. We report here the draft genome sequence of M. massiliense strain M159, isolated from the bronchial aspirate of a patient who had a pulmonary infection. This strain showed genotypic and in vitro resistance to a number of tetracyclines and beta-lactam antibiotics.

Genome sequence of Mycobacterium abscessus strain M152.

Mycobacterium abscessus is an environmental bacterium with increasing clinical relevance. Here, we report the annotated whole-genome sequence of M. abscessus strain M152.

Genomic analysis of Mycobacterium abscessus strain M139, which has an ambiguous subspecies taxonomic position.

Mycobacterium abscessus is a ubiquitous, rapidly growing species of nontuberculous mycobacteria that colonizes organic surfaces and is frequently associated with opportunistic infections in humans. We report here the draft genome sequence of Mycobacterium abscessus strain M139, which shows genomic features reported to be characteristic of both Mycobacterium abscessus subsp. abscessus and Mycobacterium abscessus subsp. massiliense.

Annotated genome sequence of Mycobacterium massiliense strain M154, belonging to the recently created taxon Mycobacterium abscessus subsp. bolletii comb. nov.

Mycobacterium massiliense has recently been proposed as a member of Mycobacterium abscessus subsp. bolletii comb. nov. Strain M154, a clinical isolate from the bronchoalveolar lavage fluid of a Malaysian patient presenting with lower respiratory tract infection, was subjected to shotgun DNA sequencing with the Illumina sequencing technology to obtain whole-genome sequence data for comparison with other genetically related strains within the M. abscessus species complex.