形態形成:ヒョウの斑点を方程式で
胚は、ほとんど同じ細胞のかたまりとして始まる。それなのに、最後には片手に五本の指ができ、シマウマには縞が、ヒョウには斑点ができる。はじめにすべてが一様なら、形はどこから来るのか。1952 年、発表した最後の大きな論文「形態形成の化学的基礎」で、チューリングは数学的な答えを提案した。生物学者がそれを確かめるには何十年もかかった。
二つの物質、二つの速さ
考え方は、チューリングがモルフォゲンと呼んだ二つの化学物質が組織の中を広がっていくところを想像することである。今日ではふつう次のように説明される。
- 活性化因子:自分自身の生成を促し、もう一方の生成も促す。
- 抑制因子:活性化因子を抑える。
二つは互いに反応し、拡散する。つまり、水に落としたインクのように、濃いところから薄いところへ広がる。決定的なのは、抑制因子のほうが活性化因子より速く広がることである。
ある一点に、たまたま活性化因子が少し多く現れたらどうなるかをたどってみよう。そこでは活性化因子も抑制因子もたくさんつくられる。遅い活性化因子はその場に集まったままでいる。速い抑制因子は両側へ逃げていき、近くにもう一つ活性化因子の山ができるのを妨げる。その結果、山のまわりに山のない帯ができ、その先にまた山、また帯と続く。間隔のそろった斑点や縞が、何もないところから生まれるのである。
計算
チューリングは方程式を書き、一様な状態の近くで線形化した。一直線に並んだ組織で、$u$ と $v$ を活性化因子と抑制因子のずれとすると、
$$\begin{aligned} \frac{\partial u}{\partial t}&=a\,u+b\,v+D_u\frac{\partial^2u}{\partial x^2},\\ \frac{\partial v}{\partial t}&=c\,u+d\,v+D_v\frac{\partial^2v}{\partial x^2}. \end{aligned}$$$a$、$b$、$c$、$d$ を含む項が反応、$D_u$ と $D_v$ を含む項が拡散である。波長 $\frac{2\pi}k$ の小さな波が育つか消えるかを見るために、$e^{\lambda t}\cos kx$ という形の解を試す。$q=k^2$ とおくと、行列
$$\begin{pmatrix}a-D_uq&b\\c&d-D_vq\end{pmatrix}$$が得られ、この行列の行列式が負になると波は育つ。
具体的な数をとろう。$a=1$、$b=-1$、$c=2$、$d=-1.5$。拡散がなければ($q=0$)、行列式は $1\cdot(-1.5)-(-1)\cdot2=0.5>0$、トレースは $-0.5<0$ である。一様な状態は安定で、化学反応だけでは模様はまったく生まれない。
ここで拡散を入れ、抑制因子を 10 倍速くする。$D_u=1$、$D_v=10$ とすると、
$$\det(q)=(1-q)(-1.5-10q)+2=10q^2-8.5q+0.5.$$根は $q\approx0.064$ と $q\approx0.786$ である。そのあいだでは行列式は負になり、その波長の波は育つ。いちばん速いのは中央の $q=0.425$ で、波長は $\frac{2\pi}{\sqrt{0.425}}\approx9.6$ である。ふつうは違いをならしてしまう拡散が、ここでは違いを生み出す。だからこの現象は拡散誘導不安定性と呼ばれる。
計算から動物へ
チューリングは計算の一部をマンチェスターの計算機で行い、得られた模様の一つは、ある種の動物に似た不規則なまだら模様だった。どの物質がモルフォゲンなのかを彼は知らず、当時の生物学には確かめる手段がなかった。この考えは何十年も、ほとんど忘れられていた。
最初の確認は化学から来た。1990 年、ボルドーのフランスの研究グループが、実験室の反応を使って、ゲルの中に安定したチューリング模様をつくりだした。生物学では 1995 年、日本の近藤滋と浅井理恒が、キンチャクダイの仲間の縞が方程式の予言どおりに動き、枝分かれすることを示した。それ以来、指、毛、羽、そして上あごのひだの間隔に、チューリング型のしくみが見いだされている。
練習:拡散の速さが同じなら模様はできない
$D_u=D_v=1$ で計算をやり直せ。ある $q>0$ で行列式は負になるか。
$\det(q)=(1-q)(-1.5-q)+2=q^2+0.5q+0.5$。係数はすべて正なので、すべての $q\geq0$ で $\det(q)>0$ である。どの波も育たない。速さの違いがなければ、拡散は不規則さをならすだけである。
練習:いちばん速い波はどれか
$\det(q)=10q^2-8.5q+0.5$ で、行列式が最小になる $q$ と、それに対応する波長を求めよ。
放物線 $\alpha q^2+\beta q+\gamma$ の最小は $q=-\frac{\beta}{2\alpha}=\frac{8.5}{20}=0.425$ にある。すると $k=\sqrt{0.425}\approx0.652$、波長は $\frac{2\pi}{0.652}\approx9.6$ である。長さ 100 の組織には、斑点がおよそ 10 個入る。
出典: Alan Turing, “The chemical basis of morphogenesis”, Philosophical Transactions of the Royal Society B, 237(1952), pp. 37–72 · Shigeru Kondo and Rihito Asai, “A reaction–diffusion wave on the skin of the marine angelfish Pomacanthus”, Nature(1995)· Philip Ball, “Forging patterns and making waves from biology to geology”, Philosophical Transactions B(2015)。