Workflow: Genomic regions intersection and visualization

Fetched 2023-01-04 16:11:06 GMT

Genomic regions intersection and visualization ============================================== 1. Merges intervals within each of the filtered peaks files from ChIP/ATAC experiments 2. Overlaps merged intervals and assigns the nearest genes to them

children parents
Workflow as SVG
  • Selected
  • Default Values
  • Nested Workflows
  • Tools
  • Inputs/Outputs

Inputs

ID Type Title Doc
alias String Experiment short name/Alias
promoter_dist Integer (Optional) Max distance from gene TSS (in both direction) to assign interval to promoter

Max distance from gene TSS (in both direction) overlapping which the interval will be assigned to the promoter region

upstream_dist Integer (Optional) Max distance from promoter (only in upstream direction) to assign interval to upstream

Max distance from the promoter (only in upstream direction) overlapping which the interval will be assigned to the upstream region

annotation_file File [TSV] Genome to use for the nearest gene assignment

Tab-separated annotation file

intervals_files File[] [BED] Filtered ChIP/ATAC experiment

Filtered peaks file from ChIP/ATAC experiment formatted as headerless BED [chrom start end name] file

intervals_aliases String[] Filtered ChIP/ATAC experiment

Filtered ChIP/ATAC experiment alias

Steps

ID Runs Label Doc
sort_intervals
merge_intervals
overlap_intervals
../tools/intervene.cwl (CommandLineTool)

Intervene: a tool for intersection and visualization of multiple genomic regions ================================================================================

Hardcoded parameters: --type genomic use genomic regions as inputs --save-overlaps save overlapping regions/names for all the combinations as bed/txt --project intervene for the guatanteed access to results location -o results for the guatanteed access to results location

Skipped parameters: --filenames as it seems not logical to have it when default for --names already uses file names

When run in \"pairwise\" mode: --save-overlaps and --overlap-thresh set to null --htype and --diagonal added

When run in \"venn\" mode: --colors aren't ignored

remove_empty_columns
../tools/custom-bash.cwl (CommandLineTool)

Tool to run custom script set as `script` input with arguments from `param`. Default script runs sed command over the input file and exports results to the file with the same name as input's basename

get_overlapped_between_all
../tools/get-file-by-name.cwl (ExpressionTool)

Returns file the first file from input File[], that match input regex expression

annotate_overlapped_intervals
refactore_overlapped_intervals
annotate_overlapped_between_all
../tools/iaintersect.cwl (CommandLineTool)

Tool assigns each peak obtained from MACS2 to a gene and region (upstream, promoter, exon, intron, intergenic)

`default_output_filename` function returns output filename with sufix set as `ext` argument. Function is called when either `output_filename` or `log_filename` inputs are not provided.

refactore_overlapped_between_all
../tools/custom-bash.cwl (CommandLineTool)

Tool to run custom script set as `script` input with arguments from `param`. Default script runs sed command over the input file and exports results to the file with the same name as input's basename

Outputs

ID Type Label Doc
overlapped_plot File [PNG] Intervals overlap diagram

Intervals overlap diagram

intervene_stderr_log File [Textual format] intervene stderr log

intervene stderr log

intervene_stdout_log File [Textual format] intervene stdout log

intervene stdout log

overlapped_combinations File (Optional) [Textual format] Overlapped combinations file

Overlapped combinations file

overlapped_between_all_bed File [BED] BED file with common for all inputs intervals

BED file with common for all inputs intervals

overlapped_intervals_files File[] [TSV] Overlapped intervals files with the nearest genes assigned

Overlapped intervals files with the nearest genes assigned

annotated_overlapped_between_all_file File [TSV] Annotated intervals common for all intersected files

Annotated intervals common for all intersected files

Permalink: https://w3id.org/cwl/view/git/60854b5d299df91e135e05d02f4be61f6a310fbc/workflows/intervene.cwl