As transcription factors (TFs) play a crucial role in regulating the transcription process through binding on the genome alone or in a combinatorial manner, TF enrichment analysis is an efficient and important procedure to locate the candidate functional TFs from a set of experimentally defined regulatory regions. While it is commonly accepted that structurally related TFs may have similar binding preference to sequences (i.e. motifs) and one TF may have multiple motifs, TF enrichment analysis is much more challenging than motif enrichment analysis. Here we present a R package for TF enrichment analysis which combine motif enrichment with the PECA model.
The package enrichTF is part of Bioconductor project starting from Bioc 3.9 built on R 3.6. To install the latest version of enrichTF, please check your current Bioconductor version and R version first. The latest version of R is recommended, and then you can download and install enrichTF and all its dependencies as follows:
if (!requireNamespace("BiocManager", quietly=TRUE))
install.packages("BiocManager")
BiocManager::install("enrichTF")Similar with other R packages, please load enrichTF each time before using the package.
library(enrichTF)It is quite convenient to run the default pipeline.
Users only need to provide a region list in BED format which contains 3 columns (chr, start, end). It could be the peak calling result of sequencing data like ATAC-seq etc.
All required data and software will be installed automatically.
# Provide your region list in BED format with 3 columns.
foregroundBedPath <- system.file(package = "enrichTF", "extdata","testregion.bed")
# Call the whole pipeline
PECA_TF_enrich(inputForegroundBed = foregroundBedPath, genome = "testgenome") # change"testgenome" to one of "hg19", "hg38", "mm9", 'mm10' ! "testgenome" is only a test example.## Configure bsgenome ...##
## Attaching package: 'BiocGenerics'## The following objects are masked from 'package:stats':
##
## IQR, mad, sd, var, xtabs## The following objects are masked from 'package:base':
##
## Filter, Find, Map, Position, Reduce, anyDuplicated, aperm, append,
## as.data.frame, basename, cbind, colnames, dirname, do.call,
## duplicated, eval, evalq, get, grep, grepl, intersect, is.unsorted,
## lapply, mapply, match, mget, order, paste, pmax, pmax.int, pmin,
## pmin.int, rank, rbind, rownames, sapply, setdiff, sort, table,
## tapply, union, unique, unsplit, which.max, which.min##
## Attaching package: 'S4Vectors'## The following object is masked from 'package:utils':
##
## findMatches## The following objects are masked from 'package:base':
##
## I, expand.grid, unname##
## Attaching package: 'Biostrings'## The following object is masked from 'package:base':
##
## strsplit##
## Attaching package: 'rtracklayer'## The following object is masked from 'package:BiocIO':
##
## FileForFormat## Configure bsgenome finished## Configure motifpwm ...## Configure motifpwm finished## Configure motifPWMOBJ ...## Configure motifPWMOBJ finished## Configure RE_gene_corr ...## Configure RE_gene_corr finished## Configure Enhancer_RE_gene_corr ...## Configure Enhancer_RE_gene_corr finished## Configure MotifTFTable ...## Configure MotifTFTable finished## Configure MotifWeights ...## Configure MotifWeights finished## Configure TFgeneRelMtx ...## Configure TFgeneRelMtx finished## Configure SampleName ...## Configure SampleName finished## Configure OpenRegion ...## Configure OpenRegion finished## Configure ConserveRegion ...## Configure ConserveRegion finished## Configure HOMER ...## Configure HOMER finished## Configure OrgDb ...## Loading required package: AnnotationDbi## Loading required package: Biobase## Welcome to Bioconductor
##
## Vignettes contain introductory material; view with
## 'browseVignettes()'. To cite Bioconductor, see
## 'citation("Biobase")', and for packages 'citation("pkgname")'.## ## Configure OrgDb finished## >>>>>>==========================================================## Step Name:pipe_UnzipAndMergeBed## All Parameters for This Step:## |Input:## | bedInput:## | "/tmp/RtmpkNLw6Y/Rinst260ef6353e1bc2/enrichTF/extdata/testregion.bed"## |Output:## | bedOutput:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_00_pipe_UnzipAndMergeBed/testregion.bed"## |Other Parameters:## __________________________________________## Begin to check if it is finished...## 2023-10-24 17:07:04.156297## New step. Start processing data:## processing file:## source:/tmp/RtmpkNLw6Y/Rinst260ef6353e1bc2/enrichTF/extdata/testregion.bed## destination:/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_00_pipe_UnzipAndMergeBed/testregion.bed## 2023-10-24 17:07:04.168808## processing finished## [1] "2023-10-24 17:07:04 EDT"
## [1] "2023-10-24 17:07:04 EDT"
## [1] "2023-10-24 17:07:04 EDT"
## [1] "2023-10-24 17:07:04 EDT"
## [1] "2023-10-24 17:07:04 EDT"
## [1] "2023-10-24 17:07:04 EDT"## All results have been saved.## ==========================================================<<<<<<## ## >>>>>>==========================================================## Step Name:pipe_GenBackground## All Parameters for This Step:## |Input:## | inputForegroundBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_00_pipe_UnzipAndMergeBed/testregion.bed"## |Output:## | outputForegroundBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_01_pipe_GenBackground/testregion.foreground.bed"## | outputBackgroundBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_01_pipe_GenBackground/testregion.background.bed"## | outputRegionBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_01_pipe_GenBackground/testregion.allregion.bed"## |Other Parameters:## | bsgenome:## | An object of BSgenome## | regionLen:## | 1000## | sampleNumb:## | 0## __________________________________________## Begin to check if it is finished...## 2023-10-24 17:07:04.444789## New step. Start processing data:## 2023-10-24 17:07:04.886819## processing finished## [1] "2023-10-24 17:07:04 EDT"
## [1] "2023-10-24 17:07:04 EDT"
## [1] "2023-10-24 17:07:04 EDT"
## [1] "2023-10-24 17:07:04 EDT"
## [1] "2023-10-24 17:07:04 EDT"
## [1] "2023-10-24 17:07:04 EDT"## All results have been saved.## ==========================================================<<<<<<## ## >>>>>>==========================================================## Step Name:pipe_RegionConnectTargetGene## All Parameters for This Step:## |Input:## | inputForegroundBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_01_pipe_GenBackground/testregion.foreground.bed"## | inputBackgroundBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_01_pipe_GenBackground/testregion.background.bed"## |Output:## | outputForegroundBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_02_pipe_RegionConnectTargetGene/testregion.gene.foreground.bed"## | ouputForgroundGeneTxt:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_02_pipe_RegionConnectTargetGene/testregion.gene.foreground.txt"## | outputBackgroundBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_02_pipe_RegionConnectTargetGene/testregion.gene.background.bed"## | regularGeneCorrBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/refdir/testgenome/RE_gene_corr.bed"## | enhancerRegularGeneCorrBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/refdir/testgenome/Enhancer_RE_gene_corr.bed"## |Other Parameters:## __________________________________________## Begin to check if it is finished...## 2023-10-24 17:07:04.918338## New step. Start processing data:## 2023-10-24 17:07:05.797411## processing finished## [1] "2023-10-24 17:07:05 EDT"
## [1] "2023-10-24 17:07:05 EDT"
## [1] "2023-10-24 17:07:05 EDT"
## [1] "2023-10-24 17:07:05 EDT"
## [1] "2023-10-24 17:07:05 EDT"
## [1] "2023-10-24 17:07:05 EDT"## All results have been saved.## ==========================================================<<<<<<## ## >>>>>>==========================================================## Step Name:pipe_FindMotifsInRegions## All Parameters for This Step:## |Input:## | inputRegionBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_01_pipe_GenBackground/testregion.allregion.bed"## |Output:## | outputRegionMotifBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_03_pipe_FindMotifsInRegions/testregion.region.motif.bed"## | outputMotifBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_03_pipe_FindMotifsInRegions/testregion.motif.bed"## |Other Parameters:## | motifRc:## | "integrate"## | pwmObj:## | An object of PWMatrixList## | genome:## | "hg19"## | threads:## | 2## __________________________________________## Begin to check if it is finished...## 2023-10-24 17:07:05.836348## New step. Start processing data:## 2023-10-24 17:07:06.837759## processing finished## [1] "2023-10-24 17:07:06 EDT"
## [1] "2023-10-24 17:07:06 EDT"
## [1] "2023-10-24 17:07:06 EDT"
## [1] "2023-10-24 17:07:06 EDT"
## [1] "2023-10-24 17:07:06 EDT"
## [1] "2023-10-24 17:07:06 EDT"## All results have been saved.## ==========================================================<<<<<<## ## >>>>>>==========================================================## Step Name:pipe_TFsEnrichInRegions## All Parameters for This Step:## |Input:## | inputRegionBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_01_pipe_GenBackground/testregion.allregion.bed"## | inputRegionMotifBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_03_pipe_FindMotifsInRegions/testregion.region.motif.bed"## | inputForegroundGeneBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_02_pipe_RegionConnectTargetGene/testregion.gene.foreground.bed"## | inputBackgroundGeneBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_02_pipe_RegionConnectTargetGene/testregion.gene.background.bed"## | inputMotifWeights:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/refdir/testgenome/MotifWeights.RData"## | inputTFgeneRelMtx:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/refdir/testgenome/TFgeneRelMtx.RData"## | inputMotifTFTable:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/refdir/testgenome/MotifTFTable.RData"## |Output:## | outputTFsEnrichTxt:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_04_pipe_TFsEnrichInRegions/testregion.PECA_TF_enrich.txt"## |Other Parameters:## __________________________________________## Begin to check if it is finished...## 2023-10-24 17:07:06.874742## New step. Start processing data:## 2023-10-24 17:07:29.545844## processing finished## [1] "2023-10-24 17:07:29 EDT"
## [1] "2023-10-24 17:07:29 EDT"
## [1] "2023-10-24 17:07:29 EDT"
## [1] "2023-10-24 17:07:29 EDT"
## [1] "2023-10-24 17:07:29 EDT"
## [1] "2023-10-24 17:07:29 EDT"## All results have been saved.## ==========================================================<<<<<<## In this work, we define four kinds of relations which will be introduced in detail. And we will show how to test on enriched TFs rather than enriched motifs.
When users input a region list in BED format, the relations between regions and one gene can be obtained from two ways:
The regions which are overlapped with the functional region of this gene. This type of relation can be obtained from gene reference data.
The regions are enhancers or other distal regulatory elements of this gene. This type of relation can be derived from 3D genome sequencing experiments such as Hi-C and HiChIP.
As we have already integrated and organized the potential relations between regions and genes from both gene reference and experiment profiles, users can use this source to obtain region-gene relations for their input region list quite easily.
Moreover, users can also customize potential relations from their own work.
The correlation scores between genes and TFs come from our previous approach named PECA. These scores are scaled in the range [-1,1]. Users can also customize this relation matrix based on their own experiment data.
As we know, one TF may relate to multiple motifs. Here, we manually annotated more than 700 motifs with their corresponding TFs.
We scan all annotated motifs on each input region to check if there are some motifs located in this region. Then, the motif-region relation table can be generated.
Based on the foreground region list from users, this package can randomly sample the same number of regions from genome as background.
Foreground and background region lists can both be connected with genes in Region-gene Relation Table, so both of them will get a gene list. We name them as the foreground gene list and background gene list, respectively.
For each TF, we will perform t-test and obtain p-values for TF-Gene scores of foreground vs those of background. These scores come from Gene-TF Relation Matrix.
Each TF may be correlated with multiple motifs. For each motif, this package will do the following Fisher’s exact test to test if one motif is enriched in foreground regions.
| # with TF’s motif | # without TF’s motif | |
|---|---|---|
| # Foreground regions | ||
| # Background regions |
Finally, the package will select the motif with the most significant p-value to represent this TF and then compute p-values for all of the TFs.
Similar to previous section, the package will perform Fisher’s exact test with TF-Gene relation matrix. Different thresholds ([-1,1] with interval 0.1) for TF-Gene relation scores are used to filter the connected genes with foreground and background. In this way, we select the most significant p-value to represent this TF and calculate the p-values for all of TFs.
Finally, false discover rate correction is applied on the results of Fisher’s exact test with threshold. We rank all of the TFs based on their enrichment significance.
There are 4 steps to customize the workflow. You can follow the order below:
printMap()The function of genBackground is to generate random background regions based on the input regions and to set the length of input sequence regions. We try to select background sequence regions that match the GC-content distribution of the input sequence regions. The default length and number of background regions are 1000 and 10,000, respectively. The length of sequence regions used for motif finding is important. Here, we set 1000 as default value, which means for each region, we get sequences from -500 to +500 relative from center. For example,
setGenome("testgenome")## Configure bsgenome ...## Configure bsgenome finished## Configure motifpwm ...## Configure motifpwm finished## Configure motifPWMOBJ ...## Configure motifPWMOBJ finished## Configure RE_gene_corr ...## Configure RE_gene_corr finished## Configure Enhancer_RE_gene_corr ...## Configure Enhancer_RE_gene_corr finished## Configure MotifTFTable ...## Configure MotifTFTable finished## Configure MotifWeights ...## Configure MotifWeights finished## Configure TFgeneRelMtx ...## Configure TFgeneRelMtx finished## Configure SampleName ...## Configure SampleName finished## Configure OpenRegion ...## Configure OpenRegion finished## Configure ConserveRegion ...## Configure ConserveRegion finished## Configure HOMER ...## Configure HOMER finished## Configure OrgDb ...## Configure OrgDb finishedforegroundBedPath <- system.file(package = "enrichTF", "extdata","testregion.bed")
gen <- genBackground(inputForegroundBed = foregroundBedPath)## >>>>>>==========================================================## Step Name:pipe_GenBackground## All Parameters for This Step:## |Input:## | inputForegroundBed:## | "/tmp/RtmpkNLw6Y/Rinst260ef6353e1bc2/enrichTF/extdata/testregion.bed"## |Output:## | outputForegroundBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_01_pipe_GenBackground/testregion.foreground.bed"## | outputBackgroundBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_01_pipe_GenBackground/testregion.background.bed"## | outputRegionBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_01_pipe_GenBackground/testregion.allregion.bed"## |Other Parameters:## | bsgenome:## | An object of BSgenome## | regionLen:## | 1000## | sampleNumb:## | 0## __________________________________________## Begin to check if it is finished...## [1] "2023-10-24 17:07:29 EDT"
## [1] "2023-10-24 17:07:29 EDT"
## [1] "2023-10-24 17:07:29 EDT"
## [1] "2023-10-24 17:07:29 EDT"
## [1] "2023-10-24 17:07:29 EDT"
## [1] "2023-10-24 17:07:29 EDT"## 2023-10-24 17:07:29.817721## New step. Start processing data:## 2023-10-24 17:07:30.014096## processing finished## [1] "2023-10-24 17:07:30 EDT"
## [1] "2023-10-24 17:07:30 EDT"
## [1] "2023-10-24 17:07:30 EDT"
## [1] "2023-10-24 17:07:30 EDT"
## [1] "2023-10-24 17:07:30 EDT"
## [1] "2023-10-24 17:07:30 EDT"## All results have been saved.## ==========================================================<<<<<<## The function of regionConnectTargetGene is to connect foreground and background regions to their target genes, which are predicted from PECA model.
For example,
conTG <- enrichRegionConnectTargetGene(gen)## >>>>>>==========================================================## Step Name:pipe_RegionConnectTargetGene## All Parameters for This Step:## |Input:## | inputForegroundBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_01_pipe_GenBackground/testregion.foreground.bed"## | inputBackgroundBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_01_pipe_GenBackground/testregion.background.bed"## |Output:## | outputForegroundBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_02_pipe_RegionConnectTargetGene/testregion.gene.foreground.bed"## | ouputForgroundGeneTxt:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_02_pipe_RegionConnectTargetGene/testregion.gene.foreground.txt"## | outputBackgroundBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_02_pipe_RegionConnectTargetGene/testregion.gene.background.bed"## | regularGeneCorrBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/refdir/testgenome/RE_gene_corr.bed"## | enhancerRegularGeneCorrBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/refdir/testgenome/Enhancer_RE_gene_corr.bed"## |Other Parameters:## __________________________________________## Begin to check if it is finished...## [1] "2023-10-24 17:07:30 EDT"
## [1] "2023-10-24 17:07:30 EDT"
## [1] "2023-10-24 17:07:30 EDT"
## [1] "2023-10-24 17:07:30 EDT"
## [1] "2023-10-24 17:07:30 EDT"
## [1] "2023-10-24 17:07:30 EDT"## 2023-10-24 17:07:30.086425## New step. Start processing data:## 2023-10-24 17:07:30.813384## processing finished## [1] "2023-10-24 17:07:30 EDT"
## [1] "2023-10-24 17:07:30 EDT"
## [1] "2023-10-24 17:07:30 EDT"
## [1] "2023-10-24 17:07:30 EDT"
## [1] "2023-10-24 17:07:30 EDT"
## [1] "2023-10-24 17:07:30 EDT"## All results have been saved.## ==========================================================<<<<<<## The function MotifsInRegions is to scan for motif occurrences using the prepared PWMs and obtain the promising candidate motifs in these regions.
For example,
findMotif <- enrichFindMotifsInRegions(gen,motifRc="integrate")## >>>>>>==========================================================## Step Name:pipe_FindMotifsInRegions## All Parameters for This Step:## |Input:## | inputRegionBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_01_pipe_GenBackground/testregion.allregion.bed"## |Output:## | outputRegionMotifBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_03_pipe_FindMotifsInRegions/testregion.region.motif.bed"## | outputMotifBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_03_pipe_FindMotifsInRegions/testregion.motif.bed"## |Other Parameters:## | motifRc:## | "integrate"## | pwmObj:## | An object of PWMatrixList## | genome:## | "hg19"## | threads:## | 2## __________________________________________## Begin to check if it is finished...## [1] "2023-10-24 17:07:30 EDT"
## [1] "2023-10-24 17:07:30 EDT"
## [1] "2023-10-24 17:07:30 EDT"
## [1] "2023-10-24 17:07:30 EDT"
## [1] "2023-10-24 17:07:30 EDT"
## [1] "2023-10-24 17:07:30 EDT"## 2023-10-24 17:07:30.891245## New step. Start processing data:## 2023-10-24 17:07:31.904471## processing finished## [1] "2023-10-24 17:07:31 EDT"
## [1] "2023-10-24 17:07:31 EDT"
## [1] "2023-10-24 17:07:31 EDT"
## [1] "2023-10-24 17:07:31 EDT"
## [1] "2023-10-24 17:07:31 EDT"
## [1] "2023-10-24 17:07:31 EDT"## All results have been saved.## ==========================================================<<<<<<## The function of TFsEnrichInRegions is to test whether each TF is enriched on input regions. For example,
result <- enrichTFsEnrichInRegions(gen)## >>>>>>==========================================================## Step Name:pipe_TFsEnrichInRegions## All Parameters for This Step:## |Input:## | inputRegionBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_01_pipe_GenBackground/testregion.allregion.bed"## | inputRegionMotifBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_03_pipe_FindMotifsInRegions/testregion.region.motif.bed"## | inputForegroundGeneBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_02_pipe_RegionConnectTargetGene/testregion.gene.foreground.bed"## | inputBackgroundGeneBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_02_pipe_RegionConnectTargetGene/testregion.gene.background.bed"## | inputMotifWeights:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/refdir/testgenome/MotifWeights.RData"## | inputTFgeneRelMtx:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/refdir/testgenome/TFgeneRelMtx.RData"## | inputMotifTFTable:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/refdir/testgenome/MotifTFTable.RData"## |Output:## | outputTFsEnrichTxt:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_04_pipe_TFsEnrichInRegions/testregion.PECA_TF_enrich.txt"## |Other Parameters:## __________________________________________## Begin to check if it is finished...## [1] "2023-10-24 17:07:31 EDT"
## [1] "2023-10-24 17:07:31 EDT"
## [1] "2023-10-24 17:07:31 EDT"
## [1] "2023-10-24 17:07:31 EDT"
## [1] "2023-10-24 17:07:31 EDT"
## [1] "2023-10-24 17:07:31 EDT"## 2023-10-24 17:07:31.97156## New step. Start processing data:## 2023-10-24 17:07:32.013041## processing finished## [1] "2023-10-24 17:07:32 EDT"
## [1] "2023-10-24 17:07:32 EDT"
## [1] "2023-10-24 17:07:32 EDT"
## [1] "2023-10-24 17:07:32 EDT"
## [1] "2023-10-24 17:07:32 EDT"
## [1] "2023-10-24 17:07:32 EDT"## All results have been saved.## ==========================================================<<<<<<## Users can build a pipeline easily based on pipeFrame (pipeline framework)
library(magrittr)##
## Attaching package: 'magrittr'## The following object is masked from 'package:GenomicRanges':
##
## subtractsetGenome("testgenome")## Configure bsgenome ...## Configure bsgenome finished## Configure motifpwm ...## Configure motifpwm finished## Configure motifPWMOBJ ...## Configure motifPWMOBJ finished## Configure RE_gene_corr ...## Configure RE_gene_corr finished## Configure Enhancer_RE_gene_corr ...## Configure Enhancer_RE_gene_corr finished## Configure MotifTFTable ...## Configure MotifTFTable finished## Configure MotifWeights ...## Configure MotifWeights finished## Configure TFgeneRelMtx ...## Configure TFgeneRelMtx finished## Configure SampleName ...## Configure SampleName finished## Configure OpenRegion ...## Configure OpenRegion finished## Configure ConserveRegion ...## Configure ConserveRegion finished## Configure HOMER ...## Configure HOMER finished## Configure OrgDb ...## Configure OrgDb finishedforegroundBedPath <- system.file(package = "enrichTF", "extdata","testregion.bed")
result <- genBackground(inputForegroundBed = foregroundBedPath) %>%
enrichRegionConnectTargetGene%>%
enrichFindMotifsInRegions(motifRc="integrate") %>%
enrichTFsEnrichInRegions## >>>>>>==========================================================## Step Name:pipe_GenBackground## All Parameters for This Step:## |Input:## | inputForegroundBed:## | "/tmp/RtmpkNLw6Y/Rinst260ef6353e1bc2/enrichTF/extdata/testregion.bed"## |Output:## | outputForegroundBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_01_pipe_GenBackground/testregion.foreground.bed"## | outputBackgroundBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_01_pipe_GenBackground/testregion.background.bed"## | outputRegionBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_01_pipe_GenBackground/testregion.allregion.bed"## |Other Parameters:## | bsgenome:## | An object of BSgenome## | regionLen:## | 1000## | sampleNumb:## | 0## __________________________________________## Begin to check if it is finished...## [1] "2023-10-24 17:07:32 EDT"
## [1] "2023-10-24 17:07:32 EDT"
## [1] "2023-10-24 17:07:32 EDT"
## [1] "2023-10-24 17:07:32 EDT"
## [1] "2023-10-24 17:07:32 EDT"
## [1] "2023-10-24 17:07:32 EDT"## The step:`pipe_GenBackground` was finished. Nothing to do.## If you need to redo or rerun this step,please call 'clearStepCache(YourStepObject)'or remove file: /tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_01_pipe_GenBackground/pipeFrame.obj.f5c63819.rds## ==========================================================<<<<<<## ## >>>>>>==========================================================## Step Name:pipe_RegionConnectTargetGene## All Parameters for This Step:## |Input:## | inputForegroundBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_01_pipe_GenBackground/testregion.foreground.bed"## | inputBackgroundBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_01_pipe_GenBackground/testregion.background.bed"## |Output:## | outputForegroundBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_02_pipe_RegionConnectTargetGene/testregion.gene.foreground.bed"## | ouputForgroundGeneTxt:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_02_pipe_RegionConnectTargetGene/testregion.gene.foreground.txt"## | outputBackgroundBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_02_pipe_RegionConnectTargetGene/testregion.gene.background.bed"## | regularGeneCorrBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/refdir/testgenome/RE_gene_corr.bed"## | enhancerRegularGeneCorrBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/refdir/testgenome/Enhancer_RE_gene_corr.bed"## |Other Parameters:## __________________________________________## Begin to check if it is finished...## [1] "2023-10-24 17:07:32 EDT"
## [1] "2023-10-24 17:07:32 EDT"
## [1] "2023-10-24 17:07:32 EDT"
## [1] "2023-10-24 17:07:32 EDT"
## [1] "2023-10-24 17:07:32 EDT"
## [1] "2023-10-24 17:07:32 EDT"## The step:`pipe_RegionConnectTargetGene` was finished. Nothing to do.## If you need to redo or rerun this step,please call 'clearStepCache(YourStepObject)'or remove file: /tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_02_pipe_RegionConnectTargetGene/pipeFrame.obj.397a439e.rds## ==========================================================<<<<<<## ## >>>>>>==========================================================## Step Name:pipe_FindMotifsInRegions## All Parameters for This Step:## |Input:## | inputRegionBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_01_pipe_GenBackground/testregion.allregion.bed"## |Output:## | outputRegionMotifBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_03_pipe_FindMotifsInRegions/testregion.region.motif.bed"## | outputMotifBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_03_pipe_FindMotifsInRegions/testregion.motif.bed"## |Other Parameters:## | motifRc:## | "integrate"## | pwmObj:## | An object of PWMatrixList## | genome:## | "hg19"## | threads:## | 2## __________________________________________## Begin to check if it is finished...## [1] "2023-10-24 17:07:32 EDT"
## [1] "2023-10-24 17:07:32 EDT"
## [1] "2023-10-24 17:07:32 EDT"
## [1] "2023-10-24 17:07:32 EDT"
## [1] "2023-10-24 17:07:32 EDT"
## [1] "2023-10-24 17:07:32 EDT"## The step:`pipe_FindMotifsInRegions` was finished. Nothing to do.## If you need to redo or rerun this step,please call 'clearStepCache(YourStepObject)'or remove file: /tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_03_pipe_FindMotifsInRegions/pipeFrame.obj.d876654a.rds## ==========================================================<<<<<<## ## >>>>>>==========================================================## Step Name:pipe_TFsEnrichInRegions## All Parameters for This Step:## |Input:## | inputRegionBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_01_pipe_GenBackground/testregion.allregion.bed"## | inputRegionMotifBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_03_pipe_FindMotifsInRegions/testregion.region.motif.bed"## | inputForegroundGeneBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_02_pipe_RegionConnectTargetGene/testregion.gene.foreground.bed"## | inputBackgroundGeneBed:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_02_pipe_RegionConnectTargetGene/testregion.gene.background.bed"## | inputMotifWeights:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/refdir/testgenome/MotifWeights.RData"## | inputTFgeneRelMtx:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/refdir/testgenome/TFgeneRelMtx.RData"## | inputMotifTFTable:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/refdir/testgenome/MotifTFTable.RData"## |Output:## | outputTFsEnrichTxt:## | "/tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_04_pipe_TFsEnrichInRegions/testregion.PECA_TF_enrich.txt"## |Other Parameters:## __________________________________________## Begin to check if it is finished...## [1] "2023-10-24 17:07:32 EDT"
## [1] "2023-10-24 17:07:32 EDT"
## [1] "2023-10-24 17:07:32 EDT"
## [1] "2023-10-24 17:07:32 EDT"
## [1] "2023-10-24 17:07:32 EDT"
## [1] "2023-10-24 17:07:32 EDT"## The step:`pipe_TFsEnrichInRegions` was finished. Nothing to do.## If you need to redo or rerun this step,please call 'clearStepCache(YourStepObject)'or remove file: /tmp/RtmpkNLw6Y/Rbuild260ef617cb760a/enrichTF/vignettes/enrichTF-pipeline/Step_04_pipe_TFsEnrichInRegions/pipeFrame.obj.a3509dd9.rds## ==========================================================<<<<<<## Here we show an example result when applying this package on our own data.
examplefile <- system.file(package = "enrichTF", "extdata","result.example.txt")
read.table(examplefile, sep='\t', header = TRUE)%>%
knitr::kable() | TF | Motif_enrichment | Targt_gene_enrichment | P_value | FDR |
|---|---|---|---|---|
| MITF | 0.0000000 | 0.0383390 | 0e+00 | 0.00e+00 |
| TFEC | 0.0000000 | 0.4857386 | 0e+00 | 0.00e+00 |
| MLX | 0.0000000 | 0.8147766 | 0e+00 | 0.00e+00 |
| TFEB | 0.0000000 | 0.5303603 | 0e+00 | 0.00e+00 |
| TFE3 | 0.0000000 | 0.4067377 | 0e+00 | 0.00e+00 |
| MZF1 | 0.0000000 | 0.0532782 | 0e+00 | 5.00e-07 |
| USF2 | 0.0000000 | 0.5766573 | 0e+00 | 9.00e-07 |
| USF1 | 0.0000000 | 0.2185547 | 0e+00 | 9.00e-07 |
| CEBPA | 0.0000012 | 0.0046147 | 0e+00 | 1.50e-06 |
| PKNOX2 | 0.0000001 | 0.0140376 | 0e+00 | 2.30e-06 |
| FOXA3 | 0.0000071 | 0.8939957 | 0e+00 | 2.30e-06 |
| CEBPB | 0.0000012 | 0.0118270 | 0e+00 | 2.30e-06 |
| PKNOX1 | 0.0000001 | 0.9561768 | 0e+00 | 2.30e-06 |
| MEIS1 | 0.0000001 | 0.1229445 | 1e-07 | 2.40e-06 |
| TGIF1 | 0.0000001 | 0.2912322 | 1e-07 | 2.40e-06 |
| TGIF2 | 0.0000001 | 0.9898682 | 1e-07 | 2.40e-06 |
| FOXD3 | 0.0000071 | 0.8125302 | 1e-07 | 4.10e-06 |
| MEIS2 | 0.0000001 | 0.8463790 | 1e-07 | 5.00e-06 |
| DLX3 | 0.0000477 | 0.0043755 | 2e-07 | 5.90e-06 |
| MYCN | 0.0000004 | 0.4520504 | 2e-07 | 7.30e-06 |
| FOXG1 | 0.0000071 | 0.8432650 | 2e-07 | 7.40e-06 |
| ID1 | 0.0000017 | 0.0093116 | 3e-07 | 9.60e-06 |
| MAX | 0.0000017 | 0.5516440 | 3e-07 | 1.04e-05 |
| FOXC2 | 0.0000071 | 0.8538362 | 4e-07 | 1.09e-05 |
| DLX5 | 0.0000477 | 0.0017341 | 4e-07 | 1.09e-05 |
| MYC | 0.0000004 | 0.2194055 | 4e-07 | 1.11e-05 |
| MEIS3 | 0.0000001 | 0.8591956 | 4e-07 | 1.11e-05 |
| RORA | 0.0001585 | 0.0021495 | 4e-07 | 1.11e-05 |
| CREB3L2 | 0.0000017 | 0.6550453 | 5e-07 | 1.15e-05 |
sessionInfo()## R version 4.3.1 (2023-06-16)
## Platform: x86_64-pc-linux-gnu (64-bit)
## Running under: Ubuntu 22.04.3 LTS
##
## Matrix products: default
## BLAS: /home/biocbuild/bbs-3.18-bioc/R/lib/libRblas.so
## LAPACK: /usr/lib/x86_64-linux-gnu/lapack/liblapack.so.3.10.0
##
## locale:
## [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
## [3] LC_TIME=en_GB LC_COLLATE=C
## [5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_US.UTF-8
## [7] LC_PAPER=en_US.UTF-8 LC_NAME=C
## [9] LC_ADDRESS=C LC_TELEPHONE=C
## [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
##
## time zone: America/New_York
## tzcode source: system (glibc)
##
## attached base packages:
## [1] stats4 stats graphics grDevices utils datasets methods
## [8] base
##
## other attached packages:
## [1] magrittr_2.0.3 org.Hs.eg.db_3.18.0
## [3] AnnotationDbi_1.64.0 Biobase_2.62.0
## [5] BSgenome.Hsapiens.UCSC.hg19_1.4.3 BSgenome_1.70.0
## [7] rtracklayer_1.62.0 BiocIO_1.12.0
## [9] Biostrings_2.70.0 XVector_0.42.0
## [11] GenomicRanges_1.54.0 GenomeInfoDb_1.38.0
## [13] IRanges_2.36.0 S4Vectors_0.40.0
## [15] BiocGenerics_0.48.0 enrichTF_1.18.0
## [17] pipeFrame_1.18.0
##
## loaded via a namespace (and not attached):
## [1] splines_4.3.1 later_1.3.1
## [3] bitops_1.0-7 ggplotify_0.1.2
## [5] filelock_1.0.2 tibble_3.2.1
## [7] R.oo_1.25.0 polyclip_1.10-6
## [9] XML_3.99-0.14 DirichletMultinomial_1.44.0
## [11] lifecycle_1.0.3 rstatix_0.7.2
## [13] fastcluster_1.2.3 lattice_0.22-5
## [15] MASS_7.3-60 heatmap3_1.1.9
## [17] backports_1.4.1 sass_0.4.7
## [19] rmarkdown_2.25 jquerylib_0.1.4
## [21] yaml_2.3.7 httpuv_1.6.12
## [23] cowplot_1.1.1 DBI_1.1.3
## [25] CNEr_1.38.0 RColorBrewer_1.1-3
## [27] abind_1.4-5 zlibbioc_1.48.0
## [29] purrr_1.0.2 R.utils_2.12.2
## [31] ggraph_2.1.0 RCurl_1.98-1.12
## [33] yulab.utils_0.1.0 pracma_2.4.2
## [35] tweenr_2.0.2 rappdirs_0.3.3
## [37] GenomeInfoDbData_1.2.11 enrichplot_1.22.0
## [39] ggrepel_0.9.4 tidytree_0.4.5
## [41] seqLogo_1.68.0 annotate_1.80.0
## [43] codetools_0.2-19 DelayedArray_0.28.0
## [45] DOSE_3.28.0 ggforce_0.4.1
## [47] tidyselect_1.2.0 aplot_0.2.2
## [49] farver_2.1.1 viridis_0.6.4
## [51] matrixStats_1.0.0 BiocFileCache_2.10.0
## [53] GenomicAlignments_1.38.0 jsonlite_1.8.7
## [55] ellipsis_0.3.2 tidygraph_1.2.3
## [57] motifmatchr_1.24.0 tools_4.3.1
## [59] treeio_1.26.0 HPO.db_0.99.2
## [61] TFMPvalue_0.0.9 Rcpp_1.0.11
## [63] glue_1.6.2 gridExtra_2.3
## [65] SparseArray_1.2.0 xfun_0.40
## [67] qvalue_2.34.0 MatrixGenerics_1.14.0
## [69] dplyr_1.1.3 withr_2.5.1
## [71] BiocManager_1.30.22 fastmap_1.1.1
## [73] fansi_1.0.5 caTools_1.18.2
## [75] digest_0.6.33 R6_2.5.1
## [77] mime_0.12 gridGraphics_0.5-1
## [79] colorspace_2.1-0 GO.db_3.18.0
## [81] gtools_3.9.4 poweRlaw_0.70.6
## [83] RSQLite_2.3.1 R.methodsS3_1.8.2
## [85] utf8_1.2.4 tidyr_1.3.0
## [87] generics_0.1.3 data.table_1.14.8
## [89] graphlayouts_1.0.1 httr_1.4.7
## [91] htmlwidgets_1.6.2 S4Arrays_1.2.0
## [93] scatterpie_0.2.1 TFBSTools_1.40.0
## [95] pkgconfig_2.0.3 gtable_0.3.4
## [97] blob_1.2.4 shadowtext_0.1.2
## [99] clusterProfiler_4.10.0 htmltools_0.5.6.1
## [101] carData_3.0-5 fgsea_1.28.0
## [103] scales_1.2.1 png_0.1-8
## [105] ggfun_0.1.3 knitr_1.44
## [107] tzdb_0.4.0 reshape2_1.4.4
## [109] rjson_0.2.21 nlme_3.1-163
## [111] visNetwork_2.1.2 curl_5.1.0
## [113] JASPAR2018_1.1.1 cachem_1.0.8
## [115] stringr_1.5.0 BiocVersion_3.18.0
## [117] parallel_4.3.1 HDO.db_0.99.1
## [119] restfulr_0.0.15 pillar_1.9.0
## [121] grid_4.3.1 vctrs_0.6.4
## [123] ggpubr_0.6.0 promises_1.2.1
## [125] car_3.1-2 dbplyr_2.3.4
## [127] xtable_1.8-4 evaluate_0.22
## [129] readr_2.1.4 cli_3.6.1
## [131] compiler_4.3.1 Rsamtools_2.18.0
## [133] rlang_1.1.1 crayon_1.5.2
## [135] ggsignif_0.6.4 fs_1.6.3
## [137] plyr_1.8.9 stringi_1.7.12
## [139] viridisLite_0.4.2 BiocParallel_1.36.0
## [141] MPO.db_0.99.7 munsell_0.5.0
## [143] lazyeval_0.2.2 GOSemSim_2.28.0
## [145] Matrix_1.6-1.1 patchwork_1.1.3
## [147] hms_1.1.3 bit64_4.0.5
## [149] ggplot2_3.4.4 KEGGREST_1.42.0
## [151] shiny_1.7.5.1 SummarizedExperiment_1.32.0
## [153] interactiveDisplayBase_1.40.0 AnnotationHub_3.10.0
## [155] broom_1.0.5 igraph_1.5.1
## [157] memoise_2.0.1 bslib_0.5.1
## [159] ggtree_3.10.0 fastmatch_1.1-4
## [161] bit_4.0.5 gson_0.1.0
## [163] ape_5.7-1