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

Probleemstelling
Myalgische Encefalomyelitis/Chronisch Vermoeidheidssyndroom (ME/CFS) is een invaliderende, vrouw-gebiaseerde ziekte die kampt met een gebrek aan diagnostische biomarkers, effectieve behandelingen of een goed begrepen etiologie. Eerdere genetische onderzoeken zijn gehinderd door kleine, heterogene cohorten en een afhankelijkheid van enkelvoudige diagnostische codes, wat leidde tot een gebrek aan reproduceerbare bevindingen. Voorheen waren er geen genomische risicoloci voor ME/CFS die repliceerden in onafhankelijke cohorten. Bovendien blijft de rol van sociaaleconomische status (SES) en gen-door-geslacht-interacties bij het modificeren van het risico controversieel vanwege tegenstrijdige rapporten en methodologische beperkingen in eerdere studies.

Methodologie
De auteurs voerden een genoombrede associatiestudie (GWAS) uit met behulp van gegevens uit de UK Biobank (UKB) en het All of Us (AoU) Research Program.

• Fenotyperingstrategie: In plaats van te vertrouwen op een enkele diagnostische code, hanteerde de studie een multi-evidence fenotyperingsstrategie om gevallen en controles te definiëren. Casussen moesten een zelfgerapporteerde klinische diagnose hebben, een slechte/matige algemene gezondheidsscore en specifieke bevestigende antwoorden op pijnvragenlijsten met betrekking tot ME/CFS. Controles moesten een goede/uitstekende gezondheidsscore hebben en geen bewijs van ME/CFS in medische dossiers of enquêtes.
• Cohorten: Drie disjuncte populaties werden vastgesteld:
  1. Discovery (UKB1): 1.268 gevallen en 113.132 controles.
  2. Replicatie 1 (UKB2): 319 gevallen en 28.471 controles.
  3. Replicatie 2 (AoU): 371 gevallen en 33.140 controles.
• Statistisch Kader: De studie maakte gebruik van TarGene (Targeted Genomic Estimation), een semi-parametrische, dubbel-robuuste schattingsworkflow gebaseerd op Targeted Learning. Deze benadering werd gekozen om de power te maximalen en bias door model-misspecificatie en populatiestratificatie te minimaliseren, waarmee de beperkingen van conventionele parametrische GWAS-methoden worden aangepakt.
• Analyse Pipeline:
  • Discovery: Een GWAS werd uitgevoerd op 670.243 genotypeerde SNP's in UKB1, waarbij de False Discovery Rate (FDR) werd gecontroleerd op 5%.
  • Replicatie: De 176 FDR-significante varianten uit de discovery-fase werden getest in de twee replicatiecohorten.
  • Fine-mapping en Colocalisatie: Voor gerepliceerde loci werd fine-mapping uitgevoerd met behulp van SuSiE om credible sets van causale varianten te genereren. Deze werden gecolocaliseerd met expression quantitative trait loci (eQTL's) uit de GTEx v10 database.
  • Interactieanalyse: Gen-door-geslacht en gen-door-deprivatie (met behulp van de Townsend Deprivation Index en Indices of Multiple Deprivation) interacties werden getest met behulp van TarGene.
  • Validatie: Er werd een lookup uitgevoerd in de DecodeME-cohort (2.143 gevallen), een studie die strikte klinische criteria gebruikt (Canadian Consensus en IoM 2015) die Post-Exertional Malaise (PEM) vereist.

Belangrijkste Resultaten

• Discovery: De analyse identificeerde 176 varianten die significant geassocieerd waren met het risico op ME/CFS (FDR < 5%) in de UKB1 discovery-cohort.
• Replicatie: Zeven genomische risicoloci repliceerden in ten minste één van de twee disjuncte replicatiecohorten (UKB2 of AoU). Echter, geen enkele variant bereikte significantie over alle drie de cohorten na correctie voor meervoudig testen.
  • De zeven gerepliceerde varianten zijn: rs115186419 (nabij ETV5), rs73175505 (nabij CSMD1), rs261902 (in BICD1), rs117553493 (nabij CLYBL), rs72741654 (nabij RORA), rs74963073 (in GRIN2A), en rs76847656 (nabij SBK1).
• Fine-mapping en Colocalisatie:
  • Fine-mapping op het chromosoom 13 locus loste een credible set op die colocaliseert met verminderde CLYBL-expressie in het putamen (basale ganglia). Het risico-allel is geassocieerd met een lagere CLYBL-expressie (Posterior Probability van H4 = 0,994).
  • Opvallend genoeg was de CLYBL Arg259 stop-gain variant (een bekende loss-of-function variant) niet geassocieerd met het risico op ME/CFS, wat suggereert dat het signaal wordt gedreven door regulatoire varianten die de expressiedosis beïnvloeden in plaats van eiwittruncatie.
• Interacties: Geen gen-door-geslacht of gen-door-deprivatie interacties overleefden de correctie voor meervoudig testen. Er werden echter nominale signalen waargenomen voor geslachts-differentiële effecten bij rs261902 en rs76847656, en een deprivatie-differentieel effect bij rs73175505.
• DecodeME Validatie: Geen van de vier testbare gerepliceerde varianten vertoonde associatie in de DecodeME-cohort. De auteurs schrijven deze discrepantie toe aan verschillen in fenotype-definitie (DecodeME vereist PEM; deze studie niet), leeftijdsverdelingen en wervingsstrategieën (community versus biobank).

Significantie en Claims
Het artikel beweert de eerste ME/CFS studie te zijn die genetische associaties rapporteert die repliceren in onafhankelijke cohorten, waarbij zeven kandidaat loci heeft geïdentificeerd. De auteurs benadrukken dat het gebruik van multi-evidence fenotypering en het TarGene statistische kader hen in staat stelde om signalen te detecteren die eerdere single-code studies hebben gemist.

De studie concludeert bescheiden dat hoewel geen enkele variant repliceerde over alle drie de cohorten (waarschijnlijk door fenotypische heterogeniteit, populatieverschillen en de "winner's curse" in kleinere replicatiecohorten), de zeven gerepliceerde loci kandidaat-targets bieden voor toekomstig biologisch onderzoek. De colocalisatie van het chromosoom 13 signaal met de CLYBL-expressie in het putamen benadrukt een potentiële link tussen mitochondriale vitamine B12-stofwisseling en de pathofysiologie van ME/CFS. De auteurs verklaren expliciet dat deze bevindingen geen causaal bewijs vormen, maar een fundament bieden voor vervolgstudies naar de biologische mechanismen van ME/CFS. Ze waarschuwen dat het gebrek aan overlap met DecodeME suggereert dat de genetische architectuur kan verschillen afhankelijk van hoe de ziekte wordt gedefinieerd (bijv. inclusie van PEM).