La progresión de la enfermedad conocida comúnmente como hígado graso metabólico suele describirse mediante categorías histológicas que abarcan desde la esteatosis hasta la esteatohepatitis asociada a disfunción metabólica —MASH—, la fibrosis, la cirrosis y el carcinoma hepatocelular. Un nuevo estudio propone complementar esa clasificación con una trayectoria molecular continua que permita observar cómo cambian los procesos metabólicos, inflamatorios y fibróticos a medida que avanza la enfermedad.
El trabajo, titulado Decoding MASLD Progression: A Molecular Trajectory-Based Framework for Modelling Disease Dynamics, combinó transcriptómica de tejido hepático, análisis de redes moleculares, deconvolución de tipos celulares y proteómica plasmática. El resultado no es un mapa anatómico del hígado, sino una representación estadística de estados moleculares compartidos por diferentes pacientes y alineados con sus características histológicas.
Una trayectoria molecular construida con biopsias hepáticas
Para desarrollar el modelo, los investigadores analizaron inicialmente datos de secuenciación de ARN de biopsias hepáticas de 136 pacientes pertenecientes a dos cohortes previamente publicadas. Después de excluir una muestra atípica, aplicaron un método de ordenamiento pseudotemporal.
Esta técnica no implica que todos los participantes hayan sido seguidos durante años. El procedimiento compara los perfiles transcriptómicos de distintas personas y los organiza según la semejanza de sus estados moleculares, con el objetivo de reconstruir una posible secuencia poblacional de progresión.
La trayectoria mostró una fuerte concordancia con las puntuaciones histológicas de esteatosis, balonización hepatocelular, inflamación, fibrosis y actividad de la enfermedad. Los autores reportaron correlaciones de Pearson entre 0,96 y 1,0 al comparar la posición promedio de los grupos histológicos, aunque no todas las comparaciones alcanzaron significación estadística.
A partir de un modelo de bosque aleatorio, el equipo seleccionó 145 genes asociados con la trayectoria y con diferentes características histopatológicas. Esta firma fue evaluada en la cohorte EPoS, integrada por 168 pacientes con MASLD, y en el conjunto Gubra, que incluyó 26 personas sanas —con peso normal o con obesidad— y 31 pacientes con MASLD o MASH.
La evaluación también incorporó una cohorte japonesa longitudinal de 58 pacientes. En este grupo, los cambios de posición dentro de la trayectoria fueron ampliamente consistentes con la evolución de las biopsias pareadas, particularmente en los casos en los que se observó regresión histológica.
La red molecular reunió más de 7.000 nodos
Los investigadores integraron módulos de coexpresión, factores de transcripción, vías biológicas y características histológicas para construir una red regulatoria de MASLD.
La red principal reunió 7.165 nodos y presentó un enriquecimiento significativo de genes previamente asociados con la enfermedad. Como validación adicional, los autores aplicaron el mismo procedimiento a la cohorte EPoS y obtuvieron una segunda red de 6.732 nodos. Ambas compartieron 5.241 componentes, con un índice de similitud de Jaccard de 0,61.
El análisis mostró que la progresión no se caracteriza por un único patrón lineal. La actividad relacionada con el metabolismo de la glucosa aumentó de manera persistente, mientras que el metabolismo lipídico presentó una dinámica más compleja, con disminución inicial y reactivación parcial en etapas posteriores. En paralelo, la señalización inmunitaria, la hipoxia, los factores de crecimiento y la remodelación de la matriz extracelular adquirieron mayor relevancia conforme avanzó la trayectoria.
Algunas señales moleculares fueron detectadas antes del incremento global de las puntuaciones histológicas. Las vías vinculadas con la matriz extracelular, por ejemplo, se activaron en una ventana anterior al aumento posterior de la fibrosis. El hallazgo sugiere que determinadas lecturas moleculares podrían identificar procesos de progresión antes de que sean plenamente visibles mediante la evaluación histológica convencional, aunque esta posibilidad requiere confirmación prospectiva.
Cambios en la composición celular del hígado
La deconvolución de los perfiles transcriptómicos permitió estimar cómo variaba la contribución relativa de diferentes poblaciones celulares.
La señal de los hepatocitos disminuyó a lo largo de la trayectoria, mientras aumentaron las contribuciones de macrófagos, células linfoides, colangiocitos y fibroblastos. Los autores precisaron que la categoría de fibroblastos representaba principalmente células estrelladas hepáticas, relacionadas con los procesos de fibrosis.
El metabolismo lipídico se asoció fundamentalmente con los hepatocitos; las vías inmunitarias, con macrófagos, neutrófilos y células dendríticas; y la organización de la matriz extracelular, con las poblaciones de fibroblastos. Estas estimaciones reflejan tendencias relativas derivadas de transcriptómica masiva y no mediciones celulares directas de cada biopsia.
Un panel candidato de 57 biomarcadores
La segunda parte del estudio buscó determinar si la trayectoria construida con tejido hepático podía trasladarse hacia marcadores potencialmente medibles en sangre.
Los investigadores partieron de 194 genes cuya expresión había mostrado correlación entre el hígado y el plasma en un trabajo anterior. De ellos, 57 también se encontraban dentro de la red molecular de MASLD y fueron seleccionados como panel candidato.
El modelo fue entrenado para diferenciar fibrosis F0-F2 de fibrosis avanzada F3-F4. En la cohorte de desarrollo alcanzó una exactitud del 86,2 % y un área bajo la curva —AUC— de 0,769. Mediante un procedimiento de selección de variables, el equipo redujo posteriormente la firma a 15 genes, una versión potencialmente más viable para futuras investigaciones.
En la cohorte Fujiwara, el modelo de 57 genes obtuvo un AUC de 0,803, frente a 0,794 para la firma de 15 genes. En EPoS, los valores fueron de 0,86 y 0,85, respectivamente, mientras que FIB-4 alcanzó un AUC de 0,76.
Cuando el clasificador entrenado con transcriptómica se aplicó a un conjunto externo de proteómica plasmática, mantuvo un AUC de 0,83 con el panel completo y de 0,785 con la versión de 15 genes. Los resultados apoyan la posibilidad de trasladar la señal identificada en el hígado hacia mediciones proteicas en plasma, pero todavía no validan una prueba diagnóstica para uso asistencial.
¿Qué implica para la investigación clínica?
Desde la perspectiva de la investigación clínica, el modelo plantea una posible herramienta para estudiar la heterogeneidad de MASLD más allá de las categorías histológicas tradicionales.
En el futuro, una trayectoria molecular validada prospectivamente podría ayudar a estratificar participantes, desarrollar biomarcadores de seguimiento y formular hipótesis terapéuticas ajustadas a las vías predominantes en cada etapa. Sin embargo, estas son aplicaciones potenciales: el estudio no evaluó un medicamento, no incorporó un ensayo clínico prospectivo y no demostró que el panel pueda medir la respuesta a una intervención.
Los propios autores señalan que los biomarcadores todavía no están listos para aplicación clínica. El modelo tampoco establece causalidad ni permite predecir el curso individual de la enfermedad. Además, el predominio de información transversal limita las conclusiones sobre reversibilidad, efectos del tratamiento y modificadores específicos como diabetes, sexo o uso de medicamentos.
El principal aporte del estudio consiste, por tanto, en ofrecer un marco computacional para organizar la progresión molecular de MASLD y orientar nuevas investigaciones sobre estratificación, monitorización y priorización de blancos terapéuticos. La propuesta deberá avanzar hacia estudios longitudinales y prospectivos antes de considerarse una herramienta para la práctica clínica.


