Workflow: MAnorm PE - quantitative comparison of ChIP-Seq paired-end data
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 PE sample 1** * previously analyzed ChIP-Seq paired-end experiment to be used as Sample 1 **ChIP-Seq PE sample 2** * previously analyzed ChIP-Seq paired-end 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
- 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 PE sample 1 |
XLS peak file from sample 1, formatted as MACS2 output |
peak_file_second | File [xls] | ChIP-Seq PE 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 PE sample 1 |
Broad peak file from sample 1 |
broad_peak_file_second | File (Optional) [ENCODE broad peak format] | ChIP-Seq PE 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) |
|
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 |
https://w3id.org/cwl/view/git/564156a9e1cc7c3679a926c479ba3ae133b1bfd4/workflows/manorm-pe.cwl