Deeptools correlation plot: between RNAseq and H3K27ac ChIPseq
Hi, I have been trying to plot correlation between my transcriptome and epigenome (H3K27ac) data but it shows almost no correlation. Within transcriptome and epigenome samples, the correlation is fine but I am not sure deeptools is a good way to plot correlation between RNAseq and H3K27ac ChIPseq. Here is my command:
python3 multiBamSummary bins --bamfiles /*R*bam --smartLabels -o multibam.npz -p 20
python3 plotCorrelation -in multibam.npz --corMethod pearson --skipZeros --whatToPlot scatterplot -o multibam_perason.png --outFileCorMatrix pearson.tab
Link to images:
• 2,758 views
•
link
updated
by
GenoMax
•
written
by
aditi225 •
•
written
by
aditi225 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
Please use these directions: How to add images to a Biostars post
I'll do that, thanks.
Presumably you have MANY bins with no signal, have you tried excluding them?
It remains the same even when I used BED files for locus specific correlations. Thanks
You might scroll around with IGV a bit just to eye-ball whether these values seem reasonable. I haven't looked at this particular correlation myself, so I don't really know what to expect. You may also use the
--log1poption, which won't change the correlation, but will make the graphs easier to see.I would recommend to call H3K27ac peaks, get nearest/overlapping genes with bedtools nearest, and plot the peaks maximum value alongside the associated genes expression with R / Python / Excel.
Otherwise, rank genes from the most expressed to the least expressed and make a heatmap with deepTools plotting H3K27ac at these genes TSS (
computeMatrixfollowed byplotHeatmap) while keeping the expression ranked genes order (--sortRegions keep). I did that kind of analysis to validate RNA-seq / FAIRE-seq correlation (noisy old-fashioned open chromatin signal) and it worked like a charm.