Seven replicated genomic associations of myalgic encephalomyelitis/chronic fatigue syndrome: a biobank study

問題の背景
筋痛性脳脊髄炎／慢性疲労症候群（ME/CFS）は、女性に偏った発症傾向があり、診断用のバイオマーカーや効果的な治療法、そして十分に解明されていない病因を伴う、衰弱性の疾患である。これまでの遺伝学的研究は、小規模で異質なコホートや、単一の診断コードへの依存によって阻害されており、再現性のある知見の欠如を招いてきた。これまで、ME/CFSに関するゲノムリスク遺伝子座が独立したコホート間で再現されたことはなかった。さらに、社会経済的地位（SES）や遺伝子と性の相互作用がリスクを修飾する役割については、先行研究における相反する報告や方法論的な限界により、議論が分かれている。

方法論
著者らは、UK Biobank（UKB）およびAll of Us（AoU）リサーチプログラムのデータを利用して、ゲノムワイド関連解析（GWAS）を実施した。

• フェノタイピング戦略： 単一の診断コードに頼るのではなく、症例（ケース）と対照（コントロール）を定義するために、マルチエビデンス・フェノタイピング戦略を採用した。症例には、自己申告による臨床診断、不良または可の全体的な健康状態評価、およびME/CFSに関する疼痛質問票への特定の肯定的な回答が求められた。対照には、良好または優の健康状態評価、および医療記録や調査におけるME/CFSの証拠がないことが求められた。
• コホート： 以下の3つの非重複集団が確立された：
  1. ディスカバリー（UKB1）： 症例1,268例、対照113,132例。
  2. リプリケーション1（UKB2）： 症例319例、対照28,471例。
  3. リプリケーション2（AoU）： 症例371例、対照33,140例。
• 統計的枠組み： 本研究では、ターゲット学習に基づく半パラメトリックな二重に頑健な推定ワークフローであるTarGene（Targeted Genomic Estimation）を利用した。このアプローチは、従来のパラメトリックなGWAS手法の限界に対処し、モデルの誤設定や集団層別化によるバイアスを最小限に抑えつつ、検出力を最大化するために選択された。
• 解析パイプライン：
  • ディスカバリー： UKB1における670,243個のジェノタイプ化されたSNPに対してGWASを実施し、偽発見率（FDR）を5%に制御した。
  • リプリケーション： ディスカバリー段階でFDR有意であった176の変異を、2つのリプリケーション・コホートでテストした。
  • ファインマッピングおよび共局在化（Colocalisation）： 再現された遺伝子座について、SuSiEを用いてファインマッピングを行い、因果変異の信頼集合（credible sets）を生成した。これらをGTEx v10データベースの発現定量形質遺伝子座（eQTL）と共局在化した。
  • 相互作用解析： TarGeneを用いて、遺伝子と性、および遺伝子と剥奪度（Townsend剥奪指数および複数の剥奪指数を使用）の相互作用をテストした。
  • 検証： 厳格な臨床基準（カナディアン・コンセンサスおよびIoM 2015）を用い、労作後の倦怠感（PEM）を必須とするDecodeMEコホート（2,143例）においてルックアップを実施した。

主な結果

• ディスカバリー： 解析の結果、UKB1のディスカバリー・コホートにおいて、ME/CFSのリスクと有意に関連する176の変異（FDR < 5%）が特定された。
• リプリケーション： 7つのゲノムリスク遺伝子座が、2つの非重複リプリケーション・コホート（UKB2またはAoU）の少なくとも一方において再現された。しかし、多重比較補正後、3つのコホートすべてで有意性に達した単一の変異は存在しなかった。
  • 再現された7つの変異は以下の通りである：rs115186419（ETV5近傍）、rs73175505（CSMD1近傍）、rs261902（BICD1内）、rs117553493（CLYBL近傍）、rs72741654（RORA近傍）、rs74963073（GRIN2A内）、およびrs76847656（SBK1近傍）。
• ファインマッピングおよび共局在化：
  • 第13染色体上の遺伝子座において、ファインマッピングにより、被蓋（被蓋体／基底核）におけるCLYBLの発現低下と共局在する信頼集合が解明された。リスクアレルは、CLYBLの低い発現と関連している（事後確率 H4 = 0.994）。
  • 特筆すべきことに、既知の機能喪失変異であるCLYBL Arg259 ストップゲイン変異は、ME/CFSのリスクとは関連していなかった。これは、シグナルがタンパク質の切断ではなく、発現量に影響を与える調節変異によって駆動されていることを示唆している。
• 相互作用： 多重比較補正を通過した遺伝子と性の相互作用、および遺伝子と剥奪度の相互作用は存在しなかった。しかし、rs261902およびrs76847656における性差の影響、およびrs73175505における剥奪度による差に関する、名目上のシグナルが観察された。
• DecodeMEによる検証： テスト可能な4つの再現された変異のうち、いずれもDecodeMEコホートにおいて関連を示さなかった。著者らは、この相違を、表現型の定義（DecodeMEはPEMを必須とする）、年齢分布、およびリクルート戦略（コミュニティベースかバイオバンクか）の違いに帰している。

意義および主張
本論文は、独立したコホートにおいて再現される遺伝的関連を報告した最初のME/CFS研究であり、7つの候補遺伝子座を特定したと主張している。著者らは、マルチエビデンス・フェノタイピングとTarGene統計フレームワークの使用により、従来の単一コードの研究が見逃していたシグナルを検出できたことを強調している。

本研究は、控えめな結論として、（コホート間の表現型の不均一性、集団の違い、および小規模なリプリケーション・コホートにおける「勝者の呪い」により）すべてのコホートで再現された変異は存在しないものの、再現された7つの遺伝子座は将来的な生物学的調査のための候補ターゲットを提供すると述べている。第13染色体のシグナルが被蓋におけるCLYBLの発現と共局在していることは、ミトコンドリアのビタミンB12代謝とME/CFSの病態生理学との潜在的な関連を浮き彫りにしている。著者らは、これらの知見が因果関係の証明にはならないことを明言しており、ME/CFSの生物学的メカニズムへの追跡調査のための基礎を提供するものであるとしている。また、DecodeMEとの重複の欠如は、疾患がどのように定義されるか（例：PEMの包含の有無）によって、遺伝的構造が異なる可能性があることを示唆していると警告している。