orfipy 0.0.3

orfipy

Introduction

orfipy is a tool written in python/cython to extract ORFs in extremely an fast and flexible manner. Other popular ORF searching tools are OrfM and getorf. Compared to OrfM and getorf, orfipy provides the most options to fine tune ORF searches. orfipy uses multiple CPU cores and is particularly faster for data containing multiple smaller fasta sequences such as de-novo transcriptome assemblies. Please read the preprint here.

Installation

Install latest stable version

pip install orfipy

Or install via conda

conda install -c bioconda orfipy

Install the development version from source

git clone https://github.com/urmi-21/orfipy.git
cd orfipy
pip install .

or use pip

pip install git+git://github.com/urmi-21/orfipy.git

Examples

Details of orfipy algorithm are in the preprint and SI. Please go through the SI if you are interested to know differences between orfipy and other ORF finder tools and how to set orfipy parameters to match the output of other tools.

Below are some usage examples for orfipy

To see full list of options use the command:

orfipy -h

Input

orfipy, version 0.0.3 and above, supports sequences in Fasta/Fastq format (orfipy uses pyfastx). Input files can be in .gz format.

Extract ORF sequences and write ORF sequences in orfs.fa file

orfipy input.fasta --dna orfs.fa --min 10 --max 10000 --procs 4 --table 1 --outdir orfs_out

Use standard codon table but use only ATG as start codon

orfipy input.fa.gz --dna orfs.fa --start ATG

Note: Users can also provide their own translation table, as a .json file, to orfipy using --table option. Example of json file containing a valid translation table is here

See available codon tables

orfipy --show-table

Extract ORFs BED file

orfipy input.fasta --bed orfs.bed --min 50 --procs 4
or
orfipy input.fasta --min 50 --procs 4 > orfs.bed 

Extract ORFs BED12 file

Note: Add --include-stop for orfipy output to be consistent with Transdecoder.Predict output .bed file.

orfipy testseq.fa --min 100 --bed12 of.bed --partial-5 --partial-3 --include-stop

Extract ORFs peptide sequences using default translation table

orfipy input.fasta --pep orfs_peptides.fa --min 50 --procs 4

API

Users can directly import the ORF search algorithm, written in cython, in their python ecosystem.

>>> import orfipy_core 
>>> seq='ATGCATGACTAGCATCAGCATCAGCAT'
>>> for start,stop,strand,description in orfipy_core.orfs(seq,minlen=3,maxlen=1000):
...     print(start,stop,strand,description)
... 
0 9 + ID=Seq_ORF.1;ORF_type=complete;ORF_len=9;ORF_frame=1;Start:ATG;Stop:TAG

orfipy_core.orfs function can take following arguments

• seq: Required input sequence (str)
• name ['Seq'] Name (str)
• minlen [0] min length (int)
• maxlen [1000000] max length (int)
• strand ['b'] Strand to use, (b)oth, (f)wd or (r)ev (char)
• starts [['TTG','CTG','ATG']] Start codons to use (list)
• stops=['TAA','TAG','TGA'] Stop codons to use (list)
• include_stop [False] Include stop codon in ORF (bool)
• partial3 [False] Report ORFs without a stop (bool)
• partial5 [False] Report ORFs without a start (bool)
• between_stops [False] Report ORFs defined as between stops (bool)

Comparison with getorf and OrfM

Comparison of orfipy features and performance with getorf and OrfM. For details see preprint and SI

Funding

This work is funded in part by the National Science Foundation award IOS 1546858, "Orphan Genes: An Untapped Genetic Reservoir of Novel Traits". This work used the Extreme Science and Engineering Discovery Environment (XSEDE), which is supported by National Science Foundation grant number ACI-1548562 (Bridges HPC environment through allocations TG-MCB190098 and TG-MCB200123 awarded from XSEDE and HPC Consortium).