一文带你get Hi

作者：杨超月

DNA在染色体上是高度折叠的，DNA与DNA片段之间不可避免的形成了高强度的交互作用。2002年提出的3C技术，用于测定染色体特定位点与位点之间的交互作用，随后发展出4C、5C 技术，分别用于测定染色体上一点到多点和多点与多点之间的交互作用。Job Dekker在2009年开发出Hi-C 技术，实现了全基因组范围内染色体片段间相互作用的捕获。Hi-C 除了用于辅助基因组组装、对已组装的基因组进行纠错，还可以揭示基因组的三维结构特征，包括从隔室（A/B Compartments）到拓扑相关结构域（TAD），最后到染色体环（loop），目前，已在人、果蝇、酵母、拟南芥、水稻和棉花等物种成功构建基因组三维结构，并完成了对不同样本基因组三维结构的比较分析。

1.Hi-C三维构象研究内容

1）绘制全基因组的Hi-C map、Cis/Trans和IDEs等分析；
2）Compartment A/B 鉴定与分析，包括与ChIP-seq、WGBS和RNAseq等联合分析；
3）TAD鉴定与分析，包括与ChIP-seq、WGBS和RNAseq等联合分析；
4）绘制全基因组LOOP模型图，需DNase-seq或ATAC-seq，并与ChIP-seq和RNAseq联合分析；
5）多样本三维结构差异分析，包括Compartment A/B、TAD和Loop差异鉴定与分析；
6）基因/重复序列互作分析，包括与RNAseq联合分析。

一文带你get Hi

图1 染色体不同层级三维构象图
（Fraser J, et al, 2015）

2.染色体三维结构重点研究内容

A/B compartments

高等动植物基因组可以人为划分为两个compartments，称作A/B compartments，通过计算PCA，第一主成分中正值表示A compartment，负值表示B compartment。A/B compartments在基因组中一般呈间隔分布，A compartments呈松散染色质状态，基因密度高，属于转录活跃区域，而B compartments则呈压缩染色质状态，基因密度低，属于转录抑制区域。通过对多个样本A/B compartments 的分布进行对比，在全基因组水平找出A/B compartments保守的区域以及发生A/B compartments转换的区域，结合RNA-Seq对发生A/B compartments转换区域内的基因的表达量进行统计，有助于解释不同样本间的差异与染色质活性状态之间的关系。如图2：

一文带你get Hi

图2 人正常细胞和癌症细胞A/B compartments
（Wu P, et al, 2017）

例如在人的多发性骨髓瘤文章中，作者利用正常B细胞（GM12878）和两种亚型的骨髓瘤细胞（近二倍体U266和近三倍体RPMI-8226）进行Hi-C、WGS和RNA-seq测序，结果发现，相比于GM12878，U266和RPMI-8226癌症细胞中6%的区域均存在B转换到A，1%的区域均存在A转换到B，结合RNA-seq数据，发现B转换到A区域内的基因变成转录活跃状态，而A转换到B区域内的基因变成转录抑制状态，KEGG富集通路发现这些基因与多发性骨髓瘤相关，主要涉及MAPK信号转导途径、TNF信号转导途径、细胞因子和细胞因子受体相互作用途径等，进一步在染色体2q11.2-q12.142处发现一个细胞因子基因簇，几种白细胞介素IL1R1，IL1R2，IL18R1和细胞因子MAP4K4下调表达，这些基因在后续研究中重点关注。

TAD

TAD（拓扑关联结构域）是一个高度自关联的连续区域，通过明显的边界与相邻区域分离开来，形成一个独立的调控单元，内部的基因拥有共同的调控元件，存在协同表达特征。TAD边界通常具有大量的绝缘子蛋白和黏连蛋白（植物中TAD边界一般缺少绝缘蛋白，边界不明显），对于维持TAD结构及稳定性具有重要作用，不但可以指导染色体折叠成高级结构，还可以正确指导远距离转录调控，该边界发生变化会导致基因调控变得紊乱。TAD边界还与组蛋白修饰、甲基化修饰等密切相关，通常与转录活性相关的表观遗传标记富集在TAD边界，而与转录失活相关的表观遗传标记远离TAD边界。通过对多个样本TAD进行对比，在全基因组水平找出发生变化的TAD（重点关注TAD边界），结合ChIP-seq、WGBS等分析TAD的变化是否与表观遗传修饰相关，进一步利用RNA-Seq对相关基因的表达量进行统计，有助于解释不同样本间空间结构的差异与表观遗传修饰及转录调控之间的关系。如图3：

一文带你get Hi

图3 不同种棉花TAD边界及表观遗传修饰
（Wang M, et al, 2018）