Stable but confounded: brain region, not diagnosis, drives a validation-passing molecular subtype in postmortem psychiatric brain

技術要約：「安定しているが、混同されている：死後脳における検証を通過する分子サブタイプを駆動するのは、診断ではなく脳領域である」

問題提起
脳のトランスクリプトームの教師なし分子サブタイピングは、精神疾患（自閉スペクトラム症、統合失調症など）の異質性を解明するためにますます利用されている。しかし、クラスタリング・アルゴリズムは、その分割が再現可能な生物学的実態を反映しているのか、あるいは単なる技術的・解剖学的なアーティファクト（偽像）なのかに関わらず、常に何らかの分割を返してしまう。本論文は、死後脳を用いた精神医学研究における3つの深刻な危険性を特定している：

1. 小規模なサンプルサイズ（数十から数百程度）による統計的パワーの限界。
2. 被験者の重複：名目上は独立しているとされるコホート間での被験者の重なりが、偽の再現性の主張を導く。
3. 解剖学的優位性：脳領域によるトランスクリプトームの変異は、診断による変異を大幅に上回る。その結果、複数の領域をプールしたデータをクラスタリングすると、疾患のサブタイプではなく、解剖学的構造が回収されてしまう。

著者らは、現在の検証慣行が、「安定した」クラスター（再現可能な構造）と「妥当な」クラスター（疾患に関連する生物学的妥当性を持つもの）を区別することに失敗していることが多いと主張している。これは、脳領域のような交絡因子が明示的にテストされていない場合に特に顕著となる。

手法
著者らは、確立されたコントロールを組み込み、明示的な交絡テストを追加した「検証優先（validation-first）」の規律を用いて、3つの公開された死後脳コホート（計384サンプル）を再解析した。データセットは以下の通り：

• 自閉スペクトラム症 (ASD): GSE28521（マイクロアレイ、前頭葉/側頭葉/小脳）および GSE64018（RNA-seq、上側頭回）。
• 統合失調症、双極性障害、主要うつ病: GSE80655（RNA-seq、前帯状回、背外側前頭前野（DLPFC）、側坐核）。

解析パイプライン:

1. 特徴量削減: 遺伝子発現を、シングルサンプル濃縮解析（ssGSEAまたは遺伝子ごとの平均zスコア）によるパスウェイレベルのスコアに集約した。
2. クラスタリング: ベイズ情報量基準（BIC）によるコンポーネント選択を用いたガウス混合モデル（GMM）を使用してサンプルをクラスタリングした。
3. 検証コントロール: すべての候補となる分割に対し、以下の4つのゲートを適用した：
   • ラベル置換ヌルモデル: 置換されたラベルからクラスターが回収可能かどうかをテスト。
   • ランダム遺伝子セット・ヌルモデル: 分割が精査されたパスウェイに特異的なものか、あるいはランダムな遺伝子セットから生じたものかをテスト。
   • ブートストラップ安定性: 元の割り当てと再サンプリングされた割り当ての間の調整ランド指数（ARI）を測定（閾値 $\ge$ 0.80）。
   • 祖先独立性: 祖先に関連する主成分との関連性を確認。
4. 交絡テスト: 分割と、脳領域、ドナーの重複、および臨床診断との間の明示的な分割表による関連性（$\chi^2$ 検定、Cramér's V）を実施した。
5. 非循環的解析: 「二重取り（double dipping）」を避けるため、疾患シグナルは、クラスタリング自体から派生した特徴量ではなく、既知の細胞型マーカー遺伝子（例：PVALB, SST, GAD1）を用いた診断に基づくケース・コントロール対比を通じて評価した。

主な結果

1. 自閉症の「GABA作動性」サブタイプの失敗:
   GABA作動性パスウェイのダウンレギュレーションを特徴とする推定ASDサブグループは、ブートストラップ安定性検証（ARI = 0.71、閾値 0.80 未満）に失敗した。さらに、第2のコホート（GSE64018）における明白な再現性は、2つのコホート間で24名中10名の被験者が重複していたために無効化された。つまり、その一致は独立した再現ではなく、同一の脳を用いたクロスプラットフォームでの裏付けであった。

2. 「検証を通過した」統合失調症の分割は、領域分類器である:
   統合失調症コホート（GSE80655）において、$k=3$ の分割はすべての統計的ゲートを通過した（ブートストラップ ARI = 0.923、ランダム遺伝子セットおよび祖先チェックを通過）。しかし、交絡テストにより、この分割は診断ではなく、完全に脳領域によって駆動されていることが明らかになった：

• サブタイプ0と2は、主に側坐核（nAcc）であった。
• サブタイプ1は、排他的に皮質（前帯状回およびDLPFC）であった。
• サブタイプと領域の関連性は圧倒的であり（$p = 4.3 \times 10^{-26}$, Cramér's V = 0.67）、一方で診断との関連性は有意ではなかった（$p = 0.41$）。
• 解析を単一の領域（DLPFC）に限定すると、安定性が検証閾値を下回った（ARI = 0.68）。これにより、「安定性」が領域軸によるアーティファクトであったことが確認された。

3. アーティファクトの系統的な性質:
   この領域駆動型のアーティファクトは、コホート内のすべての診断（双極性、主要うつ病）において系統的であり、クラスター数（$k$）にも依存しないものであった。4つの診断をプールしたコホートは、最も強い領域構造を示したが、最も弱い安定性を示した。

4. 非循環的解析による既知の生物学の回収:
   クラスタリングの枠組みは疾患シグナルを見逃したが、従来の非循環的なケース・コントロール解析は、確立されたインターニューロンの生物学を回収した：

• ASD: 皮質におけるパルブアルブミン（PVALB）マーカーの有意な減少。
• 統合失調症: 皮質におけるソマトスタチン（SST）およびGABA合成酵素（GAD1/GAD2）の有意な減少。
• 二重解離: 有意な逆方向の相互作用により、PVALBの減少はASDでより顕著であり、SSTの減少は統合失調症でより顕著であることが確認された。

5. ベンチマークの失敗:
   随伴するキャリブレーション・ベンチマーク（$R^2 = 0.889$ と主張する適応閾値モデル）は再現に失敗した。著者らはこれを指標の欠陥に起因すると突き止めた。調整ランド指数（ARI）が、退化した入力（単一の真のクラスターを持つデータセット）に対して誤った値 1.0 を返し、相関を人工的に膨張させていた。

意義および主張
本論文は、新しい分子サブタイプを発見することを目的としているのではない。むしろ、統計的な安定性チェックを厳格に通過したとしても、その分割が既知の交絡因子（脳領域）によって完全に駆動され得ることを示す「実例研究」を提供している。

著者らは以下のように結論付けている：

• 安定性 $\neq$ 妥当性: 高いブートストラップ安定性は再現性を保証するが、その分割が意図した生物学的シグナル（疾患か解剖学か）を反映していることを保証するものではない。
• 交絡テストは必須である: 複数の領域やコホートをプールした解析を行う場合、分割を脳領域やその他の交絡因子に対してテストすることは、検証プロセスにおけるルーチンの一部としなければならない。
• 発見よりも規律を: 小規模で、複数の領域をプールしたコホートにおいては、探索的な発見は失敗しやすい。著者らは、特定のチェックリストを提唱している：ヌルコントロールの報告、交絡因子に対するテスト、被験者が重複しない再現性の確保、循環的な推論の回避、そしてコントロールに失敗した分割を（隠さず）公表することである。

本研究は、分子サブタイピングが成功する場合（領域制御に成功したBowenら、2019年を例として挙げる）はあるものの、ここで記録された失敗モード（解剖学的優位性と未検討のパイプラインのデフォルトに起因するもの）が、現在の精神医学ゲノミクスにおける重大な落とし穴であることを強調している。