The Evolutionary Transition of the RNA World to Obcells to Cellular-Based Life.

The obcell hypothesis is a proposed route for the RNA world to develop into a primitive cellular one. It posits that this transition began with the emergence of the proto-ribosome which enabled RNA to colonise the external surface of lipids by the synthesis of amphipathic peptidyl-RNAs. The obcell hypothesis also posits that the emergence of a predation-based ecosystem provided a selection mechanism for continued sophistication amongst early life forms. Here, I argue for this hypothesis owing to its significant explanatory power; it offers a rationale why a ribosome which initially was capable only of producing short non-coded peptides was advantageous and it forgoes issues related to maintaining a replicating RNA inside a lipid enclosure. I develop this model by proposing that the evolutionary selection for improved membrane anchors resulted in the emergence of primitive membrane pores which enabled obcells to gradually evolve into a cellular morphology. Moreover, I introduce a model of obcell production which advances that tRNAs developed from primers of the RNA world.

AMD1 mRNA employs ribosome stalling as a mechanism for molecular memory formation.

In addition to acting as template for protein synthesis, messenger RNA (mRNA) often contains sensory sequence elements that regulate this process. Here we report a new mechanism that limits the number of complete protein molecules that can be synthesized from a single mRNA molecule of the human AMD1 gene encoding adenosylmethionine decarboxylase 1 (AdoMetDC). A small proportion of ribosomes translating AMD1 mRNA stochastically read through the stop codon of the main coding region. These readthrough ribosomes then stall close to the next in-frame stop codon, eventually forming a ribosome queue, the length of which is proportional to the number of AdoMetDC molecules that were synthesized from the same AMD1 mRNA. Once the entire spacer region between the two stop codons is filled with queueing ribosomes, the queue impinges upon the main AMD1 coding region halting its translation. Phylogenetic analysis suggests that this mechanism is highly conserved in vertebrates and existed in their common ancestor. We propose that this mechanism is used to count and limit the number of protein molecules that can be synthesized from a single mRNA template. It could serve to safeguard from dysregulated translation that may occur owing to errors in transcription or mRNA damage.

Evaluating data integrity in ribosome footprinting datasets through modelled polysome profiles.

The assessment of transcriptome-wide ribosome binding to mRNAs is useful for studying the dynamic regulation of protein synthesis. Two methods frequently applied in eukaryotic cells that operate at different levels of resolution are polysome profiling, which reveals the distribution of ribosome loads across the transcriptome, and ribosome footprinting (also termed ribosome profiling or Ribo-Seq), which when combined with appropriate data on mRNA expression can reveal ribosome densities on individual transcripts. In this study we develop methods for relating the information content of these two methods to one another, by reconstructing theoretical polysome profiles from ribosome footprinting data. Our results validate both approaches as experimental tools. Although we show that both methods can yield highly consistent data, some published ribosome footprinting datasets give rise to reconstructed polysome profiles with non-physiological features. We trace these aberrant features to inconsistencies in RNA and Ribo-Seq data when compared to datasets yielding physiological polysome profiles, thereby demonstrating that modelled polysomes are useful for assessing global dataset properties such as its quality in a simple, visual approach. Aside from using polysome profile reconstructions on published datasets, we propose that this also provides a useful tool for validating new ribosome footprinting datasets in early stages of analyses.

Translation initiation downstream from annotated start codons in human mRNAs coevolves with the Kozak context.

Eukaryotic translation initiation involves preinitiation ribosomal complex 5'-to-3' directional probing of mRNA for codons suitable for starting protein synthesis. The recognition of codons as starts depends on the codon identity and on its immediate nucleotide context known as Kozak context. When the context is weak (i.e., nonoptimal), leaky scanning takes place during which a fraction of ribosomes continues the mRNA probing. We explored the relationship between the context of AUG codons annotated as starts of protein-coding sequences and the next AUG codon occurrence. We found that AUG codons downstream from weak starts occur in the same frame more frequently than downstream from strong starts. We suggest that evolutionary selection on in-frame AUGs downstream from weak start codons is driven by the advantage of the reduction of wasteful out-of-frame product synthesis and also by the advantage of producing multiple proteoforms from certain mRNAs. We confirmed translation initiation downstream from weak start codons using ribosome profiling data. We also tested translation of alternative start codons in 10 specific human genes using reporter constructs. In all tested cases, initiation at downstream start codons was more productive than at the annotated ones. In most cases, optimization of Kozak context did not completely abolish downstream initiation, and in the specific example of CMPK1 mRNA, the optimized start remained unproductive. Collectively, our work reveals previously uncharacterized forces shaping the evolution of protein-coding genes and points to the plurality of translation initiation and the existence of sequence features influencing start codon selection, other than Kozak context.

Streptomyces rare codon UUA: from features associated with 2 adpA related locations to candidate phage regulatory translational bypassing.

In Streptomyces species, the cell cycle involves a switch from an early and vegetative state to a later phase where secondary products including antibiotics are synthesized, aerial hyphae form and sporulation occurs. AdpA, which has two domains, activates the expression of numerous genes involved in the switch from the vegetative growth phase. The adpA mRNA of many Streptomyces species has a UUA codon in a linker region between 5' sequence encoding one domain and 3' sequence encoding its other and C-terminal domain. UUA codons are exceptionally rare in Streptomyces, and its functional cognate tRNA is not present in a fully modified and acylated form, in the early and vegetative phase of the cell cycle though it is aminoacylated later. Here, we report candidate recoding signals that may influence decoding of the linker region UUA. Additionally, a short ORF 5' of the main ORF has been identified with a GUG at, or near, its 5' end and an in-frame UUA near its 3' end. The latter is commonly 5 nucleotides 5' of the main ORF start. Ribosome profiling data show translation of that 5' region. Ten years ago, UUA-mediated translational bypassing was proposed as a sensor by a Streptomyces phage of its host's cell cycle stage and an effector of its lytic/lysogeny switch. We provide the first experimental evidence supportive of this proposal.

Unusually efficient CUG initiation of an overlapping reading frame in POLG mRNA yields novel protein POLGARF.

While near-cognate codons are frequently used for translation initiation in eukaryotes, their efficiencies are usually low (<10% compared to an AUG in optimal context). Here, we describe a rare case of highly efficient near-cognate initiation. A CUG triplet located in the 5' leader of POLG messenger RNA (mRNA) initiates almost as efficiently (∼60 to 70%) as an AUG in optimal context. This CUG directs translation of a conserved 260-triplet-long overlapping open reading frame (ORF), which we call POLGARF (POLG Alternative Reading Frame). Translation of a short upstream ORF 5' of this CUG governs the ratio between POLG (the catalytic subunit of mitochondrial DNA polymerase) and POLGARF synthesized from a single POLG mRNA. Functional investigation of POLGARF suggests a role in extracellular signaling. While unprocessed POLGARF localizes to the nucleoli together with its interacting partner C1QBP, serum stimulation results in rapid cleavage and secretion of a POLGARF C-terminal fragment. Phylogenetic analysis shows that POLGARF evolved ∼160 million y ago due to a mammalian-wide interspersed repeat (MIR) transposition into the 5' leader sequence of the mammalian POLG gene, which became fixed in placental mammals. This discovery of POLGARF unveils a previously undescribed mechanism of de novo protein-coding gene evolution.

Trips-Viz: a transcriptome browser for exploring Ribo-Seq data.

Ribosome profiling (Ribo-Seq) is a technique that allows for the isolation and sequencing of mRNA fragments protected from nuclease digestion by actively translating ribosomes. Mapping these ribosome footprints to a genome or transcriptome generates quantitative information on translated regions. To provide access to publicly available ribosome profiling data in the context of transcriptomes we developed Trips-Viz (transcriptome-wide information on protein synthesis-visualized). Trips-Viz provides a large range of graphical tools for exploring global properties of translatomes and of individual transcripts. It enables analysis of aligned footprints to evaluate datasets quality, differential gene expression detection, visual identification of upstream ORFs and alternative proteoforms. Trips-Viz is available at https://trips.ucc.ie.

GWIPS-viz: 2018 update.

The GWIPS-viz browser (http://gwips.ucc.ie/) is an on-line genome browser which is tailored for exploring ribosome profiling (Ribo-seq) data. Since its publication in 2014, GWIPS-viz provides Ribo-seq data for an additional 14 genomes bringing the current total to 23. The integration of new Ribo-seq data has been automated thereby increasing the number of available tracks to 1792, a 10-fold increase in the last three years. The increase is particularly substantial for data derived from human sources. Following user requests, we added the functionality to download these tracks in bigWig format. We also incorporated new types of data (e.g. TCP-seq) as well as auxiliary tracks from other sources that help with the interpretation of Ribo-seq data. Improvements in the visualization of the data have been carried out particularly for bacterial genomes where the Ribo-seq data are now shown in a strand specific manner. For higher eukaryotic datasets, we provide characteristics of individual datasets using the RUST program which includes the triplet periodicity, sequencing biases and relative inferred A-site dwell times. This information can be used for assessing the quality of Ribo-seq datasets. To improve the power of the signal, we aggregate Ribo-seq data from several studies into Global aggregate tracks for each genome.

Insights into the mechanisms of eukaryotic translation gained with ribosome profiling.

The development of Ribosome Profiling (RiboSeq) has revolutionized functional genomics. RiboSeq is based on capturing and sequencing of the mRNA fragments enclosed within the translating ribosome and it thereby provides a 'snapshot' of ribosome positions at the transcriptome wide level. Although the method is predominantly used for analysis of differential gene expression and discovery of novel translated ORFs, the RiboSeq data can also be a rich source of information about molecular mechanisms of polypeptide synthesis and translational control. This review will focus on how recent findings made with RiboSeq have revealed important details of the molecular mechanisms of translation in eukaryotes. These include mRNA translation sensitivity to drugs affecting translation initiation and elongation, the roles of upstream ORFs in response to stress, the dynamics of elongation and termination as well as details of intrinsic ribosome behavior on the mRNA after translation termination. As the RiboSeq method is still at a relatively early stage we will also discuss the implications of RiboSeq artifacts on data interpretation.

RiboGalaxy: A browser based platform for the alignment, analysis and visualization of ribosome profiling data.

Ribosome profiling (ribo-seq) is a technique that uses high-throughput sequencing to reveal the exact locations and densities of translating ribosomes at the entire transcriptome level. The technique has become very popular since its inception in 2009. Yet experimentalists who generate ribo-seq data often have to rely on bioinformaticians to process and analyze their data. We present RiboGalaxy ( http://ribogalaxy.ucc.ie ), a freely available Galaxy-based web server for processing and analyzing ribosome profiling data with the visualization functionality provided by GWIPS-viz ( http://gwips.ucc.ie ). RiboGalaxy offers researchers a suite of tools specifically tailored for processing ribo-seq and corresponding mRNA-seq data. Researchers can take advantage of the published workflows which reduce the multi-step alignment process to a minimum of inputs from the user. Users can then explore their own aligned data as custom tracks in GWIPS-viz and compare their ribosome profiles to existing ribo-seq tracks from published studies. In addition, users can assess the quality of their ribo-seq data, determine the strength of the triplet periodicity signal, generate meta-gene ribosome profiles as well as analyze the relative impact of mRNA sequence features on local read density. RiboGalaxy is accompanied by extensive documentation and tips for helping users. In addition we provide a forum ( http://gwips.ucc.ie/Forum ) where we encourage users to post their questions and feedback to improve the overall RiboGalaxy service.

GWIPS-viz: development of a ribo-seq genome browser.

We describe the development of GWIPS-viz (http://gwips.ucc.ie), an online genome browser for viewing ribosome profiling data. Ribosome profiling (ribo-seq) is a recently developed technique that provides genome-wide information on protein synthesis (GWIPS) in vivo. It is based on the deep sequencing of ribosome-protected messenger RNA (mRNA) fragments, which allows the ribosome density along all mRNA transcripts present in the cell to be quantified. Since its inception, ribo-seq has been carried out in a number of eukaryotic and prokaryotic organisms. Owing to the increasing interest in ribo-seq, there is a pertinent demand for a dedicated ribo-seq genome browser. GWIPS-viz is based on The University of California Santa Cruz (UCSC) Genome Browser. Ribo-seq tracks, coupled with mRNA-seq tracks, are currently available for several genomes: human, mouse, zebrafish, nematode, yeast, bacteria (Escherichia coli K12, Bacillus subtilis), human cytomegalovirus and bacteriophage lambda. Our objective is to continue incorporating published ribo-seq data sets so that the wider community can readily view ribosome profiling information from multiple studies without the need to carry out computational processing.

The effect of heterogeneous Transcription Start Sites (TSS) on the translatome: implications for the mammalian cellular phenotype.

BACKGROUND: The genetic program, as manifested as the cellular phenotype, is in large part dictated by the cell's protein composition. Since characterisation of the proteome remains technically laborious it is attractive to define the genetic expression profile using the transcriptome. However, the transcriptional landscape is complex and it is unclear as to what extent it reflects the ribosome associated mRNA population (the translatome). This is particularly pertinent for genes using multiple transcriptional start sites (TSS) generating mRNAs with heterogeneous 5' transcript leaders (5'TL). Furthermore, the relative abundance of the TSS gene variants is frequently cell-type specific. Indeed, promoter switches have been reported in pathologies such as cancer. The consequences of this 5'TL heterogeneity within the transcriptome for the translatome remain unresolved. This is not a moot point because the 5'TL plays a key role in regulating mRNA recruitment onto polysomes. RESULTS: In this article, we have characterised both the transcriptome and translatome of the MCF7 (tumoural) and MCF10A (non-tumoural) cell lines. We identified ~550 genes exhibiting differential translation efficiency (TE). In itself, this is maybe not surprising. However, by focusing on genes exhibiting TSS heterogeneity we observed distinct differential promoter usage patterns in both the transcriptome and translatome. Only a minor fraction of these genes belonged to those exhibiting differential TE. Nonetheless, reporter assays demonstrated that the TSS variants impacted on the translational readout both quantitatively (the overall amount of protein expressed) and qualitatively (the nature of the proteins expressed). CONCLUSIONS: The results point to considerable and distinct cell-specific 5'TL heterogeneity within both the transcriptome and translatome of the two cell lines analysed. This observation is in-line with the ribosome filter hypothesis which posits that the ribosomal machine can selectively filter information from within the transcriptome. As such it cautions against the simple extrapolation transcriptome → proteome. Furthermore, polysomal occupancy of specific gene 5'TL variants may also serve as novel disease biomarkers.

rRNA:mRNA pairing alters the length and the symmetry of mRNA-protected fragments in ribosome profiling experiments.

MOTIVATION: Ribosome profiling is a new technique that allows monitoring locations of translating ribosomes on mRNA at a whole transcriptome level. A recent ribosome profiling study demonstrated that internal Shine-Dalgarno (SD) sequences have a major global effect on translation rates in bacteria: ribosomes pause at SD sites in mRNA. Therefore, it is important to understand how SD sites effect mRNA movement through the ribosome and generation of ribosome footprints. RESULTS: Here, we provide evidence that in addition to pausing effect, internal SD sequences induce a caterpillar-like movement of mRNA through the ribosome cavity. Once an SD site binds to the ribosome, it remains attached to it while the ribosome decodes a few subsequent codons. This leads to asymmetric progressive elongation of ribosome footprints at the 3'-end. It is likely that internal SD sequences induce a pause not on a single, but on several adjacent codons. This finding is important for our understanding of mRNA movement through the ribosome and also should facilitate interpretation of ribosome profiling data.

Ribosome profiling reveals downregulation of UMP biosynthesis as the major early response to phage infection.

Bacteria have evolved diverse defense mechanisms to counter bacteriophage attacks. Genetic programs activated upon infection characterize phage-host molecular interactions and ultimately determine the outcome of the infection. In this study, we applied ribosome profiling to monitor protein synthesis during the early stages of sk1 bacteriophage infection in Lactococcus cremoris. Our analysis revealed major changes in gene expression within 5 minutes of sk1 infection. Notably, we observed a specific and severe downregulation of several pyr operons which encode enzymes required for uridine monophosphate biosynthesis. Consistent with previous findings, this is likely an attempt of the host to starve the phage of nucleotides it requires for propagation. We also observed a gene expression response that we expect to benefit the phage. This included the upregulation of 40 ribosome proteins that likely increased the host's translational capacity, concurrent with a downregulation of genes that promote translational fidelity (lepA and raiA). In addition to the characterization of host-phage gene expression responses, the obtained ribosome profiling data enabled us to identify two putative recoding events as well as dozens of loci currently annotated as pseudogenes that are actively translated. Furthermore, our study elucidated alterations in the dynamics of the translation process, as indicated by time-dependent changes in the metagene profile, suggesting global shifts in translation rates upon infection. Additionally, we observed consistent modifications in the ribosome profiles of individual genes, which were apparent as early as 2 minutes post-infection. The study emphasizes our ability to capture rapid alterations of gene expression during phage infection through ribosome profiling. IMPORTANCE: The ribosome profiling technology has provided invaluable insights for understanding cellular translation and eukaryotic viral infections. However, its potential for investigating host-phage interactions remains largely untapped. Here, we applied ribosome profiling to Lactococcus cremoris cultures infected with sk1, a major infectious agent in dairy fermentation processes. This revealed a profound downregulation of genes involved in pyrimidine nucleotide synthesis at an early stage of phage infection, suggesting an anti-phage program aimed at restricting nucleotide availability and, consequently, phage propagation. This is consistent with recent findings and contributes to our growing appreciation for the role of nucleotide limitation as an anti-viral strategy. In addition to capturing rapid alterations in gene expression levels, we identified translation occurring outside annotated regions, as well as signatures of non-standard translation mechanisms. The gene profiles revealed specific changes in ribosomal densities upon infection, reflecting alterations in the dynamics of the translation process.

In Vivo RNAi-Mediated eIF3m Knockdown Affects Ribosome Biogenesis and Transcription but Has Limited Impact on mRNA-Specific Translation.

Translation is an essential biological process, and dysregulation is associated with a range of diseases including ribosomopathies, diabetes, and cancer. Here, we examine translation dysregulation in vivo using RNAi to knock down the m-subunit of the translation initiation factor eIF3 in the mouse liver. Transcriptome sequencing, ribosome profiling, whole proteome, and phosphoproteome analyses show that eIF3m deficiency leads to the transcriptional response and changes in cellular translation that yield few detectable differences in the translation of particular mRNAs. The transcriptional response fell into two main categories: ribosome biogenesis (increased transcription of ribosomal proteins) and cell metabolism (alterations in lipid, amino acid, nucleic acid, and drug metabolism). Analysis of ribosome biogenesis reveals inhibition of rRNA processing, highlighting decoupling of rRNA synthesis and ribosomal protein gene transcription in response to eIF3m knockdown. Interestingly, a similar reduction in eIF3m protein levels is associated with induction of the mTOR pathway in vitro but not in vivo. Overall, this work highlights the utility of a RNAi-based in vivo approach for studying the regulation of mammalian translation in vivo.

Rocaglates Induce Gain-of-Function Alterations to eIF4A and eIF4F.

Rocaglates are a diverse family of biologically active molecules that have gained tremendous interest in recent years due to their promising activities in pre-clinical cancer studies. As a result, this family of compounds has been significantly expanded through the development of efficient synthetic schemes. However, it is unknown whether all of the members of the rocaglate family act through similar mechanisms of action. Here, we present a comprehensive study comparing the biological activities of >200 rocaglates to better understand how the presence of different chemical entities influences their biological activities. Through this, we find that most rocaglates preferentially repress the translation of mRNAs containing purine-rich 5' leaders, but certain rocaglates lack this bias in translation repression. We also uncover an aspect of rocaglate mechanism of action in which the pool of translationally active eIF4F is diminished due to the sequestration of the complex onto RNA.

Polysomes Bypass a 50-Nucleotide Coding Gap Less Efficiently Than Monosomes Due to Attenuation of a 5' mRNA Stem-Loop and Enhanced Drop-off.

Efficient translational bypassing of a 50-nt non-coding gap in a phage T4 topoisomerase subunit gene (gp60) requires several recoding signals. Here we investigate the function of the mRNA stem-loop 5' of the take-off codon, as well as the importance of ribosome loading density on the mRNA for efficient bypassing. We show that polysomes are less efficient at mediating bypassing than monosomes, both in vitro and in vivo, due to their preventing formation of a stem-loop 5' of the take-off codon and allowing greater peptidyl-tRNA drop off. A ribosome profiling analysis of phage T4-infected Escherichia coli yielded protected mRNA fragments within the normal size range derived from ribosomes stalled at the take-off codon. However, ribosomes at this position also yielded some 53-nucleotide fragments, 16 longer. These were due to protection of the nucleotides that form the 5' stem-loop. NMR shows that the 5' stem-loop is highly dynamic. The importance of different nucleotides in the 5' stem-loop is revealed by mutagenesis studies. These data highlight the significance of the 5' stem-loop for the 50-nt bypassing and further enhance appreciation of relevance of the extent of ribosome loading for recoding.

Itm2a silencing rescues lamin A mediated inhibition of 3T3-L1 adipocyte differentiation.

Dysregulation of adipose tissue metabolism is associated with multiple metabolic disorders. One such disease, known as Dunnigan-type familial partial lipodystrophy (FPLD2) is characterized by defective fat metabolism and storage. FPLD2 is caused by a specific subset of mutations in the LMNA gene. The mechanisms by which LMNA mutations lead to the adipose specific FPLD2 phenotype have yet to be determined in detail. We used RNA-Seq analysis to assess the effects of wild-type (WT) and mutant (R482W) lamin A on the expression profile of differentiating 3T3-L1 mouse preadipocytes and identified Itm2a as a gene that was upregulated at 36 h post differentiation induction in these cells. In this study we identify Itm2a as a novel modulator of adipogenesis and show that endogenous Itm2a expression is transiently downregulated during induction of 3T3-L1 differentiation. Itm2a overexpression was seen to moderately inhibit differentiation of 3T3-L1 preadipocytes while shRNA mediated knockdown of Itm2a significantly enhanced 3T3-L1 differentiation. Investigation of PPARγ levels indicate that this enhanced adipogenesis is mediated through the stabilization of the PPARγ protein at specific time points during differentiation. Finally, we demonstrate that Itm2a knockdown is sufficient to rescue the inhibitory effects of lamin A WT and R482W mutant overexpression on 3T3-L1 differentiation. This suggests that targeting of Itm2a or its related pathways, including autophagy, may have potential as a therapy for FPLD2.

Comparative survey of the relative impact of mRNA features on local ribosome profiling read density.

Ribosome profiling (Ribo-seq), a promising technology for exploring ribosome decoding rates, is characterized by the presence of infrequent high peaks in ribosome footprint density and by long alignment gaps. Here, to reduce the impact of data heterogeneity we introduce a simple normalization method, Ribo-seq Unit Step Transformation (RUST). RUST is robust and outperforms other normalization techniques in the presence of heterogeneous noise. We illustrate how RUST can be used for identifying mRNA sequence features that affect ribosome footprint densities globally. We show that a few parameters extracted with RUST are sufficient for predicting experimental densities with high accuracy. Importantly the application of RUST to 30 publicly available Ribo-seq data sets revealed a substantial variation in sequence determinants of ribosome footprint frequencies, questioning the reliability of Ribo-seq as an accurate representation of local ribosome densities without prior quality control. This emphasizes our incomplete understanding of how protocol parameters affect ribosome footprint densities.

Oxygen and glucose deprivation induces widespread alterations in mRNA translation within 20 minutes.

BACKGROUND: Oxygen and glucose metabolism play pivotal roles in many (patho)physiological conditions. In particular, oxygen and glucose deprivation (OGD) during ischemia and stroke results in extensive tissue injury and cell death. RESULTS: Using time-resolved ribosome profiling, we assess gene expression levels in a neural cell line, PC12, during the first hour of OGD. The most substantial alterations are seen to occur within the first 20 minutes of OGD. While transcription of only 100 genes is significantly altered during one hour of OGD, the translation response affects approximately 3,000 genes. This response involves reprogramming of initiation and elongation rates, as well as the stringency of start codon recognition. Genes involved in oxidative phosphorylation are most affected. Detailed analysis of ribosome profiles reveals salient alterations of ribosome densities on individual mRNAs. The mRNA-specific alterations include increased translation of upstream open reading frames, site-specific ribosome pauses, and production of alternative protein isoforms with amino-terminal extensions. Detailed analysis of ribosomal profiles also reveals six mRNAs with translated ORFs occurring downstream of annotated coding regions and two examples of dual coding mRNAs, where two protein products are translated from the same long segment of mRNA, but in two different frames. CONCLUSIONS: These findings uncover novel regulatory mechanisms of translational response to OGD in mammalian cells that are different from the classical pathways such as hypoxia inducible factor (HIF) signaling, while also revealing sophisticated organization of protein coding information in certain genes.