単一細胞および空間トランスクリプトミクスデータを用いたリガンド動態の組織レベルPDEモデルの較正

治癒中の創部で細胞はどのように情報をやり取りするか

皮膚が損傷を受けると、無数の細胞が協調して損傷を修復し、瘢痕組織を過剰に残さないように働きます。これを実現するのがリガンドと呼ばれる化学的メッセージの送受信で、これらは組織中に拡散して近傍の細胞に応答を引き起こします。本研究は、単一細胞および空間解像度の現代的な遺伝子測定を、これらの化学信号の数理モデルを較正するための強力な新しい手法に変換する方法を示しており、正常な治癒や瘢痕関連疾患の理解を深める手助けをします。

なぜ化学的メッセージが組織で重要か

組織内の細胞はめったに単独で行動しません。細胞分裂や移動、休止などのタイミングを伝える信号を絶えず交換しています。多くのこうした信号は、ある細胞が放出し別の細胞の受容体に捕捉される短距離の化学物質です。これらのメッセージは創傷治癒、炎症、線維化やがんといった疾患において中心的役割を担います。研究者は長年、これら化学物質の生成、拡散、分解を記述する反応拡散モデルと呼ばれる数式を用いてきました。しかし、実際の組織を反映するような現実的なパラメータ値（生成速度や結合強度など）をどうやって決めるかは長年の課題でした。多くの既存値は理想化された実験条件に基づいており、生体内で起きていることを正確に反映しないことがあります。

生きた組織を新たな視点で見る

近年、二つの実験技術が組織観察のあり方を変えました。単一細胞RNAシーケンシングは、何千もの個々の細胞でどの遺伝子が活性化しているかを測定し、細胞の同定やシグナル伝達能力を明らかにします。一方、空間トランスクリプトミクスは、組織切片上に格子状に配置された多数の小さなスポットごとに遺伝子発現を測定し、細胞がどこに存在するかという空間情報を保ちます。本研究では、著者たちはこれら両方のデータを、創傷後30日目に採取されたヒト皮膚サンプルから組み合わせています。この時期は治癒の「再構築（リモデリング）」段階に当たり、彼らは瘢痕形成と組織修復に影響を与えることが知られる主要なシグナル分子、トランスフォーミング成長因子ベータ（TGFβ）の三つの異なる形態に着目しています。

遺伝子からモデルパラメータへ至る段階的パイプライン

研究者らは、これら豊富な遺伝子測定をTGFβの移動と作用を記述する組織スケールのモデルに結びつける計算パイプラインを構築しました。まず既存のバイオインフォマティクスツールを用いて、主要な細胞型やサブタイプ（例えば三つの異なる線維芽細胞クラス、マクロファージ、内皮細胞など）を同定し、組織格子上の各スポットに何個ずつ存在するかを推定します。次に、リガンドを出す細胞での遺伝子発現と受容体を持つ受け手側の細胞での遺伝子発現を組み合わせて、各細胞型間がTGFβを介してどれほど強くやり取りするかを推定します。こうして推定された「相互作用強度」が、数理モデルが再現すべき実験的なターゲットになります。

モデルに組織を学習させる

次に、チームはTGFβの組織内拡散を、空間トランスクリプトミクスの配置を反映した格子上で解く偏微分方程式群で表現します。ある一組のパラメータ（TGFβの拡散速度、減衰速度、各細胞型の生成や吸収量など）を選ぶと、モデルは三種類のTGFβの濃度場と細胞型間に生じる相互作用強度を予測します。これらのパラメータを調整するために、著者らは三段階の較正戦略を用います。まず近似ベイズ法で過去の実験知見に基づいた妥当なパラメータ空間を広く探索し、遺伝子データから導出した相互作用強度とよく相関する組合せを徐々に絞り込みます。次に、それら有望な候補を勾配ベースの最適化手法に入力して、より一致するようにパラメータを微調整します。このハイブリッド手法により、モデル予測とデータ由来の相互作用強度との間でほぼ完全な線形一致（相関係数0.99）を達成しました。

較正されたモデルが示す治癒の実像

よく一致したモデルを手に、研究者らは推定されたパラメータが生物学について何を示すかを探りました。遺伝子発現データから同定された三つの線維芽細胞サブタイプは機能的に異なることが確認され、一つは強い炎症プロファイルを示し、別の一つは創傷を収縮させるのに寄与するミオファイブロブラスト様の特徴を持っていました。較正された生成率は、リモデリング期には炎症が収まりつつあることと一致して、炎症性の線維芽細胞が最も少ない量のTGFβを生成していることを示唆します。意外なことに、乳頭層（パピラリー）線維芽細胞として知られるサブタイプは、瘢痕を抑え再生的結果に関連することの多いTGFβ3を比較的多く生成しているように見えました。モデルが予測するTGFβの濃度マップは、線維芽細胞やマクロファージが集まる領域と整合しており、瘢痕形成の進み方を左右する活発なクロストークのホットスポットを示唆しています。

将来の組織と疾患に適用可能な柔軟な枠組み

この創傷治癒の具体例を越えて、本研究の主な貢献は、現代の遺伝子データを良く較正された組織レベルモデルに変換する汎用的なパイプラインです。単一細胞および空間トランスクリプトミクスを高度な統計・最適化手法と組み合わせることで、リガンドの移動速度、結合の強さ、どの細胞が信号に最も寄与しているかをより厳密に推定できます。この枠組みは他の組織や他のシグナル分子、線維化やがんのような疾患状態にも適応可能であり、分子レベルの複雑なスナップショットを組織全体の挙動を予測するモデルへと翻訳する助けとなるでしょう。

引用: Daher, A., Trucu, D. & Eftimie, R. Calibrating tissue level PDE models of ligand dynamics using single cell and spatial transcriptomics data. npj Syst Biol Appl 12, 44 (2026). https://doi.org/10.1038/s41540-026-00657-8

キーワード: 創傷治癒, 細胞シグナル伝達, 空間トランスクリプトミクス, TGFベータ, 数理モデリング