nextflow.config

/*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    mskilab-org/nf-jabba Nextflow config file
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    Default config options for all compute environments
----------------------------------------------------------------------------------------
*/

// Global default params, used in configs
params {


	// TODO nf-core: Specify your pipeline's command line flags
	// Input options (Mandatory!)
	input                      = null
	postprocess_bams           = true
    tumor_only                 = false

	// References
	genome                     = 'GATK.GRCh37'
	igenomes_base              = 's3://ngi-igenomes/igenomes'
	igenomes_ignore            = false
	mski_base                  = 's3://mskilab-pipeline'
	save_reference             = false
	build_only_index           = false                            // Only build the reference indexes
	download_cache             = false                            // Do not download annotation cache

	// Options to consider
	// Main options
	no_intervals               = false                            // Intervals will be built from the fasta file
	nucleotides_per_second     = 200000                           // Default interval size
	tools                      = null                             // No default Variant_Calling or Annotation tools
	skip_tools                 = null                             // All tools (markduplicates + baserecalibrator + QC) are used by default
	split_fastq                = 0                                // FASTQ files will not be split by default by FASTP, set to 50000000 for 50M reads/fastq

	// Modify FASTQ files (trim/split) with FASTP
	trim_fastq                 = false                            // No trimming by default
	clip_r1                    = 0
	clip_r2                    = 0
	three_prime_clip_r1        = 0
	three_prime_clip_r2        = 0
	trim_nextseq               = 0
	save_trimmed               = false
	save_split_fastqs          = false

	// Alignment
	aligner                    = 'fq2bam'                       // Default is gpu-accelerated fq2bam; bwa-mem, bwa-mem2 and dragmap can be used too
    fq2bam_mark_duplicates     = true                            // Whether fq2bam should mark duplicates, set false if not using fq2bam
    fq2bam_low_memory          = false                           // Set to true if using fq2bam with gpus that have <24GB memory
    optical_duplicate_pixel_distance = 2500                      // For computing optical duplicates, 2500 for NovaSeqX+
	save_mapped                = false                            // Mapped BAMs are saved
	save_output_as_bam         = true                            // Output files from alignment are saved as bam by default and not as cram files
	seq_center                 = null                            // No sequencing center to be written in read group CN field by aligner
	seq_platform               = null                            // Default platform written in read group PL field by aligner, null by default.

	// Structural Variant Calling
	error_rate                 = 0.01                            // Default error_rate for Svaba

	//indel_mask                 = null                            // Must provide blacklist bed file for indels based on genome to run Svaba

	// SV filtering (tumor only)
	pad_junc_filter            = 1000                            // Default Padding for SV Junction filtering

    // AMBER options
    target_bed_amber = null   // pass a bed file with target regions for running AMBER on targeted sequencing sample

    // HetPileups options
    filter_hets      = "TRUE"
    max_depth   = 1000

    // SigprofilerAssignment options
    sigprofilerassignment_cosmic_version = 3.4

	// fragCounter options
	midpoint_frag              = "TRUE"                           // If TRUE only count midpoint if FALSE then count bin footprint of every fragment interval: Default=TRUE
	windowsize_frag            = 1000                            // Window / bin size : Default=200 (but dryclean uses 1000 binsize)
	minmapq_frag               = 60                              // Minimal map quality : Default = 1
	paired_frag                = "TRUE"                           // Is the dataset paired : Default = TRUE
	exome_frag                 = "FALSE"	                        // Use exons as bins instead of fixed window : Default = FALSE

    // Dryclean options
    center_dryclean                    = "TRUE"
    cbs_dryclean                         = "FALSE"
    cnsignif_dryclean                    = 0.00001
    wholeGenome_dryclean                 = "TRUE"
    blacklist_dryclean                   = "FALSE"
    blacklist_path_dryclean              = "NA"
    germline_filter_dryclean             = "FALSE"
    germline_file_dryclean               = "NA"
    human_dryclean                       = "TRUE"
    field_dryclean                       = "reads"

	// ASCAT options
	field_ascat                         = "foreground"
	hets_thresh_ascat                   = 0.2
	penalty_ascat                       = 70
	gc_correct_ascat                    = "TRUE"
	rebin_width_ascat                   = 50000
	from_maf_ascat                      = "FALSE"

    // CBS options
    cnsignif_cbs                        = 0.01
    field_cbs                           = "foreground"
    name_cbs                            = "tumor"

	// JaBbA options
    blacklist_junctions_jabba = "NULL"
    geno_jabba = "FALSE"
    indel_jabba = "exclude"
    tfield_jabba = "tier"
    iter_jabba = 2
    rescue_window_jabba = 10000
    rescue_all_jabba = "TRUE"
    nudgebalanced_jabba = "TRUE"
    edgenudge_jabba = 0.1
    strict_jabba = "FALSE"
    allin_jabba = "FALSE"
    field_jabba = "foreground"
    maxna_jabba = 0.9
    ploidy_jabba = "NA"
    purity_jabba = "NA"
    pp_method_jabba = "ppgrid"
    cnsignif_jabba = 0.00001
    slack_jabba = 20
    linear_jabba = "FALSE"
    tilim_jabba = 7200
    epgap_jabba = 0.000001
    fix_thres_jabba = -1
    lp_jabba = "TRUE"
    ism_jabba = "TRUE"
    filter_loose_jabba = "FALSE"
    gurobi_jabba = "FALSE"
    verbose_jabba = "TRUE"

    // Allelic CN (Non-Integer Balance)
	field_non_integer_balance = "foreground"
	hets_thresh_non_integer_balance = 0.2
	overwrite_non_integer_balance = "TRUE"
	lambda_non_integer_balance = 20
	allin_non_integer_balance = "TRUE"
	fix_thresh_non_integer_balance = 10
	nodebounds_non_integer_balance = "TRUE"
	ism_non_integer_balance = "FALSE"
	epgap_non_integer_balance = 0.000001
	tilim_non_integer_balance = 6000
	gurobi_non_integer_balance = "FALSE"
	pad_non_integer_balance = 101

    // Allelic CN (LP-Phased Balance)
	lambda_lp_phased_balance = 100
	cnloh_lp_phased_balance = "TRUE"
	major_lp_phased_balance = "TRUE"
	allin_lp_phased_balance = "TRUE"
	marginal_lp_phased_balance = "TRUE"
	from_maf_lp_phased_balance = "FALSE"
	ism_lp_phased_balance = "FALSE"
	epgap_lp_phased_balance = 0.001
	hets_thresh_lp_phased_balance = 0.2
	min_bins_lp_phased_balance = 3
	min_width_lp_phased_balance = 0
	trelim_lp_phased_balance = 32000
	reward_lp_phased_balance = 10
	nodefileind_lp_phased_balance = 3
	tilim_lp_phased_balance = 6000

    // HRDetect
    hrdetect_mask = null

    // Variant Calling
    only_paired_variant_calling   = false // if true, skips germline variant calling for normal-paired samples
	ascat_ploidy                  = null  // default value for ASCAT
    ascat_min_base_qual           = 20    // default value for ASCAT
    ascat_min_counts              = 10    // default value for ASCAT
    ascat_min_map_qual            = 35    // default value for ASCAT
    ascat_purity                  = null  // default value for ASCAT
    cf_ploidy                     = "2"   // default value for Control-FREEC
    cf_coeff                      = 0.05  // default value for Control-FREEC
    cf_contamination              = 0     // default value for Control-FREEC
    cf_contamination_adjustment   = false // by default we are not using this in Control-FREEC
    cf_mincov                     = 0     // ControlFreec default values
    cf_minqual                    = 0     // ControlFreec default values
    cf_window                     = null  // by default we are not using this in Control-FREEC
    cnvkit_reference              = null  // by default the reference is build from the fasta file
    concatenate_vcfs              = false // by default we don't concatenate the germline-vcf-files
    ignore_soft_clipped_bases     = false // no --dont-use-soft-clipped-bases for GATK Mutect2
    wes                           = false // Set to true, if data is exome/targeted sequencing data. Used to use correct models in various variant callers
    joint_germline                = false // g.vcf & joint germline calling are not run by default if HaplotypeCaller is selected
    joint_mutect2                = false // if true, enables patient-wise multi-sample somatic variant calling
    sentieon_haplotyper_emit_mode = "variant" // default value for Sentieon haplotyper

    // Annotation
    dbnsfp                    = null // No dbnsfp processed file
    dbnsfp_consequence        = null // No default consequence for dbnsfp plugin
    dbnsfp_fields             = "rs_dbSNP,HGVSc_VEP,HGVSp_VEP,1000Gp3_EAS_AF,1000Gp3_AMR_AF,LRT_score,GERP++_RS,gnomAD_exomes_AF" // Default fields for dbnsfp plugin
    dbnsfp_tbi                = null // No dbnsfp processed file index
    outdir_cache              = null // No default outdir cache
    snpeff_cache              = 's3://annotation-cache/snpeff_cache/'
    spliceai_indel            = null // No spliceai_indel file
    spliceai_indel_tbi        = null // No spliceai_indel file index
    spliceai_snv              = null // No spliceai_snv file
    spliceai_snv_tbi          = null // No spliceai_snv file index
    use_annotation_cache_keys = false
    vep_cache                 = null //'s3://annotation-cache/vep_cache/' but not using it
    vep_custom_args           = "--everything --filter_common --per_gene --total_length --offline --format vcf" // Default arguments for VEP
    vep_dbnsfp                = null // dbnsfp plugin disabled within VEP
    vep_include_fasta         = false // Don't use fasta file for annotation with VEP
    vep_loftee                = null // loftee plugin disabled within VEP
    vep_out_format            = "vcf"
    vep_spliceai              = null // spliceai plugin disabled within VEP
    vep_spliceregion          = null // spliceregion plugin disabled within VEP


	// MultiQC options
	multiqc_config             = null
	multiqc_title              = null
	multiqc_logo               = null
	max_multiqc_email_size     = '25.MB'
	multiqc_methods_description = null


	// Boilerplate options
	outdir                     = null
	publish_dir_mode           = 'copy'
	email                      = null
	email_on_fail              = null
	plaintext_email            = false
	monochrome_logs            = false
	hook_url                   = null
	help                       = false
	version                    = false

	// Config options
	config_profile_name        = null
	config_profile_description = null
	custom_config_version      = 'master'
	custom_config_base         = "https://raw.githubusercontent.com/nf-core/configs/${params.custom_config_version}"
	config_profile_contact     = null
	config_profile_url         = null

	// Max resource options
	// Defaults only, expecting to be overwritten
	max_memory                 = '256.GB'
	max_cpus                   = 16
	max_time                   = '360.h'
    max_accelerator            = 8

	// Schema validation default options
	validationFailUnrecognisedParams = false
	validationLenientMode            = false
	validationSchemaIgnoreParams     = 'genomes, cf_ploidy'
	validationShowHiddenParams       = false
	validate_params                  = false

}

// Load base.config by default for all pipelines
includeConfig 'conf/base.config'

// Load nf-core custom profiles from different Institutions
try {
	includeConfig "${params.custom_config_base}/nfcore_custom.config"
} catch (Exception e) {
	System.err.println("WARNING: Could not load nf-core/config profiles: ${params.custom_config_base}/nfcore_custom.config")
}

// Load mskilab-org/nf-jabba custom profiles from different institutions.
// Warning: Uncomment only if a pipeline-specific instititutional config already exists on nf-core/configs!
// try {
//   includeConfig "${params.custom_config_base}/pipeline/nfjabba.config"
// } catch (Exception e) {
//   System.err.println("WARNING: Could not load nf-core/config/nfjabba profiles: ${params.custom_config_base}/pipeline/nfjabba.config")
// }

profiles {
    chr21_test {
        process {
            withName: 'GRIDSS_GRIDSS' {
                cpus            = { check_max( 4 * task.attempt, 'cpus' ) }
                memory          = { check_max( 32.GB * task.attempt, 'memory' ) }
                time            = { check_max( 2.h  * task.attempt, 'time'    ) }
            }
            withName: 'GRIDSS_SOMATIC' {
                cpus            = { check_max( 2 * task.attempt, 'cpus' ) }
                memory          = { check_max( 8.GB * task.attempt, 'memory' ) }
                time            = { check_max( 1.h  * task.attempt, 'time'    ) }
            }
            withName: 'SAGE_SOMATIC' {
                cpus            = { check_max( 20 * task.attempt, 'cpus' ) }
                memory          = { check_max( 16.GB * task.attempt, 'memory' ) }
                time            = { check_max( 1.h  * task.attempt, 'time'    ) }
            }
            withName: 'SAGE_GERMLINE' {
                cpus            = { check_max( 8 * task.attempt, 'cpus' ) }
                memory          = { check_max( 16.GB * task.attempt, 'memory' ) }
                time            = { check_max( 1.h  * task.attempt, 'time'    ) }
            }
            withName: 'JABBA' {
                cpus            = { check_max( 4 * task.attempt, 'cpus' ) }
                memory          = { check_max( 16.GB * task.attempt, 'memory' ) }
                time            = { check_max( 1.h  * task.attempt, 'time'    ) }
            }
            withName: 'FUSIONS' {
                cpus            = { check_max( 4 * task.attempt, 'cpus' ) }
                memory          = { check_max( 16.GB * task.attempt, 'memory' ) }
                time            = { check_max( 1.h  * task.attempt, 'time'    ) }
            }
        }
    }
	debug {
		dumpHashes             = true
		process.beforeScript   = 'echo $HOSTNAME'
        	cleanup                = false
    	}
	conda {
		conda.enabled          = true
		docker.enabled         = false
		singularity.enabled    = false
		podman.enabled         = false
		shifter.enabled        = false
		charliecloud.enabled   = false
		apptainer.enabled      = false
	}
	mamba {
		conda.enabled          = true
		conda.useMamba         = true
		docker.enabled         = false
		singularity.enabled    = false
		podman.enabled         = false
		shifter.enabled        = false
		charliecloud.enabled   = false
		apptainer.enabled      = false
	}
	docker {
		docker.enabled         = true
		docker.userEmulation   = true
		conda.enabled          = false
		singularity.enabled    = false
		podman.enabled         = false
		shifter.enabled        = false
		charliecloud.enabled   = false
		apptainer.enabled      = false
	}
	arm {
		docker.runOptions = '-u $(id -u):$(id -g) --platform=linux/amd64'
	}
	singularity {
		singularity.enabled    = true
		singularity.autoMounts = true
        singularity.runOptions = '--env CPLEX_DIR=/opt/cplex'
		conda.enabled          = false
		docker.enabled         = false
		podman.enabled         = false
		shifter.enabled        = false
		charliecloud.enabled   = false
		apptainer.enabled      = false
	}
	podman {
		podman.enabled         = true
		conda.enabled          = false
		docker.enabled         = false
		singularity.enabled    = false
		shifter.enabled        = false
		charliecloud.enabled   = false
		apptainer.enabled      = false
	}
	shifter {
		shifter.enabled        = true
		conda.enabled          = false
		docker.enabled         = false
		singularity.enabled    = false
		podman.enabled         = false
		charliecloud.enabled   = false
		apptainer.enabled      = false
	}
	charliecloud {
		charliecloud.enabled   = true
		conda.enabled          = false
		docker.enabled         = false
		singularity.enabled    = false
		podman.enabled         = false
		shifter.enabled        = false
		apptainer.enabled      = false
	}
	apptainer {
		apptainer.enabled      = true
		conda.enabled          = false
		docker.enabled         = false
		singularity.enabled    = false
		podman.enabled         = false
		shifter.enabled        = false
		charliecloud.enabled   = false
	}
	gitpod {
		executor.name          = 'local'
		executor.cpus          = 16
		executor.memory        = 60.GB
	}

	//basic test config files
	test      { includeConfig 'conf/test.config'      }
	test_full { includeConfig 'conf/test_full.config' }

}

// Set default registry for Apptainer, Docker, Podman and Singularity independent of -profile
// Will not be used unless Apptainer / Docker / Podman / Singularity are enabled
// Set to your registry if you have a mirror of containers
apptainer.registry   = 'quay.io'
docker.registry      = 'quay.io'
podman.registry      = 'quay.io'
singularity.registry = 'quay.io'

// Nextflow plugins
plugins {
	id 'nf-validation' // Validation of pipeline parameters and creation of an input channel from a sample sheet
}

// Load igenomes.config if required
if (!params.igenomes_ignore) {
	includeConfig 'conf/igenomes.config'
} else {
	params.genomes = [:]
}
// Export these variables to prevent local Python/R libraries from conflicting with those in the container
// The JULIA depot path has been adjusted to a fixed path `/usr/local/share/julia` that needs to be used for packages in the container.
// See https://apeltzer.github.io/post/03-julia-lang-nextflow/ for details on that. Once we have a common agreement on where to keep Julia packages, this is adjustable.

env {
	PYTHONNOUSERSITE = 1
	R_PROFILE_USER   = "/.Rprofile"
	R_ENVIRON_USER   = "/.Renviron"
	JULIA_DEPOT_PATH = "/usr/local/share/julia"
}

// Capture exit codes from upstream processes when piping
process.shell = ['/bin/bash', '-euo', 'pipefail']

def trace_timestamp = new java.util.Date().format( 'yyyy-MM-dd_HH-mm-ss')
timeline {
	enabled = true
	file    = "${params.outdir}/pipeline_info/execution_timeline_${trace_timestamp}.html"
}
report {
	enabled = true
	file    = "${params.outdir}/pipeline_info/execution_report_${trace_timestamp}.html"
}
trace {
	enabled = true
	file    = "${params.outdir}/pipeline_info/execution_trace_${trace_timestamp}.txt"
}
dag {
	enabled = true
	file    = "${params.outdir}/pipeline_info/pipeline_dag_${trace_timestamp}.html"
}

manifest {
	name            = 'mskilab-org/nf-jabba'
	author          = """Tanubrata Dey and Shihab Dider"""
	homePage        = 'https://github.com/mskilab-org/nf-casereports'
	description     = """Clinical Core Pipeline for MSkiLab"""
	mainScript      = 'main.nf'
	nextflowVersion = '!>=23.04.0'
	version         = '1.0dev'
	doi             = ''
}

// Load modules.config for DSL2 module specific options
includeConfig 'conf/modules.config'


// prepare reference, pre-alignment steps
includeConfig 'conf/modules/prepare_cache.config'
includeConfig 'conf/modules/prepare_genome.config'
includeConfig 'conf/modules/prepare_intervals.config'

// Alignment configurations
includeConfig 'conf/modules/aligner.config'
includeConfig 'conf/modules/alignment_to_fastq.config'
includeConfig 'conf/modules/markduplicates.config'
includeConfig 'conf/modules/prepare_recalibration.config'
includeConfig 'conf/modules/recalibrate.config'
includeConfig 'conf/modules/trimming.config'

// AMBER configurations
includeConfig 'conf/modules/amber.config'

// SV configurations
includeConfig 'conf/modules/structural_variants.config'

// fragCounter configurations
includeConfig 'conf/modules/fragcounter.config'

// SNV calling configurations (sage and strelka)
includeConfig 'conf/modules/snv_calling.config'

// Variant annotation configurations (snpeff and pave)
includeConfig 'conf/modules/variant_annotation.config'

// Signatures configurations (snpeff and pave)
includeConfig 'conf/modules/sigprofilerassignment.config'

// Dryclean configurations
includeConfig 'conf/modules/dryclean.config'

// PURPLE configurations
includeConfig 'conf/modules/purple.config'

// CBS configurations
includeConfig 'conf/modules/cbs.config'

// JaBbA configurations
includeConfig 'conf/modules/jabba.config'

// Events configurations
includeConfig 'conf/modules/events.config'

// Fusions configurations
includeConfig 'conf/modules/fusions.config'

// Allelic CN configurations
includeConfig 'conf/modules/allelic_cn.config'

// HRDetect configurations
includeConfig 'conf/modules/hrdetect.config'

// Function to ensure that resource requirements don't go beyond
// a maximum limit
def check_max(obj, type) {
    if (type == 'memory') {
        try {
            if (obj.compareTo(params.max_memory as nextflow.util.MemoryUnit) == 1)
                return params.max_memory as nextflow.util.MemoryUnit
            else
                return obj
        } catch (all) {
            println "   ### ERROR ###   Max memory '${params.max_memory}' is not valid! Using default value: $obj"
            return obj
        }
    } else if (type == 'time') {
        try {
            if (obj.compareTo(params.max_time as nextflow.util.Duration) == 1)
                return params.max_time as nextflow.util.Duration
            else
                return obj
        } catch (all) {
            println "   ### ERROR ###   Max time '${params.max_time}' is not valid! Using default value: $obj"
            return obj
        }
    } else if (type == 'cpus') {
        try {
            return Math.min( obj, params.max_cpus as int )
        } catch (all) {
            println "   ### ERROR ###   Max cpus '${params.max_cpus}' is not valid! Using default value: $obj"
            return obj
        }
    } else if (type == 'accelerator') {
        try {
            return Math.min( obj, params.max_accelerator as int )
        } catch (all) {
            println "   ### ERROR ###   Max gpu '${params.max_accelerator}' is not valid! Using default value: $obj"
            return obj
        }
    }
}