Using target to predict combined binding

Mahmoud Ahmed

2023-04-25

Overview

In this document, we extend the BETA algorithm to predict cooperative and/or competitive binding events of two factors at the same region of interest. This method uses the same theory described in vignette('target') with one modification. Instead of using an individual statistics from a factor perturbation experiment, we use two separate statistics from two comparable factor perturbation experiments to define a regulatory interaction (\(RI\)). This regulatory interaction terms are used to calculate rank products of the regions of interest. In addition, we illustrate the use of this method with a simulated dataset and an example of real-world data.

Theory

This method is an extension of the BETA basic algorithm. Briefly, the rank product of a regions of interest (\(RP_g\)) is the product of two terms; the rank of a statistics from factor perturbation experiment (\(R_{ge}\)) and the rank of the regulatory potential (\(R_{gb}\)). The regulatory potential of region of interest (\(S_g\)) is the sum of the peaks within a predefined distance:

\[ S_g = \sum_{i=1}^k e^{-(0.5+4\Delta_i)} \quad\text{and}\quad RP_g = \frac{R_{gb}\times R_{ge}}{n^2} \] Where \(p\) is a peak in \(\{1,...,k\}\) with a distance, \(\Delta\), from the center of the region of interest.

Two determine the relation of two factors \(x\) and \(y\) on a common peak near a region of interest, we define a new term; the regulatory interaction (\(RI\)) as the product of two signed statistics from comparable perturbation experiments. The rank of this term is used to calculate a rank product (\(PR_g\)) for each region of interest as described above

\[ RI_{g} = x_{ge}\times y_{ge} \quad\text{and}\quad RP_g = \frac{R_{gb}\times RI_{ge}}{n^2} \]

This term would represent the interaction magnitude assuming a linear relation between the two factor. The sign of the term would define the direction of the relation were positive means cooperative and negative means competitive.

Example

The target package contain two simulated datasets. sim_peaks is random peaks with random distances from the transcripts of chromosome 1 of the mm10 mouse genome. sim_transcripts is the same transcripts with random singed statistics assigned to each. In the following two examples, we introduce changes in these statistics to simulate conditions where two factors are working cooperatively or competitively on the same transcripts.

# load libraries
library(target)
# load data
data("sim_peaks")
data("sim_transcripts")

To help visualize these cases, a plotting function plot_profiles was constructed to introduce the changes change in the statistics of the transcripts near the n number of peaks. The source code for the function is available in inst/extdata/plot-profiles.R which we source to use here. The output of the function is a series of plots to visualize the statistics of the two factors before and after introducing the changes, the peaks distances and scores and the predicted functions of the factors individually and combined.

# source the plotting function
source(system.file('extdata', 'plot-profiles.R', package = 'target'))

Cooperative factors example

The first two inputs to the plotting function is the simulated peaks and transcripts. We chose to introduce positive changes to the statistics of the transcripts with the top 5000 nearby peaks of the two factors.

# simulate and plot cooperative factors
plot_profiles(sim_peaks,
              sim_transcripts,
              n = 5000,
              change = c(3, 3))

The changes introduced above are illustrated in the right upper quadrant of the scatter plot. The predicted functions of the two factors are similar, as shown by distribution function of the regulatory potential of their targets. Finally, when the targets are predicted based on the two statistics combined, the sign of the statistics product determines the direction of the factor interactions. Here, more higher ranking transcripts had positive/red/ cooperative change associated with the two factors.

Competitive factors example

Similar to the example above, we chose to introduce change to the statistics of the transcripts with nearby peaks for the two factors. To simulate the competitive conditions of the two factors, the changes have opposite signs.

# simulate and plot competitve factors
plot_profiles(sim_peaks,
              sim_transcripts,
              n = 5000,
              change = c(3, -3))

By contrast, the changes associated with each factor are opposite in sign, as shown in the lower right quadrant. The predicted functions of the individual factors are also the opposite. The predicted combined effect of the two factors is negative/green/competitive.

References

Wang S, Sun H, Ma J, et al. Target analysis by integration of transcriptome and ChIP-seq data with BETA. Nat Protoc. 2013;8(12):2502–2515. doi:10.1038/nprot.2013.150

sessionInfo()
#> R version 4.3.0 RC (2023-04-13 r84269)
#> Platform: x86_64-pc-linux-gnu (64-bit)
#> Running under: Ubuntu 22.04.2 LTS
#> 
#> Matrix products: default
#> BLAS:   /home/biocbuild/bbs-3.17-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] stats     graphics  grDevices utils     datasets  methods   base     
#> 
#> other attached packages:
#> [1] target_1.14.0
#> 
#> loaded via a namespace (and not attached):
#>  [1] cli_3.6.1               knitr_1.42              rlang_1.1.0            
#>  [4] xfun_0.39               highr_0.10              promises_1.2.0.1       
#>  [7] shiny_1.7.4             xtable_1.8-4            jsonlite_1.8.4         
#> [10] S4Vectors_0.38.0        RCurl_1.98-1.12         htmltools_0.5.5        
#> [13] httpuv_1.6.9            sass_0.4.5              stats4_4.3.0           
#> [16] rmarkdown_2.21          evaluate_0.20           jquerylib_0.1.4        
#> [19] ellipsis_0.3.2          bitops_1.0-7            fastmap_1.1.1          
#> [22] lifecycle_1.0.3         yaml_2.3.7              IRanges_2.34.0         
#> [25] GenomeInfoDb_1.36.0     compiler_4.3.0          Rcpp_1.0.10            
#> [28] XVector_0.40.0          later_1.3.0             digest_0.6.31          
#> [31] R6_2.5.1                GenomeInfoDbData_1.2.10 magrittr_2.0.3         
#> [34] GenomicRanges_1.52.0    bslib_0.4.2             tools_4.3.0            
#> [37] mime_0.12               matrixStats_0.63.0      zlibbioc_1.46.0        
#> [40] BiocGenerics_0.46.0     cachem_1.0.7