Question: (Closed) Compare Super Enhancers in normal vs breast cancer cells
0
gravatar for aropri
9 weeks ago by
aropri70
aropri70 wrote:

I did ChIP seq using H3K27ac antibody to quantify super enhancer regions in normal vs breast cancer cells. I ran ROSE algorithm on normal mammary epithelial cells and breast cancer cells to get a list of super enhancer regions in each cell line. I have a list of super enhancer regions for both and I want to compare super enhancers between the two outputs I got. I want to compare the signal of the super enhancer regions that were classified by ROSE in normal and breast cancer cells and see calculate fold change of those regions to see if they are increasing or decreasing in intensity. I have a bed files of super enhancer regions for normal cells (MCF10A) and same for breast cancer cells (MCF10A-CA1). I know what I have to do. Briefly, collapse A.SE.bed and B.SE.bed into AllSEs.bed using bedtools merge. Quantify coverage of All SEs.bed using A.bam and B.bam and bedtools intersect. Normalize the read coverage to the sequencing depth, and calculate the fold-change between them. The problem I was having is how can I collapse or merge two or three bed files into one file that will contain all the super enhancer regions classified and quantify coverage? Any help is appreciated. Thank you so much!

chip-seq • 100 views
ADD COMMENTlink written 9 weeks ago by aropri70

Hello aropri!

Questions similar to yours can already be found at:

We have closed your question to allow us to keep similar content in the same thread.

If you disagree with this please tell us why in a reply below. We'll be happy to talk about it.

Cheers!

ADD REPLYlink written 9 weeks ago by ATpoint44k

Please follow up in the previous question.

ADD REPLYlink written 9 weeks ago by ATpoint44k
Please log in to add an answer.
The thread is closed. No new answers may be added.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2456 users visited in the last hour
_