Integrative bioinformatics prioritises exosome-associated Epstein- Barr virus microRNA targets in gastric cancer

技術要約：統合バイオインフォマティクスによる、胃癌におけるエクソソーム関連エプスタイン・バーウイルス（EBV）マイクロRNA標的の優先順位付け

問題提起
エプスタイン・バーウイルス（EBV）関連胃癌（EBVaGC）は、顕著な免疫細胞浸潤と免疫チェックポイントリガンドの過剰発現を特徴とするが、このサブタイプにおける免疫調節および免疫逃避を駆動するメカニズムは完全には解明されていない。EBVコードのBamHI-A rightward transcript (BART) マイクロRNA（miRNA）は宿主遺伝子の発現を制御することが知られているが、どの特定のmiRNAがエクソソーム中に濃縮されているか、またそれらがどのように自然免疫シグナル伝達経路、特にToll様受容体（TLR）/NF-κB/インターフェロン調節因子（IRF）軸と相互作用するかは不明である。既存の研究は、エクソソームへの濃縮を考慮せずにmiRNA標的をマッピングしているか、あるいは実験的に検証された宿主との相互作用に結びつけることなくエクソソーム内成分をプロファイリングしているかのいずれかである。エクソソームによるmiRNA濃縮、CLASHによって検証されたmiRNA-標的結合、ネットワーク・トポロジー、および腫瘍コホートによる文脈化を統合して、候補となる調節相互作用を優先順位付けするための透明性の高い計算ワークフローが欠如している。

手法
著者らは、公開されているマルチオミクスデータセットを統合し、エクソソームに関連するEBV BART miRNA-ヒト転写産物間の相互作用を優先順位付けするための、5段階のCLASH誘導型多層計算フレームワークを開発した。

1. エクソソームmiRNA濃縮解析： SNU719（EBV陽性）およびAGS（EBV陰性）胃癌細胞とその由来エクソソームを用いた小分子RNA-seqデータ（GSE227017）を使用し、著者らは差分的濃縮解析を行った。edgeR TMM正規化を用い、エクソソームに有意に濃縮されたEBV miRNA（グループA：FDR < 0.05, log2FC > 1）および中程度に濃縮されたもの（グループB：log2FC > 0.8, 0.05 < FDR < 0.15）を特定した。
2. CLASH経路フィルタリング： 直接的なArgonaute結合を伴うmiRNA-標的結合の生化学的証拠を提供するSNU719特異的CLASHデータ（GSE147228）を解析した。相互作用は3つの生物学的レプリケート間で再現性があるものに絞り込まれ、TLR/NF-κB/IRF自然免疫経路を表す精選された38個の遺伝子セットと交差させた。
3. ネットワーク構築と摂動解析： CLASH miRNA-標的エッジとSTRINGタンパク質間相互作用（PPI）データを統合した不均一ネットワークを構築した。トポロジー指標（媒介中心性、ノード強度）およびin silico摂動解析（単一ノード除去）を用いて、miRNAsを「重み付きTLR到達可能性」（miRNAからTLR経路コンポーネントへの連結性）に対する影響に基づいてランク付けした。
4. トランスクリプトーム発現フィルタリング： SNU719 RNA-seqデータ（GSE239994）を使用して、優先された宿主標的転写産物の発現を確認し、特定されたmiRNA-標的ペアが生物学的に存在する転写産物を対象としていることを保証した。
5. 腫瘍コホートの文脈化： TCGA-STAD（The Cancer Genome Atlas-Stomach Adenocarcinoma）データを分析し、優先された経路シグネチャーまたはEBVステータスが、バルク腫瘍のTLR経路活性または患者生存率と相関するかどうかを評価した。これは、提案されたメカニズムの境界条件として機能する。

主な結果

• エクソソーム濃縮： 3つのEBV miRNA（BART1-5p, BART2-3p, BHRF1-2-3p）がエクソソムに厳密に濃縮されていた（グループA）。その他の4つ（BART13-5p, BART17-5p, BART10-3p, BART8-3p）は中程度に濃縮されていた（グループB）。特に、BART17-5pとBART17-3pは相反する濃縮パターンを示しており、これは比例的な輸出ではなく選択的な分配を示唆している。
• 再現可能な相互作用： CLASHデータとTLR経路セットの統合により、3つのSNU719レプリケートすべてに存在する11個の高信頼かつ再現可能なmiRNA-標的相互作用が特定された。厳密に濃縮されたグループAのmiRNAの中に、精選されたTLRセットとの再現可能な直接的相互作用を示すものはなかった。
• 優先された候補： ネットワーク摂動解析により、BART17-5pが支配的なハブであることが特定され、in silico除去時にモデル化されたTLR到達可能性を37.28%減少させた。次いでBART8-3p（21.72%減少）およびBART10-3p（20.27%減少）が続いた。
• 具体的な相互作用： 解析により、複数の経路レベルにおける協調的な標的化が明らかになった：
  • BART8-3pとBART17-5pはNFKB1およびNFKB2を標的とする。
  • BART17-5pとBART10-3pはIRAK1に収束する。
  • BART8-3pとBART10-3pはIRF3を標的とする。
  • BART17-5pはさらにTAB2, IKBKG, TIRAP, およびRELAを標的とする。
• 発現検証： SNU719 RNA-seqにより、優先されたmiRNA-標的ペアのうち10個が細胞内で発現している転写産物を伴うことが確認された。3つの相互作用が「ティア1（Tier 1）」ステータス（高発現の標的かつ強いCLASHスコアを持つ）に達した：BART8-3pからNFKB1、BART17-5pからNFKB2、およびBART17-5pからTAB2への相互作用である。
• TCGA-STADの文脈化： 412例のTCGA-STAD腫瘍の分析により、EBV陽性腫瘍とEBV陰性腫瘍の間で、バルクのTLR、NF-κB、IFN応答、またはサイトカイン経路のスコアに有意な差は認められなかった。同様に、EBVステータスおよびTLR経路スコアは、全生存期間を層別化しなかった。

意義と主張
本論文は、公開されているマルチオミクスデータセットから、エクソソーム関連のウイルスmiRNA-ヒト標的相互作用を優先順位付けするための、再利用可能でモジュール式の計算フレームワークを提供することを主張している。核心的な貢献は、自然免疫シグナル伝達に関連する実験可能な候補として、BART17-5p, BART8-3p, およびBART10-3pを特定したことである。

著者らは、本研究が免疫抑制、腫瘍の発現、または進行を実証したものではないことを明示している。代わりに、本研究は、将来の実験的検証（例：レポーターアッセイ、miRNA摂動、およびエクソソーム受容細胞実験）のための、ランク付けされエビデンスに基づいた候補相互作用のセットを生成するものである。

TCGA-STADのバルク解析における負の結果は、フレームワークの失敗ではなく、生物学的な区別として解釈されている。すなわち、提案されたメカニズムは、腫瘍と受容免疫細胞の間の局所的なエクソソームを介したコミュニケーションに関与しており、それがバルク腫瘍のトランスクリプトームにおいて希釈または隠蔽されている可能性があるということである。したがって、本研究は、特定されたBART miRNAを、バルク腫瘍の転写シグネチャーではなく、局所的な細胞間コミュニケーションの候補として位置づけている。