Resumo Técnico: Determinantes Genéticos Comuns da Resposta Imune Humana Dinâmica in vivo ao Mycobacterium tuberculosis

Definição do Problema

Embora a maioria das infecções por Mycobacterium tuberculosis (Mtb) não progrida para a doença tuberculose (TB) ativa, os mecanismos que determinam esses desfechos clínicos divergentes permanecem mal compreendidos. Embora se saiba que as respostas imunes envolvendo células T CD4, macrófagos, IFN-γ e TNF são necessárias para a imunidade protetora, deficiências nessas vias específicas não explicam a maioria dos casos de doença incidente na população geral. Estudos de associação genômica ampla (GWAS) falharam em identificar loci de suscetibilidade à TB estatisticamente robustos e reprodutíveis. Além disso, as análises existentes de loci de traços de expressão quantitativa (eQTL) foram quase exclusivamente limitadas à estimulação ex vivo de populações de células imunes selecionadas, carecendo de validação dentro de respostas imunes in vivo complexas, dinâmicas e multicelulares em um contexto relevante para a doença.

Metodologia

Os autores utilizaram o teste tuberculínico (TST) como um modelo de desafio experimental in vivo padronizado para investigar os determinantes genéticos do hospedeiro nas respostas ao Mtb.

Coorte e Amostragem:

• Participantes: 267 indivíduos (134 homens, 140 mulheres) de diversas etnias, compreendendo duas coortes: aqueles com infecção latente por TB (ILTB) e pacientes com TB ativa (ATB) dentro de 28 dias do diagnóstico/início do tratamento.
• Biópsias: Perfis de sequenciamento de RNA pareados foram obtidos de biópsias de TST no Dia 2 (D2) e no Dia 7 (D7) para participantes com ILTB, e biópsias de D2 para participantes com ATB. Biópsias de injeção de solução salina controle foram utilizadas como linha de base.
• Genotipagem: A genotipagem em todo o genoma foi realizada em 267 indivíduos utilizando o Illumina Infinium Global Diversity Array-8, seguida de imputação pelo painel Haplotype Reference Consortium.

Pipeline Analítico:

1. Transcriptômica: Os dados de RNA-seq foram processados (trimagem, alinhamento ao hg19/GRCh38 e quantificação). A expressão diferencial foi analisada para definir o "transcriptoma de resposta ao TST" (genes significativamente regulados positivamente em comparação com os controles de salina).
2. Mapeamento de eQTL:
   • Cis-eQTL: Modelos lineares mistos (MatrixeQTL/lmerTest) foram usados para identificar associações entre SNPs e expressão gênica dentro de ±1 Mb do sítio de início da transcrição (TSS). As análises foram realizadas em conjuntos de dados combinados de D2/D7 de ILTB e D2 de ATB (415 amostras) para maximizar o poder.
   • Trans-eQTL: As associações foram testadas para SNPs localizados a >5 Mb do TSS do gene, focando nos principais cis-eQTLs identificados em D2 ou D7.
3. Análise Funcional e de Colocalização:
   • Aprendizado Profundo (Deep Learning): Modelos in silico foram utilizados para prever o impacto funcional dos SNPs líderes em sequências regulatórias (promotores/enhancers).
   • Colocalização: A análise de colocalização Bayesiana (coloc) foi realizada para testar se os eQTLs do TST compartilham variantes causais com os sinais de GWAS para TB e doença de Crohn.
   • Herdabilidade: A Regressão de Pontuação de Desvio de LD Estratificada (S-LDSC) foi usada para estimar a proporção de herdabilidade explicada pelos eQTLs de resposta ao TST em doenças imunológicas.
4. Análise do Repertório de Células T: Dados de sequenciamento de TCR de biópsias de D7 foram analisados para identificar metaclones de células T reativas ao Mtb (especificamente restritos por HLA-DRB1) e correlacionar sua expansão com variantes genéticas.

Principais Resultados

1. Variação Interindividual e Descoberta de Cis-eQTL:

• O estudo identificou uma substancial variação interindividual na expressão gênica dentro do transcriptoma do TST, revelando regulação independente de IFN-γ.
• 7.411 genes mostraram associações significativas de SNPs de ação cis (eGenes). Destes, 1.719 eram genes de resposta ao TST.
• O sinal de cis-eQTL mais forte foi mapeado para o ERAP2 (Aminopeptidase 2 do Retículo Endoplasmático).
• A análise de colocalização vinculou a redução da expressão do ERAP2 ao aumento do risco de TB em dados de GWAS, enquanto o aumento da expressão foi associado ao risco de doença de Crohn, destacando um compromisso (trade-off) entre suscetibilidade a doenças autoimunes e infecciosas.

2. Contexto Celular e Vias:

• A colocalização com catálogos públicos de eQTL e dados de célula única indicou que os cis-eQTLs identificados operam principalmente em linhagens de células mieloides (monócitos/macrófagos) e células T.
• Análise de Herdabilidade: Vias relacionadas à apresentação de antígenos da classe II do HLA e ativação de células T mostraram a maior herdabilidade. Especificamente, os níveis de expressão de HLA-DRB1 em TSTs de D2 foram determinados geneticamente e correlacionaram-se com a subsequente expansão de metaclones de células T reativas ao Mtb restritos por HLA-DRB1 em TSTs de D7.
• O estudo identificou tanto mecanismos específicos de antígeno (via QTLs de expressão de HLA-DRB1) quanto mecanismos independentes de antígeno (via regulação do ciclo celular) para a variação da resposta de células T.

3. Trans-eQTL e Regulação do Ciclo Celular:

• A análise de trans-eQTL revelou um programa proliferativo centrado no NCAPD3 (uma subunidade do complexo Condensina II).
• Um único cis-eQTL em NCAPD3 impulsionou a maioria das trans-associações (81,6%), que foram enriquecidas para genes de ciclo celular e mitose.
• O alelo de NCAPD3 foi significativamente associado à frequência de clones de células T expandidos em TSTs de D7, sugerindo um determinante de imunidade de células T geneticamente regulado e independente de antígeno. Essa associação foi validada em um estudo independente de eQTL de célula única de células T estimuladas.

4. Integração com GWAS:

• Os eQTLs de resposta ao TST mostraram enriquecimento significativo para a herdabilidade de várias doenças imunológicas em dados de GWAS.
• Na doença de Crohn, 17 SNPs apresentaram >75% de probabilidade posterior para colocalização com eQTLs de TST.
• Na TB, embora nenhum locus único tenha atingido significância de nível genômico em toda a metanálise, o estudo identificou polimorfismos específicos dentro do locus ERAP2 (por exemplo, rs3733904) associados à redução da expressão e aumento do risco de TB.

Significância e Alegações

O artigo afirma apresentar o primeiro estudo de eQTL de uma resposta imune humana dinâmica in vivo a um desafio experimental padronizado que modela a imunidade mediada por células ao Mtb.

• Estrutura Funcional: O estudo fornece uma estrutura funcional para interpretar os loci de suscetibilidade à TB, vinculando sinais de GWAS a mudanças específicas na expressão gênica em um contexto relevante para a doença.
• Insight Mecanístico: Identifica dois mecanismos genéticos distintos que regulam as respostas de células T:
  1. Específico de antígeno: Variação genética nos níveis de expressão de classe II do HLA influenciando a magnitude da expansão de células T reativas ao Mtb.
  2. Independente de antígeno: Polimorfismos genéticos que regulam genes do ciclo celular (ex: NCAPD3) que influenciam a proliferação celular e a expansão clonal de células T.
• Potencial de Biomarcador: As descobertas oferecem candidatos a biomarcadores para estratificação de risco de TB e desenvolvimento de vacinas, destacando determinantes genéticos específicos da magnitude e especificidade da resposta imune.
• Validação do Modelo: Os resultados validam o TST como um modelo que captura o repertório completo de mudanças moleculares encontradas no local natural do granuloma de TB pulmonar, estendendo as descobertas anteriores de modelos ex vivo para respostas in vivo multicelulares complexas.

Os autores enfatizam que, embora o estudo estenda a validação de diversos QTLs para respostas imunes complexas, a verificação experimental de efeitos causais para QTLs e genes específicos permanece como um passo necessário para o futuro. O estudo também observa limitações quanto ao tamanho da amostra para análises estratificadas por estado de doença e a ancestralidade mista da coorte, o que, embora desafiador para o poder estatístico, aumenta a generalizabilidade dos resultados.