[论文解读] Morphogenesis by coupled regulatory networks

Thimo Rohlf, Stefan Bornholdt|ArXiv.org|Jan 16, 2004

Gene Regulatory Network Analysis参考文献 6被引用 4

一句话总结

本文提出一种基于局部耦合调控网络的形态发生机制——以布尔网络和阈值网络建模——该机制可自组织形成具有预设比例的稳定空间模式，且不依赖于初始条件和系统规模。关键贡献在于证明此类网络通过分层信息处理，展现出对噪声和流动的鲁棒性，实现模式形成；在不同噪声和流动条件下，域尺寸表现出一阶和二阶相变的证据。

ABSTRACT

Based on a recently proposed non-equilibrium mechanism for spatial pattern formation [cond-mat/0312366] we study how morphogenesis can be controlled by locally coupled discrete dynamical networks, similar to gene regulation networks of cells in a developing multicellular organism. As an example we study the developmental problem of domain formation and proportion regulation in the presence of noise and cell flow. We find that networks that solve this task exhibit a hierarchical structure of information processing and are of similar complexity as developmental circuits of living cells. A further focus of this paper is a detailed study of noise-induced dynamics, which is a major ingredient of the control dynamics in the developmental network model. A master equation for domain boundary readjustments is formulated and solved for the continuum limit. Evidence for a first order phase transition in equilibrium domain size at vanishing noise is given by finite size scaling. A second order phase transition at increased cell flow is studied in a mean field approximation. Finally, we discuss potential applications.

研究动机与目标

• 通过局部耦合的离散动力学网络模拟形态发生，其类比于多细胞发育中的基因调控网络。
• 研究此类网络如何在存在噪声和细胞流动的情况下，实现无需预先存在的模式、具有稳定且尺寸不变的域比例的从头模式形成。
• 分析网络模型中噪声诱导动力学的统计力学特性，特别是边界重调整行为。
• 通过遗传算法优化与规则表变异性的统计分析，评估功能解的进化可及性。

提出的方法

• 将Rohlf & Bornholdt（2003）提出的随机细胞自动机模型改编为布尔网络与阈值网络框架，以模拟基因调控动力学。
• 采用遗传算法，系统规模（15–150个细胞）与种群规模（100）可变，演化满足四个约束的网络规则表：比例缩放、初始条件独立性、尺寸独立性与稳态性。
• 在连续极限下建立域边界重调整的主方程，以建模噪声诱导的动力学。
• 应用有限尺寸标度法，在噪声趋近于零时检测到平衡态域尺寸的一阶相变。
• 通过平均场近似研究在增加细胞流动下出现的二阶相变。
• 对80个高适应度解（Φ ≥ 0.96）进行统计分析，评估规则表相关性与进化鲁棒性。

实验结果

研究问题

• RQ1 在无可扩散形态发生素的情况下，局部耦合调控网络能否生成稳定且比例受控的空间模式？
• RQ2 噪声如何影响此类网络中域边界形成与重调整的动力学？
• RQ3 当噪声降低至零时，域尺寸的相变性质为何？
• RQ4 细胞流动如何影响模式形成机制的稳定性和转变行为？
• RQ5 通过突变与选择等进化过程，此类发育任务的功能解空间在多大程度上是可访问的？

主要发现

• 模型在优化后实现高鲁棒性的从头模式形成，域尺寸比例 α/(1−α) 与系统规模和初始条件无关，平均适应度 Φ ≈ 0.98。
• 通过有限尺寸标度法，在噪声趋近于零时发现平衡态域尺寸存在一阶相变的证据，表明模式稳定性发生急剧转变。
• 通过平均场近似观察到在增加细胞流动下存在二阶相变，提示流动诱导的模式破坏具有临界行为。
• 统计分析显示65%的规则表位置具有强相关性（C^ab = 1.0），表明存在保守的核心模块，而其他位置则表现出可变输出，提示进化约束。
• 解空间虽受限制但非平凡，核心规则与外围规则之间保持平衡，与动物发育调控模块的进化保守性相一致。
• 模型对初始条件波动与系统规模变化均表现出高度鲁棒性，支持其作为发育调控机制的生物学合理性。