图 2. DNMT3A 编辑细胞中的基因表达动态表明了一种不同于二进制的记忆形式。A 使用与 dCas9、PhlF 或 rTetR 融合的 KRAB、DNMT3A 或 TET1 作为 DNA 结合域 (DBD) 进行瞬时表观遗传编辑的概述。B 本研究开发的实验系统示意图。报告基因通过位点特异性染色体整合整合到内源性哺乳动物基因座中。哺乳动物组成型启动子 (EF1a) 驱动荧光蛋白 EBFP2 的表达。上游结合位点能够靶向募集表观遗传效应物,这些效应物与 DNA 结合蛋白 rTetR、PhlF 或 dCas9 融合。报告基因两侧是染色质绝缘体,以与其他基因隔离。 C 实验概述描述了瞬时转染到带有报告基因的细胞、基于转染水平的荧光激活细胞分选和时间过程流式细胞术测量。D 根据图 C 中显示的实验时间线,DNMT3A 编辑(DNMT3A-dCas9)报告基因的基因表达动态。显示的是 DNMT3A 编辑细胞的单细胞流式细胞术测量(EBFP2)。DNMT3A-dCas9 靶向启动子上游的 5 个靶位点,并使用乱序 gRNA 靶序列作为对照(图 SE.2 A、B、表 S3)。黄色阴影表示检测到转染标记的时间。显示的数据来自 3 个独立重复的代表性重复。E 转染 DNMT3A-dCas9 和细胞分选后 14 天进行 MeDIP-qPCR 和 ChIP-qPCR 分析,以获得高水平的转染。分析了启动子区域(表 S4 和方法)。显示的数据来自三个独立的重复。报告的是使用标准 ∆∆ C t 方法相对于活性状态的倍数变化及其平均值。误差线是平均值的标准差。DNMT3A-dCas9 靶向启动子 (gRNA) 上游的 5 个靶位点。使用乱序的 gRNA 靶序列 (gRNA NT) 作为对照。* P <0.05,** P <0.01,*** P <0.001,非配对双尾 t 检验。F 根据图 C 中显示的实验时间线的 KRAB 编辑 (PhlF-KRAB) 基因表达动态。显示的是单个细胞的报告基因 (EBFP2) 的流式细胞术测量值。黄色阴影区域表示在未应用 DAPG 期间检测到转染标记的时间。从第 6 天开始,在 PhlF-KRAB 和 PhlF 条件下应用 DAPG。每天测量不同的独立重复。显示的数据来自 3 个独立重复。G 转染 PhlF-KRAB 和高水平转染细胞分选后 6 天的 MeDIP-qPCR 和 ChIP-qPCR 分析。分析的是启动子区域。数据来自三个独立重复。显示的是相对于活性状态的标准 ∆∆ C t 方法确定的倍数变化及其平均值。误差线是平均值的标准差。* P <0.05,** P <0.01,*** P <0.001,非配对双尾 t 检验。H 当 KRAB = 0、TET1 = 0 时获得的染色质修饰回路。参见 SI 图 SM.1 C。I 上图:(CpGme, X) 对的剂量反应曲线。下图:DNMT3A 脉冲强度与 DNA 甲基化等级 (CpGme) 之间的剂量反应曲线。脉冲强度通过增加其高度来增加。参见 SI 图 SM.1 D 和 SM.3。J 系统基因表达的平稳概率分布,由 SI 表 SM.1 和 SM.4 中列出的反应表示,参数值在 SI 第 S.9.3 节中给出。K 系统在 t = 28 天后的基因表达概率分布,如图 J 所示,参数值和初始条件在 SI 第 S.9.4 节中给出。参见 SI 图 SM.1 B 和 SM.2。在图 I 和 J 中,DNMT3A 动力学被建模为随时间呈指数下降的脉冲(参见第 S.1.1 节 - SI 方程 (SM.7))。在我们的模型中,ε (ζ) 是衡量基础(招募)擦除率与每次修饰的自催化率之间比率的参数。参见 SI 图 SM.1 E 和 SM.3。
图2。DNMT3A募集后的基因表达动力学与数字记忆不一致。使用特定于特定于染色体的染色体整合的169个报告基因基因的示意图。哺乳动物170构成启动子(EF1A)驱动荧光蛋白EBFP2的表达。上游结合位点可实现靶向171的表观遗传效应子,该效应子与DNA结合蛋白RTETR融合在一起,PHLF或DCAS9。报告基因是由染色质绝缘子与其他基因分离出来的172。b实验概述,描述了瞬时转染到具有报告基因的173个细胞,基于转染水平的荧光激活的细胞分选,以及时间顺序的流量细胞仪174测量。根据面板中所示的175个实验时间表。显示的是四种不同水平的转染水平的报告基因176(EBFP2)的流量细胞仪测量值的分布。DNMT3A-DCAS9靶向启动子上游的5个目标位点,177用作炒GRNA目标序列作为对照(图se.2 a,b,表S3)。显示的数据来自来自3个独立重复的代表性178重复。d使用DNMT3A-179的流量细胞仪的单细胞基因表达测量值对应于面板C中所示的细胞(30天)。父母是指带有180个报告基因的未转染细胞。数据来自3个独立重复的代表性重复。平均值。e MedIP-QPCR和ChIP-QPCR 181分析DNMT3A-DCAS9和细胞分类后14天分析高水平的转染。分析了启动子区域182。显示的数据来自三个独立的重复。报道的是折叠变化及其平均值,使用183标准∆ ∆ c t方法相对于活性状态。错误条为S.D.DNMT3A-DCAS9的靶向位置为184至5个目标位点(GRNA)。使用炒GRNA目标序列(GRNA NT)作为对照。185 *p≤0.05,**p≤0.01,***p≤0.001,未配对的两尾t检验。根据面板中所示的实验时间线,krab抑制的基因表达动力学(PHLF-KRAB)186。所示是从四种不同水平的转染水平的187个报告基因基因(EBFP2)的流量细胞仪测量值的分布。每天测量一个独立的重复。显示的数据188来自3个独立重复。g重新激活细胞的百分比(400-10 5基因表达A.U.F.)对应于F. h Medip-QPCR面板中显示的189个细胞种群和CHIP-QPCR分析后6天对PHLF-KRAB和Cell 190排序进行了高水平的转染。分析是启动子区域的。数据来自三个独立的重复。191显示的是折叠变化,其平均值由标准∆ΔCT方法确定相对于活性状态。错误192条是S.D.平均值。p≤0.05,**p≤0.01,***p≤0.001,未配对的两尾t检验。参见SI图参见Si无花果。202i简化染色质修饰193当krab = 0,dnmt3a = 0,tet1 = 0时获得的电路图,而H3K9me3并未介导从头催化194 DNA甲基化的催化。SM.1 C. J顶图:(CPGME,H3K4ME3)对的剂量响应曲线。底部图:(DNMT3A,CPGME)对的剂量-195响应曲线。SM.1 D和SM.3。 k k的基因表达的概率分布196的系统,该系统由Si Tape Sm.1和Sm.3中列出的反应表示。 参见Si无花果。 SM.1 B和SM.2。 l概率197在t = 28天后的基因表达分布,如面板I所述获得。 在小组j和l中,将198 DNMT3A动力学建模为脉冲,随着时间的流逝会呈指数减小。 在我们的模型中,α'是通过抑制组蛋白修饰的DNA甲基化建立的归一化速率199,DNA甲基化擦除率200速率与激活组蛋白的擦除速率和激活的组蛋白修改速率之间的µ'是每个基准级别(ε')的级别(均为基础率(均))(招募)(招募)(招募)。修改。 参见SI图 SM.1 E和SM.3。SM.1 D和SM.3。k k的基因表达的概率分布196的系统,该系统由Si Tape Sm.1和Sm.3中列出的反应表示。参见Si无花果。SM.1 B和SM.2。 l概率197在t = 28天后的基因表达分布,如面板I所述获得。SM.1 B和SM.2。l概率197在t = 28天后的基因表达分布,如面板I所述获得。在小组j和l中,将198 DNMT3A动力学建模为脉冲,随着时间的流逝会呈指数减小。在我们的模型中,α'是通过抑制组蛋白修饰的DNA甲基化建立的归一化速率199,DNA甲基化擦除率200速率与激活组蛋白的擦除速率和激活的组蛋白修改速率之间的µ'是每个基准级别(ε')的级别(均为基础率(均))(招募)(招募)(招募)。修改。参见SI图SM.1 E和SM.3。SM.1 E和SM.3。