Drillon Guénola, Audit Benjamin, Argoul Françoise, Arneodo Alain
Univ Lyon, Ens de Lyon, Univ Claude Bernard Lyon 1, CNRS, Laboratoire de Physique, Lyon, F-69342, France.
LOMA, Université de Bordeaux, CNRS, UMR 5798, 51 Cours de le Libération, Talence, F-33405, France.
BMC Genomics. 2016 Jul 29;17:526. doi: 10.1186/s12864-016-2880-2.
Recently, a physical model of nucleosome formation based on sequence-dependent bending properties of the DNA double-helix has been used to reveal some enrichment of nucleosome-inhibiting energy barriers (NIEBs) nearby ubiquitous human "master" replication origins. Here we use this model to predict the existence of about 1.6 millions NIEBs over the 22 human autosomes.
We show that these high energy barriers of mean size 153 bp correspond to nucleosome-depleted regions (NDRs) in vitro, as expected, but also in vivo. On either side of these NIEBs, we observe, in vivo and in vitro, a similar compacted nucleosome ordering, suggesting an absence of chromatin remodeling. This nucleosomal ordering strongly correlates with oscillations of the GC content as well as with the interspecies and intraspecies mutation profiles along these regions. Comparison of these divergence rates reveals the existence of both positive and negative selections linked to nucleosome positioning around these intrinsic NDRs. Overall, these NIEBs and neighboring nucleosomes cover 37.5 % of the human genome where nucleosome occupancy is stably encoded in the DNA sequence. These 1 kb-sized regions of intrinsic nucleosome positioning are equally found in GC-rich and GC-poor isochores, in early and late replicating regions, in intergenic and genic regions but not at gene promoters.
The source of selection pressure on the NIEBs has yet to be resolved in future work. One possible scenario is that these widely distributed chromatin patterns have been selected in human to impair the condensation of the nucleosomal array into the 30 nm chromatin fiber, so as to facilitate the epigenetic regulation of nuclear functions in a cell-type-specific manner.
最近,基于DNA双螺旋序列依赖性弯曲特性的核小体形成物理模型已被用于揭示普遍存在的人类“主要”复制起点附近核小体抑制能垒(NIEB)的一些富集情况。在此,我们使用该模型预测22条人类常染色体上约160万个NIEB的存在。
我们表明,这些平均大小为153 bp的高能垒对应于体外,但同样也对应于体内的核小体缺失区域(NDR)。在这些NIEB的两侧,我们在体内和体外都观察到类似的紧密核小体排列,这表明不存在染色质重塑。这种核小体排列与GC含量的振荡以及沿这些区域的种间和种内突变谱密切相关。对这些分歧率的比较揭示了与这些内在NDR周围核小体定位相关的正选择和负选择的存在。总体而言,这些NIEB和相邻的核小体覆盖了人类基因组的37.5%,其中核小体占有率在DNA序列中稳定编码。这些1 kb大小的内在核小体定位区域在富含GC和GC含量低的等臂染色体、早期和晚期复制区域、基因间区域和基因区域中均有发现,但在基因启动子处未发现。
NIEB上选择压力的来源在未来的工作中仍有待解决。一种可能的情况是,这些广泛分布的染色质模式在人类中被选择,以削弱核小体阵列凝聚成30 nm染色质纤维,从而以细胞类型特异性方式促进核功能的表观遗传调控。