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

Planteamiento del Problema
La Encefalomielitis Miálgica/Síndrome de Fatiga Crónica (EM/SFC) es una enfermedad debilitante, con sesgo femenino, que carece de biomarcadores diagnósticos, tratamientos eficaces o una etiología bien comprendida. Las investigaciones genéticas previas se han visto obstaculizadas por cohortes pequeñas y heterogéneas, así como por la dependencia de códigos diagnósticos de una sola línea, lo que ha derivado en una falta de hallazgos reproducibles. Previamente, ningún locus de riesgo genómico para la EM/SFC se había replicado en cohortes independientes. Además, el papel del estatus socioeconómico (SES) y las interacciones gen-por-sexo en la modificación del riesgo sigue siendo contencioso debido a informes contradictorios y limitaciones metodológicas en estudios previos.

Metodología
Los autores llevaron a cabo un estudio de asociación de genoma completo (GWAS) utilizando datos del UK Biobank (UKB) y del Programa de Investigación All of Us (AoU).

• Estrategia de Fenotipado: En lugar de depender de un único código diagnóstico, el estudio empleó una estrategia de fenotipado de evidencia múltiple para definir casos y controles. Para los casos, se requería tener diagnósticos clínicos autoinformados, calificaciones de salud general pobre/regular, y respuestas afirmativas específicas en cuestionarios de dolor relacionadas con la EM/SFC. Los controles debían tener calificaciones de salud buena/excelente y no presentar evidencia de EM/SFC en sus registros médicos o encuestas.
• Cohortes: Se establecieron tres poblaciones disjuntas:
  1. Descubrimiento (UKB1): 1,268 casos y 113,132 controles.
  2. Replicación 1 (UKB2): 319 casos y 28,471 controles.
  3. Replicación 2 (AoU): 371 casos y 33,140 controles.
• Marco Estadístico: El estudio utilizó TarGene (Estimación Genómica Dirigida), un flujo de trabajo de estimación semiparamétrico y doblemente robusto basado en el Aprendizaje Dirigido (Targeted Learning). Se eligió este enfoque para maximizar el poder y minimizar el sesgo por especificación errónea del modelo y estratificación de la población, abordando las limitaciones de los métodos convencionales de GWAS paramétricos.
• Pipeline de Análisis:
  • Descubrimiento: Se realizó un GWAS sobre 670,243 SNPs genotipados en UKB1, controlando la Tasa de Falso Descubrimiento (FDR) al 5%.
  • Replicación: Se probaron las 176 variantes significativas de FDR en las dos cohortes de replicación.
  • Mapeo Fino y Colocalización: Para los loci replicados, se realizó un mapeo fino utilizando SuSiE para generar conjuntos creíbles de variantes causales. Estos se colocalizaron con los QTL de expresión (eQTL) de la base de datos GTEx v10.
  • Análisis de Interacción: Las interacciones gen-por-sexo y gen-por-privación (utilizando el Índice de Privación de Townsend e Índices de Privación Múltiple) se probaron mediante TarGene.
  • Validación: Se realizó una búsqueda en la cohorte DecodeME (2,143 casos), un estudio que utiliza criterios clínicos estrictos (Consenso Canadiense e IoM 2015) que requiere Malestar Pos-Esfuerzo (PEM).

Resultos Clave

• Descubrimiento: El análisis identificó 176 variantes significativamente asociadas con el riesgo de EM/SFC (FDR < 5%) en la cohorte de descubrimiento UKB1.
• Replicación: Siete loci de riesgo genómico se replicaron en al menos una de las dos cohortes de replicación disjuntas (UKB2 o AoU). Sin embargo, ninguna variante única alcanzó la significancia en las tres cohortes tras la corrección por pruebas múltiples.
  • Las siete variantes replicadas son: rs115186419 (cerca de ETV5), rs73175505 (cerca de CSMD1), rs261902 (en BICD1), rs117553493 (cerca de CLYBL), rs72741654 (cerca de RORA), rs74963073 (en GRIN2A) y rs76847656 (cerca de SBK1).
• Mapeo Fino y Colocalización:
  • En el locus del cromosoma 13, el mapeo fino resolvió un conjunto creíble que se colocaliza con una expresión reducida de CLYBL en el putamen (ganglios basales). El alelo de riesgo está asociado con una menor expresión de CLYBL (Probabilidad Posterior de H4 = 0.994).
  • Notablemente, la variante de pérdida de función (stop-gain) Arg259 de CLYBL no estuvo asociada con el riesgo de EM/SFC, lo que sugiere que la señal es impulsada por variantes regulatorias que afectan la dosis de expresión en lugar de la truncación de la proteína.
• Interacciones: Ninguna interacción gen-por-sexo o gen-por-privación sobrevivió a la corrección por pruebas múltiples. Sin embargo, se observaron señales nominales para efectos diferenciales de sexo en rs261902 y rs76847656, y un efecto diferencial de privación en rs73175505.
• Validación con DecodeME: Ninguna de las cuatro variantes replicables probables mostró asociación en la cohorte DecodeME. Los autores atribuyen esta discrepancia a diferencias en la definición del fenotipo (DecodeME requiere PEM; este estudio no), distribuciones de edad y estrategias de reclutamiento (comunidad frente a biobanco).

Significancia y Reivindicaciones
El artículo afirma ser el primer estudio de EM/SFC que reporta asociaciones genéticas que se replican en cohortes independientes, identificando siete loci candidatos. Los autores enfatizan que el uso de fenotipado de evidencia múltiple y el marco estadístico TarGene permitieron la detección de señales que los estudios previos de código único pasaron por alto.

El estudio concluye modestamente que, si bien ninguna variante se replicó en las tres cohortes (probablemente debido a la heterogeneidad fenotípica, diferencias poblacionales y la "maldición del ganador" en cohortes de replicación más pequeñas), los siete loci replicados proporcionan objetivos candidatos para la investigación biológica futura. La colocalización de la señal del cromosoma 13 con la expresión de CLYBL en el putamen resalta un posible vínculo entre el metabolismo de la vitamina B12 mitocondrial y la patofisiología de la EM/SFC. Los autores declaran explícitamente que estos hallazgos no constituyen una prueba causal, sino que ofrecen una base para estudios de seguimiento sobre los mecanismos biológicos de la EM/SFC. Advierten que la falta de superposición con DecodeME sugiere que la arquitectura genética puede diferir según cómo se defina la enfermedad (por ejemplo, la inclusión de PEM).