稳定、记忆与修复16 分钟

细胞怎样协调生长和分裂,维持合适的大小?

维持细胞大小,需要把每一代的偏差传回生长与周期推进过程。某些细胞能让较小的出生个体多等一段时间,再开始复制DNA;其他细胞还会调节生长速率或后续周期阶段。本篇结合2025年活体干细胞的RB家族干预与局部消融、2020年Rb稀释实验,以及2026年浓度阈值理论,用最小模型说明:抑制物被稀释本身还不够,出生时继承多少、下一轮怎样复位,同样决定能否纠正大小偏差。公开四鼠数据重绘则展示,长得更快和改变入S大小是两个不同的观测。

细胞一直在长,为什么大小没有越分越乱

同一种正在增殖的细胞,并不是每时每刻都一样大。刚出生的细胞通常较小,随后生长、复制DNA,最后分裂。因此,研究“大小是否稳定”,不能只把一张照片里最大的细胞和最小的细胞放在一起比较;要看处于相同周期阶段的细胞,尤其要看偏差如何跨世代变化。1 3

假设每个细胞都按相同的相对速率增长,到固定时间就分成两半。若这段时间刚好足够体积翻倍,那么出生时偏大的细胞,下一代仍然偏大;出生时偏小的细胞,下一代也仍然偏小。没有额外噪声时,偏差只是保留下来;一旦生长和分裂有随机变化,偏差便可能一代代累积。

所以,稳定的大小分布需要某种纠偏关系,而不仅是生长与分裂的平均数量刚好对上。 纠偏可以发生在生长速率上,也可以发生在周期推进的时机上,还可能分布在多个阶段。

哺乳动物细胞周期中的G1期,是上一轮分裂之后、复制DNA之前的阶段;进入S期意味着开始复制DNA,随后还要经过G2和M期才能完成下一次分裂。本篇首先关注G1/S这道关口,然后再问它怎样影响整个世代的大小。

环境让它先长到,还是直接让它先过关

在成年小鼠的皮肤中,研究者可以连续几天跟踪同一批表皮干细胞。2025年的一项研究提出一个具体问题:附近细胞离开或死亡后,周围细胞更容易进入周期,究竟是因为环境直接放宽了入S条件,还是因为它们长得更快,更早达到原来的条件?1

研究者用双光子激光消融局部细胞,比较距离消融位置20微米以内的邻近细胞与更远的细胞。通过三维细胞核成像和周期荧光标记,分别估计生长率,以及进入S期时的核体积。

Xie等2025年原始Fig.5,完整转载。A与E是实验和机制示意,B为真实显微时间序列,绿色圈为被消融细胞,红色圈为一个邻近细胞,标尺10微米;C、D为四只小鼠的生长与入S核大小比较。原图未裁切、改色或加注,依据CC BY-NC-ND 4.0作非商业学习展示。原作者图中“主要经生长率影响”的机制图需结合正文中的适用条件理解。
Xie等2025年原始Fig.5,完整转载。A与E是实验和机制示意,B为真实显微时间序列,绿色圈为被消融细胞,红色圈为一个邻近细胞,标尺10微米;C、D为四只小鼠的生长与入S核大小比较。原图未裁切、改色或加注,依据CC BY-NC-ND 4.0作非商业学习展示。原作者图中“主要经生长率影响”的机制图需结合正文中的适用条件理解。

图像来源:原图与说明;CC BY-NC-ND 4.0许可。

我们从公开数据表中重新提取了这两类测量。为避免把不同小鼠的体积标定混为同一把尺,以下只看每只小鼠内部,邻近组中位数除以非邻近组中位数。样本标识沿用数据表。1 5

小鼠标识 拟合生长率中位数比 入S核体积中位数比
F1 1.998 0.972
M5 1.120 1.021
Mnonclip 1.574 0.987
Mclip 1.448 1.046
本库按公开Fig.5C、5D数据重绘的同鼠中位数比。横线1表示两组中位数相等;每个点代表一只小鼠内的描述性比较。右图纵轴特意放大到0.9至1.1,不应与左图凭视觉高度比较效应大小。这不是新的显著性或等效性检验,很多细胞记录不等于很多独立动物。原始表有一组记录数与原图标注不同,未补造,详见来源记录。
本库按公开Fig.5C、5D数据重绘的同鼠中位数比。横线1表示两组中位数相等;每个点代表一只小鼠内的描述性比较。右图纵轴特意放大到0.9至1.1,不应与左图凭视觉高度比较效应大小。这不是新的显著性或等效性检验,很多细胞记录不等于很多独立动物。原始表有一组记录数与原图标注不同,未补造,详见来源记录。

四组的生长率中位数都更高,但差异强弱不一样,原论文只在四只小鼠中的三只检出显著增长率差异。相较之下,入S核体积的中位数比约落在0.97至1.05之间,原论文没有检出相应的显著大小偏移。

入S大小只来自观察窗内确实进入S期、且可明确判定事件的一部分细胞。因此,两类指标的可用细胞并不完全相同,也不能把整组生长率的变化直接当成同一批细胞逐个提前了多少小时。

这支持一条有边界的因果图:在这种局部扰动下,环境可以加快生长,使细胞更早进入与自身大小相联系的周期关口。 没检出入S大小变化,不等于证明它精确不变;这项实验也没有排除任何营养、激素或其他环境信号都可能改变周期条件。

研究还用细胞和邻域的几何特征建立预测模型,发现细胞体积是所测特征中最有力的周期状态预测量。不过,预测能力反映的是已测变量与周期的关联,不能把没有纳入模型的生化信号一并排除。对照干预比“哪一列预测最强”多提供了一层证据。

大小信号,可以藏在分子数量与体积的比例里

细胞没有必要先拿一把分子尺量出自身直径,再把数字交给周期机器。一个蛋白的数量与体积增长不同步,就可能产生大小信息。

设某种周期抑制物的分子数量为$I$,它所在的有效体积为$V$,则浓度为

$$ c=\frac{I}{V}. $$

如果$I$增长得比$V$慢,细胞越长,浓度就越低;如果数量与体积严格同比增加,浓度则不变。两者都可以看到“蛋白总量更多”,但对细胞周期的作用可能不同。

Rb及其家族成员,是这类机制的重要候选。它们能抑制E2F相关转录程序,限制细胞进入DNA复制阶段;与此同时,周期蛋白与激酶又能通过磷酸化等方式改变Rb的抑制作用。因此,数量、浓度和活性状态都要分别考虑,不能只用一个总荧光值代表全部功能。2

2020年的研究在一种端粒酶永生化的人乳腺上皮细胞系中,对内源Rb作荧光标记并追踪单细胞。G1阶段细胞核体积和总蛋白报告物增加,而Rb数量在早期数小时内相对稳定;Rb相对于总蛋白的比例随生长下降。它的行为与同时期会积累、随后被降解的p21不同。2

这一步证明的是一种分子动态:有些抑制物会被生长稀释,有些不会。 要证明这种稀释真的参与大小控制,还要改变它,再看周期与大小的关系是否改变。

研究者不仅删除Rb、调高Rb表达,也用持续表达的外源Rb使其浓度不再随细胞大小正常下降。后者可以延长平均G1时长,却没有恢复G1时长对出生大小的依赖。这比单纯看到“Rb多,细胞变大”更接近区分性证据:某个因子能拖慢周期,与它能把大小偏差传入周期,仍是两件事。

这些关键培养实验包含低剂量CDK4/6抑制条件,用来保留对Rb作用敏感的周期状态。结论应连同该条件理解,不能写成所有未经干预的培养细胞都以同样强度依赖Rb稀释。原文也把稀释与磷酸化看作可以共同起作用的途径。2

活体遗传实验支持了哪一步

2025年的皮肤研究先单独删除了Rb1,发现大小控制变化不大。随后在缺少家族成员Rbl1的背景下,再删除Rb1。双缺失后,出生时较小的细胞不再像对照那样明显延长G1,出生核大小与G1时长的相关也明显减弱。1

这里要把基因组合说清楚:对照是Rbl1单缺失,实验组是Rb1与Rbl1双缺失,不能缩写成“删除Rb1就破坏了活体大小控制”。Rbl1编码的p107可以补偿Rb1的部分功能;在培养的鼠细胞中,降低Rb1后,Rbl1表达上升,而且其浓度仍与大小有关。

这提供了RB家族参与控制的因果线索,但还有两个重要边界。

第一,这一组局部双缺失实验是在一只小鼠的两侧耳朵内作处理与对照。57个对照细胞和49个双缺失细胞,并不等于57只和49只独立动物。它有同体配对的优点,也有动物层面重复有限的限制。较晚时双缺失组织还出现细胞死亡和DNA损伤,不能把这些后果都归因为大小控制这一项功能。

第二,双缺失组的G1很短,而成像间隔为12小时。研究者认为该采样足以比较G1时长与出生大小的联系,却不足以精确估计该组入S时的大小。更直接的Rb/Rbl1浓度随大小变化,主要是在培养的MEF和NIH3T3细胞中测得,尚未在这些活体皮肤细胞内把分子浓度轨迹与入S事件完整同步测清。作者自己也保留了这个缺口。

因此,当前较有依据的是“RB家族通路有助于把大小与周期推进耦合起来”。“在活体中,具体是哪一项浓度变化、多大变化足以改变入S概率”,仍需进一步测量。

稀释本身还不够,出生时带来多少也很重要

用一个小模型把这个逻辑走通。它保留G1生长和抑制物浓度阈值,假定调控分子可及体积与所跟踪的大小保持固定比例,但不模拟完整的Rb反应网络。

设G1中体积按$V(t)=V_b e^{gt}$增长,$V_b$为出生大小,$g$为相对生长率。G1中抑制物数量$I_b$保持不变,浓度降到$c_*$时允许进入S期,于是

$$ V_{G1/S}=\frac{I_b}{c_*},\qquad t_{G1}=\frac1g\ln\frac{I_b}{c_*V_b}. $$

如果不同出生大小的细胞继承近似相同的$I_b$,它们会在相近大小进入S期。小细胞最初浓度更高,需要长得更多、等待更久;大细胞则较早达到条件。

但如果$I_b$刚好与$V_b$成正比,所有细胞的出生浓度就相同,也只需要相同的增长倍数来稀释到阈值。稀释确实发生了,却没有校正出生大小的差别。 因而不能只测到“生长时浓度下降”,就断言大小稳态已经得到解释。

本库把两种情形以及一个中间情形写成

$$ I_b=c_*V_c\left(\frac{V_b}{V_*}\right)^\beta. $$

$V_*$是无噪声下选定的目标出生大小;$\beta=0$表示继承数量与出生大小无关,$\beta=1$表示同比增加,$\beta=0.5$表示部分依赖。$V_c$是控制阈值尺度的常数。这些是示意设定,不是对Rb实验拟合出的参数。

本库抑制物稀释模型。左图显示出生大小与入S大小的关系:相同的浓度阈值,配上不同的抑制物继承规律,会产生不同大小控制。右图显示所需G1时长;绿色虚线只把G1生长率翻倍,入S体积不变,所需时间减半。坐标均已归一化,不对应实测小时数或细胞体积。
本库抑制物稀释模型。左图显示出生大小与入S大小的关系:相同的浓度阈值,配上不同的抑制物继承规律,会产生不同大小控制。右图显示所需G1时长;绿色虚线只把G1生长率翻倍,入S体积不变,所需时间减半。坐标均已归一化,不对应实测小时数或细胞体积。

2020年研究为继承规律提供了实验线索:不同大小的细胞在S/G2/M中合成的Rb量,随体积变化得较弱;Rb还比一个自由核定位荧光蛋白更均等地分配给两个大小可能不同的子细胞。Rb与染色体的关联,是作者提出的一个解释。它使“出生时的分子份额不简单按体积分配”成为可测的机制环节,而不是一句抽象的大小感知。2

单次过关怎样变成代际纠偏

进一步假定,G1以后体积再增大固定倍数$G=1.2$,然后对半分裂。选择$V_c=2V_*/G$后,下一代出生体积满足

$$ \ln\frac{V_{n+1}}{V_*} =\beta\ln\frac{V_n}{V_*}+\eta_n. $$

$\eta_n$表示每代继承数量的小幅随机波动。这个式子把反馈说得很直接:上一代的对数大小偏差,只以$\beta$倍传给下一代,再加上新的波动。

若$\beta=0$,无噪声时下一代就回到目标;若$\beta=0.5$,偏差逐代减弱;若$\beta=1$,旧偏差完整留下,新的噪声又继续累积。这里的$\beta=0.5$是对数偏差的部分纠正,并不是全范围严格相同增量的adder。

左图从目标大小的0.7倍和1.3倍出发,比较无噪声下的后代;实线与虚线分别表示大、小初态。右图为每组20000条模型谱系、120代的对数大小标准差,圆点为解析结果。具有纠偏的两组保持有限宽度,β=1组不断变宽。这是模型世代,不是小鼠的实际谱系追踪。
左图从目标大小的0.7倍和1.3倍出发,比较无噪声下的后代;实线与虚线分别表示大、小初态。右图为每组20000条模型谱系、120代的对数大小标准差,圆点为解析结果。具有纠偏的两组保持有限宽度,β=1组不断变宽。这是模型世代,不是小鼠的实际谱系追踪。

本例给每代加入独立、均匀分布在$-0.12$至$0.12$之间的对数波动,其方差为0.0048。120代后,三组对数大小方差分别约0.00481、0.00636和0.580,对应解析值0.0048、0.0064和0.576。无反馈那组的方差按世代数累积,而不是仅仅换到另一个稳定目标。

代码检查了所有模拟细胞的G1增长时间均为正,没有通过强制缩小细胞或裁剪负值来制造稳态。脚本位于research/calculations/cell_size_control.py,另附实验数据提取表和计算记录。

这个模型仍有明确的条件:它给定了每代$I_b$的规律,固定了后半周期增长倍数,并假定对称分裂。它没有从蛋白合成、降解和分配反应中自行生成全部规律。如果环境同时改变后半周期或分裂比例,目标出生大小也可能改变。仅凭入S大小相似,不能保证整个周期的所有大小统计都不变。

下一轮怎样重新开始,是机制的一部分

2026年的一项浓度控制理论强调了一个基础问题:如果分子数量和体积在分裂时都按同一比例分配,浓度就不会因为分裂自动清零。一个“达到浓度阈值就立刻触发”的模型,必须解释下一代为什么不会停在同一触发状态,以及如何准备下一次事件。4

该研究构造了多阶段方案:不同阶段切换分子的生产和降解,某一阶段完成后,系统进入另一套动态条件。这样,浓度阈值可以与反应状态的复位共同组成连续运行的周期。这是明确模型中的一种实现,不意味着所有真实细胞都必须遵循同一组方程。

这也解释了为什么前面要讲Rb在不同阶段的合成与继承,而不能只留一句“长大就把抑制物稀释了”。对于本库的小模型,出生数量规律是一项输入条件;在真实细胞中,它必须由上一轮的分子生产、降解、分配和阶段切换产生。单次过关与反复运行,是两个都需要讲通的问题。

不能把一条斜线当作完整机制

大小控制研究常画“出生大小”与“之后新增多少”的关系。但新增量本身就是末端大小减去出生大小,两个坐标共享同一次出生测量。

本库做了一个纯统计反例:设所有模拟细胞真实都增加相同体积,出生体积与新增量本来无关;给出生和末端测量加入独立噪声。由于出生误差在横轴中被加进去、在新增量中又被减掉,就会制造负相关。本例得到斜率约$-0.503$,与按误差方差推导的$-0.5$一致。这不是对论文误差大小的估计,也不否定它的结果;它说明单看负斜率不能辨认具体机制。

2025年研究也专门检查了时间采样变粗造成的偏差。12小时采样对于较短的双缺失周期并不理想,作者因此优先比较出生大小与G1时长,而没有宣称准确恢复了双缺失细胞的G1新增体积。1

测量的“大小”也需要限定。活体研究的大部分分析用细胞核体积作为全细胞体积的代理,只有某些完整几何分析直接分割整个细胞。核体积、全细胞体积、干质量和含水量不能互相直接替代;该研究也无法根据预测表现明确区分,实际被周期机制读取的究竟是核体积还是全细胞体积。

因此,更可靠的路线是把连续成长轨迹、分子数量与活性、遗传或环境干预,以及后代是否纠偏放在一起看。一个平均大小变了,既可能是目标范围改变,也可能是纠偏能力变弱,两者需要不同证据。

G1的大小依赖较弱,也可能在别处纠偏

2018年的另一组研究采用荧光排除法:细胞排开不能进入它的荧光探针,排开的空间可以用于估计细胞体积。它直接跟踪多个培养细胞类型的完整增长轨迹,发现不少对象的全周期表现接近“每轮增加近似固定的体积”,称为adder或近adder。3

固定增量并不要求偏大的细胞绝对比偏小的少长。在理想对半分裂模型中,若每轮都增加$\Delta$,则$V_{n+1}=(V_n+\Delta)/2$,偏差就会逐代缩小。较大细胞虽然也加了同样多的体积,相对自身的增幅却较小。

但相似的adder外观可以由不同过程拼出来。该研究发现,不同细胞对G1时长、S/G2时长和相对生长率的调节并不相同。把HeLa细胞先暂时阻滞、使其变大,再释放后,较大细胞的G1时长已接近下限,仍可通过与大小有关的生长率变化表现出纠偏。分子层面究竟怎样调节增长率,该研究没有完整解决。

这组培养结果与活体皮肤结果不必互相取代。测量方法、细胞类型、周期速度与培养环境都不同。它们共同说明,G1/S是一个重要控制点,但不能由这里的强弱直接判定整个细胞是否具有大小稳态。

当前能讲清楚的答案

细胞大小能够稳定,是因为生长、周期推进和分裂之间存在使偏差减弱的耦合。对某些细胞,生长改变抑制物与其他调控组分的浓度关系,使较小的细胞多等一段时间;对另一些条件,后续阶段或生长率调节也会参与。不同路径可以得到相似的统计外观。

2025年的活体研究支持RB家族和细胞自身大小与G1/S相耦合,也通过局部消融说明外界可以主要改变“长到那里的速度”。但具体的活体分子浓度轨迹、动物层面的重复,以及其他环境干预能否移动目标范围,仍未全部解决。

最值得继续追的是:把同一个活体细胞的核与胞质大小、RB家族浓度及活性、周期事件同时跟踪,再通过保留蛋白但改变其随大小变化的方式作干预。这样才能把“这个通路很重要”继续推进到“这个通路如何读取大小、怎样在每轮重新准备,并把偏差逐代压小”。

来源与依据

  1. S1 原始活体追踪、类器官、RB家族与局部细胞消融研究

    The G1/S transition in mammalian stem cells in vivo is autonomously regulated by cell size

    Shicong Xie, Shuyuan Zhang, Gustavo de Medeiros, Prisca Liberali, Jan M. Skotheim · 2025

    Nature Communications16,9071;2025-10-13发表版PDF/HTML,SI和Source Data · 相关全文已核

    Results各节、Fig3与Fig5;Methods分割/跟踪、核体积代理、CreER诱导、统计建模、消融;SI S7/S9/S10/S20说明;Source Data Figure5C和5D逐条提取;原始Fig5已查看

    除Fig4外主要用核体积,12h采样不足以准确恢复DKO入S大小,改用时长关系。Rb1单敲观察在诱导4-5周后;急性耳部双敲是在1只3月龄雄鼠两耳比较Rbl1单缺失与Rb1/Rbl1双缺失,不能按57/49细胞当独立动物。Fig3E文字把flox基因误写Rbl1,Methods和正文为Rb1-flox,本文按实际基因型记。双敲约4.5日后死亡/损伤不归因于大小一项功能;活体浓度稀释未直接测清,分子浓度证据在MEF/NIH3T3培养。消融有4鼠749生长记录、273个可用入S体积记录;F1非邻近组表中23、原图标27,未补齐。图原许可CC BY-NC-ND4.0,完整未改用于非商业学习;数据重绘为描述性计算。

  2. S2 原始内源标记、蛋白缩放/分配和大小耦合干预

    Cell growth dilutes the cell cycle inhibitor Rb to trigger cell division

    Evgeny Zatulovskiy, Shuyuan Zhang, Daniel F. Berenson, Benjamin R. Topacio, Jan M. Skotheim · 2020

    Science369,466-471;PMC7489475作者公开全文 · 相关全文已核

    主文Fig1-4说明与相关结果:HMEC-hTERT1内源标记、G1数量/体积、mitotic amount和姐妹分配、Rb缺失及恒定表达;Fig4说明的CDK4/6抑制条件

    内源Rb早期G1数量近恒,核体积/总蛋白增长;Rb与p21动态不同。强大小耦合实验G-L含100nM palbociclib、F45nM,不能冒充所有未干预细胞。恒定表达延长G1但不恢复大小依赖为区分性证据;未将29h半寿命或三倍入S率等未核SI全部方法的精确量作为本篇承重数字。SI三条通路均未得到可用PDF,只核正文与图说明,不称全部补充已读。

  3. S3 原始全周期体积追踪、控制扰动与大小稳态框架

    Size control in mammalian cells involves modulation of both growth rate and cell cycle duration

    Clotilde Cadart, Sylvain Monnier, Jacopo Grilli, Pablo J. Saez, Nishit Srivastava, Rafaele Attia, Emmanuel Terriac, Buzz Baum, Marco Cosentino-Lagomarsino, Matthieu Piel · 2018

    Nature Communications发表版;PMC6095894全文 · 相关全文已核

    Results FXm、near-adder、G1/S-G2及大HeLa生长率;Discussion;Methods FXm、培养控制和关键测量时点

    FXm测排开探针的体积,部分HeLa有轻度探针摄入;birth定义为分裂后40min、mitosis为分裂前60min,避开膨胀,不按瞬时严格体积守恒理解数据。部分初代细胞转装置后非稳态,未纳入最终统一框架。近adder是表型,不自动唯一分子机制;Roscovitine释放的大HeLa例子说明弱G1时长调节不等于全部控制消失。未取得其单独SI,也未重拟合全库增长轨迹。

  4. S4 原始浓度阈值和多阶段周期理论

    General mechanism for concentration-based cell size control

    Motasem ElGamel, Lucas Ribaudo, Andrew Mugler · 2026

    Physical Review Research8,033012;2026-07-06期刊发表版,APS harvest接口 · 相关全文已核

    §I-III:分裂不改变共同分配的浓度、局部跨代映射、给定多阶段生产/降解方程;§IV开头统计与控制类型的区分

    本文只使用阈值需要说明下一轮准备/复位以及作者多阶段构造,未独立复现其全参数扫描或酵母比较。未援引c必须在所有时间全局等于F(V)的广义必要性证明,也不宣称多阶段是所有浓度模型唯一可能实现。本库beta继承模型另行定义,不冒充本论文动力学复现。

  5. S5 原研究公开分析代码

    Image analysis and statistical modeling tools for autonomous cell size control

    Shicong Xie / Skotheim lab · 2025

    master提交30d662597259a3752c682db068d8c92c7ef8dd85固定快照 · 相关全文已核

    semiauto_tracking_segmentation/5__load_ablation_regions.py的按细胞/地区汇总和组别;tissue_dynamics_model/3__classify_g1s.py的keep_first_S及train_test_split;README

    消融分析用每细胞汇总df,而非把每时点当细胞;同鼠跨地区/实验仍相关。代码对Mclip和Mnonclip入S体积共同乘0.7,未核到该系数的充分标定解释,故本篇只比同鼠两组比值,不作跨鼠绝对标定推断。分类器保留多个G1帧后按行随机划分且去掉cellID,不能据其AUC自称按新细胞/新动物外部验证;本篇不重报预测性能作为因果证明。代码仅查看,未执行其私有路径或pickle读取。

最后更新 2026-09-09

来源