Lysosomal polygenic risk score in Parkinson's disease across populations.

技術要約：多民族集団におけるパーキンソン病のリソソーム・ポリジェニック・リスク・スコア

問題提起
パーキンソン病（PD）は、遺伝的、環境的、および加齢要因が絡み合う複雑な病因を持つ、臨床的に不均一な神経変性疾患である。ゲノムワイド関連解析（GWAS）によって多数のリスク変異が特定されているものの、それらが説明できる遺伝率はごく一部である。リソソーム機能不全は、特にGBA1遺伝子に関連して、PDにおける既知の病理学的メカニズムである。しかし、広範なリソソーム経路内における共通の遺伝的変異（特定の病原性GBA1変異を除く）が、PDのリスクや臨床表現型にどのように寄与しているかは、特に多様な祖先集団において、まだ十分に解明されていない。既存のポリジェニック・リスク・スコア（PRS）は、祖先間の転移性に欠けることが多く、リソソーム機能に特化した経路別PRSモデルについても、大規模かつ多民族のコホートを用いた包括的な評価は行われていない。

方法論
著者らはリソソーム・ポリジェニック・リスク・スコア（lyso-PRS）を構築し、Global Parkinson's Genetics Program (GP2) Release 11のデータを用いてその性能を評価した。

• コホート: 本研究では、10の異なる祖先（ヨーロッパ系、アフリカ系アメリカ人、アフリカ系、アシュケナージ系ユダヤ人、ラテン系/先住民アメリカ系、中央アジア系、東アジア系、中東系、南アジア系、および複雑な混合履歴を持つ集団）にわたる32,054例の特発性PD（iPD）症例と18,188例の健康対照群（HC）を分析した。
• 遺伝データ: ジェノタイプはNeuroboosterアレイから得られ、Genoolsを用いてインピュテーションおよび品質管理が行われた。
• 表現型の定義:
  • iPD: 確立されたPD遺伝子（GBA1, LRRK2, SNCAなど）における既知の病原性変異を持つ個人を除外して定義。
  • GBA1-PD: 病原性または疑陽性のGBA1変異を持つ個人。
  • 非発症キャリア（NMC）: GBA1変異を持つが、PDの診断を受けていない個人。
• PRSの構築:
  • リソソーム遺伝子はMolecular Signatures Database (MSigDB) を用いて定義され、519個のユニークな遺伝子が得られた。残りのリソソーム経路のポリジェニックな負荷を分離するため、GBA1およびLRRK2は遺伝子セットから除外された。
  • アレル重みは、2025年のヨーロッパ系PD GWASサマリー統計量から導出された。
  • マイナーアレル頻度（MAF）が0.01未満の変異は除外した。フィルタリング後、217,053個のSNPが保持された。
  • スコアは、連鎖不平衡（LD）クランピング（$r^2 = 0.1$, 250kb）を用いたPRSice-2を使用して算出された。
• 統計解析:
  • パフォーマンスは、AUC、オッズ比（OR）、およびバランス精度（balanced accuracy）を介して評価された。
  • 臨床アウトカム（発症年齢、UPDRSスコア、認知機能）との関連は、性別、年齢、および主成分で調整した回帰モデルを用いてテストされた。
  • ヨーロッパ系コホートは、トレーニングセット（European-GP2）と独立した検証セット（European-TUEPAC）に分割された。非ヨーロッパ系コホートは、サンプルサイズの制限によりデフォルトのp値閾値を利用した。

主な結果

1. 祖先別のパフォーマンス:

   • lyso-PRSは、ヨーロッパ系集団において最も強力な予測性能を示した。European-GP2トレーニングコホートにおけるAUCは0.602 ($p = 2.02 \times 10^{-172}$)、European-TUEPAC検証コホートにおけるAUCは0.632 ($p = 5.34 \times 10^{-20}$) であった。
   • 非ヨーロッパ系集団における性能は著しく低かった（AUC範囲: 0.51–0.59）。最も低い性能はアフリカ系集団で見られた（AUC = 0.509）。
   • 非ヨーロッパ系グループにおいて、lyso-PRSスコアは一般にHCよりもiPDで高かったが、アフリカ系アメリカ人およびアフリカ系コホートでは有意性に達しなかった。

2. ヨーロッパ系集団におけるPDサブタイプの識別:

   • lyso-PRSは、GBA1-PDを健康対照群（AUC = 0.621）および非発症キャリア（AUC = 0.621）の両方から効果的に区別した。
   • しかし、このスコアはGBA1-PDと特発性PD（iPD）の間の識別能が弱かった（AUC = 0.515）。これは、スコアがGBA1変異の状態を特異的に捉えるのではなく、一般的なPDリスクを捉えていることを示唆している。

3. 臨床的関連（ヨーロッパ系コホート）:

   • 発症年齢: 高いlyso-PRSは、European-GP2コホートにおいて早期の発症に関連していた（$\beta = -0.21$, $p = 0.022$）。ただし、効果の方向性は一致しているものの、European-TUEPACコホートでは再現されなかった（$p = 0.14$）。
   • 運動症状: 高いlyso-PRSは、トレーニングコホートにおいて、日常生活における運動経験（UPDRS-II）および運動の重症度（UPDRS-III）の悪化と関連していたが、これらの知見は検証コホートでは再現されなかった。
   • 運動合併症: 高いlyso-PRSとUPDRS-IV（運動合併症）との間に、トレーニングコホート（$\beta = 0.32$, $p = 0.0088$）および検証コホート（$\beta = 0.22$, $p = 0.045$）の両方において一貫した正の関連が見られた。
   • 認知: 検証コホートにおいて、高いlyso-PRSは認知機能の低下と関連していた（$\beta = -0.17$, $p = 3.58 \times 10^{-3}$）が、トレーニングコホートでは有意な関連は見られなかった。
   • 性別: いずれのコホートにおいても、lyso-PRSと性別の間に有意な関連は見られなかった。

意義および主張
著者らは、本研究が10の集団におけるリソソーム・ポリジェニック・リスクの最大かつ初の包括的な評価であることを主張している。主な意義は、以下のことを示した点にある：

• 経路特異的リスク: （GBA1変異とは独立した）リソソーム経路における共通の遺伝的変異は、PDの感受性に寄与し、GBA1-PDと非発症キャリアを区別することでリスクを層別化できる。
• 加算的リスクモデル: これらの知見は、GBA1変異とより広範なリソソーム・ポリジェニック負荷が、独立してPDリスクを高めるという加算的な効果を示唆している。
• 臨床的関連性: リソソーム・ポリジェニック負荷は、特定の臨床的特徴、特に運動合併症（UPDRS-IV）に影響を与える可能性があり、リソソーム機能不全と治療関連の運動問題との関連性を結びつける可能性がある。
• 祖先の限界: 本研究は、PRSの性能における祖先間の重大な格差を浮き彫りにした。これは、ヨーロッパ系のGWASサマリー統計量への依存および連鎖不接パターンの違いに起因するとされる。著者らは、世界的にポリジェニック・リスク予測の精度を向上させるために、集団特異的なGWASデータの緊急の必要性を強調している。

結論として、lyso-PRSは生物学的な層別化やリソソーム標的療法のための高リスク患者の特定において有望ではあるものの、現在の有用性はヨーロッパ系集団に限定されており、臨床的関連性の検証および集団間の転移性を向上させるためには、さらなる大規模かつ多様な研究が必要である。