《Nature Genetics》：Enhancer control of promoter activity and variability via frequency modulation of clustered transcriptional bursts
编辑推荐：
在哺乳动物细胞中，基因表达由增强子控制，这些增强子通常分散在启动子周围的大型顺式调控区域（cis-regulatory landscape）中。增强子如何决定其靶基因的转录，以及这如何取决于它们在顺式调控区域内的相对位置，仍不清楚。在这里，研究人员使用活细胞成
在哺乳动物细胞中，基因表达由增强子控制，这些增强子通常分散在启动子周围的大型顺式调控区域（cis-regulatory landscape）中。增强子如何决定其靶基因的转录，以及这如何取决于它们在顺式调控区域内的相对位置，仍不清楚。在这里，研究人员使用活细胞成像（live-cell imaging）来追踪同一个增强子控制下的启动子活性，但该增强子插入在具有最小复杂性的简化调控区域内的不同位置。结合数学建模，这揭示了启动子的RNA产生发生在转录爆发（transcriptional bursts）的簇中，增强子位置控制簇出现的频率。这导致当增强子在基因组上靠近启动子时，爆发比在大的基因组距离时更频繁且更均匀地在细胞中发生。数学建模进一步表明，增强子调节启动子从其基础转录状态转变为爆发簇变得更频繁的体系（regime）的能力。研究人员的结果揭示了哺乳动物启动子运作尚未探索的模式，并表明增强子在顺式调控区域内的位置关键地控制单细胞中转录输出的时间性和变异性。
在真核生物中，增强子通过调控转录爆发动力学来控制基因表达，但其在顺式调控区域（cis-regulatory landscape）中基因组位置如何影响启动子活性和细胞间变异性，长期未被直接解析。此前基于群体平均的染色体构象捕获和RNA定量研究暗示增强子通过三维（3D）染色质折叠实现物理邻近来传递调控信息，但这些偶遇与RNA合成在单细胞中的时空关系仍不明确，尤其在哺乳动物细胞中。单细胞内的增强子和启动子事件本质上是概率性的，转录因子和辅因子的随机结合/解离导致RNA聚合酶II（Pol II）复合物的不连续装载和许可，产生间断的转录爆发。现有的简单二状态（On/Off）或三状态模型难以解释增强子如何调节爆发频率及细胞间变异性。为隔离增强子-启动子基因组距离的影响并排除其他调控变量（如序列同一性、生化兼容性），研究人员在小鼠胚胎干细胞（mESCs）中采用简化的中性拓扑关联结构域（TAD）系统，将Sox2启动子及其同源Sox2控制区（SCR）增强子以不同距离插入，通过活细胞成像和数学建模系统研究转录爆发动力学。论文发表在《Nature Genetics》。
主要关键技术方法（≤250字）：研究人员基于piggyBac转座子系统，在mouse mESCs的中性TAD（位于染色体15，无其他活性增强子、启动子或CTCF环）中整合了由Sox2启动子驱动、含24个MS2茎环重复序列（标记于EGFP转录本的3'UTR）的报告基因，并通过PBase转座酶介导SCR增强子随机整合至距启动子不同基因组距离（5–232 kb及另一染色体作为对照），生成14个克隆系。稳定表达HaloTag标记的MS2外壳蛋白（MCP-HaloTag）以结合MS2茎环进行活细胞成像。采用倾斜照明显微镜以30秒间隔采集5小时z-stack，并结合Capture-3C测量增强子-启动子接触概率。通过数学建模（双体系二状态模型）分析爆发参数，利用Kaplan–Meier估计器校正删失数据。
研究结果：
Enhancer–promoter genomic distance affects interburst duration：通过分析5个克隆系的活细胞成像数据（去除首小时以避免荧光衰退），发现爆发大小（integrated MS2 intensity）和幅度（peak intensity）的生存分布在各个距离间无显著差异（95%置信区间内重叠），增强子距离仅轻微影响幅度。爆发持续时间（burst duration）的中位生存时间约1.5分钟，且与距离无关。相反，爆发间期间隔（interburst duration）随基因组距离增加而延长：中位值从近距离（5–30 kb）的约15分钟升至232 kb处的约30分钟，间隔超过4小时的概率从5 kb处的2%升至232 kb处的30%。这表明增强子-启动子距离主要调控爆发频率，而非爆发时长或大小。
Burst frequency and its variability across cells scale nonlinearly with contact probability：平均爆发频率（每小时爆发次数）与之前报道的转录输出和接触概率之间的非线性关系一致。近距离（5 kb）时平均约7次/4小时，远距离（149 kb）时约1次/4小时。单细胞爆发频率的变异系数（coefficient of variation）随着接触概率增加而降低，说明增强子在基因组上更靠近启动子时，细胞间爆发的时序变异性更小。
Transcriptional bursts occur in clusters：爆发间期生存分布并非单指数形式，而是最佳拟合为三个指数成分的加权和，时间尺度分别为3–7分钟、24–52分钟和>164分钟。连续两个爆发间期呈正相关（短间隔后跟随短间隔，长间隔后跟随长间隔），表明启动子并非运行在单一爆发频率下，而是在至少两个体系（regime）之间切换：一个高频体系（短间隔连续出现）和一个低频体系（长间隔连续出现）。爆发常以簇形式出现，即多个连续爆发被较长不活跃期分隔。
Modeling the multistate kinetics of promoter operation：为解释上述现象，研究人员建立了双体系二状态模型：每个体系含有一个On状态（RNA合成）和一个Off状态（无转录），但高频体系的Off→On转换速率（On rate）远高于低频体系。通过同时拟合单个克隆（5 kb距离）的爆发时长、间隔时长、稳态爆发概率以及连续间隔的相关性，模型成功重现了所有四项观测量。最佳拟合参数给出高频体系平均间隔21分钟、低频体系平均间隔4小时，高频体系平均持续4.7小时，约40%细胞处于高频体系。对于其他距离（29、149、232 kb），仅允许低频到高频的转换速率（low-to-high transition rate）变化，而固定其他参数，模型同样能够复现间隔生存分布和连续间隔的相关性趋势。该转换速率随接触概率（因而随距离减小）呈非线性增加。结论为增强子通过调节启动子从低频到高频体系的转换速率来控制爆发簇出现的频率，而每个簇内的爆发频率（即高频体系内的On rate）不受距离影响。
讨论总结：研究人员通过仅改变增强子在无结构复杂性中性区域内的位置，隔离出基因组距离对转录动力学的直接影响。结果证实增强子主要调控爆发频率而不改变爆发时长，并首次揭示爆发成簇出现和连续间隔正相关，这与简单二状态或三状态模型不相容。双体系模型很好地解释了这些现象，且增强子仅调控低频到高频的转换速率，暗示增强子-启动子相互作用可能通过暂时提高启动子附近转录因子或辅因子的可用性来调节转换。研究人员估算从低频到高频的转换平均每7小时（5 kb）至151小时（232 kb）发生一次，远低于增强子-启动子物理接触的频率（可能每几十秒一次），提示大多数接触可能无效，仅少数达到最小持续时间或由环挤压(cohesin)驱动的稀有但长久的接触才能触发体系转换。此外，近距离增强子产生更频繁、更均匀的爆发，而远距离增加细胞间变异性，表明顺式调节区域的大规模组织影响转录程序的精确性和稳健性。
研究结论（翻译）：总之，研究人员的数据揭示了远端增强子控制单细胞转录及其变异性的定量原理，并为未来研究从机制上解决增强子-启动子相互作用与转录爆发之间的相互作用提供了定量框架。
