Comparative genomics of Mycobacterium mucogenicum and Mycobacterium neoaurum clade members emphasizing tRNA and non-coding RNA

The discovery of a catalytic RNA within RNase P and its legacy

Sidney Altman's discovery of the processing of one RNA by another RNA that acts like an enzyme was revolutionary in biology and the basis for his sharing the 1989 Nobel Prize in Chemistry with Thomas Cech. These breakthrough findings support the key role of RNA in molecular evolution, where replicating RNAs (and similar chemical derivatives) either with or without peptides functioned in protocells during the early stages of life on Earth, an era referred to as the RNA world. Here, we cover the historical background highlighting the work of Altman and his colleagues and the subsequent efforts of other researchers to understand the biological function of RNase P and its catalytic RNA subunit and to employ it as a tool to downregulate gene expression. We primarily discuss bacterial RNase P-related studies but acknowledge that many groups have significantly contributed to our understanding of archaeal and eukaryotic RNase P, as reviewed in this special issue and elsewhere.

RIP-seq reveals RNAs that interact with RNA polymerase and primary sigma factors in bacteria

Bacteria have evolved structured RNAs that can associate with RNA polymerase (RNAP). Two of them have been known so far-6S RNA and Ms1 RNA but it is unclear if any other types of RNAs binding to RNAP exist in bacteria. To identify all RNAs interacting with RNAP and the primary σ factors, we have established and performed native RIP-seq in Bacillus subtilis, Corynebacterium glutamicum, Streptomyces coelicolor, Mycobacterium smegmatis and the pathogenic Mycobacterium tuberculosis. Besides known 6S RNAs in B. subtilis and Ms1 in M. smegmatis, we detected MTS2823, a homologue of Ms1, on RNAP in M. tuberculosis. In C. glutamicum, we discovered novel types of structured RNAs that associate with RNAP. Furthermore, we identified other species-specific RNAs including full-length mRNAs, revealing a previously unknown landscape of RNAs interacting with the bacterial transcription machinery.

Importance of residue 248 in Escherichia coli RNase P RNA mediated cleavage

tRNA genes are transcribed as precursors and RNase P generates the matured 5' end of tRNAs. It has been suggested that residue - 1 (the residue immediately 5' of the scissile bond) in the pre-tRNA interacts with the well-conserved bacterial RNase P RNA (RPR) residue A248 (Escherichia coli numbering). The way A248 interacts with residue - 1 is not clear. To gain insight into the role of A248, we analyzed cleavage as a function of A248 substitutions and N-1 nucleobase identity by using pre-tRNA and three model substrates. Our findings are consistent with a model where the structural topology of the active site varies and depends on the identity of the nucleobases at, and in proximity to, the cleavage site and their potential to interact. This leads to positioning of Mg2+ that activates the water that acts as the nucleophile resulting in efficient and correct cleavage. We propose that in addition to be involved in anchoring the substrate the role of A248 is to exclude bulk water from access to the amino acid acceptor stem, thereby preventing non-specific hydrolysis of the pre-tRNA. Finally, base stacking is discussed as a way to protect functionally important base-pairing interactions from non-specific hydrolysis, thereby ensuring high fidelity during RNA processing and the decoding of mRNA.

Comparative genome analysis of mycobacteria focusing on tRNA and non-coding RNA

BACKGROUND

The Mycobacterium genus encompasses at least 192 named species, many of which cause severe diseases such as tuberculosis. Non-tuberculosis mycobacteria (NTM) can also infect humans and animals. Some are of emerging concern because they show high resistance to commonly used antibiotics while others are used and evaluated in bioremediation or included in anticancer vaccines.

RESULTS

We provide the genome sequences for 114 mycobacterial type strains and together with 130 available mycobacterial genomes we generated a phylogenetic tree based on 387 core genes and supported by average nucleotide identity (ANI) data. The 244 genome sequences cover most of the species constituting the Mycobacterium genus. The genome sizes ranged from 3.2 to 8.1 Mb with an average of 5.7 Mb, and we identified 14 new plasmids. Moreover, mycobacterial genomes consisted of phage-like sequences ranging between 0 and 4.64% dependent on mycobacteria while the number of IS elements varied between 1 and 290. Our data also revealed that, depending on the mycobacteria, the number of tRNA and non-coding (nc) RNA genes differ and that their positions on the chromosome varied. We identified a conserved core set of 12 ncRNAs, 43 tRNAs and 18 aminoacyl-tRNA synthetases among mycobacteria.

CONCLUSIONS

Phages, IS elements, tRNA and ncRNAs appear to have contributed to the evolution of the Mycobacterium genus where several tRNA and ncRNA genes have been horizontally transferred. On the basis of our phylogenetic analysis, we identified several isolates of unnamed species as new mycobacterial species or strains of known mycobacteria. The predicted number of coding sequences correlates with genome size while the number of tRNA, rRNA and ncRNA genes does not. Together these findings expand our insight into the evolution of the Mycobacterium genus and as such they establish a platform to understand mycobacterial pathogenicity, their evolution, antibiotic resistance/tolerance as well as the function and evolution of ncRNA among mycobacteria.

Ms1 RNA Interacts With the RNA Polymerase Core in Streptomyces coelicolor and Was Identified in Majority of Actinobacteria Using a Linguistic Gene Synteny Search

Bacteria employ small non-coding RNAs (sRNAs) to regulate gene expression. Ms1 is an sRNA that binds to the RNA polymerase (RNAP) core and affects the intracellular level of this essential enzyme. Ms1 is structurally related to 6S RNA that binds to a different form of RNAP, the holoenzyme bearing the primary sigma factor. 6S RNAs are widespread in the bacterial kingdom except for the industrially and medicinally important Actinobacteria. While Ms1 RNA was identified in Mycobacterium, it is not clear whether Ms1 RNA is present also in other Actinobacteria species. Here, using a computational search based on secondary structure similarities combined with a linguistic gene synteny approach, we identified Ms1 RNA in Streptomyces. In S. coelicolor, Ms1 RNA overlaps with the previously annotated scr3559 sRNA with an unknown function. We experimentally confirmed that Ms1 RNA/scr3559 associates with the RNAP core without the primary sigma factor HrdB in vivo. Subsequently, we applied the computational approach to other Actinobacteria and identified Ms1 RNA candidates in 824 Actinobacteria species, revealing Ms1 RNA as a widespread class of RNAP binding sRNAs, and demonstrating the ability of our multifactorial computational approach to identify weakly conserved sRNAs in evolutionarily distant genomes.

CRISPR-Cas12a assisted precise genome editing of Mycolicibacterium neoaurum

Efficient and convenient genetic manipulation of mycobacteria, important microorganisms in human healthcare and the pharmaceutical industry, is limited. In this study, using a model strain Mycolicibacterium neoaurum ATCC 25795, the classical bacterium for the production of valuable steroidal pharmaceuticals, a genome editing system employing CRISPR-Cas12a to achieve efficient and precise genetic manipulation has been developed. Targeted genome mutations could be easily achieved by the CRISPR-Cas12a system without exogenous donor templates, assisted by innate non-homologous end-joining (NHEJ). CRISPR-Cas12a enabled rapid one-step genomic DNA fragment deletions of 1 kb, 5 kb, 10 kb, 15 kb, 20 kb and 24 kb with efficiencies of 70 %, 30 %, 30 %, 20 %, 20 % and 10 %, respectively. Combined with the pNIL/pGOAL system, CRISPR-Cas12a successfully integrated the gene of interest into the targeted genomic site by single crossover and double crossovers with efficiencies of 100 % and 9 %, respectively, using a two-plasmid system. The robust CRISPR systems developed demonstrated strong potential for precise genome editing in M. neoaurum, including targeted deletion of DNA sequences of various lengths and integration of targeted genes into desired sites in the genome.

6S-Like scr3559 RNA Affects Development and Antibiotic Production in Streptomyces coelicolor

Regulatory RNAs control a number of physiological processes in bacterial cells. Here we report on a 6S-like RNA transcript (scr3559) that affects both development and antibiotic production in Streptomyces coelicolor. Its expression is enhanced during the transition to stationary phase. Strains that over-expressed the scr3559 gene region exhibited a shortened exponential growth phase in comparison with a control strain; accelerated aerial mycelium formation and spore maturation; alongside an elevated production of actinorhodin and undecylprodigiosin. These observations were supported by LC-MS analyses of other produced metabolites, including: germicidins, desferrioxamines, and coelimycin. A subsequent microarray differential analysis revealed increased expression of genes associated with the described morphological and physiological changes. Structural and functional similarities between the scr3559 transcript and 6S RNA, and its possible employment in regulating secondary metabolite production are discussed.

Mycobacterium mucogenicum Infection in a Patient with an Open Fracture: A Case Report

Mycobacterium mucogenicum is a nontuberculous mycobacterium that is ubiquitous in nature. However, M. mucogenicum infection in patients with orthopedic trauma is rarely reported in the literature. Herein, we describe a 48 year old male Han Chinese patient whose right leg was squeezed by agricultural machinery, resulting in open tibial fractures. Postoperative antimicrobial treatment was administered because the wound had been contaminated by soil. However, no long-term wound closure occurred, and a culture of the wound exudation tested positive for M. mucogenicum. We established the clinical treatment plan according to the characteristics and drug sensitivity test results of M. mucogenicum, and the patient was discharged uneventfully. Increasingly, more reports of infection caused by nontuberculous mycobacteria are being published; however, to our knowledge, this is the first report of an orthopedic infection caused by M. mucogenicum. Because the treatment process of M. mucogenicum infection is long and complex, isolation and identification of M. mucogenicum are of great significance to effective clinical treatment.

Pulmonary disease caused by a newly identified mycobacterium: Mycolicibacterium toneyamachuris: a case report

BACKGROUND

Non-tuberculous mycobacterial pulmonary disease (NTM-PD) is becoming a significant health burden. Recent advances in analysis techniques have allowed the accurate identification of previously unknown NTM species. Here, we report a case of NTM-PD caused by a newly identified mycobacteria in an immunocompetent patient.

CASE PRESENTATION

A 44-year-old woman was referred to our hospital due to the frequent aggravation of her chronic respiratory symptoms, with NTM-PD-compatible computed tomography findings. Unidentified mycobacterium was repeatedly isolated from respiratory specimens and we diagnosed her as NTM-PD of unidentified mycobacterium. Subsequent whole-genome analysis revealed that the unidentified mycobacterium was a novel mycobacterium genetically close to Mycolicibacterium mucogenicum. We started combination therapy with clarithromycin, moxifloxacin, amikacin, and imipenem/cilastatin, referring to drug sensitivity test results and observed its effect on M. mucogenicum infection. Her symptoms and radiological findings improved significantly.

CONCLUSION

We report a case of NTM-PD caused by a newly identified mycobacteria, Mycolicibacterium toneyamachuris, genetically close to M. mucogenicum. This pathogenic mycobacterium showed different characteristics from M. mucogenicum about clinical presentation and drug sensitivity. The clinical application of genomic sequencing will advance the identification and classification of pathogenic NTM species, and enhance our understanding of mycobacterial diseases.

The rare lncRNA GOLLD is widespread and structurally conserved among Mycobacterium tRNA arrays

Noncoding RNA (ncRNA) genes produce transcripts involved in a wide range of functions, including catalytic and regulatory functions. Besides, some transcripts have highly complex structures that may impact their activities. Among the largest bacterial ncRNAs, there is the rare GOLLD RNA, which is associated with tRNA genes and supposed to be chromosome- and phage-encoded in specialized groups of bacteria, including those from Lactobacillales and Actinomycetales orders. The only GOLLD structure was inferred from a variety of sequences, including many marine metagenomes. To explore GOLLD RNA in bacterial genomes, we mined the GOLLD gene in thousands of Mycobacterium and virus genomes using Infernal software. We identified this gene in 350 mycobacteria, including megaplasmids, and 39 bacteriophages, mainly in the genomic context of tRNA arrays. Mycobacterium GOLLD genes presented a high diversity and were distributed in three phylogenetic groups: (i) Mycobacterium exclusive; (ii) Mycobacterium and mycobacteriophages; and (iii) mycobacteriophage exclusive. We also determined the GOLLD secondary structure of each group using R2 R software based on GOLLD alignments generated by Infernal software. All GOLLD groups displayed a 3' half conserved structure, including utter E-loops pseudoknots substructures, also shared by non-Mycobacterium GOLLD while the 5' half motif was different among the groups. Here, we showed that the lncRNA GOLLD is widespread in Mycobacterium within tRNA arrays and corroborated the previously predicted GOLLD secondary structure.

Insight into the biology of Mycobacterium mucogenicum and Mycobacterium neoaurum clade members

Nontuberculous mycobacteria, NTM, are of growing concern and among these members of the Mycobacterium mucogenicum (Mmuc) and Mycobacterium neoaurum (Mneo) clades can cause infections in humans and they are resistant to first-line anti-tuberculosis drugs. They can be isolated from different ecological niches such as soil, tap water and ground water. Mycobacteria, such as Mmuc and Mneo, are classified as rapid growing mycobacteria, RGM, while the most familiar, Mycobacterium tuberculosis, belongs to the slow growing mycobacteria, SGM. Modern “omics” approaches have provided new insights into our understanding of the biology and evolution of this group of bacteria. Here we present comparative genomics data for seventeen NTM of which sixteen belong to the Mmuc- and Mneo-clades. Focusing on virulence genes, including genes encoding sigma/anti-sigma factors, serine threonine protein kinases (STPK), type VII (ESX genes) secretion systems and mammalian cell entry (Mce) factors we provide insight into their presence as well as phylogenetic relationship in the case of the sigma/anti-sigma factors and STPKs. Our data further suggest that these NTM lack ESX-5 and Mce2 genes, which are known to affect virulence. In this context, Mmuc- and Mneo-clade members lack several of the genes in the glycopeptidolipid (GLP) locus, which have roles in colony morphotype appearance and virulence. For the M. mucogenicum type strain, Mmuc^T, we provide RNASeq data focusing on mRNA levels for sigma factors, STPK, ESX proteins and Mce proteins. These data are discussed and compared to in particular the SGM and fish pathogen Mycobacterium marinum. Finally, we provide insight into as to why members of the Mmuc- and Mneo-clades show resistance to rifampin and isoniazid, and why Mmuc^T forms a rough colony morphotype.