Question: Strand information when doing peak annotation with HOMER
gravatar for jiab
10 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 • 614 views
ADD COMMENTlink modified 10 months ago by Lars Juhl Jensen11k • written 10 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 10 months ago • written 10 months ago by jiab20
gravatar for Lars Juhl Jensen
10 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 10 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 10 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: 928 users visited in the last hour