Project description:Ovarian cancer (OC) is a female reproductive system tumor. RNA modifications play key roles in gene expression regulation. The growing evidence demonstrates that RNA methylation is critical for various biological functions, and that its dysregulation is related to the progression of cancer in human. OC samples were classified into different subtypes (Clusters 1 and 2) based on various RNA-modification regulatory genes (RRGs) in the process of RNA modifications (m1A, m6A, m6Am, m5C, m7G, ac4C, m3C, and Ψ) by nonnegative matrix factorization method (NMF). Based on differently expressed RRGs (DERRGs) between clusters, a pathologically specific RNA-modification regulatory gene signature was constructed with Lasso regression. Kaplan-Meier analysis and receiver operating characteristic (ROC) curves were used to evaluate the prognostic ability of the identified model. The correlations of clinicopathological features, immune subtypes, immune scores, immune cells, and tumor mutation burden (TMB) were also estimated between different NMF clusters and riskscore groups. In this study, 59 RRGs in the process of RNA modifications (m1A, m6A, m6Am, m5C, m7G, ac4C, m3C, and Ψ) were obtained from TCGA database. These RRGs were interactional, and sample clusters based on these regulators were significantly correlated with survival rate, clinical characteristics (involving survival status and pathologic stage), drug sensibility, and immune microenvironment. Furthermore, Lasso regression based on these 21 DERRGs between clusters 1 and 2 constructed a four-DERRG signature (ALYREF, ZC3H13, WTAP, and METTL1). Based on this signature, 307 OC patients were classified into high- and low-risk groups based on median value of riskscores from lasso regression. This identified signature was significantly associated with overall survival, radiation therapy, age, clinical stage, cancer status, and immune cells (involving CD4+ memory resting T cells, plasma cells, and Macrophages M1) of ovarian cancer patients. Further, GSEA revealed that multiple biological behaviors were significantly enriched in different groups. OC patients were classified into two subtypes per these RRGs. This study identified four-DERRG signature (ALYREF, ZC3H13, WTAP, and METTL1) in OC, which was an independent prognostic model for patient stratification, prognostic evaluation, and prediction of response to immunotherapy in ovarian cancer by classifying OC patients into high- and low-risk groups.

Project description:BackgroundLung adenocarcinoma (LUAD) is a major subtype of lung cancer and closely associated with poor prognosis. N6-methyladenosine (m6A), one of the most predominant modifications in mRNAs, is found to participate in tumorigenesis. However, the potential function of m6A RNA methylation in the tumor immune microenvironment is still murky.MethodsThe gene expression profile cohort and its corresponding clinical data of LUAD patients were downloaded from TCGA database and GEO database. Based on the expression of 21 m6A regulators, we identified two distinct subgroups by consensus clustering. The single-sample gene-set enrichment analysis (ssGSEA) algorithm was conducted to quantify the relative abundance of the fraction of 28 immune cell types. The prognostic model was constructed by Lasso Cox regression. Survival analysis and receiver operating characteristic (ROC) curves were used to evaluate the prognostic model.ResultConsensus classification separated the patients into two clusters (clusters 1 and 2). Those patients in cluster 1 showed a better prognosis and were related to higher immune scores and more immune cell infiltration. Subsequently, 457 differentially expressed genes (DEGs) between the two clusters were identified, and then a seven-gene prognostic model was constricted. The survival analysis showed poor prognosis in patients with high-risk score. The ROC curve confirmed the predictive accuracy of this prognostic risk signature. Besides, further analysis indicated that there were significant differences between the high-risk and low-risk groups in stages, status, clustering subtypes, and immunoscore. Low-risk group was related to higher immune score, more immune cell infiltration, and lower clinical stages. Moreover, multivariate analysis revealed that this prognostic model might be a powerful prognostic predictor for LUAD. Ultimately, the efficacy of this prognostic model was successfully validated in several external cohorts (GSE30219, GSE50081 and GSE72094).ConclusionOur study provides a robust signature for predicting patients' prognosis, which might be helpful for therapeutic strategies discovery of LUAD.

Project description:BackgroundFerroptosis is a newly discovered form of cell death characterized by iron-dependent lipid peroxidation. This study aims to investigate the potential correlation between ferroptosis and the prognosis of lung adenocarcinoma (LUAD).MethodsRNA-seq data were collected from the LUAD dataset of The Cancer Genome Atlas (TCGA) database. Based on ferroptosis-related genes, differentially expressed genes (DEGs) between LUAD and paracancerous specimens were identified. The univariate Cox regression analysis was performed to screen key genes associated with the prognosis of LUAD. LUAD patients were divided into the training set and validation set. Then, we screened out key genes and built a prognostic prediction model involving 5 genes using the least absolute shrinkage and selection operator (LASSO) regression with tenfold cross-validation and the multivariate Cox regression analysis. After dividing LUAD patients based on the median level of risk score as cut-off value, the generated prognostic prediction model was validated in the validation set. Moreover, we analyzed the somatic mutations, and estimated the scores of immune infiltration in the high-risk and low-risk groups. Functional enrichment analysis of DEGs was performed as well.ResultsHigh-risk scores indicated the worse prognosis of LUAD. The maximum area under curve (AUC) of the training set and the validation set in this study was 0.7 and 0.69, respectively. Moreover, we integrated the age, gender, and tumor stage to construct the composite nomogram. The charts indicated that the AUC of LUAD cases with the survival time of 1, 3 and 5 years was 0.698, 0.71 and 0.73, respectively. In addition, the mutation frequency of LUAD patients in the high-risk group was significantly higher than that in the low-risk group. Simultaneously, DEGs were mainly enriched in ferroptosis-related pathways by analyzing the functional results.ConclusionsThis study constructs a novel LUAD prognosis prediction model involving 5 ferroptosis-related genes, which can be used as a promising tool for decision-making of clinical therapeutic strategies of LUAD.

Project description:The magnitude of the immune response is closely associated with clinical outcome in patients with cancer. However, finding potential therapeutic targets for lung cancer in the immune system remains challenging. Here, we constructed a vital immune-prognosis genes (VIPGs) based cluster of lung adenocarcinoma (LUAD) from IMMPORT databases and The Cancer Genome Atlas. A transcription factor regulatory network for the VIPGs was also established. The tumor microenvironment of LUAD was analyzed using the ESTIMATE (Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data) algorithm and single-sample Gene Set Enrichment Analysis. The immune checkpoints and genomic alterations were explored in the different immune clusters. We identified 15 VIPGs for patients with LUAD and clustered the patients into low-immunity and high-immunity subtypes. The immune score, stromal score and ESTIMATE score were significantly higher in the high-immunity subtype, whereas tumor purity was higher in the low-immunity subtype. In addition, the immune checkpoints cytotoxic T lymphocyte associate protein-4(CTLA4), programmed cell death protein-1 and programmed death-ligand were elevated in the low-immunity subtype. The genomic results also showed that the tumor mutation burden was higher in the high-immunity subtype. Finally, Gene Set Enrichment Analysis showed that several immune-related gene sets, including interleukin-2/STAT5 signaling, inflammatory response, interleukin-6/Janus kinase(JAK)/signal transducer and activator of transcription 3 (STAT3) signaling, interferon-gamma response and allograft rejection, were elevated in the high-immunity subtype. Finally, high-immunity patients exhibited greater overall and disease-specific survival outcome compared with low-immunity patients (log rank P = 0.013 and P = 0.0097). Altogether, here we have identified 15 immune-prognosis genes and a potential immune subtype for patients with LUAD, which may provide new insights into the prognosis and treatment of LUAD.

Project description:Evidence is increasingly indicating that circular RNAs (circRNAs) are closely involved in tumorigenesis and cancer progression. However, the function and application of circRNAs in lung adenocarcinoma (LUAD) are still unknown. In this study, we constructed a circRNA-associated competitive endogenous RNA (ceRNA) network to investigate the regulatory mechanism of LUAD procession and further constructed a prognostic signature to predict overall survival for LUAD patients. Differentially expressed circRNAs (DEcircRNAs), differentially expressed miRNAs (DEmiRNAs), and differentially expressed mRNAs (DEmRNAs) were selected to construct the ceRNA network. Based on the TargetScan prediction tool and Pearson correlation coefficient, we constructed a circRNA-associated ceRNA network including 11 DEcircRNAs, 8 DEmiRNAs, and 49 DEmRNAs. GO and KEGG enrichment indicated that the ceRNA network might be involved in the regulation of GTPase activity and endothelial cell differentiation. After removing the discrete points, a PPI network containing 12 DEmRNAs was constructed. Univariate Cox regression analysis showed that three DEmRNAs were significantly associated with overall survival. Therefore, we constructed a three-gene prognostic signature for LUAD patients using the LASSO method in the TCGA-LUAD training cohort. By applying the signature, patients could be categorized into the high-risk or low-risk subgroups with significant survival differences (HR: 1.62, 95% CI: 1.12-2.35, log-rank p = 0.009). The prognostic performance was confirmed in an independent GEO cohort (GSE42127, HR: 2.59, 95% CI: 1.32-5.10, log-rank p = 0.004). Multivariate Cox regression analysis proved that the three-gene signature was an independent prognostic factor. Combining the three-gene signature with clinical characters, a nomogram was constructed. The primary and external verification C-indexes were 0.717 and 0.716, respectively. The calibration curves for the probability of 3- and 5-year OS showed significant agreement between nomogram predictions and actual observations. Our findings provided a deeper understanding of the circRNA-associated ceRNA regulatory mechanism in LUAD pathogenesis and further constructed a useful prognostic signature to guide personalized treatment of LUAD patients.

Project description:Lung adenocarcinoma (LUAD) is one of the most widespread and fatal types of lung cancer. Oxidative stress, resulting from an imbalance in the production and accumulation of reactive oxygen species (ROS), is considered a promising therapeutic target for cancer treatment. Currently, immune checkpoint blockade (ICB) therapy is being explored as a potentially effective treatment for early-stage LUAD. In this research, we aim to identify distinct subtypes of LUAD patients by investigating genes associated with oxidative stress and immunotherapy. Additionally, we aim to propose subtype-specific therapeutic strategies. We conducted a thorough search of the Gene Expression Omnibus (GEO) datasets. From this search, we pinpointed datasets that contained both expression data and survival information. We selected genes associated with oxidative stress and immunotherapy using keyword searches on GeneCards. We then combined expression data of LUAD samples from both The Cancer Genome Atlas (TCGA) and 11 GEO datasets, forming a unified dataset. This dataset was subsequently divided into two subsets, Dataset_Training and Dataset_Testing, using a random bifurcation method, with each subset containing 50% of the data. We applied consensus clustering (CC) analysis to identify distinct LUAD subtypes within the Dataset_Training. Molecular variances associated with oxidative stress levels, the tumor microenvironment (TME), and immune checkpoint genes (ICGs) were then investigated among these subtypes. Employing feature selection combined with machine learning techniques, we constructed models that achieved the highest accuracy levels. We validated the identified subtypes and models from Dataset_Training using Dataset_Testing. A hub gene with the highest importance values in the machine learning model was identified. We then utilized virtual screening to discover potential compounds targeting this hub gene. In the unified dataset, we integrated 2,154 LUAD samples from TCGA-LUAD and 11 GEO datasets. We specifically selected 1,311 genes associated with immune and oxidative stress processes. The expression data of these genes were then employed for subtype identification through CC analysis. Within Dataset_Training, two distinct subtypes emerged, each marked by different levels of immune and oxidative stress pathway values. Consequently, we named these as the OX+ and IM+ subtypes. Notably, the OX+ subtype showed increased oxidative stress levels, correlating with a worse prognosis than the IM+ subtype. Conversely, the IM+ subtype demonstrated enhanced levels of immune pathways, immune cells, and ICGs compared to the OX+ subtype. We reconfirmed these findings in Dataset_Testing. Through gene selection, we identified an optimal combination of 12 genes for predicting LUAD subtypes: ACP1, AURKA, BIRC5, CYC1, GSTP1, HSPD1, HSPE1, MDH2, MRPL13, NDUFS1, SNRPD1, and SORD. Out of the four machine learning models we tested, the support vector machine (SVM) stood out, achieving the highest area under the curve (AUC) of 0.86 and an accuracy of 0.78 on Dataset_Testing. We focused on HSPE1, which was designated as the hub gene due to its paramount importance in the SVM model, and computed the docking structures for four compounds: ZINC3978005 (Dihydroergotamine), ZINC52955754 (Ergotamine), ZINC150588351 (Elbasvir), and ZINC242548690 (Digoxin). Our study identified two subtypes of LUAD patients based on oxidative stress and immunotherapy-related genes. Our findings provided subtype-specific therapeutic strategies.

Dataset Information

A novel nomogram model for lung adenocarcinoma subtypes based on RNA-modification regulatory genes

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets