Calibrated ribosome profiling assesses the dynamics of ribosomal flux on transcripts.

Ribosome profiling, which is based on deep sequencing of ribosome footprints, has served as a powerful tool for elucidating the regulatory mechanism of protein synthesis. However, the current method has substantial issues: contamination by rRNAs and the lack of appropriate methods to measure ribosome numbers in transcripts. Here, we overcome these hurdles through the development of "Ribo-FilterOut", which is based on the separation of footprints from ribosome subunits by ultrafiltration, and "Ribo-Calibration", which relies on external spike-ins of stoichiometrically defined mRNA-ribosome complexes. A combination of these approaches estimates the number of ribosomes on a transcript, the translation initiation rate, and the overall number of translation events before its decay, all in a genome-wide manner. Moreover, our method reveals the allocation of ribosomes under heat shock stress, during aging, and across cell types. Our strategy of modified ribosome profiling measures kinetic and stoichiometric parameters of cellular translation across the transcriptome.

Ribosome Profiling and RNA Sequencing Reveal Translation and Transcription Regulation under Acute Heat Stress in Rainbow Trout (Oncorhynchus mykiss, Walbaum, 1792) Liver.

Rainbow trout (Oncorhynchus mykiss, Walbaum, 1792) is an important economic cold-water fish that is susceptible to heat stress. To date, the heat stress response in rainbow trout is more widely understood at the transcriptional level, while little research has been conducted at the translational level. To reveal the translational regulation of heat stress in rainbow trout, in this study, we performed a ribosome profiling assay of rainbow trout liver under normal and heat stress conditions. Comparative analysis of the RNA-seq data with the ribosome profiling data showed that the folding changes in gene expression at the transcriptional level are moderately correlated with those at the translational level. In total, 1213 genes were significantly altered at the translational level. However, only 32.8% of the genes were common between both levels, demonstrating that heat stress is coordinated across both transcriptional and translational levels. Moreover, 809 genes exhibited significant differences in translational efficiency (TE), with the TE of these genes being considerably affected by factors such as the GC content, coding sequence length, and upstream open reading frame (uORF) presence. In addition, 3468 potential uORFs in 2676 genes were identified, which can potentially affect the TE of the main open reading frames. In this study, Ribo-seq and RNA-seq were used for the first time to elucidate the coordinated regulation of transcription and translation in rainbow trout under heat stress. These findings are expected to contribute novel data and theoretical insights to the international literature on the thermal stress response in fish.

Rp3: Ribosome profiling-assisted proteogenomics improves coverage and confidence during microprotein discovery.

There has been a dramatic increase in the identification of non-canonical translation and a significant expansion of the protein-coding genome. Among the strategies used to identify unannotated small Open Reading Frames (smORFs) that encode microproteins, Ribosome profiling (Ribo-Seq) is the gold standard for the annotation of novel coding sequences by reporting on smORF translation. In Ribo-Seq, ribosome-protected footprints (RPFs) that map to multiple genomic sites are removed since they cannot be unambiguously assigned to a specific genomic location. Furthermore, RPFs necessarily result in short (25-34 nucleotides) reads, increasing the chance of multi-mapping alignments, such that smORFs residing in these regions cannot be identified by Ribo-Seq. Moreover, it has been challenging to identify protein evidence for Ribo-Seq. To solve this, we developed Rp3, a pipeline that integrates proteogenomics and Ribosome profiling to provide unambiguous evidence for a subset of microproteins missed by current Ribo-Seq pipelines. Here, we show that Rp3 maximizes proteomics detection and confidence of microprotein-encoding smORFs.

Comparison of software packages for detecting unannotated translated small open reading frames by Ribo-seq.

Accurate and comprehensive annotation of microprotein-coding small open reading frames (smORFs) is critical to our understanding of normal physiology and disease. Empirical identification of translated smORFs is carried out primarily using ribosome profiling (Ribo-seq). While effective, published Ribo-seq datasets can vary drastically in quality and different analysis tools are frequently employed. Here, we examine the impact of these factors on identifying translated smORFs. We compared five commonly used software tools that assess open reading frame translation from Ribo-seq (RibORFv0.1, RibORFv1.0, RiboCode, ORFquant, and Ribo-TISH) and found surprisingly low agreement across all tools. Only ~2% of smORFs were called translated by all five tools, and ~15% by three or more tools when assessing the same high-resolution Ribo-seq dataset. For larger annotated genes, the same analysis showed ~74% agreement across all five tools. We also found that some tools are strongly biased against low-resolution Ribo-seq data, while others are more tolerant. Analyzing Ribo-seq coverage revealed that smORFs detected by more than one tool tend to have higher translation levels and higher fractions of in-frame reads, consistent with what was observed for annotated genes. Together these results support employing multiple tools to identify the most confident microprotein-coding smORFs and choosing the tools based on the quality of the dataset and the planned downstream characterization experiments of the predicted smORFs.

RiboGraph: an interactive visualization system for ribosome profiling data at read length resolution.

MOTIVATION: Ribosome profiling is a widely-used technique for measuring ribosome occupancy at nucleotide resolution. However, the need to analyze this data at nucleotide resolution introduces unique challenges in data visualization and analyses. RESULTS: In this study, we introduce RiboGraph, a dedicated visualization tool designed to work with .ribo files, a specialized and efficient format for ribosome occupancy data. Unlike existing solutions that rely on large alignment files and time-consuming preprocessing steps, RiboGraph operates on a purpose designed compact file type. This efficiency allows for interactive, real-time visualization at ribosome-protected fragment length resolution. By providing an integrated toolset, RiboGraph empowers researchers to conduct comprehensive visual analysis of ribosome occupancy data. AVAILABILITY AND IMPLEMENTATION: Source code, step-by-step installation instructions and links to documentation are available on GitHub: https://github.com/ribosomeprofiling/ribograph. On the same page, we provide test files and a step-by-step tutorial highlighting the key features of RiboGraph.

Identification of genes associated with sperm storage capacity in hens at different times after insemination by RNA-seq and Ribo-seq.

BACKGROUND: Sperm storage capacity (SSC) determines the duration of fertility in hens and is an important reproduction trait that cannot be ignored in production. Currently, the genetic mechanism of SSC is still unclear in hens. Therefore, to explore the genetic basis of SSC, we analyzed the uterus-vagina junction (UVJ) of hens with different SSC at different times after insemination by RNA-seq and Ribo-seq. RESULTS: Our results showed that 589, 596, and 527 differentially expressed genes (DEGs), 730, 783, and 324 differentially translated genes (DTGs), and 804, 625, and 467 differential translation efficiency genes (DTEGs) were detected on the 5th, 10th, and 15th days after insemination, respectively. In transcription levels, we found that the differences of SSC at different times after insemination were mainly reflected in the transmission of information between cells, the composition of intercellular adhesion complexes, the regulation of ion channels, the regulation of cellular physiological activities, the composition of cells, and the composition of cell membranes. In translation efficiency (TE) levels, the differences of SSC were mainly related to the physiological and metabolic activities in the cell, the composition of the organelle membrane, the physiological activities of oxidation, cell components, and cell growth processes. According to pathway analysis, SSC was related to neuroactive ligand-receptor interaction, histidine metabolism, and PPAR signaling pathway at the transcriptional level and glutathione metabolism, oxidative phosphorylation, calcium signaling pathway, cell adhesion molecules, galactose metabolism, and Wnt signaling pathway at the TE level. We screened candidate genes affecting SSC at transcriptional levels (COL4A4, MUC6, MCHR2, TACR1, AVPR1A, COL1A1, HK2, RB1, VIPR2, HMGCS2) and TE levels(COL4A4, MUC6, CYCS, NDUFA13, CYTB, RRM2, CAMK4, HRH2, LCT, GCK, GALT). Among them, COL4A4 and MUC6 were the key candidate genes differing in transcription, translation, and translation efficiency. CONCLUSIONS: Our study used the combined analysis of RNA-seq and Ribo-seq for the first time to investigate the SSC and reveal the physiological processes associated with SSC. The key candidate genes affecting SSC were screened, and the theoretical basis was provided for the analysis of the molecular regulation mechanism of SSC.

Impact of translational regulation on diel expression revealed by time-series ribosome profiling in Arabidopsis.

Plants have developed the ability to adjust to the day/night cycle through the expression of diel genes, which allow them to effectively respond to environmental changes and optimise their growth and development. Diel oscillations also have substantial implications in many physiological processes, including photosynthesis, floral development, and environmental stress responses. The expression of diel genes is regulated by a combination of the circadian clock and responses to environmental cues, such as light and temperature. A great deal of information is available on the transcriptional regulation of diel gene expression. However, the extent to which translational regulation is involved in controlling diel changes in expression is not yet clear. To investigate the impact of translational regulation on diel expression, we conducted Ribo-seq and RNA-seq analyses on a time-series sample of Arabidopsis shoots cultivated under a 12 h light/dark cycle. Our results showed that translational regulation is involved in about 71% of the genes exhibiting diel changes in mRNA abundance or translational activity, including clock genes, many of which are subject to both translational and transcriptional control. They also revealed that the diel expression of glycosylation and ion-transporter-related genes is mainly established through translational regulation. The expression of several diel genes likely subject to translational regulation through upstream open-reading frames was also determined.

Ribosome Profiling and Mass Spectrometry Reveal Widespread Mitochondrial Translation Defects in a Striatal Cell Model of Huntington Disease.

Huntington disease (HD) is caused by an expanded polyglutamine mutation in huntingtin (mHTT) that promotes prominent atrophy in the striatum and subsequent psychiatric, cognitive deficits, and choreiform movements. Multiple lines of evidence point to an association between HD and aberrant striatal mitochondrial functions; however, the present knowledge about whether (or how) mitochondrial mRNA translation is differentially regulated in HD remains unclear. We found that protein synthesis is diminished in HD mitochondria compared to healthy control striatal cell models. We utilized ribosome profiling (Ribo-Seq) to analyze detailed snapshots of ribosome occupancy of the mitochondrial mRNA transcripts in control and HD striatal cell models. The Ribo-Seq data revealed almost unaltered ribosome occupancy on the nuclear-encoded mitochondrial transcripts involved in oxidative phosphorylation (SDHA, Ndufv1, Timm23, Tomm5, Mrps22) in HD cells. By contrast, ribosome occupancy was dramatically increased for mitochondrially encoded oxidative phosphorylation mRNAs (mt-Nd1, mt-Nd2, mt-Nd4, mt-Nd4l, mt-Nd5, mt-Nd6, mt-Co1, mt-Cytb, and mt-ATP8). We also applied tandem mass tag-based mass spectrometry identification of mitochondrial proteins to derive correlations between ribosome occupancy and actual mature mitochondrial protein products. We found many mitochondrial transcripts with comparable or higher ribosome occupancy, but diminished mitochondrial protein products, in HD. Thus, our study provides the first evidence of a widespread dichotomous effect on ribosome occupancy and protein abundance of mitochondria-related genes in HD.

Widespread stable noncanonical peptides identified by integrated analyses of ribosome profiling and ORF features.

Studies have revealed dozens of functional peptides in putative 'noncoding' regions and raised the question of how many proteins are encoded by noncanonical open reading frames (ORFs). Here, we comprehensively annotate genome-wide translated ORFs across five eukaryotes (human, mouse, zebrafish, worm, and yeast) by analyzing ribosome profiling data. We develop a logistic regression model named PepScore based on ORF features (expected length, encoded domain, and conservation) to calculate the probability that the encoded peptide is stable in humans. Systematic ectopic expression validates PepScore and shows that stable complex-associating microproteins can be encoded in 5'/3' untranslated regions and overlapping coding regions of mRNAs besides annotated noncoding RNAs. Stable noncanonical proteins follow conventional rules and localize to different subcellular compartments. Inhibition of proteasomal/lysosomal degradation pathways can stabilize some peptides especially those with moderate PepScores, but cannot rescue the expression of short ones with low PepScores suggesting they are directly degraded by cellular proteases. The majority of human noncanonical peptides with high PepScores show longer lengths but low conservation across species/mammals, and hundreds contain trait-associated genetic variants. Our study presents a statistical framework to identify stable noncanonical peptides in the genome and provides a valuable resource for functional characterization of noncanonical translation during development and disease.

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.

Ribosome Profiling Methods Adapted to the Study of RNA-Dependent Translation Regulation in Staphylococcus aureus.

Noncoding RNAs, including regulatory RNAs (sRNAs), are instrumental in regulating gene expression in pathogenic bacteria, allowing them to adapt to various stresses encountered in their host environments. Staphylococcus aureus is a well-studied model for RNA-mediated regulation of virulence and pathogenicity, with sRNAs playing significant roles in shaping S. aureus interactions with human and animal hosts. By modulating the translation and/or stability of target mRNAs, sRNAs regulate the synthesis of virulence factors and regulatory proteins required for pathogenesis. Moreover, perturbation of the levels of RNA modifications in two other classes of noncoding RNAs, rRNAs, and tRNAs, has been proposed to contribute to stress adaptation. However, the study of how these various factors affect translation regulation has often been restricted to specific genes, using in vivo reporters and/or in vitro translation systems. Genome-wide sequencing approaches offer novel perspectives for studying RNA-dependent regulation. In particular, ribosome profiling methods provide a powerful resource for characterizing the overall landscape of translational regulation, contributing to a better understanding of S. aureus physiopathology. Here, we describe protocols that we have adapted to perform ribosome profiling in S. aureus.

Utilizing high-resolution ribosome profiling for the global investigation of gene expression in Chlamydomonas.

Ribosome profiling (Ribo-seq) is a powerful method for the deep analysis of translation mechanisms and regulatory circuits during gene expression. Extraction and sequencing of ribosome-protected fragments (RPFs) and parallel RNA-seq yields genome-wide insight into translational dynamics and post-transcriptional control of gene expression. Here, we provide details on the Ribo-seq method and the subsequent analysis with the unicellular model alga Chlamydomonas reinhardtii (Chlamydomonas) for generating high-resolution data covering more than 10 000 different transcripts. Detailed analysis of the ribosomal offsets on transcripts uncovers presumable transition states during translocation of elongating ribosomes within the 5' and 3' sections of transcripts and characteristics of eukaryotic translation termination, which are fundamentally distinct for chloroplast translation. In chloroplasts, a heterogeneous RPF size distribution along the coding sequence indicates specific regulatory phases during protein synthesis. For example, local accumulation of small RPFs correlates with local slowdown of psbA translation, possibly uncovering an uncharacterized regulatory step during PsbA/D1 synthesis. Further analyses of RPF distribution along specific cytosolic transcripts revealed characteristic patterns of translation elongation exemplified for the major light-harvesting complex proteins, LHCs. By providing high-quality datasets for all subcellular genomes and attaching our data to the Chlamydomonas reference genome, we aim to make ribosome profiles easily accessible for the broad research community. The data can be browsed without advanced bioinformatic background knowledge for translation output levels of specific genes and their splice variants and for monitoring genome annotation.

Tumour circular RNAs elicit anti-tumour immunity by encoding cryptic peptides.

Emerging data have shown that previously defined noncoding genomes might encode peptides that bind human leukocyte antigen (HLA) as cryptic antigens to stimulate adaptive immunity1,2. However, the significance and mechanisms of action of cryptic antigens in anti-tumour immunity remain unclear. Here mass spectrometry of the HLA class I (HLA-I) peptidome coupled with ribosome sequencing of human breast cancer samples identified HLA-I-binding cryptic antigenic peptides that were noncanonically translated by a tumour-specific circular RNA (circRNA): circFAM53B. The cryptic peptides efficiently primed naive CD4+ and CD8+ T cells in an antigen-specific manner and induced anti-tumour immunity. Clinically, the expression of circFAM53B and its encoded peptides was associated with substantial infiltration of antigen-specific CD8+ T cells and better survival in patients with breast cancer and patients with melanoma. Mechanistically, circFAM53B-encoded peptides had strong binding affinity to both HLA-I and HLA-II molecules. In vivo, administration of vaccines consisting of tumour-specific circRNA or its encoded peptides in mice bearing breast cancer tumours or melanoma induced enhanced infiltration of tumour-antigen-specific cytotoxic T cells, which led to effective tumour control. Overall, our findings reveal that noncanonical translation of circRNAs can drive efficient anti-tumour immunity, which suggests that vaccination exploiting tumour-specific circRNAs may serve as an immunotherapeutic strategy against malignant tumours.

Improved super-resolution ribosome profiling reveals prevalent translation of upstream ORFs and small ORFs in Arabidopsis.

A crucial step in functional genomics is identifying actively translated ORFs and linking them to biological functions. The challenge lies in identifying short ORFs, as their identification is greatly influenced by data quality and depth. Here, we improved the coverage of super-resolution Ribo-seq in Arabidopsis (Arabidopsis thaliana), revealing uncharacterized translation events for nuclear, chloroplastic, and mitochondrial genes. Assisted by a transcriptome assembly, we identified 7,751 unconventional translation events, comprising 6,996 upstream ORFs (uORFs) and 209 downstream ORFs on annotated protein-coding genes, as well as 546 ORFs in presumed noncoding RNAs. Proteomic data confirmed the production of stable proteins from some of these unannotated translation events. We present evidence of active translation from primary transcripts of trans-acting small interfering RNAs (TAS1-4) and microRNAs (pri-MIR163 and pri-MIR169) and periodic ribosome stalling supporting cotranslational decay. Additionally, we developed a method for identifying extremely short uORFs, including 370 minimum uORFs (AUG-stop), and 2,921 tiny uORFs (2 to 10 amino acids) and 681 uORFs that overlap with each other. Remarkably, these short uORFs exhibit strong translational repression as do longer uORFs. We also systematically discovered 594 uORFs regulated by alternative splicing, suggesting widespread isoform-specific translational control. Finally, these prevalent uORFs are associated with numerous important pathways. In summary, our improved Arabidopsis translational landscape provides valuable resources to study gene expression regulation.

Decoding the Heterogeneity and Specialized Function of Translation Machinery Through Ribosome Profiling in Yeast Mutants of Initiation Factors.

The nuanced heterogeneity and specialized functions of translation machinery are increasingly recognized as crucial for precise translational regulation. Here, high-throughput ribosomal profiling (ribo-seq) is used to analyze the specialized roles of eukaryotic initiation factors (eIFs) in the budding yeast. By examining changes in ribosomal distribution across the genome resulting from knockouts of eIF4A, eIF4B, eIF4G1, CAF20, or EAP1, or knockdowns of eIF1, eIF1A, eIF4E, or PAB1, two distinct initiation-factor groups, the "looping" and "scanning" groups are discerned, based on similarities in the ribosomal landscapes their perturbation induced. The study delves into the cis-regulatory sequence features of genes influenced predominantly by each group, revealing that genes more dependent on the looping-group factors generally have shorter transcripts and poly(A) tails. In contrast, genes more dependent on the scanning-group factors often possess upstream open reading frames and exhibit a higher GC content in their 5' untranslated regions. From the ribosomal RNA fragments identified in the ribo-seq data, ribosomal heterogeneity associated with perturbation of specific initiation factors is further identified, suggesting their potential roles in regulating ribosomal components. Collectively, the study illuminates the complexity of translational regulation driven by heterogeneity and specialized functions of translation machinery, presenting potential approaches for targeted gene translation manipulation.

Ribosome Profiling of Plants.

Translation is a key step in control of gene expression, yet most analyses of global responses to a stimulus focus on transcription and the transcriptome. For RNA viruses in particular, which have no DNA-templated transcriptional control, control of viral and host translation is crucial. Here, we describe the method of ribosome profiling (ribo-seq) in plants, applied to virus infection. Ribo-seq is a deep sequencing technique that reveals the translatome by presenting a snapshot of the positions and relative amounts of translating ribosomes on all mRNAs in the cell. In contrast to RNA-seq, a crude cell extract is first digested with ribonuclease to degrade all mRNA not protected by a translating 80S ribosome. The resulting ribosome-protected fragments (RPFs) are deep sequenced. The number of reads mapping to a specific mRNA compared to the standard RNA-seq reads reveals the translational efficiency of that mRNA. Moreover, the precise positions of ribosome pause sites, previously unknown translatable open reading frames, and noncanonical translation events can be characterized quantitatively using ribo-seq. As this technique requires meticulous technique, here we present detailed step-by-step instructions for cell lysate preparation by flash freezing of samples, nuclease digestion of cell lysate, monosome collection by sucrose cushion ultracentrifugation, size-selective RNA extraction and rRNA depletion, library preparation for sequencing and finally quality control of sequenced data. These experimental methods apply to many plant systems, with minor nuclease digestion modifications depending on the plant tissue and species. This protocol should be valuable for studies of plant virus gene expression, and the global translational response to virus infection, or any other biotic or abiotic stress, by the host plant.

Ribosome profiling reveals the fine-tuned response of Escherichia coli to mild and severe acid stress.

IMPORTANCE: Bacteria react very differently to survive in acidic environments, such as the human gastrointestinal tract. Escherichia coli is one of the extremely acid-resistant bacteria and has a variety of acid-defense mechanisms. Here, we provide the first genome-wide overview of the adaptations of E. coli K-12 to mild and severe acid stress at both the transcriptional and translational levels. Using ribosome profiling and RNA sequencing, we uncover novel adaptations to different degrees of acidity, including previously hidden stress-induced small proteins and novel key transcription factors for acid defense, and report mRNAs with pH-dependent differential translation efficiency. In addition, we distinguish between acid-specific adaptations and general stress response mechanisms using denoising autoencoders. This workflow represents a powerful approach that takes advantage of next-generation sequencing techniques and machine learning to systematically analyze bacterial stress responses.

Discovering microproteins: making the most of ribosome profiling data.

Building a reference set of protein-coding open reading frames (ORFs) has revolutionized biological process discovery and understanding. Traditionally, gene models have been confirmed using cDNA sequencing and encoded translated regions inferred using sequence-based detection of start and stop combinations longer than 100 amino-acids to prevent false positives. This has led to small ORFs (smORFs) and their encoded proteins left un-annotated. Ribo-seq allows deciphering translated regions from untranslated irrespective of the length. In this review, we describe the power of Ribo-seq data in detection of smORFs while discussing the major challenge posed by data-quality, -depth and -sparseness in identifying the start and end of smORF translation. In particular, we outline smORF cataloguing efforts in humans and the large differences that have arisen due to variation in data, methods and assumptions. Although current versions of smORF reference sets can already be used as a powerful tool for hypothesis generation, we recommend that future editions should consider these data limitations and adopt unified processing for the community to establish a canonical catalogue of translated smORFs.

Streamlined and sensitive mono- and di-ribosome profiling in yeast and human cells.

Ribosome profiling has unveiled diverse regulation and perturbations of translation through a transcriptome-wide survey of ribosome occupancy, read out by sequencing of ribosome-protected messenger RNA fragments. Generation of ribosome footprints and their conversion into sequencing libraries is technically demanding and sensitive to biases that distort the representation of physiological ribosome occupancy. We address these challenges by producing ribosome footprints with P1 nuclease rather than RNase I and replacing RNA ligation with ordered two-template relay, a single-tube protocol for sequencing library preparation that incorporates adaptors by reverse transcription. Our streamlined approach reduced sequence bias and enhanced enrichment of ribosome footprints relative to ribosomal RNA. Furthermore, P1 nuclease preserved distinct juxtaposed ribosome complexes informative about yeast and human ribosome fates during translation initiation, stalling and termination. Our optimized methods for mRNA footprint generation and capture provide a richer translatome profile with low input and fewer technical challenges.