Question: How to classify chIP-seq reads into gene bins?
gravatar for Zee_S
22 months ago by
Zee_S50 wrote:

Hello everyone!

I am very new to RNA-seq data analysis. Your help on how to accomplish the below task is highly appreciated! I will explain my problem as simple as I can.

1) Based on FeatureCounts output, I have computed the TPM score for my set of genes of interest. Here, I defined feature as "exon".

2) Next, based on TPM scores, I classified my genes into TPM bins (0, 0.1, 1, 10, 100, 1000, >1000). I got the empirical number of genes per each of these bins.

3) Independent of this, I have a Bowtie output in bam format for a chip-seq dataset. I used FeatureCounts once again to assign the reads in this bam file to exons+introns of each gene (I ran feature counts twice, where each time the feature was defined as either "exon" or "intron").

4) Now, I want to normalize these raw exon counts and raw intron counts per gene from the chIP seq dataset in such a way that will allow me to see to which gene bins (from step2) the chip data falls into. The idea is to see how my chip enrichment regions overlap with the genes of various TPM scores.

Can I calculate a "TPM score" for the chip-seq dataset and use this as a proxy for comparing it to the actual TPM scores of the genes? this is probably not the right way.

Thanks a lot in advance for you advice!

ADD COMMENTlink modified 21 months ago by Devon Ryan93k • written 22 months ago by Zee_S50
gravatar for Devon Ryan
21 months ago by
Devon Ryan93k
Freiburg, Germany
Devon Ryan93k wrote:

The simplest way to go about this would be to make normalized (to 1x sequencing depth) bigWig files of each sample with bamCoverage (from deepTools). You can then create BED files for each of your gene groups (the various TPM bins) and use those and the bigWig files with computeMatrix. Finally, either plot a heatmap (plotHeatmap) or a profile (plotProfile) to see how the signals differ in the various groups. You could alternatively use the original BAM files and the aforementioned BED files with plotEnrichment.

ADD COMMENTlink written 21 months ago by Devon Ryan93k
Please log in to add an answer.


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