No Large Absolute Shift in the Mitochondrial Chaperonin–Import Axis in Autism Dorsolateral Prefrontal Cortex

技術要約：自閉スペクトラム症の背外側前頭前野におけるミトコンドリア・シャペロニン–インポート軸における大きな絶対的シフトの不在

問題提起
自閉スペクトラム症（ASD）における死後脳皮質トランスクリプトミクスは、ニューロン／シナプスプログラムの減少と、ミトコンドリア酸化リン酸化（OXPHOS）モジュールの減衰を特徴とするシグネチャーを確立している。しかし、既存の文献はこれらのミトコンドリアシグナルを広範に集約しており、「タンパク質品質管理およびインポート部門」（シャペロニン、マトリックスプロテアーゼ、アンフォールド・プロテイン・レスポンス（UPRmt）エフェクター、およびトランスロカーゼ）を具体的に分離して特定することは稀である。さらに、この分野は差異の「存在」をテストすることに偏っており、特定の転写物やモジュールが実質的に変化していないことを決定するための等価性テスト（equivalence testing）を採用することはほとんどない。本研究は、ASDの背外側前頭前野（DLPFC）において、ミトコンドリアのシャペロニン–インポート軸がシフトしているかどうかを検証することで、これらのギャップに対処する。具体的には、この部門が報告されているOXPHOSのダウンレギュレーションから解離している可能性があるという仮説を検証する。

方法論
本研究は、公開データセットGSE102741（ASD 13例、コントロール39例の死後脳DLPFCサンプル）を用い、2つの異なる分析レイヤーによる仮説先行型アプローチを採用した。

1. 主要な単一遺伝子解析:

   • ターゲット: HSPD1（ミトコンドリア・シャペロニンHSP60をコード）を事前に指定された主要ターゲットとした。HSPE1（共シャペロニンHSP10）およびTOMM20（インポート受容体）を二次的なプローブとした。
   • データ: 著者が公開したlog2(RPKM + 1)発現行列を分析した。
   • モデル: 事前に指定された単一の線形モデルを用い、年齢、RNA完全性指数（RIN）、性別、人種、エキソンマッピング率、および9つのトランスクリプトーム主成分（PC）を調整して、ASD–コントロール間の差を推定した。
   • 推論: HC3ヘテロスケダスティック一貫標準誤差、100,000回の反復によるフリードマン–レーン・パーミュテーション・テスト、および二標本一側テスト（TOST）を用いた。等価領域は、±0.30 log2単位（約±19%から±23%の変化）と保守的に設定した。
   • 二次的／探索的: 二次的な遺伝子にはベンジャミニ–ホフバーグ補正を適用し、探索的な対比（例：シャペロニンとインポート受容体の差）にはホルム補正を用いた。

2. モジュールレベルの再処理および解析:

   • 再処理: 生のリードをrecount3（GENCODE v26）を用いて整数カウントとして再定量し、カウントデータに対してRPKMよりも優れた平均–分散フレームワーク（TMM–voom）を適用した。
   • 遺伝子セット: 「シャペロニン–インポート」モジュール（タンパク質インポート、シャペロン、安定性、およびUPRmtをカバーするMitoCarta3.0由来の86遺伝子）と、「酸化リン酸化」モジュール（169遺伝子）を定義した。
   • テスト: 競合的（camera）および自己完結的（fry）遺伝子セットテストを適用し、モジュールレベルの等価性境界（±0.30 log2）も併せて適用した。
   • 外部チェック: データの一致性を評価するために、独立したコホート（GSE236761、n=10）を用いた限定的な再解析を行ったが、パワー不足であることも注記している。

主な結果

• 主要な単一遺伝子解析の結果:

  • HSPD1: 調整後の効果量は−0.097 log2単位（95% CI: −0.266 to 0.071; P = 0.250）であり、RPKM+1スケールで約6.5%の減少を示した。これは統計的に有意ではなかった。
  • 等価性: 90%信頼区間（−0.238 to 0.043）は、±0.30の等価性境界内に完全に収まった（TOST P = 0.010）。これにより、中程度またはそれ以上の絶対的なシフト（>19%の減少または>23%の増加と定義）を統計的に除外した。
  • 二次的な遺伝子: HSPE1はHSPD1と同様の傾向を示した（非有意）。TOMM20は名目上高かったが（P = 0.043）、多重性の補正には耐えられなかった。探索的な対比（シャペロニン minus TOMM20）は負の値を示したが、ホルム補正後には有意性に達しなかった（P = 0.072）。
  • 頑健性: 感度分析（PCの除外、pHの追加、男性への限定、leave-one-out）においても、一貫して負の非有意な推定値が得られた。

• モジュールレベルの解析結果:

  • シャペロニン–インポート・モジュール: 平均的な調整効果は−0.037 log2（90% CI: −0.262 to 0.189）であった。これは等価領域内に収まり（TOST P = 0.028）、中程度またはそれ以上の絶対的なシフトがないことが確認された。
  • 酸化リン酸化モジュール: 平均効果は−0.054 log2であったが、より広い90% CI（−0.456 to 0.348）を示し、等価性境界には収まらなかった（TOST P = 0.155）。これは、データと矛盾しない中程度の減少が可能であることを示唆している。
  • 競合的ランキング: 両モジュールとも、競合テスト（camera P < 2×10⁻⁴）においてトランスクリプトームの他部分よりもわずかに低くランク付けされたが、自己完結的テスト（*fry* *P* > 0.08）では有意性に達しなかった。
  • 外部コホート: 独立したコホート（GSE236761）では、シャペロニン軸はゼロ付近に中心があったが、等価性を確立するにはパワー不足であった。

意義および主張
本論文は、ASDのDLPFCにおけるミトコンドリア・シャペロニン–インポート軸に関する「有界な負の回答（bounded negative answer）」を主張している。その主要な貢献は、広範なOXPHOSモジュールがダウンレギュレートされていると報告されている一方で、特定のタンパク質折り畳みおよびインポート機構（HSPD1、HSPE1、および86遺伝子のモジュールによってサンプリングされたもの）は、統計的に検出可能な中程度またはそれ以上の絶対的なシフトを示さないことを示した点にある。

本研究は、経路のダウンレギュレーションは、すべてのミトコンドリア転写物がシフトしていることを意味するわけではないと論じている。折り畳みおよびインポート部門は比較的安定しており、境界内に留まっているように見えるが、エネルギー生産部門は無界（unbounded）である。これは、ASDにおけるミトコンドリアのタンパク質品質管理と酸化リン酸化の間の調節における解離を示唆している。

著者は、mRNAの等価性がタンパク質量、局在、または活性の等価性を証明するものではないことを明示しており、細胞特異的な、あるいは翻訳後調節による小さな効果が依然として存在する可能性についても言及している。本研究は、サンプルサイズが小さいこと（ASD 13例）、バルク組織を用いていること、および事前の登録が行われていないことによる制限があるが、等価性テストやパーミュテーション・チェック、および生カウントの再定量化によって分析の強固さが強化されている。結論として、本論文は、このデータセットにおいて、この特定の軸における中程度またはそれ以上の絶対的なシフトは支持されないと述べている。