Geometry and dynamical morphology of growing bacterial colonies

技術要約：成長する細菌コロニーの幾何学と動的形態学

問題提起
細菌コロニーは、自由エネルギーの最小化ではなく、継続的なエネルギー散逸によって駆動される非平衡パターン形成のアクセシブルなモデルとして機能し、多様な形態を示す。既存の文献では、静的な記述子（最終形状、粗さ、フラクタル次元など）や、栄養輸送および細胞間相互作用を組み込んだメカニズムモデルを用いてコロニー形態を広範に特徴付けてきたが、決定的な空白が残されている。すなわち、視覚的に異なるコロニーの形状が、制御された逸脱の下での共通の幾何学的枠組みから生じたものなのか、それとも根本的に異なる非平衡成長ダイナミクスを反映しているのかについては、不明なままである。さらに、静的な指標は形態発達の動的な経路を覆い隠してしまい、類似した最終形態が異なる成長履歴から生じたのかを区別したり、成長中の過渡的な再編成を特定したりすることができない。本研究は、非平衡成長がいつ単一スケールの幾何学的制約に従い、いつそこから逸脱するかを決定するための、時間分解能を持つ幾何学的枠組みの必要性に対処するものである。

手法
著者らは「幾何学第一」のアプローチを採用し、形態を静的な結果としてではなく、時間分解能を持つ顕微鏡データから再構成された動的な観測量として扱う。

• データソース: 5種類の Bacillus subtilis 株（100, 102, 106, 108, および NCIB3610）のタイムラプス顕微鏡ビデオを、文献[34]から取得した。これらは数日間にわたり6分間隔で取得されたものである。
• 画像処理: オープンソースのPythonフレームワークであるPyPetana（OpenCVをベースとする）を用い、適応型閾値処理と形態学的演算を通じて、生のビデオをバイナリマスクへとセグメンテーションした。コロニーの境界は輪郭検出アルゴリズムを用いて抽出された。
• 幾何学的観測量: 各フレームに対して、一連の時間依存的な幾何学的記述子を算出する：
  • 面積 ($A(t)$) および 周囲長 ($P(t)$): バルクの蓄積と境界の組織化を追跡するため。
  • 形状記述子: 円形度 ($C = 4\pi A/P^2$) および コンパクトネス ($\xi = P^2/A$) を用い、全域的な対称性と境界の複雑さを定量化する。
  • フラクタル次元: 境界に対するボックスカウンティング法による有効境界フラクタル次元 ($D_b$) と、充填されたマスクに対するバルク（質量）フラクタル次元 ($D_f$) を推定する。
• スケーリング解析: 解析の中核は、面積対周囲長の時系列軌跡を構築し、スケーリング関係 $A \sim P^\alpha$ をテストすることである。指数 $\alpha$ は、成長が単一の有効な幾何学的長さスケールによって支配されているかどうかの診断指標となる（2次元のコンパクトな成長では $\alpha \approx 2$）。

主な結果
解析により、5つの株の間で、視覚的な形態と基礎となる幾何学的スケーリング挙動との間の顕著な分離が明らかになった：

1. 堅牢なコンパクト・スケーリング: 株100, 108, および NCIB3610は、視覚的に異なる形態と内部構造を示すが、全成長軌跡にわたって堅牢な単一の冪乗則関係（$A \sim P^\alpha$ で $\alpha \approx 2$）に従う。これは、視覚的な不均一性にもかかわらず、それらの成長が一貫して単一の有効な幾何学的長さスケールによって支配されていることを示している。
2. 過渡的なスケーリングの崩壊: 株106および102は、特定の形態変化に伴う単一スケール挙動からの著しい逸脱を示す：
   • 株106: 2つのスケーリング領域間の明確なクロスオーバーを示す（初期段階では $\alpha \approx 1.62$、後期段階では $\alpha \approx 1.42$）。単一指数のスケーリングが失敗する中間領域は、周囲長の局所的なピークおよび円形度の急激な低下に対応している。この期間は、過渡的な境界再編成（高度に波打った構造）と一致しているが、バルク面積の成長 ($A(t)$) は単調かつ持続している。
   • 株102: 急速な境界平滑化（アメーバ状または葉状から滑らかな形態へ）に関連する、初期の曲線的で非スケーリングな軌跡を示し、その後、後期の安定した冪乗則レジーム（$\alpha \approx 1.9$）の出現を伴う。
3. バルクと境界のデカップリング: スケーリングの崩壊を示す株において、境界に敏感な記述子（$C, \xi, D_b$）は相関した極値または急速な緩和を示す一方で、バルク質量フラクタル次元 ($D_f$) はより緩やかに飽和へと向かう。これは、大幅な境界の再編成が、バルクの空間充填効率を乱すことなく起こり得ることを示している。

主な貢献

• 時間分解能を持つ幾何学的枠組み: 本研究は、コロニーの成長を連続的な動的プロセスとして分析するための手法を確立し、持続的な成長と過渡的な幾何学的再編成を区別する。
• 視覚的および幾何学的不均一性のデカップリング: 本研究は、視覚的に異なる形態が必ずしも根本的に異なる成長ダイナミクスを意味しないこと、逆に、連続的な成長を維持しながらも、顕著な幾何学的スケーリングの崩壊が起こり得ることを示している。
• マルチスケール成長モードの特定: 単一指数の面積－周囲長スケーリングからの逸脱を特定することにより、著者らは、特定の微視的メカニズムや3次元構造の仮定を必要とせずに、成長が単一の特性長さスケールからデカップリングする（例：過渡的な境界再編成などの追加の幾何学的自由度の活性化を示す）瞬間を検出する方法を提供した。

意義と主張
本論文は、時間分解能を持つ幾何学が、生きているシステムにおける非平衡成長が単一スケールの幾何学的制約から逸脱する時期を特定するための、粗視化された枠組みとして機能することを主張している。著者らは、単一スケール挙動からの逸脱（具体的には $A \sim P^\alpha$ スケーリングの崩壊）は、成長の停止やノイズではなく、固有の動的再編成を反映していると論じている。このアプローチにより、単一の幾何学的スケールによって支配される連続的な成長と、バルクの蓄積と境界の進化がデカップリングするレジームを系統的に区別することが可能になる。研究は、グローバルな幾何学的制約（コンパクトなスケーリングなど）が持続し得る一方で、それらは動的に進化する非自己相似的な形態と共存し得るため、非平衡成長を完全に特徴付けるには複数の補完的な幾何学的記述子が必要であることを結論づけている。著者らは、内部形態（例：密度勾配や内部成長フロント）を取り入れるようにこの枠組みを拡張することが、内部の組織化と界面のダイナミクスの結合を理解するための自然な次のステップであることを示唆している。