Abstract:A computer program for calculating clusters of binding sites of various transcription factors (TFs) according to the genomic coordinates of the ChIP-seq (Chromatin ImmunoPrecipitation-sequencing) profile peaks is developed. The statistical features of the distribution of the transcription factors’ binding sites (TFBSs) in the mouse genome, obtained with the help of ChIP-seq experiments in embryonic stem cells, are considered. Clusters of sites containing at least four binding sites of various TFs in the mouse genome are determined and their localization relative to the regulatory regions of the genes is described. Two types of colocalization of the sites are confirmed: clusters containing binding sites of factors Oct4, Nanog, and Sox2 located in the distal regions and clusters with n-Myc and c-Myc binding sites located mainly in the promoter regions of mouse genes. Analysis of the new ChIP-seq data on the binding of TFs Nr5a2, Tbx3, Cep, SRF, and USF1 in the same cell type confirmed the differentiation of clusters of the TFBSs into two types: those containing pluripotency regulator binding sites (Oct4, Nanog, and Sox2) and those not containing them. A computer program for the statistical processing of the data on the location of the sites in the genes is developed; it uses the experimental data on site localization obtained by ChIP-seq methods in mouse and human genomes. With the help of this program, the localization patterns of the binding sites of various TFs are detected. The distances between the closest binding sites of the TF groups Oct4, Nanog, and Sox2 and the binding sites of other factors in site clusters that serve as a basis for the analysis of the joint binding of protein complexes to DNA are calculated. The fraction of the presence of the known nucleotide motifs of TFBSs in the genomic regions of ChIP-seq is calculated. The weight matrices for such nucleotide motifs are recalculated. The correlation between the presence of motifs and the ChIP-seq binding intensity is shown. The programs implementing the computerized methods for assessing the clustering of binding sites of various TFs for new ChIP-seq data are available upon request from the authors.

Identifying ChIP-seq enrichment using MACS

Using MACS to Identify Peaks from ChIP‐Seq Data

Model-based Analysis of ChIP-Seq (MACS)

Diffreps: Detecting Differential Chromatin Modification Sites from ChIP-seq Data with Biological Replicates.

Computer Analysis of Co-Localization of Transcription Factor Binding Sites in Genome by ChIP-seq Data

An Integrated Software System for Analyzing ChIP-chip and ChIP-seq Data

QChIPat: a quantitative method to identify distinct binding patterns for two biological ChIP-seq samples in different experimental conditions

RACS: rapid analysis of ChIP-Seq data for contig based genomes

Computer Analysis of Colocalization of the TFs’ Binding Sites in the Genome According to the ChIP-seq Data

An Integrated Pipeline For The Genome-Wide Analysis Of Transcription Factor Binding Sites From Chip-Seq

Qb-13006-xsl 54..70

MAPS: Model-Based Analysis of Long-Range Chromatin Interactions from PLAC-seq and HiChIP Experiments

PICS: Probabilistic Inference for ChIP-seq

Using CisGenome to Analyze ChIP‐chip and ChIP‐seq Data

Model-based analysis of chromatin interactions from dCas9-Based CAPTURE-3C-seq

Analyzing ChIP-seq Data Based on Multiple Knowledge Sources for Histone Modification.

Unified Analysis of Multiple ChIP-Seq Datasets

CistromeFinder for ChIP-seq and DNase-seq Data Reuse

scAMACE: Model-based approach to the joint analysis of single-cell data on chromatin accessibility, gene expression and methylation

ChIP-chip versus ChIP-seq: Lessons for experimental design and data analysis

Visual annotations and a supervised learning approach for evaluating and calibrating ChIP-seq peak detectors