skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: RiboDiPA: a novel tool for differential pattern analysis in Ribo-seq data
Abstract 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.  more » « less
Award ID(s):
1764421
PAR ID:
10251805
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Nucleic Acids Research
Volume:
48
Issue:
21
ISSN:
0305-1048
Page Range / eLocation ID:
12016 to 12029
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Ellermeier, Craig D (Ed.)
    ABSTRACT Oxidative stress induces a wide range of cellular damage, often causing disease and cell death. While many organisms are susceptible to the effects of oxidative stress, haloarchaea have adapted to be highly resistant. Several aspects of the haloarchaeal oxidative stress response have been characterized; however, little is known about the impacts of oxidative stress at the translation level. Using the model archaeonHaloferax volcanii, we performed RNA-seq and ribosome profiling (Ribo-seq) to characterize the global translation landscape during oxidative stress. We identified 281 genes with differential translation efficiency (TE). Downregulated genes were enriched in ribosomal and translation proteins, in addition to peroxidases and genes involved in the TCA cycle. We also identified 42 small noncoding RNAs (sRNAs) with ribosome occupancy. Size distributions of ribosome footprints revealed distinct patterns for coding and noncoding genes, with 12 sRNAs matching the pattern of coding genes, and mass spectrometry confirming the presence of seven small proteins encoded by these sRNAs. However, the majority of sRNAs with ribosome occupancy had no evidence of coding potential. Of these ribosome-associated sRNAs, 12 had differential ribosome occupancy or TE during oxidative stress, suggesting that they may play a regulatory role during the oxidative stress response. Our findings on ribosomal regulation during oxidative stress, coupled with potential roles for ribosome-associated noncoding sRNAs and sRNA-derived small proteins inH. volcanii, revealed additional regulatory layers and underscored the multifaceted architecture of stress-responsive regulatory networks.IMPORTANCEArchaea are found in diverse environments, including as members of the human microbiome, and are known to play essential ecological roles in major geochemical cycles. The study of archaeal biology has expanded our understanding of the evolution of eukaryotes, uncovered novel biological systems, and revealed new opportunities for applications in biotechnology and bioremediation. Many archaeal systems, however, remain poorly characterized. UsingHaloferax volcaniias a model, we investigated the global translation landscape during oxidative stress. Our findings expand current knowledge of translational regulation in archaea and further illustrate the complexity of stress-responsive gene regulation. 
    more » « less
  2. Abstract BackgroundRibosome profiling, also known as Ribo-seq, is a powerful technique to study genome-wide mRNA translation. It reveals the precise positions and quantification of ribosomes on mRNAs through deep sequencing of ribosome footprints. We previously optimized the resolution of this technique in plants. However, several key reagents in our original method have been discontinued, and thus, there is an urgent need to establish an alternative protocol. ResultsHere we describe a step-by-step protocol that combines our optimized ribosome footprinting in plants with available custom library construction methods established in yeast and bacteria. We tested this protocol in 7-day-old Arabidopsis seedlings and evaluated the quality of the sequencing data regarding ribosome footprint length, mapped genomic features, and the periodic properties corresponding to actively translating ribosomes through open resource bioinformatic tools. We successfully generated high-quality Ribo-seq data comparable with our original method. ConclusionsWe established a custom library construction method for super-resolution Ribo-seq in Arabidopsis. The experimental protocol and bioinformatic pipeline should be readily applicable to other plant tissues and species. 
    more » « less
  3. null (Ed.)
    Abstract Background Translation is a fundamental process in gene expression. Ribosome profiling is a method that enables the study of transcriptome-wide translation. A fundamental, technical challenge in analyzing Ribo-Seq data is identifying the A-site location on ribosome-protected mRNA fragments. Identification of the A-site is essential as it is at this location on the ribosome where a codon is translated into an amino acid. Incorrect assignment of a read to the A-site can lead to lower signal-to-noise ratio and loss of correlations necessary to understand the molecular factors influencing translation. Therefore, an easy-to-use and accurate analysis tool is needed to accurately identify the A-site locations. Results We present RiboA, a web application that identifies the most accurate A-site location on a ribosome-protected mRNA fragment and generates the A-site read density profiles. It uses an Integer Programming method that reflects the biological fact that the A-site of actively translating ribosomes is generally located between the second codon and stop codon of a transcript, and utilizes a wide range of mRNA fragment sizes in and around the coding sequence (CDS). The web application is containerized with Docker, and it can be easily ported across platforms. Conclusions The Integer Programming method that RiboA utilizes is the most accurate in identifying the A-site on Ribo-Seq mRNA fragments compared to other methods. RiboA makes it easier for the community to use this method via a user-friendly and portable web application. In addition, RiboA supports reproducible analyses by tracking all the input datasets and parameters, and it provides enhanced visualization to facilitate scientific exploration. RiboA is available as a web service at https://a-site.vmhost.psu.edu/ . The code is publicly available at https://github.com/obrien-lab/aip_web_docker under the MIT license. 
    more » « less
  4. Visualizing Ribo-seq and other sequencing data within genes of interest is a powerful approach to studying gene expression, but its application is limited by a lack of robust tools. Here, we introduce ggRibo, a user-friendly R package for visualizing individual gene expression, integrating Ribo-seq, RNA-seq, and other genome-wide datasets with flexible scaling options. ggRibo visualizes 3-nucleotide periodicity, a hallmark of translating ribosomes, within a gene-structure context, including introns and untranslated regions, enabling the study of novel ORFs, translation of different isoforms, and mechanisms of translational regulation. ggRibo can plot multiple Ribo-seq/RNA-seq datasets from different conditions for comparison. It also contains functions for plotting single-transcript view, reading-frame decomposition, and RNA-seq coverage alone. Importantly, ggRibo supports the visualization of other omics datasets that could also be presented with single-nucleotide resolution, such as RNA degradome, transcription start sites, translation initiation sites, and epitranscriptomic modifications. We demonstrate its utility with examples of upstream ORFs, downstream ORFs, nested ORFs, and differential isoform translation in humans,Arabidopsis, tomato, and rice. We also provide examples of multiomic comparisons that reveal insights that connect the transcriptome, translatome, and degradome. In summary, ggRibo is an advanced single-gene viewer that offers a valuable resource for studying gene expression regulation through its intuitive and flexible platform. 
    more » « less
  5. Abstract A crucial step in functional genomics is identifying actively translated open reading frames (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 non-coding RNAs. Proteomics data confirmed the production of stable proteins from some of these unannotated translation events. We present evidence of active translation from primary transcripts of tasiRNAs (TAS1–4) and microRNAs (pri-MIR163, pri-MIR169), and periodic ribosome stalling supporting co-translational decay. Additionally, we developed a method for identifying extremely short uORFs, including 370 minimum uORFs (AUG-stop), and 2,921 tiny uORFs (2–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. 
    more » « less