eQTL — No interpretar la asociación de expresión como causalidad
Por qué es importante saberlo
Los GWAS de enfermedades identifican loci asociados, pero no revelan directamente a qué gen o proceso regulatorio afecta esa señal. El análisis de eQTL (locus cuantitativo del rasgo de expresión) mide si las diferencias de genotipo y las diferencias de expresión génica entre individuos varían conjuntamente, lo que ayuda a restringir las posibles vías regulatorias.
Sin embargo, el hecho de que haya asociación no significa que la variante cause directamente un cambio en la expresión. El desequilibrio de ligamiento, la composición celular y los factores ambientales pueden generar patrones similares. Los eQTL constituyen una capa de evidencia para definir objetivos de validación posterior, no una conclusión causal.
Conceptos clave
Por lo general, se denomina cis-eQTL a las asociaciones con genes cercanos y trans-eQTL a las asociadas con loci distantes o en otros cromosomas; sin embargo, los criterios de distancia y las definiciones analíticas deben especificarse según cada estudio. Si se miden uniones por empalme (splice junctions) o isoformas en lugar de la expresión total, se pueden identificar sQTL. La elección del rasgo molecular como unidad determina parte de la interpretación de los resultados.
El análisis introduce genotipos, expresión normalizada y covariables de ascendencia y técnica en un modelo de regresión, y prueba numerosas parejas variante-gen. Aunque se corrige el contraste múltiple y se ajustan los factores ocultos, una corrección excesiva puede eliminar señales biológicas reales. Con tamaños muestrales pequeños y frecuencias alélicas bajas, resulta difícil detectar efectos pequeños.
Contexto tisular y celular
La regulación génica varía según el tejido, el tipo celular, la etapa de desarrollo, las estimulaciones y el estado de enfermedad. Un eQTL evidente en hígado podría no observarse en sangre, y existen eQTL de respuesta que solo aparecen tras la estimulación inmune. Recursos multi-tisulares como GTEx proporcionan una base para comparar esta dependencia contextual.
La expresión en tejido bulk es un promedio de múltiples tipos celulares. Si el genotipo y la proporción celular varían conjuntamente, los efectos de la composición celular pueden parecerse a la regulación génica, y las señales de tipos celulares raros pueden quedar enmascaradas por el promedio. El enfoque de célula única aumenta la resolución, pero introduce otras limitaciones como la dispersión (sparsity), los lotes experimentales y los tamaños muestrales pequeños por célula.
Cuando se superpone con GWAS
Que la variante líder del GWAS y la variante líder del eQTL sean iguales o estén cerca no permite confirmar que una única variante causal cause ambos rasgos. Diferentes variantes dentro de un bloque de desequilibrio de ligamiento pueden afectar independientemente la enfermedad y la expresión. Primero, se deben evaluar las señales independientes y el conjunto creíble mediante análisis condicional y fine-mapping.
La colocalización compara probabilísticamente si dos patrones de asociación son compatibles con una señal causal compartida. El resultado depende de la cobertura de variantes, la referencia de desequilibrio de ligamiento, los supuestos sobre el número de variantes causales y las probabilidades previas. Incluso una probabilidad posterior alta no demuestra de forma definitiva el gen diana ni el mecanismo de la enfermedad.
Pequeño ejemplo
Supongamos que un locus de GWAS de una enfermedad autoinmunitaria se asocia con una menor expresión del gen A en células inmunitarias. Si ambas señales colocalizan en el tipo celular pertinente y la expresión específica de alelo y la anotación de potenciadores apoyan la misma dirección, el gen A se convierte en un buen candidato. Aun así, para reforzar una afirmación causal deben reproducirse la función inmunitaria y el fenotipo relacionado con la enfermedad mediante perturbación CRISPR o en una cohorte independiente.
Puntos que pueden generar confusión
- No se garantiza que la variante principal sea la variante causal.
- Tampoco se garantiza que el gen más cercano sea el gen diana.
- La ausencia de un eQTL no significa que no exista un efecto regulador, sino que no fue detectado en el tejido, estado y muestra medidos.
- La colocalización no sustituye los experimentos de contacto físico ni de rescate funcional.
Limitaciones
Si la ascendencia y la representatividad tisular de un recurso de eQTL son limitadas, resulta difícil generalizar a otras poblaciones y estados. El método de medición de la expresión, la anotación de referencia y la versión de la canalización también pueden cambiar los resultados. Al utilizar resultados públicos deben registrarse el alelo de efecto, la dirección, el tejido, el tamaño de la muestra y la versión del análisis.
Lectura complementaria
El paso desde un conjunto creíble de GWAS hasta un mecanismo regulador se relaciona con el concepto GWAS no codificante y contacto de cromatina; la accesibilidad de la cromatina se conecta con el concepto ATAC-seq.