Desarrollo y Validación de una Firma Pronóstica Basada en Reguladores de Metilación de ARN m6A en Adenocarcinoma de Pulmón
El cáncer de pulmón primario sigue siendo el tumor maligno más diagnosticado y la principal causa de muertes relacionadas con cáncer a nivel mundial. Entre los cánceres de pulmón primarios, el cáncer de pulmón de células no pequeñas (CPCNP) representa del 80% al 85% de los casos, y el adenocarcinoma de pulmón (LUAD) constituye aproximadamente del 30% al 35% de todos los casos de cáncer pulmonar primario. La recurrencia del LUAD ha aumentado en muchos países, lo que subraya la necesidad de mejorar las herramientas pronósticas y las estrategias terapéuticas. La metilación de N6-adenosina (m6A), una modificación interna prevalente en el ARN mensajero, ha surgido como un factor crítico en el desarrollo y progresión del cáncer. Este estudio tuvo como objetivo desarrollar y validar una firma pronóstica basada en reguladores de metilación de ARN m6A para pacientes con LUAD.
Papel de los Reguladores de m6A en la Tumorigénesis
La metilación m6A regula la expresión génica en procesos biológicos como la diferenciación de células madre, el ritmo circadiano, la embriogénesis y la respuesta al daño del ADN. En el cáncer, esta modificación influye en la estabilidad, traducción y degradación del ARN mediante tres grupos de proteínas:
- Metiltransferasas (escritoras): WTAP, KIAA1429, RBM15, RBM15B, ZC3H13, METTL3 y METTL14.
- Proteínas de unión (lectoras): IGF2BP1-3, YTHDF1-3, YTHDC1-2, HNRNPC y RBMX.
- Desmetilasas (borradoras): FTO y ALKBH5.
Construcción de la Firma Pronóstica
Se analizaron datos de la cohorte TCGA-LUAD y cuatro conjuntos de datos GEO, incluyendo 500 casos de TCGA y 891 de GEO. Los perfiles de expresión de reguladores de m6A mostraron diferencias significativas entre tejido tumoral y normal. Un análisis de red de interacción proteína-proteína (PPI) identificó seis genes centrales: RBM15, WTAP, METTL3, KIAA1429, METTL14 y YTHDC1. La regresión de Cox univariada seleccionó ocho genes asociados con la supervivencia global (OS), y un modelo de Cox con penalización Lasso generó una puntuación de riesgo basada en IGF2BP3, HNRNPA2B1 y HNRNPC.
Validación del Modelo
Los pacientes se estratificaron en grupos de bajo y alto riesgo según la mediana de la puntuación. El análisis de Kaplan-Meier mostró una OS significativamente peor en el grupo de alto riesgo (p <0,05). El área bajo la curva (AUC) para la firma m6A fue de 0,628, 0,607 y 0,603 a 2, 3 y 5 años, respectivamente, en TCGA-LUAD. Los análisis multivariados confirmaron que el estadio T y la puntuación de riesgo son factores pronósticos independientes (p <0,05).
Correlación Clínica y Validación Externa
Los niveles de IGF2BP3, HNRNPC y HNRNPA2B1 se correlacionaron negativamente con la OS. La validación en cuatro cohortes GEO mostró una OS superior en el grupo de bajo riesgo en tres de ellas (p <0,05), con AUC >0,6. Un modelo integrado que combina la firma m6A y la estadificación TNM demostró mayor precisión pronóstica que cada componente individual.
Análisis Funcional
Los análisis de enriquecimiento funcional (KEGG, GO) y GSEA revelaron vías biológicas asociadas con la firma m6A, incluyendo proliferación celular y respuesta al daño del ADN. La red PPI destacó interacciones clave entre los reguladores de m6A y genes relacionados con el cáncer.
Conclusión
Este estudio desarrolló una firma pronóstica robusta basada en reguladores de m6A para LUAD, validada en múltiples cohortes. La expresión de IGF2BP3, HNRNPA2B1 y HNRNPC se asoció con progresión tumoral y pronóstico adverso. El modelo integrado con estadificación TNM ofrece una herramienta valiosa para estratificar pacientes y guiar terapias personalizadas. Futuras investigaciones deberán explorar el potencial terapéutico de estos reguladores en LUAD.