Dataset Information

Identification and Validation of a Novel RNA-Binding Protein-Related Gene-Based Prognostic Model for Multiple Myeloma.

ABSTRACT:

Background

Multiple myeloma (MM) is a malignant hematopoietic disease that is usually incurable. RNA-binding proteins (RBPs) are involved in the development of many tumors, but their prognostic significance has not been systematically described in MM. Here, we developed a prognostic signature based on eight RBP-related genes to distinguish MM cohorts with different prognoses.

Method

After screening the differentially expressed RBPs, univariate Cox regression was performed to evaluate the prognostic relevance of each gene using The Cancer Genome Atlas (TCGA)-Multiple Myeloma Research Foundation (MMRF) dataset. Lasso and stepwise Cox regressions were used to establish a risk prediction model through the training set, and they were validated in three Gene Expression Omnibus (GEO) datasets. We developed a signature based on eight RBP-related genes, which could classify MM patients into high- and low-score groups. The predictive ability was evaluated using bioinformatics methods. Gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment, and gene set enrichment analyses were performed to identify potentially significant biological processes (BPs) in MM.

Result

The prognostic signature performed well in the TCGA-MMRF dataset. The signature includes eight hub genes: HNRNPC, RPLP2, SNRPB, EXOSC8, RARS2, MRPS31, ZC3H6, and DROSHA. Kaplan-Meier survival curves showed that the prognosis of the risk status showed significant differences. A nomogram was constructed with age; B2M, LDH, and ALB levels; and risk status as prognostic parameters. Receiver operating characteristic (ROC) curve, C-index, calibration analysis, and decision curve analysis (DCA) showed that the risk module and nomogram performed well in 1, 3, 5, and 7-year overall survival (OS). Functional analysis suggested that the spliceosome pathway may be a major pathway by which RBPs are involved in myeloma development. Moreover, our signature can improve on the R-International Staging System (ISS)/ISS scoring system (especially for stage II), which may have guiding significance for the future.

Conclusion

We constructed and verified the 8-RBP signature, which can effectively predict the prognosis of myeloma patients, and suggested that RBPs are promising biomarkers for MM.

SUBMITTER: Wang W

PROVIDER: S-EPMC8107400 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BackgroundMultiple myeloma (MM) is an incurable and relapse-prone disease with apparently prognostic heterogeneity. At present, the risk stratification of myeloma is still incomplete. Pyroptosis, a type of programmed cell death, has been shown to regulate tumor growth and may have potential prognostic value. However, the role of pyroptosis-related genes (PRGs) in MM remains undetermined. The aims of this study were to identify potential prognostic biomarkers and to construct a predictive model related to PRGs.MethodsSequencing and clinical data were obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. Non-negative matrix factorization (NMF) was performed to identify molecular subtype screening. LASSO regression was used to screen for prognostic markers, and then a risk score model was constructed. The Maxstat package was utilized to calculate the optimal cutoff value, according to which patients were divided into a high-risk group and a low-risk group, and the survival curves were plotted using the Kaplan-Meier (K-M) method. Nomograms and calibration curves were established using the rms package.ResultsA total of 33 PRGs were extracted from the TCGA database underlying which 4 MM molecular subtypes were defined. Patients in cluster 1 had poorer survival than those in cluster 2 (p = 0.035). A total of 9 PRGs were screened out as prognostic markers, and the predictive ability of the 9-gene risk score for 3-year survival was best (AUC = 0.658). Patients in the high-risk group had worse survival than those in the low-risk group (p < 0.001), which was consistent with the results verified by the GSE2658 dataset. The nomogram constructed by gender, age, International Staging System (ISS) stage, and risk score had the best prognostic predictive performance with a c-index of 0.721.ConclusionOur model could enhance the predictive ability of ISS staging and give a reference for clinical decision-making. The new, prognostic, and pyroptosis-related markers screened out by us may facilitate the development of novel risk stratification for MM.Clinical trial registrationNot applicable.

Project description:Background: Ferroptosis is a newly discovered form of regulated cell death with distinct properties and recognizing functions involved in physical conditions or various diseases, including cancers. However, the relationship between gliomas and ferroptosis-related lncRNAs (FRLs) remains unclear. Methods: We collected a total of 1850 samples from The Cancer Genome Atlas (TCGA) and Genotype Tissue Expression (GTEX) databases, including 698 tumor and 1,152 normal samples. A list of ferroptosis-related genes was downloaded from the Ferrdb website. Differentially expressed FRLs (DEFRLS) were analyzed using the “limma” package in R software. Subsequently, prognosis-related FRLs were obtained by univariate Cox analysis. Finally, a prognostic model based on the 3 FRLs was constructed using Cox regression analysis with the least absolute shrinkage and selection operator (LASSO) algorithm. The prognostic power of the model was assessed using receiver operating characteristic (ROC) curve analysis and Kaplan-Meier (K-M) survival curve analysis. In addition, we further explored the relationship of the immune landscape and somatic mutations to prognostic model characteristics. Finally, we validated the function of LINC01426 in vitro. Results: We successfully constructed a 3-FRLs signature and classified glioma patients into high-risk and low-risk groups based on the risk score calculated from this signature. Compared with traditional clinicopathological features [age, sex, grade, isocitrate dehydrogenase (IDH) status], the prognostic accuracy of this model is more stable and stronger. Additionally, the model had stable predictive power for overall survival over a 5-year period. In addition, we found significant differences between the two groups in cellular immunity, the numbers of many immune cells, including NK cells, CD4+, CD8+ T-cells, and macrophages, and the expression of many immune-related genes. Finally, the two groups were also significantly different at the level of somatic mutations, especially in glioma prognosis-related genes such as IDH1 and ATRX, with lower mutation rates in the high-risk group leading to poorer prognosis. Finally, we found that the ferroptosis process of glioma cells was inhibited after knocking down the expression of LINC01426. Conclusion: The proposed 3-FRL signature is a promising biomarker for predicting prognostic features in glioma patients.

Project description:BackgroundColorectal cancer (CRC) is one of the most commonly diagnosed cancers in the world. Mitophagy is associated with tumorigenesis and development of malignancy. However, the specific role of mitophagy has yet not been systematically explored in CRC.MethodsThe RNA-sequencing dataset of CRC from The Cancer Genome Atlas (TCGA) and microarray data of gene expression profiles of CRC from Gene Expression Omnibus (GEO) were downloaded. Mitophagy-related gene sets were obtained from the Pathway Unification database. The package "limma" was used for differential gene expression analysis. Kaplan-Meier (KM) survival analyses were utilized to evaluate the prognostic value of the mitophagy regulators. Single-sample gene set enrichment analysis (ssGSEA) was used to estimate the infiltrating immune cells and the activity of immune response. The ConsensusClusterPlus algorithm was used to determine mitophagy-related subtypes. Principal component analysis (PCA) was used to create composite measurement of mitophagy scores. The R packages "survminer" and "ReGlot" were used to plot the nomogram and calibration curves.ResultsIntegrated analysis of the GEO and TCGA databases revealed some common differentially expressed genes (DEGs) in CRC. MFN2, UBB, PINK1, and PRKN were significantly downregulated in CRC samples as compared to normal samples, and other genes were significantly upregulated in CRC samples. KM survival analyses showed that high expression of ATG12 and MAP1LC3B predicted a poor prognosis, whereas high expression of TOMM22 and TOMM40 predicted a better prognosis. Mitophagy showed significant correlation with immune-related pathways in CRC samples. We identified 2 distinct CRC subtypes with different mitophagy accumulation, of which subtype B had better prognosis and immune activity. The mitophagy score may be employed as a new and efficient clinical predictor in conjunction with other clinical indicators to predict the prognosis of CRC patients.ConclusionsWe systematically investigated the CRC heterogeneity with reference to mitophagy based on bioinformatics analyses, and the findings of this study might provide some guidance for future research into potential biomarkers for diagnosis and prognosis prediction of CRC patients.

Project description:BackgroundThe role of N7-methyladenosine (m7G)-related genes in the progression and prognosis of gastric cancer (GC) remains unclear. This study aimed to explore prognostic biomarkers for GC based on m7G methylation regulators and to construct a prognostic risk model.MethodsRNA sequencing profiles with corresponding clinicopathological information associated with GC of which the histological type was stomach adenocarcinoma (STAD) were obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO), respectively. A total of 29 m7G regulators were extracted from previous studies. According to the expression similarity of m7G regulators, the GC samples obtained from TCGA were further classified into 2 clusters demonstrating different overall survival (OS) rates and genetic heterogeneity, and the differentially expressed genes (DEGs) between these 2 clusters were defined as m7G-related genes. Univariate regression analysis and regression analysis were then used to obtain the prognostic m7G-related genes. The samples in TCGA and Genotype-Tissue Expression (GTEx) were used to verify the differential expression and prognostic value of these m7G-related genes contained in the prognostic model. Subsequently, the risk score was combined with other prognostic factors to develop a nomogram. The predictive ability of the nomogram was evaluated by the standard receiver operating characteristic (ROC) curve. Gene set enrichment analysis (GSEA) was used to identify activation pathways in both groups. Finally, the association between the prognostic model and the immune characteristics of GC were appraised.ResultsA prognostic model consisting of 11 m7G-related genes was constructed. GC patients in the high-risk group were shown to have a poor prognosis and this result was further demonstrated in each group. The risk model can be applied for patients with different clinical features. The results of GSEA showed that cell adhesion, cell junction, and focal adhesion were highly enriched in the high-risk group. In addition, we found that the expression of programmed cell death ligand 1 (PD-L1) was significantly elevated in the low-risk group, whereas programmed cell death ligand 2 (PD-L2) and tumor necrosis factor receptor superfamily member 4 (TNFRSF4) were overexpressed in the high-risk group.ConclusionsWe successfully built and verified a m7G relevant prognostic model for predicting prognosis and providing a new train of thought for improving the treatment of GC.

Project description:Previous literatures have suggested the importance of inflammatory response during lung adenocarcinoma (LUAD) development. This study aimed at exploring the inflammation-related genes and developing a prognostic signature for predicting the prognosis of LUAD. Survival‑associated inflammation-related genes were identified by univariate Cox regression analysis in the dataset of The Cancer Genome Atlas (TCGA). The least absolute shrinkage and selection operator (LASSO) penalized Cox regression model was used to derive a risk signature which is significantly negatively correlated with OS and divide samples into high-, medium- and low-risk group. Univariate and multivariate Cox analyses suggested that the level of risk group was an independent prognostic factor of the overall survival (OS). Time-dependent receiver operating characteristic (ROC) curve indicated the AUC of 1-, 3- and 5-years of the risk signature was 0.715, 0.719, 0.699 respectively. A prognostic nomogram was constructed by integrating risk group and clinical features. The independent dataset GSE30219 of Gene Expression Omnibus (GEO) was used for verification. We further explored the differences among risk groups in Gene set enrichment analysis (GSEA), tumor mutation and tumor microenvironment. Furthermore, Single Sample Gene Set Enrichment Analysis (ssGSEA) and the results of Cell-type Identification By Estimating Relative Subsets Of RNA Transcripts (CIBERSORT) suggested the status of immune cell infiltration was highly associated with risk groups. We demonstrated the prediction effect of CTLA-4 and PD-1/PD-L1 inhibitors in the low-risk group was better than that in the high-risk group using two methods of immune score include immunophenoscore from The Cancer Immunome Atlas (TCIA) and TIDE score from Tumor Immune Dysfunction and Exclusion (TIDE). In addition, partial targeted drugs and chemotherapy drugs for lung cancer had higher drug sensitivity in the high-risk group. Our findings provide a foundation for future research targeting inflammation-related genes to predictive prognosis and some reference significance for the selection of immunotherapy and drug regimen for lung adenocarcinoma.

Project description:BackgroundCuproptosis is a novel form of copper-induced cell death that targets lipoylated tricarboxylic acid (TCA) cycle proteins. However, its prognostic role in lung adenocarcinoma (LUAD) remains unclear. This study aimed to establish a cuproptosis-related prognostic signature for patients with LUAD.MethodsTranscriptome data of LUAD samples were extracted from the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. The prognostic value of cuproptosis-related genes (CRGs) was investigated using Cox regression analysis to develop a cuproptosis-related prognostic model. Kyoto Encyclopedia of Genes and Genomes (KEGG), gene ontology (GO) and gene set variation analysis (GSVA) were conducted to characterize different biological activities or pathways between high- or low-CRG groups. The expression pattern and prognostic values of CRGs were validated in 37 paired tumor-normal samples using quantitative PCR. Furthermore, in vitro experiments were performed to investigate the relationship between cuproptosis and CRG expression and to explore the function of target genes in cuproptosis.ResultsAmong the 36 CRGs, 17 genes were upregulated, and 3 genes were downregulated in LUAD. A total of 385 CRGs were identified using Pearson correlation analysis. A cuproptosis-related signature was constructed using least absolute shrinkage and selection operator (LASSO) analysis. The prognostic value of the cuproptosis-related signature was validated in six external validation cohorts and in LUAD specimens from our facility. Patients in the high-risk group based on the CRG signature score had shorter overall survival than those in the low-risk group in both the datasets and clinical specimens. In vitro experiments revealed that the expression of BARX1, GFRA3, and KHDRBS2 was upregulated after cuproptosis was induced by elesclomol-CuCL2, whereas the upregulation was suppressed on pretreatment with tetrathiomolybdate (TTM), a chelator of copper. Further, the cell proliferation assay revealed that the BARX1 and GFRA3 deficiency facilities the cuproptosis induced by elesclomol-CuCL2.ConclusionThis study established a new CRG signature that can be used to predict the OS of LUAD patients. Moreover, the knockdown of BARX1 and GFRA3 could increase the sensitivity of LUAD cells to the cuproptosis.

Project description:Dysregulated long non-coding RNAs (lncRNAs) have been found to have oncogenic and/or tumor suppressive roles in the development and progression of cancer, implying their potentials as novel independent biomarkers for cancer diagnosis and prognosis. However, the prognostic significance of expression profile-based lncRNA signature for outcome prediction in patients with multiple myeloma (MM) has not yet been investigated.LncRNA expression profiles of a large cohort of patients with MM were obtained and analyzed by repurposing the publically available microarray data. An lncRNA-focus risk score model was developed from the training dataset, and then validated in the testing and another two independent external datasets. The time-dependent receiver operating characteristic (ROC) curve was used to evaluate the prognostic performance for survival prediction. The biological function of prognostic lncRNAs was predicted using bioinformatics analysis.Four lncRNAs were identified to be significantly associated with overall survival (OS) of patients with MM in the training dataset, and were combined to develop a four-lncRNA prognostic signature to stratify patients into high-risk and low-risk groups. Patients of training dataset in the high-risk group exhibited shorter OS than those in the low-risk group (HR?=?2.718, 95 % CI?=?1.937-3.815, p <0.001). The similar prognostic values of four-lncRNA signature were observed in the testing dataset, entire GSE24080 dataset and another two independent external datasets. Multivariate Cox regression and stratified analysis showed that the prognostic power of four-lncRNA signature was independent of clinical features, including serum beta 2-microglobulin (S?2M), serum albumin (ALB) and lactate dehydrogenase (LDH). ROC analysis also demonstrated the better performance for predicting 3-year OS. Functional enrichment analysis suggested that these four lncRNAs may be involved in known genetic and epigenetic events linked to MM.Our results demonstrated potential application of lncRNAs as novel independent biomarkers for diagnosis and prognosis in MM. These lncRNA biomarkers may contribute to the understanding of underlying molecular basis of MM.

Project description:BackgroundMultiple myeloma (MM), the second most hematological malignancy, has high incidence and remains incurable till now. The pathogenesis of MM is poorly understood. This study aimed to identify novel prognostic model for MM on gene expression profiles.MethodsGene expression datas of MM (GSE6477, GSE136337) were downloaded from Gene Expression Omnibus (GEO) database. The differentially expressed genes (DEGs) in GSE6477 between case samples and normal control samples were screened by the limma package. Meanwhile, enrichment analysis was conducted, and a protein-protein interaction (PPI) network of these DEGs was established by STRING and cytoscape software. Co-expression modules of genes were built by Weighted Correlation Network Analysis (WGCNA). Key genes were identified both from hub genes and the DEGs. Univariate and multivariate Cox congression were performed to screen independent prognostic genes to construct a predictive model. The predictive power of the model was evaluated by Kaplan-Meier curve and time-dependent receiver operating characteristic (ROC) curves. Finally, univariate and multivariate Cox regression analyse were used to investigate whether the prognostic model could be independent of other clinical parameters.ResultsGSE6477, including 101 case and 15 normal control, were screened as the datasets. A total of 178 DEGs were identified, including 59 up-regulated and 119 down-regulated genes. In WGCNA analysis, module black and module purple were the most relevant modules with cancer traits, and 92 hub genes in these two modules were selected for further analysis. Next, 47 genes were chosen both from the DEGs and hub genes as key genes. Three genes (LYVE1, RNASE1, and RNASE2) were finally screened by univariate and multivariate Cox regression analyses and used to construct a risk model. In addition, the three-gene prognostic model revealed independent and accurate prognostic capacity in relation to other clinical parameters for MM patients.ConclusionIn summary, we identified and constructed a three-gene-based prognostic model that could be used to predict overall survival of MM patients.