Integrative Mendelian Randomization and Single-Cell Transcriptomic Analyses Reveal Ferroptosis-Related Protective Genes and Tumor Microenvironmental Remodeling in Colorectal Cancer

技術要約：大腸がんにおけるメンデルランダム化解析と単一細胞トランスクリプトーム解析の統合的研究

問題提起
大腸がん（CRC）は、再発、転移、および治療反応における顕著な異質性を特徴とする、がん関連死亡の主要な原因であり続けている。フェロトーシス（鉄依存的な、脂質過酸化によって駆動される制御された細胞死の一形態）は、腫瘍の進行および微小環境の再構築における重要な因子として認識されつつあるが、既存の大腸がん研究は、主にトランスクリプトームのシグネチャーや予後モデルに依存してきた。フェロトーシス関連制御因子の遺伝学的妥当性、およびその細胞型特異的な微小環境的特徴を定義する上で、大きなギャップが存在する。現在のバルク・トランスクリプトーム解析の手法では、腫瘍発生のダウンストリームの結果としての変化と、因果関係を持つ寄与因子を容易に区別することができない。また、これらの制御因子が機能する特定の細胞コンパートメントも十分に定義されていない。

方法論
本研究では、バルク・トランスクリプトーム解析、メンデルランダム化（MR）、および単一細胞RNAシーケンシング（scRNA-seq）を組み合わせたマルチオミクス統合フレームワークを用い、フェロトーシスに関連する保護的な遺伝子を同定および特性評価した。

1. データ取得および発現差異解析： GEOデータベース（発見コホートとしてGSE44861、検証コホートとしてGSE44076およびGSE25070）から得られたバルク・トランスクリプトームデータを解析し、CRC組織と正常組織間の発現変動遺伝子（DEG）を特定した。これらのDEGをFerrDbデータベースから抽出した564個のフェロトーシス関連遺伝子（FRG）のキュレーションリストと交差させ、70個の差発現フェロトーシス関連遺伝子を得た。
2. タンパク質相互作用（PPI）： STRINGデータベースを用いてPPIネットワークを構築し、70個の候補遺伝子の中からハブ遺伝子を特定した。
3. メンデルランダム化（MR）： eQTLGenコンソーシアムのcis-eQTLデータを操作変数とし、CRCのGWAS要約統計量（ebi-a-GCST90018808）をアウトカムとして、二サンプルMRフレームワークを適用した。このアプローチにより、交絡因子や逆の因果関係を厳密に制御しつつ、CRCリスクと遺伝的に予測された因果関係を持つ遺伝子を優先的に抽出した。頑健性を確保するため、感度分析（不均一性、水平多面発現性、leave-one-out分析）を実施した。
4. 診断モデリング： 最優先された遺伝子を用いて診断用ノモグラムを構築し、ROC曲線、キャリブレーションプロット、および決定曲線分析（DCA）を通じて評価した。
5. 機能および免疫解析： 遺伝子セット濃縮解析（GSEA）、免疫浸潤解析（ssGSEA）、および競合内在性RNA（ceRNA）ネットワークの構築を行い、生物学的経路および調節メカニズムを探索した。
6. 単一細胞トランスクリプトミクス： 53,253個の細胞を含むscRNA-seqデータ（GSE132465）を解析し、細胞型特異的な発現パターンを定義した。血管内皮細胞に対しては擬時間軌跡解析（Monocle3）を適用し、細胞間コミュニケーション（CellChat）を用いて微小環境の再構築を推定した。

主な貢献および結果

• 保護的遺伝子の同定： 本研究は、AKR1C2、MAPK3、MAP3K11、およびCDKN1Aを、CRCリスクの低下と遺伝的に関連する4つのフェロトーシス関連遺伝子として特定した。MR解析により、これら遺伝子の遺伝的に予測された高発現が、CRCオッズの有意な低下と関連していることが示された（例：AKR1C2 OR=0.876、CDKN1A OR=0.852）。これらの関連性は、水平多面発現性の証拠がなく、感度分析においても頑健であった。
• 検証および診断的有用性： これら4つの遺伝子は、発見コホートおよび2つの独立した検証コホートにおいて一貫してCRC組織で発現低下していた。4遺伝子による診断用ノモグラムは、良好な識別能（AUC = 0.863）、較正能、および個別の遺伝子と比較した臨床的な純利益を示した。
• 機能的濃縮： GSEAにより、これらの遺伝子が細胞周期制御（E2Fターゲット、G2Mチェックポイント）、MYC駆動型増殖、炎症シグナル（TNFα/NF-κB）、および上皮間質転換（EMT）に関連していることが明らかになった。特に、これらの遺伝子は増殖経路に対して負の濃縮を示しており、制御不能な細胞増殖を抑制する役割を示唆している。
• 免疫微小環境の再構築： 免疫浸潤解析により、CRCにおける広範な再構築が明らかになった。これは樹状細胞、NK細胞、およびインターフェロン応答におけるシグネチャーの変化として現れた。これら4つの候補遺伝子は特定の免疫シグネチャーと有意な相関を示しており、特にMAPK3とMAP3K11はより広範な免疫関連性を示し、免疫状態の再プログラミングにおける役割を示唆した。また、ceRNAネットワークは、hsa-miR-199a-5pとAKR1C2/MAP3K11を介した調節軸を予測した。
• 単一細胞解像度および血管内皮の再構築： scRNA-seq解析により、細胞型特異的な異常が明らかになった。
  • AKR1C2は免疫集団では希薄であったが、線維芽細胞および血管内皮細胞において変化が見られた。
  • MAPK3およびMAP3K11は、特に上皮および血管内皮コンパートメントにおいて広範な発現を示した。
  • CDKN1Aは、最も広範な細胞型依存的な異常を示した。
  • 血管内皮の軌跡： 擬時間解析は、正常に富食された初期の血管内皮状態から、腫瘍に富食された後期の血管内皮状態への遷移を示した。CDKN1Aの発現は、腫瘍由来の後期血管内皮細胞において増加しており、これらの遺伝子が腫瘍関連の血管内皮再構築に結びついていることを示した。
  • 細胞間コミュニケーション： 腫瘍組織では、特に上皮、ストロマ、免疫、および血管内皮細胞間の強化されたクロストークを伴う、より密なコミュニケーションネットワークが観察された。

意義および主張
著者らは、本統合的研究が、AKR1C2、MAPK3、MAP3K11、およびCDKN1Aを大腸がんにおけるフェロトーシス関連の保護的遺伝子として特定することに成功したと主張している。本知見の意義は以下の点にある：

1. 因果推論： 相関関係を超え、MRを用いることで、これら遺伝子の発現とCRCリスクの低下との間の遺伝的結合を確立した。これは、これらの遺伝子の発現低下が、二次的な影響ではなく、生物学的に意味のある「保護プログラムの喪失」であることを示唆している。
2. メカニズム的洞察： これらの遺伝子は、フェロトーシス制御、細胞周期制御、炎症シグナル、および血管内皮再構築の交差点に位置付けられている。本研究は、これらの遺伝子が単なる孤立したフェロトーシスマーカーとしてではなく、酸化ストレス、増殖抑制、および微小環境のクロストークを繋ぐ、より広範な調節軸の一部として機能することを提案している。
3. 細胞コンテキスト： 単一細胞解析を用いることで、これらの遺伝子が異常を示す特定の細胞コンパートメント（上皮、免疫、ストロマ、および血管内皮）を定義し、CRCの進行における血管内皮再構築の役割を浮き彫りにした。
4. トランスレーショナルな可能性： この4遺伝子パネルは、CRCのリスク層別化および腫瘍微小環境の特性評価のための、生物学的に裏付けられたシグネチャーを提供する。

著者らは、自らの知見が公開データセットおよび計算的推論に基づいていることを認め、控えめな姿勢を維持している。彼らは、ceRNAの仮説や特定のメカニズム（例：hsa-miR-199a-5p軸）については、フェロトーシス感受性および血管内皮再構築の調節における役割を確定するために、独立したコホートおよび機能的アッセイによる実験的検証が必要であると述べている。