A transport-only null model for apparent heterogeneity in diffusively dosed organoid arrays

技術要約：オルガノイド・アレイにおける見かけ上の不均一性のための輸送限定ヌルモデル

問題提起
患者由来の肝細胞癌オルガノイドは、腫瘍生物学の研究にますます利用されているが、名目上は同一の投与プロトコルに従っているにもかかわらず、オルガノイド間で大幅な表現型の変異を示すことがよくある。この変異は、生物学的不均一性や測定誤差に起因するとされることが多いが、しばしば見落とされる第三のメカニズムが存在する。それは、空間的な輸送による細胞内曝露の不平等である。共有された拡散プールにおいて、オルガノイドの位置（投与源やチャンバー境界、および隣接するオルガノイドとの相対的な位置）は、各コンパートメントに到達し進入する物質の量を決定する。したがって、空間的に不均一なデリバリーは、すべての生物学的パラメータが同一であっても、不均一な表現型を生じさせ得る。本論文は、幾何学的形状とデリバリー・プロトコルのみによってどの程度の個体間変異が生じ得るかを定量化するための、利用可能な「輸送限定ヌルモデル（transport-only null）」の欠如に対処するものである。

手法
著者らは、バルク拡散、表面キネティクス、および細胞内状態ダイナミクスを結合させた「輸送から表現型へのヌルモデル」を構築している。

1. バルク輸送と表面キネティクス： モデルは、細胞外領域を、$N$ 個の小さく十分に離れたコンパートメント（オルガノイド）を含む有界なチャンバーとして扱う。これは、バルク拡散を、遭遇依存的な吸着、可逆的な表面滞留、および生産的な内部移行と結合させている。

   • 更新理論（Renewal Theory）： モデルは、繰り返される吸着、滞留、脱離、および再吸着のサイクルを記述するために更新方程式を採用している。モデルは、「リセット（脱離した分子がその位置を忘れ、十分に混合された流体へと再進入する）」と「継続的な探索（空間的記憶が保持される）」を区別している。
   • 境界条件： 表面交換は、静的なシンク（吸収源）としてではなく、非マルコフ的な遭遇ダイナミクスを扱うために境界局所時間形式法を取り入れた、周波数依存的な反応性フィルターとしてモデル化されている。

2. 一致漸近展開による次元削減：

   • 二次元ドメインにおいて、一致漸近展開を用いることで、穿孔ドメインの問題を、オルガノイドの中心における $N \times N$ のグリーン関数系へと簡約化する。
   • この簡約化は、生産的なフラックスを、モノポールの強度および形状依存的なキャパシタンス（3Dの場合）または対数相互作用（2Dの場合）を用いて表現する。
   • この空間的簡約化の誤差は、円形ターゲットにおいて、$\varepsilon$（オルガノイド半径とチャンバーサイズの比）に対して $O(\varepsilon^2)$ であることが証明されている。

3. 細胞内状態ダイナミクス： 生産的な内部移行フラックスは、5次元の細胞状態モデル（前駆細胞、分化細胞、成熟細胞、アポトーシス/ストレス状態、および疾患スコア）を駆動する。モデルはすべてのオルガノイドにおいて生物学的パラメータが同一であることを前提としている。したがって、出力される表現型の変異は、純粋に幾何学に依存する生産的フラックスから生じる。

4. 数値検証：

   • 空間的検証： 簡約化されたモデルは、正確な半径方向のベンチマークおよび、穿孔ドメイン問題の完全な有限要素解（P2要素）に対して検証されている。検証は $N=2, 5, 10$ のターゲットに対して行われた。
   • 時間的検証： 時間領域の再構成は、数値的なラプラス反転によって検証されている。著者らは、二モーメント・ガンマ近似（two-moment Gamma surrogate）と厳密な反転を比較し、近似がマルコフ的（指数関数的）な滞留則に対しては非常に正確である一方、非マルコフ的（ヘビーテイル）な法則に対しては失敗する可能性があることを指摘している。

主要な貢献

1. 輸送限定ヌルの定式化： 本論文は、測定されたオルガノイドの中心位置と投与幾何学から、同一の生物学条件下での生産的取り込みと表現型の変異をマッピングする直接的なレシピを定式化している。これにより、観察された不均一性をテストするための、レイアウト固有のリファレンスが得られる。
2. 空間的広がりの定量化： 局所的な投与を行う2000個のランダムな10個体オルガノイド・アレイにおいて、モデルは中央値としての成熟係数変動係数（CV）を0.623、疾患スコアのCVを0.835と予測している。1因子設計監査によれば、この中央値は幾何学に応じて0.27から0.86の範囲で変化する。
3. メカニズムの解明：
   • 保存的 vs 損失的探索： 保存的なチャンバー（反射壁があり、バルク損失がない場合）では、最終的な総生産的取り込み確率は1となる。脱離はデリバリーを遅延させ取り込みを再分配するが、競合する損失チャネルが存在しない限り、最終的な総取り込み量を減少させることはない。
   • 滞留則の形状： 平均値は等しいが形状が異なる（例：ガンマ分布 vs 指数分布 vs ヘビーテイル分布）滞留則は、異なる一過性の表現型を生む。一過性は、単なる平均値ではなく、法則の形状に依存する。
   • 競合の正規化： 生の取り込み量は、競合とソースからの距離を混同してしまう。ソースの強制力を相殺し、ターゲット間の相互作用を孤立させるために、競合正規化された比率が導出されている。
     ло
   • リミックス（再混合）の識別可能性： 部分的なリミックス（リセットと継続的な探索の間の補間）は、滞留が不均一な場合にのみ、エンドポイントの取り込み分率から識別可能である。時間分解された取り込みは、滞留が均一な場合でも、一般的に情報量を持つ。
4. 設計への示唆：
   • 投与幾何学： 分散型投与（例：4つのポート）は、平均曝露量を一致させた後、単一の局所ポートと比較して、ベースラインの広がりを約5分の1に低減させる。
   • バッチサイズの限界： バッチサイズを増やしても、独立したバッチ内プロセス変異のみが平均化される。共有されたライン効果やバッチ効果が存在する場合、一般的な大規模バッチの受容確率は非退化的であり、共有効果が存在する場合、バッチサイズを大きくしても決定論的な結果に収束するとは限らない。

結果

• ベースラインの変動性： 局所ポートの場合、成熟CVの中央値は0.623である。これは、空間的に一様なソース（平均曝露量を一致させた後）では0.124に低下し、広がりが約5分の1に減少することを示している。
• 感度： 変動性は、ターゲット数、ターゲットサイズ、チャンバーのアスペクト比、および投与位置（例：コーナー vs センター）に対して高い感度を持つ。保存的な極限においては、バルク損失率にはほとんど影響されない。
• 収束性： 簡約化されたモデルは、生産的フラックスに対して $O(\varepsilon^2)$ の精度を達成する。ターゲット間の距離が約4ターゲット半径の場合、表現型のCVにおける誤差は $10^{-3}$ 未満となるが、個々のフラックス誤差は依然として高い（$0.5\% - 2\%$）。
• バッチ受容： 共有効果（ラインまたはバッチ）が存在する場合、バッチサイズが増加しても、バッチ受容の確率は0または1には収束しない。代わりに、共有効果の分布によって決定される内部限界へと収束する。

意義と主張
本論文は、校正されたデジタルツインではなく、幾何学固有のヌルモデルを提供することを主張している。その主な意義は、観察された不均一性を生物学的変異に帰する前に、測定された幾何学とデリバリー・プロトコルを定量化し、それらがどれだけの変異を説明できるかを確認するという基準（ベースライン）を確立することにある。

• モデルの予測を超える測定値が得られた場合は、輸送のみでは説明できない生物学的またはキネティックな変異が厳密に特定される。
• 予測を超えないことが、生物学的な均一性を証明するわけではない。それは単に、データが輸送に起因するノイズと区別がつかないことを示しているに過ぎない。
• 本フレームワークは、具体的な、テスト可能な設計推奨事項を提供する。すなわち、分散型投与は、見かけ上の不均一性を大幅に減少させ、バッチサイズを増やすことは、共有されたプロセス効果に対処しない限り、収穫逓減となる。

著者らは、数学的構成要素（更新理論、一致漸近展開、ディリクレ・トゥ・ノイマン簡約）は既存の文献から継承されたものであることを明示している。新規性は、それらを検証済みの「輸送から表現型へのモデル」として、肝細胞癌オルガノイド・アレイへ構成し、強制力の等価化、滞留の解釈、およびバッチ受容の限界に関する具体的な帰結を導出した点にある。