An Integrative Transcriptomic Analysis of a Stemness-Associated Gene Set on Ovarian Cancer Recurrence

Resumo Técnico: Uma Análise Transcriptômica Integrativa de um Conjunto de Genes Associados à Stemness no Câncer de Ovário Recorrente

Definição do Problema
O câncer de ovário (CO) permanece uma malignidade ginecológica altamente letal, com uma taxa de sobrevivência de cinco anos inferior a 30% em estágios avançados. Apesar das respostas iniciais favoráveis à quimioterapia baseada em platina, 70–80% dos pacientes apresentam recidiva da doença. Essa recorrência é frequentemente atribuída à persistência de células cancerígenas tronco (CSCs) resistentes à quimioterapia. Embora fatores de transcrição específicos (ex: NANOG, SOX2) e reguladores metabólicos (ex: família ALDH) sejam conhecidos por manter a stemness (capacidade de autorrenovação), há uma necessidade de identificar padrões de expressão gênica consistentes que caracterizem especificamente a transição do estágio primário para o recorrente. O estudo visa determinar se um painel definido de seis genes associados à stemness exibe alterações reprodutíveis na expressão em CO recorrente (RCO) em comparação ao CO primário (PCO), elucidando assim a reprogramação transcricional e metabólica subjacente à recorrência tumoral.

Metodologia
O estudo empregou uma abordagem integrativa combinando dados transcriptômicos públicos de larga escala com validação experimental:

1. Aquisição de Dados e Coortes:

   • Coorte TCGA: Dados de RNA-seq (normalizados por FPKM) de 421 amostras de CO primário (PCO) e 62 de CO recorrente (RCO).
   • Referência GEO: Um conjunto de dados (GSM4684836) de células tronco de câncer de ovário (OCSCs) validadas experimentalmente, derivadas de tumores ovarianos serosos de alto grau (HGSOC).
   • Coorte de Validação Clínica: Amostras de tecido fresco de 8 pacientes com PCO e 8 pacientes com RCO tratados na Faculdade de Medicina da Universidade de Istambul, coletadas antes da primeira linha de quimioterapia e após a recorrência documentada, respectivamente.

2. Painel de Genes Alvo:
   A análise focou em seis genes previamente implicados em stemness, autorrenovação e quimioresistência: ALDH1A1, ALDH1A2, ALDH1A3, NANOG, PROM1 (CD133) e SOX2.

3. Análise Computacional:

   • Pré-processamento: Os valores de FPKM foram filtrados (removendo genes com FPKM < 10 em >80% das amostras) e transformados em log₂ [log₂(FPKM+1)].
   • Expressão Diferencial: Testes estatísticos paramétricos com correção de Taxa de Descoberta Falsa (FDR) foram utilizados para comparar PCO e RCO. A significância foi definida como |log₂ fold change| > 1 e p ajustado < 0,05.
   • Análise de Correlação: Coeficientes de correlação de Pearson foram calculados para avaliar a semelhança dos transcriptomas de PCO e RCO em relação à assinatura de referência de OCSC.
   • Análise de Sobrevivência: Um escore de stemness composto (média do z-score dos seis genes) foi calculado para pacientes do TCGA. Curvas de Kaplan-Meier e regressão de riscos proporcionais de Cox foram utilizadas para avaliar a Sobrevivência Global (SG).

4. Validação Experimental:

   • O RNA total foi isolado das amostras clínicas e convertido em cDNA.
   • A expressão gênica foi quantificada via qRT-PCR utilizando GAPDH como gene de referência.
   • A expressão relativa foi calculada pelo método ΔΔCq, e a concordância entre as mudanças de dobra (fold-changes) de RNA-seq e qRT-PCR foi avaliada via correlação de Spearman.

Principais Resultados

• Expressão Diferencial nos Dados do TCGA:

  • Upregulados em RCO: SOX2 e PROM1 mostraram regulação positiva significativa em tumores recorrentes comparados aos tumores primários (ajustado p < 0,01).
  • Downregulados em RCO: A expressão de ALDH1A1 mostrou-se significativamente diminuída em RCO comparado a PCO (p < 0,01).
  • Sem Mudança Significativa: ALDH1A2, ALDH1A3 e NANOG não apresentaram expressão diferencial estatisticamente significativa na coorte TCGA.

• Correlação com a Assinatura de Célula Tronco:

  • Os tumores primários exibiram uma correlação negativa com a assinatura de OCSC (r = -0,59), indicando um perfil de expressão distinto de células tronco.
  • Os tumores recorrentes mostraram uma correlação quase neutra (r = -0,04), sugerindo uma convergência para um estado transcricional do tipo OCSC durante a recorrência.

• Análise de Sobrevivência:

  • Pacientes com escores de stemness compostos mais elevados mostraram uma tendência para pior Sobrevivência Global; contudo, essa diferença não foi estatisticamente significativa (log-rank p = 0,11).

• Validação por qRT-PCR:

  • Achados Consistentes: SOX2 e PROM1 foram significativamente aumentados na coorte clínica recorrente (p < 0,01 e p < 0,05, respectivamente), corroborando os achados do TCGA.
  • Achados Divergentes: ALDH1A1 apresentou um aumento estatisticamente significativo na coorte de qRT-PCR recorrente (p < 0,05), contrastando com a regulação negativa observada no conjunto de dados do TCGA.
  • Outros Genes: ALDH1A3 estava significativamente elevado na coorte de qRT-PCR (p < 0,05), mas não no TCGA. ALDH1A2 mostrou um aumento não significativo e NANOG mostrou um aumento moderado e não significativo em ambos os conjuntos de dados.
  • Concordância: Uma correlação positiva de fraca a moderada (Spearman r = 0,42) foi observada entre os valores de fold-change das plataformas de RNA-seq e qRT-PCR.

Significância e Alegações
Os autores alegam que a recorrência do câncer de ovário é impulsionada por uma reprogramação transcricional específica centrada na reativação de um subconjunto limitado de reguladores de pluripotência, especificamente o eixo SOX2–PROM1.

• Biomarcadores Robustos: A regulação positiva consistente de SOX2 e PROM1 em ambos os conjuntos de dados de larga escala (RNA-seq) e de pequena escala (qRT-PCR) posiciona-os como biomarcadores confiáveis para a ativação de células tronco associadas à recorrência.
• Heterogeneidade dos Reguladores Metabólicos: O estudo destaca que as isoformas de ALDH exibem regulação heterogênea (ex: ALDH1A1 mostrando tendências opostas em diferentes coortes), sugerindo que sua expressão é altamente influenciada pela exposição à quimioterapia, subtipo tumoral e fatores do microambiente, em vez de servirem como um marcador uniforme de recorrência.
• Insight Mecanístico: A mudança de uma correlação negativa (PCO) para uma correlação neutra (RCO) com a assinatura de OCSC sugere que a recidiva não surge de uma linhagem distinta, mas sim de subpopulações de células tumorais que gradualmente recuperam propriedades transcricionais semelhantes às de células tronco.
• Implicação Terapêutica: Os autores propõem que o alvo da assinatura de célula tronco transcricional SOX2–PROM1 poderia ser uma estratégia eficaz para eliminar populações resistentes e melhorar os desfechos a longo prazo, reconhecendo que os achados atuais são preliminares e requerem validação funcional adicional.

O artigo mantém um tom modesto quanto às aplicações clínicas, observando limitações como o tamanho menor da amostra de RCO no TCGA e o número limitado de amostras clínicas validadas, o que pode afetar a generalização da análise de sobrevivência e a detecção de diferenças transcricionais sutis.