第四章：短吻鳄的条纹模式形成与蛇类的色素沉着模式

阅读笔记 — J.D. Murray《数学生物学 II》第四章

4.1 短吻鳄条纹模式形成：实验

实验背景

本章研究的主角是美洲短吻鳄（Alligator mississippiensis）。在孵化过程中，温度显著影响胚胎发育和性别决定。当孵化温度为 33°C 时，胚胎倾向于发育为雄性；当温度为 30°C 时则倾向于雌性。实验人员（Deeming 和 Ferguson, 1988–1990）对不同温度下培养的鳄鱼胚胎进行了系统的形态测量。

形态测量

实验测量了胚胎在 30°C 和 33°C 下培养第 10 至 50 天的以下数据： - 胚胎总长度（total length） - 颈-臀部长度（nape–rump length） - 臀部-尾尖长度（rump–tail tip length）

回归分析表明，高温（33°C）下胚胎生长更快。在 33°C 条件下，颈-臀与臀-尾的比例达到 8:12 所需的时间为约 36.5 天，而在 30°C 条件下则需要约 46.5 天。

条纹数量与温度的关系

短吻鳄体表条纹数目与孵化温度之间存在显著关系：

温度	条纹数目（约）	尾尖颜色
33°C（雄性）	~20 条（8 条在躯干 + 12 条在尾部）	常为白色
30°C（雌性）	稍少	常为黑色

条纹数目与性别本身无关（33°C 下雄性平均 19.96±1.15 条，雌性 20.00±0.71 条），而是与孵化温度所决定的胚胎大小直接相关。

组织学证据

在第 28 阶段的胚胎皮肤样本中，黑色素细胞（melanocytes）的分布在黑白条纹区域存在明显差异： - 黑色条纹：黑色素细胞丰富，黑色素浓度高 - 白色条纹：黑色素细胞稀少，虽然存在但不大量产生黑色素

黑白条纹之间的边界极为清晰，表明条纹形成机制具有高度的空间选择性。

4.2 建模概念：条纹形成时间的确定

模式生成的基本概念

几乎所有模式生成模型的直觉基础都可以用短程激活-长程抑制（short-range activation and long-range inhibition）来解释。本章采用了一个细胞-趋化-扩散系统（cell-chemotaxis-diffusion system）作为说明性机制，其中细胞产生自己的趋化因子。

traveling wave 模式

从实验观察可见，短吻鳄胚胎的条纹模式首先出现在颈部，然后以波的形式向下传递到尾部。这一过程可以用一维的行波模式来描述。

图 4.2(a) 展示了一个典型的一维行波演化过程，该过程来自以下细胞-趋化模型系统（无量纲形式）：

\[\frac{\partial n}{\partial t} = D\frac{\partial^2 n}{\partial x^2} - \alpha\frac{\partial}{\partial x}\left(n\frac{\partial c}{\partial x}\right) \qquad(4.31)\]

\[\frac{\partial c}{\partial t} = \frac{\partial^2 c}{\partial x^2} + \frac{n}{n+1} - c\]

其中： - \(n = n(x,t)\)：细胞（黑色素细胞）密度 - \(c = c(x,t)\)：趋化因子浓度 - \(D\)：细胞的相对扩散系数 - \(\alpha\)：趋化系数

初始条件为均匀稳态 \(n=1, c=0.5\)，在 \(x=0\) 处施加一个小扰动。模式形成的波前从颈部向右移动，在身后留下规则的正弦式条纹模式。

条纹形成时间的估算

关键假设：模式生成机制所需的时间（数小时量级）远小于显著胚胎生长所需的时间。这意味着在机制激活时，胚胎的长度决定了能够容纳的波长数目，从而决定了条纹的数量。

根据实验数据，当颈-臀长度与臀-尾长度之比为 8:12 时：

\[L = \frac{\text{尾长} + \text{颈-臀长}}{\text{孵化时的条纹数}} \approx 4.7\text{–}4.9 \text{ mm} \qquad(4.32)\]

33°C 下第 32 天：\(L = (56.25 + 39.50)/20.35 = 4.71\) mm
30°C 下第 40 天：\(L = (53.60 + 37.20)/18.55 = 4.89\) mm

额外增加的 7.75 mm 长度可以解释雄性胚胎多出的约 1.64 条条纹。

4.3 短吻鳄的条纹与"阴影条纹"

主条纹的形成

在假设模式生成时间远快于胚胎生长时间的前提下，机制在胚胎背部产生规则等距的条纹模式。黑白条纹与细胞密度的高峰和低谷相对应： - 白条纹（低谷）：黑色素细胞密度不足，无法产生显著色素 - 黑条纹（高峰）：黑色素细胞密度高，产生深色色素

默认状态下胚胎应为深色色素形式（如果模式机制不激活），黑色素形式在自然界中极为罕见。

阴影条纹（Shadow Stripes）

如果放松"模式形成时间远快于胚胎生长"这一假设，则意味着模式形成发生在不断生长的胚胎表面之上。此时躯干部分变宽，在已形成的主条纹之间会形成额外的较不明显的阴影条纹，这些阴影条纹位于主条纹之间的区域。

数学模型的预测与以下观察相符： 1. 阴影条纹位于腹侧主条纹之间 2. 尾部条纹的颜色清晰度较低（因为模式从头部开始形成，越靠尾部形成越晚）

这一现象与第三章讨论的域生长对模式的影响一致。

角鱼（Angelfish）的条纹增加

有趣的是，年轻的角鱼（Pomacanthus）初始只有 3 条条纹，随生长逐渐增加新的条纹到 6 条、12 条。新条纹出现在已有条纹之间，这一过程涉及细胞趋化和反应扩散系统的耦合。

4.4 短吻鳄牙齿原基的空间模式

研究背景

牙齿原基（teeth primordia）在脊椎动物颌部的发育是另一个在动态生长中形成高度规则空间模式的例子。每个牙齿原基的形成涉及细胞协调，形成精确的时空序列。

短吻鳄的牙齿发育与人类腭裂和唇裂等先天畸形的研究高度相关，因为短吻鳄具有类似哺乳动物的颌骨和牙齿特征，且可以在外部卵中进行非侵入性和侵入性的在体研究。

实验数据

Westergaard 和 Ferguson（1986, 1987, 1990）对短吻鳄牙齿原基的启动和空间分布进行了详细的实验研究：

第一个牙齿原基（称为"牙决定因子"，dental determinant）形成于下颌前部
牙齿启动从牙决定因子向前后两个方向扩展
间质原基（interstitial primordia）在成熟邻居附近的空间中形成

前 7 个原基形成后，第 8 个原基开始非交替序列，随后又恢复到交替序列。整个孵化期（约 65–70 天）内，约有 19 个早期原基被重新吸收而不发挥功能，另有 7 个过渡组功能较短。

生长与牙齿数目的关系

实验数据表明，牙齿原基的数目与孵化时间呈指数关系：

\[N(t) = N_0 e^{rt}\]

下颌：\(N_0 = 0.0066\)，\(r = 0.3077/\text{day}\)
上颌：\(N_0 = 0.0047\)，\(r = 0.3442/\text{day}\)

在更长的发育周期中，牙齿数目遵循 Gompertz 型增长：

\[N(t) = N_1 \exp[-N_2 \exp(-rt)]\]

这表明早期颌部呈指数生长。

4.5 蛇类的色素沉着模式形成

生物学背景

蛇类（以及一般的爬行动物）展现出极为丰富和多样的色素模式。同一物种内部也经常存在极端的模式多态性。例如加州王蛇（Lampropeltis getulus californiae）就有侧向条纹和纵向条纹两种表型。

爬行动物的皮肤是体内最大的器官。色素细胞前体——色素母细胞（chromatoblasts）在发育期间从神经嵴迁移，分布于真皮中。皮肤是否形成色素斑取决于色素前体细胞是否产生色素。

蛇类模式研究的重要性

与蝴蝶翅膀模式不同，许多常见的蛇类模式（如菱形斑、复杂条纹等）难以用标准反应扩散模型直接生成，而需要考虑： - 生长的域（domain growth） - 变化的边界条件 - 空间变化的参数

本章表明，当模式形成机制在图案沉积期间域发生显著生长时，可以产生令人惊讶的多样性和复杂性的新图案。

4.6 细胞-趋化模型机制

模型方程

Murray 和 Myerscough（1991）提出了以下细胞-趋化模型来解释蛇类的色素沉着模式：

\[\frac{\partial n}{\partial t} = D_n \nabla^2 n - \alpha \nabla \cdot (n \nabla c) - r n (N - n) \qquad(4.19)\]

\[\frac{\partial c}{\partial t} = D_c \nabla^2 c + \frac{S n}{\beta + n} - \gamma c \qquad(4.20)\]

其中： - \(n = n(\mathbf{x},t)\)：色素细胞（色素母细胞）密度 - \(c = c(\mathbf{x},t)\)：趋化因子浓度 - \(D_n, D_c\)：扩散系数 - \(\alpha\)：趋化系数（衡量趋化效应的强度） - \(r\)：细胞有丝分裂率（逻辑增长形式） - \(N\)：初始均匀细胞密度 - \(S\)：细胞分泌化学物质的最大速率 - \(\gamma\)：趋化因子的自然降解速率 - \(\beta\)：与趋化因子产生相关的米歇利斯常数

无量纲化

令：

\[x^* = \left[\frac{\gamma}{D_c s}\right]^{1/2} x, \quad t^* = \frac{\gamma t}{s}, \quad n^* = \frac{n}{\beta}, \quad c^* = \frac{\gamma c}{S}\]

\[N^* = \frac{N}{\beta}, \quad D^* = \frac{D_n}{D_c}, \quad \alpha^* = \frac{\alpha S}{\gamma D_c}, \quad r^* = \frac{r \beta}{\gamma}\]

则无量纲方程为（省略星号）：

\[\frac{\partial n}{\partial t} = D \nabla^2 n - \alpha \nabla \cdot (n \nabla c) + s r n (N - n) \qquad(4.22)\]

\[\frac{\partial c}{\partial t} = \nabla^2 c + s \frac{n}{1+n} - c \qquad(4.23)\]

其中 \(s\) 是与皮肤大小相关的尺度因子。

线性稳定性分析

方程组存在正的齐次稳态：

\[n_0 = N, \quad c_0 = \frac{N}{1+N} \qquad(4.23)\]

在稳态附近线性化（令 \(n = N + u\)，\(c = c_0 + v\)，\(|u|,|v| \ll 1\)），得到：

\[\frac{\partial u}{\partial t} = D \nabla^2 u - \alpha N \nabla^2 v - s r N u \qquad(4.24)\]

\[\frac{\partial v}{\partial t} = \nabla^2 v + s \frac{u}{(1+N)^2} - v \qquad(4.25)\]

设解的形式为 \(\mathbf{u} \propto \exp[i\mathbf{k} \cdot \mathbf{x} + \lambda t]\)，代入可得关于 \(\lambda\) 的特征多项式：

\[\lambda^2 + [(D+1)k^2 + s r N + s]\lambda + [Dk^4 + \{s r N + D s - \frac{s N \alpha}{(1+N)^2}\}k^2 + r N s^2] = 0 \qquad(4.26)\]

当某些波数 \(k\) 对应的 \(\lambda > 0\) 时，均匀稳态失稳，形成空间图案。

行波解与 propagating pattern

在没有细胞有丝分裂（\(r=0\)）的情况下，系统变为：

\[\frac{\partial n}{\partial t} = D \nabla^2 n - \alpha \nabla \cdot (n \nabla c) \qquad(4.31)\]

\[\frac{\partial c}{\partial t} = \nabla^2 c + s \frac{n}{1+n} - c \qquad(4.32)\]

当初始扰动施加在域的一端时，会产生一个行波，在身后留下规则的条纹模式。色散关系为：

\[\lambda^2 + [(D+1)k^2 + s]\lambda + Dk^4 + \left\{Ds - \frac{sN\alpha}{(1+N)^2}\right\}k^2 = 0 \qquad(4.34)\]

4.7 简单与复杂的蛇类图案要素

侧向条纹

通过调整有丝分裂率 \(r\) 和趋化参数 \(\alpha\)，细胞-趋化模型可以产生各种侧向条纹图案，这与许多蛇类（如 Vermicella annulata、Lampropeltis getulus californiae）的条纹模式相符。

纵向条纹

对于某些参数组合，模型产生与长轴平行的纵向条纹，如束带蛇（Thamnophis sauritus sauritus）和四种线蛇（Elaphe quatuorlineata）所展示的纵向条纹模式。

有趣的是，侧向条纹和纵向条纹之间的转换只需要参数的微小变化，这解释了为什么加州王蛇同一物种内可以出现两种不同的条纹类型。

斑点模式

在适当的参数条件下，模型还可以产生规则的斑点图案，如喀普山蝰蛇（Bitis atropos atropos）所展示的交替半圆形斑点模式。

配对斑点

当趋化参数 \(\alpha\) 增加时，原本规则的斑点会转变为配对斑点模式，如 Leopard Snake（Elaphe situla）所展示的两两配对的斑点。

菱形图案与波形条纹

当域在横向生长时，侧向条纹的域会产生菱形图案（diamond patterns），如东部菱背响尾蛇（Crotalus adamanteus）和蹄形蛇（Coluber hippocrepis）所展示。当域变得较窄时，菱形图案转变为波形条纹，这与许多响尾蛇尾部图案变化的观察一致。

模式形成的直觉解释

低趋化响应（\(\alpha\) 小）：细胞需要较大面积的高密度区域才能产生足够的趋化因子梯度来招募其他细胞，因此倾向于形成条纹。
高趋化响应（\(\alpha\) 大）：较少的细胞就能产生足够陡峭的趋化因子梯度来招募其他细胞聚集，因此倾向于形成斑点。

公式汇总表

编号	公式名称/内容	方程	备注
(4.1)	指数增长定律	\(\frac{dL}{dt} = rL \Rightarrow L(t) = L_0 e^{rt}\)	颌部长度生长
(4.2)	生长域上的反应扩散方程	\(c_t = D_c \xi \xi + \gamma f(c)\)	含稀释效应
(4.3)	含稀释项的反应扩散	\(c_t = D_c \xi \xi + \gamma f(c) - rc\)	\(-rc\) 为稀释项
(4.5)	牙齿 patterning 完整模型	\(\begin{aligned} u_t &= \gamma(hc - u + u^2 v) + \partial_{xx} u \\ v_t &= \gamma(b - u^2 v) + d^2 \partial_{xx} v \\ c_t &= -\delta c + p^2 \partial_{xx} c \end{aligned}\)	三方程耦合系统
(4.6)–(4.8)	无量纲化后的牙齿模型	含稀释项 \(-ru\), \(-rv\), \(-rc\) 和时变扩散系数 \(e^{-2rt}\)	生长域变换
(4.11)	新牙齿源的逻辑增长	\(\frac{d\bar{c}_i}{dt} = k_1 \bar{c}_i(1 - \bar{c}_i/k_2)\)	抑制区建模
(4.19)	细胞趋化方程	\(\partial_t n = D_n \nabla^2 n - \alpha \nabla \cdot (n \nabla c) - rn(N-n)\)	色素细胞
(4.20)	趋化因子方程	\(\partial_t c = D_c \nabla^2 c + \frac{Sn}{\beta+n} - \gamma c\)	趋化因子分泌
(4.22)	无量纲细胞趋化模型	\(\partial_t n = D \nabla^2 n - \alpha \nabla \cdot(n \nabla c) + srn(N-n)\)	主 modeling equation
(4.23)	齐次稳态	\(n_0 = N,\quad c_0 = \frac{N}{1+N}\)	线性化基础
(4.26)	色散关系（完整）	\(\lambda^2 + [(D+1)k^2 + srN + s]\lambda + [Dk^4 + \{srN + Ds - \frac{sN\alpha}{(1+N)^2}\}k^2 + rNs^2] = 0\)	决定模式形成
(4.28)	等根条件（孤立模式）	\([srN + Ds - \frac{sN\alpha}{(1+N)^2}]^2 - 4Drs^2 N = 0\)	确保单一波数失稳
(4.29)	临界波数	\(k_c^2 = [s^2 r N / D]^{1/2}\)	模式特征尺度
(4.31)	无有丝分裂的行波系统	\(\partial_t n = D \nabla^2 n - \alpha \nabla \cdot(n \nabla c),\quad \partial_t c = \nabla^2 c + s\frac{n}{1+n} - c\)	traveling wave 分析
(4.34)	无有丝分裂色散关系	\(\lambda^2 + [(D+1)k^2 + s]\lambda + Dk^4 + \{Ds - \frac{sN\alpha}{(1+N)^2}\}k^2 = 0\)	行波分析
L	条纹波长（估算）	\(L = \frac{\text{尾长} + \text{颈-臀长}}{\text{条纹数}}\)	条纹间距经验公式

结论

本章以美洲短吻鳄和蛇类为案例，系统展示了在生长域上的模式形成这一核心主题。主要结论如下：

短吻鳄条纹：条纹数目由孵化温度决定的胚胎大小所控制，与性别本身无关。模式形成时间早于条纹可见时间。
阴影条纹：当模式形成过程与胚胎生长同时发生时，会产生主条纹之间的次级阴影条纹，这与域生长对模式选择的影响一致。
牙齿原基：牙齿形成的时空序列可以通过一个包含颌部生长的反应-扩散机制精确再现。该模型预测了移植物和物理屏障对牙齿序列的影响。
蛇类色素模式：细胞-趋化模型能够产生多种蛇类皮肤色素模式，包括侧向条纹、纵向条纹、斑点、配对斑点和菱形图案，且这些模式对参数变化高度敏感。
生长与模式形成的耦合：无论是短吻鳄的条纹还是蛇类的色素模式，胚胎生长在模式最终形态的决定中都起着至关重要的作用。