Question: Strand information when doing peak annotation with HOMER
gravatar for jiab
14 months ago by
jiab20 wrote:

As far as I can tell (by running tests manually), in the HOMER package [1] ignores strand information when searching the Transcription Start Site (TSS) nearest to each peak. After a close inspection of the various settings available, I couldn't identify one which changes this behaviour.

This was unexpected for two reasons:

  • HOMER actually requires input to have a Strand column [2]
  • The data from some NGS technologies are directional (e.g CAGE), so strand information should be used when assigning the peaks to transcript/gene models

Could someone please comment if I'm missing something here?




annotation ngs homer • 913 views
ADD COMMENTlink modified 14 months ago by Lars Juhl Jensen11k • written 14 months ago by jiab20

(A couple more things to bare in mind when searching for an annotation tool:)

PeakAnnotator [1] also ignores strand information (accordingly to its paper [2]) and consider this post [3] for strand-related issues with the ChIPpeakAnno [4] R packages.





ADD REPLYlink modified 14 months ago • written 14 months ago by jiab20
gravatar for Lars Juhl Jensen
14 months ago by
Copenhagen, Denmark
Lars Juhl Jensen11k wrote:

According to the documentation on the second page that you link to yourself, the behavior depends on the -strand argument:

-strand <+|-|both> (Count tags on specific strands relative to peak, default: both)
ADD COMMENTlink written 14 months ago by Lars Juhl Jensen11k

Well, that still ignores per peak strand information actually contained in the input. But it means I have the option of splitting the input based on strand and running HOMER twice with the appropriate -strand setting.

ADD REPLYlink written 14 months ago by jiab20
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: 981 users visited in the last hour