只有边缘留下后代
琼脂上的菌落向外长,是因为边缘的细胞向空地分裂。靠里的细胞被邻居包围,营养也跟不上,就不再有贡献。在 Hallatschek 等人的大肠杆菌菌落里,这条生长带只有约七十微米宽,只占整个种群的一小部分。
明天前沿上的每个细胞,都是今天前沿上某个细胞的后代。前沿就像一个每一步都被重新抽样的小种群,所以即使菌落里有几十亿个细胞,那里的漂变依然很强。相邻细胞很快就有了共同祖先,前沿碎成一段段单色区域。在前沿后面,这些区域被固定下来,成为扇区。
Lab 在方格上生长一个 Eden 型菌落。每个旁边有空格的占据格以速率 1 分裂,子细胞随机填入一个空的邻格,并继承颜色。本章里没有哪个细胞比别的更适应。
来源:O. Hallatschek · O. Hallatschek and D. R. Nelson
边界比随机游走跑得快
扇区两侧的边界在前沿上相遇,扇区就消失了。如果每条边界都绕一条直线做普通随机游走,它侧向位移的均方 ⟨y2⟩ 会和生长距离 L 成正比。Hallatschek 等人分析了大肠杆菌菌落中的 92 条边界,发现它们游走得更快:⟨y2⟩ 按 L2ζ 增长,ζ=0.65±0.05,高于随机游走的 1/2。
他们的解释是:边界在局部总是垂直于菌落边缘,而边缘是粗糙的。宽度为 λ 的鼓包高度约为 λχ,让边缘倾斜 λχ−1,并能维持约 λz 的时间。倾斜把边界往侧面推。在时间 t 内起作用的鼓包宽约 t1/z,所以边界以 t(χ−1)/z 量级的速度侧向漂移,位移按 t1+(χ−1)/z 增长。
Eden 模型的粗糙度指数 χ=1/2,动力学指数 z=3/2,所以 ζ=2/3。Lab 用 80 条宽 400 的固定种子条带测得 ζ=0.656±0.045,与预测和实验都相符。不过这只是对格点模型的检验,不是对细菌的新测量。
扇区边界在距离 L 上的游走,以及由 Eden 模型前沿粗糙度预测的指数。来源:O. Hallatschek
完整例题:数直线前沿上的扇区
题目。一条宽 1000 格的直线前沿,开始时随机放两种颜色。长了 32 行之后,平均有 43.2 条扇区边界。长到 400 行时还剩多少?
标度。如果边界游走的距离量级是 Lζ,那么比这更窄的扇区都已有机会被挤掉,存活扇区的典型宽度也就是这么宽。前沿宽度固定时,边界数按 L−ζ 减少。
预测。43.2×(400/32)−2/3=43.2/5.39=8.0。如果是指数为 1/2 的普通随机游走,会剩下 43.2×12.5−1/2=12.2 条。
测量。40 条固定种子条带的平均值是 8.0±0.3 条边界。在 32 到 400 行之间,对均值取对数、对距离取对数后做直线拟合,斜率为 −0.669。每条条带单独拟合,得 −0.682±0.020,离 −1/2 约有 9 个标准误。
结论。大约八条边界,也就是周期条带上大约八个扇区。下图给出完整曲线。
直线前沿生长距离 L 之后,扇区的典型宽度与边界数。宽度为 1000 的双色 Eden 条带上,扇区边界的平均数目随生长距离的变化,取 40 个固定种子菌落的平均(两轴都取对数)。实线参考线按距离的 −2/3 次幂下降,虚线按 −1/2 次幂下降,后者是普通随机游走给出的结果。来源:O. Hallatschek · K. S. Korolev
圆形菌落会留下一些扇区
直线前沿上的粗化永不停止:长得足够远,一种颜色就会占满整个前沿。圆形菌落不同。它的周长与半径成正比地增长,而边界只按 Lζ 游走,ζ<1。过了某个交叉点,膨胀占上风,相邻边界分开的速度超过它们游走靠拢的速度,存活下来的扇区就再也不会消失。
让两种效应相平衡,就得到右边的标度关系,Ns 是存活扇区数,r0 是起始液滴的半径。边界做扩散运动时,ζ=1/2,这就是 Hallatschek 和 Nelson 在他们的模型中得到的平方根律;ζ=2/3 时变成立方根。Lab 还没有针对这一关系做过基准检验。
从半径为 r0 的起始区域长出的圆形菌落中,最终存活的扇区数。来源:O. Hallatschek · O. Hallatschek and D. R. Nelson