我们报道了在由胚胎干细胞生成的拟胚体中应用定量染色体构象捕获技术并结合高通量测序的方法。该技术可用于鉴定和定量分析在胚胎干细胞分化过程中特定基因的潜在增强子与启动子区域之间的相互作用。
我们报道了在由胚胎干细胞生成的拟胚体中应用定量染色体构象捕获技术并结合高通量测序的方法。该技术可用于鉴定和定量分析在胚胎干细胞分化过程中特定基因的潜在增强子与启动子区域之间的相互作用。
在哺乳动物发育过程中,细胞命运的决定依赖于调控网络的建立,这些网络定义了基因表达的特异性、时间和空间模式。由多能干细胞衍生的拟胚体(embryoid bodies, EBs)已成为研究三大胚层分化以及细胞命运决定过程中调控回路的常用模型。尽管已知组织特异性增强子通过与启动子相互作用在这些网络中发挥重要作用,但将增强子准确分配到其对应的靶基因仍具有挑战性。为实现这一目标,需要采用定量方法来研究发育过程中增强子-启动子之间的相互作用及其动态变化。本文中,我们改进了一种4C技术,用于在EB分化模型中鉴定增强子及其与对应启动子的相互联系。该方法采用高频切割的限制性内切酶、超声破碎,以及一种可与商业化DNA文库构建试剂盒兼容的嵌套连接介导PCR流程。随后,对4C文库进行高通量测序并进行生物信息学分析,从而检测和量化所有与选定启动子存在相互作用的序列。所获得的测序数据还可用于获取分化过程中增强子-启动子相互作用动态变化的信息。本研究所描述的EB分化模型中的技术易于实施。
在小鼠中,3.5天龄胚胎的内细胞团(ICM)含有胚胎多能干细胞。ICM在第4.5天进一步发育为上胚层,产生外胚层、中胚层和内胚层细胞,即胚胎中的三个主要胚层。尽管ICM中的多能细胞在体内仅短暂存在,但可通过建立小鼠胚胎干细胞(mESCs)在体外培养中被捕获1,2,3。mESCs能够保持未分化状态并无限增殖,但在内在和外在刺激下,也可退出多能性状态,分化为三个发育胚层的细胞2,4。有趣的是,当在悬浮的小液滴中培养时,mESCs会形成三维聚集体(即拟胚体,EBs),并进一步分化为所有三个胚层的细胞5。拟胚体形成实验是研究早期谱系特化过程的重要工具。
在谱系特化过程中,各胚层的细胞会获得特定的基因表达程序4。基因的精确时空表达受到多种顺式调控元件的调控,包括核心启动子、增强子、沉默子和绝缘子6,7,8,9。增强子是一类通常跨越数百个碱基对的调控性DNA片段,可协调组织特异性的基因表达8。转录因子及其共因子通过结合增强子并调控局部染色质结构,从而激活或沉默增强子功能8,10。目前常用的鉴定潜在增强子的技术包括全基因组染色质免疫沉淀测序(ChIP-seq)以及利用测序技术检测转座酶可及性染色质(ATAC-seq)。因此,活性增强子通常具有特定的活性组蛋白修饰标记以及局部DNA可及性的增加11,12,13,14。此外,发育相关增强子被认为需要与其对应的启动子发生物理相互作用8,9。事实上,已有研究表明,破坏增强子-启动子相互作用的增强子变异或缺失可导致发育畸形15。因此,亟需开发新的技术,以提供额外信息用于识别调控发育过程中基因表达的功能性增强子。
自从染色体构象捕获(3C)技术16发展以来,染色体相互作用位点的图谱绘制已被广泛用于评估调控元件之间的物理距离。重要的是,近年来已开发出多种高通量版本的3C技术,提供了针对染色质片段间交联、消化、连接和相互作用回收的不同策略17。其中,原位Hi-C技术已成为一种流行方法,可实现全基因组范围内3C连接产物的测序18。然而,要达到适合分析增强子-启动子相互作用的分辨率,所需测序成本较高,使得该技术在特定基因座研究中难以应用。因此,研究人员开发了其他替代方法,以更高分辨率分析靶向基因座19,20,21,22。其中一种方法称为4C,即“一个对所有”策略,能够检测与选定视角位点发生相互作用的所有序列。然而,标准4C技术的一个缺点是需要进行反向PCR,该过程会扩增不同大小的片段,倾向于扩增较小产物,从而在高通量测序后的定量分析中引入偏差。最近,一种名为UMI-4C的新版4C技术被开发出来,该方法利用唯一分子标识符(UMI),实现了定量且靶向的染色体相互作用分析,克服了上述问题23。该方法采用高频切割酶、超声破碎以及嵌套连接介导的PCR流程,从而扩增出长度分布相对均一的DNA片段。这种均一性减少了PCR扩增过程中对较短序列的偏好性偏差,并有助于高效回收空间上相连的分子/片段,实现精确计数。
本文介绍了一种改进的UMI-4C技术方案,用于在拟胚体分化过程中识别并定量谱系指导性转录因子的启动子与增强子之间的染色质互作。
1. 小鼠胚胎干细胞生成拟胚体
2. 胚状体的解离
3. 固定
4. 细胞裂解与限制性内切酶消化
5. 邻近连接与交联逆转
6. DNA 打断与片段大小筛选
7. 测序文库制备
8. 4C 染色质相互作用文库的扩增与纯化
在悬浮液滴中诱导胚胎干细胞分化后的第六天,我们获得了均一的拟胚体群体,用于后续分析(图1)。我们改进了UMI-4C方法23,以定量拟胚体中谱系特异性基因启动子区域的特定染色质相互作用24。该实验流程的示意图及各步骤代表性质量控制凝胶图谱如图2A所示。第一次质量控制用于检测MboI限制性内切酶消化效率,高效消化应产生小于3 kbp的片段(图2B)。值得注意的是,小鼠胚胎干细胞和拟胚体的染色质消化较为困难,有时会残留未完全消化的染色质。第二次质量控制在连接反应后进行,以确认大部分片段大小现已> 3 kbp(图2B)。随后,通过凝胶电泳分析超声处理后获得的染色质片段,预期片段大小为400–500 bp(图2B)。
去磷酸化并连接单端接头后,进行两轮PCR以扩增目标片段。采用巢式策略为每个位点设计一对两种引物,以提高特异性。每个目标片段分别使用两对不同的引物进行扩增,以优化PCR条件(例如,Pou5f1 位点使用引物对A和B,T 位点使用引物对C和D),扩增产物在约400 bp处呈现DNA弥散条带(图2C)。 alternatively,采用多重PCR同时扩增目标A和C(图2D),纯化后获得相似大小的片段(图2D)。用于4C文库构建的引物(针对 Pou5f1 和 T 位点)见 表6。
数据分析时,首先将原始测序读段比对至参考小鼠基因组(mm10),去除所有重复及低质量(< 20)的读段。针对每个诱饵序列,通过计算读段片段的数量获得每个限制性酶切片段的信息,从而得到原始接触图谱。随后,将感兴趣区域定义为距离诱饵序列2 kbp至250 kbp范围内的所有限制性酶切片段。为平滑图谱,通过依次合并相邻的限制性酶切片段,逐步增加每个片段的大小,直至感兴趣区域内原始接触总数达到5%的阈值。为确保生物学重复间的整合及不同条件间的比较,我们在限制性酶切片段水平上同时纳入了斜率和随机截距。每种条件下的平均图谱及其之间的倍数变化如图3所示。在拟胚体分化过程中,多能性基因Pou5f1的增强子与启动子之间的相互作用减少,而中胚层-内胚层谱系决定性转录因子T的增强子-启动子相互作用则增加(图3),这些结果为理解这些发育相关增强子的功能提供了重要线索。

图1:小鼠胚胎干细胞及其衍生类胚体的代表性图像。 在无血清条件下培养的第0天小鼠胚胎干细胞(左侧)和均一的第6天类胚体(右侧),通过倒置显微镜观察。比例尺 = 500 μm。 请点击此处查看该图的放大版本。

图2:4C实验流程及本方案主要步骤的代表性图像。(A)定量4C的示意图。RS = 限制性酶切位点;US = 上游;DS = 下游;UP = 通用引物;D = RS与DS之间的距离理想情况下应为5-15 bp。(B)MboI消化后的染色质(I)、细胞核内连接的染色质(II)以及超声处理后的染色质(III)的示例。左侧数字表示通过每个样品的DNA ladder测定的DNA大小。(C)在两个位点进行PCR扩增的示例:Pou5f1(引物A和B)和T(引物C和D)。(D)使用引物A和C在Pou5f1和T位点进行多重PCR扩增的示例。ES = 胚胎干细胞;EB = 胚状体。请点击此处查看该图的放大版本。

图3:4C分析谱图示例。 在mESC和第6天EB中检测的位于Pou5f1和T基因启动子区域的诱饵所获得的定量4C谱图。上图显示了来自两个独立生物学重复实验的平均互作信号图谱;下图显示第6天EB相对于mESC的平均互作信号倍数变化(两个重复实验的平均值)。浅蓝色方框表示在分化过程中发生动态变化的增强子位置。该图改编自Tian等人的研究24。请点击此处查看该图的放大版本。
| 5 mL 用量 | |
| 1 M Tris-HCl, pH 8.0 | 50 µL |
| 5 M NaCl | 10 µL |
| 10% Igepal CA630 | 100 µL |
| 50x Roche complete protease inhibitors | 100 µL |
| MilliQ 水 | 4.74 mL |
表1:裂解缓冲液。
| 1000µL用量 | |
| MilliQ水 | 869 µL |
| 10X NEB T4 DNA连接酶缓冲液 | 120 µL |
| 20mg/mL 牛血清白蛋白 | 6 µL |
| 2000 U/µL T4 DNA连接酶 | 5 µL |
表2:连接反应主混合液的配制。
| 15 µL 体系 | |
| 5X Quick Ligation Reaction Buffer | 10 µL |
| NEBNext Adaptor | 3 µL |
| Quick T4 DNA ligase | 2 µL |
表3:接头连接反应。
| PCR 反应体系设置 | |
| 接头连接的文库(library-on-deads) | 10 µL |
| PCR 级水 | 20.25 µL |
| 10 µM 目标特异性引物 | 3.75 µL |
| 10 µM NEB Index 引物 | 3.75 µL |
| Herculase II 5X 缓冲液 | 10 µL |
| 10 mM dNTPs | 1.25 µL |
| Herculase II 聚合酶 | 1 µL |
| 总体积 | 50 µL |
| PCR 程序 | |
| 步骤 1:98 °C – 2 分钟 | |
| 步骤 2:98 °C – 20 秒 | |
| 步骤 3:65 °C – 30 秒 | |
| 步骤 4:72 °C – 45 秒 | |
| 步骤 5:返回步骤 2,共循环 15–18 次 | |
| 步骤 6:72 °C – 3 分钟 | |
| 步骤 7:4 °C – 保存 |
表4:4C染色质相互作用文库扩增,第一轮PCR。
| 巢式PCR反应体系配置 | |
| 第一次PCR所得DNA片段 | 10 µL |
| PCR级水 | 20.25 µL |
| 10 µM 特异性引物+P5 Illumina引物 | 3.75 µL |
| 10 µM P7 Illumina引物 | 3.75 µL |
| Herculase II 5X缓冲液 | 10 µL |
| 10 Mm dNTPs | 1.25 µL |
| Herculase II聚合酶 | 1 µL |
| 总体积 | 50 µL |
| 巢式PCR程序 | |
| 步骤1:98 °C - 2分钟 | |
| 步骤2:98 °C - 20秒 | |
| 步骤3:65 °C - 30秒 | |
| 步骤4:72 °C - 45秒 | |
| 步骤5:返回步骤2,共循环15-18次 | |
| 步骤6:72 °C - 3分钟 | |
| 步骤7:4 °C – 保温 |
表5:4C染色质相互作用文库扩增,巢式PCR。
| 名称 | 序列 (5'-3') |
| DS-Oct4-A | AATGATACGGCGACCACCGAGATCTACACTCTTTCCCTACACGACG CTCTTCCGATCTTCTTGCAAAGATAACTAAGCACCAGGCCAG |
| US-Oct4-A | TCTCTTGCAAAGATAACTAAGCACCAGGCC |
| DS-Oct4-B | AATGATACGGCGACCACCGAGATCTACACTCTTTCCCTACACGACG CTCTTCCGATCTGTGATGGGTCAGCAGGGCTGGAGCCGGGCT |
| US-Oct4-B | ACCAGGTGGGGGTGATGGGTCAGCAGGGCT |
| DS-T-C | AATGATACGGCGACCACCGAGATCTACACTCTTTCCCTACACGACG CTCTTCCGATCTCCTGGGTCCCTGCACATTCGCCAAAGGAGC |
| US-T-C | GATTACACCTGGGTCCCTGCACATTCGCCAA |
| DS-T-D | AATGATACGGCGACCACCGAGATCTACACTCTTTCCCTACACGACG CTCTTCCGATCTGGCTTTGGAGAGGTCAAGGAGACCCGGGAG |
| US-T-D | GCTGAGGCTTTGGAGAGGTCAAGGAGACC |
| UP-4C | CAAGCAGAAGACGGCATACGA |
| Adap-i1 | CAAGCAGAAGACGGCATACGAGATCGTGATGTGACTGGAGTTCAGA CGTGTGCTCTTCCGATC |
| Adap-i2 | CAAGCAGAAGACGGCATACGAGATACATCGGTGACTGGAGTTCAGA CGTGTGCTCTTCCGATC |
| Adap-i3 | CAAGCAGAAGACGGCATACGAGATGCCTAAGTGACTGGAGTTCAGA CGTGTGCTCTTCCGATC |
| Adap-i4 | CAAGCAGAAGACGGCATACGAGATTGGTCAGTGACTGGAGTTCAGA CGTGTGCTCTTCCGATC |
表6:用于4C文库构建的引物。
悬滴培养法无需额外添加生长因子或细胞因子,即可从预定数量的mESCs reproducibly生成均一的拟胚体(EB)群体5。本文描述了一种定量4C技术的实验方案,该方法改编自UMI-4C方法,用于在EB分化模型中定量分析谱系特异性转录因子的增强子-启动子相互作用。我们鉴定了在EB分化过程中以动态方式与Pou5f1和T基因启动子发生接触的染色质区域。Pou5f1在EB分化过程中表达下调,其启动子与远端增强子之间的接触频率也随之降低。相反,T基因在EB分化过程中表达上调,我们鉴定了三个增强子,其与启动子的接触频率在分化过程中下降(图3)。为验证上述结果,可进行活性组蛋白标记H3K27ac的染色质免疫沉淀(ChIP)实验24,因为该组蛋白修饰已被证实与增强子的激活相关,且在增强子失活过程中会丢失这一修饰标记11。
标准的4C技术已被广泛用于研究特定基因组位点的染色质互作谱25。然而,即使经过广泛的标准化处理,该方法仍难以进行定量解释26,27,28,原因在于PCR片段大小的异质性引入了偏差,且无法区分PCR重复序列。我们的定量4C方法在很大程度上与UMI-4C技术相似,后者通过超声破碎和嵌套连接介导的PCR步骤实现对单个分子的定量,从而克服经典4C方法的局限性23。然而,与使用唯一分子标识符(UMI)的UMI-4C不同,我们的定量4C方案基于超声破碎步骤产生的特异性DNA断裂来实现单分子定量。这使得我们的方案可与商业化的DNA文库构建试剂盒兼容,无需使用带有唯一分子标识符的引物。
本方案包含若干需要考虑的关键步骤。与经典的 4C 方法28类似,本方案的关键因素在于 3C 分子制备过程中消化和连接的效率。消化/连接效率低下会显著降低与目标片段相互作用的复杂性,从而导致分辨率下降。如前所述23,本方案的另一个关键步骤是文库扩增所用引物的设计。第二轮 PCR 反应的引物应位于被检测限制性酶切位点上游 5–15 nt 处。在 75 nt 的测序读长中,这可为捕获序列留下至少 40 nt 的长度用于比对。第一轮 PCR 反应所用引物应设计在第二轮引物的上游,且两者之间无重叠,同时均需具备足够的特异性以确保高效的 DNA 扩增。对于多重检测,引物应独立设计,目标熔解温度(Tm)为 60–65 °C。此外,与其他 3C 技术类似,定量 4C 方法的分辨率取决于本方案中所使用的限制性内切酶25。本方案采用识别 4 bp 序列的限制性内切酶 MboI。使用该酶所能达到的最大分辨率为约 500 bp,但此值高度依赖于基因座位置,且极少能够实现。另一项局限性在于,位于同一限制性酶切片段内的元件之间的相互作用无法被检测。此外,相距一个限制性酶切位点的相互作用无法与未消化的背景信号相区分。在连接前采用补平步骤可能有助于检测此类相互作用。
定量 4C 技术非常适合研究靶向基因座的染色质相互作用。然而,特异性 PCR 扩增步骤限制了可同时检测的基因座数量。增加靶向基因座数量的一种方法是采用多重 PCR,以同时扩增多个目标片段,但这需要所用引物之间的兼容性,并在实施前对每对引物进行测试。如果希望研究启动子区域染色质结构的全局性变化,则更适合采用全基因组范围的方法,例如 Hi-C、PC Hi-C 或 HiChIP29,30,31。
作者无任何利益冲突需要披露。
我们感谢 F. Le Dily、R. Stadhouders 以及 Graf 实验室成员提供的建议和讨论。G.S. 获得玛丽·斯克沃多夫斯卡-居里学者基金(H2020-MSCA-IF-2016,miRStem)资助,T.V.T. 获得胡安·德·拉·谢尔瓦博士后基金(MINECO,FJCI-2014-22946)资助。本研究获得欧洲研究理事会第七框架计划下的支持th 欧盟第七框架计划(ERC Synergy Grant 4D-Genome,协议号 609989,授予 T.G.),西班牙经济、工业与竞争力部(MEIC)对 EMBL 合作关系的支持,塞韦罗·奥乔亚卓越中心(2013–2017),以及加泰罗尼亚自治区政府 CERCA 计划。
| 姓名 | 公司 | 目录编号 | 评论 |
|---|---|---|---|
| 0.1% EmbryoMax 明胶 | EMD Millipore | ES-006-B | 细胞培养 |
| 0.25% 胰蛋白酶-EDTA | 25200072 | ||
| AMPure XP | Beckman Coulter | 10136224 | 4C/DNA 纯化 |
| B27 添加剂 | Gibco | 17504044 | 细胞培养 |
| β-巯基乙醇 | Gibco | 31350010 | 细胞培养 |
| Bioruptor Pico | Diagencode | B01060010 | 4C/超声处理 |
| BSA | NEB | B9000S | 4C |
| CHIR99021 | Selleck Chemicals | S1263 | 细胞培养 |
| CIP | NEB | M0212 | 4C |
| cOmplete 蛋白酶抑制剂混合物 | Roche | 4693116001 | 4C |
| DMEM/F12 培养基 | Gibco | 11320033 | 细胞培养 |
| dNTP | NEB | N0447S | 4C |
| ESGRO 白血病抑制因子 (LIF) | EMD Millipore | ESG1107 | 细胞培养 |
| 甲醛溶液 (37%) | Sigma | 252549-25ML | 4C |
| 甘氨酸 | Sigma | GE17-1323-01 | 4C |
| 糖原 | ThermoFischer | R0551 | 4C |
| Herculase II Fusion DNA 聚合酶 | Agilent | 600675 | 4C |
| IGEPAL CA-630 | Sigma | I3021-50ML | 4C |
| Knockout DMEM | 10829018 | ||
| L-谷氨酰胺 | Gibco | 25030081 | 细胞培养 |
| MboI | NEB | R0147M | 4C |
| MEM 非必需氨基酸 | Gibco | 11140050 | 细胞培养 |
| N2 添加剂 | Gibco | A1370701 | 细胞培养 |
| NEBNext DNA 文库构建试剂盒 | NEB | E6040 | 4C |
| NEBuffer 2.1 | NEB | B7202S | 4C/消化 |
| Neurobasal 培养基 | Gibco | 21103049 | 细胞培养 |
| PD0325901 | Selleck Chemicals | S1036 | 细胞培养 |
| 青霉素-链霉素 | Gibco | 15140122 | 细胞培养 |
| 蛋白酶 K | NEB | P8107S | 4C |
| Qubit 4 荧光仪 | ThermoFischer | Q33238 | 4C |
| Qubit dsDNA HS 检测试剂盒 | ThermoFischer | Q32851 | 4C |
| RNase A | ThermoFischer | EN0531 | 4C |
| 丙酮酸钠溶液 | Gibco | 11360070 | 细胞培养 |
| StemPro Accutase 细胞解离试剂 | Gibco | A1110501 | 细胞培养 |
| T4 DNA 连接酶反应缓冲液 | NEB | B0202S | 4C |
| T4 DNA 连接酶反应缓冲液 | NEB | M0202M | 4C |