A Spatial Risk Model Reveals a 25 μm Transcriptomic Event Horizon for Plaque-Associated Myeloid Activation in Alzheimer’s Disease

技術要約：空間的リスクモデルが、アルツハイマー病におけるプラーク関連マイエロイド活性化の25 μm転写学的イベント・ホライゾンを明らかにする

問題提起
空間的トランスクリプトミクス（SRT）は、アルツハイマー病（AD）の病理を高解像度でマッピングすることを可能にしたが、恒常的なミクログリアおよび血管周囲マクロファージ（PVM）が疾患関連状態へと移行する正確な空間的境界は未定義のままである。軌道推論やグラフニューラルネットワーク（GNN）といった既存の計算手法は、通常、細胞の活性化を緩やかで滑らかな連続体としてモデル化する。これらの手法は、不可逆的な細胞運命を決定づける可能性のある、急激で非線形な空間的状態遷移を捉えきれないことが多い。さらに、高解像度のSRTデータ（例：MERFISH）は本質的に疎（sparse）であり高次元であるため、次元削減を必要とするディープラーニング・フレームワークにおいては、希少な細胞集団や重要な非コードRNAが隠蔽される可能性がある。

手法
本研究では、シアトル・アルツハイマー病脳細胞アトラス（SEA-AD）の中側頭回MERFISHデータセットから得られた、計195,620個の複合ミクログリア-PVM細胞を解析した。著者らは、以下の主要コンポーネントを備えた「空間的リスクモデル」アーキテクチャを採用した：

• 状態の定義： 疾患関連ミクログリア（DAM）の標準的なマーカーであるCD74およびCTSSに基づく複合活性化スコアを生成した。完全共分散を用いた非教師あり学習による2成分ガウス混合モデル（GMM）を用い、スコア38.0において明確な活性化閾値を設定し、細胞を恒常的状態（クラス0）または疾患関連状態（クラス1）に二値化した。
• プラーク・プロキシ（代理指標）： MERFISHパネルに直接的なアミロイドβ（Aβ）染色が含まれていないため、最も活性化された500個のミクログリアを、炎症のエピセンター（震源地）を表す転写学的「プラーク・プロキシ・アンカー」として利用した。
• 空間的リスクモデリング： データリークを防ぐためにCD74とCTSSを除外した138遺伝子の特徴量セットを用いて、スパース性を考慮したXGBoost分類器を訓練した。その際、物理的な最近隣プラーク・プロキシ・アンカーへのユークリッド距離を、重要な空間的共変量として統合した。本モデルは、従来の生存分析における時間軸の代用として物理的距離を扱った。
• 不連続性解析： 活性化が連続的な勾配に従うのか、あるいは離散的な閾値を持つのかを判断するために、LOESS回帰モデルと区分的回帰（piecewise regression）モデルを比較した。構造的なブレークポイントを特定するために、5 μmから75 μmまでのノット（結節点）を反復探索する網羅的なノット・スイープ分析を実施し、赤池情報量基準（AIC）を最小化する構造的ブレークポイントを特定した。
• 解釈性： 2次元のShapley Additive Explanations（SHAP）を用い、遺伝子と空間的特徴の間の非線形な相乗作用をマッピングし、主要な制御ドライバーを特定した。
• 閾値前解析： 「警告ゾーン」（プラークから30–40 μm）における恒常的な細胞と、遠位の背景組織（>50 μm）との間で、微分発現解析（ウィルコクソン順位和検定）を行った。

主な結果

• モデルの性能： XGBoostモデルは、独立したテストセットに対してROC-AUC 0.896を達成した。これは、疾患関連状態への移行が、定義されたマーカーなしでも検出可能な、より広範な転写学的再編成に基づいていることを示している。
• 25 μmのイベント・ホライゾン： 25 μmに構造的ノットを持つ区分的回帰は、LOESSモデルのAIC 261,931.2に対し、261,601.3というAICを導出した。このΔAIC 329.9は、緩やかな連続体ではなく、離散的で急激な空間的閾値が存在することに対する強力な統計的支持を与えている。この閾値は、アンカー数の変動（250–1000細胞）に関わらず、安定していた（24.8–25.7 μm）。
• 制御ドライバー： SHAP解析により、LRRK1、L3MBTL4、PDE4Bが主要な正の予測因子であることが特定された。逆に、TACR1とGRIP2は負のSHAP値を示し、これらのダウンレギュレーションが活性化閾値を越えるための前提条件であることを示した。
• 閾値前バイオマーカー： プラークから30–40 μmの位置にある細胞（「警告ゾーン」）は、SNTB1、SLC24A2、PEX5Lの早期誘導を示しており、カルシウム恒常性、細胞骨格スキャフォールディング、および酸化ストレス管理に関する初期の細胞ストレス応答を示唆している。
• ベンチマーク： 標準的な拡散ベースの擬時間（pseudotime）アルゴリズムは、この鋭い空間的不連続性を回収できず、物理的に大きく離れた細胞（例：5 μm対2,000 μm）に対して同様の疾患スコアを割り当てた。これは、空間的軌道推論手法の限界を浮き彫りにしている。

意義および主張
著者らは、アルツハイマー病におけるプラーク関連の炎症ニッチを取り囲む、局在化した25 μmの転写学的イベント・ホライゾンを特定したと主張している。この知見は、神経炎症を均質な、あるいは緩やかに拡散する勾配とする従来の視点に挑戦するものである。むしろ、本研究は、ミクログリアの活性化とは、局所的なオリゴマーやサイトカインの濃度が生理学的許容範囲を超えたときに引き起こされる、急激で非線形なイベントであると断じている。

この空間的境界は、マイエロイド活性化のダイナミクスを理解するための定量的な枠組みを提供する。本研究は、正確な「イベント・ホライゾン」と、閾前ゾーンにおける早期警告バイオマーカー（SNTB1、SLC24A2、PEX5L）を特定することで、不可逆的な神経毒性変換が起こる前にミクログリアを安定化させるための治療介入の標的を提示している。著者らは、物理的距離を明示的にモデル化する彼らのスパース性考慮型XGBoostアーキテクチャが、AD病理に固有の鋭い空間的不連続性を捉える上で、従来の軌道推論に代わる優れた選択肢を提供すると結論付けている。

著者らが認める限界事項

• プロキシの使用： MERFISHパネルに直接的なAβ染色が含まれていないため、プラークの位置は転写学的プロキシを介して推定されており、緻密なコアプラークの物理的なサイズが数学的に含まれていない。
• 検証： 38.0という活性化閾値は数学的に導出されたものであり、明確に定義された反応性表現型を持つデータセットでの将来的な検証が必要である。
• パネルの制約： 対象となった140遺伝子パネルには、特定の血管周囲マクロファージマーカー（例：MRC1）が含まれておらず、結果はミクログリアとPVMの区別されたダイナミクスではなく、結合されたマイエロイド集団を反映している。
• 生体内検証： 本研究は、全トランスクリプトームプラットフォーム（例：Visium HD、Stereo-seq）や生体内モデルを用いた、これらの空間的なティッピングポイント（転換点）の将来的な検証を求めている。