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

Énoncé du problème
L'encéphalomyélite myalgique/syndrome de fatigue chronique (EM/SFC) est une maladie débilitante, à prédominance féminine, dépourvue de biomarqueurs diagnostiques, de traitements efficaces ou d'une étiologie bien comprise. Les recherches génétiques antérieures ont été entravées par des cohortes de petite taille et hétérogènes, ainsi que par le recours à des codes de diagnostic uniques, ce qui a conduit à un manque de résultats reproductibles. Aucun locus de risque génomique pour l'EM/SFC n'avait été précédemment répliqué dans des cohortes indépendantes. De plus, le rôle du statut socio-économique (SSE) et des interactions gène-par-sexe dans la modification du risque reste controversé en raison de rapports contradictoires et des limites méthodologiques des études antérieures.

Méthodologie
Les auteurs ont mené une étude d'association pangénomique (GWAS) en utilisant les données de la UK Biobank (UKB) et du programme de recherche All of Us (AoU).

• Stratégie de phénotypage : Au lieu de s'appuyer sur un code de diagnostic unique, l'étude a employé une stratégie de phénotypage multi-preuves pour définir les cas et les témoins. Les cas devaient présenter des diagnostics cliniques auto-déclarés, des évaluations de l'état de santé général « médiocre » ou « passable », et des réponses affirmatives spécifiques aux questionnaires sur la douleur concernant l'EM/SFC. Les témoins devaient présenter des évaluations de santé « bonnes » ou « excellentes » et ne présenter aucune trace d'EM/SFC dans leurs dossiers médicaux ou sondages.
• Cohortes : Trois populations disjointes ont été établies :
  1. Découverte (UKB1) : 1 268 cas et 113 132 témoins.
  2. Réplication 1 (UKB2) : 319 cas et 28 471 témoins.
  3. Réplication 2 (AoU) : 371 cas et 33 140 témoins.
• Cadre statistique : L'étude a utilisé TarGene (Targeted Genomic Estimation), un flux de travail d'estimation semi-paramétrique et doublement robuste basé sur l'apprentissage ciblé (Targeted Learning). Cette approche a été choisie pour maximiser la puissance tout en minimisant les biais de spécification incorrecte du modèle et de stratification de la population, répondant ainsi aux limites des méthodes conventionnelles de GWAS paramétriques.
• Pipeline d'analyse :
  • Découverte : Une GWAS a été réalisée sur 670 243 SNP génotypés dans UKB1, en contrôlant le taux de fausses découvertes (FDR) à 5 %.
  • Réplication : Les 176 variants significatifs au niveau du FDR identifiés lors de la phase de découverte ont été testés dans les deux cohortes de réplication.
  • Fine-mapping et colocalisation : Pour les locus répliqués, le fine-mapping a été réalisé à l'aide de SuSiE pour générer des ensembles crédibles de variants causaux. Ces derniers ont été colocalisés avec des QTL d'expression (eQTL) issus de la base de données GTEx v10.
  • Analyse d'interaction : Les interactions gène-par-sexe et gène-par-privation (utilisant l'indice de privation de Townsend et les indices de privation multiples) ont été testées via TarGene.
  • Validation : Une recherche (look-up) a été effectuée dans la cohorte DecodeME (2 143 cas), une étude utilisant des critères cliniques stricts (Consensus canadien et IoM 2015) qui exige la présence de Malaise Post-Effort (MPE).

Résultats clés

• Découverte : L'analyse a identifié 176 variants significativement associés au risque d'EM/SFC (FDR < 5 %) dans la cohorte de découverte UKB1.
• Réplication : Sept locus de risque génomique se sont répliqués dans au moins l'une des deux cohortes de réplication disjointes (UKB2 ou AoU). Cependant, aucun variant unique n'a atteint la significativité dans les trois cohortes après correction pour tests multiples.
  • Les sept variants répliqués sont : rs115186419 (proche de ETV5), rs73175505 (proche de CSMD1), rs261902 (dans BICD1), rs117553493 (proche de CLYBL), rs72741654 (proche de RORA), rs74963073 (dans GRIN2A), et rs76847656 (proche de SBK1).
• Fine-mapping et colocalisation :
  • Au niveau du locus du chromosome 13, le fine-mapping a résolu un ensemble crédible se colocalisant avec une réduction de l'expression de CLYBL dans le putamen (corps strié/basal). L'allèle de risque est associé à une expression plus faible de CLYBL (Probabilité postérieure de H4 = 0,994).
  • Notamment, le variant de perte de fonction (stop-gain) Arg259 de CLYBL n'était pas associé au risque d'EM/SFC, suggérant que le signal est piloté par des variants régulateurs affectant la dose d'expression plutôt que par une troncature protéique.
• Interactions : Aucune interaction gène-par-sexe ou gène-par-privation n'a survécu à la correction pour tests multiples. Cependant, des signaux nominaux ont été observés pour les effets différentiels selon le sexe à rs261902 et rs76847656, et un effet différentiel lié à la privation à rs73175505.
• Validation DecodeME : Aucun des quatre variants testables répliqués n'a montré d'association dans la cohorte DecodeME. Les auteurs attribuent cette divergence aux différences de définition du phénotype (DecodeME exige le MPE ; cette étude non), aux distributions d'âge et aux stratégies de recrutement (communautaire vs biobanque).

Signification et affirmations
L'article affirme être la première étude sur l'EM/SFC à rapporter des associations génétiques qui se répliquent dans des cohortes indépendantes, identifiant sept locus candidats. Les auteurs soulignent que l'utilisation d'un phénotypage multi-preuves et du cadre statistique TarGene a permis de détecter des signaux qui avaient échappé aux études précédentes basées sur un code de diagnostic unique.

L'étude conclut modestement que, bien qu'aucun variant ne se soit répliqué dans les trois cohortes (probablement en raison de l'hétérogénéité phénotypique, des différences de population et du « biais du vainqueur » dans les cohortes plus petites), les sept locus répliqués constituent des cibles candidates pour de futures investigations biologiques. La colocalisation du signal du chromosome 13 avec l'expression de CLYBL dans le putamen met en évidence un lien potentiel entre le métabolisme de la vitamine B12 mitochondriale et la pathophysiologie de l'EM/SFC. Les auteurs précisent explicitement que ces résultats ne constituent pas une preuve de causalité mais offrent une base pour des études de suivi sur les mécanismes biologiques de l'EM/SFC. Ils avertissent que l'absence de chevauchement avec DecodeME suggère que l'architecture génétique peut varier selon la manière dont la maladie est définie (par exemple, l'inclusion du MPE).