Assinatura de 7 Genes em Macrófagos Prevê Sobrevida e Resposta a Medicamentos no Câncer de Fígado

O carcinoma hepatocelular (HCC) é a forma dominante de câncer primário de fígado e ocupa o terceiro lugar globalmente em mortalidade relacionada ao câncer. A maioria dos pacientes é diagnosticada em estágios avançados, tornando a cirurgia curativa impossível, e as terapias sistêmicas como inibidores de tirosina quinase (TKIs) e inibidores de checkpoint imunológico (ICIs) enfrentam desafios significativos de resistência. Este estudo teve como objetivo descobrir as raízes celulares e moleculares da resistência aos EGFR-TKIs no HCC e converter esse conhecimento em uma ferramenta prognóstica clinicamente aplicável.

Os pesquisadores utilizaram dados de sequenciamento de RNA de célula única (scRNA-seq) de pacientes com HCC, processados com o pipeline Seurat em R. Eles aplicaram quatro algoritmos complementares de enriquecimento de conjuntos de genes (ssGSEA, AUCell, UCell, singscore) por meio do pacote irGSEA para pontuar cada tipo celular quanto à resistência aos EGFR-TKIs, usando um conjunto de genes curado do MSigDB. Os macrófagos apresentaram consistentemente as maiores pontuações de resistência em todos os quatro métodos. Esses macrófagos foram então sub-agrupados, e o sub-grupo com a maior pontuação de resistência teve seus genes marcadores extraídos. Esses marcadores foram cruzados com genes diferencialmente expressos do coorte TCGA-LIHC (tumor vs. tecido normal adjacente via GEPIA2), gerando um conjunto de genes candidatos que foi posteriormente refinado por regressão LASSO Cox para produzir uma assinatura de sete genes.

Os sete genes — SLC41A3, DCAF13, PPM1G, NDC80, FAM83D, FUCA2 e UQCRH — foram utilizados para calcular uma pontuação de risco por paciente. Nos coortes de treinamento, validação e completo do TCGA, pacientes de alto risco apresentaram sobrevida global significativamente pior do que pacientes de baixo risco pela análise de Kaplan-Meier. Esse achado foi replicado de forma independente no conjunto de dados GSE76427. Um nomograma integrando a pontuação de risco com variáveis clínicas (idade, sexo, grau, estadiamento) demonstrou excelente discriminação com AUCs de 0,816 (1 ano), 0,781 (3 anos) e 0,799 (5 anos). A regressão de Cox multivariada confirmou a pontuação de risco como fator prognóstico independente.

Análises de multi-ômica acrescentaram profundidade mecanicista. Pacientes de alto risco apresentaram maior carga mutacional tumoral (TMB), padrões distintos de variação do número de cópias e perfis diferenciais de metilação do DNA em comparação com pacientes de baixo risco. O algoritmo TIDE previu que pacientes de baixo risco responderiam melhor ao bloqueio de checkpoint imunológico, apresentando menores pontuações de disfunção e exclusão. Por outro lado, análises farmacogenômicas utilizando os bancos de dados GDSC e CTRP via oncoPredict sugeriram que pacientes de alto risco eram mais sensíveis ao Navitoclax (um inibidor de BCL-2) e ao Sorafenib. A validação em conjuntos de dados de scRNA-seq de imunoterapia (GSE206325 e Liu_HCC) demonstrou ainda expressão diferencial dos genes do modelo entre respondedores e não respondedores.

O estudo fornece uma ferramenta prognóstica biologicamente embasada e validada em múltiplos coortes, enraizada no compartimento de macrófagos tumorais imunossupressores. Ele conecta os mecanismos de resistência terapêutica diretamente à predição de sobrevida no nível do paciente e à sensibilidade ao tratamento, oferecendo um arcabouço para personalizar o manejo do HCC. As ressalvas incluem a natureza retrospectiva e computacional do estudo, a dependência de conjuntos de dados públicos sem validação clínica prospectiva e a ausência de experimentos funcionais que confirmem os papéis causais de cada gene na resistência.