SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
calculate genome bins' similarity GUZUMO Metagenomics 0 11-06-2017 05:38 AM
MetaCV: a composition-based algorithm to classify metagenomic reads and function whfwind Metagenomics 2 06-22-2016 08:47 AM
Why doesn't pbtranscript.py classify call reads of inserts for films with 2 to 5 rds cklopp Pacific Biosciences 3 03-03-2016 10:44 AM
DEXSeq questions on power and counting bins alittleboy Bioinformatics 0 07-02-2013 11:00 AM
Newbler --- classify many reads as Repeat NSTbioinformatics 454 Pyrosequencing 2 06-10-2010 03:43 AM

Reply
 
Thread Tools
Old 03-26-2018, 07:51 AM   #1
a_sen
Junior Member
 
Location: France

Join Date: Mar 2018
Posts: 1
Default How to classify chIP-seq reads into gene bins?

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!
a_sen is offline   Reply With Quote
Reply

Tags
chip seq, feature counts, gene expression, rnaseq analysis, tpm

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 08:43 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO