Single-cell and Spatial Transcriptomics Reveal Cholesterol Metabolism Heterogeneity in Pancreatic Ductal Adenocarcinoma and Identify NSDHL as a Key Regulatory Gene

技術要約：単一細胞および空間トランスクリプトミクスによる膵管腺癌におけるコレステロール代謝の不均一性の解明と、主要な調節遺伝子としてのNSDHLの特定

問題提起
膵管腺癌（PDAC）は、予後不良、早期転移、および治療抵抗性を特徴とする、最も致死的な悪性腫瘍の一つであり続けている。脂質代謝、特にコレステロール代謝における代謝再プログラミングがPDACの特徴として認識されているものの、単一細胞レベルでのコレステロール代謝の具体的な不均一性や、腫瘍進行過程におけるその動的な進化については、いまだ十分に理解されていない。既存のバルクトランスクリプトーム研究では、異なる細胞サブポピュレーション間の機能的な差異を解明したり、腫瘍細胞と腫瘍微小環境（TME）の間の代謝相互作用を記述したりすることはできない。その結果、PDACにおけるコレステロール代謝の細胞学的景観や、この代謝に基づいた悪性表現型を駆動する主要な調節遺伝子の同定に関する体系的な特性評価が欠如している。

手法
本研究では、単一細胞RNAシーケンシング（scRNA-seq）、空間トランスクリプトミクス、およびバルクトランスクリプトームデータ（TCGA、GEO、ICGC）を組み合わせたマルチオミクス統合戦略を用い、PDACにおけるコレステロール代謝の景観を解析した。

1. データ処理および細胞アノテーション： scRNA-seqデータ（GSE154778、GSE155698）はSeuratを用いて処理した。悪性管細胞は、B細胞をリファレンスとしたinferCNVによるコピー数多型（CNV）解析を通じて同定した。
2. 代謝層別化： 5つの独立したアルゴリズム（AUCell、UCell、singscore、ssGSEA、AddModuleScore）を用い、各悪性細胞におけるコレステロール代謝活性をスコア化した。細胞を、低（LCM）、中間（ICM）、高（HCM）コレステロール代謝の3つのサブグループに層別化した。
3. 多次元解析：
   • 細胞間コミュニケーション： CellChatを用い、代謝サブグループとTMEの間のリガンド・受容体相互作用を推論した。
   • 軌跡およびパスウェイ解析： Monocle3を用いて擬時間軌跡解析を行った。PROGENyおよびscMetabolismを利用して、シグナル伝達経路および代謝パスウェイの活性を評価した。CytoTRACEによりステムネス（幹細胞性）を評価した。
   • 遺伝子ネットワークおよび特徴量選択： 高次元加重遺伝子共発現ネットワーク解析（hdWGCNA）により、HCM表現型に関連する共発現モジュールを特定した。XGBoostアルゴリズムを用いて特徴量の重要度をランク付けし、コアとなる調節遺伝子を特定した。
   • 空間的検証： 空間トランスクリプトミクスデータ（GSE211895）をCottrazmを用いてデコンボリューションし、様々な悪性度（hMal、iMal、nMal）を持つ組織領域への遺伝子発現のマッピングを行った。
4. 実験的検証： 最上位にランクされた遺伝子であるNSDHLの役割を、安定発現ノックダウン細胞株を用いたin vitroアッセイ（CCK-8、コロニー形成、EdU、創傷治癒、Transwell）およびヌードマウスを用いたin vivo皮下腫瘍モデルを通じて検証した。

主な貢献および結果

• 代謝の不均一性の同定： 悪性管上皮細胞は、全細胞種の中で最も高いコレステロール代謝活性を示した。この細胞集団内において、三峰性の分布が認められ、LCM、ICM、HCMのサブグループ分類が可能となった。
• HCM表現型： HCMサブグループは、以下の独特な生物学的特性を示した：
  • 強化されたコミュニケーション： HCM細胞は、LCM/ICMグループと比較して、癌関連線維芽細胞（CAF）および血管内皮細胞との間で有意に強いコミュニケーションを示した。
  • ステムネスとパスウェイ： HCM細胞は高いステムネススコア（CytheTRACE）を有し、プロ腫瘍性パスウェイ（低酸素、TNFα、NF-κB）および代謝パスウェイ（糖新生、脂肪酸合成）の協調的な活性化を示した。
  • 動的な進化： 擬時間解析により、LCM → ICM → HCMという連続的な分化軌跡が明らかになり、コレステロール代謝が静的なマーカーではなく、悪性進行の動的なドライバーであることが示された。
• NSDHLの発見： hdWGCNAおよびXGBoostによる特徴量重要度ランキングを通じて、HCM表現型を駆動する最も重要な遺伝子としてNSDHL（NAD(P)依存性ステロイド脱水素酵素様タンパク質）が特定された。
  • 発現パターン： NSDHLはPDAC組織において有意に過剰発現しており、高悪性度の空間領域に濃縮され、患者の予後（OS、DSS、PFI、DFI）と相関していた。
  • 機能的検証： PDAC細胞株におけるNSDHLのノックダウンは、in vitroにおける細胞増殖、移動、浸潤を著しく抑制した。ヌードマウスモデルにおいて、NSDHLのノックダウンは腫瘍増殖を遅延させ、腫瘍重量を減少させた。
  • メカニズムの洞察： バーチャルノックアウト解析により、NSDHLが免疫調節および上皮恒常性に関与する下流遺伝子を制御していることが示唆された。高いNSDHL発現は、ネオアンチゲン負荷および増殖特性の増加と相関していた。

意義
本研究は、バルク組織解析を超え、特定の代謝サブポピュレーションを定義することで、単一細胞解像度でのPDACにおけるコレステロール代謝の不均一性の体系的なマップを提供した。本研究は、HCMサブグループを腫瘍微小環境におけるコアなシグナルハブとして確立し、高コレステロール代謝と、強化されたステムネス、攻撃的なシグナル、および不良な臨床経過を関連付けた。

NSDHLを主要な調節遺伝子として特定したことにより、本研究は予後評価および標的治療のための新たな理論的根拠を提示している。これらの知見は、NSDHLがPDACにおける「悪性代謝スイッチ」として機能することを示唆している。本研究は、バイオインフォマティクスと機能実験を通じてNSDHLがバイオマーカーおよび治療標的としての妥当性を検証している一方で、正確な分子メカニズムの解明および、前向き臨床コホートにおけるNSDHL阻害剤の有用性を検証するためにはさらなる研究が必要であることを認めている。本研究は、膵臓癌におけるコレステロール代謝を標的とした精密治療戦略の基礎を築くものである。