RNA-seq analysis pipeline for detection gene-fusions
nf-core/rnafusion is a bioinformatics best-practice analysis pipeline for RNA sequencing consisting of several tools designed for detecting and visualizing fusion genes. Results from up to 5 fusion callers tools are created, and are also aggregated, most notably in a pdf visualiation document, a vcf data collection file, and html and tsv reports.
On release, automated continuous integration tests run the pipeline on a full-sized dataset on the AWS cloud infrastructure. This ensures that the pipeline runs on AWS, has sensible resource allocation defaults set to run on real-world datasets, and permits the persistent storage of results to benchmark between pipeline releases and other analysis sources. The results obtained from the full-sized test can be viewed on the nf-core website.
In rnafusion the full-sized test includes reference building and fusion detection. The test dataset is taken from here.
--build_references triggers a parallel workflow to build references, which is a prerequisite to running the pipeline:
- Download ensembl fasta and gtf files
- Create STAR index
- Download Arriba references
- Download FusionCatcher references
- Download and build STAR-Fusion references
- Download Fusion-report DBs
- Input samplesheet check
- Concatenate fastq files per sample (cat)
- Reads quality control (FastQC)
- Optional trimming with fastp
- Arriba subworkflow
- STAR-fusion subworkflow
- Fusioncatcher subworkflow
- FusionCatcher fusion detection
- StringTie subworkflow
- Merge all fusions detected by the selected tools with Fusion-report
- Post-processing and analysis of data
- Present QC for raw reads (
- Compress bam files to cram with samtools view
As the reference building is computationally heavy (> 24h on HPC), it is recommended to test the pipeline with the
-stub parameter (creation of empty files):
First, build the references:
Then perform the analysis:
Please provide pipeline parameters via the CLI or Nextflow
-params-file option. Custom config files including those
provided by the
-c Nextflow option can be used to provide any configuration except for parameters;
- Conda is not currently supported; run with singularity or docker.
- Paths need to be absolute.
- GRCh38 is the only supported reference.
- Single-end reads are to be used as last-resort. Paired-end reads are recommended. FusionCatcher cannot be used with single-end reads shorter than 130 bp.
To see the results of an example test run with a full size dataset refer to the results tab on the nf-core website pipeline page. For more details about the output files and reports, please refer to the output documentation.
We thank the following people for their help in the development of this pipeline:
Contributions and Support
If you would like to contribute to this pipeline, please see the contributing guidelines.
If you use nf-core/rnafusion for your analysis, please cite it using the following doi: 10.5281/zenodo.3946477
An extensive list of references for the tools used by the pipeline can be found in the
You can cite the
nf-core publication as follows:
The nf-core framework for community-curated bioinformatics pipelines.
Philip Ewels, Alexander Peltzer, Sven Fillinger, Harshil Patel, Johannes Alneberg, Andreas Wilm, Maxime Ulysse Garcia, Paolo Di Tommaso & Sven Nahnsen.
Nat Biotechnol. 2020 Feb 13. doi: 10.1038/s41587-020-0439-x.