Identification of implications of m6A regulators and autophagy-associated genes for prognosis in ovarian cancer

Resumen Técnico: Identificación de las implicaciones de los reguladores de m6A y los genes asociados a la autofagia para el pronóstico en el cáncer de ovario

Planteamiento del problema
El cáncer de ovario sigue siendo una de las principales causas de mortalidad por cáncer ginecológico, caracterizado por una detección en estadios tardíos, alta heterogeneidad tumoral y quimiorresistencia frecuente. Si bien la modificación de ARN N6-metiladenosina (m6A) y la autofagia son reconocidas individualmente como reguladores críticos en la tumorigénesis y la muerte celular programada, la interacción específica entre los factores reguladores de m6A y los genes asociados a la autofagia en el pronóstico del cáncer de ovario aún no ha sido explorada a fondo. Existe la necesidad de modelos pronósticos robustos que integren estos mecanismos moleculares para guiar estrategias de tratamiento personalizado, incluyendo la quimioterapia y la inmunoterapia.

Metodología
El estudio empleó un enfoque multiómico combinando el análisis bioinformático con la validación experimental:

1. Adquisición y preprocesamiento de datos: Se recuperaron datos de RNA-seq y de variación del número de copias (CNV) de TCGA (379 muestras de cáncer de ovario), GTEx (88 muestras normales) y conjuntos de datos de GEO (GSE13866 para 130 muestras de cáncer de ovario; GSE119168 para perfiles de m6A; GSE184362 para scRNA-seq).
2. Análisis de expresión y mutación: Se evaluó la expresión diferencial de 23 reguladores de m6A (escritores, lectores y borradores) mediante pruebas de suma de rangos de Wilcoxon. Las frecuencias de mutación y CNV se analizaron utilizando "maftools".
3. Correlación con la muerte celular programada: Se utilizó el Análisis de Variación de Conjuntos de Genes (GSVA) para calcular las puntuaciones de 12 vías de muerte celular programada. El análisis de correlación de Pearson vinculó estas puntuaciones con la expresión de los reguladores de m6A, identificando la autofagia como la vía asociada de manera más significativa.
4. Construcción de la puntuación de riesgo:
   • Se fusionaron los genes asociados a la autofagia y los reguladores de m6A.
   • La regresión de Cox univariante identificó 36 genes relacionados con el pronóstico.
   • La cohorte TCGA se dividió (2:1) en conjuntos de entrenamiento y de prueba.
   • La regresión de Cox de Mínimo Cuadrado Absoluto de Selección (LASSO) con validación cruzada de diez pliegues seleccionó 10 genes clave para construir un modelo de puntuación de riesgo.
   • Los pacientes fueron estratificados en grupos de alto riesgo y bajo riesgo basándose en la mediana de las puntuaciones de riesgo.
5. Análisis clínico y terapéutico:
   • Se desarrolló un nomograma que incorpora la edad, el estadio y la puntuación de riesgo para predecir la supervivencia a 1, 3 y 5 años.
   • La sensibilidad a fármacos (IC50) para Docetaxel y Paclitaxel se estimó utilizando el paquete "oncoPredict".
   • El potencial de inmunoterapia se evaluó mediante la Carga Mutacional Tumoral (TMB) y las puntuaciones de Disfunción y Exclusión Inmune Tumoral (TIDE).
6. Validación experimental:
   • RT-qPCR: Se validó la expresión de genes clave en tejidos de cáncer de ovario y tejidos normales adyacentes archivados del Hospital Municipal de Qingdao.
   • scRNA-seq: Los datos de célula única (GSE184362) se procesaron mediante Seurat para anotar tipos celulares y analizar los patrones de expresión de genes clave (específicamente PLK2 y LEPR) a través de diferentes compartimentos celulares.

Resultes clave

• Panorama de los reguladores de m6A: Todos los 23 reguladores de m6A mostraron una expresión diferencial significativa en el cáncer de ovario. YTHDF1, RBM15B, YTHDF2, YTHDF3, RBM15, IGFBP2, ALKBH5, LRPPRC y ZC3H13 estaban sobreexpresados, mientras que otros estaban subexpresados. Las mutaciones fueron raras, pero se observaron alteraciones de CNV (amplificaciones y deleciones).
• Vínculo m6A-Autofagia: Se encontró una fuerte correlación positiva entre los reguladores de m6A y las puntuaciones de autofagia. 18 de los 23 reguladores de m6A mostraron una correlación positiva significativa con la autofagia.
• Modelo pronóstico: La regresión LASSO identificó 10 genes clave: FOXK1, FOXO1, HTR2B, LAMP3, LAMTOR5, LEPR, NPC1, PLEKHF1, PLK2 y RNF5.
  • El modelo de puntuación de riesgo resultante estratificó eficazmente a los pacientes en grupos de alto riesgo y bajo riesgo con resultados de supervivencia significativamente diferentes en el conjunto de entrenamiento, el conjunto de prueba y la cohorte externa GSE138866.
  • El modelo demostró precisión predictiva con valores de AUC de 0.5998 (1 año), 0.5956 (3 años) y 0.7005 (5 años).
  • El nomograma que combina factores clínicos (edad, estadio) y la puntuación de riesgo mostró una buena calibración frente a la supervivencia real.
• Implicaciones terapéuticas:
  • Quimioterapia: El grupo de alto riesgo exhibió una mayor sensibilidad al Docetaxel y al Paclitaxel.
  • Inmunoterapia: El grupo de bajo riesgo mostró un TMB y puntuaciones TIDE más altos, lo que sugiere un mayor potencial de beneficio de la inmunoterapia.
• Validación:
  • La RT-qPCR confirmó que PLEKHF1, FOXO1, LEPR, HTR2B, NPC1, PLK2 y FOXK1 estaban sobreexpresados en los tejidos tumorales, mientras que LAMP3 estaba subexpresado en comparación con los tejidos normales.
  • El análisis de scRNA-seq reveló que PLK2 y LEPR se expresaban significativamente de forma elevada en las células epiteliales dentro de los tejidos cancerosos. El enriquecimiento funcional indicó que PLK2 está asociado con la transición epitelio-mesénquima y las vías de p53, mientras que LEPR está vinculado al metabolismo de ácidos grasos y la hipoxia.

Significancia y afirmaciones
Los autores afirman que este estudio identifica con éxito una nueva firma pronóstica basada en la intersección de la regulación de m6A y la autofagia en el cáncer de ovario. Las contribuciones primarias son:

1. Utilidad pronóstica: La puntuación de riesgo de 10 genes proporciona una herramienta fiable para predecir la supervivencia global y estratificar a los pacientes en distintas categorías de riesgo.
2. Guía terapéutica: El modelo ofrece información para la medicina personalizada, sugiriendo que los pacientes de alto riesgo pueden beneficiarse más de quimioterapias específicas (Docetaxel/Paclitaxel), mientras que los pacientes de bajo riesgo pueden ser mejores candidatos para la inmunoterapia.
3. Perspectiva mecánica: El estudio destaca los roles específicos de PLK2 y LEPR en la tumorigénesis dentro del compartimento epitelial, vinculando la regulación mediada por m6A con las vías de la autofagia y el metabolismo.

El artículo concluye que la integración de las firmas de m6A y la autofagia puede mejorar la predicción del pronóstico y guiar la selección de tratamientos para pacientes con cáncer de ovario, aunque reconoce la necesidad de una mayor validación clínica.