mRNAsi-related metabolic risk score model identifies poor prognosis, immunoevasive contexture, and low chemotherapy response in colorectal cancer patients through machine learning

Integrated genomic analysis of the stemness index signature of mRNA expression predicts lung adenocarcinoma prognosis and immune landscape

mRNA expression-based stemness index (mRNAsi) has been used for prognostic assessment in various cancers, but its application in lung adenocarcinoma (LUAD) is limited, which is the focus of this study. Low mRNAsi in LUAD predicted a better prognosis. Eight genes (GNG7, EIF5A, ANLN, FKBP4, GAPDH, GNPNAT1, E2F7, CISH) associated with mRNAsi were screened to establish a risk model. The differentially expressed genes between the high and low risk groups were mainly enriched in the metabolism, cell cycle functions pathway. The low risk score group had higher immune cell scores. Patients with lower TIDE scores in the low risk group had better immunotherapy outcomes. In addition, risk score was effective in assessing drug sensitivity of LUAD. Reverse transcription-quantitative polymerase chain reaction (RT-qPCR) data showed that eight genes were differentially expressed in LUAD cell lines, and knockdown of EIF5A reduced the invasion and migration ability of LUAD cells. This study designed a risk model based on the eight mRNAsi-related genes for predicting LUAD prognosis. The model accurately predicted the prognosis and survival of LUAD patients, facilitating the assessment of the sensitivity of patients to immunotherapy and chemotherapy.

A novel molecular classification system based on the molecular feature score identifies patients sensitive to immune therapy and target therapy

Background

Hepatocellular carcinoma (HCC) is heterogeneous and refractory with multidimensional features. This study aims to investigate its molecular classifications based on multidimensional molecular features scores (FSs) and support classification-guided precision medicine.

Methods

Data of bulk RNA sequencing, single nucleotide variation, and single-cell RNA sequencing were collected. Feature scores (FSs) from hallmark pathways, regulatory cell death pathways, metabolism pathways, stemness index, immune scores, estimate scores, etc. were evaluated and screened. Then, the unsupervised clustering on the core FSs was performed and the characteristics of the resulting clusters were identified. Subsequently, machine learning algorithms were used to predict the classifications and prognoses. Additionally, the sensitivity to immune therapy and biological roles of classification-related prognostic genes were also evaluated.

Results

We identified four clusters with distinct characteristics. C1 is characterized by high TP53 mutations, immune suppression, and metabolic downregulation, with notable responsiveness to anti-PD1 therapy. C2 exhibited high tumor purity and metabolic activity, moderate TP53 mutations, and cold immunity. C3 represented an early phase with the most favorable prognosis, lower stemness and tumor mutations, upregulated stroma, and hypermetabolism. C4 represented a late phase with the poorest prognosis, highest stemness, higher TP53 mutations, cold immunity, and metabolic downregulation. We further developed practical software for prediction with good performance in the external validation. Additionally, FTCD was identified as a classification-specific prognostic gene with tumor-suppressing role and potential as a therapeutic target, particularly for C1 and C4 patients.

Conclusions

The four-layer classification scheme enhances the understanding of HCC heterogeneity, and we also provide robust predictive software for predicting classifications and prognoses. Notably, C1 is more sensitive to anti-PD1 therapies and FTCD is a promising therapeutic target, particularly for C1 and C4. These findings provide new insights into classification-guided precision medicine.

Construction of a pathomics model for predicting mRNAsi in lung adenocarcinoma and exploration of biological mechanism

Objective

This study aimed to predict the level of stemness index (mRNAsi) and survival prognosis of lung adenocarcinoma (LUAD) using pathomics model.

Methods

From The Cancer Genome Atlas (TCGA) database, 327 LUAD patients were randomly assigned to a training set (n = 229) and a validation set (n = 98) for pathomics model development and evaluation. PyRadiomics was used to extract pathomics features, followed by feature selection using the mRMR-RFE algorithm. In the training set, Gradient Boosting Machine (GBM) was utilized to establish a model for predicting mRNAsi in LUAD. The model's predictive performance was evaluated using ROC curves, calibration curves, and decision curve analysis (DCA). Prognostic analysis was conducted using Kaplan-Meier curves and cox regression. Additionally, gene enrichment analysis, tumor microenvironment analysis, and tumor mutational burden (TMB) analysis were performed to explore the biological mechanisms underlying the pathomics prediction model.

Results

Multivariable cox analysis (HR = 1.488, 95 % CI 1.012-2.187, P = 0.043) identified mRNAsi as a prognostic risk factor for LUAD. A total of 465 pathomics features were extracted from TCGA-LUAD histopathological images, and ultimately, the most representative 8 features were selected to construct the predictive model. ROC curves demonstrated the significant predictive value of the model for mRNAsi in both the training set (AUC = 0.769) and the validation set (AUC = 0.757). Calibration curves and Hosmer-Lemeshow goodness-of-fit test showed good consistency between the model's prediction of mRNAsi levels and the actual values. DCA indicated a good net benefit of the model. The prediction of mRNAsi levels by the pathomics model is represented using the pathomics score (PS). PS was strongly associated with the prognosis of LUAD (HR = 1.496, 95 % CI 1.008-2.222, P = 0.046). Signaling pathways related to DNA replication and damage repair were significantly enriched in the high PS group. Prediction of immune therapy response indicated significantly reduced Dysfunction in the high PS group (P < 0.001). The high PS group exhibited higher TMB values (P < 0.001).

Conclusions

The predictive model constructed based on pathomics features can forecast the mRNAsi and survival risk of LUAD. This model holds promise to aid clinical practitioners in identifying high-risk patients and devising more optimized treatment plans for patients by jointly employing therapeutic strategies targeting cancer stem cells (CSCs).

An individualized stemness-related signature to predict prognosis and immunotherapy responses for gastric cancer using single-cell and bulk tissue transcriptomes

BACKGROUND

Currently, many stemness-related signatures have been developed for gastric cancer (GC) to predict prognosis and immunotherapy outcomes. However, due to batch effects, these signatures cannot accurately analyze patients one by one, rendering them impractical in real clinical scenarios. Therefore, we aimed to develop an individualized and clinically applicable signature based on GC stemness.

METHODS

Malignant epithelial cells from single-cell RNA-Seq data of GC were used to identify stemness-related signature genes based on the CytoTRACE score. Using two bulk tissue datasets as training data, the enrichment scores of the signature genes were applied to classify samples into two subtypes. Then, using the identified subtypes as criteria, we developed an individualized stemness-related signature based on the within-sample relative expression orderings of genes.

RESULTS

We identified 175 stemness-related signature genes, which exhibited significantly higher AUCell scores in poorly differentiated GCs compared to differentiated GCs. In training datasets, GC samples were classified into two subtypes with significantly different survival times and genomic characteristics. Utilizing the two subtypes, an individualized signature was constructed containing 47 gene pairs. In four independent testing datasets, GC samples classified as high risk exhibited significantly shorter survival times, higher infiltration of M2 macrophages, and lower immune responses compared to low-risk samples. Moreover, the potential therapeutic targets and corresponding drugs were identified for the high-risk group, such as CD248 targeted by ontuxizumab.

CONCLUSIONS

We developed an individualized stemness-related signature, which can accurately predict the prognosis and efficacy of immunotherapy for each GC sample.

Integration of scRNA-seq and bulk RNA-seq constructs a stemness-related signature for predicting prognosis and immunotherapy responses in hepatocellular carcinoma

PURPOSE

Cancer stem cells are associated with unfavorable prognosis in hepatocellular carcinoma (HCC). However, existing stemness-related biomarkers and prognostic models are limited.

METHODS

The stemness-related signatures were derived from taking the union of the results obtained by performing WGCNA and CytoTRACE analysis at the bulk RNA-seq and scRNA-seq levels, respectively. Univariate Cox regression and the LASSO were applied for filtering prognosis-related signatures and selecting variables. Finally, ten gene signatures were identified to construct the prognostic model. We evaluated the differences in survival, genomic alternation, biological processes, and degree of immune cell infiltration in the high- and low-risk groups. pRRophetic and Tumor Immune Dysfunction and Exclusion (TIDE) algorithms were utilized to predict chemosensitivity and immunotherapy response. Human Protein Atlas (HPA) database was used to evaluate the protein expressions.

RESULTS

A stemness-related prognostic model was constructed with ten genes including YBX1, CYB5R3, CDC20, RAMP3, LDHA, MTHFS, PTRH2, SRPRB, GNA14, and CLEC3B. Kaplan-Meier and ROC curve analyses showed that the high-risk group had a worse prognosis and the AUC of the model in four datasets was greater than 0.64. Multivariate Cox regression analyses verified that the model was an independent prognostic indicator in predicting overall survival, and a nomogram was then built for clinical utility in predicting the prognosis of HCC. Additionally, chemotherapy drug sensitivity and immunotherapy response analyses revealed that the high-risk group exhibited a higher likelihood of benefiting from these treatments.

CONCLUSION

The novel stemness-related prognostic model is a promising biomarker for estimating overall survival in HCC.

Prediction of Prognosis and Chemotherapeutic Sensitivity Based on Cuproptosis-Associated lncRNAs in Cervical Squamous Cell Carcinoma and Endocervical Adenocarcinoma

Cervical cancer is the fourth most common cancer. The 5-year survival rate for metastatic cervical cancer is less than 10%. The survival time of patients with recurrent cervical cancer is approximately 13-17 months. Cuproptosis is a novel type of cell death related to mitochondrial respiration. Accumulative studies showed that long non-coding RNAs (lncRNAs) regulated cervical cancer progression. Compressive bioinformatic analysis showed that nine cuproptosis-related lncRNAs (CRLs), including C002128.2, AC002563.1, AC009237.14, AC048337.1, AC145423.1, AL117336.1, AP001542.3, ATP2A1-AS1, and LINC00426, were independently correlated with the overall survival (OS) of cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC) patients. The time-dependent area under curve value reached 0.716 at 1 year, 0.718 at 3 years, and 0.719 at 5 years. Notably, CESC patients in the low-risk group had increased immune cell infiltration and expression of several immune checkpoints, which indicated that they may benefit more from immune checkpoint blockade therapy. In addition, we also used the model for drug sensitivity analysis. Several drug sensitivities were more sensitive in high-risk patients and showed significant correlations with the risk models, such as Bortezomib_1191, Luminespib_1559, and Rapamycin_1084, suggesting that these drugs may be candidate clinical drugs for patients with a high risk of CESC. In summary, this study further explored the mechanism of CRLs in CESC and provided a more optimized prognostic model and some insights into chemotherapy of CESC.

Identification and Validation of Molecular Subtype and Prognostic Signature for Bladder Cancer Based on Neutrophil Extracellular Traps

Neutrophil extracellular traps (NETs) could promote tumor growth and distant metastases. Molecular subtypes of bladder cancer were identified with consensus cluster analysis. A NETs-related prognostic signature was constructed with LASSO cox regression analysis. As a result, we identified three subtypes of bladder cancer, which had a distinct difference in prognosis, immune microenvironment, TIDE score, mRNAsi score and IC50 score. We also developed a prognostic signature based on 5 NETs-related genes, which had a good performance in clinical outcome prediction of bladder cancer. These results may provide more data about the vital role of NETs in bladder cancer.