Abstract:Genome-wide association studies (GWAS) have linked thousands of genetic variants to various complex traits or diseases. However, most identified variants have weak individual effects, are correlated with nearby polymorphisms due to linkage disequilibrium (LD), and are located in non-coding cis-regulatory elements (CREs). These characteristics complicate the assessment of the direct impact of each variant on tissue specific gene expression and phenotype. To address this challenge, we have developed a novel algorithm that leverages polymer folding and 3D chromatin interactions to prioritize and identify putative causal variants and their target genes. From the millions of eQTL-Gene pairs identified by GTEx in human somatic tissues, we classify only ~10-20% as putative functional eQTL-Gene pairs supported by phenotypic associations confirmed through CRISPR deletion experiments. Our findings show that unlike most variants, functional eQTL-Gene pairs predominantly reside within the same topologically associating domain (TAD) and have strong associations with cell-type specific cis-regulatory elements (CREs), enriched for binding sites of tissue-specific transcription factors. Unlike most approaches that rely on linear distance or other chromatin features (histone code, accessibility), our algorithm emphasizes the importance of physical interactions and 3D chromatin folding in gene regulation, as the identified eQTL-Gene pairs are all among the small fraction of physical chromatin interactions sufficient for chromatin locus folding. Overall, our algorithm reduces false positive associations between DNA variants and genes identified by eQTL analysis and uncovers novel variant-gene pair associations. These findings suggest a mechanism where a small number of regulatory variants control tissue specific gene expression via their physical association with target genes confined within the same TAD. Our approach provides new insights into the molecular mechanisms driving GWAS phenotypes.

Integrative modeling of eQTLs and cis-regulatory elements suggest mechanisms underlying cell type specificity of eQTLs

Haplotype-aware modeling of cis-regulatory effects highlights the gaps remaining in eQTL data

A unified framework for cell-type-specific eQTLs prioritization by integrating bulk and scRNA-seq data

The Single-Cell Eqtlgen Consortium

Research progress on single-cell expression quantitative trait loci

Single-cell Eqtlgen Consortium: a Personalized Understanding of Disease

Promoter-interacting expression quantitative trait loci are enriched for functional genetic variants

A New Approach for Discovering Functional Links Connecting Non-Coding Regulatory Variants to Gene Targets

Impact of disease-associated chromatin accessibility QTLs across immune cell types and contexts

Cross-phenotype meta-analysis reveals large-scale trans-eQTLs mediating patterns of transcriptional co-regulation

Structural Basis of Differential Gene Expression at eQTLs Loci from High-Resolution Ensemble Models of 3D Single-Cell Chromatin Conformations

Hierarchical annotation of eQTLs by H-eQTL enables identification of genes with cell type-divergent regulation

Mapping enhancer-gene regulatory interactions from single-cell data

A statistical framework for joint eQTL analysis in multiple tissues

The impact of cell type and context-dependent regulatory variants on human immune traits

Epigenome-augmented eQTL-hotspots reveal genome-wide transcriptional programs in 36 human tissues

Probing the limits of cis-acting gene regulation using a model of allelic imbalance quantitative trait loci

Integrating single cell expression quantitative trait loci summary statistics to understand complex trait risk genes

Discovering single-cell eQTLs from scRNA-seq data only

A computational method for cell type-specific expression quantitative trait loci mapping using bulk RNA-seq data

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