Ribosome decision graphs for the representation of eukaryotic RNA translation complexity.

The application of ribosome profiling has revealed an unexpected abundance of translation in addition to that responsible for the synthesis of previously annotated protein-coding regions. Multiple short sequences have been found to be translated within single RNA molecules, within both annotated protein-coding and noncoding regions. The biological significance of this translation is a matter of intensive investigation. However, current schematic or annotation-based representations of mRNA translation generally do not account for the apparent multitude of translated regions within the same molecules. They also do not take into account the stochasticity of the process that allows alternative translations of the same RNA molecules by different ribosomes. There is a need for formal representations of mRNA complexity that would enable the analysis of quantitative information on translation and more accurate models for predicting the phenotypic effects of genetic variants affecting translation. To address this, we developed a conceptually novel abstraction that we term ribosome decision graphs (RDGs). RDGs represent translation as multiple ribosome paths through untranslated and translated mRNA segments. We termed the latter "translons." Nondeterministic events, such as initiation, reinitiation, selenocysteine insertion, or ribosomal frameshifting, are then represented as branching points. This representation allows for an adequate representation of eukaryotic translation complexity and focuses on locations critical for translation regulation. We show how RDGs can be used for depicting translated regions and for analyzing genetic variation and quantitative genome-wide data on translation for characterization of regulatory modulators of translation.

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.

Impact of eIF2α phosphorylation on the translational landscape of mouse embryonic stem cells.

The integrated stress response (ISR) is critical for cell survival under stress. In response to diverse environmental cues, eIF2α becomes phosphorylated, engendering a dramatic change in mRNA translation. The activation of ISR plays a pivotal role in the early embryogenesis, but the eIF2-dependent translational landscape in pluripotent embryonic stem cells (ESCs) is largely unexplored. We employ a multi-omics approach consisting of ribosome profiling, proteomics, and metabolomics in wild-type (eIF2α+/+) and phosphorylation-deficient mutant eIF2α (eIF2αA/A) mouse ESCs (mESCs) to investigate phosphorylated (p)-eIF2α-dependent translational control of naive pluripotency. We show a transient increase in p-eIF2α in the naive epiblast layer of E4.5 embryos. Absence of eIF2α phosphorylation engenders an exit from naive pluripotency following 2i (two chemical inhibitors of MEK1/2 and GSK3α/ß) withdrawal. p-eIF2α controls translation of mRNAs encoding proteins that govern pluripotency, chromatin organization, and glutathione synthesis. Thus, p-eIF2α acts as a key regulator of the naive pluripotency gene regulatory network.

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.

Ribosome Decision Graphs for the Representation of Eukaryotic RNA Translation Complexity.

The application of ribosome profiling has revealed an unexpected abundance of translation in addition to that responsible for the synthesis of previously annotated protein-coding regions. Multiple short sequences have been found to be translated within single RNA molecules, both within annotated protein-coding and non-coding regions. The biological significance of this translation is a matter of intensive investigation. However, current schematic or annotation-based representations of mRNA translation generally do not account for the apparent multitude of translated regions within the same molecules. They also do not take into account the stochasticity of the process that allows alternative translations of the same RNA molecules by different ribosomes. There is a need for formal representations of mRNA complexity that would enable the analysis of quantitative information on translation and more accurate models for predicting the phenotypic effects of genetic variants affecting translation. To address this, we developed a conceptually novel abstraction that we term Ribosome Decision Graphs (RDGs). RDGs represent translation as multiple ribosome paths through untranslated and translated mRNA segments. We termed the later 'translons'. Non-deterministic events, such as initiation, re-initiation, selenocysteine insertion or ribosomal frameshifting are then represented as branching points. This representation allows for an adequate representation of eukaryotic translation complexity and focuses on locations critical for translation regulation. We show how RDGs can be used for depicting translated regions, analysis of genetic variation and quantitative genome-wide data on translation for characterisation of regulatory modulators of translation.

RiboGalaxy: A Galaxy-based Web Platform for Ribosome Profiling Data Processing - 2023 Update.

Ribosome profiling (Ribo-Seq) captures a "snapshot" of ribosomes' locations at the entire transcriptome of a cell at sub-codon resolution providing insights into gene expression and enabling the discovery of novel translated regions. RiboGalaxy (https://ribogalaxy.genomicsdatascience.ie/), a Galaxy-based platform for processing Ribo-Seq data is a RiboSeq.Org (https://riboseq.org/) resource. RiboSeq.Org is an online gateway to a set of integrated tools for the processing and analysis of Ribo-Seq data. In this RiboGalaxy update we introduce changes to both the tools available on RiboGalaxy and to how the resource is managed on the backend. For example, in order to improve interoperability between Riboseq.Org resources, we added tools that link RiboGalaxy outputs with Trips-Viz and GWIPS-viz browsers for downstream analysis and visualisation. RiboGalaxy's backend now utilises Ansible configuration management which enhances its stability and jobs are executed within Singularity containers and are managed by Slurm, strengthening reproducibility and performance respectively.

Nontriplet feature of genetic code in Euplotes ciliates is a result of neutral evolution.

The triplet nature of the genetic code is considered a universal feature of known organisms. However, frequent stop codons at internal mRNA positions in Euplotes ciliates ultimately specify ribosomal frameshifting by one or two nucleotides depending on the context, thus posing a nontriplet feature of the genetic code of these organisms. Here, we sequenced transcriptomes of eight Euplotes species and assessed evolutionary patterns arising at frameshift sites. We show that frameshift sites are currently accumulating more rapidly by genetic drift than they are removed by weak selection. The time needed to reach the mutational equilibrium is several times longer than the age of Euplotes and is expected to occur after a several-fold increase in the frequency of frameshift sites. This suggests that Euplotes are at an early stage of the spread of frameshifting in expression of their genome. In addition, we find the net fitness burden of frameshift sites to be noncritical for the survival of Euplotes. Our results suggest that fundamental genome-wide changes such as a violation of the triplet character of genetic code can be introduced and maintained solely by neutral evolution.

Monitoring translation in all reading frames downstream of weak stop codons provides mechanistic insights into the impact of nucleotide and cellular contexts.

A stop codon entering the ribosome A-site is normally decoded by release factors that induce release of the polypeptide. Certain factors influence the efficiency of the termination which is in competition with elongation in either the same (readthrough) or an alternative (frameshifting) reading frame. To gain insight into the competition between these processes, we monitored translation in parallel from all three reading frames downstream of stop codons while changing the nucleotide context of termination sites or altering cellular conditions (polyamine levels). We found that P-site codon identity can have a major impact on the termination efficiency of the OPRL1 stop signal, whereas for the OAZ1 ORF1 stop signal, the P-site codon mainly influences the reading frame of non-terminating ribosomes. Changes to polyamine levels predominantly influence the termination efficiency of the OAZ1 ORF1 stop signal. In contrast, increasing polyamine levels stimulate readthrough of the OPRL1 stop signal by enhancing near-cognate decoding rather than by decreasing termination efficiency. Thus, by monitoring the four competing processes occurring at stop codons we were able to determine which is the most significantly affected upon perturbation. This approach may be useful for the interrogation of other recoding phenomena where alternative decoding processes compete with standard decoding.

Thousands of human non-AUG extended proteoforms lack evidence of evolutionary selection among mammals.

The synthesis of most proteins begins at AUG codons, yet a small number of non-AUG initiated proteoforms are also known. Here we analyse a large number of publicly available Ribo-seq datasets to identify novel, previously uncharacterised non-AUG proteoforms using Trips-Viz implementation of a novel algorithm for detecting translated ORFs. In parallel we analyse genomic alignment of 120 mammals to identify evidence of protein coding evolution in sequences encoding potential extensions. Unexpectedly we find that the number of non-AUG proteoforms identified with ribosome profiling data greatly exceeds those with strong phylogenetic support suggesting their recent evolution. Our study argues that the protein coding potential of human genome greatly exceeds that detectable through comparative genomics and exposes the existence of multiple proteins encoded by the same genomic loci.

Evolution of naturally arising SARS-CoV-2 defective interfering particles.

Defective interfering (DI) particles arise during virus propagation, are conditional on parental virus for replication and packaging, and interfere with viral expansion. There is much interest in developing DIs as anti-viral agents. Here we characterize DI particles that arose following serial passaging of SARS-CoV-2 at high multiplicity of infection. The prominent DIs identified have lost ~84% of the SARS-CoV-2 genome and are capable of attenuating parental viral titers. Synthetic variants of the DI genomes also interfere with infection and can be used as conditional, gene delivery vehicles. In addition, the DI genomes encode an Nsp1-10 fusion protein capable of attenuating viral replication. These results identify naturally selected defective viral genomes that emerged and stably propagated in the presence of parental virus.

Lack of evidence for ribosomal frameshifting in ATP7B mRNA decoding.

The research article describing the discovery of ribosomal frameshifting in the bacterial CopA gene also reported the occurrence of frameshifting in the expression of the human ortholog ATP7B based on assays using dual luciferase reporters. An examination of the publicly available ribosome profiling data and the phylogenetic analysis of the proposed frameshifting site cast doubt on the validity of this claim and prompted us to reexamine the evidence. We observed similar apparent frameshifting efficiencies as the original authors using the same type of vector that synthesizes both luciferases as a single polyprotein. However, we noticed anomalously low absolute luciferase activities from the N-terminal reporter that suggests interference of reporter activity or levels by the ATP7B test cassette. When we tested the same proposed ATP7B frameshifting cassette in a more recently developed reporter system in which the reporters are released without being included in a polyprotein, no frameshifting was detected above background levels.

The evolution and role of the periplasmic asparaginase Asp3 in yeast.

The study of nitrogen assimilation in yeast is of interest from genetic, evolutionary, and biotechnological perspectives. Over the course of evolution, yeasts have developed sophisticated control mechanisms to regulate nitrogen metabolism, with domesticated lineages sometimes displaying particular specialisation. The focus of this study was on assimilation of asparagine, which is a significant nutritional source for some alcoholic fermentations. We were particularly interested in ASP3, which encodes a periplasmic asparaginase and that was proposed to have been acquired relatively recently in S. cerevisiae by horizontal gene transfer. We examined 1680 S. cerevisiae genome assemblies to evaluate the distribution and evolutionary trajectory of ASP3. Our findings suggest an alternative hypothesis that ASP3 is an ancient Saccharomyces gene that has generally been lost over the course of evolution but has been retained in certain fermentative environments. As asparagine is the major nitrogen source in apple juice, we explored whether the presence of ASP3 would confer a growth advantage. Interestingly, we found that although ASP3 enhances growth when asparagine is the sole nitrogen source, the same effect is not seen in apple juice. These data indicate that growth in pure culture may not reflect the original selective environment for ASP3+ strains and highlight the role that complex regulation may play in optimising nitrogen assimilation in yeasts.

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.

Ghrelin rapidly elevates protein synthesis in vitro by employing the rpS6K-eEF2K-eEF2 signalling axis.

Activated ghrelin receptor GHS-R1α triggers cell signalling pathways that modulate energy homeostasis and biosynthetic processes. However, the effects of ghrelin on mRNA translation are unknown. Using various reporter assays, here we demonstrate a rapid elevation of protein synthesis in cells within 15-30 min upon stimulation of GHS-R1α by ghrelin. We further show that ghrelin-induced activation of translation is mediated, at least in part, through the de-phosphorylation (de-suppression) of elongation factor 2 (eEF2). The levels of eEF2 phosphorylation at Thr56 decrease due to the reduced activity of eEF2 kinase, which is inhibited via Ser366 phosphorylation by rpS6 kinases. Being stress-susceptible, the ghrelin-mediated decrease in eEF2 phosphorylation can be abolished by glucose deprivation and mitochondrial uncoupling. We believe that the observed burst of translation benefits rapid restocking of neuropeptides, which are released upon GHS-R1α activation, and represents the most time- and energy-efficient way of prompt recharging the orexigenic neuronal circuitry.

Development of a ribosome profiling protocol to study translation in Kluyveromyces marxianus.

Kluyveromyces marxianus is an interesting and important yeast because of particular traits such as thermotolerance and rapid growth, and for applications in food and industrial biotechnology. For both understanding its biology and developing bioprocesses, it is important to understand how K. marxianus responds and adapts to changing environments. For this, a full suite of omics tools to measure and compare global patterns of gene expression and protein synthesis is needed. We report here the development of a ribosome profiling method for K. marxianus, which allows codon resolution of translation on a genome-wide scale by deep sequencing of ribosome locations on mRNAs. To aid in the analysis and sharing of ribosome profiling data, we added the K. marxianus genome as well as transcriptome and ribosome profiling data to the publicly accessible GWIPS-viz and Trips-Viz browsers. Users are able to upload custom ribosome profiling and RNA-Seq data to both browsers, therefore allowing easy analysis and sharing of data. We also provide a set of step-by-step protocols for the experimental and bioinformatic methods that we developed.

Non-AUG translation initiation in mammals.

Recent proteogenomic studies revealed extensive translation outside of annotated protein coding regions, such as non-coding RNAs and untranslated regions of mRNAs. This non-canonical translation is largely due to start codon plurality within the same RNA. This plurality is often due to the failure of some scanning ribosomes to recognize potential start codons leading to initiation downstream-a process termed leaky scanning. Codons other than AUG (non-AUG) are particularly leaky due to their inefficiency. Here we discuss our current understanding of non-AUG initiation. We argue for a near-ubiquitous role of non-AUG initiation in shaping the dynamic composition of mammalian proteomes.