- Open Access
Predicting transcription factor site occupancy using DNA sequence intrinsic and cell-type specific chromatin features
© Kumar and Bucher. 2016
- Published: 11 January 2016
Understanding the mechanisms by which transcription factors (TF) are recruited to their physiological target sites is crucial for understanding gene regulation. DNA sequence intrinsic features such as predicted binding affinity are often not very effective in predicting in vivo site occupancy and in any case could not explain cell-type specific binding events. Recent reports show that chromatin accessibility, nucleosome occupancy and specific histone post-translational modifications greatly influence TF site occupancy in vivo. In this work, we use machine-learning methods to build predictive models and assess the relative importance of different sequence-intrinsic and chromatin features in the TF-to-target-site recruitment process.
Our study primarily relies on recent data published by the ENCODE consortium. Five dissimilar TFs assayed in multiple cell-types were selected as examples: CTCF, JunD, REST, GABP and USF2. We used two types of candidate target sites: (a) predicted sites obtained by scanning the whole genome with a position weight matrix, and (b) cell-type specific peak lists provided by ENCODE. Quantitative in vivo occupancy levels in different cell-types were based on ChIP-seq data for the corresponding TFs. In parallel, we computed a number of associated sequence-intrinsic and experimental features (histone modification, DNase I hypersensitivity, etc.) for each site. Machine learning algorithms were then used in a binary classification and regression framework to predict site occupancy and binding strength, for the purpose of assessing the relative importance of different contextual features.
We observed striking differences in the feature importance rankings between the five factors tested. PWM-scores were amongst the most important features only for CTCF and REST but of little value for JunD and USF2. Chromatin accessibility and active histone marks are potent predictors for all factors except REST. Structural DNA parameters, repressive and gene body associated histone marks are generally of little or no predictive value.
We define a general and extensible computational framework for analyzing the importance of various DNA-intrinsic and chromatin-associated features in determining cell-type specific TF binding to target sites. The application of our methodology to ENCODE data has led to new insights on transcription regulatory processes and may serve as example for future studies encompassing even larger datasets.
- Support Vector Machine
- Histone Mark
- Radial Basis Function Kernel
- Nucleosome Occupancy
- Position Weight Matrix
Genes are regulated by transcription factors (TF) binding to physiological target sites in the genome. TFs may bind directly to target sites through sequence-specific protein-DNA interactions, or indirectly through protein-protein interactions with other TFs . Understanding the mechanisms by which TFs are recruited to their target sites is essential for the understanding of gene regulation. For a long time, research in this area has been hampered by the lack of powerful assays to study TF binding events in vivo. This has drastically changed with the advent of the ChIP-seq technology which allows for comprehensive, genome-wide mapping of all in vivo bound sites of a given TF in a particular cell type at near base-pair resolution . What has become clear from ChIP-seq experiments is that the intrinsic binding specificity of a TF can only partly explain the in vivo site occupancy patterns, which in addition were found to be tissue-specific . The recruitment of TFs to target sites thus depends on both DNA-intrinsic properties and cell type specific covariates.
The intrinsic DNA binding specificity of a TF is commonly represented by a so-called position weight matrix (PWM) . A PWM is a N × 4 matrix whose elements define the TF’s binding preferences for the four bases of the DNA alphabet along a binding site of length N. Base preferences may either be expressed as occurrence probabilities or as (additive) binding energies. PWMs are basic binding site models with known limitations. For instance, they cannot model nearest neighbor dependencies nor can they account for variable spacing between reverse-complementary half-sites of homodimeric TFs . Nevertheless, it is generally agreed that at least some PWMs are good predictors of in vitro binding affinity of the corresponding TFs. Moreover, large collections of PWMs are available from public databases such as JASPAR . More advanced modeling techniques have been proposed for describing more accurately the binding specificity of a TF  but corresponding factor-specific models are not yet available for more than a handful of TFs.
Other DNA sequence-intrinsic contextual features have been used to reduce false positive rates in PWM-based in vivo TF binding site (TFBS) prediction, for instance DNA structural properties , . Double-stranded DNA possesses anisotropic flexibility, which determines its stability and rigidity, properties that potentially interfere with DNA-protein binding. These structural properties, which are broadly classified into (i) DNA conformation (A-DNA philicity and Z-DNA stability energy), (ii) flexibility (B-DNA twist, protein DNA twist, propeller twist and bending stiffness) and (iii) stability (duplex disruption and stability free energy, stacking energy and denaturation) are sequence dependent, and at least partly predictable from structural characteristics of dinucleotides as revealed by crystal structures of double-stranded oligonucleotides. Cross-species conservation is another sequence-derived feature that has been successfully exploited for distinguishing biologically functional TFBS (evolving under purifying selection) from non-functional ones .
Several recent studies have reported that TF binding is influenced (and thus potentially predictable) by chromatin contextual features such as DNA accessibility, nucleosome occupancy, or the presence of specific histone post-translational modifications . This is expected as these modifications were known to be associated with gene regulation processes long before. However, with the availability of high resolution ChIP-seq data for many histone marks in many cell types (e.g. ENCODE collection) it has become feasible to use these features in a systematic manner to build predictive models of in vivo TF binding. Proof of concept comes from a number of recent studies where chromatin features have successfully been exploited to predict TF site occupancy with high confidence –.
In the current work, we have combined multiple sequence-intrinsic parameters with experimental chromatin feature types and used machine-learning methods to assess their relative importance in the TF-to-target-site recruitment process, in the hope to gain new insights into transcription regulatory mechanisms. We used five diverse TFs as examples: (i) the insulator protein CTCF featuring 11 zinc finger domains, which has also been attributed diverse functions including transcriptional repression, genomic imprinting and tumor suppressor , (ii) JunD, a leucine zipper protein and member of the activator protein 1 (AP1) family, (iii) the transcriptional activator GABPA, the only obligatory multimeric TF within the Ets family , (iv) the transcriptional repressor REST, also known as neuron-restrictive silencer factor (NRSF), and (v) USF2 a member of the evolutionary conserved basic helix-loop-helix leucine zipper TF family. These factors were chosen on the one hand because they represent diverse classes of transcription factors, and on the other hand because they were extensively assayed by ChIP-Seq and other chromatin profiling techniques (e.g. DNase I hypersensitivity assays) in multiple cell types by the ENCODE consortium.
Dataset definition and download
Binding sites were represented by a single genomic position corresponding to the center of the binding regions. Two different types of TFBS collections were used in this study.
Genome-wide predicted sites: PWMs for CTCF (MA0139.1), JunD (MA0492.1), GABPA (MA0062.2), REST (MA0138.2) and USF2 (MA0526.1) were downloaded from the JASPAR database . The human genome assembly hg19 was scanned with these PWMs using in-house tool PWMScan (http://ccg.vital-it.ch/pwmtools) with a P-value threshold of 0.0001. Predicted sites with no or only marginal ChIP-seq tag counts (<3) in any of the cell lines analyzed were removed. Final numbers of sites for each factors were: CTCF (77878), JunD (65526), GABP (51926), REST (48171) and USF2 (37138). The midpoint of the motif was used as the reference position in the peak lists. (For PWMs of even-numbered length, the position immediately upstream of the midpoint was used).
ENCODE peak lists were downloaded from GEO (see Additional file 1 for sample ids). The “narrow peak” format version was used wherever available else the “broad peak” format was used. In either case, we used the midpoint of the peak region as reference point.
TFBS lists were annotated with cell-type specific and sequence-derived quantitative (numerical) features. These annotations were stored in a table with rows corresponding to binding sites and columns to annotation features. Cell-type specific experimental features were calculated as follows:
ChIP-seq and DNase-seq tag coverage: Files containing mapped sequence tags from a single experiment were taken from the Mass Genome Annotation (MGA) repository of the Eukaryotic Promoter Database EPD . Most of these datasets were originally downloaded in BAM format from the UCSC genome browser database (http://genome.ucsc.edu/ENCODE/downloads.html). Others were generated locally by mapping the sequence reads downloaded from the short read archive (SRA)  to the human genome assembly hg19 with the aid of Bowtie . For determining the tag counts, we used only the 5’ end positions of the mapped sequence reads. For non-histone ChIP-seq targets (TFs, cofactors, and PolII) these positions were shifted downstream or upstream towards the estimated fragment center by an appropriate distance (determined using cross correlation between + strand tags and – strand tags). The GEO sample accession numbers and specific shifting distances are given in Additional file 1. Multiple tags mapping exactly the same genome positions were counted only once per experiment. For transcription factors and corresponding cofactors (Rad21 for CTCF and FOSL for JunD) as well as for PolII, tags were counted in a window of ±100 bases relative to the TFBS center. A window of ±250 bases was used for DNase I hypersensitivity data. Histone modification tags were counted in windows of ±500 bases.
Shape-based evaluation of DNase-seq (DNase I hypersensitivity data, also referred to as “digital genomic foot-printing” (DGF)) profiles. We essentially followed the probabilistic partitioning protocol used by Nair et al.  for ranking and refocusing CTCF sites. This method attempts to classify the input data into typical and atypical examples. Limited shifting of individual genomic regions is allowed to optimally match the aggregate profile of typical examples. Here, DNase-seq tags were extracted from a ±500 bp region around TFBS center positions and binned in 10 bp windows. Shifting was limited to ±50 bp (11 shift states). The method returns a probability p that a given genomic region constitutes a typical example plus an optimal shifting distance s. The output variable p was used as quantitative shape-based DGF feature in our analyses.
Sequence-derived cell type-independent features:
Position weight matrix (PWM) score: For predicted TFBSs, the score was obtained as a by-product of scanning the whole-genome with the PWM. For ENCODE peaks, the PWM score of the best match to the JASPAR matrix within ±100 bp form the peak center was used. For computing the PWM scores, the base probability matrices from JASPAR were converted into log-odds matrices using the following base composition of the human genome as background frequencies: A - 29 %, C - 21 %, G - 21 %, T - 29 %.
Base (oligonucleotide) compositions: mono-, di-, tri, tetra- and penta-nucleotide frequencies were calculated for ±100 base sequences around TFBS sites. The sequences were split into upstream and downstream parts relative to the TFBS center. Oligonucleotide frequencies were then determined separately for the two parts and for the complete sequence region, and the frequencies obtained in this way were used as three separate annotation features.
Nucleosome occupancy prediction: average nucleosome occupancy per base was calculated for sequences around TFBS using the model built by Kaplan et al. (Version 3.0) . Sequences ±1000 bases long were used in order to avoid border effects. The final feature consists of the average score of ±10 bases around the TFBS center.
Structural features: These were calculated with the R script “dinucl.R“ from the DnaFVP package (http://dnafvp.sourceforge.net/) . In essence, this script assigns a score taken from a table to each dinucleotide within a DNA sequence. The average value over the ±10 bp regions (20 dinculeotide positions) relative to the TFBS center was used as quantitative structural feature. The following 10 features were used: A-DNA philicity, Z-DNA stability energy, duplex disruption free energy, duplex stability free energy, stacking energy, DNA denaturation, B-DNA twist, protein DNA twist, DNA propeller twist, and bending stiffness.
Conservation and polymorphism: A compressed version of the phastCons46way track of about 10 bp resolution was used as input , . The corresponding track file (available at ftp://ccg.vital-it.ch/mga/hg19/phastcons/) presents conservation profiles in a “counts per position” format similar to the ChIP-seq data files. The average PHASTCONS score within ±50 bases was used as conservation feature. SNP frequencies are based on dbSNP132 (~30 million SNPs) . All variants of type SNP (excluding indels) were intersected with the genomic coordinates of extended TFBS regions (±100 bp) using Annovar . The average SNP count within ±10 bases was used as final SNP feature.
Feature selection and model building
Data matrices consisting of features in columns and TF sites/peaks in rows were first normalized by centering columns on the mean and then scaled by dividing the centered values by the standard deviations of the corresponding columns.
Classifiers to discriminate strong from weak TF binding sites (top and bottom 20 %) were built and evaluated with the R functions train and trainControl from the Caret package. First, a given data table was randomly split into 90 % sites for model building and 10 % for testing. Support Vector Machines (SVM) and Random Forest (RF) were used to build a model under leave-one-group-out cross-validation (10 groups). Polynomial and radial basis function kernels were used with SVM. Area under the receiver operating curve (auROC) was computed with the R function roc from the pROC library and used as measure of accuracy. The rfe function from the caret library was used to carry out Recursive Feature Elimination (RFE) with SVM. In addition to RFE, which reduces the dimensionality of the parameter space in a way that takes into account correlations between features, an automatic selection of best performing SVM parameters (sigma was pre-estimated using sigest function form kernlab package and fixed for a range of cost C parameter which varied form 0.25 to 32) was used along with leave-one-group-out cross-validation to optimize models and reduce over-fitting.
Regression-based quantitative binding strength predictors were also built and cross validated with the R functions train and trainControl. The TF tag count coverage served as target variable. Cross-data set predictions were made with the R function predict and the relative importance of different features was extracted with varImp (all functions are part of the caret package). All sites from a given table were used in this case. SVM regression with radial basis function kernel was used as training method and a Pearson’s correlation coefficient (PCC) between measured and predicted TF tags was used as performance indicator.
In the next step, these peak lists were annotated with the ChIP-seq tag counts from the corresponding experiments. Tags were counted in a ±100 base region relative to the TFBS center (see Methods for details). These tag counts represent binding strength estimates of an individual TF in different tissues. Note that both predicted binding sites and ENCODE peaks were retrospectively annotated with tag counts. For predicted sites, the same list was annotated with different ChIP-seq data representing different cell-types. In the case of ENCODE peak list, different lists were annotated with different ChIP-seq data pertaining to the corresponding cell types.
In two-class prediction experiments, the goal was to discriminate between strong and weak binding sites. To this end, only the top and bottom 20 % of TF binding sites (in terms of tag coverage) were used. For regression-based quantitative binding strength prediction, all TFBS of a given list were considered.
Each peak list was annotated with a large number of features that could be used as input variables for building predictors. Some of these features were directly derived from genome sequences or general annotation resources and thus invariable between different cell types. These features include a PWM-based TFBS score, base and oligo-nucleotide compositions of the binding regions (comprehensively referred to as “sequence” features, see Methods for more details), DNA structural parameters (10 in total), the distance to the nearest TSS, predicted nucleosome occupancy based on the model by Kaplan et al. , SNP density, and a cross-genome conservation score. Cell-type specific features included ChIP-seq derived relative abundances of 8 histone marks, PolII and co-factors known to be associated with the TFs under investigation (Rad21 for CTCF and FOSL1 for JunD). In addition, two measures of chromatin accessibility were computed from DGF. One of these measures is count-based whereas the other is shape-based (see Methods for details). The annotated peak lists were then used as data input for various machine-learning algorithms. For class prediction experiments, we used the area under the ROC curve as performance measure, for quantitative binding strength prediction, we used a Pearson correlation coefficient (PCC).
To gain preliminary insights into the relative importance and usefulness of the various features, we focused first on CTCF, a kind of model system for ChIP-seq data analysis. Specifically, we built binary classifiers to distinguish strong from weak binding sites using different features and feature classes. Classifiers were built for both predicted TFBS and ENCODE peaks. Data for the K562 cell line were used in this computational experiment.
Interestingly, experimental features such as DGF and histone mark perform better for predicted sites whereas the TFBS score is more useful for discriminating strong from weak binding sites in the ENCODE peak lists. The weaker performance of TFBS on predicted site could be explained by assuming that some of the high-scoring PWM matches reside in closed chromatin regions and thus are simply inaccessible to CTCF. Such “false positives” would be difficult to identify without experimental data, as chromatin accessibility varies substantially between cell types. Moreover, the diversity in TFBS is reduced compared to experimentally determined in vivo binding sites, as the predicted sites are selected for PWM matches with a high score. For analogous reasons, the weaker performance of DGF and histone marks on ENCODE peak lists is not surprising. In simplified terms, in vivo binding sites are sampled from accessible and active chromatin regions, resulting in reduced variability of these features as compared to predicted sites.
The predictive power of the numerous sequence features combined (over 4000 in total, see Methods) was equal or lower than the TFBS score alone. Furthermore, using that many sequence features was computationally expensive and time taking. Therefore we replaced these features in subsequent analyses by a simple G + C content feature, which we knew from parallel studies, was positively correlated with CTCF binding. Since count and shape-based DGF measures showed very similar performance, we kept only the simpler count-based measure. All other features were kept even though some showed very modest performance. The particularly weak performance of SNPs may be explained by the generally low SNP counts (often 0) in TFBS regions.
We carried out RFE, to identify the most informative individual features and to analyze the degree of redundancy between features. RFE works by backward elimination, starting with all the features and eliminating one at a time. At each step, the feature judged to be least useful for prediction is eliminated and the overall performance of the predictor is re-evaluated by cross-validation. The results of recursive feature elimination are shown in Fig. 2b. We note that for both predicted TFBS and ENCODE peaks maximal performance (0.99 accuracy) is reached with as few as three features. The cofactor Rad21 and DGF (shape-based) are part of the top features in both cases, complemented by the PWM score for ENCODE peaks and histone modification H3K4me1 for predicted sites.
In the next experiment, we used a quantitative (regression-based) prediction approach to address the following questions: (i) predict the binding strength of TFs (ii) do the previous observations concerning CTCF generalize to other TFs? (iii) Can a predictor trained on one cell type be applied to another cell type? We thus extended our study to another transcription factor, JunD, and to another cell type, H1hESC, an embryonic stem cell derived cell line. We divided the complete feature collections into two major groups: (i) chromatin-related (cell-type specific) features and (ii) sequence-derived (cell-type independent) features. Predictors were trained on each subset individually and on all the features together. As we found again that SVM with radial kernel performed better than the other two methods tested, we used this method here and in all subsequent analyses described in this paper.
A markedly different picture emerges for JunD. We first note the almost total insignificance of the PWM-based binding score, which largely explains the previously observed global ineffectiveness of sequence-derived features in JunD binding strength prediction. Interestingly, this holds for predicted as well as for ChIP-seq defined binding sites. The intrinsically low binding specificity of JunD reflected by the low information content of the corresponding JASPAR PWM may partly explain this observation. We also tested a secondary PWM for JunD from the Jaspar database (MA0491.2) and observed a marginally better (but still poor) performance. In biological terms, this probably means that sequence-specific protein-DNA interactions play only a minor role among the molecular processes that recruit JunD to its physiological target sites. Another interesting observation is that the cofactor FOSL1 is assigned lower feature importance in predictors built from H1hESC cells as compared to K562 cells. A hypothetical explanation for this unexpected difference could be that JunD preferentially associates with members of the AP1 family other than FOSL1 in H1hESC cell. Generally, active histone marks are of much greater importance for predicting JunD binding strength than CTCF binding strength. Repressive histone marks (H3K27me3) and gene-body marks (H3k36me3 and H4K20me1) are of little or no use for binding strength predictions, independently of cell type.
On the other hand, the recruitment processes seem to be highly factor-specific. The most striking outlier among the five factors analyzed is REST whose binding strength is well predicted by its PWM score. However, all other features appear to be virtually useless for this factor. The fact that all dots lie on the main diagonal in Fig. 6a means that the prediction accuracy obtained with the PWM score cannot be increased by feeding more features to the learning algorithm even though histone modifications alone show a modest predictive value (Fig. 6b, R ~ 0.25). A dynamic interplay between REST binding and histone modifications was reported in a paper by Zheng et al. . However, most of the histone marks selected in this study (except H3K27me3 and H3K36me3) did not show significant alterations with REST binding there. Perhaps a different selection of histone modifications would have shown better predictive ability of REST binding. This suggests that the histone modifications sometimes seen in the vicinity of a REST site are the consequence rather than the cause of REST binding to DNA. The results for REST are in sharp contrast to CTCF, where chromatin accessibility (DGF) and histone marks substantially increase prediction accuracy in all cell types. The other three factors can more or less be placed on a continuous line connecting CTCF to GABP, USF2 and JunD, in this order. The main discriminatory feature between these factors is the PWM score that varies in terms of its predictive value from high for CTFC to very low for JunD.
We presented a quantitative framework to investigate and compare the role of sequence-intrinsic (tissue-invariant) and cell-type specific chromatin features in the biological processes that recruit TFs to their physiological target sites. Surprisingly distinct feature importance landscapes were observed for different TFs. The PWM score reflecting the intrinsic affinity of a TF to DNA was a good predictor of in vivo binding strength only for the two factors possessing an information-rich recognition motif (CTCF and REST). DGF came out high in terms of feature importance for all TFs studied except REST. Optimal performance was achieved with as few as three features suggesting high redundancy among the different features tested. DNA structural features, repressive and gene body associated histone marks were generally of little use in predicting in vivo binding site occupancy and strength. Our study highlights the immense value of large functional genomics data sets such as the ENCODE compendium for studying transcription regulation and illustrates the feasibility and effectiveness of pure in silico approaches in this field.
Computationally intensive tasks were performed at the Vital-IT high-performance computing center of the Swiss Institute of Bioinformatics. SK was supported by Swiss National Science Foundation (grant number CRSII3_154500).
Publication charges for this article have been funded by the Swiss Federal Institute of Technology Lausanne (EPFL). This article has been published as part of BMC Bioinformatics Volume 17 Supplement 1, 2016: Selected articles from the Fourteenth Asia Pacific Bioinformatics Conference (APBC 2016). The full contents of the supplements are available online at http://www.biomedcentral.com/bmcbioinformatics/supplements/17/S1.
- Gordan R, Hartemink AJ, Bulyk ML: Distinguishing direct versus indirect transcription factor-DNA interactions. Genome Res. 2009, 19 (11): 2090-100. 10.1101/gr.094144.109.View ArticlePubMedPubMed CentralGoogle Scholar
- Furey TS: ChIP-seq and beyond: new and improved methodologies to detect and characterize protein-DNA interactions. Nat Rev Genet. 2012, 13 (12): 840-52. 10.1038/nrg3306.View ArticlePubMedPubMed CentralGoogle Scholar
- Arvey A, Agius P, Noble WS, Leslie C: Sequence and chromatin determinants of cell-type-specific transcription factor binding. Genome Res. 2012, 22 (9): 1723-34. 10.1101/gr.127712.111.View ArticlePubMedPubMed CentralGoogle Scholar
- Stormo GD: Modeling the specificity of protein-DNA interactions. Quant Biol. 2013, 1 (2): 115-30. 10.1007/s40484-013-0012-4.View ArticlePubMedPubMed CentralGoogle Scholar
- Badis G, Berger MF, Philippakis AA, Talukder S, Gehrke AR, Jaeger SA, et al: Diversity and complexity in DNA recognition by transcription factors. Science. 2009, 324 (5935): 1720-3. 10.1126/science.1162327.View ArticlePubMedPubMed CentralGoogle Scholar
- Mathelier A, Zhao X, Zhang AW, Parcy F, Worsley-Hunt R, Arenillas DJ, et al: JASPAR 2014: an extensively expanded and updated open-access database of transcription factor binding profiles. Nucleic Acids Res. 2014, 42 (Database issue): D142-7. 10.1093/nar/gkt997.View ArticlePubMedGoogle Scholar
- Sharon E, Lubliner S, Segal E: A feature-based approach to modeling protein-DNA interactions. PLoS Comput Biol. 2008, 4 (8): e1000154-10.1371/journal.pcbi.1000154.View ArticlePubMedPubMed CentralGoogle Scholar
- Florquin K, Saeys Y, Degroeve S, Rouze P, Van de Peer Y: Large-scale structural analysis of the core promoter in mammalian and plant genomes. Nucleic Acids Res. 2005, 33 (13): 4255-64. 10.1093/nar/gki737.View ArticlePubMedPubMed CentralGoogle Scholar
- Tsai ZT, Shiu SH, Tsai HK: Contribution of sequence motif, chromatin state, and DNA structure features to predictive models of transcription factor binding in yeast. PLoS Comput Biol. 2015, 11 (8): e1004418-10.1371/journal.pcbi.1004418.View ArticlePubMedPubMed CentralGoogle Scholar
- Sinha S, Schroeder MD, Unnerstall U, Gaul U, Siggia ED: Cross-species comparison significantly improves genome-wide prediction of cis-regulatory modules in Drosophila. BMC Bioinformatics. 2004, 5: 129-10.1186/1471-2105-5-129.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang J, Zhuang J, Iyer S, Lin X, Whitfield TW, Greven MC, et al: Sequence features and chromatin structure around the genomic regions bound by 119 human transcription factors. Genome Res. 2012, 22 (9): 1798-812. 10.1101/gr.139105.112.View ArticlePubMedPubMed CentralGoogle Scholar
- Luo K, Hartemink AJ: Using DNase digestion data to accurately identify transcription factor binding sites. Pac Symp Biocomput. 2013, 80: 91-Google Scholar
- Liu L, Jin G, Zhou X: Modeling the relationship of epigenetic modifications to transcription factor binding. Nucleic Acids Res. 2015, 43 (8): 3873-85. 10.1093/nar/gkv255.View ArticlePubMedPubMed CentralGoogle Scholar
- Kahara J, Lahdesmaki H: BinDNase: a discriminatory approach for transcription factor binding prediction using DNase I hypersensitivity data. Bioinformatics. 2015, 31 (17): 2852-9. 10.1093/bioinformatics/btv294.View ArticlePubMedGoogle Scholar
- Kim TH, Abdullaev ZK, Smith AD, Ching KA, Loukinov DI, Green RD, et al: Analysis of the vertebrate insulator protein CTCF-binding sites in the human genome. Cell. 2007, 128 (6): 1231-45. 10.1016/j.cell.2006.12.048.View ArticlePubMedPubMed CentralGoogle Scholar
- Rosmarin AG, Resendes KK, Yang Z, McMillan JN, Fleming SL: GA-binding protein transcription factor: a review of GABP as an integrator of intracellular signaling and protein-protein interactions. Blood Cells Mol Dis. 2004, 32 (1): 143-54. 10.1016/j.bcmd.2003.09.005.View ArticlePubMedGoogle Scholar
- Dreos R, Ambrosini G, Perier RC, Bucher P: The eukaryotic promoter database: expansion of EPDnew and new promoter analysis tools. Nucleic Acids Res. 2015, 43 (Database issue): D92-6. 10.1093/nar/gku1111.View ArticlePubMedGoogle Scholar
- Kodama Y, Shumway M, Leinonen R, International Nucleotide Sequence Database C: The sequence read archive: explosive growth of sequencing data. Nucleic Acids Res. 2012, 40 (Database issue): D54-6. 10.1093/nar/gkr854.View ArticlePubMedGoogle Scholar
- Langmead B, Trapnell C, Pop M, Salzberg SL: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10 (3): R25-10.1186/gb-2009-10-3-r25.View ArticlePubMedPubMed CentralGoogle Scholar
- Nair NU, Kumar S, Moret BM, Bucher P: Probabilistic partitioning methods to find significant patterns in ChIP-Seq data. Bioinformatics. 2014, 30 (17): 2406-13. 10.1093/bioinformatics/btu318.View ArticlePubMedGoogle Scholar
- Kaplan N, Moore IK, Fondufe-Mittendorf Y, Gossett AJ, Tillo D, Field Y, et al: The DNA-encoded nucleosome organization of a eukaryotic genome. Nature. 2009, 458 (7236): 362-6. 10.1038/nature07667.View ArticlePubMedGoogle Scholar
- Siepel A, Bejerano G, Pedersen JS, Hinrichs AS, Hou M, Rosenbloom K, et al: Evolutionarily conserved elements in vertebrate, insect, worm, and yeast genomes. Genome Res. 2005, 15 (8): 1034-50. 10.1101/gr.3715005.View ArticlePubMedPubMed CentralGoogle Scholar
- Rosenbloom KR, Armstrong J, Barber GP, Casper J, Clawson H, Diekhans M, et al: The UCSC genome browser database: 2015 update. Nucleic Acids Res. 2015, 43 (Database issue): D670-81. 10.1093/nar/gku1177.View ArticlePubMedGoogle Scholar
- Sherry ST, Ward MH, Kholodov M, Baker J, Phan L, Smigielski EM, et al: dbSNP: the NCBI database of genetic variation. Nucleic Acids Res. 2001, 29 (1): 308-11. 10.1093/nar/29.1.308.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang K, Li M, Hakonarson H: ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010, 38 (16): e164-10.1093/nar/gkq603.View ArticlePubMedPubMed CentralGoogle Scholar
- Merkenschlager M, Odom DT: CTCF and cohesin: linking gene regulatory elements with their targets. Cell. 2013, 152 (6): 1285-97. 10.1016/j.cell.2013.02.029.View ArticlePubMedGoogle Scholar
- Parelho V, Hadjur S, Spivakov M, Leleu M, Sauer S, Gregson HC, et al: Cohesins functionally associate with CTCF on mammalian chromosome arms. Cell. 2008, 132 (3): 422-33. 10.1016/j.cell.2008.01.011.View ArticlePubMedGoogle Scholar
- Zheng D, Zhao K, Mehler MF: Profiling RE1/REST-mediated histone modifications in the human genome. Genome Biol. 2009, 10 (1): R9-10.1186/gb-2009-10-1-r9.View ArticlePubMedPubMed CentralGoogle Scholar
This article is published under license to BioMed Central Ltd. Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.