CInvestigating the Mechanism of the Co - regulation of Colorectal Cancer by Tumor - associated Fibroblasts and Disulfidptosis Based on Single - cell RNA Sequencing and Bulk RNA Sequencing

Resumen Técnico: Corregulación del Cáncer Colorrectal por los CAFs y la Disulfidptosis

Planteamiento del Problema
El cáncer colorrectal (CCR) sigue siendo una de las principales causas de mortalidad relacionada con el cáncer a nivel mundial, con una proporción significativa de pacientes diagnosticados en estadios avanzados o que experimentan recurrencia a pesar de una cirugía radical. Si bien se sabe que el microambiente tumoral (TME), particularmente los fibroblastos asociados al cáncer (CAFs), impulsa la progresión y la resistencia a la terapia, los mecanismos específicos mediante los cuales los CAFs interactúan con formas emergentes de muerte celular, como la disulfidptosis, aún no están claros. La disulfidptosis, una muerte celular asociada al metabolismo desencadenada por la acumulación de cistina y el estrés de disulfuro, representa una potencial vulnerabilidad terapéutica en el CCR; sin embargo, no se han establecido biomarcadores fiables que vinculen la actividad de los CAFs con esta vía.

Metodología
El estudio empleó un enfoque integrativo multiómico combinando datos de secuenciación de ARN de célula única (scRNA-seq) y de ARN masivo (RNA-seq):

1. Fuentes de Datos: Los datos de RNA-seq masivo se obtuvieron de The Cancer Genome Atlas (TCGA-CRC) y del conjunto de datos GEO GSE87211. Los datos de scRNA-seq (GSE132257) incluyeron 5 muestras de tejido de CCR y 5 muestras paraneoplásicas. Se curó una lista de 19 genes relacionados con la disulfidptosis (DRGs) a partir de la literatura existente.
2. Análisis de Célula Única: Los datos de scRNA-seq se procesaron utilizando el paquete Seurat para el control de calidad, la normalización y la reducción de la dimensionalidad (PCA, UMAP). La anotación de tipos celulares se realizó mediante el paquete SingleR. La comunicación célula-célula se analizó a través de CellChat para identificar interacciones ligando-receptor. El análisis de pseudotiempo se realizó utilizando Monocle para mapear las trayectorias de diferenciación de los CAFs.
3. Identificación de Genes: Se identificaron genes diferencialmente expresados (DEGs) en los CAFs a partir de scRNA-seq (DEGs1) y en los datos masivos de TCGA-CRC (DEGs2). Se aplicó el Análisis de Red de Coexpresión Génica Ponderada (WGCNA) a los datos de TCGA para identificar módulos de genes correlacionados con las puntuaciones de DRG.
4. Intersección y Modelado: El estudio intersectó los DEGs1, los DEGs2 y los genes del módulo WGCNA para identificar genes corregulados por los CAFs y los DRGs. Los genes pronósticos se cribaron mediante regresión de Cox univariante, seguida de regresión LASSO y modelado de riesgos proporcionales de Cox multivariante para construir una firma de riesgo.
5. Validación y Caracterización: El modelo de riesgo se validó en las cohortes TCGA y GSE87211 mediante curvas ROC, análisis de supervivencia de Kaplan-Meier y AUC dependiente del tiempo. Se realizaron análisis de enriquecimiento funcional (GO, KEGG), infiltración inmunitaria (CIBERSORT) y sensibilidad a fármacos (GDSC).
6. Validación Experimental: Se utilizó una cohorte independiente de 20 pacientes con CCR del Hospital Popular de la Provincia de Guizhou para validar los niveles de expresión proteica de genes clave mediante inmunohistoquímica (IHC) e inmunofluorescencia (IF) en tejidos tumorales frente a tejidos adyacentes normales.

Resultados Clave

• Interacciones Celulares: El análisis de scRNA-seq reveló que los CAFs y las células endoteliales exhiben la interacción recíproca más fuerte en el CCR, mediada principalmente por los ejes de señalización GDF15-TGFBR2 y WNT5A-FZD4. Se encontró que los CAFs se reducen significamente en las muestras de CCR en comparación con el tejido normal, mientras que las células T y los macrófagos se encuentran enriquecidos.
• Intersección de Genes: A partir de 127 DEGs específicos de CAFs, 5,535 DEGs masivos y 3,627 genes del módulo WGCNA, se identificaron 47 genes de intersección. El análisis funcional indicó que estos genes están enriquecidos en la organización de la matriz extracelular (ECM) y en las vías de señalización Wnt.
• Firma Pronóstica: El aprendizaje automático y la regresión de Cox identificaron dos genes pronósticos clave: SRPX (Proteína X con repetición de sushi) y PRELP (Proteína con repetición de leucina rica en prolina y arginina). Se construyó un modelo de puntuación de riesgo: Puntuación de Riesgo = 0.122 × exp(SRPX) + 0.153 × exp(PRELP).
• Estratificación y Supervivencia: Los pacientes estratificados en grupos de alto y bajo riesgo basados en este modelo mostraron diferencias significativas de supervivencia tanto en la cohorte TCGA como en la GSE87211 (P < 0.05). El modelo demostró un rendimiento predictivo robusto con valores de AUC > 0.6 a los 3, 5 y 7 años.
• Correlación Inmunitaria: El grupo de alto riesgo se correlacionó positivamente con la infiltración de macrófagos M0 (r = 0.36 para la puntuación de riesgo; r = 0.44 para PRELP) y negativamente con las células T de memoria CD4+ activadas.
• Trayectoria de Diferenciación: El análisis de pseudotiempo reveló tres vías distintas de diferenciación de CAFs. La expresión de SRPX y PRELP se observó aumentando gradualmente hacia las etapas de diferenciación terminal (CAF_myo, CAF_AP, CAF_adi), lo que sugiere que su regulación al alza está asociada con la maduración de los CAFs.
• Confirmación Experimental: Los análisis de IHC e IF confirmaron que SRPX y PRELP están significativamente regulados al alza en las áreas tumorales de pacientes con CCR en comparación con los tejidos adyacentes normales, encontrándose los niveles más altos de expresión en los CAFs primarios en comparación con los fibroblastos normales y las células cancerosas.
• Implicaciones Terapéuticas: Los pacientes de alto riesgo mostraron una mayor sensibilidad a quimioterapéuticos específicos, incluyendo Bortezomib y Docetaxel, según el análisis de sensibilidad a fármacos.

Significancia y Afirmaciones
Los autores afirman que este estudio proporciona la primera caracterización sistemática de la interacción entre los CAFs y la disulfidptosis en el CCR. Mediante la integración de datos transcriptómicos de célula única y masivos, el estudio construye un nuevo modelo pronóstico basado en genes relacionados con la disulfidptosis derivados de CAFs (SRPX y PRELP) que estratifica eficazmente el riesgo de supervivencia de los pacientes.

El artículo postula que SRPX y PRELP sirven como un vínculo molecular entre los estados funcionales de los CAFs y la disulfidptosis. Los autores sugieren que la regulación al alza de estos genes relacionados con la ECM en los CAFs en diferenciación terminal puede contribuir a la remodelación estromal y la evasión inmunitaria (específicamente al promover la infiltración de macrófagos M0 e inhibir las células T de memoria CD4+). Además, el estudio propone que estos genes pueden hacer que los CAFs o el microambiente tumoral sean vulnerables a la disulfidptosis, ofreciendo una nueva perspectiva para el abordaje terapéutico del compartimento estromal en el CCR. Los hallazgos se presentan como una base para comprender la heterogeneidad pronóstica e identificar posibles dianas terapéuticas, aunque los autores reconocen que los mecanismos moleculares precisos requieren mayor elucidación experimental.