STAR-ARRIBA
A subworkflow for fusion detection from RNA-seq data with arriba. The fusion calling is based on splice-aware, chimeric alignments done with STAR. STAR is used with specific parameters to ensure optimal functionality of the arriba fusion detection, for details, see the documentation.
Example
This meta-wrapper can be used by integrating the following into your workflow:
rule star_index:
input:
fasta="resources/genome.fasta",
gtf="resources/genome.gtf",
output:
directory("resources/star_genome"),
threads: 4
params:
sjdbOverhang=100,
extra="--genomeSAindexNbases 2",
log:
"logs/star_index_genome.log",
cache: True # mark as eligible for between workflow caching
wrapper:
"v3.8.0-49-g6f33607/bio/star/index"
rule star_align:
input:
# use a list for multiple fastq files for one sample
# usually technical replicates across lanes/flowcells
fq1="reads/{sample}_R1.1.fastq",
fq2="reads/{sample}_R2.1.fastq", #optional
idx="resources/star_genome",
annotation="resources/genome.gtf",
output:
# see STAR manual for additional output files
aln="star/{sample}/Aligned.out.bam",
reads_per_gene="star/{sample}/ReadsPerGene.out.tab",
log:
"logs/star/{sample}.log",
params:
# specific parameters to work well with arriba
extra=lambda wc, input: f"--quantMode GeneCounts --sjdbGTFfile {input.annotation}"
" --outSAMtype BAM Unsorted --chimSegmentMin 10 --chimOutType WithinBAM SoftClip"
" --chimJunctionOverhangMin 10 --chimScoreMin 1 --chimScoreDropMax 30 --chimScoreJunctionNonGTAG 0"
" --chimScoreSeparation 1 --alignSJstitchMismatchNmax 5 -1 5 5 --chimSegmentReadGapMax 3",
threads: 12
wrapper:
"v3.8.0-49-g6f33607/bio/star/align"
rule arriba:
input:
bam=rules.star_align.output.aln,
genome="resources/genome.fasta",
annotation="resources/genome.gtf",
# optional: # A custom tsv containing identified artifacts, such as read-through fusions of neighbouring genes.
# default blacklists are selected via blacklist parameter
# see https://arriba.readthedocs.io/en/latest/input-files/#blacklist
custom_blacklist=[],
output:
fusions="results/arriba/{sample}.fusions.tsv",
discarded="results/arriba/{sample}.fusions.discarded.tsv",
params:
# required if blacklist or known_fusions is set
genome_build="GRCh38",
default_blacklist=False,
default_known_fusions=True,
extra="",
log:
"logs/arriba/{sample}.log",
threads: 1
wrapper:
"v3.8.0-49-g6f33607/bio/arriba"
Note that input, output and log file paths can be chosen freely, as long as the dependencies between the rules remain as listed here. For additional parameters in each individual wrapper, please refer to their corresponding documentation (see links below).
When running with
snakemake --use-conda
the software dependencies will be automatically deployed into an isolated environment before execution.
Used wrappers
The following individual wrappers are used in this meta-wrapper:
Please refer to each wrapper in above list for additional configuration parameters and information about the executed code.