Newer
Older
rule remove_adapters_cutadapt:
BIOPZ-Katsantoni Maria
committed
'''
Remove adapters
'''
reads = os.path.join(
config["output_dir"],
"samples",
"{sample}",
"start",
"{sample}.fq1.fastq.gz")
BIOPZ-Katsantoni Maria
committed
reads = temp(os.path.join(
BIOPZ-Katsantoni Maria
committed
config["output_dir"],
BIOPZ-Katsantoni Maria
committed
"{sample}",
"{sample}.se.remove_adapters_mate1.fastq.gz"))
BIOPZ-Katsantoni Maria
committed
BIOPZ-Katsantoni Maria
committed
adapters_3 = lambda wildcards:
BIOPZ-Katsantoni Maria
committed
get_sample(
'fq1_3p',
search_id='index',
search_value=wildcards.sample),
BIOPZ-Katsantoni Maria
committed
adapters_5 = lambda wildcards:
BIOPZ-Katsantoni Maria
committed
get_sample(
'fq1_5p',
search_id='index',
search_value=wildcards.sample)
singularity:
BIOPZ-Katsantoni Maria
committed
BIOPZ-Katsantoni Maria
committed
BIOPZ-Katsantoni Maria
committed
stderr = os.path.join(
config["log_dir"],
BIOPZ-Katsantoni Maria
committed
"{sample}",
BIOPZ-Katsantoni Maria
committed
stdout = os.path.join(
config["log_dir"],
BIOPZ-Katsantoni Maria
committed
"{sample}",
shell:
"(cutadapt \
-j {threads} \
-m 10 \
-a {params.adapters_3} \
-g {params.adapters_5} \
-o {output.reads} \
BIOPZ-Katsantoni Maria
committed
1> {log.stdout} 2> {log.stderr}"
rule remove_polya_cutadapt:
BIOPZ-Katsantoni Maria
committed
'''
Remove ployA tails
'''
BIOPZ-Katsantoni Maria
committed
reads = os.path.join(
config["output_dir"],
BIOPZ-Katsantoni Maria
committed
"{sample}",
BIOPZ-Katsantoni Maria
committed
reads = temp(os.path.join(
BIOPZ-Katsantoni Maria
committed
config["output_dir"],
BIOPZ-Katsantoni Maria
committed
"{sample}",
"{sample}.se.remove_polya_mate1.fastq.gz"))
BIOPZ-Katsantoni Maria
committed
BIOPZ-Katsantoni Maria
committed
polya_3 = lambda wildcards:
BIOPZ-Katsantoni Maria
committed
get_sample(
'fq1_polya_3p',
search_id='index',
search_value=wildcards.sample),
BIOPZ-Katsantoni Maria
committed
get_sample(
'fq1_polya_5p',
search_id='index',
search_value=wildcards.sample)
BIOPZ-Katsantoni Maria
committed
singularity:
BIOPZ-Katsantoni Maria
committed
BIOPZ-Katsantoni Maria
committed
BIOPZ-Katsantoni Maria
committed
stderr = os.path.join(
config["log_dir"],
BIOPZ-Katsantoni Maria
committed
"{sample}",
BIOPZ-Katsantoni Maria
committed
stdout = os.path.join(
config["log_dir"],
BIOPZ-Katsantoni Maria
committed
"{sample}",
BIOPZ-Katsantoni Maria
committed
shell:
"(cutadapt \
-j {threads} \
-O 1 \
-m 10 \
-a {params.polya_3} \
-o {output.reads} \
BIOPZ-Katsantoni Maria
committed
1> {log.stdout} 2> {log.stderr}"
rule map_genome_star:
BIOPZ-Katsantoni Maria
committed
'''
Map to genome using STAR
'''
input:
index = lambda wildcards:
os.path.join(
config["star_indexes"],
BIOPZ-Katsantoni Maria
committed
get_sample('organism', search_id='index', search_value=wildcards.sample),
get_sample('index_size', search_id='index', search_value=wildcards.sample),
BIOPZ-Katsantoni Maria
committed
"STAR_index",
"chrNameLength.txt"),
reads = os.path.join(
config["output_dir"],
BIOPZ-Katsantoni Maria
committed
"{sample}",
BIOPZ-Katsantoni Maria
committed
BIOPZ-Katsantoni Maria
committed
bam = os.path.join(
config["output_dir"],
BIOPZ-Katsantoni Maria
committed
"{sample}",
"map_genome",
"{sample}.se.Aligned.sortedByCoord.out.bam"),
BIOPZ-Katsantoni Maria
committed
logfile = os.path.join(
config["output_dir"],
BIOPZ-Katsantoni Maria
committed
"{sample}",
"map_genome",
BIOPZ-Katsantoni Maria
committed
params:
sample_id = "{sample}",
index = lambda wildcards:
os.path.abspath(os.path.join(
config["star_indexes"],
BIOPZ-Katsantoni Maria
committed
get_sample('organism', search_id='index', search_value=wildcards.sample),
get_sample('index_size', search_id='index', search_value=wildcards.sample),
outFileNamePrefix = os.path.join(
BIOPZ-Katsantoni Maria
committed
config["output_dir"],
BIOPZ-Katsantoni Maria
committed
"{sample}",
"map_genome",
multimappers = lambda wildcards:
BIOPZ-Katsantoni Maria
committed
get_sample(
'multimappers',
search_id='index',
search_value=wildcards.sample),
soft_clip = lambda wildcards:
BIOPZ-Katsantoni Maria
committed
get_sample(
'soft_clip',
search_id='index',
search_value=wildcards.sample),
pass_mode = lambda wildcards:
BIOPZ-Katsantoni Maria
committed
get_sample(
'pass_mode',
search_id='index',
search_value=wildcards.sample)
BIOPZ-Katsantoni Maria
committed
singularity:
BIOPZ-Katsantoni Maria
committed
threads: 12
BIOPZ-Katsantoni Maria
committed
BIOPZ-Katsantoni Maria
committed
stderr = os.path.join(
config["log_dir"],
BIOPZ-Katsantoni Maria
committed
"{sample}",
BIOPZ-Katsantoni Maria
committed
shell:
"(STAR \
-- twopassMode {params.pass_mode} \
--runThreadN {threads} \
--genomeDir {params.index} \
--readFilesIn {input.reads} \
--readFilesCommand zcat \
--outFilterMultimapNmax {params.multimappers} \
--outFileNamePrefix {params.outFileNamePrefix} \
--outSAMattributes All \
--outStd BAM_SortedByCoordinate \
--outSAMtype BAM SortedByCoordinate \
--outFilterType BySJout \
--outSAMattrRGline ID:rnaseq_pipeline SM:{params.sample_id} \
BIOPZ-Katsantoni Maria
committed
--alignEndsType {params.soft_clip} > {output.bam};) \
2> {log.stderr}"
rule quantification_salmon:
BIOPZ-Katsantoni Maria
committed
'''
Quantification at transcript and gene level using Salmon
'''
input:
reads = os.path.join(
BIOPZ-Katsantoni Maria
committed
config["output_dir"],
BIOPZ-Katsantoni Maria
committed
"{sample}",
index = lambda wildcards:
os.path.join(
config["salmon_indexes"],
BIOPZ-Katsantoni Maria
committed
get_sample(
'organism',
search_id='index',
search_value=wildcards.sample),
get_sample(
'kmer',
search_id='index',
search_value=wildcards.sample),
"salmon.idx"),
BIOPZ-Katsantoni Maria
committed
gtf = lambda wildcards:
os.path.abspath(get_sample(
BIOPZ-Katsantoni Maria
committed
'gtf',
search_id='index',
search_value=wildcards.sample))
BIOPZ-Katsantoni Maria
committed
output:
gn_estimates = os.path.join(
BIOPZ-Katsantoni Maria
committed
config["output_dir"],
BIOPZ-Katsantoni Maria
committed
"{sample}",
"quant.genes.sf"),
tr_estimates = os.path.join(
BIOPZ-Katsantoni Maria
committed
config["output_dir"],
BIOPZ-Katsantoni Maria
committed
"{sample}",
"quant.sf"),
meta_info = os.path.join(
config["output_dir"],
"samples",
"{sample}",
"{sample}.salmon.se",
"aux_info",
"meta_info.json"),
flenDist = os.path.join(
config["output_dir"],
"samples",
"{sample}",
"{sample}.salmon.se",
"libParams",
"flenDist.txt")
shadow: "minimal"
BIOPZ-Katsantoni Maria
committed
params:
output_dir = os.path.join(
BIOPZ-Katsantoni Maria
committed
config["output_dir"],
BIOPZ-Katsantoni Maria
committed
"{sample}",
libType = lambda wildcards:
BIOPZ-Katsantoni Maria
committed
get_sample(
'libtype',
search_id='index',
search_value=wildcards.sample),
BIOPZ-Katsantoni Maria
committed
fraglen = lambda wildcards:
BIOPZ-Katsantoni Maria
committed
get_sample(
'mean',
search_id='index',
search_value=wildcards.sample),
BIOPZ-Katsantoni Maria
committed
fragsd = lambda wildcards:
BIOPZ-Katsantoni Maria
committed
get_sample(
'sd',
search_id='index',
search_value=wildcards.sample)
BIOPZ-Katsantoni Maria
committed
stderr = os.path.join(
config["log_dir"],
BIOPZ-Katsantoni Maria
committed
"{sample}",
BIOPZ-Katsantoni Maria
committed
stdout = os.path.join(
config["log_dir"],
BIOPZ-Katsantoni Maria
committed
"{sample}",
BIOPZ-Katsantoni Maria
committed
threads: 12
singularity:
"docker://zavolab/salmon:1.1.0-slim"
BIOPZ-Katsantoni Maria
committed
shell:
"(salmon quant \
--libType {params.libType} \
--seqBias \
--validateMappings \
--threads {threads} \
BIOPZ-Katsantoni Maria
committed
--fldMean {params.fraglen} \
--fldSD {params.fragsd} \
--index {input.index} \
--geneMap {input.gtf} \
--unmatedReads {input.reads} \
BIOPZ-Katsantoni Maria
committed
-o {params.output_dir};) \
1> {log.stdout} 2> {log.stderr}"
BIOPZ-Gypas Foivos
committed
rule genome_quantification_kallisto:
BIOPZ-Katsantoni Maria
committed
'''
Quantification at transcript and gene level using Kallisto
'''
input:
reads = os.path.join(
BIOPZ-Katsantoni Maria
committed
config["output_dir"],
BIOPZ-Katsantoni Maria
committed
"{sample}",
index = lambda wildcards:
os.path.join(
config["kallisto_indexes"],
BIOPZ-Katsantoni Maria
committed
get_sample(
'organism',
search_id='index',
search_value=wildcards.sample),
"kallisto.idx")
BIOPZ-Katsantoni Maria
committed
output:
pseudoalignment = os.path.join(
BIOPZ-Katsantoni Maria
committed
config["output_dir"],
BIOPZ-Katsantoni Maria
committed
"{sample}",
"quant_kallisto",
"{sample}.se.kallisto.pseudo.sam"),
abundances = os.path.join(
config["output_dir"],
"samples",
"{sample}",
"quant_kallisto",
"abundance.h5")
shadow: "minimal"
BIOPZ-Katsantoni Maria
committed
params:
output_dir = os.path.join(
BIOPZ-Katsantoni Maria
committed
config["output_dir"],
BIOPZ-Katsantoni Maria
committed
"{sample}",
"quant_kallisto"),
fraglen = lambda wildcards:
BIOPZ-Katsantoni Maria
committed
get_sample(
'mean',
search_id='index',
search_value=wildcards.sample),
BIOPZ-Katsantoni Maria
committed
fragsd = lambda wildcards:
BIOPZ-Katsantoni Maria
committed
get_sample(
'sd',
search_id='index',
search_value=wildcards.sample),
BIOPZ-Katsantoni Maria
committed
directionality = lambda wildcards:
BIOPZ-Katsantoni Maria
committed
get_sample(
'kallisto_directionality',
search_id='index',
search_value=wildcards.sample)
BIOPZ-Katsantoni Maria
committed
threads: 8
BIOPZ-Katsantoni Maria
committed
stderr = os.path.join(
config["log_dir"],
BIOPZ-Katsantoni Maria
committed
"{sample}",
"genome_quantification_kallisto.se.stderr.log")
BIOPZ-Katsantoni Maria
committed
singularity:
BIOPZ-Katsantoni Maria
committed
shell:
"(kallisto quant \
-i {input.index} \
-o {params.output_dir} \
--single \
-l {params.fraglen} \
-s {params.fragsd} \
--pseudobam \
BIOPZ-Katsantoni Maria
committed
{input.reads} > {output.pseudoalignment};) \