How it is define in terms if accessibility. For example in expression data when I compare A vs B either A or B as reference we define that genes are getting differentially expressed to A or B with our choice of cotrol.
Similarly in ATAC seq data when I'm comparing two stages how do i define the stage specific peaks? Such as set of peaks are only for A and other for B?
So far what i have done is I get a consensus peaks then from those peaks I get read-count and from those read count I do run deseq2 on those accessible data.
Can I do the same as I do for RNA seq data ,and compare peaks of A vs B using deseq2 and say if a peak is differentially accessible then it can be considered as that condition specific peaks.
There is other way of finding out condition or stage specific peaks ? "Identification of tissue-specific chromatin accessible regions
We used a strategy described previously based on the Shannon entropy to compute a tissue specificity index for each peak4,36,37. Specifically, for each peak, we defined its relative accessibility in a tissue type i as Ri = Ei/ΣE, where Ei is the RPM value for the peak in the tissue i, ΣE is the sum of RPM values in all tissues, and N is the total number of tissues. The entropy score for each peak across tissues can be defined as H = −1 * sum(Ri * log2Ri) (1 < i < N), where the value of H ranges between 0 to log2(N). An entropy score close to zero indicates the accessibility of this peak is highly tissue-specific, while an entropy score close to log2(N) indicates that this peak is ubiquitously accessible38. Based on the distribution of entropy scores, peaks with score less than 3.5 were selected as tissue-restricted peaks."
The above is one of the way they saying tissue specific peaks from this paper
Any suggestion would be really appreciated