第六章：间充质形态发生 (Mesenchymal Morphogenesis)
阅读笔记
6.1 引言、动机与背景生物学
本章探讨了形态发生学中的机械理论方法。与传统的化学形态因子理论不同，机械理论认为模式形成与形态发生是同时进行的单一过程。细胞与胚胎组织的形态运动持续相互作用，产生可观察的空间模式。这一方法的核心优势在于其模型基于可测量量，如细胞密度、力、组织变形等，使研究更易进行实验验证。
两种胚胎细胞类型
书中描述了两类早期胚胎细胞：
-
间充质细胞（Mesenchymal cells）：也称为成纤维细胞或真皮细胞。这类细胞具有独立运动能力，通过丝状伪足（filopodia）或片状伪足（lamellapodia）抓住粘附位点并牵拉自身移动。细胞还能分泌细胞外基质（ECM）纤维材料。
-
表皮细胞（Epidermal cells）：通常不移动，而是紧密排列成片层结构。空间模式通过细胞变形而非移动来实现。
机械理论的一个关键优势是其自我纠正能力——胚胎能够调整许多外部干扰，这与预模式存在后再进行形态发生的开放回路系统不同。
6.2 间充质形态发生的机械模型
细胞运动的影响因素
间充质细胞的运动受多种因素影响：
- 对流（Convection）：细胞可被动地随变形基质移动
- 趋化作用（Chemotaxis）：化学梯度可引导细胞运动
- 接触引导（Contact guidance）：基质提供优选运动方向
- 接触抑制（Contact inhibition）：高细胞密度抑制运动
- 粘附趋性（Haptotaxis）：细胞沿粘附梯度移动
- 扩散（Diffusion）：细胞随机运动
- 电趋性（Galvanotaxis）：电场提供运动优选方向
基本机械模型
模型基于两个关键实验事实： 1. 细胞在由纤维细胞外基质（ECM）组成的组织基质中迁移 2. 细胞能产生强大的牵拉力
模型包含三个方程：
- 细胞守恒方程
- 细胞与ECM之间的力平衡方程
- ECM守恒方程
设 \(n(\mathbf{r}, t)\) 为细胞密度，\(\rho(\mathbf{r}, t)\) 为ECM密度，\(\mathbf{u}(\mathbf{r}, t)\) 为ECM位移向量。
细胞守恒方程
一般形式为：
其中 \(\mathbf{J}\) 为细胞通量，\(M\) 为有丝分裂率。
对流通量
随机扩散
经典的Fickian扩散贡献 \(-\nabla^2 n\) 项。长程扩散效应通过非局部扩散建模：
粘附趋性（Haptotaxis）
细胞沿ECM密度梯度移动：
完整的细胞守恒方程为：
细胞-基质力学相互作用方程
在低Reynolds数条件下（胚胎运动的特征），可忽略惯性效应，力平衡方程为：
应力张量包含ECM贡献和细胞牵拉贡献：
细胞牵拉力建模为：
ECM守恒方程
在忽略分泌效应后（S=0），方程简化为：
无量纲化
使用特征长度 \(L\) 和时间尺度 \(T\)，初始矩阵密度 \(\rho_0\)，得到无量纲参数群，减少了14个参数至12个。
6.3 线性分析、色散关系与模式形成潜力
线性稳定性分析
均匀稳态解为：
线性化后寻求形式为 \(\exp[\sigma t + i\mathbf{k} \cdot \mathbf{r}]\) 的解，得到色散关系：
其中：
模式形成条件
空间非均匀解存在的必要条件是细胞牵拉力 \(\tau > 0\)。当参数确保 \(b(k^2) < 0\) 和/或 \(c(k^2) < 0\) 时，空间结构解存在。
模式形成机制直觉理解
- 当细胞牵拉力小于临界值时，任何细胞密度扰动都会被平滑消除
- 当牵拉力超过临界值时，细胞产生的力大于ECM的弹性阻力，空间不均匀性开始形成
- 进一步增加牵拉力产生更复杂的模式
6.4 产生复杂色散关系的简单机械模型
通过设置不同参数为零，可推导出多种简化模型。关键发现包括：
模型 (i)：无扩散、无趋化、无有丝分裂
当 \(D_1 = D_2 = a_1 = a_2 = r = 0\) 时：
色散关系简化为：
空间模式形成的临界牵拉力为：
关键色散关系类型
- 有限范围不稳定模式：只有特定波数范围不稳定
- 无限范围不稳定模式：所有波数都可能不稳定
- 复数 \(\sigma\)：产生行波解
表6.1：具有有限不稳定波数范围的机械模型参数条件
| 条件 | \(D_1\) | \(D_2\) | \(a_1\) | \(a_2\) | \(r\) | \(\mu\) | \(s\) | \(\gamma\) | \(\lambda\) | \(\tau\) 条件 |
|---|---|---|---|---|---|---|---|---|---|---|
| (a) | ○ | ○ | ○ | ○ | ○ | ● | ● | ● | ○ | \(\tau > \{1 + \gamma s + \cdots\}/2\) |
| (b) | ● | ○ | ○ | ○ | ○ | ● | ○ | ● | ○ | \(1 > \tau > 1/2\) |
注：● = 正非零，○ = 零
6.5 羽胚的周期模式
生物学背景
脊椎动物皮肤由两层组成：上皮表皮覆盖较厚的间充质真皮。羽毛形成的第一步是真皮凝聚物（papillae）的形成，随后是表皮的 placode 结构。
在鸡胚中，羽毛原基在受精后约6天变得可见。真皮凝聚物主要由细胞迁移形成，局部增殖起次要作用。
模式形成机制
Murray-Oster机械理论应用于羽毛胚芽形成的方案：
- 第一阶段：沿背正中线形成细胞柱，然后分裂成一行papillae
- 第二阶段：侧行papillae依次形成，与前一行的papillae交错排列
- 第三阶段：形成六角形周期阵列
模型预测与验证
数值模拟表明： - 细胞聚集与ECM密度变化同相 - ECM位移与前两者反相 - 初始细胞列形成的预应力应变场诱导相邻行列的papillae在半个波长处交错排列
关键预测
模型预测当总细胞数 \(N\) 减小时，间距增加，这与Davidson的实验观察一致。
6.6 肢体形态发生中的软骨凝聚与形态发生规则
软骨模式形成
脊椎动物肢体是研究最广泛的发育系统之一。软骨模式形成涉及软骨细胞（一种间充质细胞）的聚集。
分叉序列
随着肢芽生长（从远端开始），软骨模式依次形成：
- 单凝聚（如肱骨）→ 2. 双凝聚（如桡骨和尺骨）→ 3. 多凝聚（如指骨）
几何与尺度作为分叉参数
机械模型表明： - 初始圆形但有椭圆偏置的横截面 - 细胞凝聚产生的强中心方向应力 - 径向应力使横截面变得更椭圆 - 椭圆度变化诱导二次分叉为两个凝聚
形态发生规则
三种基本的细胞凝聚类型：
- 焦点凝聚 (F)：孤立的细胞聚集
- 分叉分叉 (B)：Y形分叉
- 节段分叉 (S)：沿轴线的重复分段
这些基本元素可以组合构建任何肢体软骨模式。
嫁接实验
将一小块组织从肢芽嫁接到另一位置（ZPA，极化活性区）会导致双肢形成。模型预测：若嫁接后的肢芽在几何上被限制于单肢对应的尺度，则无法形成双肢分叉序列。
6.7 胚胎指纹形成
皮肤纹理的生物学
指纹图案（皮肤纹理学）研究历史悠久。表皮脊纹在妊娠约第三个月开始发育，最终模式取决于不对称程度和垫块形成：对称垫产生螺旋，不对称垫产生环形，其他形式的垫块发育导致弧形。
模式形成的两阶段过程
- 脊纹形成
- 脊纹弯曲成螺旋、环形等
机械模型
基于细胞在细胞外基质上的机械运动：
细胞-基质力平衡和ECM守恒方程补充完整。
色散关系
包含矩阵分泌的模型产生三次多项式：
奇点拓扑
指纹图案具有独特的拓扑奇点，不同于向量场。Penrose表明这些是张量特征（如应力或应变）的表现。
| 奇点类型 | 索引 |
|---|---|
| 中心 | +1 |
| 鞍点 | -1 |
| 脊奇点 | +1/2 |
| 三辐脊奇点 | -1/2 |
关键结论
由于最终模式依赖于初始条件（包含随机因素），没有两个人——甚至同卵双胞胎——可以有完全相同的指纹。
公式汇总表
核心方程
| 方程编号 | 名称 | 公式 |
|---|---|---|
| (6.1) | 细胞守恒方程（一般形式） | \(\frac{\partial n}{\partial t} = -\nabla \cdot \mathbf{J} + M\) |
| (6.2) | 对流通量 | \(\mathbf{J}_c = n\frac{\partial \mathbf{u}}{\partial t}\) |
| (6.5) | 随机扩散通量 | \(\mathbf{J}_D = -D_1 \nabla n + D_2 \nabla(\nabla^2 n)\) |
| (6.6) | 粘附趋性通量 | \(\mathbf{J}_h = n(a_1 \nabla\rho - a_2 \nabla^3\rho)\) |
| (6.7) | 完整细胞守恒方程 | 包含对流、扩散、趋化、有丝分裂各项 |
| (6.10) | 力平衡方程 | \(\nabla \cdot \boldsymbol{\sigma} + \rho\mathbf{F} = 0\) |
| (6.14) | 细胞牵拉力 | \(\tau(n) = \frac{\tau_0 n}{1 + \lambda n^2}\) |
| (6.20) | ECM守恒方程 | \(\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \frac{\partial \mathbf{u}}{\partial t}) = S\) |
| (6.24) | ECM守恒（S=0） | \(\rho_t + \nabla \cdot (\rho \mathbf{u}_t) = 0\) |
| (6.31) | 色散关系多项式 | \(\sigma[\mu k^2 \sigma^2 + b(k^2)\sigma + c(k^2)] = 0\) |
| (6.36) | 简化模型色散关系 | \(\sigma(k^2) = -\frac{b(k^2)}{\mu k^2}\) |
| (6.39) | 临界牵拉力 | \(\tau_c = \frac{1}{2}(1+\lambda)^2[\cdots]\) |
| (6.63) | 指纹模型细胞方程 | \(n_t = D_1\nabla^2 n - \nabla \cdot (n\mathbf{u}_t) - a_1\nabla \cdot n\nabla\rho + rn(N-n)\) |
| (6.66) | 指纹模型色散关系 | \(a(k^2)\sigma^3 + b(k^2)\sigma^2 + c(k^2)\sigma + d(k^2) = 0\) |
模型参数说明
| 参数 | 含义 | 单位/典型值 |
|---|---|---|
| \(n\) | 细胞密度 | 细胞数/单位体积 |
| \(\rho\) | ECM密度 | 质量/单位体积 |
| \(\mathbf{u}\) | ECM位移向量 | 长度 |
| \(D_1\) | Fickian扩散系数 | 长度²/时间 |
| \(D_2\) | 长程扩散系数 | 长度⁴/时间 |
| \(a_1\) | 粘附趋性系数 | 长度²/时间 |
| \(a_2\) | 长程粘附趋性系数 | 长度⁴/时间 |
| \(\tau\) | 细胞牵拉力参数 | 力/细胞 |
| \(\lambda\) | 接触抑制参数 | - |
| \(\gamma\) | 非局部细胞-ECM相互作用参数 | 长度² |
| \(s\) | 基质附着弹性参数 | 力/长度³ |
| \(\mu_1, \mu_2\) | 粘滞系数 | 力·时间/长度² |
| \(E\) | 杨氏模量 | 力/长度² |
| \(\nu\) | 泊松比 | - |
| \(r\) | 细胞增殖率 | 1/时间 |
| \(N\) | 最大细胞密度 | 细胞数/单位体积 |
无量纲参数群
关键概念总结
-
机械理论的核心思想：模式形成与形态发生同时进行，细胞牵拉力是模式形成的主要驱动力
-
色散关系的意义：决定哪些波数模式会线性不稳定，从而决定最终空间模式的特征波长
-
临界牵拉力 \(\tau_c\)：当细胞牵拉力超过临界值时，均匀稳态变得空间不稳定，开始形成模式
-
形态发生规则：如软骨模式中的F、B、S基本元素组合，可构建各种肢体软骨模式
-
鲁棒性与耦合：组织间相互作用通过减少可能的模式选择来增强形态发生的鲁棒性
-
生物学可测试性：机械模型的参数原则上都是可测量的，这使得理论更容易通过实验验证
延伸阅读建议
- 第3章关于反应-扩散系统与机械模型的比较
- 第10章关于各向异性弹性效应
- 第8章关于血管生成中的应用
- 第10章关于伤口愈合中的应用
本笔记基于 J.D. Murray《Mathematical Biology II》第三章内容整理