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

Problema
A Encefalomielite Miálgica/Síndrome de Fadiga Crônica (EM/FC) é uma doença debilitante, com viés feminino, que carece de biomarcadores diagnósticos, tratamentos eficazes ou uma etiologia bem compreendida. Pesquisas genéticas anteriores foram dificultadas por coortes pequenas e heterogêneas e pela dependência de códigos diagnósticos de linha única, levando à falta de achados reprodutíveis. Nenhum locus de risco genômico para EM/FC havia sido replicado anteriormente em coortes independentes. Além disso, o papel do status socioeconômico (SSE) e das interações gene-por-sexo na modificação do risco permanece contencioso devido a relatórios conflitantes e limitações metodológicas em estudos anteriores.

Metodologia
Os autores conduziram um estudo de associação de genoma completo (GWAS) utilizando dados do UK Biobank (UKB) e do Programa de Pesquisa All of Us (AoU).

• Estratégia de Fenotipagem: Em vez de depender de um único código diagnóstico, o estudo empregou uma estratégia de fenotipagem de evidência múltipla para definir casos e controles. Para os casos, exigiu-se diagnósticos clínicos autorrelatados, classificações de saúde geral pobre/regular e respostas afirmativas específicas em questionários de dor em relação à EM/FC. Os controles exigiam classificações de saúde boa/excelente e nenhuma evidência de EM/FC em registros médicos ou pesquisas.
• Coortes: Três populações disjuntas foram estabelecidas:
  1. Descoberta (UKB1): 1.268 casos e 113.132 controles.
  2. Replicação 1 (UKB2): 319 casos e 28.471 controles.
  3. Replicação 2 (AoU): 371 casos e 33.140 controles.
• Estrutura Estatística: O estudo utilizou o TarGene (Estimativa Genômica Direcionada), um fluxo de trabalho de estimativa semiparamétrica e duplamente robusto baseado em Aprendizado Direcionado (Targeted Learning). Essa abordagem foi escolhida para maximizar o poder enquanto minimiza o viés de erro de especificação de modelo e estratificação populacional, abordando as limitações dos métodos convencionais de GWAS paramétricos.
• Pipeline de Análise:
  • Descoberta: Um GWAS foi realizado em 670.243 SNPs genotipados no UKB1, controlando a Taxa de Descoberta Falsa (FDR) em 5%.
  • Replicação: As 176 variantes significativamente associadas na fase de descoberta foram testadas nas duas coortes de replicação.
  • Mapeamento Fino (Fine-mapping) e Colocalização: Para os loci replicados, o mapeamento fino foi realizado usando SuSiE para gerar conjuntos de crença de variantes causais. Estes foram colocalizados com QTLs de expressão (eQTLs) do banco de dados GTEx v10.
  • Análise de Interação: Interações gene-por-sexo e gene-por-privação (usando o Índice de Privação de Townsend e Índices de Múltipla Privação) foram testadas usando TarGene.
  • Validação: Uma busca foi realizada na coorte DecodeME (2.143 casos), um estudo que utiliza critérios clínicos rigorosos (Consenso Canadense e IoM 2015) que exige Mal-estar Pós-Esforço (PEM).

Principais Resultados

• Descoberta: A análise identificou 176 variantes significativamente associadas ao risco de EM/FC (FDR < 5%) na coorte de descoberta UKB1.
• Replicação: Sete loci de risco genômico replicaram-se em pelo menos uma das duas coortes de replicação disjuntas (UKB2 ou AoU). No entanto, nenhuma variante única alcançou significância em todas as três coortes após a correção para múltiplos testes.
  • As sete variantes replicadas são: rs115186419 (próximo a ETV5), rs73175505 (próximo a CSMD1), rs261902 (em BICD1), rs117553493 (próximo a CLYBL), rs72741654 (próximo a RORA), rs74963073 (em GRIN2A) e rs76847656 (próximo a SBK1).
• Mapeamento Fino e Colocalização:
  • No locus do cromossomo 13, o mapeamento fino resolveu um conjunto de crença colocalizado com a redução da expressão de CLYBL no putâmen (corpo estriado/núcleos da base). O alelo de risco está associado à menor expressão de CLYBL (Probabilidade Posterior de H4 = 0,994).
  • Notavelmente, a variante de parada de ganho Arg259 de CLYBL (uma conhecida variante de perda de função) não foi associada ao risco de EM/FC, sugerindo que o sinal é impulsionado por variantes regulatórias que afetam a dosagem da expressão, em vez de truncamento proteico.
• Interações: Nenhuma interação gene-por-sexo ou gene-por-privação sobreviveu à correção para múltiplos testes. No entanto, sinais nominais foram observados para efeitos diferenciais de sexo em rs261902 e rs76847656, e um efeito diferencial de privação em rs73175505.
• Validação DecodeME: Nenhuma das quatro variantes testáveis replicadas mostrou associação na coorte DecodeME. Os autores atribuem essa discrepância às diferenças na definição do fenótipo (DecodeME exige PEM; este estudo não), distribuições de idade e estratégias de recrutamento (comunidade vs. biobanco).

Significância e Alegações
O artigo afirma ser o primeiro estudo de EM/FC a relatar associações genéticas que se replicam em coortes independentes, identificando sete loci candidatos. Os autores enfatizam que o uso de fenotipagem de evidência múltipla e a estrutura estatística TarGene permitiram a detecção de sinais que estudos anteriores de código único perderam.

O estudo conclui modestamente que, embora nenhum variante tenha se replicado em todas as três coortes (provavelmente devido à heterogeneidade fenotípica, diferenças populacionais e o "viés do vencedor" em coortes menores), os sete loci replicados fornecem alvos candidatos para investigações biológicas futuras. A colocalização do sinal do cromossomo 13 com a expressão de CLYBL no putâmen destaca uma potencial ligação entre o metabolismo da vitamina B12 mitocondrial e a fisiopatologia da EM/FC. Os autores declaram explicitamente que esses achados não constituem prova causal, mas oferecem uma base para estudos de acompanhamento sobre os mecanismos biológicos da EM/FC. Eles alertam que a falta de sobreposição com o DecodeME sugere que a arquitetura genética pode diferir dependendo de como a doença é definida (ex: inclusão de PEM).