Resumen Técnico: Determinantes Genéticos Comunes de la Respuesta Inmune Dinámica in vivo Humana a Mycobacterium tuberculosis

Planteamiento del Problema

Aunque la mayoría de las infecciones por Mycobacterium tuberculosis (Mtb) no progresan a la enfermedad de tuberculosis (TB) activa, los mecanismos que determinan estos resultados clínicos divergentes aún se comprenden poco. Si bien se sabe que las respuestas inmunitarias que involucran células T CD4, macrófagos, IFN-γ y TNF son necesarias para una inmunidad protectora, las deficiencias en estas vías específicas no explican la mayoría de los casos de enfermedad incidente en la población general. Los estudios de asociación de genoma completo (GWAS) no han logrado identificar loci de susceptibilidad a la TB estadísticamente robustos y reproducibles. Además, los análisis existentes de loci de rasgos de expresión cuantitativa (eQTL) se han limitado casi exclusivamente a la estimulación ex vivo de poblaciones de células inmunitarias seleccionadas, careciendo de validación dentro de respuestas inmunitarias in vivo complejas, dinámicas y multicelulares en un contexto relevante para la enfermedad.

Metodología

Los autores utilizaron la prueba de la tuberculina (TST) como un modelo de desafío experimental in vivo estandarizado para investigar los determinantes genéticos de la respuesta inmunitaria a Mtb.

Cohorte y Muestreo:

• Participantes: 267 individuos (134 hombres, 140 mujeres) de diversas etnias, que comprenden dos cohortes: aquellos con infección latente por TB (ILTB) y pacientes con TB activa (ATB) dentro de los 28 días posteriores al diagnóstico/inicio del tratamiento.
• Biopsias: Se obtuvieron perfiles de secuenciación de ARN (RNA-seq) pareados de las biopsias de la TST en el Día 2 (D2) y el Día 7 (D7) para los participantes con ILTB, y biopsias de D2 para los participantes con ATB. Las biopsias de inyección de solución salina de control se utilizaron como línea base.
• Genotipado: Se realizó el genotipado de genoma completo en 267 individuos utilizando el array Illumina Infinium Global Diversity Array-8, seguido de imputación al panel Haplotype Reference Consortium.

Pipeline Analítico:

1. Transcriptómica: Los datos de RNA-seq fueron procesados (recortados, alineados a hg19/GRCh38 y cuantificados). La expresión diferencial se analizó para definir el "transcriptoma de la respuesta a la TST" (genes significativamente sobreexpresados en comparación con los controles de solución salina).
2. Mapeo de eQTL:
   • Cis-eQTL: Se utilizaron modelos lineales mixtos (MatrixeQTL/lmerTest) para identificar asociaciones entre SNPs y la expresión génica dentro de ±1 Mb del sitio de inicio de la transcripción (TSS). Los análisis se realizaron en los conjuntos de datos combinados de D2/D7 de ILTB y D2 de ATB (415 muestras) para maximizar la potencia.
   • Trans-eQTL: Se probaron las asociaciones para SNPs ubicados a >5 Mb del TSS del gen, centrándose en los cis-eQTL líderes identificados en D2 o D7.
3. Análisis Funcional y de Colocalización:
   • Aprendizaje Profundo (Deep Learning): Se utilizaron modelos in silico para predecir el impacto funcional de los SNPs líderes en las secuencias reguladoras (promotores/potenciadores).
   • Colocalización: Se realizó un análisis de colocalización bayesiana (coloc) para probar si los eQTL de la TST comparten variantes causales con las señales de GWAS para TB y la enfermedad de Crohn.
   • Heredabilidad: Se utilizó la regresión de puntuación de desequilibrio de ligamiento estratificada (S-LDSC) para estimar la proporción de heredabilidad explicada por los eQTL de la respuesta a la TST en diversas enfermedades inmunológicas.
4. Análisis del Repertorio de Células T: Los datos de secuenciación de TCR de las biopsias de D7 se analizaron para identificar metaclones de células T reactivas a Mtb (específicamente restringidos por HLA-DRB1) y correlacionar su expansión con variantes genéticas.

Resultados Clave

1. Variación Interindividual y Descubrimiento de Cis-eQTL:

• El estudio identificó una sustancial variación interindividual en la expresión génica dentro del transcriptoma de la TST, revelando una regulación independiente de IFN-γ.
• 7.411 genes mostraron asociaciones significativas de SNPs de acción cis (eGenes). De estos, 1.719 fueron genes de respuesta a la TST.
• La señal de cis-eQTL más fuerte se mapeó a ERAP2 (Aminopeptidasa 2 del Retículo Endoplásmico).
• El análisis de colocalización vinculó la reducción de la expresión de ERAP2 con un aumento del riesgo de TB en los datos de GWAS, mientras que el aumento de la expresión se asoció con el riesgo de la enfermedad de Crohn, resaltando un compromiso (trade-off) entre la susceptibilidad a enfermedades autoinmunes e infecciosas.

2. Contexto Celular y Vías:

• La colocalización con catálogos de eQTL públicos y datos de célula única indicó que los cis-eQTL identificados operan principalmente en los linajes de células mieloides (monocitos/macrófagos) y células T.
• Análisis de Heredabilidad: Las vías relacionadas con la presentación de antígenos de la clase II de HLA y la activación de células T mostraron la mayor heredabilidad. Específicamente, los niveles de expresión de HLA-DRB1 en las TST de D2 están determinados genéticamente y se correlacionan con la expansión subsiguiente de los metaclones de células T reactivas a Mtb restringidos por HLA-DRB1 en la TST de D7.
• El estudio identificó tanto mecanismos específicos de antígeno (vía eQTL de expresión de HLA-DRB1) como mecanismos independientes de antígeno (vía regulación del ciclo celular) para la variación de la respuesta de las células T.

3. Trans-eQTL y Regulación del Ciclo Celular:

• El análisis de trans-eQTL reveló un programa proliferativo centrado en NCAPD3 (una subunidad del complejo Condensina II).
• Un único cis-eQTL en NCAPD3 impulsó la mayoría de las trans-asociaciones (81,6%), las cuales estuvieron enriquecidas para genes del ciclo celular y la mitosis.
• El alelo de NCAPD3 se asoció significativamente con la frecuencia de los clones de células T expandidos en las TST de D7, lo que sugiere un determinante de la inmunidad de las células T regulado genéticamente e independiente del antígeno. Esta asociación fue validada en un estudio independiente de eQTL de célula única de células T estimuladas.

4. Integración con GWAS:

• Los eQTL de la respuesta a la TST mostraron un enriquecimiento significativo para la heredabilidad de diversas enfermedades inmunológicas en los datos de GWAS.
• En la enfermedad de Crohn, 17 SNPs mostraron una probabilidad posterior >75% de colocalización con los eQTL de la TST.
• En la TB, aunque ningún locus individual alcanzó la significancia a nivel de genoma completo en el metaanálisis, el estudio identificó polimorfismos específicos dentro del locus ERAP2 (por ejemplo, rs3733904) asociados con la reducción de la expresión y el aumento del riesgo de TB.

Significancia y Reivindicaciones

El artículo afirma presentar el primer análisis de eQTL de una respuesta inmunitaria humana dinámica in vivo ante un desafío experimental estandarizado que modela la inmunidad celular a Mtb.

• Marco Funcional: El estudio proporciona un marco funcional para interpretar los loci de susceptibilidad a la TB al vincular las señales de GWAS con cambios específicos en la expresión génica en un contexto relevante para la enfermedad.
• Información Mecanística: Identifica dos mecanismos genéticos distintos que regulan las respuestas de las células T:
  1. Específico de antígeno: La variación genética en los niveles de expresión de la clase II de HLA influye en la magnitud de la expansión de las células T reactivas a Mtb.
  2. Independiente de antígeno: Los polimorfismos genéticos que regulan los genes del ciclo celular (por ejemplo, NCAPD3) influyen en la proliferación celular y la expansión clonal de las células T.
• Potencial de Biomarcador: Los hallazgos ofrecen biomarcadores candidatos para la estratificación del riesgo de TB y el desarrollo de vacunas, al resaltar determinantes genéticos específicos de la magnitud y especificidad de la respuesta inmunitaria.
• Validación del Modelo: Los resultados validan la TST como un modelo que captura el repertorio completo de cambios moleculares encontrados en el sitio natural del granuloma de TB pulmonar, extendiendo los hallazgos previos de los modelos ex vivo a las respuestas in vivo multicelulares complejas.

Los autores enfatizan que, si bien el estudio extiende la validación de diversos QTL a respuestas inmunitarias complejas, la verificación experimental de los efectos causales de QTLs y genes específicos sigue siendo un paso necesario para el futuro. El estudio también señala las limitaciones en cuanto al tamaño de la muestra para los análisis estratificados por estado de la enfermedad y la ascendencia mixta de la cohorte, lo cual, aunque supone un desafío para la potencia estadística, mejora la generalizabilidad de los resultados.