- Open Access
RiboMiner: a toolset for mining multi-dimensional features of the translatome with ribosome profiling data
BMC Bioinformatics volume 21, Article number: 340 (2020)
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.
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.
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.
Based on deep sequencing of the ribosome-protected mRNA fragments, ribosome profiling enables genome-wide investigations of translation with sub-codon resolution . In the past decade, ribosome profiling has been widely used for studies of translation under various contexts in human and almost all the major model organisms such as bacteria, C. elegans, yeast, and mouse, et al. [2, 3]. In these studies, the most popular applications of ribosome profiling data include quantifications of translation efficiency (TE) [4, 5], annotations of open reading frames (ORFs) [4, 6], meta-gene analysis of the ribosome distribution patterns [7, 8], and identification of translation initiation sites and pausing regions [3, 7].
A large variety of algorithms, software, and online resources have been developed to accommodate the ever-growing needs for data processing and analysis of ribosome profiling as well as presentation and interpretation of the results [9,10,11]. For example, RiboCode , ORF-RATER , Ribowave , ORFscore , RiboTaper , and Ribo-TISH  were designed for annotation of active ORFs. Xtail , anota , Riborex , Babel , and RiboDiff  were developed for quantification of TE changes. In addition to these highly focused analyses, more customized down-stream analyses of ribosome profiling data are of great value as well for better understanding of the translatome and the translation regulation. As summarized in Fig. 1, these pipelines implemented a broad range of data analyses, from pre-processing to down-stream information mining [17, 23,24,25,26,27,28,29,30,31,32,33,34,35,36]. For example, RiboPip is focused on the pre-processing of ribosome profiling and RNA-seq data; mQC  and Ribo-seQC  were mainly designed for quality control of the data; RiboPlot  and Shoelaces  can be used for quality control and visualization of the data; PausePred  was specifically designed to identify pausing motifs; RiboTools  is a Galaxy toolbox with functions including detection of translational ambiguities and identification of readthrough events; PROTEOFORMER  is a python package for automatic processing of ribosome profiling data, which includes SNP calling, ORF assembly and TIS identification; Plastid  is used for differential translation analysis and metagene analysis; RiboProfiling , riboSeqR , and systemPipeR  are R packages for tasks such as data pre-processing, quality control, TE calculation and ORF annotation, et al. These tools together cover many key procedures for ribosome profiling data analysis and presentation. However, as shown in Fig. 1, new types of analyses, many of which are related to the complicated machinery of translation regulation and only emerged in recent studies, have not been implemented by the current methods.
Here, we present RiboMiner, a python toolset for ribosome profiling data analysis, in response to the emerging needs for deeper mining of the hidden information about translation and its regulation, such as co-translation events, metagene analyses of the ribosome footprint density, and the potential regulatory factors of ribosome distribution, e.g., codon usage, tRNA gene copy numbers, and properties of nascent amino acids (Fig. 1). Characterization of these multi-dimensional features, under specific experimental or physiological conditions, provides closer views on the shift of translation landscapes with more details, which could be informative for pursuing the underlying machinery of translation regulation.
Overview of RiboMiner
The pipeline of RiboMiner is composed of 4 major function modules (Fig. 2): 1) Quality Control (QC), designed for assessment of ribosome profiling data quality with a multitude of benchmarks including 3-nt periodicity, distribution of read lengths, reads distribution in non-coding genome/transcriptome, etc. 2) Metagene Analysis (MA) for global distribution patterns of ribosome footprints on pooled transcripts of the full or any given subset of the transcriptome. This analysis helps identifying potential ribosome stalling events that take place at the global scale or just for subsets of the transcripts. Via direct comparisons between different conditions, metagene analysis of the footprints shed lights on how and at which stage these conditions perturb translation. 3) Feature Analysis (FA) for mining of various features that are enriched in predefined gene sets, e.g., the subset of genes with ribosome stalling events detected by MA as introduced above. The features being tested here include ribosome footprint densities on different amino acid (AA) and tri-amino acid (tri-AA) motifs, local and global codon adaptation index, local and global tRNA adaptation index, AA hydrophobicity and charge, etc. By identifying such features that are enriched in selected genes, this module of functions could provide insights into the machinery of translation landscape shift, from a multitude of different perspectives. 4) Enrichment Analysis (EA), designed for direct comparison of the ribosome footprint distributions across two experiments. For instance, this function can be used for data analysis of selective ribosome profiling to help the identification of co-translation regulation.
The python package of RiboMiner is freely available at https://github.com/xryanglab/RiboMiner and https://pypi.org/project/RiboMiner. We also offer a Docker image for RiboMiner at https://hub.docker.com/r/yanglab/ribocode_ribominer. The RiboMiner pipeline, including the testing data used in the present study, is also available as a Gene Container Service (GCS) on the Huawei Cloud.
Inputs of RiboMiner
As previously described, the common practice of raw data pre-processing before downstream analysis of ribosome profiling data includes quality control of sequencing, adapter trimming, read mapping and read counting (Fig. 2). The counts files could be used for differential expression analysis or quantification of differential translation efficiency with specialized tools such as DESeq2  or Xtail . The BAM files generated by read mapping can be used for TIS identification [23, 26] or ORF annotation [12, 13, 15] with specialized tools such as RiboCode  and PROTEOFORMER . There are three categories of input files for RiboMiner (Fig. 2): 1) genome sequences in FASTA format and an annotation file in GTF format, which would be used for annotations of the RNA transcripts and the protein coding sequences; 2) RPF sequences in FASTA format, after pre-processing such as adapter trimming and quality filtering; 3) Two BAM files generated by mapping of the RPF reads to the genome and the transcriptome. These BAM files can be supplied by the users or generated by the script we have offered, which can be found on the GitHub page of RiboMiner. See supplementary file for a detailed tutorial of RiboMiner, including preparation of input files and all the down-stream analyses.
Alternative splicing of the eukaryotic genes generates multiple transcript isoforms for each gene . To avoid ambiguous alignments of the sequencing reads, the longest transcript of each protein coding gene would be used for the following analyses. The function OutputTranscriptInfo performs this task and generates annotation files containing all the selected transcripts of the protein coding genes. GetProteinCodingSequence returns the sequences of these transcripts, the protein coding sequences, and the amino acid sequences, whereas GetUTRSequences extracts the UTR sequences of these transcripts specifically. It should be noted that all the functions above are dependent on the transcript annotation file generated by prepare_transcripts of RiboCode  that our group developed before.
Four functions were designed for quality control of the ribosome footprints. Periodicity, which is adapted from metaplots of RiboCode , is used for assessment of the 3-nt periodicity and identification of P-sites of the ribosome footprints. It reports the distributions of RPFs aligned by their 5′ end in relative to the start and stop codons, which is done for the reads of each specific length or for all the reads combined. RiboDensityOfDiffFrames returns the read densities of each reading frame. LengthDistribution or ReadsLengthOfSpecificRegions provides the length distributions of all the ribosome footprints or the ones from specific regions such as CDS, 5’UTR, and 3’UTR. StatisticReadsOnDNAsContam counts the RPF reads mapped to introns and intergenic regions of the genome, which are potentially DNA contaminations or other non-ribosome-footprint fragments and could be indicative to the data quality in general.
Metagene analysis aligns the transcripts of all the genes or a pre-defined gene set by their start codons and quantifies the relative read densities at each nucleotide or codon. Such analysis is particularly useful for identifying the potential global ribosome stalling sites under certain experimental or physiological conditions. MetageneAnalysisForTheWholeRegions can be used for calculating the read densities along the transcripts for an overall view, which would be helpful for testing whether the ribosome distributions are biasedly enriched. The function PolarityCalculation then evaluates the ribosome distribution bias for each gene and returns an overall distribution of such bias for all the genes. Finally, MetageneAnalysis was developed to zoom in and study the footprint densities in any particular region of the transcripts, including the UTR regions. This is particularly useful for allocating the ribosome stalling regions. All the results from the functions above can be readily presented as figures (see Supplementary file).
A number of factors have been shown to be involved in the regulation of translation initiation and elongation, such as poly-proline motifs [8, 39, 40], codon usage , tRNA gene copy numbers , amino acids with positive charges and high isoelectric points (pI) , etc. RiboMiner provides a series of functions for mining of such hidden features that are related to ribosome occupancy, with a goal of providing valuable insights into the molecular machinery of translation regulation. RiboDensityAtEachKindAAOrCodon calculates the ribosome footprint density at each amino acid (AA) or codon to show the differences under the experimental or physiological conditions. Furthermore, PausingScore was developed to quantify the ribosome density at each tri-AA motif and identify the motifs with enriched ribosome occupancy. Next, RiboDensityAroundTripleteAAMotifs can be used for computing the ribosome occupancy around the P/E site of the tri-AA motifs identified by the functions above.
tRNA adaptation index (tAI) and codon adaptation index (cAI) have been found to be potentially influential to translation, especially during elongation [41,42,43]. RiboMiner thereby provides two functions tAI and cAI to calculate the global tAI and cAI values as well as the local tAI and cAI values at each position along the transcripts for a specified organism. Note that although the weights of tAI are not exactly the same in all species, most of them are highly correlated to the weights fitted in yeast . Finally, hydropathyCharge calculates the hydrophobicity and charge indexes of amino acids encoded by each codon along the transcripts.
Previously, the ratio of RPF read counts in 5’UTR to CDS was reported to be negatively correlated with the translation efficiency (TE) . RiboMiner includes a function, RPFdist, to calculate these ratios, which would be potentially informative for evaluating the translation efficiencies in cases when the RNA-seq data in parallel with the ribosome profiling data is not available. See Supplementary file for more details about usages of these functions above.
Direct comparison between different ribosome profiling data could reveal the translation landscape shifts in details. For example, selective ribosome profiling (SeRP) is a powerful tool for studying the interaction of molecular chaperones and their potential targeting factors in the process of translation elongation . It also reveals the co-translation events among different subunits of protein complexes . Detailed analysis of the data from SeRP has been carefully done [46, 48], but these pipelines are not available in any of the current data analysis programs. Thus, RiboMiner incorporated a specially designed pipeline for mining of ribosome footprint enrichments with selective ribosome profiling data, which is quite different from the normal ribosome profiling data. This pipeline is composed of at least three steps: First, RiboDensityAtEachPosition calculates the ribosome density at each position for each transcript; Second, EnrichmentAnalysis performs an enrichment analysis by calculating the ribosome density ratio at each position in one ribosome profiling data over another one, e.g., ribosome profiling with IP of a specific protein and the normal ribosome profiling data. Finally, PlotEnrichmentRatio or EnrichmentAnalysisForSingleTrans generates plots of the results. In cases of multiple replicates for each condition, enrichmentMeanDensity can be used to calculate the mean ribosome density at each position ahead of the function EnrichmentAnalysis. It is worth noting that although this module was designed for SeRP data analysis, it could also be used for read density comparison between two normal ribosome profiling data under different conditions.
To showcase the application of RiboMiner for mining of insights into the translation regulation machineries, we used two published datasets, one from the study of eIF5A in translation regulation in yeast (GSE89704)  and the other from the study of co-translation of the aminoacyl-tRNA-synthetase complex (GSE116570) , which is specifically used for the Enrichment Analysis.
Benchmarks for the quality of a ribosome profiling dataset include 3-nt periodicity, appropriate distribution of the footprint lengths (usually ~ 28–30 nt), and enrichment in coding regions. As for GSE89704 (SRR5008135 for example), RiboMiner exhibited strong 3-nt periodicity for the reads with a specific length or for all the reads combined and aligned by their P-site positions (Fig. 3a). The lengths of ribosome footprints are usually 28 ~ 30 nt . However, the footprint length distribution of SRR5008135, generated by RiboMiner, showed that although the main peak was indeed around 28–30 nt, there is a small peak at 19 nt (Fig. 3b), which seems quite abnormal. This was further addressed by genome mapping of the RPF reads by RiboMiner, which showed that a small proportion of the reads were mapped to the intergenic regions and introns, suggesting some level of potential DNA contamination (Fig. 3b, c). After removal of these reads, the footprints strongly unified around 28–30 nt. Finally, frame analysis of the footprints showed that most reads were enriched in the first reading frame (Fig. 3d). Taken together, the results above indicate generally good quality of the data and identified the source of potential contamination, which could be easily eliminated by size selection of the raw reads.
RiboMiner was used for metagene analysis of the ribosome footprint distributions with a similar procedure as previously described . As shown in Fig. 4a, RiboMiner reproduced the dramatically changed pattern of ribosome occupancy upon knock-down of eIF5A, suggesting strong stalling at the early stage of elongation in the first 100 ~ 150 codons (Fig. 4a, c). Distributions of the polarity scores, generated by RiboMiner, confirmed the significant shift of the ribosome footprints towards the 5′ ends (Fig. 4b). In addition to the metagene analysis along the whole transcript or the CDS region, RiboMiner could also be used for metagene analysis for the UTR regions (Fig. 4d). Finally, it is worth noting that the metagene analysis with RiboMiner can be done at the global scale or just for a subset of the genes.
Mining of features related to translation regulation
Although the global metagene analysis as shown above indicates strong ribosome stalling upon eIF5A knock-down (Fig. 4a, c), it is possible that such pattern was mainly contributed by a subset of the genes. To specifically identify the genes with ribosome footprints enriched in the first 100 codons upon eIF5A knock-down, we compared the read densities on the first 100 codons for each transcript between the control and si-eIF5A samples. Two thousand nine hundred fifty-four genes showed up-regulated ribosome densities (“up-regulated genes”) with the ratio of si-eIF5A/si-Ctrl > = 1.5. This list of genes could then be used for mining of common features that are potentially related to the ribosome stalling during early elongation.
For example, RiboMiner can calculate ribosome density at each codon or amino acid (AA) for a set of genes (Fig. 5a, b) and further evaluate the change of ribosome density under different conditions (Fig. 5c, d). The results showed that read densities on proline (P) and aspartic acid (D) were largely increased after eIF5A-KD (si-eIF5A), suggesting that these amino acids may have contributed to the ribosome stalling observed above.
Next, RiboMiner was also used for calculating the ribosome densities on tri-AA motifs. Results showed that the motifs of poly-proline (PPP) and poly-aspartic acid (DDD) were significantly enriched by ribosome upon eIF5A knock-down (Fig. 6a), whereas the ribosome density on poly-lysine was decreased, but with no statistical significance (Fig. 6b). These results are consistent to the original reports in  (Fig. 6c). Besides, although the charges of amino acids have been reported to be negatively correlated to the speed of translation elongation , we found no difference in the charge of amino acids among the different gene sets (Fig. 6d).
Interestingly, analysis with RiboMiner revealed some more novel features that have not been reported in the original study. For example, we found that both the local and global cAI and tAI values of the genes with ribosome enrichment during early elongation (“up-regulated gene”) tend to be much smaller than those of the other transcripts used as control (“unblocked genes” and “down-regulated genes”) and so did the local cAI and tAI values (Fig. 7a, b, d). This suggests that the “up-regulated genes” have more sub-optimal codons, which then potentially led to slower moving of the ribosomes along the transcripts. In addition, the hydrophobicity of the amino acids encoded by these “up-regulated genes” are also much smaller, which we suspect may be related to ribosome stalling as well (Fig. 7c).
In summary, based on the results of RiboMiner, we reproduced the ribosome stalling event upon eIF5A knock-down as reported by the original study. Such stalling was mainly contributed by a subset of genes, and further mining of multiple features of these genes indicated that tri-AA motifs, codon usage, tRNA copy numbers, and amino acid hydrophobicity may be related to the translation dysregulations due to eIF5A knock-down. Many of these new insights have not been reported by the original study and therefore could worth further investigations.
Enrichment analysis with selective ribosome profiling data
We used RiboMiner to revisit the selective ribosome profiling data in the study of the assembly of a hetero-trimeric complex, the multi-aminoacyl-tRNA synthetase . This complex is composed of three major subunits, i.e., the essential methionyl- and glutamyl-tRNA synthetases MetRS and GluRS (encoded by MES1 and GUS1, respectively) and the Arc1p cofactor (encoded by ARC1) regulating the catalytic activities and subcellular distributions of the complex . RiboMiner regenerated the distributions of ribosome footprints that are enriched in the selective ribosome profiling (Fig. 8). This confirms the main conclusion of the original study, i.e., proteins MetRS and GluRS could co-translate with each other, and both participate in the translation of Arc1p starting at a specific position .
Running time of different modules of RiboMiner
Last, the running time of each module of RiboMiner with the data of 4 different model organisms are provided in Table 1. We did not consider the pre-filtering and data preparation steps. In general, the time mainly depends on the total reads, i.e., the sequencing depth. All jobs were run on a Linux server with Intel R Xeon R CPUs at 2.40GHz, 64 G memory. Note that the table reports the total running time of all the functions in each module. The scripts and results for evaluation of the running time are available at https://github.com/xryanglab/RiboMiner/tree/master/RuningTimeTest.
RiboMiner is a python toolset for deep mining of multi-dimensional features of the translatomes with ribosome profiling data. A multitude of functions incorporated in RiboMiner are useful for quality control of ribosome profiling data, metagene analysis for detection of translation dysregulations such as ribosome stalling, mining of various features related to the translation dysregulations, and exploration of selective ribosome profiling for fine maps of translation regulation such as co-translation. Applications of RiboMiner on two published datasets did not only reproduce the main results reported before, but also generated novel insights into the translation regulation processes. RiboMiner provides quantitative data for visualization as well as statistical analyses. In summary, here we present RiboMiner as a complementary toolset to the current methods, to facilitate the comprehensive and thorough dissections of the translation landscapes as well as the translation regulations with the technique of ribosome profiling.
Availability and requirements
Project name: RiboMiner
Project home page: https://github.com/xryanglab/RiboMiner
Operating system(s): Platform independent
Programming language: Python
Other requirements: Python (version > = 3.6)
License: GNU GPL
Any restrictions to use by non-academics: Licence needed
Availability of data and materials
RiboMiner is freely available at https://github.com/xryanglab/RiboMiner and https://pypi.org/project/RiboMiner. Detailed instructions and the information of all the sequencing data used in the current manuscript are available at https://github.com/xryanglab/RiboMiner/blob/master/Implementation.md. We also offer a Docker image for RiboMiner at https://hub.docker.com/r/yanglab/ribocode_ribominer. The RiboMiner pipeline, including the testing data used in the present study, is also available as a Gene Container Service (GCS) on the Huawei Cloud.
Open reading frames
Translation initiation site
Ribosome profiling footprint
tRNA adaptation index
Codon adaptation index
Selective ribosome profiling
Ingolia NT, Ghaemmaghami S, Newman JRS, Weissman JS. Genome-wide analysis in vivo of translation with nucleotide resolution using ribosome profiling. Science. 2009;324(5924):218–23.
Ingolia NT. Ribosome profiling: new views of translation, from single codons to genome scale. Nat Rev Genet. 2014;15(3):205–13.
Ingolia NT, Hussmann JA, Weissman JS. Ribosome Profiling: Global Views of Translation. CSH PERSPECT BIOL. 2019;11(5):a032698. https://doi.org/10.1101/cshperspect.a032698.
Zou Q, Xiao Z, Huang R, Wang X, Wang X, Zhao H, Yang X. Survey of the translation shifts in hepatocellular carcinoma with ribosome profiling. Theranostics. 2019;9(14):4141–55.
Lin Z, Hsu PJ, Xing X, Fang J, Lu Z, Zou Q, Zhang KJ, Zhang X, Zhou Y, Zhang T, et al. Mettl3−/Mettl14-mediated mRNA N (6)-methyladenosine modulates murine spermatogenesis. Cell Res. 2017;27(10):1216–30.
Jackson R, Kroehling L, Khitun A, Bailis W, Jarret A, York AG, Khan OM, Brewer JR, Skadow MH, Duizer C, et al. The translation of non-canonical open reading frames controls mucosal immunity. Nature. 2018;564(7736):434–8.
Lin Y, Li F, Huang L, Polte C, Duan H, Fang J, Sun L, Xing X, Tian G, Cheng Y, et al. eIF3 associates with 80S ribosomes to promote translation elongation, mitochondrial homeostasis, and muscle health. Mol Cell. 2020. https://doi.org/10.1016/j.molcel.2020.06.003.
Schuller AP, Wu CC, Dever TE, Buskirk AR, Green R. eIF5A functions globally in translation elongation and termination. Mol Cell. 2017;66(2):194–205.e195.
Wang H, Wang Y, Xie Z. Computational resources for ribosome profiling: from database to web server and software. Brief Bioinform. 2019;20(1):144–55.
Andreev DE, O'Connor PBF, Loughran G, Dmitriev SE, Baranov PV, Shatsky IN. Insights into the mechanisms of eukaryotic translation gained with ribosome profiling. Nucleic Acids Res. 2017;45(2):513–26.
Wang Y, Zhang H, Lu J. Recent advances in ribosome profiling for deciphering translational regulation. Meth. 2020;176:46–54. https://doi.org/10.1016/j.ymeth.2019.05.011.
Xiao Z, Huang R, Xing X, Chen Y, Deng H, Yang X. De novo annotation and characterization of the translatome with ribosome profiling data. Nucleic Acids Res. 2018;46(10):e61.
Fields AP, Rodriguez EH, Jovanovic M. A regression-based analysis of ribosome-profiling data reveals a conserved complexity to mammalian translation. Mol Cell. 2015;60(5):11.
Xu ZY, Hu L, Shi BB, Geng SS, Xu LC, Wang D, Lu ZJ. Ribosome elongating footprints denoised by wavelet transform comprehensively characterize dynamic cellular translation events. Nucleic Acids Res. 2018;46(18):e109.
Bazzini AA, Johnstone TG, Christiano R, Mackowiak SD, Obermayer B, Fleming ES, Vejnar CE, Lee MT, Rajewsky N, Walther TC, et al. Identification of small ORFs in vertebrates using ribosome footprinting and evolutionary conservation. EMBO J. 2014;33(9):981–93.
Calviello L, Mukherjee N, Wyler E, Zauber H, Hirsekorn A, Selbach M, Landthaler M, Obermayer B, Ohler U. Detecting actively translated open reading frames in ribosome profiling data. Nat Methods. 2016;13(2):165–70.
Zhang P, He DD, Xu Y, Hou JK, Pan BF, Wang YF, Liu T, Davis CM, Ehli EA, Tan L, et al. Genome-wide identification and differential analysis of translational initiation. Nat Commun. 2017;8(1):1749.
Xiao Z, Zou Q, Liu Y, Yang X. Genome-wide assessment of differential translations with ribosome profiling data. Nat Commun. 2016;7:11194.
Larsson O, Sonenberg N, Nadon R. anota: analysis of differential translation in genome-wide studies. Bioinformatics. 2011;27(10):1440–1.
Li W, Wang W, Uren PJ, Penalva LOF, Smith AD. Riborex: fast and flexible identification of differential translation from Ribo-seq data. Bioinformatics. 2017;33(11):1735–7.
Olshen AB, Hsieh AC, Stumpf CR, Olshen RA, Ruggero D, Taylor BS. Assessing gene-level translational control from ribosome profiling. Bioinformatics. 2013;29(23):2995–3002.
Zhong Y, Karaletsos T, Drewe P, Sreedharan VT, Kuo D, Singh K, Wendel HG, Ratsch G. RiboDiff: detecting changes of mRNA translation efficiency from ribosome footprints. Bioinformatics. 2017;33(1):139–41.
Crappe J, Ndah E, Koch A, Steyaert S, Gawron D, De Keulenaer S, De Meester E, De Meyer T, Van Criekinge W, Van Damme P, et al. PROTEOFORMER: deep proteome coverage through ribosome profiling and MS integration. Nucleic Acids Res. 2015;43(5):e29.
Legendre R, Baudin-Baillieu A, Hatin I, Namy O. RiboTools: a galaxy toolbox for qualitative ribosome profiling analysis. Bioinformatics. 2015;31(15):2586–8.
Chung BY, Hardcastle TJ, Jones JD, Irigoyen N, Firth AE, Baulcombe DC, Brierley I. The use of duplex-specific nuclease in ribosome profiling and a user-friendly software package for Ribo-seq data analysis. RNA. 2015;21(10):1731–45.
de Klerk E, Fokkema IF, Thiadens KA, Goeman JJ, Palmblad M, den Dunnen JT, von Lindern M, ‘t Hoen PA. Assessing the translational landscape of myogenic differentiation by ribosome profiling. Nucleic Acids Res. 2015;43(9):4408–28.
Dunn JG, Weissman JS. Plastid: nucleotide-resolution analysis of next-generation sequencing and genomics data. BMC Genomics. 2016;17(1):958.
Michel AM, Mullan JPA, Velayudhan V, O'Connor PBF, Donohue CA, Baranov PV. RiboGalaxy: a browser based platform for the alignment, analysis and visualization of ribosome profiling data. RNA Biol. 2016;13(3):316–9.
Backman TWH, Girke T. systemPipeR: NGS workflow and report generation environment. BMC Bioinformatics. 2016;17:388.
Popa A, Lebrigand K, Paquet A, Nottet N, Robbe-Sermesant K, Waldmann R, Barbry P. RiboProfiling: a Bioconductor package for standard Ribo-seq pipeline processing. F1000Res. 2016;5:1309.
Wang H, Kingsford C, McManus CJ. Using the Ribodeblur pipeline to recover A-sites from yeast ribosome profiling data. Methods. 2018;137:67–70.
Lauria F, Tebaldi T, Bernabo P, Groen EJN, Gillingwater TH, Viero G. riboWaltz: optimization of ribosome P-site positioning in ribosome profiling data. PLoS Comput Biol. 2018;14(8):e1006169.
Kumari R, Michel AM, Baranov PV. PausePred and Rfeet: webtools for inferring ribosome pauses and visualizing footprint density from ribosome profiling data. RNA. 2018;24(10):1297–304.
Birkeland A, Chyzynska K, Valen E. Shoelaces: an interactive tool for ribosome profiling processing and visualization. BMC Genomics. 2018;19(1):543.
Verbruggen S, Menschaert G. mQC: A post-mapping data exploration tool for ribosome profiling. Comput Methods Programs Biomed. 2019;181:104806. https://doi.org/10.1016/j.cmpb.2018.10.018.
Calviello L, Sydow D, Harnett D, Ohler U. Ribo-seQC: Comprehensive analysis of cytoplasmic and organellar ribosome profiling data. Cold Spring Harbor Laboratory Press; 2019. https://doi.org/10.1101/601468.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.
Black DL. Mechanisms of alternative pre-messenger RNA splicing. Annu Rev Biochem. 2003;72:291–336.
Doerfel LK, Wohlgemuth I, Kothe C, Peske F, Urlaub H, Rodnina MV. EF-P is essential for rapid synthesis of proteins containing consecutive Proline residues. Science. 2013;339(6115):85–8.
Woolstenhulme CJ, Guydosh NR, Green R, Buskirk AR. High-precision analysis of translational pausing by ribosome profiling in Bacteria lacking EFP. Cell Rep. 2015;11(1):13–21.
Zhang S, Hu H, Zhou J, He X, Jiang T, Zeng J. Analysis of ribosome stalling and translation elongation dynamics by deep learning. Cell Syst. 2017;5(3):212–20 e216.
Tuller T, Carmi A, Vestsigian K, Navon S, Dorfan Y, Zaborske J, Pan T, Dahan O, Furman I, Pilpel Y. An evolutionarily conserved mechanism for controlling the efficiency of protein translation. Cell. 2010;141(2):344–54.
Riba A, Di Nanni N, Mittal N, Arhne E, Schmidt A, Zavolan M. Protein synthesis rates and ribosome occupancies reveal determinants of translation elongation rates. Proc Natl Acad Sci U S A. 2019;116(30):15023–32.
Sabi R, Daniel RV, Tuller T. stAI (calc): tRNA adaptation index calculator based on species-specific weights. Bioinformatics. 2017;33(4):589–91.
Murat P, Marsico G, Herdy B, Ghanbarian AT, Portella G, Balasubramanian S. RNA G-quadruplexes at upstream open reading frames cause DHX36- and DHX9-dependent translation of human mRNAs. Genome Biol. 2018;19(1):229.
Becker AH, Oh E, Weissman JS, Kramer G, Bukau B. Selective ribosome profiling as a tool for studying the interaction of chaperones and targeting factors with nascent polypeptide chains and ribosomes. Nat Protoc. 2013;8(11):2212–39.
Shiber A, Doring K, Friedrich U, Klann K, Merker D, Zedan M, Tippmann F, Kramer G, Bukau B. Cotranslational assembly of protein complexes in eukaryotes revealed by ribosome profiling. Nature. 2018;561(7722):268–72.
Galmozzi CV, Merker D, Friedrich UA, Döring K, Kramer G. Selective ribosome profiling to study interactions of translating ribosomes in yeast. Nat Protoc. 2019;14(8):2279–317.
The authors wish to acknowledge the supports from the core facilities of Genome Sequencing & Analysis and High-Performance Computing of the National Protein Science Facility (Beijing).
This work was supported by the National key research and development program, Precision Medicine Project (2016YFC0906001), the Tsinghua University Spring Breeze Fund, the Tsinghua University Initiative Scientific Research Program (2019Z06QCX01), and the National Natural Science Foundation of China (81972912 and 31671381). The funding bodies did not play any roles in the design of the study, in the collection, analysis, interpretation of data, or in writing of the manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Li, F., Xing, X., Xiao, Z. et al. RiboMiner: a toolset for mining multi-dimensional features of the translatome with ribosome profiling data. BMC Bioinformatics 21, 340 (2020). https://doi.org/10.1186/s12859-020-03670-8