Systematic assessment of sequencing depth requirements for Hi-C-derived metrics
Abstract
Background: Hi-C experiments produce genome-wide chromatin contact maps from which structural features can be quantified across multiple genomic scales, ranging from megabase-scale compartments to kilobase-scale boundaries and chromatin loops. Although high-resolution analyses commonly rely on hundreds of millions of sequenced read pairs, the minimum sequencing depth required for different classes of Hi-C-derived features has not been systematically established. We therefore sought to determine these depth requirements systematically.
Results: We performed progressive random subsampling of twelve Hi-C libraries representing multiple cell types and experimental protocols. Loop-density and loop-size inference from the P(s) log-derivative was evaluated across all libraries, whereas compartment, insulation, and boundary analyses were performed on a subset of nine libraries with comparable full-depth coverage. Rather than defining a “sufficient depth” as an arbitrary fraction of the full-depth value, we introduced a biologically motivated criterion based on the variability of the metric between independent biological replicates. The required sequencing depth for each metric was defined as the point at which the subsampling scatter first reached the variability observed between independent biological replicates, representing the accuracy that additional sequencing cannot improve upon. Using this criterion, the genome-average loop-density estimate reached its biological floor at approximately 7 million read pairs, the genome-average loop-size estimate at approximately 20 million, insulation scores at approximately 80 million, and compartment eigenvectors at 100-kilobase resolution at approximately 110 million read pairs.
Conclusions: Different Hi-C-derived metrics require substantially different sequencing depths to achieve biologically meaningful accuracy. For all metrics, depth-induced variability fell below biological replicate variability well before full sequencing depth was reached. These thresholds apply to genome-averaged characteristics, not local metrics of specific loci. They therefore provide practical guidance for genome-wide experimental design and sequencing budget allocation, suggesting that increasing the number of biological replicates is likely to improve the inference of global biophysical features more effectively than sequencing individual libraries to greater depth.
Related
- FHIR Data Mapping Quality Assessment: A Scoping Review and Expert-Informed Requirements for Future Pre-Mapping Frameworks.
- Systematic evaluation of integration methods and parameters on single-cell RNA-sequencing biological insights: a case study on cattle embryos
- Refining the TCGA Framework: The Pivotal Role of Targeted Next-Generation Sequencing in Resolving Diagnostic Ambiguities in Endometrial Cancer—Systematic Review
- Systematic quantification and removal of host DNA contamination in 16S rRNA gene sequencing
- Controllable Cascade Aggregation Ion Programmed Nanomachines and Simple Isolation Enable Tumor‐Derived Small Extracellular Vesicles Detection for Enhanced Breast Cancer Assessment
BibTeX
@article{granovsky2026systematic,
title = {{Systematic assessment of sequencing depth requirements for Hi-C-derived metrics}},
author = {Alexander Granovsky and Kirill Polovnikov},
journal = {BMC Genomics},
year = {2026},
month = oct,
doi = {10.1186/s12864-026-13455-7},
url = {https://link.springer.com/article/10.1186/s12864-026-13455-7},
note = {CC BY 4.0},
}
RIS
TY - JOUR
TI - Systematic assessment of sequencing depth requirements for Hi-C-derived metrics
AU - Alexander Granovsky
AU - Kirill Polovnikov
PY - 2026
DA - 2026/10/15
JO - BMC Genomics
DO - 10.1186/s12864-026-13455-7
UR - https://link.springer.com/article/10.1186/s12864-026-13455-7
AB - Background: Hi-C experiments produce genome-wide chromatin contact maps from which structural features can be quantified across multiple genomic scales, ranging from megabase-scale compartments to kilobase-scale boundaries and chromatin loops. Although high-resolution analyses commonly rely on hundreds of millions of sequenced read pairs, the minimum sequencing depth required for different classes of Hi-C-derived features has not been systematically established. We therefore sought to determine these depth requirements systematically. Results: We performed progressive random subsampling of twelve Hi-C libraries representing multiple cell types and experimental protocols. Loop-density and loop-size inference from the P(s) log-derivative was evaluated across all libraries, whereas compartment, insulation, and boundary analyses were performed on a subset of nine libraries with comparable full-depth coverage. Rather than defining a “sufficient depth” as an arbitrary fraction of the full-depth value, we introduced a biologically motivated criterion based on the variability of the metric between independent biological replicates. The required sequencing depth for each metric was defined as the point at which the subsampling scatter first reached the variability observed between independent biological replicates, representing the accuracy that additional sequencing cannot improve upon. Using this criterion, the genome-average loop-density estimate reached its biological floor at approximately 7 million read pairs, the genome-average loop-size estimate at approximately 20 million, insulation scores at approximately 80 million, and compartment eigenvectors at 100-kilobase resolution at approximately 110 million read pairs. Conclusions: Different Hi-C-derived metrics require substantially different sequencing depths to achieve biologically meaningful accuracy. For all metrics, depth-induced variability fell below biological replicate variability well before full sequencing depth was reached. These thresholds apply to genome-averaged characteristics, not local metrics of specific loci. They therefore provide practical guidance for genome-wide experimental design and sequencing budget allocation, suggesting that increasing the number of biological replicates is likely to improve the inference of global biophysical features more effectively than sequencing individual libraries to greater depth.
SN - 1471-2164
ER -