Leaderless genes in bacteria: clue to the evolution of translation initiation mechanisms in prokaryotes

Diversity and Evolution of Mitochondrial Translation Apparatus

The evolution of mitochondria has proceeded independently in different eukaryotic lines, which is reflected in the diversity of mitochondrial genomes and mechanisms of their expression in eukaryotic species. Mitochondria have lost most of bacterial ancestor genes by transferring them to the nucleus or eliminating them. However, mitochondria of almost all eukaryotic cells still retain relatively small genomes, as well as their replication, transcription, and translation apparatuses. The dependence on the nuclear genome, specific features of mitochondrial transcripts, and synthesis of highly hydrophobic membrane proteins in the mitochondria have led to significant changes in the translation apparatus inherited from the bacterial ancestor, which retained the basic structure necessary for protein synthesis but became more specialized and labile. In this review, we discuss specific properties of translation initiation in the mitochondria and how the evolution of mitochondria affected the functions of main factors initiating protein biosynthesis in these organelles.

Comprehensive quantitative modeling of translation efficiency in a genome-reduced bacterium

Translation efficiency has been mainly studied by ribosome profiling, which only provides an incomplete picture of translation kinetics. Here, we integrated the absolute quantifications of tRNAs, mRNAs, RNA half-lives, proteins, and protein half-lives with ribosome densities and derived the initiation and elongation rates for 475 genes (67% of all genes), 73 with high precision, in the bacterium Mycoplasma pneumoniae (Mpn). We found that, although the initiation rate varied over 160-fold among genes, most of the known factors had little impact on translation efficiency. Local codon elongation rates could not be fully explained by the adaptation to tRNA abundances, which varied over 100-fold among tRNA isoacceptors. We provide a comprehensive quantitative view of translation efficiency, which suggests the existence of unidentified mechanisms of translational regulation in Mpn.

Ultradeep characterisation of translational sequence determinants refutes rare-codon hypothesis and unveils quadruplet base pairing of initiator tRNA and transcript

Translation is a key determinant of gene expression and an important biotechnological engineering target. In bacteria, 5'-untranslated region (5'-UTR) and coding sequence (CDS) are well-known mRNA parts controlling translation and thus cellular protein levels. However, the complex interaction of 5'-UTR and CDS has so far only been studied for few sequences leading to non-generalisable and partly contradictory conclusions. Herein, we systematically assess the dynamic translation from over 1.2 million 5'-UTR-CDS pairs in Escherichia coli to investigate their collective effect using a new method for ultradeep sequence-function mapping. This allows us to disentangle and precisely quantify effects of various sequence determinants of translation. We find that 5'-UTR and CDS individually account for 53% and 20% of variance in translation, respectively, and show conclusively that, contrary to a common hypothesis, tRNA abundance does not explain expression changes between CDSs with different synonymous codons. Moreover, the obtained large-scale data provide clear experimental evidence for a base-pairing interaction between initiator tRNA and mRNA beyond the anticodon-codon interaction, an effect that is often masked for individual sequences and therefore inaccessible to low-throughput approaches. Our study highlights the indispensability of ultradeep sequence-function mapping to accurately determine the contribution of parts and phenomena involved in gene regulation.

CvkR is a MerR-type transcriptional repressor of class 2 type V-K CRISPR-associated transposase systems

Certain CRISPR-Cas elements integrate into Tn7-like transposons, forming CRISPR-associated transposon (CAST) systems. How the activity of these systems is controlled in situ has remained largely unknown. Here we characterize the MerR-type transcriptional regulator Alr3614 that is encoded by one of the CAST (AnCAST) system genes in the genome of cyanobacterium Anabaena sp. PCC 7120. We identify a number of Alr3614 homologs across cyanobacteria and suggest naming these regulators CvkR for Cas V-K repressors. Alr3614/CvkR is translated from leaderless mRNA and represses the AnCAST core modules cas12k and tnsB directly, and indirectly the abundance of the tracr-CRISPR RNA. We identify a widely conserved CvkR binding motif 5'-AnnACATnATGTnnT-3'. Crystal structure of CvkR at 1.6 Å resolution reveals that it comprises distinct dimerization and potential effector-binding domains and that it assembles into a homodimer, representing a discrete structural subfamily of MerR regulators. CvkR repressors are at the core of a widely conserved regulatory mechanism that controls type V-K CAST systems.

A Copper-Responsive Two-Component System Governs Lipoprotein Remodeling in Listeria monocytogenes

Bacterial lipoproteins are membrane-associated proteins with a characteristic acylated N-terminal cysteine residue anchoring C-terminal globular domains to the membrane surface. While all lipoproteins are modified with acyl chains, the number, length, and position can vary depending on host. The acylation pattern also alters ligand recognition by the Toll-like receptor 2 (TLR2) protein family, a signaling system that is central to bacterial surveillance and innate immunity. In select Listeria monocytogenes isolates carrying certain plasmids, copper exposure converts the lipoprotein chemotype into a weak TLR2 ligand through expression of the enzyme lipoprotein intramolecular acyltransferase (Lit). In this study, we identify the response regulator (CopR) from a heavy metal-sensing two-component system as the transcription factor that integrates external copper levels with lipoprotein structural modifications. We show that phosphorylated CopR controls the expression of three distinct transcripts within the plasmid cassette encoding Lit2, prolipoprotein diacylglyceryl transferase (Lgt2), putative copper resistance determinants, and itself (the CopRS two-component system). CopR recognizes a direct repeat half-site consensus motif (TCTACACA) separated by 3 bp that overlaps the -35 promoter element. Target gene expression and lipoprotein conversion were not observed in the absence of the response regulator, indicating that CopR phosphorylation is the dominant mechanism of regulation. IMPORTANCE Copper is a frontline antimicrobial used to limit bacterial growth in multiple settings. Here, we demonstrate how the response regulator CopR from a plasmid-borne two-component system in the opportunistic pathogen L. monocytogenes directly induces lipoprotein remodeling in tandem with copper resistance genes due to extracellular copper stress. Activation of CopR by phosphorylation converts the lipoprotein chemotype from a high- to low-immunostimulatory TLR2 ligand. The two-component system-mediated coregulation of copper resistance determinants, in tandem with lipoprotein biosynthesis demonstrated here in L. monocytogenes, may be a common feature of transmissible copper resistance cassettes found in other Firmicutes.

Initiation at AUGUG and GUGUG sequences can lead to translation of overlapping reading frames in E. coli

During initiation, the ribosome is tasked to efficiently recognize open reading frames (ORFs) for accurate and fast translation of mRNAs. A critical step is start codon recognition, which is modulated by initiation factors, mRNA structure, a Shine Dalgarno (SD) sequence and the start codon itself. Within the Escherichia coli genome, we identified more than 50 annotated initiation sites harboring AUGUG or GUGUG sequence motifs that provide two canonical start codons, AUG and GUG, in immediate proximity. As these sites may challenge start codon recognition, we studied if and how the ribosome is accurately guided to the designated ORF, with a special focus on the SD sequence as well as adenine at the fourth coding sequence position (A4). By in vitro and in vivo experiments, we characterized key requirements for unambiguous start codon recognition, but also discovered initiation sites that lead to the translation of both overlapping reading frames. Our findings corroborate the existence of an ambiguous translation initiation mechanism, implicating a multitude of so far unrecognized ORFs and translation products in bacteria.

Flow-Seq Method: Features and Application in Bacterial Translation Studies

The Flow-seq method is based on using reporter construct libraries, where a certain element regulating the gene expression of fluorescent reporter proteins is represented in many thousands of variants. Reporter construct libraries are introduced into cells, sorted according to their fluorescence level, and then subjected to next-generation sequencing. Therefore, it turns out to be possible to identify patterns that determine the expression efficiency, based on tens and hundreds of thousands of reporter constructs in one experiment. This method has become common in evaluating the efficiency of protein synthesis simultaneously by multiple mRNA variants. However, its potential is not confined to this area. In the presented review, a comparative analysis of the Flow-seq method and other alternative approaches used for translation efficiency evaluation of mRNA was carried out; the features of its application and the results obtained by Flow-seq were also considered.

Regulation of Leaderless mRNA Translation in Bacteria

In bacteria, the translation of genetic information can begin through at least three different mechanisms: canonical or Shine-Dalgarno-led initiation, readthrough or 70S scanning initiation, or leaderless initiation. Here, we discuss the main features and regulation of the last, which is characterized mainly by the ability of 70S ribosomal particles to bind to AUG located at or near the 5′ end of mRNAs to initiate translation. These leaderless mRNAs (lmRNAs) are rare in enterobacteria, such as Escherichia coli, but are common in other bacteria, such as Mycobacterium tuberculosis and Deinococcus deserti, where they may represent more than 20% and even up to 60% of the genes. Given that lmRNAs are devoid of a 5′ untranslated region and the Shine-Dalgarno sequence located within it, the mechanism of translation regulation must depend on molecular strategies that are different from what has been observed in the Shine-Dalgarno-led translation. Diverse regulatory mechanisms have been proposed, including the processing of ribosomal RNA and changes in the abundance of translation factors, but all of them produce global changes in the initiation of lmRNA translation. Thus, further research will be required to understand how the initiation of the translation of particular lmRNA genes is regulated.

Identification and Cross-Characterisation of Artificial Promoters and 5′ Untranslated Regions in Vibrio natriegens

Vibrio natriegens has recently gained attention as a novel fast-growing bacterium in synthetic biology applications. Currently, a limited set of genetic elements optimised for Escherichia coli are used in V. natriegens due to the lack of DNA parts characterised in this novel host. In this study, we report the identification and cross-characterisation of artificial promoters and 5′ untranslated regions (artificial regulatory sequence, ARES) that lead to production of fluorescent proteins with a wide-range of expression levels. We identify and cross-characterise 52 constructs in V. natriegens and E. coli. Furthermore, we report the DNA sequence and motif analysis of the ARESs using various algorithms. With this study, we expand the pool of characterised genetic DNA parts that can be used for different biotechnological applications using V. natriegens as a host microorganism.

Prokaryotic Genome Annotation

In the last decade, the high-throughput and relatively low cost of short-read sequencing technologies have revolutionized prokaryotic genomics. This has led to an exponential increase in the number of bacterial and archaeal genome sequences available, as well as corresponding increase of genome assembly and annotation tools developed. Together, these hardware and software technologies have given scientists unprecedented options to study their chosen microbial systems without the need for large teams of bioinformaticists or supercomputing facilities. While these analysis tools largely fall into only a few categories, each may have different requirements, caveats and file formats, and some may be rarely updated or even abandoned. And so, despite the apparent ease in sequencing and analyzing a prokaryotic genome, it is no wonder that the budding genomicist may quickly find oneself overwhelmed. Here, we aim to provide the reader with an overview of genome annotation and its most important considerations, as well as an easy-to-follow protocol to get started with annotating a prokaryotic genome.

Importance of the 5' regulatory region to bacterial synthetic biology applications

The field of synthetic biology is evolving at a fast pace. It is advancing beyond single-gene alterations in single hosts to the logical design of complex circuits and the development of integrated synthetic genomes. Recent breakthroughs in deep learning, which is increasingly used in de novo assembly of DNA components with predictable effects, are also aiding the discipline. Despite advances in computing, the field is still reliant on the availability of pre-characterized DNA parts, whether natural or synthetic, to regulate gene expression in bacteria and make valuable compounds. In this review, we discuss the different bacterial synthetic biology methodologies employed in the creation of 5' regulatory regions - promoters, untranslated regions and 5'-end of coding sequences. We summarize methodologies and discuss their significance for each of the functional DNA components, and highlight the key advances made in bacterial engineering by concentrating on their flaws and strengths. We end the review by outlining the issues that the discipline may face in the near future.

Translation of a Leaderless Reporter Is Robust During Exponential Growth and Well Sustained During Stress Conditions in Mycobacterium tuberculosis

Mycobacterium tuberculosis expresses a large number of leaderless mRNA transcripts; these lack the 5' leader region, which usually contains the Shine-Dalgarno sequence required for translation initiation in bacteria. In M. tuberculosis, transcripts encoding proteins like toxin-antitoxin systems are predominantly leaderless and the overall ratio of leaderless to Shine-Dalgarno transcripts significantly increases during growth arrest, suggesting that leaderless translation might be important during persistence in the host. However, whether these two types of transcripts are translated with differing efficiencies during optimal growth conditions and during stress conditions that induce growth arrest, is unclear. Here, we have used the desA1 (Rv0824c) and desA2 (Rv1094) gene pair as representative for Shine-Dalgarno and leaderless transcripts in M. tuberculosis respectively; and used them to construct bioluminescent reporter strains. We detect robust leaderless translation during exponential in vitro growth, and we show that leaderless translation is more stable than Shine-Dalgarno translation during adaptation to stress conditions. These changes are independent from transcription, as transcription levels did not significantly change following quantitative real-time PCR analysis. Upon entrance into nutrient starvation and after nitric oxide exposure, leaderless translation is significantly less affected by the stress than Shine-Dalgarno translation. Similarly, during the early stages of infection of macrophages, the levels of leaderless translation are transiently more stable than those of Shine-Dalgarno translation. These results suggest that leaderless translation may offer an advantage in the physiology of M. tuberculosis. Identification of the molecular mechanisms underlying this translational regulation may provide insights into persistent infection.

Characterization of the DdrD protein from the extremely radioresistant bacterium Deinococcus radiodurans

Here, we report the in vitro and in vivo characterization of the DdrD protein from the extraordinary stress-resistant bacterium, D. radiodurans. DdrD is one of the most highly induced proteins following cellular irradiation or desiccation. We confirm that DdrD belongs to the Radiation Desiccation Response (RDR) regulon protein family whose expression is regulated by the IrrE/DdrO proteins after DNA damage. We show that DdrD is a DNA binding protein that binds to single-stranded DNA In vitro, but not to duplex DNA unless it has a 5' single-stranded extension. In vivo, we observed no significant effect of the absence of DdrD on the survival of D. radiodurans cells after exposure to γ-rays or UV irradiation in different genetic contexts. However, genome reassembly is affected in a ∆ddrD mutant when cells recover from irradiation in the absence of nutrients. Thus, DdrD likely contributes to genome reconstitution after irradiation, but only under starvation conditions. Lastly, we show that the absence of the DdrD protein partially restores the frequency of plasmid transformation of a ∆ddrB mutant, suggesting that DdrD could also be involved in biological processes other than the response to DNA damage.

Autorepressor PsrA is required for optimal Legionella pneumophila growth in Acanthamoeba castellanii protozoa

Legionella pneumophila possesses a unique intracellular lifecycle featuring distinct morphological stages that include replicative forms and transmissive cyst forms. Expression of genes associated with virulence traits and cyst morphogenesis is concomitant, and governed by a complex stringent response based-regulatory network and the stationary phase sigma factor RpoS. In Pseudomonas spp., rpoS expression is controlled by the autorepressor PsrA, and orthologs of PsrA and RpoS are required for cyst formation in Azotobacter. Here we report that the L. pneumophila psrA ortholog, expressed as a leaderless monocistronic transcript, is also an autorepressor, but is not a regulator of rpoS expression. Further, the binding site sequence recognized by L. pneumophila PsrA is different from that of Pseudomonas PsrA, suggesting a repertoire of target genes unique to L. pneumophila. While PsrA was dispensable for growth in human U937-derived macrophages, lack of PsrA affected bacterial intracellular growth in Acanthamoeba castellanii protozoa, but also increased the quantity of poly-3-hydroxybutyrate (PHB) inclusions in matured transmissive cysts. Interestingly, overexpression of PsrA increased the size and bacterial load of the replicative vacuole in both host cell types. Taken together, we report that PsrA is a host-specific requirement for optimal temporal progression of L. pneumophila intracellular lifecycle in A. castellanii.

Lysinibacillus sphaericus III(3)7 and Plasmid Vector pMK4: New Challenges in Cloning Platforms

The acquisition and especially the maintenance of a plasmid usually brings a fitness cost that reduces the reproductive rate of the bacterial host; for strains like Lysinibacillus sphaericus III(3)7, which possesses important environmental properties, this alteration along with morphological changes and reduced sporulation rates may exert a negative effect on metabolic studies using plasmids as cloning platforms. The aim of this study is to approach the metabolic behavior of pMK4-bearing cells of L. sphaericus III(3)7 through the use of bioinformatic and in vitro analyses. An incompatibility model between the pMK4 vector and a predicted megaplasmid, pBsph, inside III(3)7 cells was constructed based on an incA region. Additionally, in vitro long-term plasmid stability was not found in plasmid-bearing cells. Alignments between replicons, mobile genetic elements and RNA-RNA interactions were assessed, pairwise alignment visualization, graphic models and morphological changes were evaluated by SEM. Metabolite analysis was done through HPLC coupled to a Q-TOF 6545, and electrospray ionization was used, finally, Aedes aegypti and Culex quinquefasciatus larvae were used for larvicidal activity assessment. Results found, a decreased growth rate, spore formation reduction and morphological changes, which supported the idea of metabolic cost exerted by pMK4. An incompatibility between pMK4 and pBsph appears to take place inside L. sphaericus III(3)7 cells, however, further in vitro studies are needed to confirm it.

Chemotranscriptomic Profiling Defines Drug-Specific Signatures of the Glycopeptide Antibiotics Dalbavancin, Vancomycin and Chlorobiphenyl-Vancomycin in a VanB-Type-Resistant Streptomycete

Dalbavancin, vancomycin and chlorobiphenyl-vancomycin share a high degree of structural similarity and the same primary mode of drug action. All inhibit bacterial cell wall biosynthesis through complexation with intermediates in peptidoglycan biosynthesis mediated via interaction with peptidyl-d-alanyl-d-alanine (d-Ala-d-Ala) residues present at the termini of the intermediates. VanB-type glycopeptide resistance in bacteria encodes an inducible reprogramming of bacterial cell wall biosynthesis that generates precursors terminating with d-alanyl-d-lactate (d-Ala-d-Lac). This system in Streptomyces coelicolor confers protection against the natural product vancomycin but not dalbavancin or chlorobiphenyl-vancomycin, which are semi-synthetic derivatives and fail to sufficiently activate the inducible VanB-type sensory response. We used transcriptome profiling by RNAseq to identify the gene expression signatures elucidated in S. coelicolor in response to the three different glycopeptide compounds. An integrated comparison of the results defines both the contribution of the VanB resistance system to the control of changes in gene transcription and the impact at the transcriptional level of the structural diversity present in the glycopeptide antibiotics used. Dalbavancin induces markedly more extensive changes in the expression of genes required for transport processes, RNA methylation, haem biosynthesis and the biosynthesis of the amino acids arginine and glutamine. Chlorobiphenyl-vancomycin exhibits specific effects on tryptophan and calcium-dependent antibiotic biosynthesis and has a stronger repressive effect on translation. Vancomycin predictably has a uniquely strong effect on the genes controlled by the VanB resistance system and also impacts metal ion homeostasis and leucine biosynthesis. Leaderless gene transcription is disfavoured in the core transcriptional up- and down-regulation taking place in response to all the glycopeptide antibiotics, while HrdB-dependent transcripts are favoured in the down-regulated group. This study illustrates the biological impact of peripheral changes to glycopeptide antibiotic structure and could inform the design of future semi-synthetic glycopeptide derivatives.

Integrative characterization of the near-minimal bacterium Mesoplasma florum

The near-minimal bacterium Mesoplasma florum is an interesting model for synthetic genomics and systems biology due to its small genome (~ 800 kb), fast growth rate, and lack of pathogenic potential. However, fundamental aspects of its biology remain largely unexplored. Here, we report a broad yet remarkably detailed characterization of M. florum by combining a wide variety of experimental approaches. We investigated several physical and physiological parameters of this bacterium, including cell size, growth kinetics, and biomass composition of the cell. We also performed the first genome-wide analysis of its transcriptome and proteome, notably revealing a conserved promoter motif, the organization of transcription units, and the transcription and protein expression levels of all protein-coding sequences. We converted gene transcription and expression levels into absolute molecular abundances using biomass quantification results, generating an unprecedented view of the M. florum cellular composition and functions. These characterization efforts provide a strong experimental foundation for the development of a genome-scale model for M. florum and will guide future genome engineering endeavors in this simple organism.

Lost and Found: Re-searching and Re-scoring Proteomics Data Aids Genome Annotation and Improves Proteome Coverage

Delineation of open reading frames (ORFs) causes persistent inconsistencies in prokaryote genome annotation. We demonstrate that by advanced (re)analysis of omics data, a higher proteome coverage and sensitive detection of unannotated ORFs can be achieved, which can be exploited for conditional bacterial genome (re)annotation, which is especially relevant in view of annotating the wealth of sequenced prokaryotic genomes obtained in recent years.

Prokaryotic genome annotation is heavily dependent on automated gene annotation pipelines that are prone to propagate errors and underestimate genome complexity. We describe an optimized proteogenomic workflow that uses ribosome profiling (ribo-seq) and proteomic data for Salmonella enterica serovar Typhimurium to identify unannotated proteins or alternative protein forms. This data analysis encompasses the searching of cofragmenting peptides and postprocessing with extended peptide-to-spectrum quality features, including comparison to predicted fragment ion intensities. When this strategy is applied, an enhanced proteome depth is achieved, as well as greater confidence for unannotated peptide hits. We demonstrate the general applicability of our pipeline by reanalyzing public Deinococcus radiodurans data sets. Taken together, our results show that systematic reanalysis using available prokaryotic (proteome) data sets holds great promise to assist in experimentally based genome annotation.

IMPORTANCE Delineation of open reading frames (ORFs) causes persistent inconsistencies in prokaryote genome annotation. We demonstrate that by advanced (re)analysis of omics data, a higher proteome coverage and sensitive detection of unannotated ORFs can be achieved, which can be exploited for conditional bacterial genome (re)annotation, which is especially relevant in view of annotating the wealth of sequenced prokaryotic genomes obtained in recent years.

The Transcriptomic Landscape of Cupriavidus metallidurans CH34 Acutely Exposed to Copper

Bacteria are increasingly used for biotechnological applications such as bioremediation, biorecovery, bioproduction, and biosensing. The development of strains suited for such applications requires a thorough understanding of their behavior, with a key role for their transcriptomic landscape. We present a thorough analysis of the transcriptome of Cupriavidus metallidurans CH34 cells acutely exposed to copper by tagRNA-sequencing. C. metallidurans CH34 is a model organism for metal resistance, and its potential as a biosensor and candidate for metal bioremediation has been demonstrated in multiple studies. Several metabolic pathways were impacted by Cu exposure, and a broad spectrum of metal resistance mechanisms, not limited to copper-specific clusters, was overexpressed. In addition, several gene clusters involved in the oxidative stress response and the cysteine-sulfur metabolism were induced. In total, 7500 transcription start sites (TSSs) were annotated and classified with respect to their location relative to coding sequences (CDSs). Predicted TSSs were used to re-annotate 182 CDSs. The TSSs of 2422 CDSs were detected, and consensus promotor logos were derived. Interestingly, many leaderless messenger RNAs (mRNAs) were found. In addition, many mRNAs were transcribed from multiple alternative TSSs. We observed pervasive intragenic TSSs both in sense and antisense to CDSs. Antisense transcripts were enriched near the 5' end of mRNAs, indicating a functional role in post-transcriptional regulation. In total, 578 TSSs were detected in intergenic regions, of which 35 were identified as putative small regulatory RNAs. Finally, we provide a detailed analysis of the main copper resistance clusters in CH34, which include many intragenic and antisense transcripts. These results clearly highlight the ubiquity of noncoding transcripts in the CH34 transcriptome, many of which are putatively involved in the regulation of metal resistance.

Polycysteine-encoding leaderless short ORFs function as cysteine-responsive attenuators of operonic gene expression in mycobacteria

Genome-wide transcriptomic analyses have revealed abundant expressed short open reading frames (ORFs) in bacteria. Whether these short ORFs, or the small proteins they encode, are functional remains an open question. One quarter of mycobacterial mRNAs are leaderless, beginning with a 5'-AUG or GUG initiation codon. Leaderless mRNAs often encode unannotated short ORFs as the first gene of a polycistronic transcript. Here, we show that polycysteine-encoding leaderless short ORFs function as cysteine-responsive attenuators of operonic gene expression. Detailed mutational analysis shows that one polycysteine short ORF controls expression of the downstream genes. Our data indicate that ribosomes stalled in the polycysteine tract block mRNA structures that otherwise sequester the ribosome-binding site of the 3'gene. We assessed endogenous proteomic responses to cysteine limitation in Mycobacterium smegmatis using mass spectrometry. Six cysteine metabolic loci having unannotated polycysteine-encoding leaderless short ORF architectures responded to cysteine limitation, revealing widespread cysteine-responsive attenuation in mycobacteria. Individual leaderless short ORFs confer independent operon-level control, while their shared dependence on cysteine ensures a collective response mediated by ribosome pausing. We propose the term ribulon to classify ribosome-directed regulons. Regulon-level coordination by ribosomes on sensory short ORFs illustrates one utility of the many unannotated short ORFs expressed in bacterial genomes.

Loading messenger RNAs onto ribosomes in human mitochondria: lessons learned from a bacterial toxin

Mitochondria are peculiar organelles because their function depends on genetic information that is present in two genomes: nuclear and mitochondrial. The expression of mitochondrially encoded information requires dedicated machinery. Many efforts have been made to identify this machinery and describe its relevant mechanisms. Recently, Bruni et al. reported a cellular model that they established to investigate the pathway for loading messenger RNAs onto ribosomes in human mitochondria. Their study revealed a role for monosome formation in the stability of mitochondrial mRNAs. Comment on: https://doi.org/10.1111/febs.15342.

Primary transcriptome and translatome analysis determines transcriptional and translational regulatory elements encoded in the Streptomyces clavuligerus genome

Determining transcriptional and translational regulatory elements in GC-rich Streptomyces genomes is essential to elucidating the complex regulatory networks that govern secondary metabolite biosynthetic gene cluster (BGC) expression. However, information about such regulatory elements has been limited for Streptomyces genomes. To address this limitation, a high-quality genome sequence of β-lactam antibiotic-producing Streptomyces clavuligerus ATCC 27 064 is completed, which contains 7163 newly annotated genes. This provides a fundamental reference genome sequence to integrate multiple genome-scale data types, including dRNA-Seq, RNA-Seq and ribosome profiling. Data integration results in the precise determination of 2659 transcription start sites which reveal transcriptional and translational regulatory elements, including −10 and −35 promoter components specific to sigma (σ) factors, and 5′-untranslated region as a determinant for translation efficiency regulation. Particularly, sequence analysis of a wide diversity of the −35 components enables us to predict potential σ-factor regulons, along with various spacer lengths between the −10 and −35 elements. At last, the primary transcriptome landscape of the β-lactam biosynthetic pathway is analyzed, suggesting temporal changes in metabolism for the synthesis of secondary metabolites driven by transcriptional regulation. This comprehensive genetic information provides a versatile genetic resource for rational engineering of secondary metabolite BGCs in Streptomyces.

A TetR-Family Protein (CAETHG_0459) Activates Transcription From a New Promoter Motif Associated With Essential Genes for Autotrophic Growth in Acetogens

Acetogens can fix carbon (CO or CO2) into acetyl-CoA via the Wood-Ljungdahl pathway (WLP) that also makes them attractive cell factories for the production of fuels and chemicals from waste feedstocks. Although most biochemical details of the WLP are well understood and systems-level characterization of acetogen metabolism has recently improved, key transcriptional features such as promoter motifs and transcriptional regulators are still unknown in acetogens. Here, we use differential RNA-sequencing to identify a previously undescribed promoter motif associated with essential genes for autotrophic growth of the model-acetogen Clostridium autoethanogenum. RNA polymerase was shown to bind to the new promoter motif using a DNA-binding protein assay and proteomics enabled the discovery of four candidates to potentially function directly in control of transcription of the WLP and other key genes of C1 fixation metabolism. Next, in vivo experiments showed that a TetR-family transcriptional regulator (CAETHG_0459) and the housekeeping sigma factor (σA) activate expression of a reporter protein (GFP) in-frame with the new promoter motif from a fusion vector in Escherichia coli. Lastly, a protein-protein interaction assay with the RNA polymerase (RNAP) shows that CAETHG_0459 directly binds to the RNAP. Together, the data presented here advance the fundamental understanding of transcriptional regulation of C1 fixation in acetogens and provide a strategy for improving the performance of gas-fermenting bacteria by genetic engineering.

DeepRibo: a neural network for precise gene annotation of prokaryotes by combining ribosome profiling signal and binding site patterns

Annotation of gene expression in prokaryotes often finds itself corrected due to small variations of the annotated gene regions observed between different (sub)-species. It has become apparent that traditional sequence alignment algorithms, used for the curation of genomes, are not able to map the full complexity of the genomic landscape. We present DeepRibo, a novel neural network utilizing features extracted from ribosome profiling information and binding site sequence patterns that shows to be a precise tool for the delineation and annotation of expressed genes in prokaryotes. The neural network combines recurrent memory cells and convolutional layers, adapting the information gained from both the high-throughput ribosome profiling data and ribosome binding translation initiation sequence region into one model. DeepRibo is designed as a single model trained on a variety of ribosome profiling experiments, used for the identification of open reading frames in prokaryotes without a priori knowledge of the translational landscape. Through extensive validation of the model trained on various sets of data, multiple species sequence similarity, mass spectrometry and Edman degradation verified proteins, the effectiveness of DeepRibo is highlighted.

The influence of sigma factors and ribosomal recognition elements on heterologous expression of cyanobacterial gene clusters in Escherichia coli

Cyanobacterial natural products offer new possibilities for drugs and lead compounds but many factors can inhibit the production of sufficient yields for pharmaceutical processes. While Escherichia coli and Streptomyces sp. have been used as heterologous expression hosts to produce cyanobacterial natural products, they have not met with resounding success largely due to their inability to recognize cyanobacterial promoter regions. Recent work has shown that the filamentous freshwater cyanobacterium Anabaena sp. strain PCC 7120 recognizes various cyanobacterial promoter regions and can produce lyngbyatoxin A from the native promoter. Introduction of Anabaena sigma factors into E. coli might allow the native transcriptional machinery to recognize cyanobacterial promoters. Here, all 12 Anabaena sigma factors were expressed in E. coli and subsets were found to initiate transcription from several cyanobacterial promoters based on transcriptional fusions to the chloramphenicol acetyltransferase (CAT) reporter. Expression of individual Anabaena sigma factors in E. coli did not result in lyngbyatoxin A production from its native cyanobacterial gene cluster, possibly hindered by deficiencies in recognition of cyanobacterial ribosomal binding sites by native E. coli translational machinery. This represents an important step toward engineering E. coli into a general heterologous expression host for cyanobacterial biosynthetic gene cluster expression.

Translational regulation in mycobacteria and its implications for pathogenicity

Protein synthesis is a fundamental requirement of all cells for survival and replication. To date, vast numbers of genetic and biochemical studies have been performed to address the mechanisms of translation and its regulation in Escherichia coli, but only a limited number of studies have investigated these processes in other bacteria, particularly in slow growing bacteria like Mycobacterium tuberculosis, the causative agent of human tuberculosis. In this Review, we highlight important differences in the translational machinery of M. tuberculosis compared with E. coli, specifically the presence of two additional proteins and subunit stabilizing elements such as the B9 bridge. We also consider the role of leaderless translation in the ability of M. tuberculosis to establish latent infection and look at the experimental evidence that translational regulatory mechanisms operate in mycobacteria during stress adaptation, particularly focussing on differences in toxin-antitoxin systems between E. coli and M. tuberculosis and on the role of tuneable translational fidelity in conferring phenotypic antibiotic resistance. Finally, we consider the implications of these differences in the context of the biological adaptation of M. tuberculosis and discuss how these regulatory mechanisms could aid in the development of novel therapeutics for tuberculosis.

Noncanonical features and modifications on the 5'-end of bacterial sRNAs and mRNAs

Although many eukaryotic transcripts contain cap structures, it has been long thought that bacterial RNAs do not carry any special modifications on their 5'-ends. In bacteria, primary transcripts are produced by transcription initiated with a nucleoside triphosphate and are therefore triphosphorylated on 5'-ends. Some transcripts are then processed by nucleases that yield monophosphorylated RNAs for specific cellular activities. Many primary transcripts are also converted to monophosphorylated species by removal of the terminal pyrophosphate for 5'-end-dependent degradation. Recent studies surprisingly revealed an expanded repertoire of chemical groups on 5'-ends of bacterial RNAs. In addition to mono- and triphosphorylated moieties, some mRNAs and sRNAs contain cap-like structures and diphosphates on their 5'-ends. Although incorporation and removal of these groups have become better understood in recent years, the physiological significance of these modifications remain obscure. This review highlights recent studies aimed at identification and elucidation of novel modifications on the 5'-ends of bacterial RNAs and discusses possible physiological applications of the modified RNAs. This article is categorized under: RNA Turnover and Surveillance > Regulation of RNA Stability RNA Structure and Dynamics > RNA Structure, Dynamics, and Chemistry RNA Processing > Capping and 5' End Modifications.

Translation in Prokaryotes

This review summarizes our current understanding of translation in prokaryotes, focusing on the mechanistic and structural aspects of each phase of translation: initiation, elongation, termination, and ribosome recycling. The assembly of the initiation complex provides multiple checkpoints for messenger RNA (mRNA) and start-site selection. Correct codon-anticodon interaction during the decoding phase of elongation results in major conformational changes of the small ribosomal subunit and shapes the reaction pathway of guanosine triphosphate (GTP) hydrolysis. The ribosome orchestrates proton transfer during peptide bond formation, but requires the help of elongation factor P (EF-P) when two or more consecutive Pro residues are to be incorporated. Understanding the choreography of transfer RNA (tRNA) and mRNA movements during translocation helps to place the available structures of translocation intermediates onto the time axis of the reaction pathway. The nascent protein begins to fold cotranslationally, in the constrained space of the polypeptide exit tunnel of the ribosome. When a stop codon is reached at the end of the coding sequence, the ribosome, assisted by termination factors, hydrolyzes the ester bond of the peptidyl-tRNA, thereby releasing the nascent protein. Following termination, the ribosome is dissociated into subunits and recycled into another round of initiation. At each step of translation, the ribosome undergoes dynamic fluctuations between different conformation states. The aim of this article is to show the link between ribosome structure, dynamics, and function.

Leaderless mRNAs in the Spotlight: Ancient but Not Outdated!

Previously, leaderless mRNAs (lmRNAs) were perceived to make up only a minor fraction of the transcriptome in bacteria. However, advancements in RNA sequencing technology are uncovering vast numbers of lmRNAs, particularly in archaea, Actinobacteria , and extremophiles and thus underline their significance in cellular physiology and regulation. Due to the absence of conventional ribosome binding signals, lmRNA translation initiation is distinct from canonical mRNAs and can therefore be differentially regulated. The ribosome’s inherent ability to bind a 5′-terminal AUG can stabilize and protect the lmRNA from degradation or allow ribosomal loading for downstream initiation events. As a result, lmRNAs remain translationally competent during a variety of physiological conditions, allowing them to contribute to multiple regulatory mechanisms. Furthermore, the abundance of lmRNAs can increase during adverse conditions through the upregulation of lmRNA transcription from alternative promoters or by the generation of lmRNAs from canonical mRNAs cleaved by an endonucleolytic toxin. In these ways, lmRNA translation can continue during stress and contribute to regulation, illustrating their importance in the cell. Due to their presence in all domains of life and their ability to be translated by heterologous hosts, lmRNAs appear further to represent ancestral transcripts that might allow us to study the evolution of the ribosome and the translational process.

Translation and Translational Control in Dinoflagellates

Dinoflagellates are unicellular protists that feature a multitude of unusual nuclear features, including large genomes, packaging of DNA without histones, and multiple gene copies organized as tandem gene arrays. Furthermore, all dinoflagellate mRNAs experience trans-splicing with a common 22-nucleotide splice leader (SL) sequence. These features challenge some of the concepts and assumptions about the regulation of gene expression derived from work on model eukaryotes such as yeasts and mammals. Translational control in the dinoflagellates, based on extensive study of circadian bioluminescence and by more recent microarray and transcriptome analyses, is now understood to be a crucial element in regulating gene expression. A picture of the translation machinery of dinoflagellates is emerging from the recent availability of transcriptomes of multiple dinoflagellate species and the first complete genome sequences. The components comprising the translational control toolkit of dinoflagellates are beginning to take shape and are outlined here.

PKS-NRPS Enzymology and Structural Biology: Considerations in Protein Production

The structural diversity and complexity of marine natural products have made them a rich and productive source of new bioactive molecules for drug development. The identification of these new compounds has led to extensive study of the protein constituents of the biosynthetic pathways from the producing microbes. Essential processes in the dissection of biosynthesis have been the elucidation of catalytic functions and the determination of 3D structures for enzymes of the polyketide synthases and nonribosomal peptide synthetases that carry out individual reactions. The size and complexity of these proteins present numerous difficulties in the process of going from gene to structure. Here, we review the problems that may be encountered at the various steps of this process and discuss some of the solutions devised in our and other labs for the cloning, production, purification, and structure solution of complex proteins using Escherichia coli as a heterologous host.

Transcriptome sequencing of the human pathogen Corynebacterium diphtheriae NCTC 13129 provides detailed insights into its transcriptional landscape and into DtxR-mediated transcriptional regulation

Background

The human pathogen Corynebacterium diphtheriae is the causative agent of diphtheria. In the 1990s a large diphtheria outbreak in Eastern Europe was caused by the strain C. diphtheriae NCTC 13129. Although the genome was sequenced more than a decade ago, not much is known about its transcriptome. Our aim was to use transcriptome sequencing (RNA-Seq) to close this knowledge gap and gain insights into the transcriptional landscape of a C. diphtheriae tox⁺ strain.

Results

We applied two different RNA-Seq techniques, one to retrieve 5′-ends of primary transcripts and the other to characterize the whole transcriptional landscape in order to gain insights into various features of the C. diphtheriae NCTC 13129 transcriptome. By examining the data we identified 1656 transcription start sites (TSS), of which 1202 were assigned to genes and 454 to putative novel transcripts. By using the TSS data promoter regions recognized by the housekeeping sigma factor σ^A and its motifs were analyzed in detail, revealing a well conserved −10 but an only weakly conserved −35 motif, respectively. Furthermore, with the TSS data 5’-UTR lengths were explored. The observed 5’-UTRs range from zero length (leaderless transcripts), which make up 20% of all genes, up to over 450 nt long leaders, which may harbor regulatory functions. The C. diphtheriae transcriptome consists of 471 operons which are further divided into 167 sub-operon structures. In a differential expression analysis approach, we discovered that genetic disruption of the iron-sensing transcription regulator DtxR, which controls expression of diphtheria toxin (DT), causes a strong influence on general gene expression. Nearly 15% of the genome is differentially transcribed, indicating that DtxR might have other regulatory functions in addition to regulation of iron metabolism and DT. Furthermore, our findings shed light on the transcriptional landscape of the DT encoding gene tox and present evidence for two tox antisense RNAs, which point to a new way of transcriptional regulation of toxin production.

Conclusions

This study presents extensive insights into the transcriptome of C. diphtheriae and provides a basis for future studies regarding gene characterization, transcriptional regulatory networks, and regulation of the tox gene in particular.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4481-8) contains supplementary material, which is available to authorized users.

Transcriptome Landscape of Mycobacterium smegmatis

The non-pathogenic bacterium Mycobacterium smegmatis mc²155 has been widely used as a model organism in mycobacterial research, yet a detailed study about its transcription landscape remains to be established. Here we report the transcriptome, expression profiles and transcriptional structures through growth-phase-dependent RNA sequencing (RNA-seq) as well as other related experiments. We found: (1) 2,139 transcriptional start sites (TSSs) in the genome-wide scale, of which eight samples were randomly selected and further verified by 5′-RACE; (2) 2,233 independent monocistronic or polycistronic mRNAs in the transcriptome within the operon/sub-operon structures which are classified into five groups; (3) 47.50% (1016/2139) genes were transcribed into leaderless mRNAs, with the TSSs of 41.3% (883/2139) mRNAs overlapping with the first base of the annotated start codon. Initial amino acids of MSMEG_4921 and MSMEG_6422 proteins were identified by Edman degradation, indicating the presence of distinctive widespread leaderless features in M. smegmatis mc²155. (4) 150 genes with potentially wrong structural annotation, of which 124 proposed genes have been corrected; (5) eight highly active promoters, with their activities further determined by β-galactosidase assays. These data integrated the transcriptional landscape to genome information of model organism mc²155 and lay a solid foundation for further works in Mycobacterium.

Diversity of Translation Initiation Mechanisms across Bacterial Species Is Driven by Environmental Conditions and Growth Demands

The Shine-Dalgarno (SD) sequence motif is frequently found upstream of protein coding genes and is thought to be the dominant mechanism of translation initiation used by bacteria. Experimental studies have shown that the SD sequence facilitates start codon recognition and enhances translation initiation by directly interacting with the highly conserved anti-SD sequence on the 30S ribosomal subunit. However, the proportion of SD-led genes within a genome varies across species and the factors governing this variation in translation initiation mechanisms remain largely unknown. Here, we conduct a phylogenetically informed analysis and find that species capable of rapid growth contain a higher proportion of SD-led genes throughout their genomes. We show that SD sequence utilization covaries with a suite of genomic features that are important for efficient translation initiation and elongation. In addition to these endogenous genomic factors, we further show that exogenous environmental factors may influence the evolution of translation initiation mechanisms by finding that thermophilic species contain significantly more SD-led genes than mesophiles. Our results demonstrate that variation in translation initiation mechanisms across bacterial species is predictable and is a consequence of differential life-history strategies related to maximum growth rate and environmental-specific constraints.

Leveraging genome-wide datasets to quantify the functional role of the anti-Shine-Dalgarno sequence in regulating translation efficiency

Studies dating back to the 1970s established that sequence complementarity between the anti-Shine–Dalgarno (aSD) sequence on prokaryotic ribosomes and the 5′ untranslated region of mRNAs helps to facilitate translation initiation. The optimal location of aSD sequence binding relative to the start codon, the full extents of the aSD sequence and the functional form of the relationship between aSD sequence complementarity and translation efficiency have not been fully resolved. Here, we investigate these relationships by leveraging the sequence diversity of endogenous genes and recently available genome-wide estimates of translation efficiency. We show that—after accounting for predicted mRNA structure—aSD sequence complementarity increases the translation of endogenous mRNAs by roughly 50%. Further, we observe that this relationship is nonlinear, with translation efficiency maximized for mRNAs with intermediate levels of aSD sequence complementarity. The mechanistic insights that we observe are highly robust: we find nearly identical results in multiple datasets spanning three distantly related bacteria. Further, we verify our main conclusions by re-analysing a controlled experimental dataset.

Detailed transcriptome analysis of the plant growth promoting Paenibacillus riograndensis SBR5 by using RNA-seq technology

BACKGROUND

The plant growth promoting rhizobacterium Paenibacillus riograndensis SBR5 is a promising candidate to serve as crop inoculant. Despite its potential in providing environmental and economic benefits, the species P. riograndensis is poorly characterized. Here, we performed for the first time a detailed transcriptome analysis of P. riograndensis SBR5 using RNA-seq technology.

RESULTS

RNA was isolated from P. riograndensis SBR5 cultivated under 15 different growth conditions and combined together in order to analyze an RNA pool representing a large set of expressed genes. The resultant total RNA was used to generate 2 different libraries, one enriched in 5'-ends of the primary transcripts and the other representing the whole transcriptome. Both libraries were sequenced and analyzed to identify the conserved sequences of ribosome biding sites and translation start motifs, and to elucidate operon structures present in the transcriptome of P. riograndensis. Sequence analysis of the library enriched in 5'-ends of the primary transcripts was used to identify 1082 transcription start sites (TSS) belonging to novel transcripts and allowed us to determine a promoter consensus sequence and regulatory sequences in 5' untranslated regions including riboswitches. A putative thiamine pyrophosphate dependent riboswitch upstream of the thiamine biosynthesis gene thiC was characterized by translational fusion to a fluorescent reporter gene and shown to function in P. riograndensis SBR5.

CONCLUSIONS

Our RNA-seq analysis provides insight into the P. riograndensis SBR5 transcriptome at the systems level and will be a valuable basis for differential RNA-seq analysis of this bacterium.

Comparative genomic analysis of translation initiation mechanisms for genes lacking the Shine-Dalgarno sequence in prokaryotes

In prokaryotes, translation initiation is believed to occur through an interaction between the 3΄ tail of a 16S rRNA and a corresponding Shine–Dalgarno (SD) sequence in the 5΄ untranslated region (UTR) of an mRNA. However, some genes lack SD sequences (non-SD genes), and the fraction of non-SD genes in a genome varies depending on the prokaryotic species. To elucidate non-SD translation initiation mechanisms in prokaryotes from an evolutionary perspective, we statistically examined the nucleotide frequencies around the initiation codons in non-SD genes from 260 prokaryotes (235 bacteria and 25 archaea). We identified distinct nucleotide frequency biases upstream of the initiation codon in bacteria and archaea, likely because of the presence of leaderless mRNAs lacking a 5΄ UTR. Moreover, we observed overall similarities in the nucleotide patterns between upstream and downstream regions of the initiation codon in all examined phyla. Symmetric nucleotide frequency biases might facilitate translation initiation by preventing the formation of secondary structures around the initiation codon. These features are more prominent in species’ genomes that harbor large fractions of non-SD sequences, suggesting that a reduced stability around the initiation codon is important for efficient translation initiation in prokaryotes.

Selection on start codons in prokaryotes and potential compensatory nucleotide substitutions

Reconstruction of the evolution of start codons in 36 groups of closely related bacterial and archaeal genomes reveals purifying selection affecting AUG codons. The AUG starts are replaced by GUG and especially UUG significantly less frequently than expected under the neutral expectation derived from the frequencies of the respective nucleotide triplet substitutions in non-coding regions and in 4-fold degenerate sites. Thus, AUG is the optimal start codon that is actively maintained by purifying selection. However, purifying selection on start codons is significantly weaker than the selection on the same codons in coding sequences, although the switches between the codons result in conservative amino acid substitutions. The only exception is the AUG to UUG switch that is strongly selected against among start codons. Selection on start codons is most pronounced in evolutionarily conserved, highly expressed genes. Mutation of the start codon to a sub-optimal form (GUG or UUG) tends to be compensated by mutations in the Shine-Dalgarno sequence towards a stronger translation initiation signal. Together, all these findings indicate that in prokaryotes, translation start signals are subject to weak but significant selection for maximization of initiation rate and, consequently, protein production.

Discovery of numerous novel small genes in the intergenic regions of the Escherichia coli O157:H7 Sakai genome

In the past, short protein-coding genes were often disregarded by genome annotation pipelines. Transcriptome sequencing (RNAseq) signals outside of annotated genes have usually been interpreted to indicate either ncRNA or pervasive transcription. Therefore, in addition to the transcriptome, the translatome (RIBOseq) of the enteric pathogen Escherichia coli O157:H7 strain Sakai was determined at two optimal growth conditions and a severe stress condition combining low temperature and high osmotic pressure. All intergenic open reading frames potentially encoding a protein of ≥ 30 amino acids were investigated with regard to coverage by transcription and translation signals and their translatability expressed by the ribosomal coverage value. This led to discovery of 465 unique, putative novel genes not yet annotated in this E. coli strain, which are evenly distributed over both DNA strands of the genome. For 255 of the novel genes, annotated homologs in other bacteria were found, and a machine-learning algorithm, trained on small protein-coding E. coli genes, predicted that 89% of these translated open reading frames represent bona fide genes. The remaining 210 putative novel genes without annotated homologs were compared to the 255 novel genes with homologs and to 250 short annotated genes of this E. coli strain. All three groups turned out to be similar with respect to their translatability distribution, fractions of differentially regulated genes, secondary structure composition, and the distribution of evolutionary constraint, suggesting that both novel groups represent legitimate genes. However, the machine-learning algorithm only recognized a small fraction of the 210 genes without annotated homologs. It is possible that these genes represent a novel group of genes, which have unusual features dissimilar to the genes of the machine-learning algorithm training set.

Evolutionary Origins of Two-Barrel RNA Polymerases and Site-Specific Transcription Initiation

Evolution-related multisubunit RNA polymerases (RNAPs) carry out RNA synthesis in all domains life. Although their catalytic cores and fundamental mechanisms of transcription elongation are conserved, the initiation stage of the transcription cycle differs substantially in bacteria, archaea, and eukaryotes in terms of the requirements for accessory factors and details of the molecular mechanisms. This review focuses on recent insights into the evolution of the transcription apparatus with regard to (a) the surprisingly pervasive double-Ψ β-barrel active-site configuration among different nucleic acid polymerase families, (b) the origin and phylogenetic distribution of TBP, TFB, and TFE transcription factors, and

Noncanonical Translation Initiation Comes of Age

The canonical translation initiation mechanism involves base pairing between the mRNA and 16S rRNA. However, a variety of identified mechanisms deviate from this conventional route. Beck and Janssen (J Bacteriol 199:e00091-17, 2017, https://doi.org/10.1128/JB.00091-17) have recently described another noncanonical mode of translation initiation. Here, we describe how this process differs from previously reported mechanisms, with the hope that it will foster increased awareness of the diversity of regulatory mechanisms that await discovery.

The alr-groEL1 operon in Mycobacterium tuberculosis: an interplay of multiple regulatory elements

Threonylcarbamoyladenosine is a universally conserved essential modification of tRNA that ensures translational fidelity in cellular milieu. TsaD, TsaB and TsaE are identified as tRNA-A₃₇-threonylcarbamoyl (t⁶A)-transferase enzymes that have been reconstituted in vitro, in few bacteria recently. However, transcriptional organization and regulation of these genes are not known in any of these organisms. This study describes the intricate architecture of a complex multicistronic alr-groEL1 operon, harboring essential genes, namely tsaD, tsaB, tsaE, groES, groEL1, and alr (required for cell wall synthesis), and rimI encoding an N-α- acetyltransferase in Mycobacterium tuberculosis. Using northern blotting, RT-PCR and in vivo fluorescence assays, genes alr to groEL1 were found to constitute an ~6.3 kb heptacistronic operon with multiple internal promoters and an I-shaped intrinsic hairpin-like cis-regulatory element. A strong promoter PtsaD within the coding sequence of rimI gene is identified in M. tuberculosis, in addition. The study further proposes an amendment in the known bicistronic groESL1 operon annotation by providing evidence that groESL1 is co-transcribed as sub-operon of alr-groEL1 operon. The architecture of alr-groEL1 operon, conservation of the genetic context and a mosaic transcriptional profile displayed under various stress conditions convincingly suggest the involvement of this operon in stress adaptation in M. tuberculosis.

Differentiation of ncRNAs from small mRNAs in Escherichia coli O157:H7 EDL933 (EHEC) by combined RNAseq and RIBOseq - ryhB encodes the regulatory RNA RyhB and a peptide, RyhP

Background

While NGS allows rapid global detection of transcripts, it remains difficult to distinguish ncRNAs from short mRNAs. To detect potentially translated RNAs, we developed an improved protocol for bacterial ribosomal footprinting (RIBOseq). This allowed distinguishing ncRNA from mRNA in EHEC. A high ratio of ribosomal footprints per transcript (ribosomal coverage value, RCV) is expected to indicate a translated RNA, while a low RCV should point to a non-translated RNA.

Results

Based on their low RCV, 150 novel non-translated EHEC transcripts were identified as putative ncRNAs, representing both antisense and intergenic transcripts, 74 of which had expressed homologs in E. coli MG1655. Bioinformatics analysis predicted statistically significant target regulons for 15 of the intergenic transcripts; experimental analysis revealed 4-fold or higher differential expression of 46 novel ncRNA in different growth media. Out of 329 annotated EHEC ncRNAs, 52 showed an RCV similar to protein-coding genes, of those, 16 had RIBOseq patterns matching annotated genes in other enterobacteriaceae, and 11 seem to possess a Shine-Dalgarno sequence, suggesting that such ncRNAs may encode small proteins instead of being solely non-coding. To support that the RIBOseq signals are reflecting translation, we tested the ribosomal-footprint covered ORF of ryhB and found a phenotype for the encoded peptide in iron-limiting condition.

Conclusion

Determination of the RCV is a useful approach for a rapid first-step differentiation between bacterial ncRNAs and small mRNAs. Further, many known ncRNAs may encode proteins as well.

Electronic supplementary material

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

Four translation initiation pathways employed by the leaderless mRNA in eukaryotes

mRNAs lacking 5′ untranslated regions (leaderless mRNAs) are molecular relics of an ancient translation initiation pathway. Nevertheless, they still represent a significant portion of transcriptome in some taxons, including a number of eukaryotic species. In bacteria and archaea, the leaderless mRNAs can bind non-dissociated 70 S ribosomes and initiate translation without protein initiation factors involved. Here we use the Fleeting mRNA Transfection technique (FLERT) to show that translation of a leaderless reporter mRNA is resistant to conditions when eIF2 and eIF4F, two key eukaryotic translation initiation factors, are inactivated in mammalian cells. We report an unconventional translation initiation pathway utilized by the leaderless mRNA in vitro, in addition to the previously described 80S-, eIF2-, or eIF2D-mediated modes. This mechanism is a bacterial-like eIF5B/IF2-assisted initiation that has only been reported for hepatitis C virus-like internal ribosome entry sites (IRESs). Therefore, the leaderless mRNA is able to take any of four different translation initiation pathways in eukaryotes.

The role of mRNA structure in bacterial translational regulation

The characteristics of bacterial messenger RNAs (mRNAs) that influence translation efficiency provide many convenient handles for regulation of gene expression, especially when coupled with the processes of transcription termination and mRNA degradation. An mRNA's structure, especially near the site of initiation, has profound consequences for how readily it is translated. This property allows bacterial gene expression to be altered by changes to mRNA structure induced by temperature, or interactions with a wide variety of cellular components including small molecules, other RNAs (such as sRNAs and tRNAs), and RNA-binding proteins. This review discusses the links between mRNA structure and translation efficiency, and how mRNA structure is manipulated by conditions and signals within the cell to regulate gene expression. The range of RNA regulators discussed follows a continuum from very complex tertiary structures such as riboswitch aptamers and ribosomal protein-binding sites to thermosensors and mRNA:sRNA interactions that involve only base-pairing interactions. Furthermore, the high degrees of diversity observed for both mRNA structures and the mechanisms by which inhibition of translation occur have significant consequences for understanding the evolution of bacterial translational regulation. WIREs RNA 2017, 8:e1370. doi: 10.1002/wrna.1370 For further resources related to this article, please visit the WIREs website.

The dynamic transcriptional and translational landscape of the model antibiotic producer Streptomyces coelicolor A3(2)

Individual Streptomyces species have the genetic potential to produce a diverse array of natural products of commercial, medical and veterinary interest. However, these products are often not detectable under laboratory culture conditions. To harness their full biosynthetic potential, it is important to develop a detailed understanding of the regulatory networks that orchestrate their metabolism. Here we integrate nucleotide resolution genome-scale measurements of the transcriptome and translatome of Streptomyces coelicolor, the model antibiotic-producing actinomycete. Our systematic study determines 3,570 transcription start sites and identifies 230 small RNAs and a considerable proportion (∼21%) of leaderless mRNAs; this enables deduction of genome-wide promoter architecture. Ribosome profiling reveals that the translation efficiency of secondary metabolic genes is negatively correlated with transcription and that several key antibiotic regulatory genes are translationally induced at transition growth phase. These findings might facilitate the design of new approaches to antibiotic discovery and development.

In silico analysis of 5'-UTRs highlights the prevalence of Shine-Dalgarno and leaderless-dependent mechanisms of translation initiation in bacteria and archaea, respectively

In prokaryotes, a heterogeneous set of protein translation initiation mechanisms such as Shine-Dalgarno (SD) sequence-dependent, SD sequence-independent or ribosomal protein S1 mediated and leaderless transcript-dependent exists. To estimate the distribution of coding sequences employing a particular translation initiation mechanism, a total of 107 prokaryotic genomes were analysed using in silico approaches. Analysis of 5'-untranslated regions (UTRs) of genes reveals the existence of three types of mRNAs described as transcripts with and without SD motif and leaderless transcripts. Our results indicate that although all the three types of translation initiation mechanisms are widespread among prokaryotes, the number of SD-dependent genes in bacteria is higher than that of archaea. In contrast, archaea contain a significantly higher number of leaderless genes than SD-led genes. The correlation analysis between genome size and SD-led & leaderless genes suggests that the SD-led genes are decreasing (increasing) with genome size in bacteria (archaea). However, the leaderless genes are increasing (decreasing) in bacteria (archaea) with genome size. Moreover, an analysis of the start-codon biasness confirms that among ATG, GTG and TTG codons, ATG is indeed the most preferred codon at the translation initiation site in most of the coding sequences. In leaderless genes, however, the codons GTG and TTG are also observed at the translation initiation site in some species contradicting earlier studies which suggested the usage of only ATG codon. Henceforth, the conventional mechanism of translation initiation cannot be generalized as an exclusive way of initiating the process of protein biosynthesis in prokaryotes.

A novel enrichment strategy reveals unprecedented number of novel transcription start sites at single base resolution in a model prokaryote and the gut microbiome

Background

The initiating nucleotide found at the 5’ end of primary transcripts has a distinctive triphosphorylated end that distinguishes these transcripts from all other RNA species. Recognizing this distinction is key to deconvoluting the primary transcriptome from the plethora of processed transcripts that confound analysis of the transcriptome. The currently available methods do not use targeted enrichment for the 5′end of primary transcripts, but rather attempt to deplete non-targeted RNA.

Results

We developed a method, Cappable-seq, for directly enriching for the 5' end of primary transcripts and enabling determination of transcription start sites at single base resolution. This is achieved by enzymatically modifying the 5′ triphosphorylated end of RNA with a selectable tag. We first applied Cappable-seq to E. coli, achieving up to 50 fold enrichment of primary transcripts and identifying an unprecedented 16539 transcription start sites (TSS) genome-wide at single base resolution. We also applied Cappable-seq to a mouse cecum sample and identified TSS in a microbiome.

Conclusions

Cappable-seq allows for the first time the capture of the 5′ end of primary transcripts. This enables a unique robust TSS determination in bacteria and microbiomes.  In addition to and beyond TSS determination, Cappable-seq depletes ribosomal RNA and reduces the complexity of the transcriptome to a single quantifiable tag per transcript enabling digital profiling of gene expression in any microbiome.

Electronic supplementary material

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