Search

High stemness index scores are associated with poor survival in patients with lung cancer. Studies on the mRNA expression-based stemness index (mRNAsi) are typically conducted using tumor tissues; however, mRNAsi-related expression signatures based on cell-free RNA (cfRNA) are yet to be comprehensively investigated. The present study aimed to elucidate the gene expression profiles of tumor stemness in lung cancer tissues and corresponding cfRNAs in blood, and to assess their links with immune infiltration. Tumor tissue, paracancerous tissue, peripheral blood and lymph node samples were collected from patients with stage I-III non-small cell lung cancer and RNA sequencing was performed. The TCGAbiolinks package was used to calculate the mRNAsi for each of these four types of sample. Weighted gene co-expression network analysis and differentially expressed gene analyses were performed to investigate mRNAsi-related genes, and pathway enrichment analysis was performed using the Kyoto Encyclopedia of Genes and Genomes (KEGG) orthology-based annotation system. In addition, the STAR-Fusion tool was used to detect fusion variants, and CIBERSORT was used to analyze the correlations of stemness signatures in tissues and blood with immune cell infiltration. The mRNAsi values in peripheral blood and lymph nodes were found to be higher than those in cancer tissues. 'Hematopoietic cell lineage' was the only KEGG pathway enriched in mRNAsi-related genes in both lung cancer tissues and peripheral blood. In addition, the protein tyrosine phosphatase receptor type C associated protein gene was the only gene commonly associated with the mRNAsi in these two types of sample. The expression of mRNAsi-related genes was increased in the dendritic and Treg cells in tumor tissues, but was elevated in Treg and CD8 cells in the blood. In conclusion, cfRNAs in the blood exhibit unique stemness signatures that have potential for use in the diagnosis of lung cancer. [ABSTRACT FROM AUTHOR]

BackgroundCancer stem cells are characterized by self-renewal, clonal tumor initiation capacity, and treatment resistance, which play essential roles in the tumor progression of prostate cancer (PCa). In this study, we aim to explore the features of cancer stemness and characterize the expression of stem cell-related genes for PCa.MethodsWe downloaded RNA-seq data and related clinical information from The Cancer Genome Atlas (TCGA) database. The mRNA stemness index (mRNAsi) was analyzed for various clinical features, overall survival (OS), and disease-free survival (DFS), and a weighted gene co-expression network analysis (WGCNA) was performed to identify crucial gene modules and key genes, which may play a role in CSCs. The key gene functions were verified using multiple databases, including the TCGA and Gene Expression Omnibus database (GEO). Next, we explored the potential function of the modules and genes obtained using WGCNA using an enrichment analysis. Finally, we performed in vitro experiments for further verification.ResultsWe found that mRNAsi were higher in PCa tissues than in normal tissues, and the mRNAsi were closely related to the clinical characteristics of PCa. A total of 16 key genes associated with the mRNAsi scores were identified by WGCNA analysis, including NCAPG, NEK2, DLGAP5, CENPA, CENPF, TPX2, GTSE1, KIF4A, NEIL3, CDC25C, UBE2C, CDCA5, MELK, SKA3, NUF2, and BIRC5. These genes were explicitly highly expressed in PCa across TCGA cancers and were validated in 3 independent GEO PCa datasets. The functional annotations of the key genes were linked with the cell proliferation processes. NUF2 may be a potential biomarker for PCa. In vitro experiments showed that knockdown NUF2 reduced the proliferation and migration of PCa cells.ConclusionThe 16 key genes identified in this study significantly correlate with PCa stem cell characteristics and showed prognosis-oriented effects in PCa patients. Further, the NUF2 gene may be used as a drug target for treating PCa.

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

Abstract Lung adenocarcinoma (LUAD) is one of the sole causes of death in lung cancer patients. This study combined with single-cell RNA-seq analysis to identify tumor stem-related prognostic models to predict the prognosis of lung adenocarcinoma, chemotherapy agents, and immunotherapy efficacy. mRNA expression-based stemness index (mRNAsi) was determined by One Class Linear Regression (OCLR). Differentially expressed genes (DEGs) were detected by limma package. Single-cell RNA-seq analysis in GSE123902 dataset was performed using Seurat package. Weighted Co-Expression Network Analysis (WGCNA) was built by rms package. Cell differentiation ability was determined by CytoTRACE. Cell communication analysis was performed by CellCall and CellChat package. Prognosis model was constructed by 10 machine learning and 101 combinations. Drug predictive analysis was conducted by pRRophetic package. Immune microenvironment landscape was determined by ESTIMATE, MCP-Counter, ssGSEA analysis. Tumor samples have higher mRNAsi, and the high mRNAsi group presents a worse prognosis. Turquoise module was highly correlated with mRNAsi in TCGA-LUAD dataset. scRNA analysis showed that 22 epithelial cell clusters were obtained, and higher CSCs malignant epithelial cells have more complex cellular communication with other cells and presented dedifferentiation phenomenon. Cellular senescence and Hippo signaling pathway are the major difference pathways between high- and low CSCs malignant epithelial cells. The pseudo-temporal analysis shows that cluster1, 2, high CSC epithelial cells, are concentrated at the end of the differentiation trajectory. Finally, 13 genes were obtained by intersecting genes in turquoise module, Top200 genes in hdWGCNA, DEGs in high- and low- mRNAsi group as well as DEGs in tumor samples vs. normal group. Among 101 prognostic models, average c-index (0.71) was highest in CoxBoost + RSF model. The high-risk group samples had immunosuppressive status, higher tumor malignancy and low benefit from immunotherapy. This work found that malignant tumors and malignant epithelial cells have high CSC characteristics, and identified a model that could predict the prognosis, immune microenvironment, and immunotherapy of LUAD, based on CSC-related genes. These results provided reference value for the clinical diagnosis and treatment of LUAD.

Objective: An mRNA expression-based stemness index (mRNAsi) has been developed to characterize cancer stemness. However, the predictive value of mRNAsi-based signature in therapeutic resistance and immunotherapy in thyroid cancer (THCA) remains unclarified. This study evaluated and validated the role of mRNAsi in drug sensitivity, its relationship between mRNAsi and THCA clinical features and immunity based on bioinformatics. Methods: Based on transcriptome data of THCA patients from the Tumor Genome Atlas Project (TCGA) database, and expression data of multifunctional stem cell samples from the Progenitor Cell Biology Consortium (PCBC) databases, mRNAsi was calculated by the '' one class logistic regression (OCLR)'' method, Molecular subtypes of TCGA-THCA samples were identified with mRNAsi-related genes using ConsensusClusterPlus method. The gene mutation, clinical characteristics, immune characteristics, TIDE and drug sensitivity were compared among molecular subtypes. A prognostic model was designed with Lasso cox method. Modulation of malignant phenotype of THCA cell lines by model characterization genes is validated by CCK-8, flow cytometry. DNA methylation disorder in promoter region was analyzed between risk groups. The model was validated for survival in the internal Test dataset, while TCGA pan-cancer and immunotherapy datasets were further employed to validate the performance of this model. Results: We obtained a total of 78 stem cell samples, each containing the expression profile of 8087 mRNA genes. Based on mRNAsi, THCA was divided into 3 subtypes. Subtype C2 had the poorest prognosis and highest immune score, while subtype C3 had the best prognosis, lowest mRANsi and highest TIDE score. Patients in subtype C2 showed higher sensitivity to Cisplatin, Erlotinib, Paclitaxel, and Lapatinib. The prognostic signature was generated using 5 mRNAsi-related genes, which could predict prognosis for THCA. qRT-PCR results showed that the expression of 5 genes were various in Hth7 and KTC-1 cells, and inhibition CELSR3 expression increased percentage of apoptosis in Hth7 and KTC-1 cells. mRNAsi related DNA methylation sites were mainly enriched in tumor related pathways. Good performance of this model was validated in Test dataset, pan-cancer and immunotherapy datasets. Conclusion: This study identified three subtypes for classification and developed a prognostic model with mRNAsi-related genes, which provided great potential for prognosis and immunotherapy prediction.

Background: Lung adenocarcinoma (LUAD) is the most prevalent subtype of lung cancer, and comprehending its molecular mechanisms is pivotal for advancing treatment efficacy. This study aims to explore the prognostic and functional significance of base excision repair (BER)-related long non-coding RNAs (BERLncs) in LUAD. Methods: A risk score model for BERLncs was developed using the least absolute shrinkage and selection operator regression and Cox regression analysis. Model validation and prognostic evaluation were performed using Kaplan–Meier and receiver-operating characteristic curve analyses. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses were conducted to elucidate the potential biological functions of BERLncs. Comparative analyses were carried out to investigate disparities in tumor mutation burden (TMB), immune infiltration, tumor immune dysfunction and exclusion (TIDE) score, chemosensitivity, and immune checkpoint gene expression between the two risk groups. Results: A predictive risk score model comprising 19 BERLncs was successfully developed. Patients were divided into high-risk and low-risk groups according to the median risk score. The high-risk subgroup exhibited significantly inferior overall survival. Functional enrichment analysis revealed pathways associated with lung cancer development, notably the neuroactive ligand–receptor interaction pathway. High-risk patients demonstrated elevated TMB, diminished TIDE scores, and an immunosuppressive tumor microenvironment, while low-risk patients displayed potential benefits from immunotherapy. Additionally, the risk model identified potential anticancer agents. Conclusion: The risk score model based on BERLncs shows promise as a prognostic biomarker for LUAD patients, providing valuable insights for clinical decision-making, therapeutic strategies, and understanding of underlying biological mechanisms. [ABSTRACT FROM AUTHOR]

Due to the high heterogeneity, lung adenocarcinoma (LUAD) cannot be distinguished into precise molecular subtypes, thereby resulting in poor therapeutic effect and low 5-year survival rate clinically. Although the tumor stemness score (mRNAsi) has been shown to accurately characterize the similarity index of cancer stem cells (CSCs), whether mRNAsi can serve as an effective molecular typing tool for LUAD isn't reported to date. In this study, we first demonstrate that mRNAsi is significantly correlated with the prognosis and disease degree of LUAD patients, i.e., the higher the mRNAsi, the worse the prognosis and the higher the disease degree. Second, we identify 449 mRNAsi-related genes based on both weighted gene co-expression network analysis (WGCNA) and univariate regression analysis. Third, our results display that 449 mRNAsi-related genes can accurately distinguish the LUAD patients into two molecular subtypes: ms-H subtype (with high mRNAsi) and ms-L subtype (with low mRNAsi), particularly the ms-H subtype has a worse prognosis. Remarkably, significant differences in clinical characteristics, immune microenvironment, and somatic mutation exist between the two molecular subtypes, which might lead to the poorer prognosis of the ms-H subtype patients than that of the ms-L subtype ones. Finally, we establish a prognostic model containing 8 mRNAsi-related genes, which can effectively predict the survival rate of LUAD patients. Taken together, our work provides the first molecular subtype related to mRNAsi in LUAD, and reveals that these two molecular subtypes, the prognostic model and marker genes may have important clinical value for effectively monitoring and treating LUAD patients. [ABSTRACT FROM AUTHOR]

Abstract Background N6-methyladenosine (m6A) modification is involved in tumorigenesis and progression as well as closely correlated with stem cell differentiation and pluripotency. Moreover, tumor progression includes the acquisition of stemness characteristics and accumulating loss of differentiation phenotype. Therefore, we integrated m6A modification and stemness indicator mRNAsi to classify patients and predict prognosis for LGG. Methods We performed consensus clustering, weighted gene co-expression network analysis, and least absolute shrinkage and selection operator Cox regression analysis to identify an m6A regulation- and mRNAsi-related prognostic index (MRMRPI). Based on this prognostic index, we also explored the differences in immune microenvironments between high- and low-risk populations. Next, immunotherapy responses were also predicted. Moreover, single-cell RNA sequencing data was further used to verify the expression of these genes in MRMRPI. At last, the tumor-promoting and tumor-associated macrophage polarization roles of TIMP1 in LGG were validated by in vitro experiments. Results Ten genes (DGCR10, CYP2E1, CSMD3, HOXB3, CABP4, AVIL, PTCRA, TIMP1, CLEC18A, and SAMD9) were identified to construct the MRMRPI, which was able to successfully classify patients into high- and low-risk group. Significant differences in prognosis, immune microenvironment, and immunotherapy responses were found between distinct groups. A nomogram integrating the MRMRPI and other prognostic factors were also developed to accurately predict prognosis. Moreover, in vitro experiments illustrated that inhibition of TIMP1 could inhibit the proliferation, migration, and invasion of LGG cells and also inhibit the polarization of tumor-associated macrophages. Conclusion These findings provide novel insights into understanding the interactions of m6A methylation regulation and tumor stemness on LGG development and contribute to guiding more precise immunotherapy strategies.

Background: Increasing evidence suggests that hepatocellular carcinoma (HCC) stem cells (LCSCs) play an essential part in HCC recurrence, metastasis, and chemotherapy and radiotherapy resistance. Multiple studies have demonstrated that stemness-related genes facilitate the progression of tumors. However, the mechanism by which stemness-related genes contribute to HCC is not well understood. Here, we aim to construct a stemness-related score (SRscores) model for deeper analysis of stemness-related genes, assisting with the prognosis and individualized treatment of HCC patients.Further, we found that the gene LPCAT1 was highly expressed in tumor tissues by immunohistochemistry, and sphere-forming assay revealed that knockdown of LPCAT1 inhibited the sphere-forming ability of hepatocellular carcinoma cells. Methods: We used the TCGA-LIHC dataset to screen stemness-related genes of HCC from the MSigDB database. Prognosis, tumor microenvironment, immunological checkpoints, tumor immune dysfunction, rejection, treatment sensitivity, and putative biological pathways were examined. Random forest created the SRscores model. The anti-PD-1/anti-CTLA4 immunotherapy, tumor mutational burden, medication sensitivity, and cancer stem cell index were compared between the high- and low-risk score groups. We also examined risk scores for different cell types using single-cell RNA sequencing data and correlated transcription factor activity in cancer stem cells with SRscores genes. Finally, we tested core marker expression and biological functions. Results: Patients can be divided into two subtypes (Cluster1 and Cluster2) based on the TCGA-LIHC dataset's identification of 11 stemness-related genes. Additionally, a SRscores was developed based on subtypes. Cluster2 and the group with the lowest SRscores had superior survival and immunotherapy response than Cluster1 and the group with the highest SRscores. The group with a high SRscores was significantly more enriched in classical tumor pathways than the group with a low SRscores. Multiple transcription factors and SRscores genes are correlated. The core gene LPCAT1 is highly expressed in rat liver cancer tissues and promotes tumor cell sphere formation. Conclusion: A SRscores model can be utilized to predict the prognosis of HCC patients as well as their response to immunotherapy. [ABSTRACT FROM AUTHOR]

Hepatocellular carcinoma (HCC) is a malignant lethal tumor and both cancer stem cells (CSCs) and metabolism reprogramming have been proven to play indispensable roles in HCC. This study aimed to reveal the connection between metabolism reprogramming and the stemness characteristics of HCC, established a new gene signature related to stemness and metabolism and utilized it to assess HCC prognosis and immunotherapy response. The clinical information and gene expression profiles (GEPs) of 478 HCC patients came from the Gene Expression Omnibus (GEO) and the Cancer Genome Atlas (TCGA). The one-class logistic regression (OCLR) algorithm was employed to calculate the messenger ribonucleic acid expression-based stemness index (mRNAsi), a new stemness index quantifying stemness features. Differentially expressed analyses were done between high- and low-mRNAsi groups and 74 differentially expressed metabolism-related genes (DEMRGs) were identified with the help of metabolism-related gene sets from Molecular Signatures Database (MSigDB). After integrated analysis, a risk score model based on the three most efficient prognostic DEMRGs, including Recombinant Phosphofructokinase Platelet (PFKP), phosphodiesterase 2A (PDE2A) and UDP-glucuronosyltransferase 1A5 (UGT1A5) was constructed and HCC patients were divided into high-risk and low-risk groups. Significant differences were found in pathway enrichment, immune cell infiltration patterns, and gene alterations between the two groups. High-risk group patients tended to have worse clinical outcomes and were more likely to respond to immunotherapy. A stemness-metabolism-related model composed of gender, age, the risk score model and tumor-node- metastasis (TNM) staging was generated and showed great discrimination and strong ability in predicting HCC prognosis and immunotherapy response. [ABSTRACT FROM AUTHOR]

Background: Stem cells play an important role in acute myeloid leukemia (AML). However, their precise effect onAML tumorigenesis and progression remains unclear. Methods: The present study aimed to characterize stem cell-related gene expression and identify stemness biomarker genes in AML. We calculated the stemness index (mRNAsi) based on transcription data using the one-class logistic regression (OCLR) algorithm for patients in the training set. According to the mRNAsi score, we performed consensus clustering and identified two stemness subgroups. Eight stemness-related genes were identified as stemness biomarkers through gene selection by three machine learning methods. Results: We found that patients in stemness subgroup I had a poor prognosis and benefited from nilotinib, MK-2206 and axitinib treatment. In addition, the mutation profiles of these two stemness subgroups were different, which suggested that patients in different subgroups had different biological processes. There was a strong significant negative correlation between mRNAsi and the immune score (r= -0.43, p<0.001). Furthermore, we identified eight stemness-related genes that have potential to be biomarkers, including SLC43A2, CYBB, CFP, GRN, CST3, TIMP1, CFD and IGLL1. These genes, except IGLL1, had a negative correlation with mRNAsi. SLC43A2 is expected to be a potential stemness-related biomarker in AML. Conclusion: Overall, we established a novel stemness classification using the mRNAsi score and eight stemness-related genes that may be biomarkers. Clinical decision-making should be guided by this new signature in prospective studies. [ABSTRACT FROM AUTHOR]

Highlights IL-6, FPR2, RLN3 had higher degree in PPI network. The patients in cluster 1 had better prognosis. Total 10 immune checkpoint genes showed higher expression in cluster 1. TMB value in cluster 2 was higher than that in cluster 1.

Simple Summary: In this study; four immune-related predictive biomarkers for LGG were identified and proven to be IRGs. The development of more efficient immunotherapy techniques was facilitated by the creation of a prognostic signature to evaluate and forecast the prognosis of LGG patients. The immune-related prognostic score was then used to build a nomogram that was subsequently used to more accurately predict the prognosis of LGG patients Background: Low-grade gliomas (LGGs), which are the second most common intracranial tumor, are diagnosed in seven out of one million people, tending to develop in younger people. Tumor stem cells and immune cells are important in the development of tumorigenesis. However, research on prognostic factors linked to the immune microenvironment and stem cells in LGG patients is limited. We critically need accurate related tools for assessing the risk of LGG patients. Methods: In this study, we aimed to identify immune-related genes (IRGs) in LGG based on the mRNAsi score. We employed differentially expressed gene (DEG) methods and weighted correlation network analysis (WGCNA). The risk signature was then further established using a lasso Cox regression analysis and a multivariate Cox analysis. Next, we used immunohistochemical sections (HPA) and a survival analysis to identify the hub genes. A nomogram was built to assess the prognosis of patients based on their clinical information and risk scores and was validated using a DCA curve, among other methods. Results: Four hub genes were obtained: C3AR1 (HR = 0.98, p < 0.001), MSR1 (HR = 1.02, p < 0.001), SLC11A1 (HR = 1.01, p < 0.01), and IL-10 (HR = 1.01, p < 0.001). For LGG patients, we created an immune-related prognostic signature (IPS) based on mRNAsi for estimating risk scores; different risk groups showed significantly different survival rates (p = 3.3 × 10−16). Then, via an evaluation of the IRG-related signature, we created a nomogram for predicting LGG survival probability. Conclusion: The outcome suggests that, when predicting the prognosis of LGG patients, our nomogram was more effective than the IPS. In this study, four immune-related predictive biomarkers for LGG were identified and proven to be IRGs. Therefore, the development of efficient immunotherapy techniques can be facilitated by the creation of the IPS. [ABSTRACT FROM AUTHOR]

Background: N6-methyladenosine (m6A) modification is involved in tumorigenesis and progression as well as closely correlated with stem cell differentiation and pluripotency. Moreover, tumor progression includes the acquisition of stemness characteristics and accumulating loss of differentiation phenotype. Therefore, we integrated m6A modification and stemness indicator mRNAsi to classify patients and predict prognosis for LGG. Methods: We performed consensus clustering, weighted gene co-expression network analysis, and least absolute shrinkage and selection operator Cox regression analysis to identify an m6A regulation- and mRNAsi-related prognostic index (MRMRPI). Based on this prognostic index, we also explored the differences in immune microenvironments between high- and low-risk populations. Next, immunotherapy responses were also predicted. Moreover, single-cell RNA sequencing data was further used to verify the expression of these genes in MRMRPI. At last, the tumor-promoting and tumor-associated macrophage polarization roles of TIMP1 in LGG were validated by in vitro experiments. Results: Ten genes (DGCR10, CYP2E1, CSMD3, HOXB3, CABP4, AVIL, PTCRA, TIMP1, CLEC18A, and SAMD9) were identified to construct the MRMRPI, which was able to successfully classify patients into high- and low-risk group. Significant differences in prognosis, immune microenvironment, and immunotherapy responses were found between distinct groups. A nomogram integrating the MRMRPI and other prognostic factors were also developed to accurately predict prognosis. Moreover, in vitro experiments illustrated that inhibition of TIMP1 could inhibit the proliferation, migration, and invasion of LGG cells and also inhibit the polarization of tumor-associated macrophages. Conclusion: These findings provide novel insights into understanding the interactions of m6A methylation regulation and tumor stemness on LGG development and contribute to guiding more precise immunotherapy strategies. [ABSTRACT FROM AUTHOR]

In our study, 49 key genes significantly associated with renal cell carcinoma (RCC) stemness were obtained. Next, we developed a molecular prognostic signature associated with stemness features of pan-RCC. The difference in overall survival (OS) between the high- and low-risk groups was statistically significant (p <.05). The area under the receiver operating characteristic curve for 1-year OS, 5-year OS, and 10-year OS was 0.759, 0.712, and 0.918, respectively. The results of validation in The Cancer Genome Atlas cohort and International Cancer Genome Consortium cohort revealed the predictive capability of this signature. Furthermore, we selected three genes and further validation showed that these three hub genes were potential hub biomarkers for pan-RCC stemness features. [ABSTRACT FROM AUTHOR]

Metastasis of osteosarcoma is an important adverse factor affecting patients' survival, and cancer stemness is the crucial cause of distant metastasis. Capsaicin, the main component of pepper, has been proven in our previous work to inhibit osteosarcoma proliferation and enhance its drug sensitivity to cisplatin at low concentrations. This study aims to further explore the anti-osteosarcoma effect of capsaicin at low concentrations (100 μ M, 24 h) on stemness and metastasis. The stemness of human osteosarcoma (HOS) cells was decreased significantly by capsaicin treatment. Additionally, the capsaicin treatment's inhibition of cancer stem cells (CSCs) was dose-dependent on both sphere formation and sphere size. Meanwhile, capsaicin inhibited invasion and migration, which might be associated with 25 metastasis-related genes. SOX2 and EZH2 were the most two relevant stemness factors for capsaicin's dose-dependent inhibition of osteosarcoma. The mRNAsi score of HOS stemness inhibited by capsaicin was strongly correlated with most metastasis-related genes of osteosarcoma. Capsaicin downregulated six metastasis-promoting genes and up-regulated three metastasis-inhibiting genes, which significantly affected the overall survival and/or disease-free survival of patients. In addition, the CSC re-adhesion scratch assay demonstrated that capsaicin inhibited the migration ability of osteosarcoma by inhibiting its stemness. Overall, capsaicin exerts a significant inhibitory effect on the stemness expression and metastatic ability of osteosarcoma. Moreover, it can inhibit the migratory ability of osteosarcoma by suppressing its stemness via downregulating SOX2 and EZH2. Therefore, capsaicin is expected to be a potential drug against osteosarcoma metastasis due to its ability to inhibit cancer stemness. [ABSTRACT FROM AUTHOR]

Objective: Cancer stemness and M2 macrophages are intimately linked to the prognosis of lung adenocarcinoma (LUAD). For this reason, this investigation sought to identify the key genes relevant to cancer stemness and M2 macrophages, explore the relationship between these genes and clinical characteristics, and determine the potential mechanism. Methods: LUAD transcriptomic data was analyzed from The Cancer Genome Atlas (TCGA) as well as the Gene Expression Omnibus databases. Differential expression analysis was performed to discern abnormally expressed genes between LUAD and control samples in TCGA cohort. The Cell type Identification by Estimating Relative Subsets of RNA Transcripts (CIBERSORT) algorithm was applied to determine varyingly infiltrated immune cells in LUAD compared with the control samples in TCGA cohort. Weighted correlation network analysis (WGCNA) was performed to identify genes associated with mRNA expression-based stemness index (mRNAsi) and M2 macrophages. Least absolute shrinkage and selection operator (LASSO), RandomForest (RF) and support vector machine-recursive feature elimination (SVM-RFE) machine learning methods were conducted to detect gene signatures. Global survival evaluation (Kaplan-Meier curve) was applied to investigate the relationship between gene signatures and the survival time of LUAD patients. Receiver operating characteristic (ROC) curves were produced to define biomarkers relevant to diagnosis. Gene Set Enrichment Analysis (GSEA) was performed to probe the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways related to diagnostic biomarkers. The public single-cell dataset of LUAD (GSE123902) was used to investigate the expression differences of diagnostic biomarkers in various cell types in LUAD. Real-time quantitative PCR (qRT-PCR) was performed to confirm key genes in lung adenocarcinoma cells. Results: A total of 5,410 differentialy expressed genes (DEGs) as well as 15 differentially infiltrated immune cells were identified between LUAD and control sepcimens in TCGA cohort. Thirty-seven DEGs were associated with both M2 macrophages and mRNAsi according to the WGCNA analysis. Sixteen common gene signatures were obtained using three diverse algorithms. CBFA2T3, DENND3 and FCAMR were correlated to overall and disease-free survival of LUAD patients. ROC curves revealed that CBFA2T3 and DENND3 expression accurately classified LUAD and control samples. The results of single cell related analysis showed that two diagnostic biomarkers expressions were differed between the different tissue sources in M2-like macrophages. QRT-PCR demonstrated the mRNA expressions of CBFA2T3 and DENND3 were upregulated in lung adenocarcinoma cells A549 and H2122. Conclusion: Our study identified CBFA2T3 and DENND3 as key genes associated with mRNAsi and M2 macrophages in LUAD and investigated the potential molecular mechanisms underlying this relationship.

Background: Recent studies have shown that the mRNA expression-based stemness index (mRNAsi) can accurately quantify the similarity of cancer cells to stem cells, and mRNAsi-related genes are used as biomarkers for cancer. However, mRNAsi-driven tumor heterogeneity is rarely investigated, especially whether mRNAsi can distinguish hepatocellular carcinoma (HCC) into different molecular subtypes is still largely unknown. Methods: Using OCLR machine learning algorithm, weighted gene co-expression network analysis, consistent unsupervised clustering, survival analysis and multivariate cox regression etc. to identify biomarkers and molecular subtypes related to tumor stemness in HCC. Results: We firstly demonstrate that the high mRNAsi is significantly associated with the poor survival and high disease grades in HCC. Secondly, we identify 212 mRNAsi-related genes that can divide HCC into three molecular subtypes: low cancer stemness cell phenotype (CSCP-L), moderate cancer stemness cell phenotype (CSCP-M) and high cancer stemness cell phenotype (CSCP-H), especially over-activated ribosomes, spliceosomes and nucleotide metabolism lead to the worst prognosis for the CSCP-H subtype patients, while activated amino acids, fatty acids and complement systems result in the best prognosis for the CSCP-L subtype. Thirdly, we find that three CSCP subtypes have different mutation characteristics, immune microenvironment and immune checkpoint expression, which may cause the differential prognosis for three subtypes. Finally, we identify 10 robust mRNAsi-related biomarkers that can effectively predict the survival of HCC patients. Conclusions: These novel cancer stemness-related CSCP subtypes and biomarkers in this study will be of great clinical significance for the diagnosis, prognosis and targeted therapy of HCC patients.

Background: Tumor stem cells (TSCs) have been widely reported to play a critical role in tumor progression and metastasis. We explored the role of tumor stemness in intrahepatic cholangiocarcinoma (iCCA) and established a prognostic risk model related to tumor stemness for prognosis prediction and clinical treatment guidance in iCCA patients.Materials and Methods: The expression profiles of iCCA samples (E-MTAB-6389 and GSE107943 cohorts) were used in the study. One-class logistic regression algorithm calculated the mRNA stemness index (mRNAsi). The mRNAsi-related genes were used as a basis for the identification of mRNAsi-related molecular subtypes through consensus clustering. The immune characteristics and biological pathways of different subtypes were assessed. The mRNAsi-related risk model was constructed with differentially expressed genes (DEGs) between subtypes.Results: The patients with high mRNAsi had longer overall survival than that with low mRNAsi. Two subtypes were identified with that C2 had higher mRNAsi and better prognosis than C1. Tumor-related pathways such as TGF-β and epithelial-mesenchymal transition (EMT) were activated in C1. C1 had higher enrichment of cancer-associated fibroblasts and tumor-associated macrophages, as well as higher immune response and angiogenesis score than C2. We screened a total 98 prognostic DEGs between C1 and C2. Based on the prognostic DEGs, we constructed a risk model containing three genes (ANO1, CD109, and CTNND2) that could divide iCCA samples into high- and low-risk groups. The two groups had distinct prognosis and immune characteristics. Notably, the risk score was negatively associated with mRNAsi (R = −0.53). High-risk group had higher enrichment score of T cell inflamed GEP, INF-γ, and cytolytic activity, and lower score of estimated IC50 of 5-fluorouracil and cisplatin than low-risk group.Conclusions: This study clarified the important role of tumor stemness in iCCA and developed an mRNAsi-related risk model for predicting the prognosis and supporting the clinical treatment in iCCA patients. The three genes (ANO1, CD109, and CTNND2) may serve as potential targets for iCCA treatment.

Background: To comprehensively analyze the stemness characteristics related to prognosis and the immune microenvironment in lung adenocarcinoma (LUAD).Methods: The OCLR machine learning method was used to calculate the stemness index (mRNAsi) of the LUAD samples. DEGs common between the low mRNAsi, normal, and high mRNAsi groups were screened and the immune-stemness genes were obtained. Then the PPI network was created and enrichment analyses were performed. Moreover, different subtypes based on immune-stemness genes associated with prognosis were identified, and the relationships between LUAD stemness and TIME variables were systematically analyzed, followed by TMB analysis.Results: Patients in the high mRNAsi groups with poor prognosis were screened along with 144 immune-stemness genes. IL-6, FPR2, and RLN3 showed a higher degree in the PPI network. A total of 26 immune-stemness genes associated with prognosis were screened. Two clusters were obtained (cluster 1 and cluster 2). Survival analysis revealed that patients in cluster 2 had a poor prognosis. A total of 12 immune cell subpopulations exhibited significant differences between cluster 1 and cluster 2 (P < 0.05). A total of 10 immune checkpoint genes exhibited significantly higher expression in cluster 1 (P < 0.05) than in cluster 2. Further, the TMB value in cluster 2 was higher than that in cluster 1 (P < 0.05).Conclusion: Immune-stemness genes, including L-6, FPR2, and RLN3, might play significant roles in LUAD development via cytokine-cytokine receptor interaction, neuroactive ligand‒receptor interaction, and the JAK‒STAT pathway. Immune-stemness genes were related to tumor-infiltrating immune cells, TMB, and expression of immune checkpoint gene. [ABSTRACT FROM AUTHOR]

Colorectal cancer (CRC) is one of the most fatal cancers of the digestive system. Although cancer stem cells and metabolic reprogramming have an important effect on tumor progression and drug resistance, their combined effect on CRC prognosis remains unclear. Therefore, we generated a 21-gene mRNA stemness index-related metabolic risk score model, which was examined in The Cancer Genome Atlas and Gene Expression Omnibus databases (1323 patients) and validated using the Zhongshan Hospital cohort (200 patients). The high-risk group showed more immune infiltrations; higher levels of immunosuppressive checkpoints, such as CD274, tumor mutation burden, and resistance to chemotherapeutics; potentially better response to immune therapy; worse prognosis; and advanced stage of tumor node metastasis than the low-risk group. The combination of risk score and clinical characteristics was effective in predicting overall survival. Zhongshan cohort validated that high-risk score group correlated with malignant progression, worse prognosis, inferior adjuvant chemotherapy responsiveness of CRC, and shaped an immunoevasive contexture. This tool may provide a more accurate risk stratification in CRC and screening of patients with CRC responsive to immunotherapy. [ABSTRACT FROM AUTHOR]

Background: Pancreatic adenocarcinoma (PAAD) is one of the most aggressive and fatal gastrointestinal malignancies with high morbidity and mortality worldwide. Accumulating evidence has revealed the clinical significance of the interaction between the hypoxic microenvironment and cancer stemness in pancreatic cancer progression and therapies. This study aims to identify a hypoxia-stemness index-related gene signature for risk stratification and prognosis prediction in PAAD. Methods: The mRNA expression-based stemness index (mRNAsi) data of PAAD samples from The Cancer Genome Atlas (TCGA) database were calculated based on the one-class logistic regression (OCLR)machine learning algorithm. Univariate Cox regression and LASSO regression analyses were then performed to establish a hypoxia-mRNAsi-related gene signature, and its prognostic performance was verified in both the TCGA-PAAD and GSE62452 corhorts by Kaplan-Meier and receiver operating characteristic (ROC) analyses. Additionally, we further validated the expression levels of signature genes using the TCGA, GTEx and HPA databases as well as qPCR experiments. Moreover, we constructed a prognostic nomogram incorporating the eight-gene signature and traditional clinical factors and analyzed the correlations of the risk score with immune infiltrates and immune checkpoint genes. Results: The mRNAsi values of PAAD samples were significantly higher than those of normal samples (p < 0.001), and PAAD patients with high mRNAsi values exhibited worse overall survival (OS). A novel prognostic risk model was successfully constructed based on the eight-gene signature comprising JMJD6, NDST1, ENO3, LDHA, TES, ANKZF1, CITED, and SIAH2, which could accurately predict the 1-, 3-, and 5-year OS of PAAD patients in both the training and external validation datasets. Additionally, the eight-gene signature could distinguish PAAD samples from normal samples and stratify PAAD patients into low- and high-risk groups with distinct OS. The risk score was closely correlated with immune cell infiltration patterns and immune checkpoint molecules. Moreover, calibration analysis showed the excellent predictive ability of the nomogram incorporating the eight-gene signature and traditional clinical factors. Conclusion: We developed a hypoxia-stemness-related prognostic signature that reliably predicts the OS of PAAD. Our findings may aid in the risk stratification and individual treatment of PAAD patients. [ABSTRACT FROM AUTHOR]

Background: Lung adenocarcinoma (LUAD) is one of the most lethal malignancies and is currently lacking in effective biomarkers to assist in diagnosis and therapy. The aim of this study is to investigate hub genes and develop a risk signature for predicting prognosis of LUAD patients. Methods: RNA-sequencing data and relevant clinical data were downloaded from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database. Weighted gene co-expression network analysis (WGCNA) was performed to identify hub genes associated with mRNA expression-based stemness indices (mRNAsi) in TCGA. We utilized LASSO Cox regression to assemble our predictive model. To validate our predictive model, me applied it to an external cohort. Results: mRNAsi index was significantly associated with the tissue type of LUAD, and high mRNAsi scores may have a protective influence on survival outcomes seen in LUAD patients. WGCNA indicated that the turquoise module was significantly correlated with the mRNAsi. We identified a 9-gene signature (CENPW, MCM2, STIL, RACGAP1, ASPM, KIF14, ANLN, CDCA8, and PLK1) from the turquoise module that could effectively identify a high-risk subset of these patients. Using the Kaplan-Meier survival curve, as well as the time-dependent receiver operating characteristic (tdROC) analysis, we determined that this gene signature had a strong predictive ability (AUC = 0.716). By combining the 9-gene signature with clinicopathological features, we were able to design a predictive nomogram. Finally, we additionally validated the 9-gene signature using two external cohorts from GEO and the model proved to be of high value. Conclusion: Our study shows that the 9-gene mRNAsi-related signature can predict the prognosis of LUAD patient and contribute to decisions in the treatment and prevention of LUAD patients. [ABSTRACT FROM AUTHOR]

Background: Lung adenocarcinoma (LUAD) is one of the most lethal malignancies and is currently lacking in effective biomarkers to assist in diagnosis and therapy. The aim of this study is to investigate hub genes and develop a risk signature for predicting prognosis of LUAD patients. Methods: RNA-sequencing data and relevant clinical data were downloaded from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database. Weighted gene co-expression network analysis (WGCNA) was performed to identify hub genes associated with mRNA expression-based stemness indices (mRNAsi) in TCGA. We utilized LASSO Cox regression to assemble our predictive model. To validate our predictive model, me applied it to an external cohort. Results: mRNAsi index was significantly associated with the tissue type of LUAD, and high mRNAsi scores may have a protective influence on survival outcomes seen in LUAD patients. WGCNA indicated that the turquoise module was significantly correlated with the mRNAsi. We identified a 9-gene signature (CENPW, MCM2, STIL, RACGAP1, ASPM, KIF14, ANLN, CDCA8, and PLK1) from the turquoise module that could effectively identify a high-risk subset of these patients. Using the Kaplan-Meier survival curve, as well as the time-dependent receiver operating characteristic (tdROC) analysis, we determined that this gene signature had a strong predictive ability (AUC = 0.716). By combining the 9-gene signature with clinicopathological features, we were able to design a predictive nomogram. Finally, we additionally validated the 9-gene signature using two external cohorts from GEO and the model proved to be of high value. Conclusion: Our study shows that the 9-gene mRNAsi-related signature can predict the prognosis of LUAD patient and contribute to decisions in the treatment and prevention of LUAD patients.

In this paper, high-grade serous ovarian cancer (HGSOC) is studied, which is the most common histological subtype of ovarian cancer. We use a new analytical procedure to combine the bulk RNA-Seq sample for ovarian cancer, mRNA expression-based stemness index (mRNAsi), and single-cell data for ovarian cancer. Through integrating bulk RNA-Seq sample of cancer samples from TCGA, UCSC Xena and single-cell RNA-Seq (scRNA-Seq) data of HGSOC from GEO, and performing a series of computational analyses on them, we identify stemness markers and survival-related markers, explore stem cell populations in ovarian cancer, and provide potential treatment recommendation. As a result, 171 key genes for capturing stem cell characteristics are screened and one vital cancer stem cell subpopulation is identified. Through further analysis of these key genes and cancer stem cell subpopulation, more critical genes can be obtained as LCP2, FCGR3A, COL1A1, COL1A2, MT-CYB, CCT5, and PAPPA, are closely associated with ovarian cancer. So these genes have the potential to be used as prognostic biomarkers for ovarian cancer. [ABSTRACT FROM AUTHOR]

Abstract Background Cancer stem cells (CSCs) play an important role in drug resistance, recurrence, and metastasis of tumors. Considering the heterogeneity of tumors, this study aimed to explore the key genes regulating stem cells in intestinal-type and diffuse-type gastric cancer. Methods RNA-seq data and related clinical information were downloaded from The Cancer Genome Atlas (TCGA). WGCNA was used to clustered differentially expressed genes with similar expression profiles to form modules. Furtherly, based on the mRNA expression-based stemness index (mRNAsi), significant modules and key genes were identified. Next, the expression of key genes was further verified by the Oncomine database. Results MRNAsi scores of GC were significantly higher than that of normal tissue. Additionally, mRNAsi scores of intestinal-type GC (IGC) were significantly higher than that of diffuse-type GC (DGC). WGCNA showed that the blue module of IGC and the brown module of DGC were both the most significantly associated with mRNAsi. We screened out 16 and 43 key genes for IGC and DGC and found that these genes were closely related, respectively. Functional analysis showed the relationship between the key genes confirmed in the Oncomine database and the fate of cells. Conclusions In this study, 16 and 43 genes related to the characteristics of CSCs were identified in IGC and DGC, respectively. These genes were both associated with cell cycle, which could serve as therapeutic targets for the inhibition of stem cells from both types of GC.

Abstract Background With the gradual unveiling of tumour heterogeneity, cancer stem cells (CSCs) are now being considered the initial component of tumour initiation. However, the mechanisms of the growth and maintenance of breast cancer (BRCA) stem cells are still unknown. Methods To explore the crucial genes modulating BRCA stemness characteristics, we combined the gene expression value and mRNA expression-based stemness index (mRNAsi) of samples from The Cancer Genome Atlas (TCGA), and the mRNAsi was corrected using the tumour purity (corrected mRNAsi). mRNAsi and corrected mRNAsi were analysed and showed a close relationship with BRCA clinical characteristics, including tumour depth, pathological staging and survival status. Next, weighted gene co-expression network analysis (WGCNA) was applied to distinguish crucial gene modules and key genes. A series of functional analyses and expression validation of key genes were conducted using multiple databases, including Oncomine, Gene Expression Omnibus (GEO) and Gene Expression Profiling Integrative Analysis (GEPIA). Results This study found that mRNAsi and corrected mRNAsi scores were higher in BRCA tissues than that in normal tissues, and both of them increased with tumour stage. Higher corrected mRNAsi scores showed worse overall survival outcomes. We screened 3 modules and 32 key genes, and those key genes were found to be strongly correlated with each other. Functional analysis revealed that the key genes were related to cell fate decision events such as the cell cycle, cellular senescence, chromosome segregation and mitotic nuclear division. Among 32 key genes, we identified 12 genes that strongly correlated with BRCA survival. Conclusions Thirty-two genes were found to be closely related to BRCA stem cell characteristics; among them, 12 genes showed prognosis-oriented effects in BRCA patients. The most significant signalling pathway related to stemness in BRCA was the cell cycle pathway, which may support new ideas for screening therapeutic targets to inhibit BRCA stem characteristics. These findings may highlight some therapeutic targets for inhibiting BRCA stem cells.

In this paper, high-grade serous ovarian cancer (HGSOC) is studied, which is the most common histological subtype of ovarian cancer. We use a new analytical procedure to combine the bulk RNA-Seq sample for ovarian cancer, mRNA expression-based stemness index (mRNAsi), and single-cell data for ovarian cancer. Through integrating bulk RNA-Seq sample of cancer samples from TCGA, UCSC Xena and single-cell RNA-Seq (scRNA-Seq) data of HGSOC from GEO, and performing a series of computational analyses on them, we identify stemness markers and survival-related markers, explore stem cell populations in ovarian cancer, and provide potential treatment recommendation. As a result, 171 key genes for capturing stem cell characteristics are screened and one vital cancer stem cell subpopulation is identified. Through further analysis of these key genes and cancer stem cell subpopulation, more critical genes can be obtained as LCP2, FCGR3A, COL1A1, COL1A2, MT-CYB, CCT5, and PAPPA, are closely associated with ovarian cancer. So these genes have the potential to be used as prognostic biomarkers for ovarian cancer.

High stemness index scores are associated with poor survival in patients with lung cancer. Studies on the mRNA expression-based stemness index (mRNAsi) are typically conducted using tumor tissues; however, mRNAsi-related expression signatures based on cell-free RNA (cfRNA) are yet to be comprehensively investigated. The present study aimed to elucidate the gene expression profiles of tumor stemness in lung cancer tissues and corresponding cfRNAs in blood, and to assess their links with immune infiltration. Tumor tissue, paracancerous tissue, peripheral blood and lymph node samples were collected from patients with stage I-III non-small cell lung cancer and RNA sequencing was performed. The TCGAbiolinks package was used to calculate the mRNAsi for each of these four types of sample. Weighted gene co-expression network analysis and differentially expressed gene analyses were performed to investigate mRNAsi-related genes, and pathway enrichment analysis was performed using the Kyoto Encyclopedia of Genes and Genomes (KEGG) orthology-based annotation system. In addition, the STAR-Fusion tool was used to detect fusion variants, and CIBERSORT was used to analyze the correlations of stemness signatures in tissues and blood with immune cell infiltration. The mRNAsi values in peripheral blood and lymph nodes were found to be higher than those in cancer tissues. 'Hematopoietic cell lineage' was the only KEGG pathway enriched in mRNAsi-related genes in both lung cancer tissues and peripheral blood. In addition, the protein tyrosine phosphatase receptor type C associated protein gene was the only gene commonly associated with the mRNAsi in these two types of sample. The expression of mRNAsi-related genes was increased in the dendritic and Treg cells in tumor tissues, but was elevated in Treg and CD8 cells in the blood. In conclusion, cfRNAs in the blood exhibit unique stemness signatures that have potential for use in the diagnosis of lung cancer., Competing Interests: ZW, JL, and RW are employees of Berry Oncology Corporation, and some of the materials used in the study were provided by Berry Oncology Corporation. BL is an employee of Liaoning Kanghui Biotechnology Corporation, and assays and analyses in this study were provided free of charge by Liaoning Kanghui Biotechnology Corporation. The other authors declare that they have no competing interests., (Copyright: © 2024 Yan et al.)

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)., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2024 The Authors. Published by Elsevier Ltd.)

Background: The precise function of Tolloid Like 2 (TLL2) remains uncertain within the context of Lung Adenocarcinoma (LUAD)., Objective: The primary objective of this investigation was to conduct a thorough analysis., Methods: To assess its diagnostic utility, data from The Cancer Genome Atlas (TCGA) database were used to assess TLL2 expression in pan-cancer and LUAD. The study has also investigated the correlation between TLL2 expression levels and LUAD symptoms and prognosis. Furthermore, the study has explored possible regulatory networks involving TLL2, including its association with immune infiltration, tumor stemness index (mRNAsi), and drug sensitivity in LUAD. We have explored TLL2 expression in single-cell sequencing of LUAD and the genomic variation and clinical significance of TLL2 in LUAD. The expression of TLL2 has been validated in GSE87340 and cell lines by quantitative Real-time PCR (qRT-PCR)., Results: An abnormal expression of TLL2 has been found in pan-cancer and LUAD. In LUAD patients, elevated levels of TLL2 were significantly related to the T stage (p = 0.046) and the pathological stage (p = 0.016). The expression of TLL2 in patients with LUAD was significantly associated with poorer Overall Survival (OS) (p < 0.001). The expression of TLL2 was determined to be an independent predictor of poorer OS (p = 0.042). TLL2 was associated with ribosome, neuroactive ligand-receptor interaction, allograft rejection, ECM receptor interaction, asthma, porphyrin and chlorophyll metabolism, focal adhesion, pentose and glucuronate inter-conversions, and ascorbate and aldarate metabolism. The expression of TLL2 in LUAD was correlated with immune infiltration and mRNAsi. The expression of TLL2 was significantly and negatively correlated with TAK-715, XMD13-2, STF-62247, OSI-930, and EHT-1864 in LUAD. The TLL2 gene was up- -regulated in multiple individual LUAD cells. LUAD patients with altered TLL2 had a shorter PFS as opposed to those with unaltered TLL2. The expression of TLL2 was significantly increased in LUAD cells., Conclusion: For patients with LUAD, TLL2 may serve as an immunotherapeutic target and a useful prognosis biomarker., (Copyright© Bentham Science Publishers; For any queries, please email at epub@benthamscience.net.)

Glioblastoma (GBM) is the most common glial tumour and has extremely poor prognosis. GBM stem-like cells drive tumorigenesis and progression. However, a systematic assessment of stemness indices and their association with immunological properties in GBM is lacking. We collected 874 GBM samples from four GBM cohorts (TCGA, CGGA, GSE4412, and GSE13041) and calculated the mRNA expression-based stemness indices (mRNAsi) and corrected mRNAsi (c_mRNAsi, mRNAsi/tumour purity) with OCLR algorithm. Then, mRNAsi/c_mRNAsi were used to quantify the stemness traits that correlated significantly with prognosis. Additionally, confounding variables were identified. We used discrimination, calibration, and model improvement capability to evaluate the established models. Finally, the CIBERSORTx algorithm and ssGSEA were implemented for functional analysis. Patients with high mRNAsi/c_mRNAsi GBM showed better prognosis among the four GBM cohorts. After identifying the confounding variables, c_mRNAsi still maintained its prognostic value. Model evaluation showed that the c_mRNAsi-based model performed well. Patients with high c_mRNAsi exhibited significant immune suppression. Moreover, c_mRNAsi correlated negatively with infiltrating levels of immune-related cells. In addition, ssGSEA revealed that immune-related pathways were generally activated in patients with high c_mRNAsi. We comprehensively evaluated GBM stemness indices based on large cohorts and established a c_mRNAsi-based classifier for prognosis prediction.

S-palmitoylation, catalyzed by a family of 23 zinc finger Asp-His-His-Cys (DHHC) domain-containing (ZDHHC) protein acyltransferases localized on the cell membrane. However, stemness genes modulated by ZDHHCs in lung adenocarcinoma (LUAD) remain to be defined. Previously, we have constructed a network of cancer stem cell genes, including INCENP, based on mRNA stemness indices (mRNAsi) of LUAD. INCENP has the function of a chromosomal passenger complex locating to centromeres, which is performed by the conserved region of its N-terminal domain. INCENP protein with a deletion of the first non-conserved 26 amino acid sequence failed to target centromeres. However, the exact function of the deleted sequence has not been elucidated. To identify novel cancer stem cell-relevant palmitoylated proteins and responsible ZDHHC enzymes in LUAD, we analyzed multi-omics data obtained from the database of The Cancer Genome Atlas (TCGA), Gene Expression Omnibus (GEO), Clinical Proteomic Tumor Analysis Consortium (CPTAC), and the Human Protein Atlas (HPA). ZDHHC5 is distinguished from the ZDHHC family for being up-regulated in mRNA and protein levels and associated with malignant prognosis. ZDHHC5 was positively associated with INCENP, and the correlation score increased with LUAD stages. CSS-Palm results showed Cys15 was the S-palmitoylation site of INCENP. Interestingly, Cys15 locates in the 1–26 aa sequence of INCENP, and is a conserved site across species. As INCENP is a nuclear protein, we predicted that the nuclear localization signal of ZDHHC5 was specific to the importin αβ pathway, and the result of immunofluorescence proves that ZDHHC5 is located in the nucleoplasm, in addition to the plasma membrane. Therefore, our study indicates the S-palmitoylation of INCENP mediated by ZDHHC5 as a potential mechanism of S-palmitoylation to modulate CSCs in LUAD.

Tumor progression includes the obtainment of progenitor and stem cell-like features and the gradual loss of a differentiated phenotype. Stemness was defined as the potential for differentiation and self-renewal from the cell of origin. Previous studies have confirmed the effective application of stemness in a number of malignancies. However, the mechanisms underlying the growth and maintenance of multiple myeloma (MM) stem cells remain unclear. We calculated the stemness index for samples of MM by utilizing a novel one-class logistic regression (OCLR) machine learning algorithm and found that mRNA expression-based stemness index (mRNAsi) was an independent prognostic factor of MM. Based on the same cutoff value, mRNAsi could stratify MM patients into low and high groups with different outcomes. We identified 127 stemness-related signatures using weighted gene co-expression network analysis (WGCNA) and differential expression analysis. Functional annotation and pathway enrichment analysis indicated that these genes were mainly involved in the cell cycle, cell differentiation, and DNA replication and repair. Using the molecular complex detection (MCODE) algorithm, we identified 34 pivotal signatures. Meanwhile, we conducted unsupervised clustering and classified the MM cohorts into three MM stemness (MMS) clusters with distinct prognoses. Samples in MMS-cluster3 possessed the highest stemness fractions and the worst prognosis. Additionally, we applied the ESTIMATE algorithm to infer differential immune infiltration among the three MMS clusters. The immune core and stromal score were significantly lower in MMS-cluster3 than in the other clusters, supporting the negative relation between stemness and anticancer immunity. Finally, we proposed a prognostic nomogram that allows for individualized assessment of the 3- and 5-year overall survival (OS) probabilities among patients with MM. Our study comprehensively assessed the MM stemness index based on large cohorts and built a 34-gene based classifier for predicting prognosis and potential strategies for stemness treatment.

Glioblastoma (GBM) is the most malignant and lethal intracranial tumor, with extremely limited treatment options. Immunotherapy has been widely studied in GBM, but none can significantly prolong the overall survival (OS) of patients without selection. Considering that GBM cancer stem cells (CSCs) play a non-negligible role in tumorigenesis and chemoradiotherapy resistance, we proposed a novel stemness-based classification of GBM and screened out certain population more responsive to immunotherapy. The one-class logistic regression algorithm was used to calculate the stemness index (mRNAsi) of 518 GBM patients from The Cancer Genome Atlas (TCGA) database based on transcriptomics of GBM and pluripotent stem cells. Based on their stemness signature, GBM patients were divided into two subtypes via consensus clustering, and patients in Stemness Subtype I presented significantly better OS but poorer progression-free survival than Stemness Subtype II. Genomic variations revealed patients in Stemness Subtype I had higher somatic mutation loads and copy number alteration burdens. Additionally, two stemness subtypes had distinct tumor immune microenvironment patterns. Tumor Immune Dysfunction and Exclusion and subclass mapping analysis further demonstrated patients in Stemness Subtype I were more likely to respond to immunotherapy, especially anti-PD1 treatment. The pRRophetic algorithm also indicated patients in Stemness Subtype I were more resistant to temozolomide therapy. Finally, multiple machine learning algorithms were used to develop a 7-gene Stemness Subtype Predictor, which were further validated in two external independent GBM cohorts. This novel stemness-based classification could provide a promising prognostic predictor for GBM and may guide physicians in selecting potential responders for preferential use of immunotherapy. [ABSTRACT FROM AUTHOR]

Tumor progression includes the obtainment of progenitor and stem cell-like features and the gradual loss of a differentiated phenotype. Stemness was defined as the potential for differentiation and self-renewal from the cell of origin. Previous studies have confirmed the effective application of stemness in a number of malignancies. However, the mechanisms underlying the growth and maintenance of multiple myeloma (MM) stem cells remain unclear. We calculated the stemness index for samples of MM by utilizing a novel one-class logistic regression (OCLR) machine learning algorithm and found that mRNA expression-based stemness index (mRNAsi) was an independent prognostic factor of MM. Based on the same cutoff value, mRNAsi could stratify MM patients into low and high groups with different outcomes. We identified 127 stemness-related signatures using weighted gene co-expression network analysis (WGCNA) and differential expression analysis. Functional annotation and pathway enrichment analysis indicated that these genes were mainly involved in the cell cycle, cell differentiation, and DNA replication and repair. Using the molecular complex detection (MCODE) algorithm, we identified 34 pivotal signatures. Meanwhile, we conducted unsupervised clustering and classified the MM cohorts into three MM stemness (MMS) clusters with distinct prognoses. Samples in MMS-cluster3 possessed the highest stemness fractions and the worst prognosis. Additionally, we applied the ESTIMATE algorithm to infer differential immune infiltration among the three MMS clusters. The immune core and stromal score were significantly lower in MMS-cluster3 than in the other clusters, supporting the negative relation between stemness and anticancer immunity. Finally, we proposed a prognostic nomogram that allows for individualized assessment of the 3- and 5-year overall survival (OS) probabilities among patients with MM. Our study comprehensively assessed the MM stemness index based on large cohorts and built a 34-gene based classifier for predicting prognosis and potential strategies for stemness treatment. [ABSTRACT FROM AUTHOR]

Background: The relationship between long noncoding RNAs (lncRNAs) and the mRNA stemness index (mRNAsi) in colorectal cancer (CRC) is still unclear. Materials & methods: The mRNAsi, mRNAsi-related lncRNAs and their clinical significance were analyzed by bioinformatic approaches in The Cancer Genome Atlas (TCGA)-COREAD dataset. Results: mRNAsi was negatively related to pathological features but positively related to overall survival and recurrence-free survival in CRC. A five mRNAsi-related lncRNAs prognostic signature was further developed and showed independent prognostic factors related to overall survival in CRC patients, due to the five mRNAsi-related lncRNAs involved in several pathways of the cancer stem cells and malignant cancer cell phenotypes. Conclusion: The present study highlights the potential roles of mRNAsi-related lncRNAs as alternative prognostic markers. [ABSTRACT FROM AUTHOR]

Background: Stem cells play an important role in hepatocellular carcinoma (HCC). However, their precise effect on HCC tumorigenesis and progression remains unclear. The present study aimed to characterize stem cell-related gene expression in HCC. Methods: The mRNA expression-based stemness index (mRNAsi) was used to analyze the clinical characteristics and prognosis of HCC patients. The weighted gene co-expression network analysis (WGCNA) was used to construct a gene co-expression network of 374 HCC patients. Finally, six genes were used to construct the prognosis signature. Results: HCC patients had a higher mRNAsi score than healthy people, suggesting poor prognosis. Two gene modules highly related to mRNAsi were identified. Multivariate Cox analysis was carried out to establish a Cox proportional risk regression model. The risk score for each patient was the sum of the product of each gene expression and its coefficient. Survival analysis suggested that the low-risk group had a significantly better prognosis. Conclusions: The established six-gene signature was able to predict patient prognosis accurately. This new signature should be verified in prospective studies in order to determine patient prognosis in clinical decision-making. [ABSTRACT FROM AUTHOR]

Prostate cancer stemness (PCS) cells have been reported to drive tumor progression, recurrence and drug resistance. However, there is lacking systematical assessment of stemlike indices and associations with immunological properties in prostate adenocarcinoma (PRAD). We thus collected 7 PRAD cohorts with 1465 men and calculated the stemlike indices for each sample using one-class logistic regression machine learning algorithm. We selected the mRNAsi to quantify the stemlike indices that correlated significantly with prognosis and accordingly identified 21 PCS-related CpG loci and 13 pivotal signature. The 13-gene based PCS model possessed high predictive significance for progression-free survival (PFS) that was trained and validated in 7 independent cohorts. Meanwhile, we conducted consensus clustering and classified the total cohorts into 5 PCS clusters with distinct outcomes. Samples in PCScluster5 possessed the highest stemness fractions and suffered from the worst prognosis. Additionally, we implemented the CIBERSORT algorithm to infer the differential abundance across 5 PCS clusters. The activated immune cells (CD8+ T cell and dendritic cells) infiltrated significantly less in PCScluster5 than other clusters, supporting the negative regulations between stemlike indices and anticancer immunity. High mRNAsi was also found to be associated with up-regulation of immunosuppressive checkpoints, like PDL1. Lastly, we used the Connectivity Map (CMap) resource to screen potential compounds for targeting PRAD stemness, including the top hits of cell cycle inhibitor and FOXM1 inhibitor. Taken together, our study comprehensively evaluated the PRAD stemlike indices based on large cohorts and established a 13-gene based classifier for predicting prognosis or potential strategies for stemness treatment. [ABSTRACT FROM AUTHOR]

BackgroundBladder urothelial carcinoma (BLCA) is one of the most common urinary system malignancies with a high metastasis rate. Cancer stem cells (CSCs) play an important role in the occurrence and progression of BLCA, however, its roles in bone metastasis and the prognostic stemness biomarkers have not been identified in BLCA.MethodIn order to identify the roles of CSC in the tumorigenesis, bone metastasis and prognosis of BLCA, the RNA sequencing data of patients with BLCA were retrieved from The Cancer Genome Atlas (TCGA) databases. The mRNA expression-based stemness index (mRNAsi) and the differential expressed genes (DEGs) were evaluated and identified. The associations between mRNAsi and the tumorigenesis, bone metastasis, clinical stage and overall survival (OS) were also established. The key prognostic stemness-related genes (PSRGs) were screened by Lasso regression, and based on them, the predict model was constructed. Its accuracy was tested by the area under the curve (AUC) of the receiver operator characteristic (ROC) curve and the risk score. Additionally, in order to explore the key regulatory network, the relationship among differentially expressing TFs, PSRGs, and absolute quantification of 50 hallmarks of cancer were also identified by Pearson correlation analysis. To verify the identified key TFs and PSRGs, their expression levels were identified by our clinical samples via immunohistochemistry (IHC).ResultsA total of 8,647 DEGs were identified between 411 primary BLCAs and 19 normal solid tissue samples. According to the clinical stage, mRNAsi and bone metastasis of BLCA, 2,383 stage-related DEGs, 3,680 stemness-related DEGs and 716 bone metastasis-associated DEGs were uncovered, respectively. Additionally, compared with normal tissue, mRNAsi was significantly upregulated in the primary BLCA and also associated with the prognosis (P = 0.016), bone metastasis (P < 0.001) and AJCC clinical stage (P < 0.001) of BLCA patients. A total of 20 PSRGs were further screened by Lasso regression, and based on them, we constructed the predict model with a relatively high accuracy (AUC: 0.699). Moreover, we found two key TFs (EPO, ARID3A), four key PRSGs (CACNA1E, LINC01356, CGA and SSX3) and five key hallmarks of cancer gene sets (DNA repair, myc targets, E2F targets, mTORC1 signaling and unfolded protein response) in the regulatory network. The tissue microarray of BLCA and BLCA bone metastasis also revealed high expression of the key TFs (EPO, ARID3A) and PRSGs (SSX3) in BLCA.ConclusionOur study identifies mRNAsi as a reliable index in predicting the tumorigenesis, bone metastasis and prognosis of patients with BLCA and provides a well-applied model for predicting the OS for patients with BLCA based on 20 PSRGs. Besides, we also identified the regulatory network between key PSRGs and cancer gene sets in mediating the BLCA bone metastasis.

Background: Bladder urothelial carcinoma (BLCA) is one of the most common urinary system malignancies with a high metastasis rate. Cancer stem cells (CSCs) play an important role in the occurrence and progression of BLCA, however, its roles in bone metastasis and the prognostic stemness biomarkers have not been identified in BLCA. Method: In order to identify the roles of CSC in the tumorigenesis, bone metastasis and prognosis of BLCA, the RNA sequencing data of patients with BLCA were retrieved from The Cancer Genome Atlas (TCGA) databases. The mRNA expression-based stemness index (mRNAsi) and the differential expressed genes (DEGs) were evaluated and identified. The associations between mRNAsi and the tumorigenesis, bone metastasis, clinical stage and overall survival (OS) were also established. The key prognostic stemness-related genes (PSRGs) were screened by Lasso regression, and based on them, the predict model was constructed. Its accuracy was tested by the area under the curve (AUC) of the receiver operator characteristic (ROC) curve and the risk score. Additionally, in order to explore the key regulatory network, the relationship among differentially expressing TFs, PSRGs, and absolute quantification of 50 hallmarks of cancer were also identified by Pearson correlation analysis. To verify the identified key TFs and PSRGs, their expression levels were identified by our clinical samples via immunohistochemistry (IHC). Results: A total of 8,647 DEGs were identified between 411 primary BLCAs and 19 normal solid tissue samples. According to the clinical stage, mRNAsi and bone metastasis of BLCA, 2,383 stage-related DEGs, 3,680 stemness-related DEGs and 716 bone metastasis-associated DEGs were uncovered, respectively. Additionally, compared with normal tissue, mRNAsi was significantly upregulated in the primary BLCA and also associated with the prognosis (P = 0.016), bone metastasis (P < 0.001) and AJCC clinical stage (P < 0.001) of BLCA patients. A total of 20 PSRGs were further screened by Lasso regression, and based on them, we constructed the predict model with a relatively high accuracy (AUC: 0.699). Moreover, we found two key TFs (EPO, ARID3A), four key PRSGs (CACNA1E, LINC01356, CGA and SSX3) and five key hallmarks of cancer gene sets (DNA repair, myc targets, E2F targets, mTORC1 signaling and unfolded protein response) in the regulatory network. The tissue microarray of BLCA and BLCA bone metastasis also revealed high expression of the key TFs (EPO, ARID3A) and PRSGs (SSX3) in BLCA. Conclusion: Our study identifies mRNAsi as a reliable index in predicting the tumorigenesis, bone metastasis and prognosis of patients with BLCA and provides a well-applied model for predicting the OS for patients with BLCA based on 20 PSRGs. Besides, we also identified the regulatory network between key PSRGs and cancer gene sets in mediating the BLCA bone metastasis. [ABSTRACT FROM AUTHOR]

Objective: An mRNA expression-based stemness index (mRNAsi) has been developed to characterize cancer stemness. However, the predictive value of mRNAsi-based signature in therapeutic resistance and immunotherapy in thyroid cancer (THCA) remains unclarified. This study evaluated and validated the role of mRNAsi in drug sensitivity, its relationship between mRNAsi and THCA clinical features and immunity based on bioinformatics., Methods: Based on transcriptome data of THCA patients from the Tumor Genome Atlas Project (TCGA) database, and expression data of multifunctional stem cell samples from the Progenitor Cell Biology Consortium (PCBC) databases, mRNAsi was calculated by the " one class logistic regression (OCLR)" method, Molecular subtypes of TCGA-THCA samples were identified with mRNAsi-related genes using ConsensusClusterPlus method. The gene mutation, clinical characteristics, immune characteristics, TIDE and drug sensitivity were compared among molecular subtypes. A prognostic model was designed with Lasso cox method. Modulation of malignant phenotype of THCA cell lines by model characterization genes is validated by CCK-8, flow cytometry. DNA methylation disorder in promoter region was analyzed between risk groups. The model was validated for survival in the internal Test dataset, while TCGA pan-cancer and immunotherapy datasets were further employed to validate the performance of this model., Results: We obtained a total of 78 stem cell samples, each containing the expression profile of 8087 mRNA genes. Based on mRNAsi, THCA was divided into 3 subtypes. Subtype C2 had the poorest prognosis and highest immune score, while subtype C3 had the best prognosis, lowest mRANsi and highest TIDE score. Patients in subtype C2 showed higher sensitivity to Cisplatin, Erlotinib, Paclitaxel, and Lapatinib. The prognostic signature was generated using 5 mRNAsi-related genes, which could predict prognosis for THCA. qRT-PCR results showed that the expression of 5 genes were various in Hth7 and KTC-1 cells, and inhibition CELSR3 expression increased percentage of apoptosis in Hth7 and KTC-1 cells. mRNAsi related DNA methylation sites were mainly enriched in tumor related pathways. Good performance of this model was validated in Test dataset, pan-cancer and immunotherapy datasets., Conclusion: This study identified three subtypes for classification and developed a prognostic model with mRNAsi-related genes, which provided great potential for prognosis and immunotherapy prediction., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2024 The Authors.)

Background: CESC is the second most commonly diagnosed gynecological malignancy. Given the pivotal involvement of metabolism-related genes (MRGs) in the etiology of multiple tumors, our investigation aims to devise a prognostic risk signature rooted in cancer stemness and metabolism., Methods: The stemness index based on mRNA expression (mRNAsi) of samples from the TCGA dataset was computed using the One-class logistic regression (OCLR) algorithm. Furthermore, potential metabolism-related genes related to mRNAsi were identified through weighted gene co-expression network analysis (WGCNA). We construct a stemness-related metabolic gene signature through shrinkage estimation and univariate analysis, thereby calculating the corresponding risk scores. Moreover, we selected corresponding DEGs between groups with high- and low-risk score and conducted routine bioinformatic analyses. Furthermore, we validated the expression of four hub genes at the protein level through immunohistochemistry (IHC) in samples obtained from our patient cohort., Results: According to the findings, it was found that six genes-AKR1B10, GNA15, ALDH1B1, PLOD2, LPCAT1, and GPX8- were differentially expressed in both TCGA-CSEC and GEO datasets among 23 differentially expressed metabolism-related genes (DEMRGs). mRNAsi exhibited a notable association with the extent of key oncogene mutation. The results showed that the AUC values for forecasting survival at 1, 3, and 5 years are 0.715, 0.689, and 0.748, individually. We observed a notable association between the risk score and different immune cell populations, along with enrichment in crucial signaling pathways in CESC. Four genes differentially expressed between different risk score groups were validated by IHC to be highly expressed in the CESC samples at the protein level., Conclusion: The current investigation indicated that a 3-gene signature based on stemness-related metabolic and 4 hub genes with differential expression between high and low-risk score subgroups may serve as valuable prognostic markers and potential therapeutic targets in CESC.

BackgroundWith characteristic self-renewal and multipotent differentiation, cancer stem cells (CSCs) have a crucial influence on the metastasis, relapse and drug resistance of gastric cancer (GC). However, the genes that participates in the stemness of GC stem cells have not been identified.MethodsThe mRNA expression-based stemness index (mRNAsi) was analyzed with differential expressions in GC. The weighted gene co-expression network analysis (WGCNA) was utilized to build a co-expression network targeting differentially expressed genes (DEG) and discover mRNAsi-related modules and genes. We assessed the association between the key genes at both the transcription and protein level. Gene Expression Omnibus (GEO) database was used to validate the expression levels of the key genes. The risk model was established according to the least absolute shrinkage and selection operator (LASSO) Cox regression analysis. Furthermore, we determined the prognostic value of the model by employing Kaplan-Meier (KM) plus multivariate Cox analysis.ResultsGC tissues exhibited a substantially higher mRNAsi relative to the healthy non-tumor tissues. Based on WGCNA, 17 key genes (ARHGAP11A, BUB1, BUB1B, C1orf112, CENPF, KIF14, KIF15, KIF18B, KIF4A, NCAPH, PLK4, RACGAP1, RAD54L, SGO2, TPX2, TTK, and XRCC2) were identified. These key genes were clearly overexpressed in GC and validated in the GEO database. The protein-protein interaction (PPI) network as assessed by STRING indicated that the key genes were tightly connected. After LASSO analysis, a nine-gene risk model (BUB1B, NCAPH, KIF15, RAD54L, KIF18B, KIF4A, TTK, SGO2, C1orf112) was constructed. The overall survival in the high-risk group was relatively poor. The area under curve (AUC) of risk score was higher compared to that of clinicopathological characteristics. According to the multivariate Cox analysis, the nine-gene risk model was a predictor of disease outcomes in GC patients (HR, 7.606; 95% CI, 3.037–19.051; P < 0.001). We constructed a prognostic nomogram with well−fitted calibration curve based on risk score and clinical data.ConclusionThe 17 mRNAsi-related key genes identified in this study could be potential treatment targets in GC treatment, considering that they can inhibit the stemness properties. The nine-gene risk model can be employed to predict the disease outcomes of the patients.

Background: Cancer stem cells (CSCs) play an important role in drug resistance, recurrence, and metastasis of tumors. Considering the heterogeneity of tumors, this study aimed to explore the key genes regulating stem cells in intestinal-type and diffuse-type gastric cancer.Methods: RNA-seq data and related clinical information were downloaded from The Cancer Genome Atlas (TCGA). WGCNA was used to clustered differentially expressed genes with similar expression profiles to form modules. Furtherly, based on the mRNA expression-based stemness index (mRNAsi), significant modules and key genes were identified. Next, the expression of key genes was further verified by the Oncomine database.Results: MRNAsi scores of GC were significantly higher than that of normal tissue. Additionally, mRNAsi scores of intestinal-type GC (IGC) were significantly higher than that of diffuse-type GC (DGC). WGCNA showed that the blue module of IGC and the brown module of DGC were both the most significantly associated with mRNAsi. We screened out 16 and 43 key genes for IGC and DGC and found that these genes were closely related, respectively. Functional analysis showed the relationship between the key genes confirmed in the Oncomine database and the fate of cells.Conclusions: In this study, 16 and 43 genes related to the characteristics of CSCs were identified in IGC and DGC, respectively. These genes were both associated with cell cycle, which could serve as therapeutic targets for the inhibition of stem cells from both types of GC. [ABSTRACT FROM AUTHOR]

Background: Head and neck squamous cell carcinoma (HNSCC) is a type of invasive malignancy and the seventh most common cancer in the worldwide. Cancer stem cells (CSCs) are self‐renewal cells in tumors and can produce heterogeneous tumor cells, which play an important role in the development of HNSCC. In our research, we aimed to identify genes related to the CSCs characteristics in HNSCC. Methods: Messenger RNA (mRNA) expression‐based stiffness index (mRNAsi) can be used as a quantitative characterization of CSCs. We used one‐class logistic regression machine learning algorithm (OCLR) to calculate the mRNAsi and investigate the relationship between mRNAsi and clinical characteristics of HNSCC. Then, a weighted gene co‐expression network analysis (WGCNA) and protein‐protein interaction (PPI) network was constructed to screen hub genes related to mRNAsi of HNSCC. Results: The results indicated that the score of mRNAsi in HNSCC tissues is higher than in paracancer tissues, while the mRNAsi was not statistically correlated with the prognosis and clinical characteristics of HNSCC. Six positive and six negative hub genes related to mRNAsi of HNSCC were selected, which may act as therapeutic targets for inhibiting CSCs within HNSCC. Conclusions: In conclusion, our research selected 12 hub genes related to mRNAsi of HNSCC through weighted gene co‐expression network analysis. These genes may become therapeutic targets to inhibit the CSCs of HNSCC in the future. [ABSTRACT FROM AUTHOR]

Background: Accounting for tumor heterogeneity, cancer stem cells (CSC) are involved in tumor metastasis, relapse, and drug resistance. Genes regulating CSC characteristics in lung adenocarcinoma (LUAD) were explored and validated in this study. Methods: The mRNA stemness index (mRNAsi) of more than 500 LUAD cases from The Cancer Genome Atlas database were calculated using a one-class logistic regression machine learning algorithm based on the mRNA expression of pluripotent stem cells and their differentiated progeny. mRNAsi-related key genes were identified by weighted correlation network analysis. The expression levels and prognostic roles of key genes were analyzed in Oncomine, PrognoScan, and Kaplan–Meier plotter databases, and validated using data from our center. Results: The mRNAsi was significantly higher in LUAD compared with normal lung tissues. LUAD patients of advanced stage exhibited a higher mRNAsi and worse overall survival (OS). Eight key genes were identified: heat shock 70 kDa protein 4 (HSPA4), cell division cycle associated 7 (CDCA7), cell division cycle 20 (CDC20), cyclin-dependent kinase 1 (CDK1), CAP-GLY domain containing linker protein 1 (CLIP1), cyclin B1 (CCNB1), H2A histone family, member X (H2AFX), and Bloom syndrome, RecQ helicase-like (BLM). These genes were differentially expressed in various types of malignancies and validated in the LUAD cases. LUAD patients with low expression of CDC20, CDK1, CCNB1, H2AFX, or BLM had a significantly better OS, whereas OS was reduced for patients with low expression of CLIP1. In addition, the expression of CDCA7 did not significantly impact the OS of LUAD patients. The protein–protein interaction networks evaluated by STRING demonstrated strong relationships between these key genes, which were validated in our cases. Conclusions: The mRNAsi was significantly higher in LUAD compared with normal samples. Eight mRNAsi-related key genes were associated with prognosis and the cell cycle, and were strongly correlated with each other and differentially expressed in tumor and normal samples. We provide a new strategy for exploring stemness-related genes in LUAD cases. [ABSTRACT FROM AUTHOR]

Background: With the gradual unveiling of tumour heterogeneity, cancer stem cells (CSCs) are now being considered the initial component of tumour initiation. However, the mechanisms of the growth and maintenance of breast cancer (BRCA) stem cells are still unknown.Methods: To explore the crucial genes modulating BRCA stemness characteristics, we combined the gene expression value and mRNA expression-based stemness index (mRNAsi) of samples from The Cancer Genome Atlas (TCGA), and the mRNAsi was corrected using the tumour purity (corrected mRNAsi). mRNAsi and corrected mRNAsi were analysed and showed a close relationship with BRCA clinical characteristics, including tumour depth, pathological staging and survival status. Next, weighted gene co-expression network analysis (WGCNA) was applied to distinguish crucial gene modules and key genes. A series of functional analyses and expression validation of key genes were conducted using multiple databases, including Oncomine, Gene Expression Omnibus (GEO) and Gene Expression Profiling Integrative Analysis (GEPIA).Results: This study found that mRNAsi and corrected mRNAsi scores were higher in BRCA tissues than that in normal tissues, and both of them increased with tumour stage. Higher corrected mRNAsi scores showed worse overall survival outcomes. We screened 3 modules and 32 key genes, and those key genes were found to be strongly correlated with each other. Functional analysis revealed that the key genes were related to cell fate decision events such as the cell cycle, cellular senescence, chromosome segregation and mitotic nuclear division. Among 32 key genes, we identified 12 genes that strongly correlated with BRCA survival.Conclusions: Thirty-two genes were found to be closely related to BRCA stem cell characteristics; among them, 12 genes showed prognosis-oriented effects in BRCA patients. The most significant signalling pathway related to stemness in BRCA was the cell cycle pathway, which may support new ideas for screening therapeutic targets to inhibit BRCA stem characteristics. These findings may highlight some therapeutic targets for inhibiting BRCA stem cells. [ABSTRACT FROM AUTHOR]

Objective: Cancer stem cells are self-renewal cells in tumors and can produce heterogeneous tumor cells, which play an important role in the development of lung squamous cell carcinoma (LSCC). In our research, we aimed to explore the expression of genes related to LSCC stem cells.Methods: We downloaded the RNAseq data, the pathological and prognostic profiles of LSCC cases from the public database TCGA. The mRNA expression-based stiffness index (mRNAsi) of LSCC was calculated and the prognostic value of mRNAsi was discussed. Then, we constructed a weighted gene co-expression network analysis (WGCNA) to screen key genes related to mRNAsi of LSCC.Results: MRNAsi is an independent prognostic factor in LSCC. We screened 5 key genes (BUB1, BIRC5, CCNB2, KIF15 and SPAG5) related to mRNAsi of LSCC based on WGCNA. The key genes were highly expressed in the tumor samples compared to the normal samples. In addition, there is a strong interaction between proteins of these key genes and a strong co-expression relationship at the transcriptional level.Conclusions: To conclude, mRNAsi play an important role in LSCC. Five key genes (BUB1, BIRC5, CCNB2, KIF15 and SPAG5) related to mRNAsi were screened, which may act as therapeutic targets for inhibiting the stem cell characteristics of LSCC. [ABSTRACT FROM AUTHOR]