Dataset Information

Identification of prognostic stemness biomarkers in colon adenocarcinoma drug resistance.

ABSTRACT:

Background

Colon adenocarcinoma (COAD) is one of the leading causes of death worldwide. Cancer stem cells (CSCs) are vital for COAD chemoresistance and recurrence, however little is known about stem cell-related biomarkers in drug resistance and COAD prognosis prediction.

Methods

To uncover the roles of CSC in COAD tumorigenesis, chemoresistance, and prognosis, we retrieved COAD patients' RNAseq data from TCGA (The Cancer Genome Atlas). We further performed analysis of differentially expressed genes (DEGs) and mRNA expression-based stemness index (mRNAsi) to identify stemness-related COAD biomarkers. We then evaluated the roles of mRNAsi in tumorigenesis, clinical-stage, overall survival (OS), and chemoresistance. Afterward, we used identified prognostic stemness-related genes (PSRGs) to construct a prediction model. After constructing the prediction model, we used elastic Net regression and area under the curve (AUC) to explore the prediction value of PSRGs based on risk scores and the receiver operator characteristic (ROC) curve. To elucidate the underlying interconnected systems, we examined relationships between the levels of TFs, PSRGs, and 50 cancer hallmarks by a Pearson correlation analysis.

Results

Twelve thousand one hundred eight DEGs were identified by comparing 456 primary COADs and 41 normal solid tissue samples. Furthermore, we identified 4351 clinical stage-related DEGs, 16,516 stemness-associated DEGs, and 54 chemoresistance-related DEGs from cancer stages: mRNAsi, and COAD chemoresistance. Compared to normal tissue samples, mRNAsi in COAD patients were marked on an elevation and involved in prognosis (p = 0.027), stemness-related DEGs based on chemoresistance (OR = 3.28, p ≤ 0.001) and AJCC clinical stage relating (OR = 4.02, p ≤ 0.001) to COAD patients. The prediction model of prognosis were constructed using the 6 PSRGs with high accuracy (AUC: 0.659). The model identified universal correlation between NRIP2 and FDFT1 (key PRSGs), and some cancer related transcription factors (TFs) and trademarks of cancer gene were in the regulatory network.

Conclusion

We found that mRNAsi is a reliable predictive biomarker of tumorigenesis and COAD prognosis. Our established prediction model of COAD chemoresistance, which includes the six PSRGs, is effective, as the model provides promising therapeutic targets in the COAD.

SUBMITTER: Li Z

PROVIDER: S-EPMC9261069 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description: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.

Project description:Colorectal cancer is one of the commoner digestive tract malignant tumor types, and its incidence and mortality rate are high. Accumulating evidence indicates that long?chain non?coding RNAs (lncRNAs) and protein?coding RNAs interact with each other by competing with the same micro(mi)RNA response element (MREs) and serve an important role in the regulation of gene expression in a variety of tumor types. However, the regulatory mechanism and prognostic role of lncRNA?mediated competing endogenous (ce)RNA networks in colon cancer have yet to be elucidated. The expression profiles of mRNAs, lncRNAs and miRNAs from 471 colon cancer and 41 paracancerous tissue samples were downloaded from The Cancer Genome Atlas database. A lncRNA?miRNA?mRNA ceRNA network in colon cancer was constructed and comprised 17 hub lncRNAs, 87 hub miRNA and 144 hub mRNAs. The topological properties of the network were analyzed, and the random walk algorithm was used to identify the nodes significantly associated with colon cancer. Survival analysis using the UALCAN database indicated that 2/17 lncRNAs identified [metastasis?associated lung adenocarcinoma transcript (MALAT1) and maternally expressed gene 3 (MEG3)] and 5/144 mRNAs [FES upstream region (FURIN), nuclear factor of activated T?cells 5 (NFAT5), RNA Binding Motif Protein 12B (RBM12B), Ras related GTP binding A (RRAGA) and WD repeat domain phosphoinositide?interacting protein 2 (WIPI2)] were significantly associated with the overall survival of patients with colon cancer, and may therefore be used as potential prognostic biomarkers of colon cancer. According to extracted lncRNA?miRNA?mRNA interaction pairs, the GSE26334 dataset was used to confirm that the lncRNA MALAT1/miR?129?5p/NFAT5 axis may represent a novel regulatory mechanism concerning the progression of colon cancer. The clusterProfiler package was used to analyze Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways in colon cancer. Finally, drugs that significantly interact with the core genes identified in colon cancer were predicted using a hypergeometric test. Of these, fostamatinib was identified to be a targeted drug for colon cancer therapy. The present findings provide a novel perspective for improved understanding of the lncRNA?associated ceRNA network and may facilitate the development of novel targeted therapeutics in colon cancer.

Project description: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.