DiBiG
ICBR BioinformaticsPowered by Actor, v1.0

ChIPseq - Alignment and peak finding

Title: NSD2_CS_RCHACV_070218_Mohammad
Project: https://basespace.illumina.com/run/96022930/UTX_CS_ARP_HA_KARPAS_061218_DAPHNE
Started on: 8/29/2018 11:00:27
Hostname: login1.ufhpc
Run directory: /ufrc/riva/ariva/collab/licht/runs/NSD2-E1099K-Project/NSD2_CS_RCHACV_070218_Mohammad/NSD2_CS_RCHACV_070218_Mohammad
Configuration chipseq.conf
1. Input data
The following table summarizes the samples, conditions, and contrasts in this analysis. A readset is either a single fastq file or a pair of fastq files (for paired-end sequencing).

CategoryData
Summary of input data
Contrasts:2C-EZH2 vs. 9B-EZH2, 2C-H2AK119ub vs. 9B-H2AK119ub, 2C-H3K27ac vs. 9B-H3K27ac, 2C-H3K27me3 vs. 9B-H3K27me3, 2C-H3K36me2 vs. 9B-H3K36me2, 2C-NSD2 vs. 9B-NSD2
Conditions:2C-EZH2, 2C-H2AK119ub, 2C-H3K27ac, 2C-H3K27me3, 2C-H3K36me2, 2C-NSD2, 9B-EZH2, 9B-H2AK119ub, 9B-H3K27ac, 9B-H3K27me3, 9B-H3K36me2, 9B-NSD2
2C-EZH2RCHACV-2C-EZH2 (4 readsets)
2C-H2AK119ubRCHACV-2C-H2AK119ub (4 readsets)
2C-H3K27acRCHACV-2C-H3K27ac (4 readsets)
2C-H3K27me3RCHACV-2C-H3K27me3 (4 readsets)
2C-H3K36me2RCHACV-2C-H3K36me2 (4 readsets)
2C-NSD2RCHACV-2C-NSD2 (4 readsets)
9B-EZH2RCHACV-9B-EZH2 (4 readsets)
9B-H2AK119ubRCHACV-9B-H2AK119ub (4 readsets)
9B-H3K27acRCHACV-9B-H3K27ac (4 readsets)
9B-H3K27me3RCHACV-9B-H3K27me3 (4 readsets)
9B-H3K36me2RCHACV-9B-H3K36me2 (4 readsets)
9B-NSD2RCHACV-9B-NSD2 (4 readsets)
Table 1. Summary of input data

2. Trimming and quality control
The input sequences were trimmed using trimmomatic. Quality control was performed before and after trimming using FastQC. The following table provides links to the quality control reports before and after trimming, as well as the number of reads in the trimmed files.

SampleReadsetReads before trimQC before trimReads after trimQC after trim% Retained
RCHACV-2C-EZH2RCHACV-2C-EZH2_r111,686,727RCHACV-2C-EZH2_S0_L001_R1_001
RCHACV-2C-EZH2_S0_L001_R2_001
11,116,569RCHACV-2C-EZH2_S0_L001_R1_001.trim.paired
RCHACV-2C-EZH2_S0_L001_R2_001.trim.paired
95.12%
RCHACV-2C-EZH2_r211,983,230RCHACV-2C-EZH2_S0_L002_R1_001
RCHACV-2C-EZH2_S0_L002_R2_001
11,385,306RCHACV-2C-EZH2_S0_L002_R1_001.trim.paired
RCHACV-2C-EZH2_S0_L002_R2_001.trim.paired
95.01%
RCHACV-2C-EZH2_r311,642,324RCHACV-2C-EZH2_S0_L003_R1_001
RCHACV-2C-EZH2_S0_L003_R2_001
11,127,542RCHACV-2C-EZH2_S0_L003_R1_001.trim.paired
RCHACV-2C-EZH2_S0_L003_R2_001.trim.paired
95.58%
RCHACV-2C-EZH2_r412,066,534RCHACV-2C-EZH2_S0_L004_R1_001
RCHACV-2C-EZH2_S0_L004_R2_001
11,462,389RCHACV-2C-EZH2_S0_L004_R1_001.trim.paired
RCHACV-2C-EZH2_S0_L004_R2_001.trim.paired
94.99%
RCHACV-2C-H2AK119ubRCHACV-2C-H2AK119ub_r111,129,932RCHACV-2C-H2AK119ub_S0_L001_R1_001
RCHACV-2C-H2AK119ub_S0_L001_R2_001
10,589,134RCHACV-2C-H2AK119ub_S0_L001_R1_001.trim.paired
RCHACV-2C-H2AK119ub_S0_L001_R2_001.trim.paired
95.14%
RCHACV-2C-H2AK119ub_r211,449,589RCHACV-2C-H2AK119ub_S0_L002_R1_001
RCHACV-2C-H2AK119ub_S0_L002_R2_001
10,887,687RCHACV-2C-H2AK119ub_S0_L002_R1_001.trim.paired
RCHACV-2C-H2AK119ub_S0_L002_R2_001.trim.paired
95.09%
RCHACV-2C-H2AK119ub_r311,102,220RCHACV-2C-H2AK119ub_S0_L003_R1_001
RCHACV-2C-H2AK119ub_S0_L003_R2_001
10,617,877RCHACV-2C-H2AK119ub_S0_L003_R1_001.trim.paired
RCHACV-2C-H2AK119ub_S0_L003_R2_001.trim.paired
95.64%
RCHACV-2C-H2AK119ub_r411,556,088RCHACV-2C-H2AK119ub_S0_L004_R1_001
RCHACV-2C-H2AK119ub_S0_L004_R2_001
10,981,933RCHACV-2C-H2AK119ub_S0_L004_R1_001.trim.paired
RCHACV-2C-H2AK119ub_S0_L004_R2_001.trim.paired
95.03%
RCHACV-2C-H3K27acRCHACV-2C-H3K27ac_r17,474,397RCHACV-2C-H3K27ac_S0_L001_R1_001
RCHACV-2C-H3K27ac_S0_L001_R2_001
7,140,689RCHACV-2C-H3K27ac_S0_L001_R1_001.trim.paired
RCHACV-2C-H3K27ac_S0_L001_R2_001.trim.paired
95.54%
RCHACV-2C-H3K27ac_r27,633,201RCHACV-2C-H3K27ac_S0_L002_R1_001
RCHACV-2C-H3K27ac_S0_L002_R2_001
7,288,641RCHACV-2C-H3K27ac_S0_L002_R1_001.trim.paired
RCHACV-2C-H3K27ac_S0_L002_R2_001.trim.paired
95.49%
RCHACV-2C-H3K27ac_r37,423,877RCHACV-2C-H3K27ac_S0_L003_R1_001
RCHACV-2C-H3K27ac_S0_L003_R2_001
7,123,365RCHACV-2C-H3K27ac_S0_L003_R1_001.trim.paired
RCHACV-2C-H3K27ac_S0_L003_R2_001.trim.paired
95.95%
RCHACV-2C-H3K27ac_r47,666,536RCHACV-2C-H3K27ac_S0_L004_R1_001
RCHACV-2C-H3K27ac_S0_L004_R2_001
7,320,015RCHACV-2C-H3K27ac_S0_L004_R1_001.trim.paired
RCHACV-2C-H3K27ac_S0_L004_R2_001.trim.paired
95.48%
RCHACV-2C-H3K27me3RCHACV-2C-H3K27me3_r110,073,957RCHACV-2C-H3K27me3_S0_L001_R1_001
RCHACV-2C-H3K27me3_S0_L001_R2_001
9,604,658RCHACV-2C-H3K27me3_S0_L001_R1_001.trim.paired
RCHACV-2C-H3K27me3_S0_L001_R2_001.trim.paired
95.34%
RCHACV-2C-H3K27me3_r210,302,342RCHACV-2C-H3K27me3_S0_L002_R1_001
RCHACV-2C-H3K27me3_S0_L002_R2_001
9,814,367RCHACV-2C-H3K27me3_S0_L002_R1_001.trim.paired
RCHACV-2C-H3K27me3_S0_L002_R2_001.trim.paired
95.26%
RCHACV-2C-H3K27me3_r310,013,431RCHACV-2C-H3K27me3_S0_L003_R1_001
RCHACV-2C-H3K27me3_S0_L003_R2_001
9,590,492RCHACV-2C-H3K27me3_S0_L003_R1_001.trim.paired
RCHACV-2C-H3K27me3_S0_L003_R2_001.trim.paired
95.78%
RCHACV-2C-H3K27me3_r410,455,813RCHACV-2C-H3K27me3_S0_L004_R1_001
RCHACV-2C-H3K27me3_S0_L004_R2_001
9,946,969RCHACV-2C-H3K27me3_S0_L004_R1_001.trim.paired
RCHACV-2C-H3K27me3_S0_L004_R2_001.trim.paired
95.13%
RCHACV-2C-H3K36me2RCHACV-2C-H3K36me2_r17,765,767RCHACV-2C-H3K36me2_S0_L001_R1_001
RCHACV-2C-H3K36me2_S0_L001_R2_001
7,414,115RCHACV-2C-H3K36me2_S0_L001_R1_001.trim.paired
RCHACV-2C-H3K36me2_S0_L001_R2_001.trim.paired
95.47%
RCHACV-2C-H3K36me2_r27,979,659RCHACV-2C-H3K36me2_S0_L002_R1_001
RCHACV-2C-H3K36me2_S0_L002_R2_001
7,618,550RCHACV-2C-H3K36me2_S0_L002_R1_001.trim.paired
RCHACV-2C-H3K36me2_S0_L002_R2_001.trim.paired
95.47%
RCHACV-2C-H3K36me2_r37,727,634RCHACV-2C-H3K36me2_S0_L003_R1_001
RCHACV-2C-H3K36me2_S0_L003_R2_001
7,410,445RCHACV-2C-H3K36me2_S0_L003_R1_001.trim.paired
RCHACV-2C-H3K36me2_S0_L003_R2_001.trim.paired
95.90%
RCHACV-2C-H3K36me2_r48,046,533RCHACV-2C-H3K36me2_S0_L004_R1_001
RCHACV-2C-H3K36me2_S0_L004_R2_001
7,677,944RCHACV-2C-H3K36me2_S0_L004_R1_001.trim.paired
RCHACV-2C-H3K36me2_S0_L004_R2_001.trim.paired
95.42%
RCHACV-2C-NSD2RCHACV-2C-NSD2_r18,947,748RCHACV-2C-NSD2_S0_L001_R1_001
RCHACV-2C-NSD2_S0_L001_R2_001
8,544,109RCHACV-2C-NSD2_S0_L001_R1_001.trim.paired
RCHACV-2C-NSD2_S0_L001_R2_001.trim.paired
95.49%
RCHACV-2C-NSD2_r29,107,494RCHACV-2C-NSD2_S0_L002_R1_001
RCHACV-2C-NSD2_S0_L002_R2_001
8,705,655RCHACV-2C-NSD2_S0_L002_R1_001.trim.paired
RCHACV-2C-NSD2_S0_L002_R2_001.trim.paired
95.59%
RCHACV-2C-NSD2_r38,891,732RCHACV-2C-NSD2_S0_L003_R1_001
RCHACV-2C-NSD2_S0_L003_R2_001
8,532,075RCHACV-2C-NSD2_S0_L003_R1_001.trim.paired
RCHACV-2C-NSD2_S0_L003_R2_001.trim.paired
95.96%
RCHACV-2C-NSD2_r49,143,621RCHACV-2C-NSD2_S0_L004_R1_001
RCHACV-2C-NSD2_S0_L004_R2_001
8,735,662RCHACV-2C-NSD2_S0_L004_R1_001.trim.paired
RCHACV-2C-NSD2_S0_L004_R2_001.trim.paired
95.54%
RCHACV-9B-EZH2RCHACV-9B-EZH2_r17,150,569RCHACV-9B-EZH2_S0_L001_R1_001
RCHACV-9B-EZH2_S0_L001_R2_001
6,795,885RCHACV-9B-EZH2_S0_L001_R1_001.trim.paired
RCHACV-9B-EZH2_S0_L001_R2_001.trim.paired
95.04%
RCHACV-9B-EZH2_r27,346,761RCHACV-9B-EZH2_S0_L002_R1_001
RCHACV-9B-EZH2_S0_L002_R2_001
6,988,604RCHACV-9B-EZH2_S0_L002_R1_001.trim.paired
RCHACV-9B-EZH2_S0_L002_R2_001.trim.paired
95.12%
RCHACV-9B-EZH2_r37,145,735RCHACV-9B-EZH2_S0_L003_R1_001
RCHACV-9B-EZH2_S0_L003_R2_001
6,825,728RCHACV-9B-EZH2_S0_L003_R1_001.trim.paired
RCHACV-9B-EZH2_S0_L003_R2_001.trim.paired
95.52%
RCHACV-9B-EZH2_r47,372,435RCHACV-9B-EZH2_S0_L004_R1_001
RCHACV-9B-EZH2_S0_L004_R2_001
7,011,792RCHACV-9B-EZH2_S0_L004_R1_001.trim.paired
RCHACV-9B-EZH2_S0_L004_R2_001.trim.paired
95.11%
RCHACV-9B-H2AK119ubRCHACV-9B-H2AK119ub_r110,761,036RCHACV-9B-H2AK119ub_S0_L001_R1_001
RCHACV-9B-H2AK119ub_S0_L001_R2_001
10,179,322RCHACV-9B-H2AK119ub_S0_L001_R1_001.trim.paired
RCHACV-9B-H2AK119ub_S0_L001_R2_001.trim.paired
94.59%
RCHACV-9B-H2AK119ub_r211,162,986RCHACV-9B-H2AK119ub_S0_L002_R1_001
RCHACV-9B-H2AK119ub_S0_L002_R2_001
10,558,232RCHACV-9B-H2AK119ub_S0_L002_R1_001.trim.paired
RCHACV-9B-H2AK119ub_S0_L002_R2_001.trim.paired
94.58%
RCHACV-9B-H2AK119ub_r310,791,781RCHACV-9B-H2AK119ub_S0_L003_R1_001
RCHACV-9B-H2AK119ub_S0_L003_R2_001
10,266,409RCHACV-9B-H2AK119ub_S0_L003_R1_001.trim.paired
RCHACV-9B-H2AK119ub_S0_L003_R2_001.trim.paired
95.13%
RCHACV-9B-H2AK119ub_r411,250,736RCHACV-9B-H2AK119ub_S0_L004_R1_001
RCHACV-9B-H2AK119ub_S0_L004_R2_001
10,625,105RCHACV-9B-H2AK119ub_S0_L004_R1_001.trim.paired
RCHACV-9B-H2AK119ub_S0_L004_R2_001.trim.paired
94.44%
RCHACV-9B-H3K27acRCHACV-9B-H3K27ac_r18,010,614RCHACV-9B-H3K27ac_S0_L001_R1_001
RCHACV-9B-H3K27ac_S0_L001_R2_001
7,654,929RCHACV-9B-H3K27ac_S0_L001_R1_001.trim.paired
RCHACV-9B-H3K27ac_S0_L001_R2_001.trim.paired
95.56%
RCHACV-9B-H3K27ac_r28,166,314RCHACV-9B-H3K27ac_S0_L002_R1_001
RCHACV-9B-H3K27ac_S0_L002_R2_001
7,789,946RCHACV-9B-H3K27ac_S0_L002_R1_001.trim.paired
RCHACV-9B-H3K27ac_S0_L002_R2_001.trim.paired
95.39%
RCHACV-9B-H3K27ac_r37,950,796RCHACV-9B-H3K27ac_S0_L003_R1_001
RCHACV-9B-H3K27ac_S0_L003_R2_001
7,626,533RCHACV-9B-H3K27ac_S0_L003_R1_001.trim.paired
RCHACV-9B-H3K27ac_S0_L003_R2_001.trim.paired
95.92%
RCHACV-9B-H3K27ac_r48,248,803RCHACV-9B-H3K27ac_S0_L004_R1_001
RCHACV-9B-H3K27ac_S0_L004_R2_001
7,869,517RCHACV-9B-H3K27ac_S0_L004_R1_001.trim.paired
RCHACV-9B-H3K27ac_S0_L004_R2_001.trim.paired
95.40%
RCHACV-9B-H3K27me3RCHACV-9B-H3K27me3_r18,787,859RCHACV-9B-H3K27me3_S0_L001_R1_001
RCHACV-9B-H3K27me3_S0_L001_R2_001
8,357,344RCHACV-9B-H3K27me3_S0_L001_R1_001.trim.paired
RCHACV-9B-H3K27me3_S0_L001_R2_001.trim.paired
95.10%
RCHACV-9B-H3K27me3_r28,750,966RCHACV-9B-H3K27me3_S0_L002_R1_001
RCHACV-9B-H3K27me3_S0_L002_R2_001
8,306,023RCHACV-9B-H3K27me3_S0_L002_R1_001.trim.paired
RCHACV-9B-H3K27me3_S0_L002_R2_001.trim.paired
94.92%
RCHACV-9B-H3K27me3_r38,519,822RCHACV-9B-H3K27me3_S0_L003_R1_001
RCHACV-9B-H3K27me3_S0_L003_R2_001
8,136,672RCHACV-9B-H3K27me3_S0_L003_R1_001.trim.paired
RCHACV-9B-H3K27me3_S0_L003_R2_001.trim.paired
95.50%
RCHACV-9B-H3K27me3_r48,850,690RCHACV-9B-H3K27me3_S0_L004_R1_001
RCHACV-9B-H3K27me3_S0_L004_R2_001
8,399,248RCHACV-9B-H3K27me3_S0_L004_R1_001.trim.paired
RCHACV-9B-H3K27me3_S0_L004_R2_001.trim.paired
94.90%
RCHACV-9B-H3K36me2RCHACV-9B-H3K36me2_r16,959,779RCHACV-9B-H3K36me2_S0_L001_R1_001
RCHACV-9B-H3K36me2_S0_L001_R2_001
6,616,531RCHACV-9B-H3K36me2_S0_L001_R1_001.trim.paired
RCHACV-9B-H3K36me2_S0_L001_R2_001.trim.paired
95.07%
RCHACV-9B-H3K36me2_r27,159,359RCHACV-9B-H3K36me2_S0_L002_R1_001
RCHACV-9B-H3K36me2_S0_L002_R2_001
6,803,214RCHACV-9B-H3K36me2_S0_L002_R1_001.trim.paired
RCHACV-9B-H3K36me2_S0_L002_R2_001.trim.paired
95.03%
RCHACV-9B-H3K36me2_r36,954,174RCHACV-9B-H3K36me2_S0_L003_R1_001
RCHACV-9B-H3K36me2_S0_L003_R2_001
6,644,734RCHACV-9B-H3K36me2_S0_L003_R1_001.trim.paired
RCHACV-9B-H3K36me2_S0_L003_R2_001.trim.paired
95.55%
RCHACV-9B-H3K36me2_r47,213,391RCHACV-9B-H3K36me2_S0_L004_R1_001
RCHACV-9B-H3K36me2_S0_L004_R2_001
6,850,475RCHACV-9B-H3K36me2_S0_L004_R1_001.trim.paired
RCHACV-9B-H3K36me2_S0_L004_R2_001.trim.paired
94.97%
RCHACV-9B-NSD2RCHACV-9B-NSD2_r19,390,114RCHACV-9B-NSD2_S0_L001_R1_001
RCHACV-9B-NSD2_S0_L001_R2_001
8,977,382RCHACV-9B-NSD2_S0_L001_R1_001.trim.paired
RCHACV-9B-NSD2_S0_L001_R2_001.trim.paired
95.60%
RCHACV-9B-NSD2_r29,533,900RCHACV-9B-NSD2_S0_L002_R1_001
RCHACV-9B-NSD2_S0_L002_R2_001
9,109,524RCHACV-9B-NSD2_S0_L002_R1_001.trim.paired
RCHACV-9B-NSD2_S0_L002_R2_001.trim.paired
95.55%
RCHACV-9B-NSD2_r39,328,941RCHACV-9B-NSD2_S0_L003_R1_001
RCHACV-9B-NSD2_S0_L003_R2_001
8,958,012RCHACV-9B-NSD2_S0_L003_R1_001.trim.paired
RCHACV-9B-NSD2_S0_L003_R2_001.trim.paired
96.02%
RCHACV-9B-NSD2_r49,602,825RCHACV-9B-NSD2_S0_L004_R1_001
RCHACV-9B-NSD2_S0_L004_R2_001
9,173,005RCHACV-9B-NSD2_S0_L004_R1_001.trim.paired
RCHACV-9B-NSD2_S0_L004_R2_001.trim.paired
95.52%
Table 2. Number of reads in input files and links to QC reports.

The following two tables report the number of reads before and after QC in each sample and in each condition.

SampleReads before QCReads after QC% Retained
RCHACV-2C-EZH247,378,81545,091,80695.17%
RCHACV-2C-H2AK119ub45,237,82943,076,63195.22%
RCHACV-2C-H3K27ac30,198,01128,872,71095.61%
RCHACV-2C-H3K27me340,845,54338,956,48695.38%
RCHACV-2C-H3K36me231,519,59330,121,05495.56%
RCHACV-2C-NSD236,090,59534,517,50195.64%
RCHACV-9B-EZH229,015,50027,622,00995.20%
RCHACV-9B-H2AK119ub43,966,53941,629,06894.68%
RCHACV-9B-H3K27ac32,376,52730,940,92595.57%
RCHACV-9B-H3K27me334,909,33733,199,28795.10%
RCHACV-9B-H3K36me228,286,70326,914,95495.15%
RCHACV-9B-NSD237,855,78036,217,92395.67%
Table 3. Number of reads in each sample before and after QC.



ConditionReads before QCReads after QC% Retained
2C-EZH247,378,81545,091,80695.17%
2C-H2AK119ub45,237,82943,076,63195.22%
2C-H3K27ac30,198,01128,872,71095.61%
2C-H3K27me340,845,54338,956,48695.38%
2C-H3K36me231,519,59330,121,05495.56%
2C-NSD236,090,59534,517,50195.64%
9B-EZH229,015,50027,622,00995.20%
9B-H2AK119ub43,966,53941,629,06894.68%
9B-H3K27ac32,376,52730,940,92595.57%
9B-H3K27me334,909,33733,199,28795.10%
9B-H3K36me228,286,70326,914,95495.15%
9B-NSD237,855,78036,217,92395.67%
Table 4. Number of reads in each condition before and after QC.

3. Alignment to genome
The input sequences were aligned to the genome using Bowtie 2.3.4.2. The following table reports the number of aligned reads for each sample. The WIG files can be uploaded to the UCSC Genome Browser as custom tracks.

SampleTotal readsAligned readsConcordant alignment rateBowtie2 report
RCHACV-2C-EZH245,091,80624,226,98153.73%RCHACV-2C-EZH2.bt2stats.html
RCHACV-2C-H2AK119ub43,076,63123,731,97855.09%RCHACV-2C-H2AK119ub.bt2stats.html
RCHACV-2C-H3K27ac28,872,71014,393,76149.85%RCHACV-2C-H3K27ac.bt2stats.html
RCHACV-2C-H3K27me338,956,48617,256,95344.30%RCHACV-2C-H3K27me3.bt2stats.html
RCHACV-2C-H3K36me230,121,05416,659,39855.31%RCHACV-2C-H3K36me2.bt2stats.html
RCHACV-2C-NSD234,517,50118,784,99954.42%RCHACV-2C-NSD2.bt2stats.html
RCHACV-9B-EZH227,622,00916,167,07058.53%RCHACV-9B-EZH2.bt2stats.html
RCHACV-9B-H2AK119ub41,629,06828,199,83067.74%RCHACV-9B-H2AK119ub.bt2stats.html
RCHACV-9B-H3K27ac30,940,92512,105,50639.12%RCHACV-9B-H3K27ac.bt2stats.html
RCHACV-9B-H3K27me333,199,28710,739,47732.35%RCHACV-9B-H3K27me3.bt2stats.html
RCHACV-9B-H3K36me226,914,95416,013,11359.50%RCHACV-9B-H3K36me2.bt2stats.html
RCHACV-9B-NSD236,217,92315,329,19542.32%RCHACV-9B-NSD2.bt2stats.html
Table 5. Number of alignments to genome.

4. Genome coverage
The following table reports the overall and effective genome coverage in each sample. The Total nt column reports the total number of nucleotides sequenced, i.e. the number of aligned reads times the length of each read. Coverage is this number divided by the size of the genome. Effective bp reports the number of bases in the genome having coverage greater than 5, and the Effective Perc column shows what percentage this is of the genome size. Note that, especially in the case of RNA-seq, the effective genome size may be much smaller than the full size. Eff Coverage is the average coverage over the effectively covered fraction of the genome.

NameTotal ntCoverageEffective bpEffective PercEff Coverage
RCHACV-2C-EZH2298,967,2150.1028,902,0440.90%10.34
RCHACV-2C-H2AK119ub477,661,5870.1557,666,5141.90%8.28
RCHACV-2C-H3K27ac261,058,6240.0828,331,5800.90%9.21
RCHACV-2C-H3K27me3424,097,5810.1430,447,4921.00%13.93
RCHACV-2C-H3K36me2196,419,7310.0622,256,2290.70%8.83
RCHACV-2C-NSD2263,725,0740.0929,353,8511.00%8.98
RCHACV-9B-EZH2195,978,7770.0618,599,9240.60%10.54
RCHACV-9B-H2AK119ub585,385,2770.1971,836,4452.30%8.15
RCHACV-9B-H3K27ac199,797,4780.0621,796,5070.70%9.17
RCHACV-9B-H3K27me3372,071,4800.1219,437,1350.60%19.14
RCHACV-9B-H3K36me2140,910,2840.0514,832,6510.50%9.50
RCHACV-9B-NSD2175,585,5360.0618,894,6540.60%9.29
Table 6. Genome coverage by sample.

The following table reports the overall and effective genome coverage in each condition.

NameTotal ntCoverageEffective bpEffective PercEff Coverage
2C-EZH2298,967,2150.1028,902,0440.90%10.34
2C-H2AK119ub477,661,5870.1557,666,5141.90%8.28
2C-H3K27ac261,058,6240.0828,331,5800.90%9.21
2C-H3K27me3424,097,5810.1430,447,4921.00%13.93
2C-H3K36me2196,419,7310.0622,256,2290.70%8.83
2C-NSD2263,725,0740.0929,353,8511.00%8.98
9B-EZH2195,978,7770.0618,599,9240.60%10.54
9B-H2AK119ub585,385,2770.1971,836,4452.30%8.15
9B-H3K27ac199,797,4780.0621,796,5070.70%9.17
9B-H3K27me3372,071,4800.1219,437,1350.60%19.14
9B-H3K36me2140,910,2840.0514,832,6510.50%9.50
9B-NSD2175,585,5360.0618,894,6540.60%9.29
Table 7. Genome coverage by condition

File: NSD2_CS_RCHACV_070218_Mohammad.sample.cov.xlsx
Size: 30.20 kB
Description: Per-chromosome coverage data, by sample.

File: NSD2_CS_RCHACV_070218_Mohammad.cond.cov.xlsx
Size: 30.18 kB
Description: Per-chromosome coverage data, by condition.

5. Peak detection
Peak detection was performed using MACS version 2.1.1.20160309. The following table provides links to the Pileup, narrowPeaks, and Summits files for each condition. All files are in bedGraph format.

ConditionPileupPeaksSummits
2C-EZH22C-EZH2.macs/2C-EZH2.bedGraph2C-EZH2.macs/2C-EZH2.npeaks.bedGraph2C-EZH2.macs/2C-EZH2.summits.bedGraph
2C-H2AK119ub2C-H2AK119ub.macs/2C-H2AK119ub.bedGraph2C-H2AK119ub.macs/2C-H2AK119ub.npeaks.bedGraph2C-H2AK119ub.macs/2C-H2AK119ub.summits.bedGraph
2C-H3K27ac2C-H3K27ac.macs/2C-H3K27ac.bedGraph2C-H3K27ac.macs/2C-H3K27ac.npeaks.bedGraph2C-H3K27ac.macs/2C-H3K27ac.summits.bedGraph
2C-H3K27me32C-H3K27me3.macs/2C-H3K27me3.bedGraph2C-H3K27me3.macs/2C-H3K27me3.npeaks.bedGraph2C-H3K27me3.macs/2C-H3K27me3.summits.bedGraph
2C-H3K36me22C-H3K36me2.macs/2C-H3K36me2.bedGraph2C-H3K36me2.macs/2C-H3K36me2.npeaks.bedGraph2C-H3K36me2.macs/2C-H3K36me2.summits.bedGraph
2C-NSD22C-NSD2.macs/2C-NSD2.bedGraph2C-NSD2.macs/2C-NSD2.npeaks.bedGraph2C-NSD2.macs/2C-NSD2.summits.bedGraph
9B-EZH29B-EZH2.macs/9B-EZH2.bedGraph9B-EZH2.macs/9B-EZH2.npeaks.bedGraph9B-EZH2.macs/9B-EZH2.summits.bedGraph
9B-H2AK119ub9B-H2AK119ub.macs/9B-H2AK119ub.bedGraph9B-H2AK119ub.macs/9B-H2AK119ub.npeaks.bedGraph9B-H2AK119ub.macs/9B-H2AK119ub.summits.bedGraph
9B-H3K27ac9B-H3K27ac.macs/9B-H3K27ac.bedGraph9B-H3K27ac.macs/9B-H3K27ac.npeaks.bedGraph9B-H3K27ac.macs/9B-H3K27ac.summits.bedGraph
9B-H3K27me39B-H3K27me3.macs/9B-H3K27me3.bedGraph9B-H3K27me3.macs/9B-H3K27me3.npeaks.bedGraph9B-H3K27me3.macs/9B-H3K27me3.summits.bedGraph
9B-H3K36me29B-H3K36me2.macs/9B-H3K36me2.bedGraph9B-H3K36me2.macs/9B-H3K36me2.npeaks.bedGraph9B-H3K36me2.macs/9B-H3K36me2.summits.bedGraph
9B-NSD29B-NSD2.macs/9B-NSD2.bedGraph9B-NSD2.macs/9B-NSD2.npeaks.bedGraph9B-NSD2.macs/9B-NSD2.summits.bedGraph
Table 8. Results of peak detection with MACS.

The following tables shows the number of peaks found for each condition, and their classification.

ConditionTotalUpstreamExonCodingExonIntronDownstreamIntergenic
2C-EZH22358920.09%7.25%2.95%26.95%7.47%35.30%
2C-H2AK119ub224979.36%3.51%0.81%31.07%7.57%47.68%
2C-H3K27ac6609917.29%7.48%1.78%43.44%7.79%22.23%
2C-H3K27me35122620.28%7.16%2.54%36.91%8.61%24.50%
2C-H3K36me2104954.73%1.53%0.26%25.68%5.53%62.28%
2C-NSD273524.58%1.31%0.19%9.47%4.56%79.90%
9B-EZH22991821.54%8.18%3.53%33.04%9.40%24.31%
9B-H2AK119ub2403611.38%4.07%1.29%30.38%8.65%44.23%
9B-H3K27ac6042519.21%8.63%2.28%42.45%8.33%19.10%
9B-H3K27me36611519.63%7.13%2.76%38.36%9.09%23.03%
9B-H3K36me283885.65%2.10%0.41%18.97%6.59%66.29%
9B-NSD252116.01%1.65%0.23%9.92%6.56%75.63%
Table 9. Classification of peaks in genome regions

6. Peak finding
Peak finding was performed using HOMER. Three detection types were performed: Peaks, Regions, and SuperEnhancers. Please see the HOMER documentation for details. The following table provides links to the bedGraph files for each detection type in each condition, and to an Excel file containing the results of each detection type split by genomic region (e.g. Promoters, Exons, Introns, Intergenic). The Combined file contains the results of each detection type in all conditions (one sheet per condition).

ConditionCountbedGraphClassified
Peaks
2C-EZH2182C-EZH2.tags.d/peaks.bedGraph2C-EZH2-peaks-distr.xlsx
2C-H2AK119ub9382C-H2AK119ub.tags.d/peaks.bedGraph2C-H2AK119ub-peaks-distr.xlsx
2C-H3K27ac1,4982C-H3K27ac.tags.d/peaks.bedGraph2C-H3K27ac-peaks-distr.xlsx
2C-H3K27me31222C-H3K27me3.tags.d/peaks.bedGraph2C-H3K27me3-peaks-distr.xlsx
2C-H3K36me21,8942C-H3K36me2.tags.d/peaks.bedGraph2C-H3K36me2-peaks-distr.xlsx
2C-NSD212C-NSD2.tags.d/peaks.bedGraph2C-NSD2-peaks-distr.xlsx
9B-EZH27889B-EZH2.tags.d/peaks.bedGraph9B-EZH2-peaks-distr.xlsx
9B-H2AK119ub09B-H2AK119ub.tags.d/peaks.bedGraph9B-H2AK119ub-peaks-distr.xlsx
9B-H3K27ac2,3439B-H3K27ac.tags.d/peaks.bedGraph9B-H3K27ac-peaks-distr.xlsx
9B-H3K27me31,3779B-H3K27me3.tags.d/peaks.bedGraph9B-H3K27me3-peaks-distr.xlsx
9B-H3K36me27839B-H3K36me2.tags.d/peaks.bedGraph9B-H3K36me2-peaks-distr.xlsx
9B-NSD22789B-NSD2.tags.d/peaks.bedGraph9B-NSD2-peaks-distr.xlsx
Combined10,040peaks.xlsx 
Regions
2C-EZH2262C-EZH2.tags.d/regions.bedGraph2C-EZH2-regions-distr.xlsx
2C-H2AK119ub572C-H2AK119ub.tags.d/regions.bedGraph2C-H2AK119ub-regions-distr.xlsx
2C-H3K27ac5632C-H3K27ac.tags.d/regions.bedGraph2C-H3K27ac-regions-distr.xlsx
2C-H3K27me3572C-H3K27me3.tags.d/regions.bedGraph2C-H3K27me3-regions-distr.xlsx
2C-H3K36me21,2072C-H3K36me2.tags.d/regions.bedGraph2C-H3K36me2-regions-distr.xlsx
2C-NSD212C-NSD2.tags.d/regions.bedGraph2C-NSD2-regions-distr.xlsx
9B-EZH23569B-EZH2.tags.d/regions.bedGraph9B-EZH2-regions-distr.xlsx
9B-H2AK119ub969B-H2AK119ub.tags.d/regions.bedGraph9B-H2AK119ub-regions-distr.xlsx
9B-H3K27ac1809B-H3K27ac.tags.d/regions.bedGraph9B-H3K27ac-regions-distr.xlsx
9B-H3K27me3469B-H3K27me3.tags.d/regions.bedGraph9B-H3K27me3-regions-distr.xlsx
9B-H3K36me27339B-H3K36me2.tags.d/regions.bedGraph9B-H3K36me2-regions-distr.xlsx
9B-NSD23039B-NSD2.tags.d/regions.bedGraph9B-NSD2-regions-distr.xlsx
Combined3,625regions.xlsx 
Enhancers
2C-EZH262C-EZH2.tags.d/superEnhancers.bedGraph2C-EZH2-enhancers-distr.xlsx
2C-H2AK119ub822C-H2AK119ub.tags.d/superEnhancers.bedGraph2C-H2AK119ub-enhancers-distr.xlsx
2C-H3K27ac942C-H3K27ac.tags.d/superEnhancers.bedGraph2C-H3K27ac-enhancers-distr.xlsx
2C-H3K27me3102C-H3K27me3.tags.d/superEnhancers.bedGraph2C-H3K27me3-enhancers-distr.xlsx
2C-H3K36me2762C-H3K36me2.tags.d/superEnhancers.bedGraph2C-H3K36me2-enhancers-distr.xlsx
2C-NSD212C-NSD2.tags.d/superEnhancers.bedGraph2C-NSD2-enhancers-distr.xlsx
9B-EZH2399B-EZH2.tags.d/superEnhancers.bedGraph9B-EZH2-enhancers-distr.xlsx
9B-H2AK119ub09B-H2AK119ub.tags.d/superEnhancers.bedGraph9B-H2AK119ub-enhancers-distr.xlsx
9B-H3K27ac1889B-H3K27ac.tags.d/superEnhancers.bedGraph9B-H3K27ac-enhancers-distr.xlsx
9B-H3K27me3679B-H3K27me3.tags.d/superEnhancers.bedGraph9B-H3K27me3-enhancers-distr.xlsx
9B-H3K36me2409B-H3K36me2.tags.d/superEnhancers.bedGraph9B-H3K36me2-enhancers-distr.xlsx
9B-NSD2139B-NSD2.tags.d/superEnhancers.bedGraph9B-NSD2-enhancers-distr.xlsx
Combined616enhancers.xlsx 
Table 10. Results of peak detection with HOMER.

The following plots show the distribution of peak locations in the different conditions.

(png format, 37.91 kB)
Distribution of peak locations


(png format, 37.66 kB)
Distribution of region locations


(png format, 37.95 kB)
Distribution of enhancer locations

7. Differential peak finding
Differential peak finding was performed using the HOMER getDifferentialPeaks command. The following files contain the peaks showing significant differences in each contrast.

TestControlTest peaksControl peaksDiffPeaks
2C-EZH29B-EZH207882C-EZH2.vs.9B-EZH2.diffPeaks.xlsx
2C-H2AK119ub9B-H2AK119ub002C-H2AK119ub.vs.9B-H2AK119ub.diffPeaks.xlsx
2C-H3K27ac9B-H3K27ac38510202C-H3K27ac.vs.9B-H3K27ac.diffPeaks.xlsx
2C-H3K27me39B-H3K27me32311942C-H3K27me3.vs.9B-H3K27me3.diffPeaks.xlsx
2C-H3K36me29B-H3K36me24872842C-H3K36me2.vs.9B-H3K36me2.diffPeaks.xlsx
2C-NSD29B-NSD202782C-NSD2.vs.9B-NSD2.diffPeaks.xlsx
Table 11. Results of differential peak detection with HOMER.

8. Motif finding
ChIP-Seq peaks were analyzed with the HOMER findMotifs function. The links in the following table lead to the motif finding report for each condition.

ConditionKnownde novo
2C-EZH22C-EZH2.motifs/knownResults.html2C-EZH2.motifs/homerResults.html
2C-H2AK119ub2C-H2AK119ub.motifs/knownResults.html2C-H2AK119ub.motifs/homerResults.html
2C-H3K27ac2C-H3K27ac.motifs/knownResults.html2C-H3K27ac.motifs/homerResults.html
2C-H3K27me32C-H3K27me3.motifs/knownResults.html2C-H3K27me3.motifs/homerResults.html
2C-H3K36me22C-H3K36me2.motifs/knownResults.html2C-H3K36me2.motifs/homerResults.html
2C-NSD22C-NSD2.motifs/knownResults.html2C-NSD2.motifs/homerResults.html
9B-EZH29B-EZH2.motifs/knownResults.html9B-EZH2.motifs/homerResults.html
9B-H2AK119ub9B-H2AK119ub.motifs/knownResults.html9B-H2AK119ub.motifs/homerResults.html
9B-H3K27ac9B-H3K27ac.motifs/knownResults.html9B-H3K27ac.motifs/homerResults.html
9B-H3K27me39B-H3K27me3.motifs/knownResults.html9B-H3K27me3.motifs/homerResults.html
9B-H3K36me29B-H3K36me2.motifs/knownResults.html9B-H3K36me2.motifs/homerResults.html
9B-NSD29B-NSD2.motifs/knownResults.html9B-NSD2.motifs/homerResults.html
Table 12. Results of motif finding with HOMER.

9. Insert size distribution
Insert size distribution was analyzed using the Picard CollectInsertSizeMetrics tool. The following PDF files contain the plot of insert size distribution in each condition.
File: 2C-EZH2.insert-size-metrics.pdf
Size: 8.52 kB
Description: Insert size distribution in 2C-EZH2

File: 2C-H2AK119ub.insert-size-metrics.pdf
Size: 8.14 kB
Description: Insert size distribution in 2C-H2AK119ub

File: 2C-H3K27ac.insert-size-metrics.pdf
Size: 7.64 kB
Description: Insert size distribution in 2C-H3K27ac

File: 2C-H3K27me3.insert-size-metrics.pdf
Size: 8.19 kB
Description: Insert size distribution in 2C-H3K27me3

File: 2C-H3K36me2.insert-size-metrics.pdf
Size: 8.31 kB
Description: Insert size distribution in 2C-H3K36me2

File: 2C-NSD2.insert-size-metrics.pdf
Size: 8.15 kB
Description: Insert size distribution in 2C-NSD2

File: 9B-EZH2.insert-size-metrics.pdf
Size: 8.65 kB
Description: Insert size distribution in 9B-EZH2

File: 9B-H2AK119ub.insert-size-metrics.pdf
Size: 10.21 kB
Description: Insert size distribution in 9B-H2AK119ub

File: 9B-H3K27ac.insert-size-metrics.pdf
Size: 6.82 kB
Description: Insert size distribution in 9B-H3K27ac

File: 9B-H3K27me3.insert-size-metrics.pdf
Size: 6.88 kB
Description: Insert size distribution in 9B-H3K27me3

File: 9B-H3K36me2.insert-size-metrics.pdf
Size: 8.60 kB
Description: Insert size distribution in 9B-H3K36me2

File: 9B-NSD2.insert-size-metrics.pdf
Size: 7.02 kB
Description: Insert size distribution in 9B-NSD2

10. MultiQC report
MultiQC is a general Quality Control tool for a large number of bioinformatics pipeline. The report on this analysis (generated using MultiQC version 1.1) is available here:

MultiQC report
11. UCSC hub

UCSC Genome Browser: use the previous link to display the data tracks automatically, or copy the the URL http://licht:licht@lichtlab.cancer.ufl.edu/reports/NSD2//NSD2_CS_RCHACV_070218_Mohammad/hub/hub.txt and paste it into the "My Hubs" form in this page.

WashU EpiGenome Browser: use the previous link to display the data tracks automatically, or copy the following URL into the "Datahub by URL Link" field: http://licht:licht@lichtlab.cancer.ufl.edu/reports/NSD2//NSD2_CS_RCHACV_070218_Mohammad/hub/hub.json.



Completed: 8-29-2018@11:01
© 2018, A. Riva, University of Florida.