merging Grange result from runFIMO
Entering edit mode
3 months ago
berlin • 0

I have three different groups of gene sets (Control, Treat1 , Treat 2), and for TFBS prediction of each group, I have tried to run FIMO(= Find Individual Motif Occurrences ) for each group from MEME suit. For individual transcription factors, HOCOMOCO data was used with their official TF name.

With result of FIMO, I have acquired multiple grange objects. As there are 763 TF from HOCOMOCO, I will have around 3763 each promotor seq.

So my question will be,

  1. Would be ideal to merge all the dataset from each group ?
  2. If yes, how should I merge this multiple grange object?
  3. For interpretation of fimo result, would it be wise to compare motif occurence?
  4. Is there any better way to compare TFBS prediction from different groups of gene set, only with their name. (not with chip data)
  • targeted sequences are was upstream 5k downstream 3k from TSS

Example of code for fimo :

  dummyseq <- getSeq(BSgenome.Hsapiens.UCSC.hg38, found_PR)
  writeXStringSet(dummyseq, file="DummySeq.fasta")

  foo_dummy<- MotifDb::MotifDb %>%  # Query the database for the HOCOMOCO motif using it's gene name
  MotifDb::query(HOCOMOCO$`Transcription factor`[1]) %>%  # Convert from motifdb format to universalmotif format
  universalmotif::convert_motifs() %>%  
 .[[1]]         # The result is a list, to simplify the object, return it as a single universalmotif

  Foo_dummy_Fimo_result_CT <- runFimo(DummySeq.fasta, foo_dummy)

  class(Foo_dummy_Fimo_result_CT) # grange 

  length(Foo_dummy_Fimo_result_CT) # 1428 
grange TFBS Motif JASPAR HOCOMOCO • 190 views
Entering edit mode

In the end, tried to convert grange data to dataframe with annoGR2DF function and rbind them. It works!


Login before adding your answer.

Traffic: 886 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6