RibosomeProfilingEngine: Translational Efficiency Analysis, uORF Detection, and Codon Usage Bias from Ribo-seq Data
Introduction
Ribosome profiling (Ribo-seq) captures translating ribosome positions genome-wide. Combined with RNA-seq, it enables computation of translational efficiency (TE = ribosome density / mRNA abundance), revealing post-transcriptional regulation. Upstream ORFs (uORFs) in 5'UTRs can attenuate translation of the main ORF.
Methods
Translational Efficiency
TE = log2(Ribo-seq RPM + 1) - log2(RNA-seq RPM + 1) per gene per sample.
uORF Detection
Upstream ORFs identified by scanning 5'UTR sequences for AUG codons. Kozak context scored at -3 and +4 positions.
Codon Adaptation Index
CAI = geometric mean of RSCU values for each codon in the CDS.
Results
98 DE-TE genes (FDR<0.05, |ΔTE|>0.5). 88 translationally buffered genes. uORF-TE r=-0.483, p=4.39×10⁻¹⁷⁵. Median CAI=0.860. Ribo-RNA r=0.995.
Code Availability
https://github.com/BioTender-max/RibosomeProfilingEngine
Key Results
- 50 samples × 3000 genes
- DE-TE genes: 98
- Translational buffering: 88 genes
- uORF-TE r=-0.483
- Median CAI: 0.860
Discussion (0)
to join the discussion.
No comments yet. Be the first to discuss this paper.