[BioC] HTSeq-Count

Pickl, Julia j.pickl at dkfz-heidelberg.de
Mon Sep 8 19:10:38 CEST 2014

Hi Steve,

thank you very much for your help.

Which tool for DE did you use?
I used edger, however I´ve read that edgeR and DESeq2 might be overstringent for RIP-Seq (f.e. RIP-Seeker package Paper, Supplement).

Best wishes,

-----Ursprüngliche Nachricht-----
Von: mailinglist.honeypot at gmail.com [mailto:mailinglist.honeypot at gmail.com] Im Auftrag von Steve Lianoglou
Gesendet: Freitag, 5. September 2014 18:35
An: Julia [guest]
Cc: bioconductor at r-project.org list; Pickl, Julia
Betreff: Re: [BioC] HTSeq-Count

Hi Julia,

On Fri, Sep 5, 2014 at 7:41 AM, Julia [guest] <guest at bioconductor.org> wrote:
> Hi all,
> I am new to the field of seq and performed a RIP-Seq experiment using HTSeq count as counter.
> I get now the following (using union, but doesn´t look better for interesection_strict):
> __no_feature    1503377
> __ambiguous     490772
> __too_low_aQual 0
> __not_aligned   0
> __alignment_not_unique  5277314
> When I sum up counts for all genes, I get 3227845.
> The number for __no_feature, __ambiguous, __alignment_not_unique look very high.
> Does somebody have an idea for that?

While I haven't worked with RIP-seq data myself, I do have some experience with HITS-CLIP and PAR-CLIP, which I believe are quite similar (at least in principle) -- these experiments must incredibly difficult to pull off, however, because I'd say most of these types of datasets that came my way were notoriously/incredibly noisy.

This is just to say the problem you are seeing may not be an informatics problem, and could (quite possibly) be an experimental one.


Steve Lianoglou
Computational Biologist

More information about the Bioconductor mailing list