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