peakC: a flexible, non-parametric peak calling package for 4C and Capture-C data

High-resolution methods such as 4C and Capture-C enable the study of chromatin loops such as those formed between promoters and enhancers or CTCF/cohesin binding sites. An important aspect of 4C/CapC analyses is the identification of robust peaks in the data for the identification of chromatin loops. Here we present an R package for the analysis of 4C/CapC data. We generated 4C data for 10 viewpoints in 2 tissues in triplicate to test our methods. We developed a non-parametric peak caller based on rank-products. Sampling analysis shows that not read depth but template quality is the most important determinant of success in 4C experiments. By performing peak calling on single experiments we show that the peak calling results are similar to the replicate experiments, but that false positive rates are significantly reduced by performing replicates. 4C-seq was performed with either 3 or 6 replicates on a range of viewpoints in ES 14 embryonic stem cells (IB10) and e14.5 fetal liver cells