Abstract:Methods that analyze single-cell paired RNA-seq and ATAC-seq multiome data have shown great promise in linking regulatory elements to genes. However, existing methods differ in their modeling assumptions and approaches to account for biological and technical noise-leading to low concordance in their linking scores-and do not capture the effects of genomic distance. We propose pgBoost, an integrative modeling framework that trains a non-linear combination of existing linking strategies (including genomic distance) on fine-mapped eQTL data to assign a probabilistic score to each candidate SNP-gene link. We applied pgBoost to single-cell multiome data from 85k cells representing 6 major immune/blood cell types. pgBoost attained higher enrichment for fine-mapped eSNP-eGene pairs (e.g. 21x at distance >10kb) than existing methods (1.2-10x; p-value for difference = 5e-13 vs. distance-based method and < 4e-35 for each other method), with larger improvements at larger distances (e.g. 35x vs. 0.89-6.6x at distance >100kb; p-value for difference < 0.002 vs. each other method). pgBoost also outperformed existing methods in enrichment for CRISPR-validated links (e.g. 4.8x vs. 1.6-4.1x at distance >10kb; p-value for difference = 0.25 vs. distance-based method and < 2e-5 for each other method), with larger improvements at larger distances (e.g. 15x vs. 1.6-2.5x at distance >100kb; p-value for difference < 0.009 for each other method). Similar improvements in enrichment were observed for links derived from Activity-By-Contact (ABC) scores and GWAS data. We further determined that restricting pgBoost to features from a focal cell type improved the identification of SNP-gene links relevant to that cell type. We highlight several examples where pgBoost linked fine-mapped GWAS variants to experimentally validated or biologically plausible target genes that were not implicated by other methods. In conclusion, a non-linear combination of linking strategies, including genomic distance, improves power to identify target genes underlying GWAS associations.

Enhlink infers distal and context-specific enhancer–promoter linkages

Linking enhancers to genes

Single-cell multi-ome regression models identify functional and disease-associated enhancers and enable chromatin potential analysis

Demystifying the secret mission of enhancers: linking distal regulatory elements to target genes

Global View Of Enhancer-Promoter Interactome In Human Cells

Enhancer–promoter interactions are encoded by complex genomic signatures on looping chromatin

EXTRA-seq: a genome-integrated extended massively parallel reporter assay to quantify enhancer-promoter communication

Predicting Enhancers with Deep Convolutional Neural Networks

iEnhancer-ELM: improve enhancer identification by extracting position-related multiscale contextual information based on enhancer language models

Linking regulatory variants to target genes by integrating single-cell multiome methods and genomic distance

Predicting unrecognized enhancer-mediated genome topology by an ensemble machine learning model

DeepEnhancer: Predicting Enhancers by Convolutional Neural Networks.

iEnhancer-XG: interpretable sequence-based enhancers and their strength predictor

Supervised learning of enhancer-promoter specificity based on genome-wide perturbation studies highlights areas for improvement in learning

Enhancer-MDLF: a novel deep learning framework for identifying cell-specific enhancers

CelLink: integrating single-cell multi-omics data with weak feature linkage and imbalanced cell populations

Shared and distinct molecular effects of regulatory genetic variants provide insight into mechanisms of distal enhancer-promoter communication

Supervised enhancer prediction with epigenetic pattern recognition and targeted validation

A comparison of experimental assays and analytical methods for genome-wide identification of active enhancers

Exploiting sequence-based features for predicting enhancer–promoter interactions

Genomic context sensitizes regulatory elements to genetic disruption