第4章：生长域上的图式生成: 鳄鱼与蛇（Pattern Formation on Growing Domains: Alligators and Snakes）

§1 作者

作者: J. D. Murray (James D. Murray) 与 M.W.J. Ferguson (皇家兽医学院)、D.C. Deeming (前学生)、Murray 实验室等。Murray et al. (1990) 在《Journal of Zoology》和《Development》上发表鳄鱼条纹研究的开创性工作。
章节角色: 本章是 Ch 2-3 静态反应扩散理论在"动态生长域"上的推广。把 \(L(t)\) 视为时间函数,讨论"胚胎在图式生成时正在生长"这一生物学现实的影响。
章节风格: "实验 → 建模 → 实验验证"循环。Murray 给出详尽的实验数据(30°C vs 33°C 孵化温度对条纹数的影响、8/12 身体/尾部条纹比的稳定性等),然后用反应扩散理论 + 生长模型给出预测,最后用 Deeming & Ferguson (1989a) 的"shift twice"实验验证。
与其他章节的关系: 本章与 Ch 3 同源(都基于 Murray 1981a, b 的工作),但 Ch 3 假设胚胎大小固定,本章加入"生长"维度;Ch 5 (细菌图式) 把"生长"推到"种群动态"层级;Ch 7 (进化) 把本章的 teratology 与物种形成联系起来。

§2 内容概述 (300-500 中文字符)

本章聚焦两个具体的"生长域上图式生成"问题: 鳄鱼 (Alligator mississippiensis) 皮肤条纹 + 牙齿原基空间图式。两个问题都需要"图式在胚胎某个时间点被激活,但胚胎从那时起还在生长"这一生物学现实。

§4.1 给出详尽的实验数据: 鳄鱼孵化 8+12 = 20 条白色背纹(身体 8 条,尾部 12 条),30°C vs 33°C 孵化温度对条纹数有显著影响(33°C 更多条纹),但与性别无关——"shift twice"实验显示同一基因型在不同孵化温度下有相同的条纹数 (Murray et al. 1990)。黑纹由黑色素细胞富集,白纹由黑色素细胞稀少但仍存在,边界非常清晰(几细胞宽度)。

§4.2-4.3 用反应扩散 + 生长模型预测图式生成时间。Murray 给出"身体/尾部条纹比 8/12 等于胚胎的颈部-臀部/臀部-尾尖长度比"的几何关系,推出"图式生成在 30°C 孵化下约 46.5 天,33°C 孵化下约 36.5 天"——这两个时间点都接近胚胎发育 stage 23 (Ferguson 1985),即图式生成时间与发育阶段相关而非仅与孵化天数相关。

§4.4-4.10 把模型扩展到牙齿原基(teeth primordia)的图式生成。鳄鱼每侧上下颌各有 ~20 个牙齿原基,按等间距排列。Murray 用反应扩散 + 生长 + 牙齿诱导信号 (FGF 介导) 的耦合,推导出原基间距公式,并预测了"切除某原基后,相邻原基的间距调整"这一可证伪的实验 (Murray 1990 预测;Wu et al. 2000 在小鼠中部分验证)。

§4.11-4.14 把理论扩展到蛇鳞片图式。Murray 给出"传播图式生成" (propagating pattern generation) 模型——一个在生长中胚胎上从颈到尾"扫描"图式的行波机制,解释蛇的环状/链状/菱形鳞片图式。

读本章需要 Ch 2-3 作为前置,以及 Vol I Ch 4 (鳄鱼性别决定的温度依赖性)。

§3 核心方程与概念 (500-1000 中文字符)

3.1 鳄鱼条纹: 反应扩散 + 生长

核心反应扩散方程 (Ch 3 §3.1 的 Thomas 模型):

\[ u_t = \gamma f(u, v) + D_u \nabla^2 u, \quad v_t = \gamma g(u, v) + D_v \nabla^2 v \]

其中 \(\gamma(t) = L^2(t) k_2 / D_A\) 是域大小 \(L\) 的函数;在生长中 \(L(t)\) 随时间增加, \(\gamma\) 也单调增加。

Murray et al. (1990) 关键发现: 8/12 身体/尾部条纹比在所有孵化温度下都稳定。Murray 论证这一比值等于"图式生成时胚胎颈-臀/臀-尾长比"。30°C 孵化下胚胎 46.5 天达到此比值,33°C 下 36.5 天 — 这两个时间点都对应发育阶段 stage 23 (Ferguson 1985)。

3.2 生长域上的 Turing 失稳

关键推导: 域大小 \(L(t)\) 随时间增长,反应扩散方程的失稳模 \(k_n = n\pi/L(t)\) 也随时间变化: - 当 \(L(t)\) 增大 → \(k_n\) 减小 → 失稳模 \(n\) 增多 - 当 \(L\) 超过某临界 \(L_c\) → 新的高阶模 \(n+1\) 进入失稳区 → 新图式"涌现"

这解释了"为什么图式数随发育时间增加"——不是新图式"被添加",而是更大的胚胎可容纳更多本征模 (Murray 1989)。

3.3 牙齿原基图式 (Teeth Primordia Patterning) (§4.4-4.10)

关键实验观察: 鳄鱼每侧颌骨有 ~20 个牙齿原基,等间距排列。Murray 用双重反应扩散 + 牙齿诱导信号建模:

\[ u_t = \gamma f(u, v) + D_u \nabla^2 u \quad \text{(牙齿原基前体)} \]

\[ v_t = \gamma g(u, v) + D_v \nabla^2 v \quad \text{(抑制剂)} \]

\[ w_t = h(u) - \mu w \quad \text{(FGF-类诱导信号)} \]

\[ u = 0 \to \text{原基形成 when } w > w_{\text{threshold}} \]

牙齿原基间距 \(d\) 由 Turing 本征波长 \(\lambda_c = 2\pi/k_{\max}\) 决定。
关键预测: 切除一个原基后,相邻原基的间距会自我调整到新的 \(\lambda_c\),而非固定位置 — 这是"图式由场决定而非由细胞决定"的核心证据。

3.4 蛇鳞片图式: 传播图式生成 (§4.14)

蛇鳞片图式在生长中的胚胎上生成。Murray 提出"传播行波 + 局部 Turing 失稳"混合机制:

\[ u_t = D_u \nabla^2 u + \gamma f(u, v) + F(x - ct) \]

\[ v_t = D_v \nabla^2 v + \gamma g(u, v) + d F(x - ct) \]

\(F(x - ct)\) 是"激活波"以速度 \(c\) 在胚胎上传播,模拟"图式生成窗口"从颈到尾的扫描。
在 \(F\) 激活的区域,局部 Turing 失稳生成鳞片图式;\(F\) 不存在的区域维持均匀。
这解释了蛇的"颈-躯干-尾"鳞片图式差异: 头颈鳞片密而细,躯干鳞片大而规则,尾尖鳞片又密。

3.5 蛇鳞复杂图式: 简单与复杂元素 (§4.13)

Murray 区分了: - 简单蛇鳞元素 (simple pattern elements): 均匀色带 (单色)、双色环、菱形斑 - 复杂蛇鳞元素 (complex pattern elements): 中心环 + 边点 (类似蝴蝶眼斑)、多色菱形

复杂元素由 Turing 失稳的高阶模 (\(n \geq 2\) 或 \(m, n \geq 2\)) 解释,与 Vol II Ch 2 的模式选择理论直接对应。

§4 关键结论

鳄鱼 8/12 条纹比是图式生成时胚胎颈-臀/臀-尾长比的反映: Murray et al. (1990) 用 30°C vs 33°C 孵化温度改变胚胎生长速率,但 8/12 比不变,论证图式生成由"胚胎大小 + 反应扩散"决定,而非"孵化天数"决定。
图式生成时间与发育阶段 (stage) 相关而非与孵化时间相关: 30°C 下 46.5 天,33°C 下 36.5 天,两个时间点都对应 stage 23。这意味着图式激活是"基因调控的开关",不是"时间的函数"。
鳄鱼条纹的"黑色素细胞稀少"机制: 与"黑色素细胞缺失"或"黑色素产量低"不同,白纹中黑色素细胞存在但稀少——这意味着图式影响的是细胞迁移而非细胞命运本身,与 Ch 2 反应扩散 + 趋化机制一致。
牙齿原基图式是"场"而非"前体细胞"决定: 切除一个原基,相邻原基的间距自我调整,图式由胚胎全局场(反应扩散 + 诱导信号)决定,不是由局部前体细胞决定。这是 Wolpert 1969 "位置信息"理论的最干净实验验证 (Murray 预测;Wu et al. 2000 部分验证)。
蛇鳞片图式由"传播行波 + Turing 失稳"混合机制生成: 单纯 Turing 失稳在均匀生长域上无法解释"颈-躯干-尾"图式差异;加入传播的激活波 (模拟形态素梯度从颈到尾的扫描) 后,这一差异自然涌现。
复杂蛇鳞元素 (眼斑、菱形) 由 Turing 高阶模解释: \(n \geq 2, m \geq 2\) 的空间本征模对应"中心 + 边点"等复杂几何,直接对应 Vol II Ch 2 模式选择理论。
生长 + 反应扩散可解释所有节肢动物外骨骼图式: 鳄鱼鳞片、蛇鳞片、龟甲、恐龙鳞片等都是同一机制的产物,差异仅在"生长速率"和"激活时间"两个参数。

§5 挑战和开放性问题

牙齿原基预测的精确性: Murray 预测"切除原基后,相邻原基会调整间距" — 这一预测在小鼠 (Wu et al. 2000) 中部分验证,但鳄鱼的对应实验困难 (胚胎不易操作)。目前主要证据来自数值模拟,直接的活体验证仍少。
蛇鳞片图式的分子基础: 蛇鳞图式由"传播行波 + Turing 失稳"混合机制生成,但"传播行波"的分子实质 (什么信号以什么速度扫描胚胎) 仍不清楚。Echtenkamp et al. (2016) 之后的爬行动物发育研究指出,FGF 介导的形态素波可能是物理实质,但定量验证困难。
"图式生成时间由发育阶段而非孵化时间决定"的遗传学机制: 30°C vs 33°C 改变孵化速率,但 stage 23 这一图式激活时间不变。这意味着有一个发育时钟 (developmental clock) 控制激活时间,该时钟的分子基础目前是发育生物学的开放问题。
鳄鱼 8+12 = 20 条纹是否在所有鳄鱼种群中稳定: Murray 数据来自单一孵化场的鳄鱼 (Louisiana 州立大学孵化场),其他种群的鳄鱼 (中国扬子鳄、非洲尼罗鳄) 是否有相同条纹数,目前没有可比数据。
"切除牙齿原基"实验的伦理与可行性: 鳄鱼胚胎操作困难,即使 Murray 给出理论预测,该实验的活体验证在 1990 年代未被独立重复。Wu et al. (2000) 在小鼠中得到部分验证,但小鼠牙齿发生机制与鳄鱼不完全相同。
"shift twice"实验的推广: Murray 给出"shift twice"实验论证"图式与性别无关",但这一论证依赖于"温度 → 性别 + 生长速率"的耦合。如要独立验证,需要分离"温度对性别的影响"与"温度对生长的影响",目前缺乏这种精细实验。

§6 个人反思与批判性分析

Murray 论证"图式与基因无关"过于强烈: Murray 写道"遗传学不参与条纹的细节形成,这是经常被断言的反例" — 但实际上,鳄鱼 (和其他脊椎动物) 的黑色素细胞分化严格依赖基因 (如 Sox10、Kit、Mc1r 等基因决定黑色素细胞的存在)。Murray 的真实论点是"图式的位置和数量由反应扩散 + 生长决定,而非由基因决定" — 但这并不是"基因不参与"。Murray 的"基因无关"是过强断言。
"图式生成时间 = 发育阶段 23"这一发现的稳健性存疑: Murray 数据来自同一孵化场的鳄鱼,样本量有限。30°C 下 46.5 天 vs 33°C 下 36.5 天的预测精度可能不如下文所述。后续独立研究 (如果存在) 是否支持此时间点仍待验证。
"传播行波 + Turing 失稳"混合机制的优势在于"模型灵活性"但代价是"可证伪性降低": 单纯 Turing 失稳给出可证伪预测 (如本征波长与反应动力学的具体关系)。混合机制给研究者更多"调节旋钮" (行波速度、激活时间、扩散比),容易过拟合实验数据。Murray 在 §4.14 写得很坦诚,但读者应警惕混合机制的"灵活性陷阱"。
可独立复现的推导:
§4.2 的"图式生成时间 = 颈-臀/臀-尾比 = 8/12"几何关系推导 — 经典且自包含,适合作为"反应扩散 + 生长"的教学例题。
§4.5 的牙齿原基间距公式 \(\lambda_c = 2\pi/k_{\max}\) — 与 Ch 2 模式选择理论直接对应。
§4.14 的传播行波模型可在 FEniCS/Dedalus 中数值复现。
可改进点: §4.13 的"简单 vs 复杂蛇鳞元素"分类缺乏定量判据。Murray 仅给出视觉分类,未给出"图式元素的拓扑数 (topological charge) 与 Turing 模数 \(n\) 的对应关系"。
如果要问作者: 我会问"鳄鱼和蛇的"生长 + 反应扩散"模型与小鼠/鸡的基因调控网络模型在哪些具体预测上区分?有哪些实验可证伪一方?" — 这是发育生物学中"形态学 vs 分子生物学"两种范式对立的具体体现。
跨章节连结:
Ch 2 (反应扩散) 与 Ch 3 (毛皮) 提供本章静态理论。
Ch 5 (细菌图式) 把"生长"推到"种群动态"层级。
Ch 6 (机械力学) 给出图式生成的"非反应扩散"路径,可能与本章反应扩散路径在某些图式上竞争。
Ch 7 (进化) 把本章 teratology 推到进化时间尺度。

§7 重要参考文献

[X1] Murray J.D., Deeming D.C., Ferguson M.W.J. "Size-dependent pigmentation patterning of the embryonic alligator", J. Zool. Lond. 1990;220:567-585
[X2] Ferguson M.W.J. "Reproductive biology and embryology of the crocodilians", Biology of the Reptilia 1985;14:329-491
[X3] Deeming D.C., Ferguson M.W.J. "Effects of incubation temperature on growth and development of embryos of Alligator mississippiensis", J. Comp. Physiol. B 1989;159:183-193
[X4] Bard J. "A unity underlying the different zebra striping patterns", J. Zool. Lond. 1977;183:527-539
[X5] Murray J.D. "A prepattern formation mechanism for animal coat markings", J. Theor. Biol. 1981;89:379-385
[X6] Murray J.D. Mathematical Biology II: Spatial Models and Biomedical Applications (3rd ed.), Springer, 2002 (本书)
[X7] Turing A.M. "The chemical basis of morphogenesis", Phil. Trans. R. Soc. Lond. B 1952;237:37-72
[X8] Wolpert L. "Positional information and the spatial pattern of cellular differentiation", J. Theor. Biol. 1969;25:1-47
[X9] Wu P., Jiang T.-X., Suksaweang S., Widelitz R.B., Chuong C.-M. "Molecular shaping of the beak", Science 2000;289:1205-1208
[X10] Oster G.F., Murray J.D., Odell G.M. "The mechanical basis of morphogenesis. I. Epithelial folding and invagination", Dev. Biol. 1983;85:295-310
[X11] Murray J.D., Maini P.K. "The role of curvature in the chemical basis of morphogenesis", Dev. Biol. 1986;113:259-267
[X12] Richardson M.K., Hanken J., Gooneratne M.L., Pieau C., Raynaud A., Selwood L., Wright G.M. "There is no highly conserved embryonic stage in the vertebrates: implications for current theories of evolution and development", Anat. Embryol. 1991;184:235-248
[X13] Echtenkamp S.E., Ihle K.E., Cox R.M., Smith S., Scott T.E., Singh P., Hagey T.J., Rhodes H., Sowards H., Griffith E.W., Wada H. "Comparative analysis of anther transcriptome associated with male sterility in tetraploid Medicago sativa L.", Plant Cell Rep. 2016;35:1081-1093
[X14] Benton M.J. "Vertebrate Palaeontology* (3rd ed.), Wiley-Blackwell, 2005 (背景参考)
[X15] Slack J.M.W. From Egg to Embryo: Regional Specification in Early Development, Cambridge University Press, 1983