Genomewide comparison and novel ncRNAs of Aquificales

Multiple structural flavors of RNase P in precursor tRNA processing

The precursor transfer RNAs (pre-tRNAs) require extensive processing to generate mature tRNAs possessing proper fold, structural stability, and functionality required to sustain cellular viability. The road to tRNA maturation follows an ordered process: 5'-processing, 3'-processing, modifications at specific sites, if any, and 3'-CCA addition before aminoacylation and recruitment to the cellular protein synthesis machinery. Ribonuclease P (RNase P) is a universally conserved endonuclease in all domains of life, performing the hydrolysis of pre-tRNA sequences at the 5' end by the removal of phosphodiester linkages between nucleotides at position -1 and +1. Except for an archaeal species: Nanoarchaeum equitans where tRNAs are transcribed from leaderless-position +1, RNase P is indispensable for life and displays fundamental variations in terms of enzyme subunit composition, mechanism of substrate recognition and active site architecture, utilizing in all cases a two metal ion-mediated conserved catalytic reaction. While the canonical RNA-based ribonucleoprotein RNase P has been well-known to occur in bacteria, archaea, and eukaryotes, the occurrence of RNA-free protein-only RNase P in eukaryotes and RNA-free homologs of Aquifex RNase P in prokaryotes has been discovered more recently. This review aims to provide a comprehensive overview of structural diversity displayed by various RNA-based and RNA-free RNase P holoenzymes towards harnessing critical RNA-protein and protein-protein interactions in achieving conserved pre-tRNA processing functionality. Furthermore, alternate roles and functional interchangeability of RNase P are discussed in the context of its employability in several clinical and biotechnological applications. This article is categorized under: RNA Processing > tRNA Processing RNA Evolution and Genomics > RNA and Ribonucleoprotein Evolution RNA Interactions with Proteins and Other Molecules > RNA-Protein Complexes.

Clustered Regularly Interspaced Short Palindromic Repeat/CRISPR-Associated Protein and Its Utility All at Sea: Status, Challenges, and Prospects

Initially discovered over 35 years ago in the bacterium Escherichia coli as a defense system against invasion of viral (or other exogenous) DNA into the genome, CRISPR/Cas has ushered in a new era of functional genetics and served as a versatile genetic tool in all branches of life science. CRISPR/Cas has revolutionized the methodology of gene knockout with simplicity and rapidity, but it is also powerful for gene knock-in and gene modification. In the field of marine biology and ecology, this tool has been instrumental in the functional characterization of 'dark' genes and the documentation of the functional differentiation of gene paralogs. Powerful as it is, challenges exist that have hindered the advances in functional genetics in some important lineages. This review examines the status of applications of CRISPR/Cas in marine research and assesses the prospect of quickly expanding the deployment of this powerful tool to address the myriad fundamental marine biology and biological oceanography questions.

Characterization of RNA-based and protein-only RNases P from bacteria encoding both enzyme types

A small group of bacteria encode two types of RNase P, the classical ribonucleoprotein (RNP) RNase P as well as the protein-only RNase P HARP (homolog of Aquifex RNase P). We characterized the dual RNase P activities of five bacteria that belong to three different phyla. All five bacterial species encode functional RNA (gene rnpB) and protein (gene rnpA) subunits of RNP RNase P, but only the HARP of the thermophile Thermodesulfatator indicus (phylum Thermodesulfobacteria) was found to have robust tRNA 5'-end maturation activity in vitro and in vivo in an Escherichia coli RNase P depletion strain. These findings suggest that both types of RNase P are able to contribute to the essential tRNA 5'-end maturation activity in T. indicus, thus resembling the predicted evolutionary transition state in the progenitor of the Aquificaceae before the loss of rnpA and rnpB genes in this family of bacteria. Remarkably, T. indicus RNase P RNA is transcribed with a P12 expansion segment that is posttranscriptionally excised in vivo, such that the major fraction of the RNA is fragmented and thereby truncated by ∼70 nt in the native T. indicus host as well as in the E. coli complementation strain. Replacing the native P12 element of T. indicus RNase P RNA with the short P12 helix of Thermotoga maritima RNase P RNA abolished fragmentation, but simultaneously impaired complementation efficiency in E. coli cells, suggesting that intracellular fragmentation and truncation of T. indicus RNase P RNA may be beneficial to RNA folding and/or enzymatic activity.

Insights into 6S RNA in lactic acid bacteria (LAB)

BACKGROUND

6S RNA is a regulator of cellular transcription that tunes the metabolism of cells. This small non-coding RNA is found in nearly all bacteria and among the most abundant transcripts. Lactic acid bacteria (LAB) constitute a group of microorganisms with strong biotechnological relevance, often exploited as starter cultures for industrial products through fermentation. Some strains are used as probiotics while others represent potential pathogens. Occasional reports of 6S RNA within this group already indicate striking metabolic implications. A conceivable idea is that LAB with 6S RNA defects may metabolize nutrients faster, as inferred from studies of Echerichia coli. This may accelerate fermentation processes with the potential to reduce production costs. Similarly, elevated levels of secondary metabolites might be produced. Evidence for this possibility comes from preliminary findings regarding the production of surfactin in Bacillus subtilis, which has functions similar to those of bacteriocins. The prerequisite for its potential biotechnological utility is a general characterization of 6S RNA in LAB.

RESULTS

We provide a genomic annotation of 6S RNA throughout the Lactobacillales order. It laid the foundation for a bioinformatic characterization of common 6S RNA features. This covers secondary structures, synteny, phylogeny, and product RNA start sites. The canonical 6S RNA structure is formed by a central bulge flanked by helical arms and a template site for product RNA synthesis. 6S RNA exhibits strong syntenic conservation. It is usually flanked by the replication-associated recombination protein A and the universal stress protein A. A catabolite responsive element was identified in over a third of all 6S RNA genes. It is known to modulate gene expression based on the available carbon sources. The presence of antisense transcripts could not be verified as a general trait of LAB 6S RNAs.

CONCLUSIONS

Despite a large number of species and the heterogeneity of LAB, the stress regulator 6S RNA is well-conserved both from a structural as well as a syntenic perspective. This is the first approach to describe 6S RNAs and short 6S RNA-derived transcripts beyond a single species, spanning a large taxonomic group covering multiple families. It yields universal insights into this regulator and complements the findings derived from other bacterial model organisms.

Heuristic algorithms for best match graph editing

BACKGROUND

Best match graphs (BMGs) are a class of colored digraphs that naturally appear in mathematical phylogenetics as a representation of the pairwise most closely related genes among multiple species. An arc connects a gene x with a gene y from another species (vertex color) Y whenever it is one of the phylogenetically closest relatives of x. BMGs can be approximated with the help of similarity measures between gene sequences, albeit not without errors. Empirical estimates thus will usually violate the theoretical properties of BMGs. The corresponding graph editing problem can be used to guide error correction for best match data. Since the arc set modification problems for BMGs are NP-complete, efficient heuristics are needed if BMGs are to be used for the practical analysis of biological sequence data.

RESULTS

Since BMGs have a characterization in terms of consistency of a certain set of rooted triples (binary trees on three vertices) defined on the set of genes, we consider heuristics that operate on triple sets. As an alternative, we show that there is a close connection to a set partitioning problem that leads to a class of top-down recursive algorithms that are similar to Aho's supertree algorithm and give rise to BMG editing algorithms that are consistent in the sense that they leave BMGs invariant. Extensive benchmarking shows that community detection algorithms for the partitioning steps perform best for BMG editing.

CONCLUSION

Noisy BMG data can be corrected with sufficient accuracy and efficiency to make BMGs an attractive alternative to classical phylogenetic methods.

Minimal protein-only RNase P structure reveals insights into tRNA precursor recognition and catalysis

Ribonuclease P (RNase P) is an endoribonuclease that catalyzes the processing of the 5' leader sequence of precursor tRNA (pre-tRNA). Ribonucleoprotein RNase P and protein-only RNase P (PRORP) in eukaryotes have been extensively studied, but the mechanism by which a prokaryotic nuclease recognizes and cleaves pre-tRNA is unclear. To gain insights into this mechanism, we studied homologs of Aquifex RNase P (HARPs), thought to be enzymes of approximately 23 kDa comprising only this nuclease domain. We determined the cryo-EM structure of Aq880, the first identified HARP enzyme. The structure unexpectedly revealed that Aq880 consists of both the nuclease and protruding helical (PrH) domains. Aq880 monomers assemble into a dimer via the PrH domain. Six dimers form a dodecamer with a left-handed one-turn superhelical structure. The structure also revealed that the active site of Aq880 is analogous to that of eukaryotic PRORPs. The pre-tRNA docking model demonstrated that 5' processing of pre-tRNAs is achieved by two adjacent dimers within the dodecamer. One dimer is responsible for catalysis, and the PrH domains of the other dimer are responsible for pre-tRNA elbow recognition. Our study suggests that HARPs measure an invariant distance from the pre-tRNA elbow to cleave the 5' leader sequence, which is analogous to the mechanism of eukaryotic PRORPs and the ribonucleoprotein RNase P. Collectively, these findings shed light on how different types of RNase P enzymes utilize the same pre-tRNA processing.

Microbe Profile: Aquifex aeolicus: an extreme heat-loving bacterium that feeds on gases and inorganic chemicals

The bacterium ‘ Aquifex aeolicus ’ is the model organism for the deeply rooted phylum Aquificae . This ‘water-maker’ is an H2-oxidizing microaerophile that flourishes in extremely hot marine habitats, and it also thrives on the sulphur compounds commonly found in volcanic environments. ‘ A. aeolicus ’ has hyper-stable proteins and a fully sequenced genome, with some of its essential metabolic pathways deciphered (including energy conservation). Many of its proteins have also been characterized (especially structurally), including many of the enzymes involved in replication, transcription, RNA processing and cell envelope biosynthesis. Enzymes that are of promise for biotechnological applications have been widely investigated in this species. ‘ A. aeolicus ’ has also added to our understanding of the origins of life and evolution.

6S RNA, a Global Regulator of Transcription

6S RNA is a small RNA regulator of RNA polymerase (RNAP) that is present broadly throughout the bacterial kingdom. Initial functional studies in Escherichia coli revealed that 6S RNA forms a complex with RNAP resulting in regulation of transcription, and cells lacking 6S RNA have altered survival phenotypes. The last decade has focused on deepening the understanding of several aspects of 6S RNA activity, including (i) addressing questions of how broadly conserved 6S RNAs are in diverse organisms through continued identification and initial characterization of divergent 6S RNAs; (ii) the nature of the 6S RNA-RNAP interaction through examination of variant proteins and mutant RNAs, cross-linking approaches, and ultimately a cryo-electron microscopic structure; (iii) the physiological consequences of 6S RNA function through identification of the 6S RNA regulon and promoter features that determine 6S RNA sensitivity; and (iv) the mechanism and cellular impact of 6S RNA-directed synthesis of product RNAs (i.e., pRNA synthesis). Much has been learned about this unusual RNA, its mechanism of action, and how it is regulated; yet much still remains to be investigated, especially regarding potential differences in behavior of 6S RNAs in diverse bacteria.

Minimal and RNA-free RNase P in Aquifex aeolicus

RNase P is an essential tRNA-processing enzyme in all domains of life. We identified an unknown type of protein-only RNase P in the hyperthermophilic bacterium Aquifex aeolicus: Without an RNA subunit and the smallest of its kind, the 23-kDa polypeptide comprises a metallonuclease domain only. The protein has RNase P activity in vitro and rescued the growth of Escherichia coli and Saccharomyces cerevisiae strains with inactivations of their more complex and larger endogenous ribonucleoprotein RNase P. Homologs of Aquifex RNase P (HARP) were identified in many Archaea and some Bacteria, of which all Archaea and most Bacteria also encode an RNA-based RNase P; activity of both RNase P forms from the same bacterium or archaeon could be verified in two selected cases. Bioinformatic analyses suggest that A. aeolicus and related Aquificaceae likely acquired HARP by horizontal gene transfer from an archaeon.

Dissemination of 6S RNA among bacteria

6S RNA is a highly abundant small non-coding RNA widely spread among diverse bacterial groups. By competing with DNA promoters for binding to RNA polymerase (RNAP), the RNA regulates transcription on a global scale. RNAP produces small product RNAs derived from 6S RNA as template, which rearranges the 6S RNA structure leading to dissociation of 6S RNA:RNAP complexes. Although 6S RNA has been experimentally analysed in detail for some species, such as Escherichia coli and Bacillus subtilis, and was computationally predicted in many diverse bacteria, a complete and up-to-date overview of the distribution among all bacteria is missing. In this study we searched with new methods for 6S RNA genes in all currently available bacterial genomes. We ended up with a set of 1,750 6S RNA genes, of which 1,367 are novel and bona fide, distributed among 1,610 bacteria, and had a few tentative candidates among the remaining 510 assembled bacterial genomes accessible. We were able to confirm two tentative candidates by Northern blot analysis. We extended 6S RNA genes of the Flavobacteriia significantly in length compared to the present Rfam entry. We describe multiple homologs of 6S RNAs (including split 6S RNA genes) and performed a detailed synteny analysis.

Comprehensive insights in the Mycobacterium avium subsp. paratuberculosis genome using new WGS data of sheep strain JIII-386 from Germany

Mycobacterium avium (M. a.) subsp. paratuberculosis (MAP)—the etiologic agent of Johne’s disease—affects cattle, sheep, and other ruminants worldwide. To decipher phenotypic differences among sheep and cattle strains (belonging to MAP-S [Type-I/III], respectively, MAP-C [Type-II]), comparative genome analysis needs data from diverse isolates originating from different geographic regions of the world. This study presents the so far best assembled genome of a MAP-S-strain: Sheep isolate JIII-386 from Germany. One newly sequenced cattle isolate (JII-1961, Germany), four published MAP strains of MAP-C and MAP-S from the United States and Australia, and M. a. subsp. hominissuis (MAH) strain 104 were used for assembly improvement and comparisons. All genomes were annotated by BacProt and results compared with NCBI (National Center for Biotechnology Information) annotation. Corresponding protein-coding sequences (CDSs) were detected, but also CDSs that were exclusively determined by either NCBI or BacProt. A new Shine–Dalgarno sequence motif (5′-AGCTGG-3′) was extracted. Novel CDSs including PE-PGRS family protein genes and about 80 noncoding RNAs exhibiting high sequence conservation are presented. Previously found genetic differences between MAP-types are partially revised. Four of ten assumed MAP-S-specific large sequence polymorphism regions (LSP^Ss) are still present in MAP-C strains; new LSP^Ss were identified. Independently of the regional origin of the strains, the number of individual CDSs and single nucleotide variants confirms the strong similarity of MAP-C strains and shows higher diversity among MAP-S strains. This study gives ambiguous results regarding the hypothesis that MAP-S is the evolutionary intermediate between MAH and MAP-C, but it clearly shows a higher similarity of MAP to MAH than to Mycobacterium intracellulare.

Phylogenomics with paralogs

Phylogenomics heavily relies on well-curated sequence data sets that comprise, for each gene, exclusively 1:1 orthologos. Paralogs are treated as a dangerous nuisance that has to be detected and removed. We show here that this severe restriction of the data sets is not necessary. Building upon recent advances in mathematical phylogenetics, we demonstrate that gene duplications convey meaningful phylogenetic information and allow the inference of plausible phylogenetic trees, provided orthologs and paralogs can be distinguished with a degree of certainty. Starting from tree-free estimates of orthology, cograph editing can sufficiently reduce the noise to find correct event-annotated gene trees. The information of gene trees can then directly be translated into constraints on the species trees. Although the resolution is very poor for individual gene families, we show that genome-wide data sets are sufficient to generate fully resolved phylogenetic trees, even in the presence of horizontal gene transfer.