A tribute to Sidney Altman, one of the architects of modern RNA biology.

This special issue of JBC pays tribute to Sidney Altman, whose discovery of a catalytic role for RNA, a breakthrough made independently by Thomas Cech, overturned the long-held dogma that only proteins can serve as catalysts in biological systems. The discovery of RNA catalysis galvanized biologists to think expansively in new directions and has given rise to a remarkable RNAissance in science and medicine. The collection of articles begins with the story of the discovery of RNase P and builds up to the emerging picture of an unexpectedly vast repertoire of RNase P variants in the three domains of life, including insights derived from recent high-resolution structures on how RNAs, ribonucleoproteins, or protein scaffolds can be used variably to generate an active site for catalyzing the same RNA processing reaction. The series of articles ends with a discussion of more recently discovered endonucleases (Argonautes, Cas), whose parallels with RNase P underscore recurring themes in diverse biological contexts.

Suppression of the E. coli rnpA49 conditionally lethal phenotype by different compensatory mutations.

RNase P is an essential enzyme found across all domains of life that is responsible for the 5'-end maturation of precursor tRNAs. For decades, numerous studies have sought to elucidate the mechanisms and biochemistry governing RNase P function. However, much remains unknown about the regulation of RNase P expression, the turnover and degradation of the enzyme, and the mechanisms underlying the phenotypes and complementation of specific RNase P mutations, especially in the model bacterium, Escherichia coli. In E. coli, the temperature-sensitive rnpA49 mutation in the protein subunit of RNase P has arguably been one of the most well-studied mutations for examining the enzyme's activity in vivo. Here, we report for the first time naturally-occurring temperature-resistant suppressor mutations of E. coli strains carrying the rnpA49 allele. We find that rnpA49 strains can partially compensate the temperature-sensitive defect via gene amplifications of either RNase P subunit (rnpA49 or rnpB) or by the acquisition of loss-of-function mutations in Lon protease or RNase R. Our results agree with previous plasmid overexpression and gene deletion complementation studies, and importantly suggest the involvement of Lon protease in the degradation and/or regulatory pathway(s) of the mutant protein subunit of RNase P. This work offers novel insights into the behavior and complementation of the rnpA49 allele in vivo and provides direction for follow-up studies regarding RNase P regulation and turnover in E. coli.

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.

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.

Intracellular localization of the mycobacterial stressosome complex.

Microorganisms survive stresses by alternating the expression of genes suitable for surviving the immediate and present danger and eventually adapt to new conditions. Many bacteria have evolved a multiprotein "molecular machinery" designated the "Stressosome" that integrates different stress signals and activates alternative sigma factors for appropriate downstream responses. We and others have identified orthologs of some of the Bacillus subtilis stressosome components, RsbR, RsbS, RsbT and RsbUVW in several mycobacteria and we have previously reported mutual interactions among the stressosome components RsbR, RsbS, RsbT and RsbUVW from Mycobacterium marinum. Here we provide evidence that "STAS" domains of both RsbR and RsbS are important for establishing the interaction and thus critical for stressosome assembly. Fluorescence microscopy further suggested co-localization of RsbR and RsbS in multiprotein complexes visible as co-localized fluorescent foci distributed at scattered locations in the M. marinum cytoplasm; the number, intensity and distribution of such foci changed in cells under stressed conditions. Finally, we provide bioinformatics data that 17 (of 244) mycobacteria, which lack the RsbRST genes, carry homologs of Bacillus cereus genes rsbK and rsbM indicating the existence of alternative σF activation pathways among mycobacteria.

Author Correction: Extensive genomic diversity among Mycobacterium marinum strains revealed by whole genome sequencing.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

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, MmucT, 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 MmucT forms a rough colony morphotype.

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

BACKGROUND: Mycobacteria occupy various ecological niches and can be isolated from soil, tap water and ground water. Several cause diseases in humans and animals. To get deeper insight into our understanding of mycobacterial evolution focusing on tRNA and non-coding (nc)RNA, we conducted a comparative genome analysis of Mycobacterium mucogenicum (Mmuc) and Mycobacterium neoaurum (Mneo) clade members. RESULTS: Genome sizes for Mmuc- and Mneo-clade members vary between 5.4 and 6.5 Mbps with the complete MmucT (type strain) genome encompassing 6.1 Mbp. The number of tRNA genes range between 46 and 79 (including one pseudo tRNA gene) with 39 tRNA genes common among the members of these clades, while additional tRNA genes were probably acquired through horizontal gene transfer. Selected tRNAs and ncRNAs (RNase P RNA, tmRNA, 4.5S RNA, Ms1 RNA and 6C RNA) are expressed, and the levels for several of these are higher in stationary phase compared to exponentially growing cells. The rare tRNAIleTAT isoacceptor and two for mycobacteria novel ncRNAs: the Lactobacillales-derived GOLLD RNA and a homolog to the antisense Salmonella typhimurium phage Sar RNA, were shown to be present and expressed in certain Mmuc-clade members. CONCLUSIONS: Phages, IS elements, horizontally transferred tRNA gene clusters, and phage-derived ncRNAs appears to have influenced the evolution of the Mmuc- and Mneo-clades. While the number of predicted coding sequences correlates with genome size, the number of tRNA coding genes does not. The majority of the tRNA genes in mycobacteria are transcribed mainly from single genes and the levels of certain ncRNAs, including RNase P RNA (essential for the processing of tRNAs), are higher at stationary phase compared to exponentially growing cells. We provide supporting evidence that Ms1 RNA represents a mycobacterial 6S RNA variant. The evolutionary routes for the ncRNAs RNase P RNA, tmRNA and Ms1 RNA are different from that of the core genes.

Gene miaA for post-transcriptional modification of tRNAXXA is important for morphological and metabolic differentiation in Streptomyces.

Members of actinobacterial genus Streptomyces possess a sophisticated life cycle and are the deepest source of bioactive secondary metabolites. Although morphogenesis and secondary metabolism are subject to transcriptional co-regulation, streptomycetes employ an additional mechanism to initiate the aforementioned processes. This mechanism is based on delayed translation of rare leucyl codon UUA by the only cognate tRNALeu UAA (encoded by bldA). The bldA-based genetic switch is an extensively documented example of translational regulation in Streptomyces. Yet, after five decades since the discovery of bldA, factors that shape its function and peculiar conditionality remained elusive. Here we address the hypothesis that post-transcriptional tRNA modifications play a role in tRNA-based mechanisms of translational control in Streptomyces. Particularly, we studied two Streptomyces albus J1074 genes, XNR_1074 (miaA) and XNR_1078 (miaB), encoding tRNA (adenosine(37)-N6)-dimethylallyltransferase and tRNA (N6-isopentenyl adenosine(37)-C2)-methylthiotransferase respectively. These enzymes produce, in a sequential manner, a hypermodified ms2 i6 A37 residue in most of the A36-A37-containing tRNAs. We show that miaB and especially miaA null mutant of S. albus possess altered morphogenesis and secondary metabolism. We provide genetic evidence that miaA deficiency impacts translational level of gene expression, most likely through impaired decoding of codons UXX and UUA in particular.

Extended insight into the Mycobacterium chelonae-abscessus complex through whole genome sequencing of Mycobacterium salmoniphilum outbreak and Mycobacterium salmoniphilum-like strains.

Members of the Mycobacterium chelonae-abscessus complex (MCAC) are close to the mycobacterial ancestor and includes both human, animal and fish pathogens. We present the genomes of 14 members of this complex: the complete genomes of Mycobacterium salmoniphilum and Mycobacterium chelonae type strains, seven M. salmoniphilum isolates, and five M. salmoniphilum-like strains including strains isolated during an outbreak in an animal facility at Uppsala University. Average nucleotide identity (ANI) analysis and core gene phylogeny revealed that the M. salmoniphilum-like strains are variants of the human pathogen Mycobacterium franklinii and phylogenetically close to Mycobacterium abscessus. Our data further suggested that M. salmoniphilum separates into three branches named group I, II and III with the M. salmoniphilum type strain belonging to group II. Among predicted virulence factors, the presence of phospholipase C (plcC), which is a major virulence factor that makes M. abscessus highly cytotoxic to mouse macrophages, and that M. franklinii originally was isolated from infected humans make it plausible that the outbreak in the animal facility was caused by a M. salmoniphilum-like strain. Interestingly, M. salmoniphilum-like was isolated from tap water suggesting that it can be present in the environment. Moreover, we predicted the presence of mutational hotspots in the M. salmoniphilum isolates and 26% of these hotspots overlap with genes categorized as having roles in virulence, disease and defense. We also provide data about key genes involved in transcription and translation such as sigma factor, ribosomal protein and tRNA genes.

Extensive genomic diversity among Mycobacterium marinum strains revealed by whole genome sequencing.

Mycobacterium marinum is the causative agent for the tuberculosis-like disease mycobacteriosis in fish and skin lesions in humans. Ubiquitous in its geographical distribution, M. marinum is known to occupy diverse fish as hosts. However, information about its genomic diversity is limited. Here, we provide the genome sequences for 15 M. marinum strains isolated from infected humans and fish. Comparative genomic analysis of these and four available genomes of the M. marinum strains M, E11, MB2 and Europe reveal high genomic diversity among the strains, leading to the conclusion that M. marinum should be divided into two different clusters, the "M"- and the "Aronson"-type. We suggest that these two clusters should be considered to represent two M. marinum subspecies. Our data also show that the M. marinum pan-genome for both groups is open and expanding and we provide data showing high number of mutational hotspots in M. marinum relative to other mycobacteria such as Mycobacterium tuberculosis. This high genomic diversity might be related to the ability of M. marinum to occupy different ecological niches.

Critical domain interactions for type A RNase P RNA catalysis with and without the specificity domain.

The natural trans-acting ribozyme RNase P RNA (RPR) is composed of two domains in which the catalytic (C-) domain mediates cleavage of various substrates. The C-domain alone, after removal of the second specificity (S-) domain, catalyzes this reaction as well, albeit with reduced efficiency. Here we provide experimental evidence indicating that efficient cleavage mediated by the Escherichia coli C-domain (Eco CP RPR) with and without the C5 protein likely depends on an interaction referred to as the "P6-mimic". Moreover, the P18 helix connects the C- and S-domains between its loop and the P8 helix in the S-domain (the P8/ P18-interaction). In contrast to the "P6-mimic", the presence of P18 does not contribute to the catalytic performance by the C-domain lacking the S-domain in cleavage of an all ribo model hairpin loop substrate while deletion or disruption of the P8/ P18-interaction in full-size RPR lowers the catalytic efficiency in cleavage of the same model hairpin loop substrate in keeping with previously reported data using precursor tRNAs. Consistent with that P18 is not required for cleavage mediated by the C-domain we show that the archaeal Pyrococcus furiosus RPR C-domain, which lacks the P18 helix, is catalytically active in trans without the S-domain and any protein. Our data also suggest that the S-domain has a larger impact on catalysis for E. coli RPR compared to P. furiosus RPR. Finally, we provide data indicating that the absence of the S-domain and P18, or the P8/ P18-interaction in full-length RPR influences the charge distribution near the cleavage site in the RPR-substrate complex to a small but reproducible extent.

Inhibition of Bacterial RNase P RNA by Phenothiazine Derivatives.

There is a need to identify novel scaffolds and targets to develop new antibiotics. Methylene blue is a phenothiazine derivative, and it has been shown to possess anti-malarial and anti-trypanosomal activities. Here, we show that different phenothiazine derivatives and pyronine G inhibited the activities of three structurally different bacterial RNase P RNAs (RPRs), including that from Mycobacterium tuberculosis, with Ki values in the lower µM range. Interestingly, three antipsychotic phenothiazines (chlorpromazine, thioridazine, and trifluoperazine), which are known to have antibacterial activities, also inhibited the activity of bacterial RPRs, albeit with higher Ki values than methylene blue. Phenothiazines also affected lead(II)-induced cleavage of bacterial RPR and inhibited yeast tRNA(Phe), indicating binding of these drugs to functionally important regions. Collectively, our findings provide the first experimental data showing that long, noncoding RNAs could be targeted by different phenothiazine derivatives.

Cleavage of Model Substrates by Arabidopsis thaliana PRORP1 Reveals New Insights into Its Substrate Requirements.

Two broad classes of RNase P trim the 5' leader of precursor tRNAs (pre-tRNAs): ribonucleoprotein (RNP)- and proteinaceous (PRORP)-variants. These two RNase P types, which use different scaffolds for catalysis, reflect independent evolutionary paths. While the catalytic RNA-based RNP form is present in all three domains of life, the PRORP family is restricted to eukaryotes. To obtain insights on substrate recognition by PRORPs, we examined the 5' processing ability of recombinant Arabidopsis thaliana PRORP1 (AtPRORP1) using a panel of pre-tRNASer variants and model hairpin-loop derivatives (pATSer type) that consist of the acceptor-T-stem stack and the T-/D-loop. Our data indicate the importance of the identity of N-1 (the residue immediately 5' to the cleavage site) and the N-1:N+73 base pair for cleavage rate and site selection of pre-tRNASer and pATSer. The nucleobase preferences that we observed mirror the frequency of occurrence in the complete suite of organellar pre-tRNAs in eight algae/plants that we analyzed. The importance of the T-/D-loop in pre-tRNASer for tight binding to AtPRORP1 is indicated by the 200-fold weaker binding of pATSer compared to pre-tRNASer, while the essentiality of the T-loop for cleavage is reflected by the near-complete loss of activity when a GAAA-tetraloop replaced the T-loop in pATSer. Substituting the 2'-OH at N-1 with 2'-H also resulted in no detectable cleavage, hinting at the possible role of this 2'-OH in coordinating Mg2+ ions critical for catalysis. Collectively, our results indicate similarities but also key differences in substrate recognition by the bacterial RNase P RNP and AtPRORP1: while both forms exploit the acceptor-T-stem stack and the elbow region in the pre-tRNA, the RNP form appears to require more recognition determinants for cleavage-site selection.

The Mycobacterium phlei Genome: Expectations and Surprises.

Mycobacterium phlei, a nontuberculosis mycobacterial species, was first described in 1898-1899. We present the complete genome sequence for theM. phlei CCUG21000(T)type strain and the draft genomes for four additional strains. The genome size for all five is 5.3 Mb with 69.4% Guanine-Cytosine content. This is ≈0.35 Mbp smaller than the previously reported M. phlei RIVM draft genome. The size difference is attributed partly to large bacteriophage sequence fragments in theM. phlei RIVM genome. Comparative analysis revealed the following: 1) A CRISPR system similar to Type 1E (cas3) in M. phlei RIVM; 2) genes involved in polyamine metabolism and transport (potAD,potF) that are absent in other mycobacteria, and 3) strain-specific variations in the number of σ-factor genes. Moreover,M. phlei has as many as 82 mce(mammalian cell entry) homologs and many of the horizontally acquired genes in M. phlei are present in other environmental bacteria including mycobacteria that share similar habitat. Phylogenetic analysis based on 693 Mycobacterium core genes present in all complete mycobacterial genomes suggested that its closest neighbor is Mycobacterium smegmatis JS623 and Mycobacterium rhodesiae NBB3, while it is more distant toM. smegmatis mc2 155.

Comparative Sigma Factor-mRNA Levels in Mycobacterium marinum under Stress Conditions and during Host Infection.

We have used RNASeq and qRT-PCR to study mRNA levels for all σ-factors in different Mycobacterium marinum strains under various growth and stress conditions. We also studied their levels in M. marinum from infected fish and mosquito larvae. The annotated σ-factors were expressed and transcripts varied in relation to growth and stress conditions. Some were highly abundant such as sigA, sigB, sigC, sigD, sigE and sigH while others were not. The σ-factor mRNA profiles were similar after heat stress, during infection of fish and mosquito larvae. The similarity also applies to some of the known heat shock genes such as the α-crystallin gene. Therefore, it seems probable that the physiological state of M. marinum is similar when exposed to these different conditions. Moreover, the mosquito larvae data suggest that this is the state that the fish encounter when infected, at least with respect to σ-factor mRNA levels. Comparative genomic analysis of σ-factor gene localizations in three M. marinum strains and Mycobacterium tuberculosis H37Rv revealed chromosomal rearrangements that changed the localization of especially sigA, sigB, sigD, sigE, sigF and sigJ after the divergence of these two species. This may explain the variation in species-specific expression upon exposure to different growth conditions.

Characterization of Three Mycobacterium spp. with Potential Use in Bioremediation by Genome Sequencing and Comparative Genomics.

We provide the genome sequences of the type strains of the polychlorophenol-degrading Mycobacterium chlorophenolicum (DSM43826), the degrader of chlorinated aliphatics Mycobacterium chubuense (DSM44219) and Mycobacterium obuense (DSM44075) that has been tested for use in cancer immunotherapy. The genome sizes of M. chlorophenolicum, M. chubuense, and M. obuense are 6.93, 5.95, and 5.58 Mb with GC-contents of 68.4%, 69.2%, and 67.9%, respectively. Comparative genomic analysis revealed that 3,254 genes are common and we predicted approximately 250 genes acquired through horizontal gene transfer from different sources including proteobacteria. The data also showed that the biodegrading Mycobacterium spp. NBB4, also referred to as M. chubuense NBB4, is distantly related to the M. chubuense type strain and should be considered as a separate species, we suggest it to be named Mycobacterium ethylenense NBB4. Among different categories we identified genes with potential roles in: biodegradation of aromatic compounds and copper homeostasis. These are the first nonpathogenic Mycobacterium spp. found harboring genes involved in copper homeostasis. These findings would therefore provide insight into the role of this group of Mycobacterium spp. in bioremediation as well as the evolution of copper homeostasis within the Mycobacterium genus.

Differentiation and phylogenetic relationships in Mycobacterium spp with special reference to the RNase P RNA gene rnpB.

The rnpB gene encodes for the RNA subunit of the catalytic ribonuclease RNase P and is present in all bacteria and has both conserved and highly variable sequence regions. Determination of rnpB in 35 Mycobacterium spp. showed species specific sequences for all species except the Mycobacterium tuberculosis complex (four species). High sequence variation was seen in the P3, P15 and P19 regions of suggested secondary structures of the corresponding RNase P RNA molecules. Phylogenetic analysis showed that rnpB gave similar tree topologies as 16S rRNA and hsp65 genes. A combined analysis of the three genes increased the number of nodes with significant support from 10 to 19. The results indicate that rnpB is useful for phylogenetic studies and is a possible target for identification and detection of Mycobacterium spp.