SEARCH · PubMed Health

Results for “univariate and multivariate Cox”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 recordsLinked to original sources

PubMed2025

Construction of a prognostic model for gastric cancer based on immune infiltration and microenvironment, and exploration of MEF2C gene function.

BACKGROUND: Advanced gastric cancer (GC) exhibits a high recurrence rate and a dismal prognosis. Myocyte enhancer factor 2c (MEF2C) was found to contribute to the development of various types of cancer. Therefore, our aim is to develop a prognostic model that predicts the prognosis of GC patients and initially explore the role of MEF2C in immunotherapy for GC. METHODS: Transcriptome sequence data of GC was obtained from The Cancer Genome Atlas (TCGA), the Gene Expression Omnibus (GEO) and PRJEB25780 cohort for subsequent immune infiltration analysis, immune microenvironment analysis, consensus clustering analysis and feature selection for definition and classification of gene M and N. Principal component analysis (PCA) modeling was performed based on gene M and N for the calculation of immune checkpoint inhibitor (ICI) Score. Then, a Nomogram was constructed and evaluated for predicting the prognosis of GC patients, based on univariate and multivariate Cox regression. Functional enrichment analysis was performed to initially investigate the potential biological mechanisms. Through Genomics of Drug Sensitivity in Cancer (GDSC) dataset, the estimated IC50 values of several chemotherapeutic drugs were calculated. Tumor-related transcription factors (TFs) were retrieved from the Cistrome Cancer database and utilized our model to screen these TFs, and weighted correlation network analysis (WGCNA) was performed to identify transcription factors strongly associated with immunotherapy in GC. Finally, 10 patients with advanced GC were enrolled from Sun Yat-sen University Cancer Center, including paired tumor tissues, paracancerous tissues and peritoneal metastases, for preparing sequencing library, in order to perform external validation. RESULTS: Lower ICI Score was correlated with improved prognosis in both the training and validation cohorts. First, lower mutant-allele tumor heterogeneity (MATH) was associated with lower ICI Score, and those GC patients with lower MATH and lower ICI Score had the best prognosis. Second, regardless of the T or N staging, the low ICI Score group had significantly higher overall survival (OS) compared to the high ICI Score group. For its mechanisms, consistently, for Camptothecin, Doxorubicin, Mitomycin, Docetaxel, Cisplatin, Vinblastine, Sorafenib and Paclitaxel, all of the IC50 values were significantly lower in the low ICI Score group compared to the high ICI Score group. As a result, based on univariate and multivariate Cox regression, ICI Score was considered to be an independent prognostic factor for GC. And our Nomogram showed good agreement between predicted and actual probabilities. Based on CIBERSORT deconvolution analysis, there was difference of immune cell composition found between high and low ICI Score groups, probably affecting the efficacy of immunotherapy. Then, MEF2C, a tumor-related transcription factor, was screened out by WGCNA analysis. Higher MEF2C expression is significantly correlated with a worse OS. Moreover, its higher expression is also negatively correlated with tumor mutation burden (TMB) and microsatellite instability (MSI), but positively correlated with several immunosuppressive molecules, indicating MEF2C may exert its influence on tumor development by upregulating immunosuppressive molecules. Finally, based on transcriptome sequencing data on 10 paired tumor tissues from Sun Yat-sen University Cancer Center, MEF2C expression was significantly lower in paracancerous tissues compared to tumor tissues and peritoneal metastases, and it was also lower in tumor tissues compared to peritoneal metastases, indicating a potential positive association between MEF2C expression and tumor invasiveness. CONCLUSIONS: Our prognostic model can effectively predict outcomes and facilitate stratification GC patients, offering valuable insights for clinical decision-making. The identified transcription factor MEF2C can serve as a biomarker for assessing the efficacy of immunotherapy for GC.

Humans↗

PubMed2025

CS Ratio is an immune-related prognostic biomarker for cervical cancer.

BACKGROUND: The tumor microenvironment (TME) plays a crucial role in cancer progression but its complex structure significant variability among patients present considerable challenges for research. Recent studies have demonstrated that macrophage polarization states defined by the expression levels of CXCL9 SPP1 (CS Ratio) are more prognostically relevant than traditional M1/M2 markers. The CS polarization state reflects a highly coordinated network of pro-tumor anti-tumor variables offering a simplified yet effective immune response indicator for the complex TME. The CS Ratio has been shown to correlate with the abundance of anti-tumor immune cells the gene expression programs of tumor-infiltrating cells responses to immunotherapy. Cervical cancer, one of the most common gynecological malignancies, still faces limited therapeutic options. CXCL9, a member of the CXC chemokine family, plays a critical role in immune regulation, inflammation, tumor growth, angiogenesis, and metastasis. Similarly, SPP1, a cytokine, influences immune-related pathways by regulating molecules such as interferon-γ and interleukin-12. However, no studies have systematically investigated the role of the CS Ratio in cervical cancer or its relationship with immunotherapy characteristics. Research in this area could provide critical insights into the role and clinical potential of the CS Ratio in cervical cancer and related tumors. METHODS: The expression ratio of CXCL9 to SPP1 was analyzed in cervical cancer patients using data from the Gene Expression Omnibus (GEO) database, which revealed significant differences. Data for cervical cancer patients were obtained from The Cancer Genome Atlas (TCGA) database. The optimal cutoff value for the CS Ratio was determined using the maxstat package in R, and Kaplan-Meier (KM) survival curves were constructed. Patients were categorized into High and Low groups based on the median CS Ratio. Immune scores were analyzed, and immune cell infiltration was assessed using CIBERSORT. Differences in the CS Ratio were evaluated across patients with varying pathological T stages and FIGO stages. Additionally, receiver operating characteristic (ROC) analysis was performed using the pROC package in R to calculate the area under the curve (AUC). Univariate and multivariate Cox regression analyses were performed to evaluate the potential of the CS Ratio as an independent prognostic factor in cervical cancer. A Cox regression-based nomogram integrating four key features was subsequently developed for the TCGA-CESC cohort. Nomogram performance was assessed using calibration curves and ROC analysis. RESULTS: The CS Ratio was significantly lower in cervical cancer patients compared to normal controls (P < 0.05). KM survival curves indicated that patients in the CS High group exhibited better prognoses. Immune score analysis revealed significantly higher immune scores (P < 0.05) and lower tumor purity (P < 0.05)in the CS High group compared to the Low group. CIBERSORT analysis revealed significantly higher proportions of CD8+ T cells (P < 0.05) and M1 macrophages (P < 0.05), and a significantly lower proportion of M2 macrophages (P < 0.05), in the CS High group compared to the Low group. The CS Ratio significantly decreased with advancing FIGO stage (P < 0.05). Both univariate (P < 0.05) and multivariate Cox regression analyses (P < 0.05) confirmed the CS Ratio as an independent prognostic factor. ROC analysis demonstrated that the CS Ratio had higher AUC values for predicting 1-year (AUC=0.69), 3-year (AUC=0.66), and 5-year OS (AUC=0.68) than CXCL9 or SPP1 alone. The Cox regression-based nomogram integrating four key features demonstrated predictive capability for 1-, 3-, and 5-year OS in CESC patients (Concordance Index = 0.751; 95% CI: 0.678-0.824; p = 1.50Í10-11). Significant survival differences were observed between the high-risk and low-risk groups based on the nomogram score. ROC analysis yielded high AUC values for survival prediction: 0.85 (95% CI: 0.94-0.75) at 1-year, 0.74 (95% CI:0.84-0.64) at 3-year, and 0.72 (95% CI:0.84-0.61) at 5-year. CONCLUSION: The CS Ratio may serve as a more effective prognostic biomarker for cervical cancer patients.

CXCL9↗

PubMed2025

Association Between HLA-DRB1 Serotype and HLA-DQB1 Allele Mismatches and Acute Rejection in Kidney Transplantation.

The purpose of this single-center case-control study was to investigate the association between HLA serotype mismatch (MM), compared to other HLA MM modalities, and the occurrence of acute rejection (AR) within the first year after deceased donor kidney transplantation. The study included 198 transplants in 99 pairs of recipients of kidneys from the same donor, where one recipient experienced AR and the other survived the first year without AR. Donors and recipients were typed with NGS for 11 HLA loci at high resolution. HLA MM categories included allele groups, alleles, serotypes, amino acids, EMMA, eplet and PIRCHE-II. Additionally, we investigated Cytomegalovirus LIL peptide (CMV LIL) MM. Recipients with AR presented higher frequencies of pre-transplant HLA-ABDR DSA (20.2% vs. 6.1%, p = 0.005) and CMV LIL MM (24.2% vs. 10.1%, p = 0.01). Univariate and multivariate Cox proportional hazards regression for matched-pair analyses were used to test the association between HLA MM and AR. Univariate analyses indicated significant association with DRB1 ST, HLA-DQB1 AG, HLA-DQB1 AL, EMMA C, EMMA DQB1, Eplet ABC and Eplet DQ MM. Different models were tested in multivariate analyses, all including pre-transplant HLA-ABDR DSA and CMV LIL MM. The models were compared using the Akaike Information Criterion (AIC). The best estimate for AR prediction (AIC = 97.6) was the model that included pre-transplant HLA-ABDR DSA (HR = 11.97; p = 0.003), CMV LIL MM (HR = 367.2; p < 0.001), HLA-DRB1 serotype MM (9.65; p = 0.002) and HLA-DQB1 allele MM (HR = 3.54; p = 0.033). In conclusion, this original report demonstrates an association between the HLA-DRB1 serotype MM and AR, highlighting that serotypes are clinically relevant.

Humans↗

PubMed2026

A cuproptosis-related lncRNAs-based risk signature for predicting prognosis and immune status in glioma.

BACKGROUND: Glioma is one of the most prevalent primary malignant brain tumors, characterized by poor prognosis and limited treatment options. Recent studies have identified cuproptosis, a novel copper-dependent form of regulated cell death, as a critical mechanism involved in tumor progression. However, the role of cuproptosis-related long non-coding RNAs (lncRNAs) in glioma remains not fully clarified. This study aimed to develop and validate a prognostic model based on cuproptosis-associated lncRNAs to predict patient outcomes and guide individualizing therapeutic strategies. METHODS: Transcriptomic profiles and clinical data were obtained from The Cancer Genome Atlas (TCGA), The Genotype-Tissue Expression (GTEx), and the Chinese Glioma Genome Atlas (CGGA) databases. Cuproptosis -related prognostic lncRNAs were filtered via univariate and multivariate Cox and Least absolute shrinkage and selection operator (LASSO) regression analyses, which were selected to establish a prognostic model for glioma. Samples were divided into high- and low-risk groups, and the predictive performance of the prognostic model was evaluated based on receiver operating characteristic (ROC) curves, Kaplan-Meier (K-M) survival curves, and a nomogram. In addition, immune cell infiltration, tumor mutational burden (TMB), immunophenoscore (IPS), Tumor Immune Dysfunction and Exclusion (TIDE) and drug sensitivity were analyzed. Expression levels of selected lncRNAs and proteins were validated using quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR) and Western blotting. RESULTS: An 11-lncRNA signature associated with cuproptosis was established, and the risk score derived from this model was identified as an independent prognostic factor for glioma. The model exhibited excellent predictive ability, with area under the curve (AUC) values of 0.880, 0.913, and 0.866 for 1-, 3-, and 5-year survival, respectively. Higher TMB, immune checkpoint expression, and IPS were observed in the high-risk group and no significant difference was observed in TIDE between risk groups. Drug sensitivity analysis identified TPCA-1, KIN001-135, and ispinesib mesylate as potential therapeutic agents. Expression validation in glioma cells further supported the biological relevance of the selected lncRNAs. CONCLUSIONS: This cuproptosis-related lncRNA-based signature demonstrates strong prognostic value and may serve as a promising tool for glioma risk stratification and personalized treatment selection.

Glioma↗

PubMed2026

Coagulation activation is associated with genomic-instability-related features in TP53-mutated AML and MDS: routine laboratory patterns beyond classical disseminated intravascular coagulation.

BACKGROUND: Disseminated intravascular coagulation (DIC) is a serious complication of acute myeloid leukemia (AML) associated with poor prognosis. In TP53-mutated AML and myelodysplastic syndrome (MDS), however, the classical ISTH criteria rarely identify overt DIC, although bleeding and thrombotic complications are well documented in acute leukaemia. We hypothesized that these patients exhibit a lower-grade, subclinical coagulation activation that is associated with the underlying genomic-instability-related features of TP53-mutant disease. METHODS: We retrospectively analyzed 107 consecutive patients with TP53-mutated AML (n = 52) or high-risk MDS (MDS, n = 55), median age 65 years, diagnosed and initially evaluated at our centre between 2018 and 2025. Seven routine coagulation markers and 46 co-mutated genes were evaluated for associations with overall survival (OS) using univariate and multivariable Cox regression, continuous dose-response modeling, and unsupervised k-means clustering. Internal validity was assessed by 1000 bootstrap resamples. RESULTS: Overt DIC according to ISTH criteria was rare (15%). Subclinical activation was common: 50% of patients had a D-dimer ≥1 μg/mL, 41% a fibrinogen ≥4 g/L, and 29% an INR ≥1.2. In univariate analysis, D-dimer, fibrinogen, INR, prothrombin time, and activated partial thromboplastin time were each associated with OS (HR 1.33-1.38 per SD; all p < 0.05). Complex karyotype correlated with higher D-dimer (median 1.39 vs. 0.60 μg/mL, p = 0.022) and fibrinogen (3.91 vs. 2.53 g/L, p = 0.007), while TP53 variant allele frequency (VAF) showed modest positive correlations with D-dimer (ρ = 0.21), INR (ρ = 0.27), and PT (ρ = 0.27; all p < 0.05). Clustering identified three coagulation phenotypes: Silent (51%), Thrombo-inflammatory (31%), and Consumption-like (18%), showing a graded but statistically non-significant gradient in molecular features and a stepwise decline in median OS (14, 10 and 8 months; log-rank p = 0.041). After adjustment for complex karyotype, TP53 VAF, and favorable co-mutation count, the Consumption-like phenotype was associated with a non-significant increased risk (HR 1.83, 95% CI 0.92-3.65, p = 0.084), whereas favorable co-mutation pathways remained independently protective (HR 0.56, 95% CI 0.35-0.90, p = 0.016). CONCLUSION: In TP53-mutated AML/MDS, coagulation activation intensity is associated with the degree of genomic instability. The three phenotypes may add biological resolution beyond classical DIC and cytogenetic risk groups, but represent laboratory patterns rather than validated bleeding or thrombosis prediction tools. However, after accounting for genomic features, phenotypes were not independent predictors of outcome, with complex karyotype, TP53 VAF, and favorable co-mutation count driving prognosis. Because treatment intensity and other clinical confounders were not available, these survival associations are hypothesis-generating. Coagulation profiling remains inexpensive, widely accessible, and offers a practical window into disease biology that warrants prospective validation.

TP53↗

PubMed2020

Exploration of predictive and prognostic alternative splicing signatures in lung adenocarcinoma using machine learning methods.

BACKGROUND: Alternative splicing (AS) plays critical roles in generating protein diversity and complexity. Dysregulation of AS underlies the initiation and progression of tumors. Machine learning approaches have emerged as efficient tools to identify promising biomarkers. It is meaningful to explore pivotal AS events (ASEs) to deepen understanding and improve prognostic assessments of lung adenocarcinoma (LUAD) via machine learning algorithms. METHOD: RNA sequencing data and AS data were extracted from The Cancer Genome Atlas (TCGA) database and TCGA SpliceSeq database. Using several machine learning methods, we identified 24 pairs of LUAD-related ASEs implicated in splicing switches and a random forest-based classifiers for identifying lymph node metastasis (LNM) consisting of 12 ASEs. Furthermore, we identified key prognosis-related ASEs and established a 16-ASE-based prognostic model to predict overall survival for LUAD patients using Cox regression model, random survival forest analysis, and forward selection model. Bioinformatics analyses were also applied to identify underlying mechanisms and associated upstream splicing factors (SFs). RESULTS: Each pair of ASEs was spliced from the same parent gene, and exhibited perfect inverse intrapair correlation (correlation coefficient = - 1). The 12-ASE-based classifier showed robust ability to evaluate LNM status of LUAD patients with the area under the receiver operating characteristic (ROC) curve (AUC) more than 0.7 in fivefold cross-validation. The prognostic model performed well at 1, 3, 5, and 10 years in both the training cohort and internal test cohort. Univariate and multivariate Cox regression indicated the prognostic model could be used as an independent prognostic factor for patients with LUAD. Further analysis revealed correlations between the prognostic model and American Joint Committee on Cancer stage, T stage, N stage, and living status. The splicing network constructed of survival-related SFs and ASEs depicts regulatory relationships between them. CONCLUSION: In summary, our study provides insight into LUAD researches and managements based on these AS biomarkers.

Adenocarcinoma of Lung↗

PubMed2025

Expression and prognosis of CXCL13 in uterine corpus endometrial carcinoma based on bioinformatics analysis.

OBJECTIVE: The biological significance of the chemokine ligand C-X-C motif chemokine ligand 13 (CXCL13) may play a significant role in the pathogenesis of uterine corpus endometrial carcinoma (UCEC). This study aims to identify and verify CXCL13 with predictive value for prognosis in UCEC. METHODS: CXCL13 mRNA expression differences were analyzed using R software in three independent datasets: one each from The Cancer Genome Atlas (TCGA) and two from the Gene Expression Omnibus (GEO), namely GSE17025 and GSE106191. The correlation between CXCL13 expression and prognosis was evaluated by Kaplan-Meier analysis. Univariate and multivariate Cox analyses were utilized to construct a prognostic nomogram. Tumor Immune Estimation Resource (TIMER) and the Tumor and Immune System Interaction Database (TISIDB) were employed to assess the relationship between CXCL13 and tumor immune infiltration. Coexpressed genes with CXCL13 were identified by the Spearman correlation analysis. A CXCL13 protein-protein interaction (PPI) network was constructed with the STRING website tool and hub genes were screened out. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genome (KEGG) analyses were performed with the "clusterProfiler" R package. Gene set enrichment analysis (GSEA) was used to identify underlying biological mechanisms. A drug-gene interaction network was constructed in the Comparative Toxicogenomics Database (CTD). RESULTS: High CXCL13 mRNA expression were validated in UCEC in the above three independent datasets. High CXCL13 expression was associated with favorable prognosis in UCEC. A nomogram for predicting the 1-, 3-, and 5-year survival probability in UCEC was construct based on CXCL13 expression and other clinical parameters. The use of Spearman correlation indicated certain correlation between CXCL13 and immune cells and immune checkpoint (ICP) genes. Seven hub genes were upregulated in UCEC, namely CXCL9, IFNG, CXCL10, CXCL11, GBP5, CCL18, and GZMB. The expression and prognostic relevance of CXCL9, IFNG, GBP5, and GZMB were in accordance with CXCL13. The main biological processes enriched were cytokine-cytokine receptor interaction and chemokine signaling pathway. CONCLUSIONS: The above comprehensive analyses suggest that CXCL13 may serve as a potential prognostic biomarker for UCEC, specifically for early-stage UCEC.

CXCL13↗

PubMed2026

Integrative analysis of the roles and prognostic value of RNA-binding proteins in papillary renal cell carcinoma.

RNA-binding proteins (RBPs) serve essential roles in various cancer types, but their functions in papillary renal cell carcinoma (pRCC) have not been elucidated to date. In our work, differentially expressed RBPs in pRCC were identified after acquisition of RNA-sequencing and clinical data related to pRCC from The Cancer Genome Atlas database(TCGA). Functional enrichment analysis and protein interaction network analysis, along with univariate and multivariate Cox regression analyses, were performed to uncover potential biological effects of the identified RBPs and screen the hub RBPs for pRCC prognosis. We identified 251 up-regulated and 129 down-regulated RBPs, and filtered out seven hub RBPs, namely, SRSF8, CD3EAP, HBS1L, ELAC2, MRPL34, NOP2 and IGF2BP2, for their prognostic relevance. A prognostic risk score model for overall survival of pRCC patients was constructed based on the seven hub RBPs. Further analysis showed that the low-risk group had higher survival rate than the high-risk group in both training and validation cohorts. The predictive accuracy was verified in the Human Protein Atlas database.In addition, we introduced the GSE15641 dataset from the Gene Expression Omnibus (GEO) database for independent external validation, and confirmed the expression levels of HBS1L, MRPL34 and IGF2BP2 through real-time quantitative PCR (RT-qPCR) and Western blotting (WB) using human renal tubular epithelial cell line HK-2 and human papillary renal cell carcinoma cell line Caki-2. In pRCC, CD3EAP was significantly elevated, while ELAC2, IGF2BP2, MRPL34, SRSF8 and HBS1L were significantly reduced. There was no significant difference between tumor and normal tissues in NOP2 expression. Risk score and tumor grade were independent prognostic factors associated with overall survival. In addition, we established a nomogram based on the seven prognostic RBPs to help predict overall survival at 1-3 years. In conclusion, seven differentially expressed hub RBPs were identified as potential prognostic biomarkers for pRCC. Our prognostic model might serve as a support for better treatment decision-making. Our work could provide potential new ideas for diagnosis and research on targeted drugs for pRCC.

Bioinformatics↗

PubMed2026

Integrative analysis identifies a glycosylation-related lncRNA signature associated with prognosis in kidney renal clear cell carcinoma.

BACKGROUND: Glycosylation and long non-coding RNAs (lncRNAs) play critical roles in tumor progression. However, the prognostic significance of glycosylation-related lncRNAs (GRLncs) in kidney renal clear cell carcinoma (KIRC) remains largely unclear. This study aimed to identify prognostic GRLncs and construct a predictive model for KIRC prognosis. METHODS: Transcriptomic and clinical data of KIRC patients were analyzed to identify GRLncs associated with overall survival (OS). A prognostic model was constructed based on selected GRLncs, and its predictive performance was evaluated using Kaplan-Meier (KM) survival analysis, receiver operating characteristic (ROC) curves, and univariate and multivariate Cox regression analyses. Patients were stratified into high- and low-risk groups according to the median risk score, and internal validation was performed using training and testing cohorts to assess the stability of the model. Tumor microenvironment characteristics, immune checkpoint expression, immunotherapy response, and drug sensitivity were further analyzed. In addition, the expression of three signature lncRNAs was validated by real-time quantitative polymerase chain reaction (RT-qPCR) in 10 paired KIRC tumor and adjacent normal tissues. Functional roles of selected lncRNAs were investigated using antisense oligonucleotides (ASOs)-mediated knockdown in KIRC cell lines, followed by Cell Counting Kit 8 (CCK-8), 5-ethynyl-2'-deoxyuridine (EdU) incorporation, colony formation, and migration assays. RESULTS: Five GRLncs (AC093278.2, EPB41L4A-DT, DLGAP1-AS2, AC084876.1, and AC005261.3) were identified and used to construct a prognostic model. AC093278.2 and EPB41L4A-DT were protective factors, whereas DLGAP1-AS2, AC084876.1, and AC005261.3 were risk factors. KM analysis on GRLncs-based risk score stratification revealed patients in the high-risk group had significantly poorer OS than those in the low-risk group. ROC analysis and Cox regression demonstrated that the GRLnc-based risk score served as an independent predictor of KIRC prognosis and exhibited favorable predictive performance compared with conventional clinical variables. High- and low-risk groups also exhibited distinct immune microenvironment characteristics, immune checkpoint expression patterns, and predicted drug sensitivities. RT-qPCR detected significant downregulation of protective factor-EPB41L4A-DT in KIRC tissues, while risk factors-DLGAP1-AS2 and AC084876.1 showed expression trends consistent with their predicted risk attributes. Functional experiments further revealed that knockdown of DLGAP1-AS2 and AC084876.1 suppressed proliferation and migration of KIRC cells, whereas knockdown of EPB41L4A-DT promoted these processes, supporting the biological relevance of these three signature lncRNAs. CONCLUSIONS: This study establishes a novel prognostic model based on five GRLncs that showed promising performance in The Cancer Genome Atlas (TCGA)-based analyses of KIRC. The combined clinical expression analysis and functional validation of three constituent GRLncs (DLGAP1-AS2, EPB41L4A-DT, and AC084876.1) supports the biological plausibility of the model and suggest that GRLncs may serve as potential prognostic biomarkers and therapeutic targets for KIRC.

Kidney renal clear cell carcinoma (KIRC)↗

PubMed2026

A per- and polyfluoroalkyl substances-based gene signature links prognosis to immune landscapes in thyroid cancer.

BACKGROUND: Thyroid cancer (THCA) is the most common endocrine malignancy with a rising global incidence and significant heterogeneity. Although per- and polyfluoroalkyl substances (PFAS) exposure is linked to thyroid dysfunction, the prognostic value of per- and polyfluoroalkyl substances-related genes (PFASRGs) and their role in the tumor immune microenvironment (TME) remain poorly understood. This study aims to systematically screen key PFASRGs and evaluate their prognostic value as biomarkers for THCA. METHODS: Utilizing The Cancer Genome Atlas (TCGA)-THCA transcriptomic data and PFASRGs, we constructed a prognostic model through differential expression analysis, univariate and multivariate Cox regression analyses, and the least absolute shrinkage and selection operator (LASSO). The model's robustness was validated using receiver operating characteristic (ROC) curves, Kaplan-Meier analysis, and clinical nomograms. Furthermore, the TME, immunotherapy response, and drug sensitivities were systematically evaluated. Distinct molecular landscapes were characterized by stratifying the cohort via unsupervised consensus clustering analysis. RESULTS: The eight-gene prognostic model demonstrated robust performance, with area under the curve (AUC) values exceeding 0.85 across all validation cohorts. High-risk patients exhibited significantly shorter overall survival and an "inflamed" TME characterized by high immune scores and checkpoint expression. In contrast, the therapeutic efficacy of anti-programmed death-ligand 1 (PD-L1) agents was more pronounced in the low-risk category, as evidenced by a superior objective response. Furthermore, distinct molecular subtypes and risk-specific sensitivities to targeted agents, such as sorafenib and sunitinib, were identified, highlighting the model's clinical utility for personalized treatment. CONCLUSIONS: We established a novel THCA prognostic framework based on eight PFASRGs. This model exhibits superior performance in risk stratification, effectively distinguishing cohorts with divergent clinical trajectories, unique immune microenvironment features, and varied therapeutic responses. Our findings provide a powerful predictive tool for refining prognostic evaluation and facilitating the implementation of personalized management strategies for THCA patients.

Per- and polyfluoroalkyl substances-related genes ↗

PubMed2026

Development and internal validation of a six-gene prognostic model based on galactose metabolism for overall survival in lung adenocarcinoma.

BACKGROUND: Lung cancer remains a leading cause of cancer incidence and mortality globally. Metabolic reprogramming promotes tumor progression and shapes an immunosuppressive tumor microenvironment. Galactose metabolism is involved in multiple malignancies, but its prognostic value in lung adenocarcinoma (LUAD) remains unclear. This study aimed to develop and internally validate a galactose metabolism-related multigene prognostic model for LUAD. METHODS: A retrospective prognostic model development and internal validation study was performed using RNA sequencing (RNA-seq) and clinical data from 585 LUAD patients in The Cancer Genome Atlas (TCGA). Differential expression, functional enrichment, univariate and multivariate Cox regression were applied to construct a prognostic gene signature. Internal validation was performed using bootstrap resampling. Model performance was evaluated by time-dependent receiver operating characteristic (ROC), C-index, calibration, and Kaplan-Meier analysis. Associations between the model and immune infiltration, immunotherapy responsiveness, and tumor stemness were also analyzed. RESULTS: A six-gene prognostic model (GALT, GANC, PGM1, GALM, B4GALT1, PGM2) was developed. The model showed good discrimination with 1-, 3-, and 5-year area under the curve (AUC) values of 0.719, 0.693, and 0.684, respectively. The low-risk group exhibited significantly longer survival, increased antitumor immune infiltration (CD8+ T cells, M1 macrophages, activated CD4+ memory T cells), higher expression of T cell proliferation-related genes, lower immune checkpoint expression, better predicted immunotherapy response, and lower tumor stemness compared with the high-risk group. CONCLUSIONS: We developed and internally validated a six-gene prognostic model for LUAD based on galactose metabolism. The model shows moderate prognostic performance and is associated with antitumor immunity and tumor stemness. It may be used for prognostic risk stratification and to guide personalized immunotherapy in LUAD.

Galactose metabolism↗

PubMed2026

Crucial role of telomere maintenance-related genes in survival prediction and subtype identification in colorectal cancer.

BACKGROUND: Telomere maintenance-related genes (TMRGs) are implicated in Colorectal cancer (CRC) development, but their prognostic value and clinical relevance remain insufficiently explored. This study aims to develop a TMRG-based prognostic model and elucidate its clinical utility in CRC management. METHODS: The Cancer Genome Atlas database was utilized to download RNA-seq data from 638 CRC and 51 control samples. Differential expressed genes were screened and intersected with 2086 TMRGs, resulting in the identification of 976 TMRGs. Through univariate and multivariate Cox regression analysis, a prognostic model comprising three telomere maintenance-related biomarkers (PDE1B, TFAP2B, and HSPA1A) was developed and validated using an external dataset. By integrating the model risk score with clinical features, a nomogram was constructed to predict the survival outcomes of CRC patients. Additionally, an in-depth investigation of the immuno-infiltration, functional variation and drug sensitivity analysis were performed in two risk subgroups defined by the prognostic model. Finally, the functional significance of PDE1B in CRC cell lines was investigated through MTT assays, cell colony formation assays, transwell assays and flow cytometry. RESULTS: A total of 976 DE-TMRGs were enriched in telomere/DNA replication pathways. A three-gene signature (PDE1B, TFAP2B, and HSPA1A) stratified patients into high- and low-risk groups with divergent survival (AUC >0.60, validated externally). High-risk patients had advanced N/M stages, elevated M0/M2 macrophages, reduced CD4+ memory T cells, and upregulated immune checkpoints. Nomogram integrating risk score, age, and N/M stage accurately predicted 1-/3-/5-year survival. Low-risk patients showed greater 5-fluorouracil sensitivity. PDE1B expression was significantly reduced in CRC tissues and correlated with advanced stages. Functional assays confirmed PDE1B overexpression suppressed proliferation, migration, invasion, and induced apoptosis in CRC cells. CONCLUSION: This study identifies a moderately predictive telomere maintenance-related gene signature as an independent prognostic predictor in CRC. The risk stratification model effectively discriminates patients with distinct survival patterns, tumor microenvironments, and therapeutic responses, while the integrated nomogram offers additional reference information for survival analysis, albeit with only moderate predictive accuracy. These findings indicate telomere maintenance-related gene signature could serve as a preliminary auxiliary risk stratification tool for postoperative CRC patients, PDE1B may also serve as a potential epithelial tumor-suppressor target for future preclinical studies.

PDE1B↗

PubMed2026

The Impact of Molecular Characteristics on the Efficacy of Frontline Immune Checkpoint Inhibitor Therapy in Patients with Metastatic Melanoma.

Background: Metastatic melanoma in East Asian populations is enriched for acral and mucosal subtypes and harbors a distinct molecular landscape compared with Western cutaneous melanoma. However, data on the efficacy of first-line immune checkpoint inhibitor (ICI) therapy and on the impact of molecular characteristics on treatment effect in this population remain limited. We aimed to characterize the genomic landscape of an East Asian metastatic melanoma cohort and to evaluate the predictive values of molecular markers for first-line ICI therapy. Methods: This study included 135 patients with metastatic melanoma who received first-line ICI at Samsung Medical Center between January 2022 and December 2025. Of these, 108 with paired next-generation sequencing (NGS) data were included in the molecular analysis. Survival outcomes were estimated by the Kaplan-Meier method, and the prognostic value of molecular subtype was assessed using univariable and multivariable Cox proportional hazards models. Results: Mucosal melanoma was the most common primary site (58, 43.0%), followed by acral melanoma (31, 23.0%) and cutaneous melanoma (25, 18.5%); 4 (3.0%) had uveal melanoma and 17 (12.6%) had melanoma of other or unknown primary origin. The overall objective response rate to first-line ICI was 37.8% (51 of 135), and median progression-free survival (PFS) was 6.2 months (95% confidence interval [CI] 4.0-9.6). Using the hierarchical classification, 108 patients were classified into five molecular subtypes: BRAF-altered (n = 17, 15.7%), RAS-altered (n = 17, 15.7%), KIT-altered (n = 15, 13.9%), and NF1-altered (n = 7, 6.5%), and quadruple-wild-type (n = 52, 48.1%). TMB-high status (10.2%) showed no significant association with outcomes. Molecular subtype was significantly associated with PFS (log-rank p = 0.009). After multivariable adjustment, BRAF fusion (hazard ratio [HR] 5.05, 95% CI 1.70-14.98, p = 0.003) and KIT-altered status (HR 2.91, 95% CI 1.46-5.77, p = 0.002) emerged as independent adverse prognostic factors, whereas BRAF V600 single-nucleotide variants did not differ significantly from quadruple-wild-type. Conclusions: In this East Asian metastatic melanoma cohort, BRAF fusion and KIT alterations were independent adverse prognostic factors for first-line ICI. These findings suggest that BRAF fusion and KIT-altered tumors may represent distinct subgroups that warrant further investigation.

immune checkpoint inhibitor↗

PubMed2026

Amplification-Driven S100A11 Overexpression in Hepatocellular Carcinoma Is Associated with Metabolic Reprogramming, ECM Remodelling, and Immune Evasion: A Pan-Cancer Genomic Study.

BACKGROUND: S100A11, a calcium-binding S100 family protein, is increasingly implicated in carcinogenesis, yet its molecular regulation and clinical relevance across cancers remain unclear. Hepatocellular carcinoma (HCC) carries a dismal prognosis, in part due to a lack of reliable biomarkers for risk stratification of established disease. METHODS: We conducted a pan-cancer analysis of S100A11 genomic alterations across 31 studies (10,767 samples) obtained from TCGA, encompassing copy number alterations, somatic mutations, and DNA methylation. HCC-specific analyses evaluated S100A11 expression, its potential as a diagnostic/prognostic marker, co-expression networks, and pathway enrichment using TCGA-LIHC data, with univariate and multivariate Cox regression to assess survival associations. RESULTS: S100A11 alterations were predominantly driven by copy number amplification, with the highest frequencies in hepatobiliary cancers, lung and breast cancers. Copy number amplification showed a consistent inverse relationship with promoter methylation, indicating amplification-driven transcriptional activation. In HCC, S100A11 was markedly overexpressed compared with normal liver tissue, with strong diagnostic discriminatory capacity. High S100A11 expression was significantly associated with inferior overall survival (log-rank p = 0.032; HR = 1.46, 95% CI 1.03-2.06) and remained an independent predictor of overall survival after adjustment for age, sex, and AJCC pathologic stage (HR = 1.27, 95% CI 1.01-1.60, p = 0.038). Co-expression and pathway analyses demonstrated an association between S100A11 and metabolic reprogramming, extracellular matrix remodelling, and immune dysregulation. CONCLUSIONS: These findings identify S100A11 as a candidate diagnostic and prognostic biomarker in HCC whose overexpression is associated with metabolic reprogramming, ECM remodelling, and immune dysregulation, warranting experimental validation of a mechanistic role.

ECM↗

PubMed2024

Construction and Analysis of a Mitochondrial Metabolism-Related Prognostic Model for Breast Cancer to Evaluate Survival and Immunotherapy.

As one of the most prevalent malignancies among women, breast cancer (BC) is tightly linked to metabolic dysfunction. However, the correlation between mitochondrial metabolism-related genes (MMRGs) and BC remains unclear. The training and validation datasets for BC were obtained from The Cancer Genome Atlas and Gene Expression Omnibus databases, respectively. MMRG-related data were obtained from the Molecular Signatures Database. A risk score prognostic model incorporating MMRGs was established based on univariate, LASSO, and multivariate Cox regression analyses. Independent factors affecting BC prognosis were identified through regression analysis and presented in a nomogram. Single-sample gene set enrichment analysis was employed to assess the immune levels of high-risk (HR) and low-risk (LR) groups. The sensitivity of BC patients in the two groups to common anti-tumor drugs was evaluated by utilizing the Genomics of Drug Sensitivity in Cancer database. 12 MMRGs significantly associated with survival were selected from 1234 MMRGs. A 12-gene risk score prognostic model was built. In the multivariate regression analysis incorporating classical clinical factors, the MMRG-related risk score remained an independent prognostic factor. As revealed by tumor immune microenvironment analysis, the LR group with higher survival rates had elevated immune levels. The drug sensitivity results unmasked that the LR group demonstrated higher sensitivity to Irinotecan, Nilotinib, and Oxaliplatin, while the HR group demonstrated higher sensitivity to Lapatinib. The development of MMRG characteristics provides a comprehensive understanding of mitochondrial metabolism in BC, aiding in the prediction of prognosis and tumor microenvironment, and offering promising therapeutic choices for BC patients with different MMRG risk scores.

Humans↗

PubMed2026

Development and Validation of a Prognostic Signature Based on Transcription Factors Associated with Endoplasmic Reticulum Stress in Pancreatic Adenocarcinoma.

BACKGROUND: Endoplasmic reticulum stress (ER stress) plays a crucial role in influencing the malignant behaviors of various tumors. Targeting the expression or degradation of transcription factors (TFs) offers a promising avenue for cancer treatment. However, a detailed understanding of how ER stress affects TF function and their interactions remains limited. This study aims to develop a prognostic model and identify TFs associated with ER stress in pancreatic ductal adenocarcinoma (PDAC). METHODS: We obtained gene expression profiles and corresponding clinical data from The Cancer Genome Atlas (TCGA). To develop a prognostic signature, we performed several analyses, including unsupervised clustering, enrichment analysis, immune infiltration assessment, as well as univariate, LASSO, and multivariate Cox regression analyses. Four transcription factors-STAT1, IRF6, NRF1, and RXRA-were incorporated into a risk model, which was subsequently validated using the GSE dataset. Additionally, we examined IRF6 through quantitative PCR, western blotting, flow cytometry, and immunohistochemistry in vitro using pancreatic cancer cell lines and a tissue microarray. RESULTS: The high-risk group identified by the model exhibited significant associations with immune cell infiltration and poorer survival outcomes, though there was no significant correlation with tumor purity (p = 0.19). Furthermore, IRF6 downregulation in vitro was found to inhibit pancreatic cancer cell proliferation and promote apoptosis. IRF6 depletion also increased the expression of key molecules involved in ER stress at both the transcriptional and translational levels. Immunohistochemical analysis revealed marked differences in IRF6 expression between tumor and adjacent non-tumor tissues (59.29±29.88 vs. 95.22±40.80, p<0.001). CONCLUSION: This study provides evidence that the constructed risk model can effectively predict prognosis in PDAC patients. Transcription factors related to ER stress, such as IRF6, show promise as both prognostic biomarkers and potential therapeutic targets for PDAC.

Humans↗

PubMed2025

An anti-androgen resistance-related gene signature acts as a prognostic marker and increases enzalutamide efficacy via PLK1 inhibition in prostate cancer.

BACKGROUND: Anti-androgen resistance remains a major clinical challenge in the treatment of prostate cancer (PCa), leading to disease progression and treatment failure. Despite extensive research on resistance mechanisms, a reliable prognostic model for predicting patient outcomes and guiding therapeutic strategies is still lacking. This study aimed to develop a novel gene signature related to anti-androgen resistance and evaluate its prognostic and therapeutic implications. METHODS: Anti-androgen resistance-related differentially expressed genes (ARRDEGs) were identified through transcriptomic analysis of enzalutamide- and dual enzalutamide abiraterone-resistant PCa cell lines from the GEO database. Functional enrichment analysis was performed to determine the biological roles of these genes. A prognostic gene signature was developed using univariate Cox regression, LASSO, and multivariate Cox regression models. The model was validated in independent PCa cohorts from The Cancer Genome Atlas (TCGA). Additionally, we assessed the correlation between the signature, immune infiltration, immune checkpoint expression, and drug sensitivity. The efficacy of PLK1 inhibition combined with enzalutamide was further explored using in vitro and in vivo experiments. RESULTS: We identified 304 ARRDEGs, from which three key genes (LMNB1, SSPO, and PLK1) were selected to construct a prognostic signature. This gene signature effectively stratified PCa patients into high- and low-risk groups, with the high-risk group exhibiting shorter recurrence-free survival and distinct immune characteristics. High-risk patients demonstrated elevated immune checkpoint expression (B7H3, CTLA-4, B7-1, and TIGIT), increased M2 macrophage infiltration, and enhanced sensitivity to chemotherapy and targeted therapy. Mechanistically, PLK1 inhibition potentiated the antitumor effect of enzalutamide by downregulating SLC7A11 and inducing ferroptosis, providing a potential therapeutic strategy to overcome anti-androgen resistance. CONCLUSION: We established a novel ARRDEGs-based prognostic signature that predicts PCa progression and response to chemotherapy and targeted therapy. The integration of this signature with immune profiling and drug sensitivity analysis provides a valuable tool for precision oncology in PCa. Our findings highlight the potential of PLK1 inhibition as a therapeutic strategy to enhance enzalutamide efficacy and overcome resistance.

Humans↗

PubMed2026

Posterior Segment Risk Factors for Penetrating Keratoplasty Failure.

PURPOSE: To analyze the relationship between intraoperative and postpenetrating keratoplasty (PK) posterior segment variables and PK graft survival. DESIGN: Retrospective clinical cohort study. SUBJECTS: Patients undergoing PK between May 1, 2007 and September 1, 2018 at a single tertiary center. METHODS: Chart review for PKs performed was conducted, and the first PK completed at the institution for each patient was included for analysis. Data collected included demographics, medical and ocular history, preoperative and intraoperative findings, and intraoperative and postoperative posterior segment factors (pars plana vitrectomy [PPV], endolaser, retinal detachment [RD], and vitreous hemorrhage [VH]). After univariable analysis, variables were selected for multivariable Cox regression analysis. MAIN OUTCOME MEASURE: Graft failure, defined as irreversible and visually significant corneal edema, haze, or scarring. RESULTS: Eight hundred and thirty-five eyes of 835 patients were included. Mean age was 57.1 ± 22.0 (range: 0-100) years, and mean time from PK to final follow-up or graft failure was 3.2 ± 2.9 (range: 0.01-16.1) years. Graft failure occurred in 35.0% of cases with a mean onset of 1.9 ± 2.0 (range: 0.04-11.4) years after PK. After multivariable analysis, 9 variables had significant associations with failure. Two posterior segment variables were significant: intraoperative VH at the time of PK (hazard ratio [HR] 6.6, 95% confidence interval [CI] 1.6-27.7, P = .010) and silicone oil (SO) tamponade after the PK (HR 3.2, 95% CI 1.4-7.4, P = .007). CONCLUSIONS: Graft failure is a serious complication of PK. VH at the time of the PK and SO tamponade after the PK were associated with graft failure. In complex eyes that are undergoing PK grafts and that may also require posterior segment interventions, these findings may guide patient counseling and discussion of graft prognosis.

Humans↗