Copper metabolism-related genes SCO2 and PARK7 in pulmonary tuberculosis: diagnostic value, immune infiltration and single-cell virtual knockout analysis with clinical validation

技術要約：肺結核における銅代謝関連遺伝子SCO2およびPARK7

問題提起
肺結核（PTB）は、依然として単一の感染性因子による死因の主要なものとなっている。現在の痰を用いた微生物学的検査への依存は、乏血性（paucibacillary）疾患における感度の低さに課題があり、免疫学的アッセイは活動性感染と潜在性感染の区別に苦慮している。したがって、活動性疾患の病態生理学的状態を反映する、宿主由来の血液バイオマーカーに対する未充足のニーズが存在する。さらに、銅は宿主とマイコバクテリウム・ツベルクローシス（Mycobacterium tuberculosis）の相互作用において中心的であることが知られているが（マクロファージに対する殺菌的な武器であると同時に、細菌の解毒標的でもある）、活動性PTB患者の末梢血における銅代謝関連遺伝子の具体的な発現プロファイルや、全身性の免疫不全との関係性は十分に定義されていない。

方法論
本研究では、多段階のバイオインフォマティクスおよび臨床検証パイプラインを用いた：

1. データ取得および前処理:

   • トレーニングセット: GEOデータベースから取得したGSE28623（活動性PTB患者46名、健康対照群37名）。
   • 検証セット: GSE83456（活動性PTB患者45名、健康対照群61名）。
   • シングルセルデータ: in silico ノックアウト解析のためのGSE281148（結核性肉芽腫の穿刺吸引物）。
   • 遺伝子セット: GeneCardsから取得した2,513個の銅代謝関連遺伝子。

2. 遺伝子のスクリーニングおよび選択:

   • 発現変動解析: トレーニングセットにおいて、2,240個の変動発現遺伝子（DEGs）を特定した（調整済み P < 0.05, |log2FC| > 0.5）。
   • WGCNA: 加重遺伝子共発現ネットワークを構築し、「blue」モジュールがPTBと最も強く相関していることを特定した（r = 0.85）。
   • 交差解析: DEGs、WGCNAのblueモジュール、および銅代謝遺伝子セットの三方向の交差により、15個の候補遺伝子を得た。
   • ネットワーク解析: タンパク質相互作用（PPI）ネットワークおよびMCODE解析により、8遺伝子のコアモジュールを特定した。
   • 特徴量選択: 10分割交差検証を用いたLASSO回帰により、4つの特徴量遺伝子（PIK3CA, BECN1, SCO2, PARK7）を選択した。外部検証セット（GSE83456）において、PIK3CAとBECN1は有意な診断能を示さなかったが、SCO2とPARK7のみが有意な診断性能を維持した。

3. 臨床検証:

   • 成都郫都区人民病院における、治療未経験のPTB患者20名と健康対照群20名の末梢血を用いて、リアルタイムPCR（qRT-PCR）を実施した。
   • 診断性能は、ROC曲線、ノモグラム、および決定曲線分析（DCA）を用いて評価した。

4. 機能および免疫解析:

   • 免疫浸潤: ssGSEAを用いて免疫細胞の存在量を定量化し、遺伝子発現と免疫細胞の相関を分析した。
   • 薬剤予測: NetworkAnalyst 3.0を用い、遺伝子-化合物関連ネットワークを構築した。
   • 仮想ノックアウト: シングルセルデータセット（GSE281148）に対してscTenifoldKnk法を適用し、SCO2およびPARK7の機能喪失をシミュレートし、下流の制御ネットワークの摂動を解析した。

主な結果

• 遺伝子発現: SCO2およびPARK7はいずれも活動性PTBにおいて有意にアップレギュレートされていた。
  • 外部検証（GSE83456）: SCO2（AUC 0.952）およびPARK7（AUC 0.755）は強い診断能を示した。
  • 臨床検証: qRT-PCRにより、臨床検体におけるアップレギュレーションが確認された（PARK7: 約1.38倍, P < 0.001; SCO2: 約1.82倍, P = 0.002）。
  • 結合モデル: SCO2とPARK7を組み合わせたノモグラムは、外部検証セットにおいてAUC 0.960を達成した。
• 免疫相関:
  • 両遺伝子は自然免疫細胞（マクロファージ、好中球）と正の相関を示した。
  • 両遺伝子は防御的な獲得免疫細胞（T細胞、NK細胞、セントラルメモリーT細胞）と負の相関を示した。
  • SCO2は特にTh2および細胞傷害性細胞と負の相関を示し、PARK7はNK細胞と負の相関を示した。
• 仮想ノックアウト・シミュレーション:
  • PARK7ノックアウト: 摂動は酸化ストレス、小胞体ストレス、およびアポトーシスモジュール（例：BAX, PRDX5, p53経路）に集中した。これは、PARK7が酸化損傷を緩衝し、細胞生存を調節する機能を担っていることを示唆している。
  • SCO2ノックアウト: 摂動はインターフェロンシグナル、Th1/Th17分化、および獲得免疫経路（例：STAT1, IFNG, CXCL10）に集中した。これは、SCO2が抗結核免疫調節ネットワークに不可欠であることを示唆している。
• 薬剤/化合物関連:
  • PARK7は酸化ストレス剤および重金属に関連していた。
  • SCO2は銅代謝物質およびミトコンドリア毒性物質に関連していた。
  • 著者らは、これらの関連は直接的な治療候補ではなく、毒性ストレスや代謝経路に関するメカニズム的な手がかりを提供するものであると述べている。

意義および主張
本研究は、SCO2およびPARK7を活動性PTBにおけるコアとなる銅代謝関連マーカーとして特定したと主張している。著者らは、これらの遺伝子が末梢血において異常発現しており、銅の不均衡および免疫不全に関連していると仮定している。

• 診断的価値: 本研究は、SCO2およびPARK7を、中程度から高い精度を持つ候補診断バイオマーカーとして提示しており、PTBの補助的鑑別診断のためのツールとなる可能性を示している。
• 病理学的洞察: 本知見は、ミトコンドリアの銅代謝（SCO2）および酸化ストレス調節（PARK7）を通じて、銅の恒常性不全が結核の進行に関与していることを示唆している。
• 治療的可能性: 本研究は、これら遺伝子を宿主指向型療法（HDT）の候補ターゲットとして提案しており、特にミトコンドリア機能および免疫調節における役割に注目している。
• 限界事項: 著者らは、臨床サンプル数が少ないこと（20ペア）、単一施設での設計であること、群間での年齢の不均衡があること、および直接的なin vitroまたはin vivoの介入的証拠ではなくin silicoシミュレーションに依存していることなどの限界を明示的に述べている。彼らは、仮想ノックアウトの結果は予測的な仮説であり、さらなる実験的検証が必要であることを強調している。

結論として、本論文は銅の恒常性不全を結核の病理プロセスに結びつけ、診断のための候補分子、および銅代謝を標的とした宿主指向型介入への理論的な手がかりを提供しているが、メカニズムを検証するためのより大規模な多施設共同研究および基礎研究の必要性を認めている。