Common cis-eQTLs do not explain chromosome 21 expression deviations from the ploidy expectation in Down syndrome

Resumen Técnico: Los cis-eQTL comunes no explican las desviaciones en la expresión génica del cromosoma 21 en la trisomía 21

Planteamiento del problema
El síndrome de Down (SD) es causado por la trisomía 21 (T21), la triplicación del cromosoma humano 21 (HSA21). Aunque la mayoría de los genes de HSA21 exhiben un aumento esperado de 1.5 veces en su expresión en individuos con SD en comparación con individuos euploides (disómicos, D21), un pequeño subconjunto de genes se desvía de esta expectativa, mostrando niveles de expresión más cercanos a las líneas base de los euploides. Investigaciones previas, incluyendo un estudio de base familiar por Hunter et al., sugirieron que la variación regulatoria genética común (específicamente los locus de rasgos cuantitativos de expresión en cis o cis-eQTL) podría explicar estas desviaciones en individuos específicos. Sin embargo, no estaba claro si la variación regulatoria cis común explica estas desviaciones a nivel poblacional. Este estudio investiga si las diferencias en la frecuencia o el efecto de los cis-eQTL comunes entre las poblaciones T21 y D21 explican el amortiguamiento (buffering) o la amplificación de los efectos de dosis de genes específicos.

Metodología
Los autores utilizaron el conjunto de datos transcriptómicos específicos de SD más grande hasta la fecha, derivado del Human Trisome Project (HTP), que comprende datos de RNA-seq de sangre total de 274 participantes con T21 y 80 participantes euploides (D21). Para asegurar la precisión en un contexto triploide, el estudio empleó un riguroso flujo de trabajo analítico consciente de la ploidía:

1. Genotipado consciente de la ploidía: Los flujos de trabajo estándar de genotipado germinal asumen la diploidía y asignan erróneamente la dosis alélica en HSA21. Los autores volvieron a llamar las variantes directamente de los datos de secuenciación de genoma completo (WGS) utilizando Octopus, estableciendo explícitamente la ploidía de HSA21 como 3 para las muestras T21 y 2 para las D21. Esto generó recuentos de alelos precisos (0–3 copias) para la llamada de variantes.
2. Expresión diferencial consciente de la ploidía: Utilizando DESeq2, los autores realizaron un análisis de expresión diferencial con un modelo que estableció el cambio de pliegue (fold-change) esperado para los genes de HSA21 en las muestras T21 en 1.5 (log2 1.5 ≈ 0.585), en lugar del valor predeterminado de 1.0. Las covariables incluyeron edad, sexo, IMC, fuente de la muestra y 19 proporciones de tipos celulares derivadas de citometría de masa para controlar las diferencias composicionales.
3. Identificación de desviaciones: Los genes se clasificaron como "desviados" si su cambio de pliegue corregido por ploidía difería del 1.5 en al menos un 33% (p ajustado < 0.01).
4. Pruebas de asociación de cis-eQTL: Para cada gen desviado, los autores cruzaron los cis-eQTL comunes (ventana de ±1 Mb) con el conjunto de datos de sangre total GTEx v10. Regredieron la expresión ajustada por covariables contra la dosis de genotipo (0–3 copias) dentro de la cohorte T21. La significancia se evaluó mediante pruebas de permutación a nivel de gen (1,000 permutaciones) ajustadas para pruebas múltiples.
5. Análisis comparativo: El estudio comparó la tasa de detección y los tamaños de efecto de los cis-eQTL en genes desviados frente a los genes de HSA21 que muestran el cambio de 1.5 veces esperado. También modelaron si las diferencias en la frecuencia de alelos entre los participantes de T21 y los donantes de GTEx podían explicar matemáticamente las desviaciones de expresión observadas.

Resultados clave

• Desviaciones de expresión: Tras filtrar por bajos recuentos de lecturas y alto contenido de repeticiones, y aplicar la corrección de ploidía, 13 de los 318 genes objetivo de HSA21 se desviaron significativamente del cambio esperado de 1.5 veces (8.5%). Entre los 10 genes desviados con una variante de GTEx para probar, 5 fueron superiores a lo esperado y 5 fueron inferiores a lo esperado.
• Detección de cis-eQTL: Entre los 10 genes desviados que eran testeables (GTEx tenía datos y variantes), 8 presentaron un cis-eQTL detectado en la cohorte T21. Sin embargo, 59 de los 99 genes de HSA21 con el cambio de 1.5 veces esperado también tenían cis-eQTL significativos.
• Tamaños de efecto y frecuencias: La magnitud de los efectos de los cis-eQTL en los participantes de T21 fue similar a la observada en GTEx (mediana del cambio de pliegue alélico ~0.73 para genes desviados frente a ~0.76 para genes de dosis esperada). Crucialmente, las frecuencias de alelos de las variantes probadas en los participantes de T21 no difirieron sistemáticamente de las de los donantes de GTEx (diferencia mediana −0.0005).
• Poder explicativo: El modelado matemático demostró que incluso si los participantes de T21 portaran el alelo alternativo en las tres copias o en ninguna, las diferencias de frecuencia observadas no podrían producir la magnitud de las desviaciones de expresión observadas. Las mejores variantes explicaban solo entre el 2% y el 33% de la varianza de la expresión entre los participantes de T21.
• Casos específicos: Dos genes con expresión menor a la esperada (BACE2 y OLIG2) no presentaron cis-eQTL detectados en la cohorte T21, a pesar de tener fuertes eQTL en GTEx. Sin embargo, la falta de detección no se debió a la baja potencia, y sus patrones de desviación se mantuvieron robustos en los modelos de covariables.

Significancia y afirmaciones
El artículo concluye que, si bien la variación regulatoria cis común actúa sobre los genes de HSA21 en individuos con T21 de manera similar a como actúa en individuos euploides, esta no explica qué genes específicos se desvían de la dosis de ploidía esperada a nivel poblacional.

• Población vs. Individuo: Los autores aclaran que, aunque el genotipo de un individuo puede explicar sus desviaciones específicas (como mostró el estudio de base familiar de Hunter et al.), los cis-eQTL comunes no crean una diferencia sistemática entre las poblaciones T21 y D21. Las desviaciones no son impulsadas por un cambio en la frecuencia de los alelos reguladores en la población de nacidos vivos con T21.
• Rechazo de la hipótesis genética: Los hallazgos no respaldan la hipótesis de que la variación regulatoria cis común es el mecanismo principal para el amortiguamiento o la amplificación de la dosis de genes específicos de HSA21 en el SD.
• Mecanismos alternativos: Los autores sugieren que los mecanismos que impulsan estas desviaciones probablemente residen en otros lugares, involucrando potencialmente factores no genéticos como bucles de retroalimentación autorreguladora, silenciamiento post-transcripcional por miRNAs triplicados, splicing alternativo, remodelación de la cromatina o variantes regulatorias raras no capturadas por el mapeo de eQTL comunes.

El estudio enfatiza que la mayoría de los genes de HSA21 siguen el aumento de dosis de 1.5 veces esperado cuando se analizan con métodos apropiados conscientes de la ploidía, y que las desviaciones observadas en un pequeño subconjunto de genes siguen siendo una pregunta abierta que requiere investigación más allá de la variación regulatoria cis común.