RiboFlow
RiboFlow is a Nextflow based pipeline for processing ribosome profiling data.
Installation
Requirements
First, follow the instructions in Nextflow website and install Nextflow.
The easiest way of using RiboFLow is using Docker. If using Docker is not an option, you can install the dependencies using Conda and run RiboFlow without Docker.
Docker Option
Install Docker. Here is a tutorial for Ubuntu.
All remaining dependencies come in the Docker image ceniklab/riboflow. This image is automatically pulled by RiboFlow when run with Docker (see test runs below).
Conda Option
This option has been tested on Linux systems only.
Install Conda.
All other dependencies can be installed using the environment file, environment.yaml, in this repository.
git clone https://github.com/ribosomeprofiling/riboflow.git
conda env create -f riboflow/environment.yaml
The above command will create a conda environment called ribo and install dependencies in it. To start using RiboFlow, you need to activate the ribo environment.
conda activate ribo
Test Run
For fresh installations, before running RiboFlow on actual data, it is recommended to do a test run.
Clone this repository in a new folder and change your working directory to the RiboFlow folder.
mkdir rf_test_run && cd rf_test_run
git clone https://github.com/ribosomeprofiling/riboflow.git
cd riboflow
Obtain a copy of the sample data in the working directory.
git clone https://github.com/ribosomeprofiling/rf_sample_data.git
Run Using Docker
Provide the argument -profile docker_local
to Nextflow to indicate Docker use.
nextflow RiboFlow.groovy -params-file project.yaml -profile docker_local
Run Using Conda Environment
Make sure that you have created the conda environment, called ribo, using the instructions above. Then activate the conda environment.
conda activate ribo
If the above command fails to activate the ribo environment, try
source activate ribo
Now RiboFlow is ready to run.
nextflow RiboFlow.groovy -params-file project.yaml
Output
Pipeline run may take several minutes.
When finished, the resulting files are in the ./output
folder.
Mapping statistics are compiled in a csv file called stats.csv
ls output/stats/stats.csv
Ribosome occupancy data is in a single
ribo file called all.ribo
.
ls output/ribo/all.ribo
You can use RiboR or RiboPy to work with ribo files.
Actual Run
For running RiboFlow on actual data, files must be organized and a parameters file must be prepared. You can examine the sample run above to see an example.
- Organize your data. The following files are required for RiboFlow
- Ribosome profiling sequencing data: in gzipped fastq files
- Transcriptome Reference: Bowtie2 index files
- Filter Reference: Bowtie2 index files (typically for rRNA sequences)
- Annotation: A bed file defining CDS, UTR5 and UTR3 regions.
- Transcript Lengths: A two column tsv file containing transcript lengths
-
Prepare a custom
project.yaml
file. You can use the sample fileproject.yaml
, provided in this repository, as template. -
In
project.yaml
, provide RiboFlow parameters such asclip_arguments
, alignment arguments etc. You can simply modify the arguments in the sample fileproject.yaml
in this repository. -
You can adjust the hardware and computing environment settings in Nextflow configuration file(s). For Docker option, see
configs/docker_local.config
. If you are not using Docker, seeconfigs/local.config
. -
RNA-Seq data is optional for RiboFlow. So, if you do NOT have RNA-Seq data, in the project file, set
do_rnaseq: false
If you have RNA-Seq data to be paired with ribosome profiling data, see the Advanced Features below.
- Metadata is optional for RiboFlow.. If you do NOT have metadata, in the project file, set
do_metadata: false
If you have metadata, see Advanced Features below.
- Run RiboFlow using the new parameters file
project.yaml
.
Using Docker:
nextflow RiboFlow.groovy -params-file project.yaml -profile docker_local
Without Docker:
nextflow RiboFlow.groovy -params-file project.yaml
Advanced Features
RNA-Seq Data
If you have RNA-Seq data that you want to pair with ribosome profiling experiments,
provide the paths of the RNA-Seq (gzipped) fastq files in the configuration file in
input -> metadata. See the file project.yaml
in this repository for an example.
Note that the names in defining RNA-Seq files must match the names in definig ribosome profiling data.
Also turn set the do_rnaseq flag to true, in the project file:
do_rnaseq: true
Transcript abundance data will be stored in the output ribo file.
Metadata
If you have metadata files for the ribosome profiling experiments,
provide the paths of the metadata files (in yaml format) in the configuration file in
input -> metadata. See the file project.yaml
in this repository for an example.
Note that the names in defining metadata files must match the names in definig ribosome profiling data.
Also turn set the metadata flag to true, in the project file:
do_metadata: true
Metadata will be stored in the output ribo file.
nextflow pipeline
This repository is a template and a library repository to help you build nextflow pipeline. You can fork this repository to build your own pipeline. To get the last commits from this repository into your fork use the following commands:
git remote add upstream gitlab_lbmc:pipelines/nextflow.git
git pull upstream master
If you created your .config
file before version 0.4.0
you need to run the script src/.update_config.sh
to use the latest docker, singularity and conda configuration (don't forget to check your config files afterward for typos).
Getting Started
These instructions will get you a copy of the project up and running on your local machine for development and testing purposes. See deployment for notes on how to deploy the project on a live system.
Available tools
Projects using nextflow
A list of projects using nextflow at the LBMC.
Contributing
Please read CONTRIBUTING.md for details on our code of conduct, and the process for submitting pull requests to us.
Versioning
We use SemVer for versioning. For the versions available, see the tags on this repository.
Authors
- Laurent Modolo - Initial work
See also the list of contributors who participated in this project.
License
This project is licensed under the CeCiLL License- see the LICENSE file for details