Assessing base-resolution DNA mechanics on the genome scale

Abstract

Intrinsic DNA properties such as bending play a crucial role in diverse biological systems. A recent advantage in the high-throughput method called loop-seq makes it possible to determine bendability of hundred thousand 50-bp DNA duplexes in one experiment. However, its still infeasible to assess whole sequence bendability in large genomes such as human, which needs thousands of loop-seq experiments. Here we introduce BendNet - a neural network to accurately predict the intrinsic DNA bending at base-resolution by only given DNA sequences. BendNet can increase the resolution of experimental results, and can predict DNA bendability for any new given sequences in high accuracy. We applied BendNet to the human genome and observed high-stiffness regions located at both transcriptional start sites and transcriptional end sites. Such stiffness patterns are different for coding and non-coding genes, which matches distinct nucleosome occupancy patterns. As expected, most transcription factors (TFs) bind in DNA of low bendability. In contrast, we observed an unusually high bendability within binding elements of specific TFs such as EBF1 and regulators of genome folding such as CTCF. These factors either co-bind or compete with nucleosomes to carry out their functions. More interestingly, CTCF binding regions exhibit the highest bendability than other DNA elements, implying their potential role in trapping and holding the CTCF in the exact locations to make sure CTCF as stable anchor in loop extrusion process. Our work provides a tool to assess DNA bendability for large-scale DNA sequences and expands our understanding on DNA mechanics in chromatin regulation and genome folding.