Integrated Multi-Omics and Network-Based Identification of Prognostic Cell-Cycle Hub Genes in Lung cancer

技術要約：肺癌における予後に関連する細胞周期ハブ遺伝子の統合的マルチオミクスおよびネットワークに基づく同定

問題提起
肺癌（LC）は、依然として世界における癌関連死亡の主要な原因であり、2022年には約180万人の死者が記録されている。標的療法（例：EGFR、ALK、KRAS阻害剤）や免疫療法が進展しているにもかかわらず、進行した段階での診断、腫瘍の異質性、薬剤耐性、およびグローバルに文書化された予後マーカーの欠如によって、臨床的なアウトカムはしばしば制限されている。現在の病期分類（AJCC）や臨床病理学的な限界は、正確な予後予測をしばしば不明瞭にしており、診断、治療標的、および患者の生存予測を向上させるための、新規かつ堅牢なバイオマーカーの同定が必要とされている。

方法論
本研究では、予後に関連するハブ遺伝子を同定するために、包括的な in silico マルチオミクスおよびネットワークベースのアプローチを採用した：

1. データ取得および前処理： トランスクリプトームデータを3つのGEOデータセット（GSE81089、GSE40419、GSE141569）から取得した。生のRNA-seqデータに対し、FastQCによる品質管理、Cutadaptによるアダプター除去、STARを用いたヒト参照ゲノム（GRCh38.p14）へのアライメント、およびfeatureCountsによるリードカウントを行った。
2. 発現変動解析： ComBat-seqを用いてバッチ効果を補正した。DESelQパッケージを用い、負の二項分布を適用して、差次的発現遺伝子（DEG）を特定した。遺伝子は、補正済みp値 < 0.05 かつ |log2FC| > 2 の場合に有意であると分類された。
3. 機能的濃縮解析： topGO パッケージ（weight01 アルゴリズム、フィッシャーの正確検定）を用いて、遺伝子オントロジー（GO）濃縮解析を実施した。また、MSigDB Hallmark、KEGG、および Reactome データベースにわたる経路レベルの変化を特定するために、fgsea パッケージを用いて遺伝子セット濃縮解析（GSEA）を行った。
4. ネットワーク構築およびハブ同定： STRING データベース（信頼度スコア > 0.4）を用いてタンパク質相互作用（PPI）ネットワークを構築し、Cytoscape を用いて可視化した。上位1,000個の相互作用を解析した。ハブ遺伝子は、6つのトポロジーアルゴリズム（Degree、Maximum Clique Centrality (MCC)、Maximum Neighborhood Component (MNC)、Edge Percolated Component (EPC)、Radiality、Closeness）を用いた CytoHubba プラグインを使用して特定された。これらの手法の交差部分から最終的なハブ候補が得られた。
5. マルチオミクスによる検証：
   • 遺伝学的改変： cBioPortal を通じて解析し、変異、コピー数多型（CNV）、および増幅を評価した。
   • エピジェネティックな調節： UALCAN データベース（TCGA データ）を用いて、プロモーターのメチル化レベルを調査した。
   • 免疫浸潤： ハブ遺伝子の発現と免疫細胞浸潤（B細胞、CD4+/CD8+ T細胞、好中球、マクロファージ、樹状細胞）との相関を TIMER 2.0 を用いて評価した。
   • 生存解析： GEPIA2 サーバー（カプラン＝マイヤー・プロットおよびハザード比）を用いて、予後価値を検証した。

主な結果

• 差次的発現： 解析により、6,849個の有意な DEG（アップレギュレートされたものが4,846個、ダウンレギュレートされたものが2,003個）が特定された。
• 経路の活性化： GSEA は、細胞周期に関連する経路（E2F ターゲット、G2M チェックポイント、DNA 複製、MYC ターゲット）および代謝経路（解糖系、mTORC1）の有意な活性化を明らかにした。逆に、免疫関連の経路（補体、炎症反応、TNF-α シグナル伝達、インターフェロンガンマ応答）は有意に抑制されていた。
• ハブ遺伝子の同定： 6つのトポロジー手法すべてにおいて、CDK1、CCNB1、CCNA2、TOP2A の4つのハブ遺伝子が一貫して特定された。
• 遺伝学的およびエピジェネティックなプロファイル：
  • CDK1、CCNB1、および TOP2A： 腫瘍組織においてプロモーターの低メチル化を示し、これが過剰発現と相関していた。遺伝学的改変には、ミスセンス変異、増幅、およびコピー数増加が含まれていた。
  • CCNA2： 過剰発現しているにもかかわらずプロモーターの高メチル化を示しており、他の転写または転写後調節メカニズムを介した調節を示唆している。主にミスセンス変異を示した。
• 免疫相関： 4つのハブ遺伝すべてが、免疫細胞浸潤と有意な相関を示した。特に、これらは一般的に B 細胞および CD4+/CD8+ T 細胞とは負の相関を示し、好中球とは正の相関を示しており、腫瘍免疫微小環境の調節における役割を示唆している。
• 予後価値： 生存解析により、これら4つの遺伝子の高発現が全生存期間（OS）の悪化と有意に関連していることが確認された。
  • CDK1 & CCNA2： ハザード比 (HR) = 1.9。
  • CCNB1： HR = 1.8。
  • TOP2A： HR = 1.5。
  • すべての関連において統計的有意性が認められた（Log-rank p < 0.05）。

意義および主張
著者らは、統合的な解析によって CDK1、CCNB1、CCNA2、および TOP2A が、肺癌における有望な予後バイオマーカーおよび潜在的な治療標的として特定されたと主張している。本研究は、これらの遺伝子が肺癌における細胞周期チェックポイント、染色体不安定性、および免疫微小環境のリモデリングの異常を特徴づけていると断じている。

本論文は、これらのハブ遺伝子が肺癌の進行に関する詳細な分子学的理解を提供し、革新的な診断技術および治療戦略の開発を支援する可能性があると結論付けている。ただし、著者らは、in silico の知見は堅牢であるものの、実際の患者コホートにおける生物学的優先順位と治療的可能性を確認するためには、さらなる実験的および臨床的な検証が必要であると控えめに述べている。本研究は特定の新しい薬物製剤を提案するものではないが、肺癌の増殖と転移を抑制するために、これらの指標を標的とする阻害剤の検討が可能であることを示唆している。