Continuous Topological Accessibility (RFA) as a 3D-Aware GRN Filter: Anticorrelation with eQTL and No Prediction of ΔRNA
Timofeev AndreyiD, Alexander BratishchikoviD, Alexander AnufrieviD
Abstract
We formalize Relative Forest Accessibility (RFA), a topological accessibility index computed as the solution of (I - {gamma}P )x = v0, where P is a Markov transition matrix built from Hi-C contacts and v0 is a vector of TF peaks (or ATAC peaks). RFA yields a continuous accessibility profile across all genomic bins. We test RFA on three targets. First, on promoter TF counts (GM12878 chr1, 10 kb, 811 genes), RFA loses to 1D distance (AUC 0.836 vs 0.976); this is an artifact because the target and the 1D distance look at the same window. Second, on eQTL (GTEx v10 LCL, chr1+17+19, 10 kb, 1625 genes), RFA anticorrelates with eQTL: {rho}(RFA, - log10 pmin) = -0.09, AUC(RFA) = 0.386 (below chance), AUC(-RFA) = 0.614. High RFA corresponds to fewer eQTL, and dose-response analysis confirms that AUC(RFA) decreases with eQTL strength (mean trend -0.056, 19/20 TFs). A distance-preserving shuffle of the Hi-C matrix reproduces this anticorrelation ({triangleup}AUC vs real: [-0.074, +0.076] across 5 TFs, all CIs include 0), whereas a degree-preserving rewiring destroys it ({triangleup}AUC [isin] [-0.261, -0.003], CIs exclude 0 for 5/5 TFs); the anticorrelation therefore reflects distance-decay of Hi-C contacts (local 3D accessibility), not higher-order topology. The anticorrelation replicates on chr17 ({triangleup}{rho} = -0.13 for rewired, 0/5 for shuffled) and on chr19 with an inverted control pattern, consistent with local within-diagonal heterogeneity as the carrier on short gene-dense chromosomes. Third, on {triangleup}RNA (MCF10A vs MCF7, 40 kb, ATAC peaks), RFA does not predict {triangleup}RNA ({rho} = +0.005, permutation p = 0.87, bootstrap 95% CI [-0.040, +0.053], n = 1730); the CI includes zero and does not exclude the practical boundary |{rho}| = 0.05. Within-line sanity checks are also null. We conclude that RFA anticorrelates with eQTL via local distance-decay and does not predict {triangleup}RNA. The anticorrelation is carried by distance-decay, not by loops or TADs, and may reflect expression buffering: genes with strong eQTL are less dependent on local 3D accessibility. RFA is not a universal predictor but a filter with an inverted sign on eQTL and no signal on {triangleup}RNA.
02 /Specification & Classification
Click any tag to browse matching papers- Organism
- Assay
- Research Focus
- Approach
Reproducibility ledger
Nobody has reported trying to reproduce this paper yet. Did you replicate an experiment or run the code? Your report saves the next lab a week.
Sign in to file a reproducibility report.
Related papers
- ReQTL - an allele-level measure of variation-expression genomic relationships
Liam Flinn Spurr, Nawaf Alomran, … Anélia Horvath·bioRxiv·Nov 7, 2018
- Improved Quality Metrics for Association and Reproducibility in Chromatin Accessibility Data Using Mutual Information
Cullen Roth, Vrinda Venu, … Shawn R. Starkenburg·BMC BioinformaticsPeer-reviewed·Nov 21, 2023
- Chromatin co-accessibility is highly structured, spans entire chromosomes, and mediates long range regulatory genetic effects
William W. Young Greenwald, Agnieszka D’Antonio‐Chronowska, … Kelly A. Frazer·bioRxiv·Apr 9, 2019
- Changes in chromatin accessibility are not concordant with transcriptional changes for single-factor perturbations
Karun Kiani, Eric M. Sanford, … Arjun Raj·Molecular Systems BiologyPeer-reviewed·Sep 1, 2022
- An automated ATAC-seq method reveals sequence determinants of transcription factor dose response in the open chromatin
Betty B. Liu, Masaru Shimasawa, … William James Greenleaf·bioRxiv·Jul 24, 2025
- Hierarchical Markov Random Field model captures spatial dependency in gene expression, demonstrating regulation via the 3D genome
Naihui Zhou, Iddo Friedberg, Mark S. Kaiser·bioRxiv·Dec 16, 2019