Nuleic Acids Res：三维基因组数据的动态染色质可接近性模型

 来自北京大学生命科学学院、生命科学联合中心、统计科学中心李程研究组使用基于Hi-C数据的染色质区域间的相对空间距离，利用布朗运动描述转录因子和染色质修饰蛋白在不同染色质区域之间的动态转移过程，并以此建立了染色质区域附近分子浓度变化和趋于稳定的马尔可夫模型。并指出干细胞分化过程中的一个重要因素是时空特异地调控染色质的空间组织和结构，从而促进细胞阶段和类型特异性的基因得以表达。这一研究成果公布在《Nucleic Acids Research》杂志上。

真核细胞的染色质在细胞核中折叠成复杂的三维结构，但这些结构的形成和功能还没有被完全理解。近年来，以Hi-C为代表的染色质构象捕获技术揭示了较精细的染色质三维结构单元，如A、B间隔(A/B Compartments)和拓扑结构域(Topologically Associated Domains , TADs)，它们在不同细胞类型、物种内都有保守性，是定量化地理解染色质三维结构形成和功能的关键性进展。三维基因组与基因表达调控具有密切的关联，例如，一个TAD内的基因往往被共同调控，但一个TAD内的的增强子由于被TAD边界阻隔，不会调控另一个TAD内的基因。在疾病形成和干细胞分化过程中，表达调控、表观遗传的变化经常伴随着三维基因组的结构变化。

在此项研究中，李程研究组使用基于Hi-C数据的染色质区域间的相对空间距离，利用布朗运动描述转录因子和染色质修饰蛋白在不同染色质区域之间的动态转移过程，并以此建立了染色质区域附近分子浓度变化和趋于稳定的马尔可夫模型。模型由此推导出的平衡态分布(Steady-State Distribution, SSD)表征了分子在染色质区域之间由于随机运动达到的动态平衡的浓度，是领域内首次基于三维基因组信息的动态染色质可接近性的度量。SSD不同于直观空间结构意义上的紧密程度，因为紧密或松散的染色质区域上都有可能出现高SSD值。

李程研究组进一步发现，相比传统的测量染色质可接近性的实验方法如DNase-seq和FAIRE-seq，SSD与表征染色质活跃程度的多个组蛋白修饰、以及基因表达有更强的相关性,证明了SSD作为动态可接近性度量的优势。SSD模型说明染色质修饰蛋白和转录因子的分布遵循由染色质三维结构决定的动态平衡状态，其浓度则进一步影响局部基因转录表达的水平。此研究还发现SSD可以帮助揭示干细胞分化过程中染色质结构变化的两个阶段。在胚胎干细胞经过神经元前体细胞阶段、向神经元细胞分化的过程中，与染色质结构组织相关的基因所在区间的SSD（动态可接近性）首先升高，帮助这些基因的表达，进而促进更多神经元细胞特异基因所在区间的SSD的升高和相应基因的表达。这说明干细胞分化过程中的一个重要因素是时空特异地调控染色质的空间组织和结构，从而促进细胞阶段和类型特异性的基因得以表达。

该研究得到北大－清华生命科学联合中心、科技部国家重点研发计划、国家自然科学基金重点项目的资助。李程研究员为该论文的通讯作者，生命科学中心13级博士生王轶楠为该论文第一作者，生科院14级本科生范操琦是论文第二作者，生命科学中心16级博士生郑宇轩参与了部分工作。