The adapted Activity-By-Contact model for enhancer-gene assignment and its application to single-cell data

In our work, we implemented the ABC-model and could show that one assay for measuring the openness of enhancers is sufficient. Further, we propose an adapted calculation of the ABC-score, which describes enhancer activity in a gene-specific manner without requiring any additional data. We combined our implementation of the ABC-score with an approach to quantify TF binding affinity into STARE: a framework to derive TF affinities to genes. STARE was also designed for potential application on single-cell data. You can find the code and more details here: https://github.com/schulzlab/stare We provide the data for the validation of our ABC-implementation on two CRISPR-screens. We also provide the results of our analysis of single-cell data of the human heart with STARE. All data is in hg19. Folder structure: <strong>K562_CandidateEnhancer</strong>: K562 enhancer with the 4th column for enhancer activity, one file for each activity representation that was measured. <strong>K562_ABC_Predictions: </strong>Regular ABC-scores and adapted ABC-scores for each activity measurement. The files contain all scored interactions for a 10MB window, without any cut-off. We also included the results of the implementation of the ABC-score of Fulco et al. (2019). <strong>STARE_Hocker_*</strong>: Whole STARE output for human heart single-cell data, one for regular ABC, adapted ABC and one for window-based approach. All approaches were run with a 5 MB window, the ABC-based runs with a score cut-off of 0.02. Each folder contains two subdirectories, one for the ABC-scoring and one for the Gene-TF affinity matrices. The 'ABC_output' also contains a GeneInfo file for each cell type, summarising different attributes per gene. <strong>INVOKE_Hocker_*</strong>: Folder with the input and output of INVOKE (see https://github.com/schulzlab/tepic), based on the STARE runs. CS genes stands for cell type-specific genes, defined as genes with a z-score across cell types of ≥ 2 and TPM ≥ 0.5. The INVOKE commands were as follows:<br> Rscript INVOKE.R --dataDir=&lt;TF-Gene matrix&gt; --outDir=&lt;out_path&gt; --response=Expression --regularization=E --performance=TRUE --seed=1234 Importantly, the results are based on data of the following publications: K562 predictions: Fulco, C. P. et al. (2019). Activity-by-contact model of enhancer–promoter<br> regulation from thousands of CRISPR perturbations. Nature Genetics,<br> 51(12), 1664–1669 Hi-C matrix for K562 predictions: Rao, S. et al. (2014). A 3D Map of the Human Genome at Kilobase<br> Resolution Reveals Principles of Chromatin Looping. Cell, 159(7),<br> 1665–1680 STARE and INVOKE runs: Hocker, J. D. et al. (2021). Cardiac cell type–specific gene regula-<br> tory programs and disease risk association. Science Advances, 7(20),<br> eabf1444 H3K27ac HiChIP for STARE runs: Anene-Nzelu, C. G. et al. (2020). Assigning Distal Genomic Enhancers<br> to Cardiac Disease–Causing Genes. Circulation, 142(9), 910–912