Exactly Solvable Population Model with Square-Root Growth Noise and Cell-Size Regulation

技术摘要：具有平方根增长噪声与细胞大小调节机制的可解种群模型

问题陈述
本文探讨了定量细胞生理学中一个核心且尚未解决的问题：微观层面的单细胞生长速率波动如何塑造种群层面的适应度（fitness）与大小统计特性。虽然单细胞技术已经能够对分裂统计和大小控制策略进行详细表征，但单细胞生长速率的时间统计特性仍缺乏充分的约束。现有的理论模型通常假设瞬时指数增长率 $\lambda(t)$ 与细胞大小 $s$ 相互独立。然而，从微观角度来看，生物量产生的内在波动预计会随细胞大小的平方根而缩放（由于生产噪声的广延性），这导致了一个 Feller 型扩散过程。作者旨在通过引入一个包含尺寸相关（平方根）增长噪声的精确可解种群模型，来弥补这一差距，并确定其对种群增长率、稳态大小分布以及种群规模波动的影响。

研究方法
作者分析了一个尺寸结构化种群，其中个体细胞根据代表 Feller 平方根扩散的随机微分方程（SDE）进行生长：
$$\frac{ds}{dt} = k s + \sigma \sqrt{2k \bar{s}_b s} \, \xi(t)$$
其中 $k$ 是平均增长率，$\bar{s}_b$ 是平均出生大小，$\sigma$ 是无量纲噪声强度，$\xi(t)$ 是高斯白噪声。细胞分裂与大小调节采用通用建模方式：细胞在大小 $s_d$ 时分裂（$s_d$ 可能是一个依赖于出生大小 $s_b$ 的随机变量），并分裂成两个大小分别为 $r s_d$ 和 $(1-r)s_d$ 的子细胞。

分析过程分为以下几个解析步骤：

1. 种群质量动力学： 作者利用平方根噪声的可加性，证明了总种群质量 $M_t$ 遵循与单个细胞相同的 SDE，从而实现了关于总质量和种群规模的精确解。
2. 稳态大小分布： 对于恒定分裂大小的简化情况，作者求解了相关的 Fokker-Planck 方程（一个包含漂移项和扩散项的偏微分方程），以推导出细胞大小的精确稳态密度。
3. 分支过程映射： 该模型被映射为一个 Bellman-Harris 分支过程。作者通过将代际时间分布定义为 Feller 过程的一阶到达时间（first-passage time），并利用其拉普拉斯变换来计算种群规模的各阶矩。
4. 泛化： 通过将生长动力学与分裂统计学解耦，将结果扩展到包含噪声分裂大小及通用大小调节机制（如 adder、sizer、timer 模型）的情形。

核心贡献与结果

1. 精确的种群增长率：
   主要结果是，渐近种群增长率 $\Lambda$ 恰好等于单细胞平均增长率 $k$（即 $\Lambda = k$）。这一结论对于所有噪声强度 $\sigma$ 以及所有维持大小稳态的分裂/大小调节方案均成立。

   • 对比： 这与基于大小无关的增长率噪声的模型形成了鲜明对比，在后者中，$\Lambda$ 会对噪声的方差和相关时间表现出敏感性。
   • 机制： 这种“强解耦”现象专门源于加性平方根噪声的结构，这并非随机生长模型的通用属性。

2. 稳态细胞大小分布：
   作者推导出了细胞大小稳态分布的精确解析解。

   • 卷积结构： 当 $\sigma > 0$ 时，解是通过将确定性（$\sigma=0$）解与宽度为 $\sigma^2$ 的核进行单侧指数卷积得到的。
   • 矩（Moments）： 种群平均细胞大小相对于确定性情况偏移了 $-\sigma^2$，而方差仅在 $\sigma^4$ 阶发生改变。
   • 解耦： 分布的形状仅通过谱系出生大小分布受大小调节影响，而生长噪声引入了一种与分裂无关的平滑效应。

3. 种群规模波动：
   推导出了平均重标度后的种群规模 $W = \lim_{t\to\infty} N_t / \langle N_t \rangle$ 的闭式分布。

   • 分布类型： $W$ 服从复合泊松-指数分布，该分布仅取决于生长噪声参数 $\sigma$。
   • 独立性： 该分布对于分裂大小噪声、不对称分配或特定的大小调节机制具有不变性。这些因素仅影响将总质量与细胞数量联系起来的比例因子，而不影响波动的形状。
   • 方差： 长期相对方差由 $\text{Var}(N_t) / \langle N_t \rangle^2 = 2\sigma^2$ 给出。
   • 可解性： 作者指出，除了 Yule 过程（指数代际时间）之外，这是极少数已知其重标度后种群规模完整平稳分布的非平凡分支模型之一。

意义与主张
本文将带有平方根噪声的随机指数生长识别为研究受大小控制种群中生长波动的“精确可解基准”。作者声称，其结果揭示了一系列解耦层次：

1. 种群增长率独立于生长噪声和分裂细节。
2. 快照式的细胞大小分布仅通过谱系出生大小分布受大小调节影响，并伴随特定的卷积效应。
3. 种群规模波动的形状完全由生长噪声决定，且独立于分裂机制。

作者断言，这些预测与基于大小无关增长率噪声的模型存在本质区别。因此，所推导出的特征（特别是种群增长率对噪声强度的独立性以及种群规模波动的特定形式）提供了具体的、可测试的标准，用于利用现有的及未来的单细胞数据来区分不同的微观模型。这项工作并非提出新的实验方案，而是提供了一个用于解释经验数据的理论框架。