Search
Question: Strand-aware CSAW analysis
1
gravatar for IV
2.8 years ago by
IV20
USA
IV20 wrote:

Dear all,

I've started working with the fantastic csaw package.

I was searching for options or tweaks to take strand into account, since we have some strand-specific IP library preps performed.

Is there a way to incorporate this into the analysis?

Thank you in advance,

Ioannis

***Edit

We are analyzing stranded library preps such as: MeDIP-Seq, various RBP-IP protocols such as CLIP-Seq (some are stranded), etc.

***/Edit

ADD COMMENTlink modified 2.8 years ago • written 2.8 years ago by IV20
3
gravatar for Aaron Lun
2.8 years ago by
Aaron Lun17k
Cambridge, United Kingdom
Aaron Lun17k wrote:

It depends on what you mean by "strand-aware". Currently, csaw is slightly strand-aware in that directional read extension behaves differently for reads mapped to different strands. Other than that, though, there's not much distinction made between strands.

I've have to know more about the nature of the strand specificity in your experiment, in order to give some better advice. I don't think I've ever dealt with strand-specific IP's before, given that the pulled-down DNA is double-stranded to start with (unless you've got stuff like ChIP-Exo data).

---

Edit: Given that strand-specific counting seems to be pretty useful, I've added a forward option to the readParam object. This'll tell windowCounts (or regionCounts) whether or not to extract and count reads on the forward, reverse or both strands.

I've also added a strandedCounts wrapper function to easily get both forward- and reverse-strand counts for genomic bins or regions. Word of warning, though; the downstream functions in csaw that work with GRanges are not designed to be strand-aware, so some extra work will be necessary if you want, e.g., sensible strand-specific clustering. Check out ?strandedCounts for examples.

These updates should be active in the 1.1.19 build on BioC-devel; give it a day or so to go through.

ADD COMMENTlink modified 2.8 years ago • written 2.8 years ago by Aaron Lun17k

MeDIP-seq gives strand-specific methylation measurements.

ADD REPLYlink written 2.8 years ago by Dario Strbenac1.4k

Okay, fair enough. I'll throw in some strand-specific options in the read counting functions.

ADD REPLYlink written 2.8 years ago by Aaron Lun17k

Thanks Aaron! I'll give it a go.

ADD REPLYlink written 2.8 years ago by IV20
Please log in to add an answer.

Help
Access

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