Workflow: MAnorm SE - quantitative comparison of ChIP-Seq single-read data

Fetched 2023-01-03 19:45:10 GMT

What is MAnorm? -------------- MAnorm is a robust model for quantitative comparison of ChIP-Seq data sets of TFs (transcription factors) or epigenetic modifications and you can use it for: * Normalization of two ChIP-seq samples * Quantitative comparison (differential analysis) of two ChIP-seq samples * Evaluating the overlap enrichment of the protein binding sites(peaks) * Elucidating underlying mechanisms of cell-type specific gene regulation How MAnorm works? ---------------- MAnorm uses common peaks of two samples as a reference to build the rescaling model for normalization, which is based on the empirical assumption that if a chromatin-associated protein has a substantial number of peaks shared in two conditions, the binding at these common regions will tend to be determined by similar mechanisms, and thus should exhibit similar global binding intensities across samples. The observed differences on common peaks are presumed to reflect the scaling relationship of ChIP-Seq signals between two samples, which can be applied to all peaks. What do the inputs mean? ---------------- ### General **Experiment short name/Alias** * short name for you experiment to identify among the others **ChIP-Seq SE sample 1** * previously analyzed ChIP-Seq single-read experiment to be used as Sample 1 **ChIP-Seq SE sample 2** * previously analyzed ChIP-Seq single-read experiment to be used as Sample 2 **Genome** * Reference genome to be used for gene assigning ### Advanced **Reads shift size for sample 1** * This value is used to shift reads towards 3' direction to determine the precise binding site. Set as half of the fragment length. Default 100 **Reads shift size for sample 2** * This value is used to shift reads towards 5' direction to determine the precise binding site. Set as half of the fragment length. Default 100 **M-value (log2-ratio) cutoff** * Absolute M-value (log2-ratio) cutoff to define biased (differential binding) peaks. Default: 1.0 **P-value cutoff** * P-value cutoff to define biased peaks. Default: 0.01 **Window size** * Window size to count reads and calculate read densities. 2000 is recommended for sharp histone marks like H3K4me3 and H3K27ac, and 1000 for TFs or DNase-seq. Default: 2000

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

Inputs

ID Type Title Doc
alias String Experiment short name/Alias
window_size Integer (Optional) Window size (2000 is recommended for sharp histone marks like H3K4me3 and H3K27ac)

Window size to count reads and calculate read densities. 2000 is recommended for sharp histone marks like H3K4me3 and H3K27ac, and 1000 for TFs or DNase-seq. Default: 2000

promoter_dist Integer (Optional) Promoter distance, bp

Max distance from gene TSS (in both direction) overlapping which the peak will be assigned to the promoter region. Default: 1000 bp

upstream_dist Integer (Optional) Upstream distance, bp

Max distance from the promoter (only in upstream direction) overlapping which the peak will be assigned to the upstream region. Default: 20,000 bp

bam_file_first File [BAM] BAM file from sample 1

BAM alignment file from sample 1

m_value_cutoff Float (Optional) M-value (log2-ratio) cutoff

Absolute M-value (log2-ratio) cutoff to define biased (differential binding) peaks. Default: 1.0

p_value_cutoff Float (Optional) P-value cutoff

P-value cutoff to define biased peaks. Default: 0.01

annotation_file File [TSV] Annotation file

Tab-separated annotation file

bam_file_second File [BAM] BAM file from sample 2

BAM alignment file from sample 2

peak_file_first File [xls] ChIP-Seq SE sample 1

XLS peak file from sample 1, formatted as MACS2 output

peak_file_second File [xls] ChIP-Seq SE sample 2

XLS peak file from sample 2, formatted as MACS2 output

shift_size_first Integer (Optional) Reads shift size for sample 1

Reads shift size of sample 1. This value is used to shift reads towards 3' direction to determine the precise binding site. Set as half of the fragment length. Default 100

shift_size_second Integer (Optional) Reads shift size for sample 2

Reads shift size of sample 2. This value is used to shift reads towards 5' direction to determine the precise binding site. Set as half of the fragment length. Default 100

broad_peak_file_first File (Optional) [ENCODE broad peak format] ChIP-Seq SE sample 1

Broad peak file from sample 1

broad_peak_file_second File (Optional) [ENCODE broad peak format] ChIP-Seq SE sample 2

Broad peak file from sample 2

Steps

ID Runs Label Doc
manorm
../tools/manorm.cwl (CommandLineTool)
MAnorm - quantitative comparison of ChIP-Seq data

--wa argument is skipped

assign_genes
../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.

filter_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

restore_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

Outputs

ID Type Label Doc
common_peak_file File [TSV] MAnorm common peak file with assigned genes

\"File contains nearest gene information, the M-A values and normalized read density of each peak, common peaks from two samples are merged together. Coordinates in a result file is under 1-based coordinate-system\"

a_values_wig_file File [BED] Genome track file for A-values

Genome track file for A-values

m_values_wig_file File [BED] Genome track file for M-values

Genome track file for M-values

manorm_stderr_log File [Textual format] MAnorm stderr log

MAnorm stderr log

manorm_stdout_log File [Textual format] MAnorm stdout log

MAnorm stdout log

p_values_wig_file File [BED] Genome track file for P-values

Genome track file for P-values

unbiased_peak_file File [BED] Unbiased peak file

Unbiased peak file

ma_with_P_value_plot File [PNG] MA-values with P-values plot

MA-values with P-values plot

above_m_cutoff_peak_file File [BED] Above M-value cutoff peak file

Above M-value cutoff peak file

below_m_cutoff_peak_file File [BED] Below M-value cutoff peak file

Below M-value cutoff peak file

ma_after_normalization_plot File [PNG] MA-values after normalization plot

MA-values after normalization plot

ma_before_normalization_plot File [PNG] MA-values before normalization plot

MA-values before normalization plot

read_density_on_common_peaks_plot File [PNG] Read density on common peaks plot

Read density on common peaks plot

Permalink: https://w3id.org/cwl/view/git/8a92669a566589d80fde9d151054ffc220ed4ddd/workflows/manorm-se.cwl