接触广泛丢失,表达变化却很少
黏连蛋白是一种环状蛋白复合物,负责挤出 DNA 环;去除它会在全基因组范围内抹去环和接触结构域。论文的起点是一个定量矛盾:在 HCT-116 结直肠癌细胞中用生长素诱导黏连蛋白缺失 6 小时后,尽管 3D 结构被广泛破坏,只有 146 个基因的表达变化超过 1.75 倍。如果增强子靠与启动子发生物理接触来发挥作用,为什么失去这些接触后大部分转录仍然完好?
作者的回答是:接触丢失并不均匀。它集中在长距离的增强子-启动子对上,而且只有真正丢失接触的增强子才失去调控效力。这把问题从“黏连蛋白是否重要”改写为“哪些增强子-启动子对依赖它”。
活性乘以接触预测哪些对被丢失
为了在缺失前后绘制增强子-基因对,研究使用了活性-接触(ABC)模型。ABC 把每个候选元件-基因对的得分定义为增强子活性乘以增强子-启动子接触频率,并以该基因 5 Mb 内所有其他元件的乘积做归一化。活性来自 DNase-seq 可及性和 H3K27ac ChIP-seq;接触来自匹配的 ENCODE Hi-C,分辨率 5 Kb。若某对的 ABC 得分超过 0.027,就判定为增强子互作;该阈值依据此前 K562 CRISPRi 数据设定,以达到 70% 召回率。
由于两种条件之间活性基本稳定而接触发生变化,ABC 得分的变化紧密跟随接触变化。这使 ABC 成为接触依赖性调控丢失的有用读数,但它仍是预测:模型假设增强子效应随活性乘以接触而缩放,论文用实验检验这一假设,而不是直接接受它。
接触丢失集中在 50 Kb 以外
距离依赖是最清晰的定量结果。黏连蛋白缺失后,相距超过 50 Kb 的增强子-启动子对平均丢失 62.2% 的接触频率,而相距不足 50 Kb 的对丢失 18.6%。作者把 50 Kb 视为接触频率变化约两倍的临界点,并用它区分远端与近端元件。
这种不对称重塑了预测的调控图谱。有黏连蛋白时,ABC 在 10,393 个基因上判定 28,575 个增强子-基因对;缺失后剩下 9,489 个基因上的 21,523 对。相距超过 50 Kb 的对中有 84.0% 被丢失,而不足 50 Kb 的对中只有 23.7%;互作距离的第 90 百分位从 156.3 Kb 降到 69.2 Kb。基于相同数据的逻辑回归替代模型 ENCODE-rE2G 也重现了优先丢失远端增强子的同一趋势。
距离并非全部。不同远端增强子的接触下降幅度各异,并未收敛到单一的距离衰减曲线;靠近 CTCF 位点和远离 CTCF 位点的对都丢失了接触。作者据此认为,环挤出决定长距离接触,而包括区室化在内的其他位点特异机制帮助部分接触在没有黏连蛋白时得以维持。
CRUDO 在有和无黏连蛋白时检验增强子效应
预测需要直接检验,因此作者建立了 CRUDO:在 degron 操作下对调控元件做 CRISPRi。他们把抑制候选元件的混合 CRISPRi 筛选,与在同一细胞系中用生长素诱导降解 RAD21(黏连蛋白的主要亚基)结合起来。目标基因 1 Mb 内的每个元件都被扰动,基因表达效应在有黏连蛋白和无黏连蛋白两种条件下都用 FlowFISH 读出。
该设计给出一个逐元件的指标,论文称之为黏连蛋白依赖性。如果抑制某个增强子在有黏连蛋白时使表达下降 15%,无黏连蛋白时只下降 7.7%,那么这个增强子的一部分效应依赖黏连蛋白。形式上,该指标比较增强子在两种条件下的效应,假设增强子活性与启动子以乘法方式组合;显著性通过针对同一元件的多条 gRNA 做配对 t 检验,并以 Benjamini-Hochberg 校正低于 0.05 判定。
CCND1 的预实验展示了这一逻辑。启动子和两个近端增强子在两种条件下效应相近,黏连蛋白依赖性分别为 −4.9%、−13.3% 和 2.3%。第三个远端增强子(E1)的效应从 15.3% 降到 7.7%,黏连蛋白依赖性为 49.7%,它也是三者中唯一与启动子接触发生变化的增强子。
三十个增强子-基因对,十六个依赖黏连蛋白
完整筛选覆盖了五个黏连蛋白敏感基因周围的 1,039 个候选增强子:CCND1、FAM3C、KITLG、SSFA2 和 MYC。这些基因因高表达(便于检测)以及拥有接触变化范围较广的 ABC 预测增强子而被选中。筛选鉴定出 30 个显著的增强子-基因对,其中 16 个具有黏连蛋白依赖性效应。
解读这个数字时有两个边界。第一,这 30 对来自五个位点而非全基因组,因此比例描述的是这组被选中的基因,而不是所有增强子。第二,所测效应往往很小:讨论部分报告相对基线表达的变化为 5–24%,并用 FlowFISH 和 DC-TAP-seq 两种读数验证。论文把如此温和效应的检出归因于 CRUDO 的灵敏度。
敏感增强子都是远端,但远端未必敏感
核心模式是一种不对称。该集合中每个依赖黏连蛋白的增强子都是远端,距其目标启动子超过 50 Kb。但反过来不成立:许多远端增强子在黏连蛋白缺失后仍保持调控效应。因此距离只是敏感性的前提条件,并不决定敏感性。
与敏感性对应更好的变量是接触频率本身的变化。在 30 个元件中,黏连蛋白依赖性与增强子-启动子 3D 接触的百分比变化相关,论文的回归比较支持接触变化优于单看距离:仅含距离的模型拟合不如包含接触变化的模型。在接触下降超过 50% 的 22 对中,15 对被判定为依赖黏连蛋白,2 对不依赖,4 对检验效能不足。
这是跨元件的相关性,作者也如此表述:接触变化预测效应变化,而区室化等位点特异因素决定哪些远端对保留接触。实际解读是:仅凭远端位置不能把增强子标记为黏连蛋白敏感;测得的接触下降才可以。
为何多数基因几乎不受黏连蛋白缺失影响
全基因组 ABC 比较为最初的矛盾提供了解释。大多数预测的增强子-基因互作在黏连蛋白缺失后并未出现显著接触变化,主要因为多数预测对位于接触得以维持的较短距离。如果增强子仍能接触其启动子,抑制它应当仍会降低表达,因此该基因保持对黏连蛋白不敏感。
按此解释,黏连蛋白敏感基因数量少,并不能证明 3D 接触与增强子功能无关。它反映的是只有少数增强子-启动子对的接触依赖黏连蛋白。作者把这一点与工程位点和合成激活子的早期工作联系起来——那些工作强调距离——并通过展示内源远端增强子在接触丢失程度上差异很大而加以扩展。
接触变化比距离更能预测敏感性
对于想判断某个增强子是否会对黏连蛋白扰动有反应的读者,论文的优先项很明确:看测得的增强子-启动子接触频率变化,而不是只看基因组距离。在该数据集中,超过 50 Kb 的距离是黏连蛋白依赖性的必要条件而非充分条件;接触丢失超过 50% 的 22 对分成 15 对依赖、2 对不依赖和 4 对效能不足,而不是一个干净的“是”。
范围是一个 HCT-116 细胞系中五个黏连蛋白敏感基因、生长素处理 6 小时,效应通过 CRISPRi 抑制候选元件测得。在这一范围内,研究的贡献是一条定量联系:接触丢失预测增强子效应丢失,而多数效应幅度温和,解释了为何全基因组表达变化仍然有限。
本文引用的文献
Quantitative changes in 3D contact predict changes in enhancer-promoter regulation upon cohesin depletion
黏连蛋白缺失对超过 50 Kb 的增强子-启动子接触削弱最强,CRUDO 显示丢失接触的增强子也丢失调控效应,而全基因组大多数预测对保持接触。
