Dataset Information

System analysis based on the ER stress-related genes identifies WFS1 as a novel therapy target for colon cancer.

ABSTRACT:

Background

Colon cancer (COAD) is the third-largest common malignant tumor and the fourth major cause of cancer death in the world. Endoplasmic reticulum (ER) stress has a great influence on cell growth, migration, proliferation, invasion, angiogenesis, and chemoresistance of massive tumors. Although ER stress is known to play an important role in various types of cancer, the prognostic model based on ER stress-related genes (ERSRGs) in colon cancer has not been constructed yet. In this study, we established an ERSRGs prognostic risk model to assess the survival of COAD patients.

Methods

The COAD gene expression profile and clinical information data of the training set were obtained from the GEO database (GSE40967) and the test set COAD gene expression profile and clinical informative data were downloaded from the TCGA database. The endoplasmic reticulum stress-related genes (ERSRGs) were obtained from Gene Set Enrichment Analysis (GSEA) website. Differentially expressed ERSRGs between normal samples and COAD samples were identified by R "limma" package. Based on the univariate, lasso, and multivariate Cox regression analysis, we developed an ERSRGs prognostic risk model to predict survival in COAD patients. Finally, we verified the function of WFS1 in COAD through in vitro experiments.

Results

We built a 9-gene prognostic risk model based on the univariate, lasso, and multivariate Cox regression analysis. Kaplan-Meier survival analysis and Receiver operating characteristic (ROC) curve revealed that the prognostic risk model has good predictive performance. Subsequently, we screened 60 compounds with significant differences in the estimated half-maximal inhibitory concentration (IC50) between high-risk and low-risk groups. In addition, we found that the ERSRGs prognostic risk model was related to immune cell infiltration and the expression of immune checkpoint molecules. Finally, we determined that knockdown of the expression of WFS1 inhibits the proliferation of colon cancer cells.

Conclusions

The prognostic risk model we built may help clinicians accurately predict the survival of patients with COAD. Our findings provide valuable insights into the role of ERSRGs in COAD and may provide new targets for COAD therapy.

SUBMITTER: Yang X

PROVIDER: S-EPMC9740360 | biostudies-literature | 2022 Nov

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

System analysis based on the ER stress-related genes identifies WFS1 as a novel therapy target for colon cancer.

Yang Xianguang X Zhang Chaoyang C Yan Cheng C Ma Liukai L Ma Jiahao J Meng Xiaoke X

Aging 20221128 22

<h4>Background</h4>Colon cancer (COAD) is the third-largest common malignant tumor and the fourth major cause of cancer death in the world. Endoplasmic reticulum (ER) stress has a great influence on cell growth, migration, proliferation, invasion, angiogenesis, and chemoresistance of massive tumors. Although ER stress is known to play an important role in various types of cancer, the prognostic model based on ER stress-related genes (ERSRGs) in colon cancer has not been constructed yet. In this ...[more]

PMID: 36445321

Similar Datasets

Project description:BackgroundPancreatic adenocarcinoma (PAAD) is one of the most common malignant tumors of the digestive tract. Pyroptosis is a newly discovered programmed cell death that highly correlated with the prognosis of tumors. However, the prognostic value of pyroptosis in PAAD remains unclear.MethodsA total of 178 pancreatic cancer PAAD samples and 167 normal samples were obtained from The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) databases. The "DESeq2" R package was used to identify differntially expressed pyroptosis-related genes between normal pancreatic samples and PAAD samples. The prognostic model was established in TCGA cohort based on univariate Cox and the least absolute shrinkage and selection operator (LASSO) Cox regression analyses, which was validated in test set from Gene Expression Omnibus (GEO) cohort. Univariate independent prognostic analysis and multivariate independent prognostic analysis were used to determine whether the risk score can be used as an independent prognostic factor to predict the clinicopathological features of PAAD patients. A nomogram was used to predict the survival probability of PAAD patients, which could help in clinical decision-making. The R package "pRRophetic" was applied to calculate the drug sensitivity of each samples from high- and low-risk group. Tumor immune infiltration was investigated using an ESTIMATE algorithm. Finally, the pro-tumor phenotype of GSDMC was explored in PANC-1 and CFPAC-1 cells.ResultOn the basis of univariate Cox and LASSO regression analyses, we constructed a risk model with identified five pyroptosis-related genes (IL18, CASP4, NLRP1, GSDMC, and NLRP2), which was validated in the test set. The PAAD samples were divided into high-risk and low-risk groups on the basis of the risk score's median. According to Kaplan Meier curve analysis, samples from high-risk groups had worse outcomes than those from low-risk groups. The time-dependent receiver operating characteristics (ROC) analysis revealed that the risk model could predict the prognosis of PAAD accurately. A nomogram accompanied by calibration curves was presented for predicting 1-, 2-, and 3-year survival in PAAD patients. More importantly, 4 small molecular compounds (A.443654, PD.173074, Epothilone. B, Lapatinib) were identified, which might be potential drugs for the treatment of PAAD patients. Finally, the depletion of GSDMC inhibits the proliferation, invasion, and migration of pancreatic adenocarcinoma cells.ConclusionIn this study, we developed a pyroptosis-related prognostic model based on IL18, CASP4, NLRP1, NLRP2, and GSDMC , which may be helpful for clinicians to make clinical decisions for PAAD patients and provide valuable insights for individualized treatment. Our result suggest that GSDMC may promote the proliferation and migration of PAAD cell lines. These findings may provide new insights into the roles of pyroptosis-related genes in PAAD, and offer new therapeutic targets for the treatment of PAAD.

Project description:BackgroundLiver hepatocellular carcinoma (LIHC) ranks sixth among the most common types of cancer with a high mortality rate. Cuproptosis is a newly discovered type of cell death in tumor, which is characterized by accumulation of intracellular copper leading to the aggregation of mitochondrial lipoproteins and destabilization of proteins. Thus, understanding the exact effects of cuproptosis-related genes in LIHC and determining their prognosticvalue is critical. However, the prognostic model of LIHC based on cuproptosis-related genes has not been reported.MethodsFirstly, we downloaded transcriptome data and clinical information of LIHC patients from TCGA and GEO (GSE76427), respectively. We then extracted the expression of cuproptosis-related genes and established a prognostic model by lasso cox regression analysis. Afterwards, the prediction performance of the model was evaluated by Kaplan-Meier survival analysis and receiver operating characteristic curve (ROC). Then, the prognostic model and the expression levels of the three genes were validated using the dataset from GEO. Subsequently, we divided LIHC patients into two subtypes by non-negative matrix factorization (NMF) classification and performed survival analysis. We constructed a Sankey plot linking different subtypes and prognostic models. Next, we calculate the drug sensitivity of each sample from patients in the high-risk group and low-risk group by the R package pRRophetic. Finally, we verified the function of LIPT1 in LIHC.ResultsUsing lasso cox regression analysis, we developed a prognostic risk model based on three cuproptosis-related genes (GCSH, LIPT1 and CDKN2A). Both in the training and in the test sets, the overall survival (OS) of LIHC patients in the low-risk group was significantly longer than that in the high-risk group. By performing NMF cluster, we identified two molecular subtypes of LIHC (C1 and C2), with C1 subtype having significantly longer OS and PFS than C2 subtype. The ROC analysis indicated that our model had a precisely predictive capacity for patients with LIHC. The multivariate Cox regression analysis indicated that the risk score is an independent predictor. Subsequently, we identified 71 compounds with IC50 values that differed between the high-risk and low-risk groups. Finally, we determined that knockdown of LIPT1 gene expression inhibited proliferation and invasion of hepatoma cells.ConclusionIn this study, we developed a novel prognostic model for hepatocellular carcinoma based on cuproptosis-related genes that can effectively predict the prognosis of LIHC patients. The model may be helpful for clinicians to make clinical decisions for patients with LIHC and provide valuable insights for individualized treatment. Two distinct subtypes of LIHC were identified based on cuproptosis-related genes, with different prognosis and immune characteristics. In addition, we verified that LIPT1 may promote proliferation, invasion and migration of LIHC cells. LIPT1 might be a new potential target for therapy of LIHC.

Project description:The identification and characterization of tumor suppressor genes has enhanced our understanding of the biology of cancer and enabled the development of new diagnostic and therapeutic modalities. Whereas in past decades, a handful of tumor suppressors have been slowly identified using techniques such as linkage analysis, large-scale sequencing of the cancer genome has enabled the rapid identification of a large number of genes that are mutated in cancer. However, determining which of these many genes play key roles in cancer development has proven challenging. Specifically, recent sequencing of human breast and colon cancers has revealed a large number of somatic gene mutations, but virtually all are heterozygous, occur at low frequency, and are tumor-type specific. Keywords: Microarray, Hypermethylome, DNA-hypermethylation, DAC, TSA, Epigenetic, Mutation, Colorectal Cancer, Breast Cancer The candidate hypermethylome genes were identified as described below. These genes were compared to the CAN genes identified by Sjöblom T et al (Science, 2006). Cell culture and treatment. For drug treatments, log phase MCF7, T-47D, MDA-MB231 and MDA-MB-468 cells were cultured in McCoys 5A media (Invitrogen) containing 10% BCS and 1x penicillin/streptomycin with 5M 5aza-deoxycytidine (DAC) (Sigma; stock solution: 1mM in PBS) for 96 hours, replacing media and DAC every 24 hours. Cell treatment with 300nM Trichostatin A (Sigma; stock solution: 1.5mM dissolved in ethanol) was performed for 18 hours. Control cells underwent mock treatment in parallel with addition of equal volume of PBS or ethanol without drugs. Microarray analysis. Total RNA was harvested from log phase cells using the Qiagen kit according to the manufacturers instructions, including a DNAase step. RNA was quantified using the NanoDrop ND-100 followed by quality assessment with 2100 Bioanalyzer (Agilent Technologies). RNA concentrations for individual samples were greater than 200ng/ul, with 28s/18s ratios greater than 2.2 and RNA integrity numbers of 10 (highest). Sample amplification and labeling procedures were carried out using the Low RNA Input Fluorescent Linear Amplification Kit (Agilent Technologies) according to the manufacturers instructions. The labeled cRNA was purified using the RNeasy mini kit (Qiagen) and quantified. RNA spike-in controls (Agilent Technologies) were added to RNA samples before amplification. 0.75 microgram of samples labeled with Cy3 or Cy5 were mixed with control targets (Agilent Technologies), assembled on Oligo Microarray, hybridized, and processed according to the Agilent microarray protocol. Scanning was performed with the Agilent G2565BA microarray scanner under default settings recommended by Agilent Technologies. Data analysis. All arrays were subject to quality checks recommended by the manufacturer. Images were visually inspected for artifacts and distributions of signal and background intensity of both red and green channels were examined to identify anomalous arrays. No irregularities were observed, and all arrays were retained and used. All calculations were performed using the R statistical computing platform and packages from Bioconductor bioinformatics software project. The log ratio of red signal to green signal was calculated and Loess normalization as implemented in the limma package from Bioconductor.

Dataset Information

System analysis based on the ER stress-related genes identifies WFS1 as a novel therapy target for colon cancer.

Background

Methods

Results

Conclusions

Publications

System analysis based on the ER stress-related genes identifies WFS1 as a novel therapy target for colon cancer.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets