CLIP enrichment analysis
1
0
Entering edit mode
6 weeks ago
noahhelton98 ▴ 10

Hi there, I am analyzing a dataset from another paper and am trying to recreate one their methods on CLIP Seq. I have the bam files and in their methods they claimed:

'For each RNA, we enumerated 100 nucleotide windows across the entire RNA. For each window, we calculated the enrichment by computing the number of reads overlapping the window in the protein elution sample divided by the total number of reads within the protein elution sample. We normalized this ratio by the number of reads in the input sample divided by the total number of reads in the input sample.'

I have two questions with this :

First, when they say they protein elution sample, should I be taking the total number of reads within the entire BAM file? Or just the ones reads within the RNA I am looking at?

Second, what methods are used for normalization of CLIP data? Is there an efficient way to do this in python? I created my sliding window and plots of the reads in python so I would like to stick with it to try and completely recreate this. Unfortunately, the paper didn't link any github or source code - so this is all I really have.

Thanks so much.

ChIP-Seq • 129 views
ADD COMMENT
1
Entering edit mode
6 weeks ago

Its not possible to tell from their description, but I would guess the normalization is applied separately for each RNA. That what I have done when I've done iCLIP.

There is as yet, no standard/consensus on the analysis of PAR/i/eCLIP data, and I suspect the ideal approach depends on the sitaution and goal of analysis.

If its any help, I have written a few python functions for dealing with CLIP-like data. Its not the best documented in the world, and its not published, but if you find it helpful, you can look at www.github.com/sudlab/iCLIPlib .

ADD COMMENT
0
Entering edit mode

Thanks so much I'll take a look! I appreciate it tons.

ADD REPLY

Login before adding your answer.

Traffic: 2336 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6