Artigo de investigação

Mineração em Câncer de Endométrio com base no banco de dados TCGA e construção de modelos prognósticos

DOI:

10.3791/69650

16 de janeiro de 2026

Neste artigo

Resumo

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Este protocolo delineia uma abordagem computacional para identificar genes associados ao metabolismo da nicotinamida no câncer de endométrio, utilizando conjuntos de dados TCGA e GEO. Para facilitar a descoberta de biomarcadores e aprimorar a estimativa de risco no carcinoma endometrial do corpo uterino, esta pesquisa descreve a análise de expressão distinta, aprimoramento de rotas, desenho prognóstico do modelo e procedimentos de avaliação.

Resumo

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

O câncer de endométrio (CE) está entre os tumores ginecológicos menos prevalentes no mundo, com aumento da incidência devido ao envelhecimento e à obesidade. A metástase dos gânglios linfáticos continua comum no Carcinoma Endometrial do Corpo Uterino (UCEC), exigindo novos biomarcadores prognósticos para orientar o tratamento. Nesta pesquisa, as informações do UCEC do Atlas do Genoma do Câncer (TCGA) foram analisadas, e os resultados foram validados usando o Gene Expression Omnibus (GEO). Dezoito fatores genéticos diferencialmente expressos associados ao metabolismo da nicotinamida (NMRDEGs) foram identificados. A Análise de Enriquecimento de Conjuntos Gênicos (GSEA) indicou seu papel no estresse oxidativo, hipóxia, glicólise e processos de apoptose. A regressão de Cox univariada identificou seis genes-chave (AURKA, CDKN3, FOXM1, CDKN2A, TK1 e CDK1), utilizados para avançar em uma estrutura de previsão de riscos. A análise de interação proteína-proteína (PPI) revelou genes hubs adicionais como CDK2, CCNA2, TP53 e FOXM1. Os seis genes-chave apresentaram forte valor prognóstico, e o modelo de risco do estudo pode orientar decisões clínicas. O metabolismo da nicotinamida foi encontrado significativamente ligado à progressão da CE. Este estudo oferece novas percepções sobre o papel do metabolismo da nicotinamida na CE e sugere possíveis caminhos para avanços no tratamento.

Introdução

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

O câncer de endométrio (CE) é comum em pacientesdo sexo feminino 1. A CE pode ser classificada no tipo dependente de hormônios e no tipo2 independente de hormônios. Em 2020, as atualizações das diretrizes para o manejo de pacientes com EC incluíram detecção e tipagem molecular de marcadores moleculares, influenciando significativamente o tratamento adjuvante pós-operatório e o prognóstico clínico para EC3. A patogênese da CE é multifacetada, envolvendo mutações ou perda de vários oncogenes e/ou genes supressores de tumores, juntamente com anomalias em múltiplas vias de sinalização que podem influenciar a progressão e o desfecho da doença. Buscar novos biomarcadores e identificar alvos de tratamento eficazes são fundamentais. Essa busca aprofundar a compreensão da patogênese da CE e dos potenciais biomarcadores prognósticos, estabelecendo uma base para a detecção precoce e a terapia direcionada.

A niacinamida se assemelha à niacina e desempenha um papel na absorção de gordura e na quebra do glicogênio. É essencial para o corpo humano como um componente importante da coenzima. O metabolismo da nicotinamida abrange os processos pelos quais o corpo humano absorve, transforma e utiliza a nicotinamida.

O metabolismo da niacinamida começa com a absorção. Pode ser adquirido por meio da ingestão de alimentos ou conversão de aminoácidos dentro do corpo e é absorvido pelas células epiteliais da mucosa intestinal. Depois, ele se transforma em formas mais complexas, cruciais em muitos processos bioquímicos. Além disso, a nicotinamida é convertível em outros compostos, incluindo o metabólito 1-metilnicotinamida por meio da nicotinamida adenina dinucleotídea metiltransferase (NMT). Essas transformações podem influenciar o início e a progressão de várias doenças metabólicas. Conforme relatado, a suplementação com niacinamida traz vários efeitos benéficos, incluindo a prevenção e tratamento detumores. Estudos recentes confirmaram que ele pode reduzir o nível de fosforilação da STAT3Y705 em células de câncerde pulmão 5. Um estudo descobriu que o uso de Niacinamida sozinha pode inibir tumores invasivos de membrana não melanoma em pacientes de altorisco 6. Kourtzidis et al.7 descobriram que suplementar com niacinamida pode inibir o ganho de peso em ratos, aumentar a atividade da transcetolose (TK) e da NA+K+- ATPase dos glóbulos vermelhos, além de afetar o metabolismo energético no corpo. Um estudo anterior reconheceu a nicotinamida N-metiltransferase (NNMT), uma metiltransferase, como uma questão importante na atividade dos fibroblastos associadosao câncer 8. Além disso, foi demonstrado que a nicotinamida fosforibosil transferase contribui para a resistência ao venetoclax em compartimentos estalos de leucemia mieloide críticarecaída 9. Portanto, o objetivo era identificar de forma abrangente genes relacionados ao metabolismo do NAD+ (NMRGs) ligados ao risco de CE.

Acesso restrito. Inicie sessão ou comece um teste para visualizar este conteúdo.

Protocolo

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Este estudo utilizou dados clínicos e transcriptômicos públicos e desidentificados do The Cancer Genome Atlas e do Gene Expression Omnibus. Todos os estudos contribuintes tinham aprovação prévia do conselho de revisão institucional e consentimento informado. Como apenas a análise secundária dos dados anonimizados foi realizada, não foi necessária aprovação ética adicional. Os bancos de dados e softwares utilizados estão listados na Tabela de Materiais.

1. Download de dados

A pesquisa utilizou o conjunto de dados EC (TCGA-Uterine Corpus Endometrial Carcinoma (TCGA-UCEC)), que consiste em 589 amostras, incluindo 554 amostras de tecido tumoral de pacientes UCEC (grupo UCEC) e dados de sequenciamento de 35 tecidos normais adjacentes (grupo normal). O banco de dados Xena da UCSC foi utilizado para recuperar os dados clínicoscorrespondentes 11, excluindo aqueles que não possuíam informações clínicas completas. No total, 577 amostras com dados clínicos estavam disponíveis para análise. Informações detalhadas de referência são fornecidas na Tabela 1.

Conjuntos de dados adicionais relacionados ao EC, GSE115810 e GSE6367812 foram baixados usando o pacoteGEOquery 13. O conjunto de dados GSE115810 e GSE63678 foram combinados para criar os Conjuntos de Dados Combinados para análise posterior (Tabela 2).

Genes relacionados ao metabolismo (NMRGs) do NAD+ estão situados após o GeneCardsrecord14 e a literaturarelevante 15. Usando "metabolismo da niacinamida" como termo de busca no GeneCards, identificaram 345 NMRGs com pontuações de relevância acima de 4. A combinação e remoção de duplicados dos 42 NMRGs encontrados na literatura compilou um total de 371 NMRGs (Tabela Suplementar 1). Os dados clínicos foram adquiridos como arquivos fenotípicos em tsv; os dados eram baixados no formato HTSeq-FPKM. Amostras excluídas tinham mais de 20% de seus dados clínicos faltando. O FPKM foi transformadologaritaritar-2 e convertido em TPM (transcrições por milhão). Os IDs das sondas foram mapeados para símbolos gênicos para conjuntos de dados GEO, e sondas duplicadas foram mediadas.

2. Genes diferencialmente expressos do metabolismo da nicotinamida

A pesquisa começou aplicando o conjunto Rsva16 para eliminar posses conjuntas após os conjuntos de dados GSE115810 e GSE63678, resultando em um conjunto mútuo contendo 31 amostras EC (UCEC) e 8 amostras normais adjacentes. Em seguida, a utilização do conjunto limma17 para realizar exame de expressão gênica de discrepância no conjunto de dados TCGA-UCEC.

A interseção dos DEGs da análise TCGA-UCEC com 337 NMRGs para identificar os DEGs conectados ao metabolismo da nicotinamida. Isso produziu uma lista de genes diferencialmente expressos pelo metabolismo da niacinamida (NMRDEGs), que foram ilustrados por uma ilustração de Venn. Os resultados do exame de aparência de discrepância foram ilustrados pelo pacote ggplot2 R18, enquanto um mapa de calor dos NMRDEGs foi gerado usando o conjunto de pheatmap19. A variância entre conjuntos de dados é eliminada por correção em lote usando ComBat (Bayes empírico). O modelo linear empírico de Bayes foi empregado na análise DEG. Limiares de expressão diferenciais aplicados explicitamente:

|log 2FC| ≥ 1

FDR menos de 0,05.

A lista do NMRG só se cruzava com DEGs que atendiam a ambos os requisitos. Plots de vulcões e mapas de calor feitos com ggplot2 e pheatmap.

3. Exame de função dos NMRDEGs (GO), aprimoramento do caminho (KEGG)

Os exames de aprimoramento do GO20 e KEGG21foram concluídos pelo conjunto 22 do clusterProfiler. Para ambas as análises, limiares de significância foram reconhecidos em p. ajustar< 0,05 e FDR (valor q) < 0,25. A pesquisa também integrou os valores logFC na análise de enriquecimento, representando os resultados em diagramas circulares e de cordas. Limiares para significância de enriquecimento: valor p ajustado < 0,05 FDR < 0,25 (valor q). O Cluster Profiler foi utilizado para estudos GO e KEGG. A direcionalidade genética é exibida usando gráficos de corda e círculo que incorporam dados de mudançalogarítmica 2 vezes.

4. Análise de Enriquecimento de Conjuntos Gênicos (GSEA)

O tipo de conjuntos de fatores hereditários que mais contribuíram para o fenótipo pode ser identificado usando GSEA23. Para essa análise, o conjunto de dados TCGA-UCEC foi classificado com base nos valores logFC, e um exame de melhoria foi realizado usando o pacote clusterProfiler . Restrições principais incluíam um valor semente de 2022 e 10.000 permutações. O conjunto genético MSigDB "c2.all.v2022.1.Hs.symbols.gmt" foi utilizadoem 24. As vias mais enriquecidas, incluindo genes induzidos por hipóxia de Manalo, senescência induzida por estresse oxidativo, glicólise e apoptose, foram visualizadas usando um gráfico de montanha. Antes do GSEA, os genes eram ordenados por log2 vezes de variação. 10.000 permutações foram usadas na análise. c2.all.v2022.1.Hs.symbols.gmt é a coleção do MSigDB que foi utilizada. Para reprodutibilidade, uma semente aleatória fixa (2022) foi empregada. Caminhos significativos foram aqueles com p < 0,05 e q < 0,25.

5. Construção do modelo Cox e exame prognóstico relacionado

Para determinar o valor preditivo dos genes diferencialmente expressos ligados ao metabolismo da nicotinamida (NMRDEGs) no carcinoma endometrial (UCEC), os pesquisadores utilizaram a análise de regressão de Cox univariada para primeiro classificar os fatores hereditários dos candidatos; aqueles com proporção de perigo (HR) > 1 e < de valor p de 0,1 foram considerados adequados para o quadro multivariado de riscos relativos de Cox.

Critérios para seleção univariada de Cox: p < 0,10 e HR > 1. Valores de expressãonormalizados com 2-TPM foram empregados no modelo multivariado de Cox. Uma combinação linear de coeficientes de Cox × expressão gênica é usada para determinar a pontuação de risco. Probabilidades de OS de 1, 3 e 5 anos foram empregadas na calibração de nomogramas. Valores de AUC de 1, 3 e 5 anos foram empregados em ROC dependente do tempo. A abordagem surv_cutpoint max-statistic foi usada para encontrar valores de corte de sobrevivência. Tanto as análises KM quanto as ROC usaram os mesmos limiares.

Um nomografo foi construído a partir do modelo multivariado de Cox para avaliar sua precisão ou capacidade preditiva e calcular as chances de existência total de 1, 3 e 5 anos. Arcos de padronização são empregados para avaliar a estabilidade entre os futuros resultados e os resultados reais, e a análise da curva de decisão (DCA) foi utilizada para medir a eficacia médica daestrutura 25.

Os níveis de expressão de mRNA foram determinados como valores normalizados de registros por milhão (TPM) transformados em log₂ por meio do pacote DESeq2. As TPMs levaram em conta a complexidade do sequenciamento e a medição de genes para fornecer estimativas robustas e não tendenciosas dos níveis de expressão entre amostras.

Usando os coeficientes do modelo Cox multivariado, a avaliação prognóstica de risco de cada paciente foi determinada da seguinte forma:

riskScore = Σi Coeficiente (genei) *Expressão de mRNA (genei) (1)

Arcos de existência Kaplan-Meier (KM) foram preparados para avaliar a resistência geral de agrupamentos de alto e baixo risco criados com base em determinadas classificações de risco. Arcos de características operacionais dependentes do tempo do receptor (ROC) foram produzidos na avaliação da rotina de estruturas em períodos de 1, 3 e 5anos 26,27.

Para categorizar a expressão gênica por coleções de alta e baixa expressão para estratificação de sobrevivência, utiliza-se a utilização do papel de ponto surv_cut após o pacote R do survminer. Essa função determina o maior valor de corte maximizando a estatística padronizada de classificação logarítmica, fornecendo um ponto de corte imparcial e estatisticamente ótimo.

Os valores de corte obtidos para cada gene prognóstico são indicados nas curvas ROC como linhas tracejadas. A pesquisa aplicou os mesmos limiares para todas as análises de sobrevivência e ROC.

O RCGA RNA-SEQ foi baixado no formato HTSeq-FPKM; os dados clínicos eram importados como arquivos de fenótipo TSV; O FPKM foi convertido para TPM e transformado em log₂; Conjuntos de dados GEO foram mapeados de IDs de sonda para símbolos de genes usando anotações de plataforma; sondas duplicadas foram médias para um único valor do gene; Amostras com mais de 20% de informação clínica ausente foram excluídas; ComBat (Bayes empírico) foi usado para correção em lote para conjuntos de dados GSE; PCA e boxplots foram usados para verificar se a correção por lote foi bem-sucedida. Normalização de TPM usando técnicas padrão de transformação de expressões; correção por lote usando ComBat com origem do conjunto de dados como variável por lote; expressão diferencial calculada usando modelagem linear limma (matriz de desenho tumoral vs. normal); listas geneais ranqueadas geradas a partir de alterações de fold log₂ para entrada GSEA; e regressões de Cox univariadas e multivariadas realizadas usando ferramentas de análise de sobrevivência

6. Análise de variação de conjuntos gênicos (GSVA)

OGSVA 28 foi utilizado para medir o aumento de processos entre os clusters. No conjunto de dados TCGA-UCEC, 50 vias características foram enriquecidas, com 41 apresentando alterações importantes entre as assembleias binárias. O GSVA foi usado com conjuntos genéticos característicos para obter atividade de vias por amostra; Dados de interação com proteínas STRING foram importados para o Cytoscape; o algoritmo MCC foi usado para identificar genes hub; as pontuações de risco foram calculadas como a soma dos valores de expressão gênica multiplicada por seus coeficientes de Cox; curvas ROC dependentes do tempo foram geradas usando rotinas ROC de tempo de sobrevivência. Para cada amostra, a GSVA calculou as pontuações de enriquecimento em nível de caminho. O teste de soma de ranks de Wilcoxon é usado para avaliar diferenças na atividade do caminho característico. Dos cinquenta caminhos de assinatura, quarenta e um foram significativamente diferentes (p ajustado < 0,05).

7. Sistema de interação proteína-proteína (PPI)

Um sistema PPI contendo os genes importantes (AURKA, CDKN3, FOXM1, CDKN2A, TK1 e CDK1) foi criado usando o arquivoSTRING 29 e um limiar de valor de comunicação de 0,70, indicando alta confiança. Essa rede foi criada usando o Cytoscape30, destacando interações que podem desempenhar papéis cruciais na patogênese do UCEC. O método31 da Centralidade Máxima da Clique (MCC) foi útil para classificar o fator hereditário criado nas pontuações de interação deles na rede. A sequência proteica superior 10 com as pontuações de interface mais altas foi reconhecida, incluindo CDK2, CDK4, CCNA2, CCNB1, CCNE1, CDK1, TP53 e FOXM1. Esses genes foram posteriormente analisados quanto ao seu envolvimento em processos biológicos críticos. A plataformaGeneMANIA 32 também foi usada para prever interações proteicas adicionais e fornecer um contexto mais amplo para os papéis dos genes-chave na progressão do UCEC. Limiar para a pontuação de confiança STRING: >0,70 (alta confiança). O Cytoscape exibe a rede. A técnica de Centralidade Máxima da Clique (MCC) é usada para classificar genes centrais. A classificação MCC foi usada para identificar os genes que mais interagem (CDK2, CCNA2, TP53, etc.). Previsões adicionais de interação são feitas usando o GeneMANIA.

8. Roteiro tecnológico

O fluxo de trabalho geral e os métodos empregados neste estudo são resumidos no roteiro tecnológico exibido na Figura 1. Este roteiro descreve os passos desde a aquisição de conjuntos de dados e análise de expressões diferenciais até a construção de modelos prognósticos e análises de enriquecimento.

9. Análise estatística

O processamento de dados e a estimativa estatística foram realizados usando o programa R (v4.3.0). O teste Mann-Whitney U ou teste t do Estudante Independente foi usado para comparações em dois grupos; o teste de Kruskal-Wallis foi empregado para três ou mais montagens. Os dados descritivos foram avaliados usando qui-quadrado ou teste exato de Fisher. Além disso, foram realizadas análises de correlação de Spearman e de sobrevivência de Kaplan-Meier; p < 0,05 foi considerado significativo.

Aplicação de testes estatísticos de acordo com a distribuição de dados: Dados normais usando o teste t do aluno. Teste Mann-Whitney U para dados que não são normais. Kruskal-Wallis testam para mais de três grupos. Exato de Fisher e qui-quadrado para dados categóricos. A significância estatística é definida como p < 0,05.

A confiabilidade dos dados é mantida pelos checkpoints de pré-tratamento, onde os boxplots devem mostrar variância consistente de expressão entre amostras, e os gráficos PCA devem demonstrar a ausência de clusters específicos de lote após o ajuste ComBat. Mapas de calor que demonstrem agrupamento tumor-normal e gráficos vulcânicos que ilustram claramente a regulação genética para cima e para baixo são necessários para validação da DEG. Para o modelo de projeção de Cox, os gráficos de calibração devem corresponder à sobrevivência prevista e real, os valores de AUC do ROC devem ser maiores que 0,65, e os arcos KM devem mostrar uma diferença substancial de sobrevivência. Diferentes atividades de rotas entre grupos de risco devem ser demonstradas por análise GSVA, em linha com mecanismos estabelecidos como rotas de expansão ou ciclo celular. Para verificar a resiliência da rede, nós altamente acoplados na rede PPI devem aparecer centralmente, e genes hub determinados pelo MCC devem corresponder a reguladores fisiologicamente significativos.

Acesso restrito. Inicie sessão ou comece um teste para visualizar este conteúdo.

Resultados

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Fusão da correção do conjunto de dados GEO

A pesquisa utilizou o conjunto SVA R para precisar posses em lote nos conjuntos combinados de GSE115810 e GSE63678, produzindo um conjunto unificado denominado Conjuntos de Dados Combinados. Os resultados da eliminação de resultados em grupo foram verificados usando caixas de distribuição (Figuras 2A, B) e gráficos PCA (Figura 2C, D), que dem...

Acesso restrito. Inicie sessão ou comece um teste para visualizar este conteúdo.

Discussão

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

O papel do metabolismo energético no microambiente das células cancerígenas tem atraído considerável atenção recentemente. Pesquisas mostraram que genes envolvidos no metabolismo da nicotinamida estão intimamente ligados à progressão de várias malignidades. Uma vez absorvida pelo corpo, a nicotinamida é convertida em NAD e NADP, ambos fundamentais em inúmeros processos relacionados ao câncer. Consequentemente, o metabolismo do NAD está emergindo como um objetivo terapêutico promissor par...

Acesso restrito. Inicie sessão ou comece um teste para visualizar este conteúdo.

Divulgações

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Os autores não têm nada a revelar.

Agradecimentos

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Não aplicável.

Acesso restrito. Inicie sessão ou comece um teste para visualizar este conteúdo.

Materiais

Lista de materiais utilizados neste artigo
NomeEmpresaNúmero de catálogoComentários
Pacote clusterProfilerBiocondutorv4.8.0Análise de enriquecimento GO e KEGG
CitoscapeConsórcio Cytoscapev3.9.1Visualização de redes PPI
Pacote DESeq2Biocondutorv1.40.0Normalização e análise de expressão
Banco de dados GeneCardsSuíte GeneCardsN/AFonte dos genes relacionados ao metabolismo da nicotinamida
Plataforma GeneMANIAgenemania.orgN/APrevisão da função e interação do gênio
Conjuntos de dados GEO (GSE115810, GSE63678)Omnibus de Expressão Gênica (GEO)N/AConjuntos de dados públicos de expressão gênica para EC
Pacote ggplot2CRANv3.4.2Visualização de dados
Pacote LimmaBiocondutorv3.54.0 (exemplo)Análise diferencial de expressão
Pacote pheatmapCRANv1.0.12Plotagem de mapas de calor
Software R (v4.3.0)Fundação R para Computação EstatísticaVersão 4.3.0Plataforma de análise de dados
Banco de dados STRINGCORDA (string-db.org)v11.5Análise de interação proteína-proteína
Pacote survminerCRANv0.4.9Análise de sobrevivência
Pacote SVABiocondutorv3.46.0Remoção do efeito em lote
Conjunto de dados TCGA-UCECO Atlas do Genoma do Câncer (TCGA)N/ATranscriptoma e dados clínicos de pacientes com UCEC
UCSC XenaUniversidade da Califórnia, Santa CruzN/APlataforma para download de dados clínicos/genômicos

Referências

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,
  1. Bray, F., et al. Global cancer statistics 2018: Globocan estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 68 (6), 394-424 (2018).
  2. Ni, J., et al. All human genes of the uteroglobin family are localized on chromosome 11q12.2 and form a dense cluster. Ann N Y Acad Sci. 923, 25-42 (2000).
  3. Concin, N., et al. Esgo/estro/esp guidelines for the management of patients with endometrial carcinoma. Virchows Arch. 478 (2), 153-190 (2021).
  4. Nikas, I. P., Paschou, S. A., Ryu, H. S. The role of Nicotinamide in cancer chemoprevention and therapy. Biomolecules. 10 (3), 477(2020).
  5. Wang, W., Karamanlidis, G., Tian, R. Decreased NAD activates STAT3 and integrin pathways to drive epithelial-mesenchymal transition. Mol Cell Proteomics. 17 (10), 2005-2017 (2018).
  6. Chen, A. C., et al. A phase 3 randomized trial of Nicotinamide for skin-cancer chemoprevention. N Engl J Med. 373 (17), 1618-1626 (2015).
  7. Kourtzidis, I. A., et al. Nicotinamide riboside supplementation dysregulates redox and energy metabolism in rats: Implications for exercise performance. Exp Physiol. 103 (10), 1357-1366 (2018).
  8. Eckert, M. A., et al. Proteomics reveals NNMT as a master metabolic regulator of cancer-associated fibroblasts. Nature. 569 (7758), 723-728 (2019).
  9. Jones, C. L., et al. Nicotinamide metabolism mediates resistance to venetoclax in relapsed acute myeloid leukemia stem cells. Cell Stem Cell. 27 (5), 748-764 (2020).
  10. Colaprico, A., et al. TCGAbiolinks: An R/Bioconductor package for integrative analysis of TCGA data. Nucleic Acids Res. 44 (8), e71(2016).
  11. Goldman, M. J., et al. Visualizing and interpreting cancer genomics data via the Xena platform. Nat Biotechnol. 38 (6), 675-678 (2020).
  12. Zhao, D., Liu, J., Wang, C., Zhang, Y., Chen, L. Diagnosis and prediction of endometrial carcinoma using machine learning and artificial neural networks based on public databases. Genes (Basel). 13 (6), 935(2022).
  13. Davis, S., Meltzer, P. S. GEOquery: A bridge between the Gene Expression Omnibus (GEO) and Bioconductor. Bioinformatics. 23 (14), 1846-1847 (2007).
  14. Stelzer, G., et al. The GeneCards suite: From gene data mining to disease genome sequence analyses. Curr Protoc Bioinformatics. 54, 1.30.1-1.30.33 (2016).
  15. Cui, H., Liu, Y., Jiang, X., Li, H., Zhang, L., Wang, Z. Comprehensive analysis of Nicotinamide metabolism-related signature for predicting prognosis and immunotherapy response in breast cancer. Front Immunol. 14, 1145552(2023).
  16. Leek, J. T., Johnson, W. E., Parker, H. S., Jaffe, A. E., Storey, J. D. The SVA package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics. 28 (6), 882-883 (2012).
  17. Huang, J., Wang, S., Li, Y., Liu, X., Zhang, Y., Chen, Z. Proposing a novel molecular subtyping scheme for predicting distant recurrence-free survival in breast cancer post-neoadjuvant chemotherapy with close correlation to metabolism and senescence. Front Endocrinol (Lausanne). 14, 1265520(2023).
  18. Ito, K., Murphy, D. Application of ggplot2 to pharmacometric graphics. CPT Pharmacometrics Syst Pharmacol. 2 (10), e79(2013).
  19. Li, G. M., Zhang, C. L., Rui, R. P., Sun, B., Guo, W. Bioinformatics analysis of common differential genes of coronary artery disease and ischemic cardiomyopathy. Eur Rev Med Pharmacol Sci. 22 (11), 3553-3569 (2018).
  20. Mi, H., Muruganujan, A., Ebert, D., Huang, X., Thomas, P. D. PANTHER version 14: More genomes, a new PANTHER GO-slim and improvements in enrichment analysis tools. Nucleic Acids Res. 47 (D1), D419-D426 (2019).
  21. Ogata, H., et al. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 27 (1), 29-34 (1999).
  22. Yu, G., Wang, L. G., Han, Y., He, Q. Y. ClusterProfiler: An R package for comparing biological themes among gene clusters. OMICS. 16 (5), 284-287 (2012).
  23. Subramanian, A., et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 102 (43), 15545-15550 (2005).
  24. Liberzon, A., Birger, C., Thorvaldsdóttir, H., Ghandi, M., Mesirov, J. P., Tamayo, P. The Molecular Signatures Database hallmark gene set collection. Cell Syst. 1 (6), 417-425 (2015).
  25. Tataranni, T., Piccoli, C. Dichloroacetate (DCA) and cancer: An overview towards clinical applications. Oxid Med Cell Longev. 2019, 8201079(2019).
  26. Barakat, A., Mittal, A., Ricketts, D., Rogers, B. A. Understanding survival analysis: actuarial life tables and the Kaplan-Meier plot. Br J Hosp Med (Lond). 80 (11), 642-646 (2019).
  27. Heagerty, P. J., Lumley, T., Pepe, M. S. Time-dependent ROC curves for censored survival data and a diagnostic marker. Biometrics. 56 (2), 337-344 (2000).
  28. Hänzelmann, S., Castelo, R., Guinney, J. GSVA: Gene set variation analysis for microarray and RNA-seq data. BMC Bioinformatics. 14, 7(2013).
  29. Wakonig, B., Auersperg, A. M. I., O'Hara, M. String-pulling in the Goffin's cockatoo (Cacatua goffiniana). Learn Behav. 49 (1), 124-136 (2021).
  30. Shannon, P., et al. Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res. 13 (11), 2498-2504 (2003).
  31. Boughorbel, S., Jarray, F., El-Anbari, M. Optimal classifier for imbalanced data using Matthews correlation coefficient metric. PLoS One. 12 (6), e0177678(2017).
  32. Franz, M., et al. GeneMANIA update 2018. Nucleic Acids Res. 46 (W1), W60-W64 (2018).
  33. Chen, A., Jing, W., Qiu, J., Zhang, R. Prediction of cervical cancer outcome by identifying and validating a NAD⁺ metabolism-derived gene signature. J Pers Med. 12 (12), 2031(2022).
  34. Lin, L., et al. Identification of NAD⁺ metabolism-derived gene signatures in ovarian cancer prognosis and immunotherapy. Front Genet. 13, 905238(2022).
  35. Li, H., et al. circABHD2 inhibits malignant progression of endometrial cancer by regulating the NAD⁺/NAMPT metabolism axis. Mol Biotechnol. 67 (7), 2644-2659 (2025).
  36. Hu, D., et al. Comprehensive analysis of a NAD⁺ metabolism-derived gene signature to predict the prognosis and immune landscape in endometrial cancer. Biomol Biomed. 24 (2), 346-359 (2024).
  37. Deng, B., et al. AURKA emerges as a vulnerable target for KEAP1-deficient non-small cell lung cancer by activation of asparagine synthesis. Cell Death Dis. 15 (3), 233(2024).
  38. Zhang, C., Shen, Q., Gao, M., Li, J., Pang, B. The role of cyclin-dependent kinase inhibitor 3 (CDKN3) in promoting human tumors: literature review and pan-cancer analysis. Heliyon. 10 (4), e26061(2024).
  39. Bai, X., Xiong, J., Li, L., Yu, C., Sun, C. Suppression of hypoxia-induced CAV1 autophagic degradation enhances nanoalbumin-paclitaxel transcytosis and improves therapeutic activity in pancreatic cancer. Eur J Pharmacol. 969, 176431(2024).
  40. Zhao, Z., et al. Activation of the FOXM1/ASF1B/PRDX3 axis confers hyperproliferative and antioxidative stress reactivity to gastric cancer. Cancer Lett. 589, 216796(2024).
  41. Yang, C., Plum, P. S., Gockel, I., Thieme, R. Pan-cancer analysis and in vitro validation of the oncogenic and prognostic roles of AURKA in human cancers. Front Oncol. 13, 1186101(2023).
  42. Shen, X., et al. YTHDC1-dependent m6A modification-modulated FOXM1 promotes glycolysis and tumor progression through CENPA in triple-negative breast cancer. Cancer Sci. 115 (6), 1881-1895 (2024).
  43. Jiang, J., et al. HNRNPA2B1-mediated m6A modification of FOXM1 promotes drug resistance and inhibits ferroptosis in endometrial cancer via regulation of LCN2. Funct Integr Genomics. 24 (1), 3(2023).
  44. Yue, C., et al. Autophagy-related risk signature based on CDKN2A to facilitate survival prediction of patients with endometrial cancer. J Gene Med. 26 (1), e3648(2024).
  45. Fu, Y., et al. Discovery of potential biomarkers for early diagnosis of endometrial cancer via integrating metabolomics and transcriptomics. Comput Biol Med. 173, 108327(2024).
  46. Omar, I. S., et al. Aberrant upregulation of CDK1 contributes to medroxyprogesterone acetate resistance in cancer-associated fibroblasts of the endometrium. Biochem Biophys Res Commun. 628, 133-140 (2022).
  47. Zhang, B., et al. A novel glycolysis-related gene signature for predicting the prognosis of multiple myeloma. Front Cell Dev Biol. 11, 1198949(2023).
  48. Huang, H., et al. A novel five-gene metabolism-related risk signature for predicting prognosis and immune infiltration in endometrial cancer: A TCGA data mining study. Comput Biol Med. 155, 106632(2023).
  49. Mao, X., et al. High glucose levels promote glycolysis and cholesterol synthesis via ERRα and suppress the autophagy-lysosomal pathway in endometrial cancer. Cell Death Dis. 16 (1), 182(2025).

Acesso restrito. Inicie sessão ou comece um teste para visualizar este conteúdo.

Reimpressões e permissões

Solicitar permissão para reutilizar o texto ou as figuras deste artigo JoVE

Solicitar permissão

Etiquetas

Endometrial CancerTCGA DatabasePrognostic ModelsUterine Corpus Endometrial CarcinomaNicotinamide MetabolismGene Set EnrichmentDifferential Gene ExpressionCox RegressionProtein Interaction AnalysisPrognostic Biomarkers

Artigos relacionados