ATACseqTFEA - preparebindingSites() for whole genome?
1
0
Entering edit mode
jfmay • 0
@518122bc
Last seen 12 months ago
Canada

Hi, I'm new to R and struggling to follow the guide for ATACseqTFEA found here . The guide's example is for the Drerio genome chromosome 1, in a limited Granges to scan for binding sites.

However, for my full experiment, I want to scan the whole mouse genome. My problem is that the code I'm running for this purpose takes forever to run (so far has taken ~24 hours and currently still running). Is that very normal for this type of analysis? I understand I'm working with big data, but does a more efficient solution exist? Below is the code I entered. I omitted the arguments specifying chromosome/Granges within the genome in order to include the whole thing.

library(ATACseqTFEA)
library(ATACseqQC)
library(BSgenome.Mmusculus.UCSC.mm10)

motifs <- readRDS(system.file("extdata", "PWMatrixList.rds",
                              package="ATACseqTFEA"))
bindingSites <- 
  prepareBindingSites(motifs, Mmusculus, 
                      p.cutoff = 5e-05)

Any guidance would be much appreciated.

ATACseqTFEA • 692 views
ADD COMMENT
0
Entering edit mode
Ou, Jianhong ★ 1.3k
@ou-jianhong-4539
Last seen 4 days ago
United States

You can import the results from fimo which will be much faster. In background, ATACseqTFEA is using motifmatchr::matchMotifs to predict the binding sites. If multiple TFs invovled (for example more than 600), it will take some time. To import the results from fimo, you just import the results into a GRanges object with a meta column as motif. Let me know if you still have any questions.

Jianhong.

ADD COMMENT
0
Entering edit mode

Hello Jianhong, Thanks for developing this tool. It's really helpful for the analysis of ATAC. I have the same problem as mentioned above. It will take very long time to finish prepare Binding Sites. You also mentioned that the FIMO could be used to speed up this process. Could you please put the codes here to show how to use FIMO for this step. Thanks. Weiming

ADD REPLY
0
Entering edit mode

I will work on this to write a help function for importing FIMO results.

ADD REPLY
0
Entering edit mode

I add a function named as importFimoBindingSites to help to import FIMO results. It is available at jianhong/ATACseqTFEA@devel. And will be available at next release version.

ADD REPLY

Login before adding your answer.

Traffic: 377 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

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

Powered by the version 2.3.6