Skip to main content

A RNAseq pipeline from raw reads to feature counts

Project description

This is is the rnaseq pipeline from the Sequana projet

Overview:RNASeq analysis from raw data to feature counts
Input:A set of Fastq Files and genome reference and annotation.
Output:MultiQC reports and feature Counts
Status:Production
Citation:Cokelaer et al, (2017), ‘Sequana’: a Set of Snakemake NGS pipelines, Journal of Open Source Software, 2(16), 352, JOSS DOI doi:10.21105/joss.00352

Installation

You must install Sequana first:

pip install sequana

Then, just install this package:

pip install sequana_rnaseq

Usage

sequana_pipelines_rnaseq --help
sequana_pipelines_rnaseq --input-directory DATAPATH --genome-directory --aligner STAR

This creates a directory with the pipeline and configuration file. You will then need to execute the pipeline:

cd rnaseq
sh rnaseq.sh  # for a local run

This launch a snakemake pipeline. If you are familiar with snakemake, you can retrieve the pipeline itself and its configuration files and then execute the pipeline yourself with specific parameters:

snakemake -s rnaseq.rules -c config.yaml --cores 4 --stats stats.txt

Or use sequanix interface.

Requirements

This pipelines requires the following executable(s):

  • bowtie
  • bowtie2
  • STAR
  • featureCounts
  • picard
  • multiqc

More may be needed depending on the configuration file options. For instance, you may use fastq_screen, in which case you need to install it and configure it.

https://raw.githubusercontent.com/sequana/sequana_rnaseq/master/sequana_pipelines/rnaseq/dag.png

Details

This pipeline runs rnaseq in parallel on the input fastq files (paired or not). A brief sequana summary report is also produced.

Rules and configuration details

Here is the latest documented configuration file to be used with the pipeline. Each rule used in the pipeline may have a section in the configuration file.

Warning

the RNAseQC rule is switch off and is not currently functional in version 0.9.X

Changelog

0.9: Major refactorisation. remove sartools, kraken rules. Indexing is now optional and can be set in the configuration. Configuration file is simplified with a general section to enter the genome location and aligner. Fixed rules in sequana (0.8.0) that were not up-to-date with several executables used in the pipeline including picard, fastq_screen, etc. See Sequana Changelog for details with respect to rules changes. Copying the feature counts in main directory ready to use for a differential analysis.

Project details


Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Files for sequana-rnaseq, version 0.9.3
Filename, size File type Python version Upload date Hashes
Filename, size sequana_rnaseq-0.9.3.tar.gz (85.8 kB) File type Source Python version None Upload date Hashes View hashes

Supported by

Elastic Elastic Search Pingdom Pingdom Monitoring Google Google BigQuery Sentry Sentry Error logging AWS AWS Cloud computing DataDog DataDog Monitoring Fastly Fastly CDN SignalFx SignalFx Supporter DigiCert DigiCert EV certificate StatusPage StatusPage Status page