Dataset Information

Ribo-ODDR: Oligo design pipeline for experiment-specific rRNA depletion in ribo-seq.

ABSTRACT:

Motivation

Ribosome Profiling (Ribo-seq) has revolutionized the study of RNA translation by providing information on ribosome positions across all translated RNAs with nucleotide-resolution. Yet several technical limitations restrict the sequencing depth of such experiments, the most common of which is the overabundance of rRNA fragments. Various strategies can be employed to tackle this issue, including the use of commercial rRNA depletion kits. However, as they are designed for more standardized RNAseq experiments, they may perform suboptimally in Ribo-seq. In order to overcome this, it is possible to use custom biotinylated oligos complementary to the most abundant rRNA fragments, however currently no computational framework exists to aid the design of optimal oligos.

Results

Here, we first show that a major confounding issue is that the rRNA fragments generated via Ribo-seq vary significantly with differing experimental conditions, suggesting that a "one-size-fits-all" approach may be inefficient. Therefore we developed Ribo-ODDR, an oligo design pipeline integrated with a user-friendly interface that assists in oligo selection for efficient experiment-specific rRNA depletion. Ribo-ODDR uses preliminary data to identify the most abundant rRNA fragments, and calculates the rRNA depletion efficiency of potential oligos. We experimentally show that Ribo-ODDR designed oligos outperform commercially available kits and lead to a significant increase in rRNA depletion in Ribo-seq.

Availability

Ribo-ODDR is freely accessible at https://github.com/fallerlab/Ribo-ODDR.

Supplementary information

Supplementary data are available at Bioinformatics online.

SUBMITTER: Alkan F

PROVIDER: S-EPMC8428583 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Ribosome profiling (Ribo-Seq) is a highly sensitive method to quantify ribosome occupancies along individual mRNAs on a genome-wide scale. Hereby, ribosome-protected fragments (= footprints) are generated by nuclease digestion, isolated, and sequenced together with the corresponding randomly fragmented input samples, to determine ribosome densities (RD). For library preparation, equal amounts of total RNA are used. Subsequently, all transcript fragments are subjected to linker ligation, cDNA synthesis, and PCR amplification. Importantly, the number of reads obtained for every transcript in input and footprint samples during sequencing depends on sequencing depth and library size, as well as the relative abundance of the transcript in the sample. However, the information pertaining to the absolute amount of input and footprint sequences is lost during sample preparation, hence ruling out any conclusion whether translation is generally suppressed or activated in one condition over the other. Therefore, the RD fold-changes determined for individual genes do not reflect absolute regulation, but have to be interpreted as relative to bulk mRNA translation. Here, we modified the original ribosome profiling protocol that was first established by Ingolia et al. (2009), by adding small amounts of yeast lysate to the mammalian lysates of interest as a spike-in. This allows us to not only detect changes in the RD of specific transcripts relative to each other, but also to simultaneously measure global differences in RD (normalized ribosome density values) between samples. Graphic abstract: Global changes in translation efficiency can be detected with polysome profiling, where the proportion of polysomal ribosomes is interpreted as a proxy for ribosome density (RD) on bulk mRNA. Ribo-Seq measures changes in RD of specific mRNAs relative to bulk mRNA. The addition of a yeast-lysate, as a spike-in for normalization of read counts, allows for an absolute measurement of changes in RD.

Project description:Ribosome profiling (Ribo-Seq) is a powerful method to study translation in bacteria. However, Ribo-Seq signal can be observed across RNAs that one would not expect to be bound by ribosomes. For example, Escherichia coli Ribo-Seq libraries also capture reads from most noncoding RNAs (ncRNAs). While some of these ncRNAs may overlap coding regions, this alone does not explain the majority of observed signal across ncRNAs. These fragments of ncRNAs in Ribo-Seq data pass all size selection steps of the Ribo-Seq protocol and survive hours of micrococcal nuclease (MNase) treatment. In this work, we specifically focus on Ribo-Seq signal across ncRNAs and provide evidence to suggest that RNA structure, as opposed to ribosome binding, protects them from degradation and allows them to persist in the Ribo-Seq sequencing library preparation. By inspecting these "contaminant reads" in bacterial Ribo-Seq, we show that data previously disregarded in bacterial Ribo-Seq experiments may, in fact, be used to gain partial information regarding the in vivo secondary structure of ncRNAs.IMPORTANCE Structured ncRNAs are pivotal mediators of bioregulation in bacteria, and their functions are often reliant on their specific structures. Here, we first inspect Ribo-Seq reads across noncoding regions, identifying contaminant reads in these libraries. We observe that contaminant reads in bacterial Ribo-Seq experiments that are often disregarded, in fact, strongly overlap with structured regions of ncRNAs. We then perform several bioinformatic analyses to determine why these contaminant reads may persist in Ribo-Seq libraries. Finally, we highlight some structured RNA contaminants in Ribo-Seq and support the hypothesis that structures in the RNA protect them from MNase digestion. We conclude that researchers should be cautious when interpreting Ribo-Seq signal as coding without considering signal distribution. These findings also may enable us to partially resolve RNA structures, identify novel structured RNAs, and elucidate RNA structure-function relationships in bacteria at a large scale and in vivo through the reanalysis of existing Ribo-Seq data sets.