Firma Epigenética de Seis Genes Predice la Supervivencia en Cáncer de Hígado y la Respuesta a Fármacos

El carcinoma hepatocelular (HCC) sigue siendo uno de los cánceres más letales a nivel mundial, con más de 900.000 nuevos casos anuales. A pesar de los avances en cirugía, terapia dirigida e inmunoterapia, los resultados siguen siendo deficientes debido a la heterogeneidad biológica y a la escasez de biomarcadores adecuados para la estratificación de pacientes. Este estudio aborda esa brecha construyendo un modelo pronóstico basado en la biología de la desmetilación del DNA, un mecanismo epigenético cada vez más reconocido como central en la activación de oncogenes y la evasión inmune en el cáncer de hígado.

Los investigadores recopilaron 3.743 genes relacionados con la desmetilación (DRGs, por sus siglas en inglés) de GeneCards y los intersectaron con genes diferencialmente expresados identificados en TCGA-LIHC (346 tumores) y GEO-GSE112790 (183 tumores). El Análisis de Redes de Co-expresión Génica Ponderada (WGCNA) identificó un «módulo azul» con la correlación más fuerte respecto al estado tumoral. La regresión de Cox univariante, la regularización LASSO y el modelado de Cox multivariante destilaron el conjunto de candidatos hasta obtener una firma de seis genes: CEP41, SUB1, CDC20, G6PD, VPS72 y SPINDOC. Se calculó una puntuación de riesgo continua por paciente, estratificándolos en grupos de alto y bajo riesgo según el umbral de la mediana.

El modelo demostró una sólida precisión predictiva con AUCs ≥ 0,70 en los puntos temporales de 1, 3 y 5 años en las cohortes de entrenamiento, prueba interna y validación externa (GSE76427). Los análisis de Kaplan–Meier confirmaron una supervivencia global significativamente peor en los pacientes de alto riesgo (p < 0,001). El análisis SHAP cuantificó la contribución de cada gen a las predicciones de riesgo individuales, mejorando la interpretabilidad del modelo. Un nomograma que integra la puntuación de riesgo con variables clínicas permitió la estimación individualizada de la probabilidad de supervivencia.

Los análisis de enriquecimiento funcional revelaron que los tumores de alto riesgo están enriquecidos en vías del ciclo celular, reparación del DNA y reprogramación metabólica. El perfil inmune mediante ssGSEA mostró que los pacientes de alto riesgo albergan un microambiente inmunosupresor dominado por células T reguladoras y células supresoras de origen mieloide, un patrón consistente con una respuesta deficiente a la inmunoterapia. El análisis de mutaciones somáticas reveló una mayor frecuencia de mutaciones en TP53 en los pacientes de alto riesgo. El modelado de sensibilidad farmacológica (GDSC/pRRophetic) predijo una mayor sensibilidad al inhibidor de aurora quinasa Tozasertib y al inhibidor de BCL-2 Navitoclax en los pacientes de alto riesgo, lo que sugiere vulnerabilidades terapéuticas susceptibles de intervención. G6PD, una enzima clave de la vía de las pentosas fosfato implicada en la reprogramación metabólica, fue validada mediante inmunohistoquímica como significativamente sobreexpresada en el tejido tumoral del HCC, lo que respalda su papel como factor impulsor funcional.

Aunque los hallazgos son convincentes, el modelo se deriva de conjuntos de datos retrospectivos de acceso público y requiere validación clínica prospectiva. Las predicciones de sensibilidad farmacológica son estimaciones computacionales basadas en farmacogenómica de líneas celulares y no han sido evaluadas en pacientes con HCC. No obstante, este marco integrador que vincula la biología de la desmetilación con el pronóstico, la inmunidad tumoral y la respuesta farmacológica representa un avance significativo hacia la oncología de precisión en el cáncer de hígado.