A Regression-Based Analysis of Ribosome-Profiling Data Reveals a Conserved Complexity to Mammalian Translation

Identification of unannotated coding sequences and their physiological functions

Most protein-coding sequences (CDSs) are predicted sequences based on criteria such as a size sufficient to encode a product of at least 100 amino acids and with translation starting at an AUG initiation codon. However, recent studies based on ribosome profiling and mass spectrometry have shown that several RNAs annotated as long noncoding RNAs (lncRNAs) are actually translated to generate polypeptides of fewer than 100 amino acids, and that many proteins are translated from near-cognate initiation codons such as CUG and GUG. Furthermore, studies of genetically engineered mouse models have revealed that such polypeptides and proteins contribute to diverse physiological processes. In this review, we describe the latest methods for the identification of unannotated CDSs and provide examples of their physiological functions.

Identification and analysis of smORFs in Chlamydomonas reinhardtii

Small open reading frames (smORFs) have been acknowledged as an important partner in organism functions ranging from bacteria to higher eukaryotes. However, lack of investigation of smORFs in green algae, despite their importance in ecology and evolution. We applied bioinformatic analysis, ribosome profiling, and small peptide proteomics to provide a genome-wide and high-confident smORF database in the model green alga Chlamydomonas reinhardtii. The whole genome was screened first to mine potential coding smORFs. Then conservative analysis, ribosome profiling, and proteomics data were processed to identify conserved smORFs and generate translation evidence. The combination of procedures resulted in 2014 smORFs that might exist in the C. reinhardtii genome. The expression of smORFs in Cd treatment suggested that two smORFs might participate in redox reaction, three in inorganic phosphate transport, and one in DNA repair under stress. Our study built a genome-widely database in C. reinhardtii, providing target smORFs for further research.

Complex regulation of Gephyrin splicing is a determinant of inhibitory postsynaptic diversity

Gephyrin (GPHN) regulates the clustering of postsynaptic components at inhibitory synapses and is involved in pathophysiology of neuropsychiatric disorders. Here, we uncover an extensive diversity of GPHN transcripts that are tightly controlled by splicing during mouse and human brain development. Proteomic analysis reveals at least a hundred isoforms of GPHN incorporated at inhibitory Glycine and gamma-aminobutyric acid A receptors containing synapses. They exhibit different localization and postsynaptic clustering properties, and altering the expression level of one isoform is sufficient to affect the number, size, and density of inhibitory synapses in cerebellar Purkinje cells. Furthermore, we discovered that splicing defects reported in neuropsychiatric disorders are carried by multiple alternative GPHN transcripts, demonstrating the need for a thorough analysis of the GPHN transcriptome in patients. Overall, we show that alternative splicing of GPHN is an important genetic variation to consider in neurological diseases and a determinant of the diversity of postsynaptic inhibitory synapses.

The protein gephyrin is involved in organizing synapses. Here, the authors show how different transcripts of gephyrin form and regulate inhibitory synapses.

Proteomic-driven identification of short open reading frame-encoded peptides

Accumulating evidence has shown that a large number of short open reading frames (sORFs) also have the ability to encode proteins. The discovery of sORFs opens up a new research area, leading to the identification and functional study of sORF encoded peptides (SEPs) at the omics level. Besides bioinformatics prediction and ribosomal profiling, mass spectrometry (MS) has become a significant tool as it directly detects the sequence of SEPs. Though MS-based proteomics methods have proved to be effective for qualitative and quantitative analysis of SEPs, the detection of SEPs is still a great challenge due to their low abundance and short sequence. To illustrate the progress in method development, we described and discussed the main steps of large-scale proteomics identification of SEPs, including SEP extraction and enrichment, MS detection, data processing and quality control, quantification, and function prediction and validation methods. This article is protected by copyright. All rights reserved.

RiboReport - benchmarking tools for ribosome profiling-based identification of open reading frames in bacteria

Small proteins encoded by short open reading frames (ORFs) with 50 codons or fewer are emerging as an important class of cellular macromolecules in diverse organisms. However, they often evade detection by proteomics or in silico methods. Ribosome profiling (Ribo-seq) has revealed widespread translation in genomic regions previously thought to be non-coding, driving the development of ORF detection tools using Ribo-seq data. However, only a handful of tools have been designed for bacteria, and these have not yet been systematically compared. Here, we aimed to identify tools that use Ribo-seq data to correctly determine the translational status of annotated bacterial ORFs and also discover novel translated regions with high sensitivity. To this end, we generated a large set of annotated ORFs from four diverse bacterial organisms, manually labeled for their translation status based on Ribo-seq data, which are available for future benchmarking studies. This set was used to investigate the predictive performance of seven Ribo-seq-based ORF detection tools (REPARATION_blast, DeepRibo, Ribo-TISH, PRICE, smORFer, ribotricer and SPECtre), as well as IRSOM, which uses coding potential and RNA-seq coverage only. DeepRibo and REPARATION_blast robustly predicted translated ORFs, including sORFs, with no significant difference for ORFs in close proximity to other genes versus stand-alone genes. However, no tool predicted a set of novel, experimentally verified sORFs with high sensitivity. Start codon predictions with smORFer show the value of initiation site profiling data to further improve the sensitivity of ORF prediction tools in bacteria. Overall, we find that bacterial tools perform well for sORF detection, although there is potential for improving their performance, applicability, usability and reproducibility.

Small open reading frames in plant research: from prediction to functional characterization

Gene prediction is a laborious and time-consuming task. The advancement of sequencing technologies and bioinformatics tools, coupled with accelerated rate of ribosome profiling and mass spectrometry development, have made identification of small open reading frames (sORFs) (< 100 codons) in various plant genomes possible. The past 50 years have seen sORFs being isolated from many organisms. However, to date, a comprehensive sORF annotation pipeline is as yet unavailable, hence, addressed in our review. Here, we also provide current information on classification and functions of plant sORFs and their potential applications in crop improvement programs.

The dark proteome: translation from noncanonical open reading frames

Omics-based technologies have revolutionized our understanding of the coding potential of the genome. In particular, these studies revealed widespread unannotated open reading frames (ORFs) throughout genomes and that these regions have the potential to encode novel functional (micro-)proteins and/or hold regulatory roles. However, despite their genomic prevalence, relatively few of these noncanonical ORFs have been functionally characterized, likely in part due to their under-recognition by the broader scientific community. The few that have been investigated in detail have demonstrated their essentiality in critical and divergent biological processes. As such, here we aim to discuss recent advances in understanding the diversity of noncanonical ORFs and their roles, as well as detail biologically important examples within the context of the mammalian genome.

Trinucleotide CGG Repeat Diseases: An Expanding Field of Polyglycine Proteins?

Microsatellites are repeated DNA sequences of 3–6 nucleotides highly variable in length and sequence and that have important roles in genomes regulation and evolution. However, expansion of a subset of these microsatellites over a threshold size is responsible of more than 50 human genetic diseases. Interestingly, some of these disorders are caused by expansions of similar sequences, sizes and localizations and present striking similarities in clinical manifestations and histopathological features, which suggest a common mechanism of disease. Notably, five identical CGG repeat expansions, but located in different genes, are the causes of fragile X-associated tremor/ataxia syndrome (FXTAS), neuronal intranuclear inclusion disease (NIID), oculopharyngodistal myopathy type 1 to 3 (OPDM1-3) and oculopharyngeal myopathy with leukoencephalopathy (OPML), which are neuromuscular and neurodegenerative syndromes with overlapping symptoms and similar histopathological features, notably the presence of characteristic eosinophilic ubiquitin-positive intranuclear inclusions. In this review we summarize recent finding in neuronal intranuclear inclusion disease and FXTAS, where the causing CGG expansions were found to be embedded within small upstream ORFs (uORFs), resulting in their translation into novel proteins containing a stretch of polyglycine (polyG). Importantly, expression of these polyG proteins is toxic in animal models and is sufficient to reproduce the formation of ubiquitin-positive intranuclear inclusions. These data suggest the existence of a novel class of human genetic pathology, the polyG diseases, and question whether a similar mechanism may exist in other diseases, notably in OPDM and OPML.

Micropeptides translated from putative long non-coding RNAs

Long non-coding RNAs (lncRNAs) transcribed in mammals and eukaryotes were thought to have no protein coding capability. However, recent studies have suggested that plenty of lncRNAs are mis-annotated and virtually contain coding sequences which are translated into functional peptides by ribosomal machinery, and these functional peptides are called micropeptides or small peptides. Here we review the rapidly advancing field of micropeptides translated from putative lncRNAs, describe the strategies for their identification, and elucidate their critical roles in many fundamental biological processes. We also discuss the prospects of research in micropeptides and the potential applications of micropeptides.

Unannotated proteins expand the MHC-I-restricted immunopeptidome in cancer

Tumor-associated epitopes presented on MHC-I that can activate the immune system against cancer cells are typically identified from annotated protein-coding regions of the genome, but whether peptides originating from novel or unannotated open reading frames (nuORFs) can contribute to antitumor immune responses remains unclear. Here we show that peptides originating from nuORFs detected by ribosome profiling of malignant and healthy samples can be displayed on MHC-I of cancer cells, acting as additional sources of cancer antigens. We constructed a high-confidence database of translated nuORFs across tissues (nuORFdb) and used it to detect 3,555 translated nuORFs from MHC-I immunopeptidome mass spectrometry analysis, including peptides that result from somatic mutations in nuORFs of cancer samples as well as tumor-specific nuORFs translated in melanoma, chronic lymphocytic leukemia and glioblastoma. NuORFs are an unexplored pool of MHC-I-presented, tumor-specific peptides with potential as immunotherapy targets.

Translation and emerging functions of non-coding RNAs in inflammation and immunity

Regulatory non‐coding RNAs (ncRNAs) including small non‐coding RNAs (sRNAs), long non‐coding RNAs (lncRNAs), and circular RNAs (circRNAs) have gained considerable attention in the last few years. This is mainly due to their condition‐ and tissue‐specific expression and their various modes of action, which suggests them as promising biomarkers and therapeutic targets. One important mechanism of ncRNAs to regulate gene expression is through translation of short open reading frames (sORFs). These sORFs can be located in lncRNAs, in non‐translated regions of mRNAs where upstream ORFs (uORFs) represent the majority, or in circRNAs. Regulation of their translation can function as a quick way to adapt protein production to changing cellular or environmental cues, and can either depend solely on the initiation and elongation of translation, or on the roles of the produced functional peptides. Due to the experimental challenges to pinpoint translation events and to detect the produced peptides, translational regulation through regulatory RNAs is not well studied yet. In the case of circRNAs, they have only recently started to be recognized as regulatory molecules instead of mere artifacts of RNA biosynthesis. Of the many roles described for regulatory ncRNAs, we will focus here on their regulation during inflammation and in immunity.

Small Open Reading Frames, How to Find Them and Determine Their Function

Advances in genomics and molecular biology have revealed an abundance of small open reading frames (sORFs) across all types of transcripts. While these sORFs are often assumed to be non-functional, many have been implicated in physiological functions and a significant number of sORFs have been described in human diseases. Thus, sORFs may represent a hidden repository of functional elements that could serve as therapeutic targets. Unlike protein-coding genes, it is not necessarily the encoded peptide of an sORF that enacts its function, sometimes simply the act of translating an sORF might have a regulatory role. Indeed, the most studied sORFs are located in the 5′UTRs of coding transcripts and can have a regulatory impact on the translation of the downstream protein-coding sequence. However, sORFs have also been abundantly identified in non-coding RNAs including lncRNAs, circular RNAs and ribosomal RNAs suggesting that sORFs may be diverse in function. Of the many different experimental methods used to discover sORFs, the most commonly used are ribosome profiling and mass spectrometry. These can confirm interactions between transcripts and ribosomes and the production of a peptide, respectively. Extensions to ribosome profiling, which also capture scanning ribosomes, have further made it possible to see how sORFs impact the translation initiation of mRNAs. While high-throughput techniques have made the identification of sORFs less difficult, defining their function, if any, is typically more challenging. Together, the abundance and potential function of many of these sORFs argues for the necessity of including sORFs in gene annotations and systematically characterizing these to understand their potential functional roles. In this review, we will focus on the high-throughput methods used in the detection and characterization of sORFs and discuss techniques for validation and functional characterization.

Identifying Small Open Reading Frames in Prokaryotes with Ribosome Profiling

Small proteins encoded by open reading frames (ORFs) shorter than 50 codons (small ORFs [sORFs]) are often overlooked by annotation engines and are difficult to characterize using traditional biochemical techniques. Ribosome profiling has tremendous potential to empirically improve the annotations of prokaryotic genomes. Recent improvements in ribosome profiling methods for bacterial model organisms have revealed many new sORFs in well-characterized genomes. Antibiotics that trap ribosomes just after initiation have played a key role in these developments by allowing the unambiguous identification of the start codons (and, hence, the reading frame) for novel ORFs. Here, we describe these new methods and highlight critical controls and considerations for adapting ribosome profiling to different prokaryotic species.

The Small Open Reading Frame-Encoded Peptides: Advances in Methodologies and Functional Studies

Small open reading frames (sORFs) are an important class of genes with less than 100 codons. They were historically annotated as noncoding or even junk sequences. In recent years, accumulating evidence suggests that sORFs could encode a considerable number of polypeptides, many of which play important roles in both physiology and disease pathology. However, it has been technically challenging to directly detect sORF-encoded peptides (SEPs). Here, we discuss the latest advances in methodologies for identifying SEPs with mass spectrometry, as well as the progress on functional studies of SEPs.

The translatome of neuronal cell bodies, dendrites, and axons

To form synaptic connections and store information, neurons continuously remodel their proteomes. The impressive length of dendrites and axons imposes logistical challenges to maintain synaptic proteins at locations remote from the transcription source (the nucleus). The discovery of thousands of messenger RNAs (mRNAs) near synapses suggested that neurons overcome distance and gain autonomy by producing proteins locally. It is not generally known, however, if, how, and when localized mRNAs are translated into protein. To investigate the translational landscape in neuronal subregions, we performed simultaneous RNA sequencing (RNA-seq) and ribosome sequencing (Ribo-seq) from microdissected rodent brain slices to identify and quantify the transcriptome and translatome in cell bodies (somata) as well as dendrites and axons (neuropil). Thousands of transcripts were differentially translated between somatic and synaptic regions, with many scaffold and signaling molecules displaying increased translation levels in the neuropil. Most translational changes between compartments could be accounted for by differences in RNA abundance. Pervasive translational regulation was observed in both somata and neuropil influenced by specific mRNA features (e.g., untranslated region [UTR] length, RNA-binding protein [RBP] motifs, and upstream open reading frames [uORFs]). For over 800 mRNAs, the dominant source of translation was the neuropil. We constructed a searchable and interactive database for exploring mRNA transcripts and their translation levels in the somata and neuropil [MPI Brain Research, The mRNA translation landscape in the synaptic neuropil. https://public.brain.mpg.de/dashapps/localseq/ Accessed 5 October 2021]. Overall, our findings emphasize the substantial contribution of local translation to maintaining synaptic protein levels and indicate that on-site translational control is an important mechanism to control synaptic strength.

Unraveling the hidden role of a uORF-encoded peptide as a kinase inhibitor of PKCs

Approximately 40% of human messenger RNAs (mRNAs) contain upstream open reading frames (uORFs) in their 5' untranslated regions. Some of these uORF sequences, thought to attenuate scanning ribosomes or lead to mRNA degradation, were recently shown to be translated, although the function of the encoded peptides remains unknown. Here, we show a uORF-encoded peptide that exhibits kinase inhibitory functions. This uORF, upstream of the protein kinase C-eta (PKC-η) main ORF, encodes a peptide (uPEP2) containing the typical PKC pseudosubstrate motif present in all PKCs that autoinhibits their kinase activity. We show that uPEP2 directly binds to and selectively inhibits the catalytic activity of novel PKCs but not of classical or atypical PKCs. The endogenous deletion of uORF2 or its overexpression in MCF-7 cells revealed that the endogenously translated uPEP2 reduces the protein levels of PKC-η and other novel PKCs and restricts cell proliferation. Functionally, treatment of breast cancer cells with uPEP2 diminished cell survival and their migration and synergized with chemotherapy by interfering with the response to DNA damage. Furthermore, in a xenograft of MDA-MB-231 breast cancer tumor in mice models, uPEP2 suppressed tumor progression, invasion, and metastasis. Tumor histology showed reduced proliferation, enhanced cell death, and lower protein expression levels of novel PKCs along with diminished phosphorylation of PKC substrates. Hence, our study demonstrates that uORFs may encode biologically active peptides beyond their role as translation regulators of their downstream ORFs. Together, we point to a unique function of a uORF-encoded peptide as a kinase inhibitor, pertinent to cancer therapy.

SmProt: A Reliable Repository with Comprehensive Annotation of Small Proteins Identified from Ribosome Profiling

Small proteins specifically refer to proteins consisting of less than 100 amino acids translated from small open reading frames (sORFs), which were usually missed in previous genome annotation. The significance of small proteins has been revealed in current years, along with the discovery of their diverse functions. However, systematic annotation of small proteins is still insufficient. SmProt was specially developed to provide valuable information on small proteins for scientific community. Here we present the update of SmProt, which emphasizes reliability of translated sORFs, genetic variants in translated sORFs, disease-specific sORF translation events or sequences, and remarkably increased data volume. More components such as non-ATG translation initiation, function, and new sources are also included. SmProt incorporated 638,958 unique small proteins curated from 3,165,229 primary records, which were computationally predicted from 419 ribosome profiling (Ribo-seq) datasets or collected from literature and other sources from 370 cell lines or tissues in 8 species (Homo sapiens, Mus musculus, Rattus norvegicus, Drosophila melanogaster, Danio rerio, Saccharomyces cerevisiae, Caenorhabditis elegans, and Escherichia coli). In addition, small protein families identified from human microbiomes were also collected. All datasets in SmProt are free to access, and available for browse, search, and bulk downloads at http://bigdata.ibp.ac.cn/SmProt/.

Exploring Evidence of Non-coding RNA Translation With Trips-Viz and GWIPS-Viz Browsers

Detection of translation in so-called non-coding RNA provides an opportunity for identification of novel bioactive peptides and microproteins. The main methods used for these purposes are ribosome profiling and mass spectrometry. A number of publicly available datasets already exist for a substantial number of different cell types grown under various conditions, and public data mining is an attractive strategy for identification of translation in non-coding RNAs. Since the analysis of publicly available data requires intensive data processing, several data resources have been created recently for exploring processed publicly available data, such as OpenProt, GWIPS-viz, and Trips-Viz. In this work we provide a detailed demonstration of how to use the latter two tools for exploring experimental evidence for translation of RNAs hitherto classified as non-coding. For this purpose, we use a set of transcripts with substantially different patterns of ribosome footprint distributions. We discuss how certain features of these patterns can be used as evidence for or against genuine translation. During our analysis we concluded that the MTLN mRNA, previously misannotated as lncRNA LINC00116, likely encodes only a short proteoform expressed from shorter RNA transcript variants.

Trips-Viz: an environment for the analysis of public and user-generated ribosome profiling data

Trips-Viz (https://trips.ucc.ie/) is an interactive platform for the analysis and visualization of ribosome profiling (Ribo-Seq) and shotgun RNA sequencing (RNA-seq) data. This includes publicly available and user generated data, hence Trips-Viz can be classified as a database and as a server. As a database it provides access to many processed Ribo-Seq and RNA-seq data aligned to reference transcriptomes which has been expanded considerably since its inception. Here, we focus on the server functionality of Trips-viz which also has been greatly improved. Trips-viz now enables visualisation of proteomics data from a large number of processed mass spectrometry datasets. It can be used to support translation inferred from Ribo-Seq data. Users are now able to upload a custom reference transcriptome as well as data types other than Ribo-Seq/RNA-Seq. Incorporating custom data has been streamlined with RiboGalaxy (https://ribogalaxy.ucc.ie/) integration. The other new functionality is the rapid detection of translated open reading frames (ORFs) through a simple easy to use interface. The analysis of differential expression has been also improved via integration of DESeq2 and Anota2seq in addition to a number of other improvements of existing Trips-viz features.

MetamORF: a repository of unique short open reading frames identified by both experimental and computational approaches for gene and metagene analyses

The development of high-throughput technologies revealed the existence of non-canonical short open reading frames (sORFs) on most eukaryotic ribonucleic acids. They are ubiquitous genetic elements conserved across species and suspected to be involved in numerous cellular processes. MetamORF (https://metamorf.hb.univ-amu.fr/) aims to provide a repository of unique sORFs identified in the human and mouse genomes with both experimental and computational approaches. By gathering publicly available sORF data, normalizing them and summarizing redundant information, we were able to identify a total of 1 162 675 unique sORFs. Despite the usual characterization of ORFs as short, upstream or downstream, there is currently no clear consensus regarding the definition of these categories. Thus, the data have been reprocessed using a normalized nomenclature. MetamORF enables new analyses at locus, gene, transcript and ORF levels, which should offer the possibility to address new questions regarding sORF functions in the future. The repository is available through an user-friendly web interface, allowing easy browsing, visualization, filtering over multiple criteria and export possibilities. sORFs can be searched starting from a gene, a transcript and an ORF ID, looking in a genome area or browsing the whole repository for a species. The database content has also been made available through track hubs at UCSC Genome Browser. Finally, we demonstrated an enrichment of genes harboring upstream ORFs among genes expressed in response to reticular stress. Database URL  https://metamorf.hb.univ-amu.fr/.

MS-Based HLA-II Peptidomics Combined With Multiomics Will Aid the Development of Future Immunotherapies

Immunotherapies have emerged to treat diseases by selectively modulating a patient's immune response. Although the roles of T and B cells in adaptive immunity have been well studied, it remains difficult to select targets for immunotherapeutic strategies. Because human leukocyte antigen class II (HLA-II) peptides activate CD4+ T cells and regulate B cell activation, proliferation, and differentiation, these peptide antigens represent a class of potential immunotherapy targets and biomarkers. To better understand the molecular basis of how HLA-II antigen presentation is involved in disease progression and treatment, systematic HLA-II peptidomics combined with multiomic analyses of diverse cell types in healthy and diseased states is required. For this reason, MS-based innovations that facilitate investigations into the interplay between disease pathologies and the presentation of HLA-II peptides to CD4+ T cells will aid in the development of patient-focused immunotherapies.

smORFer: a modular algorithm to detect small ORFs in prokaryotes

Emerging evidence places small proteins (≤50 amino acids) more centrally in physiological processes. Yet, their functional identification and the systematic genome annotation of their cognate small open-reading frames (smORFs) remains challenging both experimentally and computationally. Ribosome profiling or Ribo-Seq (that is a deep sequencing of ribosome-protected fragments) enables detecting of actively translated open-reading frames (ORFs) and empirical annotation of coding sequences (CDSs) using the in-register translation pattern that is characteristic for genuinely translating ribosomes. Multiple identifiers of ORFs that use the 3-nt periodicity in Ribo-Seq data sets have been successful in eukaryotic smORF annotation. They have difficulties evaluating prokaryotic genomes due to the unique architecture (e.g. polycistronic messages, overlapping ORFs, leaderless translation, non-canonical initiation etc.). Here, we present a new algorithm, smORFer, which performs with high accuracy in prokaryotic organisms in detecting putative smORFs. The unique feature of smORFer is that it uses an integrated approach and considers structural features of the genetic sequence along with in-frame translation and uses Fourier transform to convert these parameters into a measurable score to faithfully select smORFs. The algorithm is executed in a modular way, and dependent on the data available for a particular organism, different modules can be selected for smORF search.

SARS-CoV-2 uses a multipronged strategy to impede host protein synthesis

The coronavirus SARS-CoV-2 is the cause of the ongoing pandemic of COVID-19¹. Coronaviruses have developed a variety of mechanisms to repress host mRNA translation to allow the translation of viral mRNA, and concomitantly block the cellular innate immune response^2,3. Although several different proteins of SARS-CoV-2 have previously been implicated in shutting off host expression^4-7, a comprehensive picture of the effects of SARS-CoV-2 infection on cellular gene expression is lacking. Here we combine RNA sequencing, ribosome profiling and metabolic labelling of newly synthesized RNA to comprehensively define the mechanisms that are used by SARS-CoV-2 to shut off cellular protein synthesis. We show that infection leads to a global reduction in translation, but that viral transcripts are not preferentially translated. Instead, we find that infection leads to the accelerated degradation of cytosolic cellular mRNAs, which facilitates viral takeover of the mRNA pool in infected cells. We reveal that the translation of transcripts that are induced in response to infection (including innate immune genes) is impaired. We demonstrate this impairment is probably mediated by inhibition of nuclear mRNA export, which prevents newly transcribed cellular mRNA from accessing ribosomes. Overall, our results uncover a multipronged strategy that is used by SARS-CoV-2 to take over the translation machinery and to suppress host defences.

ORFLine: a bioinformatic pipeline to prioritise small open reading frames identifies candidate secreted small proteins from lymphocytes

MOTIVATION

The annotation of small open reading frames (smORFs) of less than 100 codons (<300 nucleotides) is challenging due to the large number of such sequences in the genome.

RESULTS

In this study, we developed a computational pipeline, which we have named ORFLine, that stringently identifies smORFs and classifies them according to their position within transcripts. We identified a total of 5744 unique smORFs in datasets from mouse B and T lymphocytes and systematically characterized them using ORFLine. We further searched smORFs for the presence of a signal peptide, which predicted known secreted chemokines as well as novel micropeptides. Four novel micropeptides show evidence of secretion and are therefore candidate mediators of immunoregulatory functions.

AVAILABILITY

Freely available on the web at https://github.com/boboppie/ORFLine.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The role of micropeptides in biology

Micropeptides are small polypeptides coded by small open-reading frames. Progress in computational biology and the analyses of large-scale transcriptomes and proteomes have revealed that mammalian genomes produce a large number of transcripts encoding micropeptides. Many of these have been previously annotated as long noncoding RNAs. The role of micropeptides in cellular homeostasis maintenance has been demonstrated. This review discusses different types of micropeptides as well as methods to identify them, such as computational approaches, ribosome profiling, and mass spectrometry.

uORF-seqr: A Machine Learning-Based Approach to the Identification of Upstream Open Reading Frames in Yeast

The identification of upstream open reading frames (uORFs) using ribosome profiling data is complicated by several factors such as the noise inherent to the procedure, the substantial increase in potential translation initiation sites (and false positives) when one includes non-canonical start codons, and the paucity of molecularly validated uORFs. Here we present uORF-seqr, a novel machine learning algorithm that uses ribosome profiling data, in conjunction with RNA-seq data, as well as transcript aware genome annotation files to identify statistically significant AUG and near-cognate codon uORFs.

Global mapping of translation initiation sites by TIS profiling in budding yeast

Translation initiation site (TIS) profiling allows for the genome-wide identification of TISs in vivo by exclusively capturing mRNA fragments within ribosomes that have just completed translation initiation. It leverages translation inhibitors, such as harringtonine and lactimidomycin (LTM), that preferentially capture ribosomes at start codon positions, protecting TIS-derived mRNA fragments from nuclease digestion. Here, we describe a step-by-step protocol for TIS profiling in LTM-treated budding yeast that we developed to identify TISs and open reading frames in vegetative and meiotic cells. For complete details on the use and execution of this protocol, please refer to Eisenberg et al. (2020).

Most non-canonical proteins uniquely populate the proteome or immunopeptidome

Combining RNA sequencing, ribosome profiling, and mass spectrometry, we elucidate the contribution of non-canonical translation to the proteome and major histocompatibility complex (MHC) class I immunopeptidome. Remarkably, of 14,498 proteins identified in three human B cell lymphomas, 2,503 are non-canonical proteins. Of these, 28% are novel isoforms and 72% are cryptic proteins encoded by ostensibly non-coding regions (60%) or frameshifted canonical genes (12%). Cryptic proteins are translated as efficiently as canonical proteins, have more predicted disordered residues and lower stability, and critically generate MHC-I peptides 5-fold more efficiently per translation event. Translating 5' "untranslated" regions hinders downstream translation of genes involved in transcription, translation, and antiviral responses. Novel protein isoforms show strong enrichment for signaling pathways deregulated in cancer. Only a small fraction of cryptic proteins detected in the proteome contribute to the MHC-I immunopeptidome, demonstrating the high preferential access of cryptic defective ribosomal products to the class I pathway.

Best practices to ensure robust investigation of circular RNAs: pitfalls and tips

Pre-mRNAs from thousands of eukaryotic genes can be non-canonically spliced to generate circular RNAs (circRNAs) that have covalently linked ends. Most mature circular RNAs are expressed at low levels, but some have known physiological functions and/or accumulate to higher levels than their associated linear mRNAs. These observations have sparked great interest into this class of previously underappreciated RNAs and prompted the development of new experimental approaches to study them, especially methods to measure or modulate circular RNA expression levels. Nonetheless, each of these approaches has caveats and potential pitfalls that must be controlled for when designing experiments and interpreting results. Here, we provide practical advice, tips, and suggested guidelines for performing robust identification, validation, and functional characterization of circular RNAs. Beyond promoting rigor and reproducibility, these suggestions should help bring clarity to the field, especially how circular RNAs function and whether these transcripts may sponge microRNAs/proteins or serve as templates for translation.

Revisiting sORFs: overcoming challenges to identify and characterize functional microproteins

Short ORFs (sORFs), that is, occurrences of a start and stop codon within 100 codons or less, can be found in organisms of all domains of life, outnumbering annotated protein-coding ORFs by orders of magnitude. Even though functional proteins smaller than 100 amino acids are known, the coding potential of sORFs has often been overlooked, as it is not trivial to predict and test for functionality within the large number of sORFs. Recent advances in ribosome profiling and mass spectrometry approaches, together with refined bioinformatic predictions, have enabled a huge leap forward in this field and identified thousands of likely coding sORFs. A relatively low number of small proteins or microproteins produced from these sORFs have been characterized so far on the molecular, structural, and/or mechanistic level. These however display versatile and, in some cases, essential cellular functions, allowing for the exciting possibility that many more, previously unknown small proteins might be encoded in the genome, waiting to be discovered. This review will give an overview of the steadily growing microprotein field, focusing on eukaryotic small proteins. We will discuss emerging themes in the molecular action of microproteins, as well as advances and challenges in microprotein identification and characterization.

A Mini-review of Computational Approaches to Predict Functions and Findings of Novel Micro Peptides

:

New techniques in bioinformatics and the study of the transcriptome at a wide-scale have uncovered the fact that a large part of the genome is being translated than recently perceived thoughts and research, bringing about the creation of a various quantity of RNA with proteincoding and noncoding potential. A lot of RNA particles have been considered as noncoding due to many reasons, according to developing proofs. Like many sORFs that encode many functional micro peptides have neglected due to their tiny sizes.

:

Advanced studies reveal many major biological functions of these sORFs and their encoded micro peptides in a different and wide range of species. All the achievement in the identification of these sORFs and micro peptides is due to the progressive bioinformatics and high-throughput sequencing methods. This field has pulled in more consideration due to the detection of a large number of more sORFs and micro peptides. Nowadays, COVID-19 grabs all the attention of science as it is a sudden outbreak. sORFs of COVID-19 should be revealed for new ways to understand this virus. This review discusses ongoing progress in the systems for the identification and distinguishing proof of sORFs and micro peptides.

Codon Resolution Analysis of Ribosome Profiling Data

Translation is a central biological process in living cells. Ribosome profiling approach enables assessing translation on a global, cell-wide level. Extracting versatile information from the ribosome profiling data usually requires specialized expertise for handling the sequencing data that is not available to the broad community of experimentalists. Here, we provide an easy-to-use and modifiable workflow that uses a small set of commands and enables full data analysis in a standardized way, including precise positioning of the ribosome-protected fragments, for determining codon-specific translation features. The workflow is complemented with simple step-by-step explanations and is accessible to scientists with no computational background.

RiboDiPA: a novel tool for differential pattern analysis in Ribo-seq data

Ribosome profiling, also known as Ribo-seq, has become a popular approach to investigate regulatory mechanisms of translation in a wide variety of biological contexts. Ribo-seq not only provides a measurement of translation efficiency based on the relative abundance of ribosomes bound to transcripts, but also has the capacity to reveal dynamic and local regulation at different stages of translation based on positional information of footprints across individual transcripts. While many computational tools exist for the analysis of Ribo-seq data, no method is currently available for rigorous testing of the pattern differences in ribosome footprints. In this work, we develop a novel approach together with an R package, RiboDiPA, for Differential Pattern Analysis of Ribo-seq data. RiboDiPA allows for quick identification of genes with statistically significant differences in ribosome occupancy patterns for model organisms ranging from yeast to mammals. We show that differential pattern analysis reveals information that is distinct and complimentary to existing methods that focus on translational efficiency analysis. Using both simulated Ribo-seq footprint data and three benchmark data sets, we illustrate that RiboDiPA can uncover meaningful pattern differences across multiple biological conditions on a global scale, and pinpoint characteristic ribosome occupancy patterns at single codon resolution.

Rules are made to be broken: a "simple" model organism reveals the complexity of gene regulation

Global methods for assaying translation have greatly improved our understanding of the protein-coding capacity of the genome. In particular, it is now possible to perform genome-wide and condition-specific identification of translation initiation sites through modified ribosome profiling methods that selectively capture initiating ribosomes. Here we discuss our recent study applying such an approach to meiotic and mitotic timepoints in the simple eukaryote, budding yeast, as an example of the surprising diversity of protein products-many of which are non-canonical-that can be revealed by such methods. We also highlight several key challenges in studying non-canonical protein isoforms that have precluded their prior systematic discovery. A growing body of work supports expanded use of empirical protein-coding region identification, which can help relieve some of the limitations and biases inherent to traditional genome annotation approaches. Our study also argues for the adoption of less static views of gene identity and a broader framework for considering the translational capacity of the genome.

Translation of noncoding RNAs and cancer

The human genome contains thousands of noncoding RNAs (ncRNAs), which are thought to lack open reading frames (ORFs) and cannot be translated. Some ncRNAs reportedly have important functions, including epigenetic regulation, chromatin remolding, protein modification, and RNA degradation, but the functions of most ncRNAs remain elusive. Through the application and development of ribosome profiling and sequencing technologies, an increasing number of studies have discovered the translation of ncRNAs. Although ncRNAs were initially defined as noncoding RNAs, a number of ncRNAs actually contain ORFs that are translated into peptides. Here, we summarize the available methods, tools, and databases for identifying and validating ncRNA-encoded peptides/proteins, and the recent findings regarding ncRNA-encoded small peptides/proteins in cancer are compiled and synthesized. Importantly, the role of ncRNA-encoding peptides/proteins has application prospects in cancer research, but some potential challenges remain unresolved. The aim of this review is to provide a theoretical basis that might promote the discovery of more peptides/proteins encoded by ncRNAs and aid the further development of novel diagnostic and prognostic cancer markers and therapeutic targets.

RiboToolkit: an integrated platform for analysis and annotation of ribosome profiling data to decode mRNA translation at codon resolution

Ribosome profiling (Ribo-seq) is a powerful technology for globally monitoring RNA translation; ranging from codon occupancy profiling, identification of actively translated open reading frames (ORFs), to the quantification of translational efficiency under various physiological or experimental conditions. However, analyzing and decoding translation information from Ribo-seq data is not trivial. Although there are many existing tools to analyze Ribo-seq data, most of these tools are designed for specific or limited functionalities and an easy-to-use integrated tool to analyze Ribo-seq data is lacking. Fortunately, the small size (26–34 nt) of ribosome protected fragments (RPFs) in Ribo-seq and the relatively small amount of sequencing data greatly facilitates the development of such a web platform, which is easy to manipulate for users with or without bioinformatic expertise. Thus, we developed RiboToolkit (http://rnabioinfor.tch.harvard.edu/RiboToolkit), a convenient, freely available, web-based service to centralize Ribo-seq data analyses, including data cleaning and quality evaluation, expression analysis based on RPFs, codon occupancy, translation efficiency analysis, differential translation analysis, functional annotation, translation metagene analysis, and identification of actively translated ORFs. Besides, easy-to-use web interfaces were developed to facilitate data analysis and intuitively visualize results. Thus, RiboToolkit will greatly facilitate the study of mRNA translation based on ribosome profiling.

Induction of an Alternative mRNA 5' Leader Enhances Translation of the Ciliopathy Gene Inpp5e and Resistance to Oncolytic Virus Infection

Residual cell-intrinsic innate immunity in cancer cells hampers infection with oncolytic viruses. Translational control of mRNA is an important feature of innate immunity, yet the identity of translationally regulated mRNAs functioning in host defense remains ill-defined. We report the translatomes of resistant murine "4T1" breast cancer cells infected with three of the most clinically advanced oncolytic viruses: herpes simplex virus 1, reovirus, and vaccinia virus. Common among all three infections are translationally de-repressed mRNAs, including Inpp5e, encoding an inositol 5-phosphatase that modifies lipid second messenger signaling. We find that viral infection induces the expression of an Inpp5e mRNA variant that lacks repressive upstream open reading frames (uORFs) within its 5' leader and is efficiently translated. Furthermore, we show that INPP5E contributes to antiviral immunity by altering virus attachment. These findings uncover a role for translational control through alternative 5' leader expression and assign an antiviral function to the ciliopathy gene Inpp5e.

Accurate detection of short and long active ORFs using Ribo-seq data

MOTIVATION

Ribo-seq, a technique for deep-sequencing ribosome-protected mRNA fragments, has enabled transcriptome-wide monitoring of translation in vivo. It has opened avenues for re-evaluating the coding potential of open reading frames (ORFs), including many short ORFs that were previously presumed to be non-translating. However, the detection of translating ORFs, specifically short ORFs, from Ribo-seq data, remains challenging due to its high heterogeneity and noise.

RESULTS

We present ribotricer, a method for detecting actively translating ORFs by directly leveraging the three-nucleotide periodicity of Ribo-seq data. Ribotricer demonstrates higher accuracy and robustness compared with other methods at detecting actively translating ORFs including short ORFs on multiple published datasets across species inclusive of Arabidopsis, Caenorhabditis elegans, Drosophila, human, mouse, rat, yeast and zebrafish.

AVAILABILITY AND IMPLEMENTATION

Ribotricer is available at https://github.com/smithlabcode/ribotricer. All analysis scripts and results are available at https://github.com/smithlabcode/ribotricer-results.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The coding capacity of SARS-CoV-2

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the cause of the ongoing coronavirus disease 2019 (COVID-19) pandemic¹. To understand the pathogenicity and antigenic potential of SARS-CoV-2 and to develop therapeutic tools, it is essential to profile the full repertoire of its expressed proteins. The current map of SARS-CoV-2 coding capacity is based on computational predictions and relies on homology with other coronaviruses. As the protein complement varies among coronaviruses, especially in regard to the variety of accessory proteins, it is crucial to characterize the specific range of SARS-CoV-2 proteins in an unbiased and open-ended manner. Here, using a suite of ribosome-profiling techniques^2-4, we present a high-resolution map of coding regions in the SARS-CoV-2 genome, which enables us to accurately quantify the expression of canonical viral open reading frames (ORFs) and to identify 23 unannotated viral ORFs. These ORFs include upstream ORFs that are likely to have a regulatory role, several in-frame internal ORFs within existing ORFs, resulting in N-terminally truncated products, as well as internal out-of-frame ORFs, which generate novel polypeptides. We further show that viral mRNAs are not translated more efficiently than host mRNAs; instead, virus translation dominates host translation because of the high levels of viral transcripts. Our work provides a resource that will form the basis of future functional studies.

Control of translation by eukaryotic mRNA transcript leaders-Insights from high-throughput assays and computational modeling

Eukaryotic gene expression is tightly regulated during translation of mRNA to protein. Mis-regulation of translation can lead to aberrant proteins which accumulate in cancers and cause neurodegenerative diseases. Foundational studies on model genes established fundamental roles for mRNA 5' transcript leader (TL) sequences in controlling ribosome recruitment, scanning, and initiation. TL cis-regulatory elements and their corresponding trans-acting factors control cap-dependent initiation under unstressed conditions. Under stress, cap-dependent initiation is suppressed, and specific mRNA structures and sequences promote translation of stress-responsive transcripts to remodel the proteome. In this review, we summarize current knowledge of TL functions in translation initiation. We focus on insights from high-throughput analyses of ribosome occupancy, mRNA structure, RNA Binding Protein occupancy, and massively parallel reporter assays. These data-driven approaches, coupled with computational analyses and modeling, have paved the way for a comprehensive understanding of TL functions. Finally, we will discuss areas of future research on the roles of mRNA sequences and structures in translation. This article is categorized under: Translation > Translation Mechanisms RNA Evolution and Genomics > Computational Analyses of RNA RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems.

Translation Initiation Site Profiling Reveals Widespread Synthesis of Non-AUG-Initiated Protein Isoforms in Yeast

Genomic analyses in budding yeast have helped define the foundational principles of eukaryotic gene expression. However, in the absence of empirical methods for defining coding regions, these analyses have historically excluded specific classes of possible coding regions, such as those initiating at non-AUG start codons. Here, we applied an experimental approach to globally annotate translation initiation sites in yeast and identified 149 genes with alternative N-terminally extended protein isoforms initiating from near-cognate codons upstream of annotated AUG start codons. These isoforms are produced in concert with canonical isoforms and translated with high specificity, resulting from initiation at only a small subset of possible start codons. The non-AUG initiation driving their production is enriched during meiosis and induced by low eIF5A, which is seen in this context. These findings reveal widespread production of non-canonical protein isoforms and unexpected complexity to the rules by which even a simple eukaryotic genome is decoded.

Alternative Translation Initiation Generates a Functionally Distinct Isoform of the Stress-Activated Protein Kinase MK2

Alternative translation is an important mechanism of post-transcriptional gene regulation leading to the expression of different protein isoforms originating from the same mRNA. Here, we describe an abundant long isoform of the stress/p38^MAPK-activated protein kinase MK2. This isoform is constitutively translated from an alternative CUG translation initiation start site located in the 5' UTR of its mRNA. The RNA helicase eIF4A1 is needed to ensure translation of the long and the known short isoforms of MK2, of which the molecular properties were determined. Only the short isoform phosphorylated Hsp27 in vivo, supported migration and stress-induced immediate early gene (IEG) expression. Interaction profiling revealed short-isoform-specific binding partners that were associated with migration. In contrast, the long isoform contains at least one additional phosphorylatable serine in its unique N terminus. In sum, our data reveal a longer isoform of MK2 with distinct physiological properties.

RiboMiner: a toolset for mining multi-dimensional features of the translatome with ribosome profiling data

Background

Ribosome profiling has been widely used for studies of translation under a large variety of cellular and physiological contexts. Many of these studies have greatly benefitted from a series of data-mining tools designed for dissection of the translatome from different aspects. However, as the studies of translation advance quickly, the current toolbox still falls in short, and more specialized tools are in urgent need for deeper and more efficient mining of the important and new features of the translation landscapes.

Results

Here, we present RiboMiner, a bioinformatics toolset for mining of multi-dimensional features of the translatome with ribosome profiling data. RiboMiner performs extensive quality assessment of the data and integrates a spectrum of tools for various metagene analyses of the ribosome footprints and for detailed analyses of multiple features related to translation regulation. Visualizations of all the results are available. Many of these analyses have not been provided by previous methods. RiboMiner is highly flexible, as the pipeline could be easily adapted and customized for different scopes and targets of the studies.

Conclusions

Applications of RiboMiner on two published datasets did not only reproduced the main results reported before, but also generated novel insights into the translation regulation processes. Therefore, being complementary to the current tools, RiboMiner could be a valuable resource for dissections of the translation landscapes and the translation regulations by mining the ribosome profiling data more comprehensively and with higher resolution. RiboMiner is freely available at https://github.com/xryanglab/RiboMiner and https://pypi.org/project/RiboMiner.

The small peptide world in long noncoding RNAs

Long noncoding RNAs (lncRNAs) are a group of transcripts that are longer than 200 nucleotides (nt) without coding potential. Over the past decade, tens of thousands of novel lncRNAs have been annotated in animal and plant genomes because of advanced high-throughput RNA sequencing technologies and with the aid of coding transcript classifiers. Further, a considerable number of reports have revealed the existence of stable, functional small peptides (also known as micropeptides), translated from lncRNAs. In this review, we discuss the methods of lncRNA classification, the investigations regarding their coding potential and the functional significance of the peptides they encode.

Reconsidering proteomic diversity with functional investigation of small ORFs and alternative ORFs

The discovery of functional yet non-annotated open reading frames (ORFs) throughout the genome of several species presents an unprecedented challenge in current genome annotation. These novel ORFs are shorter than annotated ones and many can be found on the same RNA, in opposition to current assumptions in annotation methodologies. Whilst the literature lacks consensus, these novel ORFs are commonly referred to as small ORFs (sORFs) or alternative ORFs (alt-ORFs). Unannotated ORFs represent an overlooked layer of complexity in the coding potential of genomes and are transforming our current vision of the nature of coding genes. In this review, we outline what constitutes a sORF or an alt-ORF and emphasize differences between both nomenclatures. We then describe complementary large-scale methods to accurately discover novel ORFs as well as yield functional insights on the novel proteins they encode. While serendipitous discoveries highlighted the functional importance of some novel ORFs, omics methods facilitate and improve their characterization to better understand physiological and pathological pathways. Functional annotation of sORFs, alt-ORFs and their corresponding microproteins will likely help fundamental and clinical research.

Computational methods for ribosome profiling data analysis

Since the introduction of the ribosome profiling technique in 2009 its popularity has greatly increased. It is widely used for the comprehensive assessment of gene expression and for studying the mechanisms of regulation at the translational level. As the number of ribosome profiling datasets being produced continues to grow, so too does the need for reliable software that can provide answers to the biological questions it can address. This review describes the computational methods and tools that have been developed to analyze ribosome profiling data at the different stages of the process. It starts with initial routine processing of raw data and follows with more specific tasks such as the identification of translated open reading frames, differential gene expression analysis, or evaluation of local or global codon decoding rates. The review pinpoints challenges associated with each step and explains the ways in which they are currently addressed. In addition it provides a comprehensive, albeit incomplete, list of publicly available software applicable to each step, which may be a beneficial starting point to those unexposed to ribosome profiling analysis. The outline of current challenges in ribosome profiling data analysis may inspire computational biologists to search for novel, potentially superior, solutions that will improve and expand the bioinformatician's toolbox for ribosome profiling data analysis. This article is characterized under: Translation > Ribosome Structure/Function RNA Evolution and Genomics > Computational Analyses of RNA Translation > Translation Mechanisms Translation > Translation Regulation.

Alternative ORFs and small ORFs: shedding light on the dark proteome

Traditional annotation of protein-encoding genes relied on assumptions, such as one open reading frame (ORF) encodes one protein and minimal lengths for translated proteins. With the serendipitous discoveries of translated ORFs encoded upstream and downstream of annotated ORFs, from alternative start sites nested within annotated ORFs and from RNAs previously considered noncoding, it is becoming clear that these initial assumptions are incorrect. The findings have led to the realization that genetic information is more densely coded and that the proteome is more complex than previously anticipated. As such, interest in the identification and characterization of the previously ignored 'dark proteome' is increasing, though we note that research in eukaryotes and bacteria has largely progressed in isolation. To bridge this gap and illustrate exciting findings emerging from studies of the dark proteome, we highlight recent advances in both eukaryotic and bacterial cells. We discuss progress in the detection of alternative ORFs as well as in the understanding of functions and the regulation of their expression and posit questions for future work.

Pervasive functional translation of noncanonical human open reading frames

Ribosome profiling has revealed pervasive but largely uncharacterized translation outside of canonical coding sequences (CDSs). In this work, we exploit a systematic CRISPR-based screening strategy to identify hundreds of noncanonical CDSs that are essential for cellular growth and whose disruption elicits specific, robust transcriptomic and phenotypic changes in human cells. Functional characterization of the encoded microproteins reveals distinct cellular localizations, specific protein binding partners, and hundreds of microproteins that are presented by the human leukocyte antigen system. We find multiple microproteins encoded in upstream open reading frames, which form stable complexes with the main, canonical protein encoded on the same messenger RNA, thereby revealing the use of functional bicistronic operons in mammals. Together, our results point to a family of functional human microproteins that play critical and diverse cellular roles.

Ribosome profiling of HEK293T cells overexpressing codon optimized coagulation factor IX

Ribosome profiling provides the opportunity to evaluate translation kinetics at codon level resolution. Here, we describe ribosome profiling data, generated from two HEK293T cell lines. The ribosome profiling data are composed of Ribo-seq (mRNA sequencing data from ribosome protected fragments) and RNA-seq data (total RNA sequencing). The two HEK293T cell lines each express a version of the F9 gene, both of which are translated into identical proteins in terms of their amino acid sequences. However, these F9 genes vary drastically in their codon usage and predicted mRNA structure. We also provide the pipeline that we used to analyze the data. Further analyzing this dataset holds great potential as it can be used i) to unveil insights into the composition and regulation of the transcriptome, ii) for comparison with other ribosome profiling datasets, iii) to measure the rate of protein synthesis across the proteome and identify differences in elongation rates, iv) to discover previously unidentified translation of peptides, v) to explore the effects of codon usage or codon context in translational kinetics and vi) to investigate cotranslational folding. Importantly, a unique feature of this dataset, compared to other available ribosome profiling data, is the presence of the F9 gene in two very distinct coding sequences.