Hi all, I have a question about sample grouping when calling peaks using MACS2:
I have data from two studies containing Input and H3K27ac samples from multiple tissues (3 Input-H3K27ac pairs for each tissue type in one study, and some samples from the other study).
What is the best way to group samples for peak calling: combine all samples from one study for peak calling, or sort them by tissue and do peak calling for each tissue within each study? Or do peak calling for each pair of samples separately?
Thank you!
1 answer
Definitely do not mix studies. Batch effects will introduce spurious calls. You could indeed group by study and tissue so you get one set of peaks, for example study1-organ1, study2-organ1, study1-organ2, study2-organ2 etc...and then assess reproducible peaks with something like the Irreproducible Discovery Rate framework. That would give you a list of peaks per organ that is consistent between studies.
I used this as guideline in the past: https://hbctraining.github.io/Intro-to-ChIPseq/lessons/07_handling-replicates-idr.html
HBC has a lot of great resources for training.
Log in to answer this question.
-
Can DiffBind Perform Peak Calling and Differential Analysis Directly from BAM Files?
written by Siqi •
-
Normalizing Peaks in ChIP seq
written by mropri
-
High Variability when Peak Calling with SEACR
written by hunter.a.ogg •
-
Controls for ChIP-seq Peak Calling
written by Pappu
-
MACS2 option for ATAC seq?
written by star
-
Peak calling using macs2?
written by star
-
Metrics for comparing Hi-C and ATAC-seq
written by iMorbit •
-
peaks overlap: ATAC & CHIP seq
written by sophialovechan •
-
Comparing numbers of ChIP-seq peaks between different sample types
written by nash.claire
-
ChIP-Seq Histone modification time span study
written by Sudhir Jadhao
