Hi all,
I called 90,112 peaks from H3K27ac CHIP-seq data using MACS2, and then identified 1,108 super enhancers (SEs) using ROSE with default parameter '-s 12500 -t 2500'. However, I found that most of SEs are located in promoter regions.
The genomic annotation of 1,108 SEs by ROSE:
The genomic annotation of 90,112 peaks by MACS2:
I am kind of interested in SEs located in distal intergenic regions. But now the number of this type of SEs is limited. Could any one guide me to solve this problem? Any suggestion will be highly appreciated. Thank you in advance.
•
written
by
Erica •
1 answer
Hi, You will need to get the tss information for your genome :
wget -q -O " http://ftp.ensembl.org/pub/grch37/current/gtf/homo_sapiens/Homo_sapiens.GRCh37.87.gtf.gz" | gunzip -c | grep -v "#" ./data/Homo_sapiens.GRCh37.87.gtf | awk '($3=="gene")' | grep protein_coding | awk '{OFS="\t"};{if($7 == "+"){start = $4} else if($7 == "-"){start = $5}};{print $1, start-1, start}' | grep -v "\." | grep -v "_" > tss.bed
With this, you can use awk to extend the TSS as much as you want (here 2000bp each way).
awk '{OFS="\t" print $1,$2-2000,$3+2000}'
Then bedtools intersect :
bedtools intersect -a yourpeakfile -b tss.be -v > filtered_tss.bed
You can check again with Chipseeker if this work
Log in to answer this question.
-
Strategy for noisy ATAC-seq: Fragment size selection for pioneer factor DA analysis
written by Abeer H •
-
Difficulties in calling peaks for CUT&RUN (IgG control with high RPKM value)
written by Basti •
-
Question about an alignment
written by barslmn •
-
Integration of Cut&Tag to ATACseq and RNAseq data - impact of a mutation
written by Alexandre Eeckhoutte •
-
ROSE called super-enhancers overlapping with promoters
written by os306 •
-
deeptools heatmap mostly empty
written by a_bis •
-
Can I use All Peaks (MACS2 -q 1) for analysis?
written by Tian •
Hi Erica, were you able to find a solution? I am trying to do the same. Also, did you use chipseeker to find the feature distribution? I wonder which output file from Rose you use to run this with SE. Thanks!!