czbiohub/rnaseq is a bioinformatics analysis pipeline used for RNA sequencing data.
The workflow processes raw data from FastQ inputs (FastQC, fastp), aligns the reads (STAR or HiSAT2), generates gene counts (htseq-count, StringTie) and performs extensive quality-control on the results (RSeQC, dupRadar, Preseq, edgeR, MultiQC). See the output documentation for more details of the results.
Additionally, the pipeline is expanded to be able to quantify transcript, exon, alternative splicing and TxRevise expressions. See optional quantification methods for details.
The pipeline is built using Nextflow, a bioinformatics workflow tool to run tasks across multiple compute infrastructures in a very portable manner. It comes with docker / singularity containers making installation trivial and results highly reproducible.
The czbiohub/rnaseq pipeline comes with documentation about the pipeline, found in the docs/ directory:
- Installation
- Pipeline configuration
- Running the pipeline (Gene expression)
- Running the pipeline (With additional quantification methods)
- Output and how to interpret the results
- Troubleshooting
The schema shown below represents the high level structure of the pipeline.
These scripts were originally written for use at the National Genomics Infrastructure, part of SciLifeLab in Stockholm, Sweden, by Phil Ewels (@ewels) and Rickard Hammarén (@Hammarn). They have since taken on a life of their own at Chan Zuckerberg Biohub where they are maintained by Olga Botvinnik.
Many thanks to other who have helped out along the way too, including (but not limited to): @Galithil, @pditommaso, @orzechoj, @apeltzer, @colindaven.