Skip to content

Latest commit

 

History

History
39 lines (30 loc) · 1.64 KB

Progenio.md

File metadata and controls

39 lines (30 loc) · 1.64 KB
Raw

�

Progenio

Progenio is an open-source proteogenomics pipeline for multistrain genome annotation, variant and Interupted Coding Sequence (ICDS) detection between strains and relative to one or more reference strains.

1. Clone the repository

git clone https://github.com/uct-cbio/proteomics-pipelines.git

2. Progenio with Singularity (Recommended)

Ensure singularity is installed, and make sure at least 2 cores and 4 GB of RAM are available.

cd proteomics-pipelines/singularity/progenio�
./create_image.sh

3. Run the example

Clone the example data:

git clone https://[email protected]/thys_potgieter/cbio-proteogenomics-tests.git

Create an interactive singularity session:

cd proteomics-pipelines/singularity/progenio
singularity shell progenio_v1.0.0.img

Run the example:

cd cbio-proteogenomics-tests/lib/testdata/proteogenomics_sample
sixframe.sh genomes/ genomes/config.yaml # create the six frame database used by the next step (genomes/sixframe/combined_prot.fasta)
Progenio.sh progenio_test.yml # analyze the data using the six frame translated database

Now adjust the config file for our own data.

4. Notes

MaxQuant protein search database entries are mapped to STOP to STOP genome six frame translated open reading frame sequences for each strain. Alligned ORFs with a evalue of 0 are considered mapped. ORFs are mapped to UniProt proteomes using BLAST, with the top hit per ORF for each proteome selected with an evalue cutoff of 0.0001 used for annotation. The mapped proteins from the Reference proteome are used for functional analysis using limma t-test and gene set enrichment analysis.