Skip to content

4.8 smORF analysis

Function

The smORF analysis module provides a complete workflow for scanning candidate small open reading frames (smORFs), filtering candidates with sequence and Kozak-context criteria, evaluating Ribo-seq translation evidence, and integrating sample-level evidence into ORF-level summary tables.

The workflow is designed for transcript-centric smORF discovery. Candidate ORFs are first generated from a genome FASTA and genePred annotation, then filtered and evaluated using P-site density tracks from Ribo-seq data.

Workflow

smorf_scanner → smorf_filter → smorf_evidence → smorf_integrate
Step Command Main purpose
1 smorf_scanner Scan transcript-centric candidate ORFs from genome FASTA and genePred annotation.
2 smorf_filter Filter scanned ORFs by length, start codon, category, strand, completeness, and Kozak context.
3 smorf_evidence Evaluate Ribo-seq P-site density, frame periodicity, coverage, pausing, and release evidence.
4 smorf_integrate Integrate long-format sample-level evidence into ORF-level matrices and summary tables.

Input overview

Input Used by Description
Genome FASTA smorf_scanner Genome sequence used to reconstruct transcript sequences.
genePred annotation smorf_scanner, optional for smorf_evidence Transcript annotation. Scanner uses it to define transcript structures; evidence uses it to recover ORF exon blocks when needed.
Scanner message table smorf_filter Candidate ORF table generated by smorf_scanner.
Filtered ORF table smorf_evidence Passed ORF table generated by smorf_filter.
P-site density files smorf_evidence Strand-specific or unstranded bedGraph/WIG files containing Ribo-seq P-site density.
Evidence table smorf_integrate Long-format ORF-by-sample table generated by smorf_evidence.

Main evidence levels

Evidence Meaning
ORF sequence Start codon, stop codon, ORF length, strand, and category.
Kozak context Start-codon context scored by annotated, built-in, PWM, or sequence-derived Kozak models.
Ribo-seq signal Total RPF signal, covered nucleotides/codons, and coverage ratio.
Periodicity Frame-specific signal distribution, especially frame-0 enrichment.
Start/stop behavior Start pausing, pre-stop pausing, and post-stop release signatures.
Coverage shape Uniform, disperse, or skewed RPF distribution across the ORF.
Multi-sample support Reproducibility and support level across multiple Ribo-seq samples.
Integration ORF-level summary, best sample, evidence labels, and frame-density matrix.

Notes

  • smorf_scanner produces many candidate ORFs. Scanner output alone should not be treated as evidence of translation.
  • smorf_filter removes candidates according to sequence features and Kozak context, but it does not use Ribo-seq evidence.
  • smorf_evidence requires either a density list, one unstranded density file, or both plus- and minus-strand density files.
  • smorf_integrate is most useful when multiple Ribo-seq samples or replicates are available.
  • Strand-specific P-site density files can be generated with rpf_Bam2bw and then supplied to smorf_evidence as bedGraph or WIG tracks.