Identification of the Signature Associated With m6A RNA Methylation Regulators and m6A-Related Genes and Construction of the Risk Score for Prognostication in Early-Stage Lung Adenocarcinoma

BackgroundN6-methyladenosine (m6A) RNA modification is vital for cancers because methylation can alter gene expression and even affect some functional modification. Our study aimed to analyze m6A RNA methylation regulators and m6A-related genes to understand the prognosis of early lung adenocarcinoma.MethodsThe relevant datasets were utilized to analyze 21 m6A RNA methylation regulators and 5,486 m6A-related genes in m6Avar. Univariate Cox regression analysis, random survival forest analysis, Kaplan–Meier analysis, Chi-square analysis, and multivariate cox analysis were carried out on the datasets, and a risk prognostic model based on three feature genes was constructed.ResultsRespectively, we treated GSE31210 (n = 226) as the training set, GSE50081 (n = 128) and TCGA data (n = 400) as the test set. By performing univariable cox regression analysis and random survival forest algorithm in the training group, 218 genes were significant and three prognosis-related genes (ZCRB1, ADH1C, and YTHDC2) were screened out, which could divide LUAD patients into low and high-risk group (P < 0.0001). The predictive efficacy of the model was confirmed in the test group GSE50081 (P = 0.0018) and the TCGA datasets (P = 0.014). Multivariable cox manifested that the three-gene signature was an independent risk factor in LUAD. Furthermore, genes in the signature were also externally validated using the online database. Moreover, YTHDC2 was the important gene in the risk score model and played a vital role in readers of m6A methylation.ConclusionThe findings of this study suggested that associated with m6A RNA methylation regulators and m6A-related genes, the three-gene signature was a reliable prognostic indicator for LUAD patients, indicating a clinical application prospect to serve as a potential therapeutic target.

Download Full-text

Identification of a ferroptosis-related gene signature (FRGS) for predicting clinical outcome in lung adenocarcinoma

PeerJ ◽

10.7717/peerj.11233 ◽

2021 ◽

Vol 9 ◽

pp. e11233

Author(s):

Sheng Wang ◽

Chunlei Wu ◽

Dehua Ma ◽

Quanteng Hu

Keyword(s):

Regression Analysis ◽

Lung Adenocarcinoma ◽

Risk Score ◽

Prognostic Value ◽

Cox Regression ◽

Gene Signature ◽

Tnm Stage ◽

Related Gene ◽

Cox Regression Analysis ◽

Prognostic Predictors

Background Lung adenocarcinoma (LUAD) is the most common pathological subtype of lung cancer. Ferroptosis, an oxidative, iron-dependent form of necrotic cell death, is highly associated with tumorigenesis and cancer progression. However, the prognostic value of ferroptosis progress in LUAD was still rarely be investigated. Methods Herein, we collected three mRNA expression profiles and 85 ferroptosis-related genes from public databases. The “limma” package was used to identify ferroptosis-related differentially expressed genes (DEGs). Univariate Cox regression analysis and LASSO regression analysis were applied to screen and develop a ferroptosis-related gene signature (FRGS) and a formula to calculate the risk score. Multivariate Cox regression analysis was implemented to determine independent prognostic predictors of overall survival (OS). The area under the receiver operating characteristic curve (AUC) and calibration plot were used to evaluate the predictive accuracy of the FRGS and nomogram. Results We developed a FRGS with five genes (CYBB, CISD1, FADD, SAT2, VDAC2). The AUC of the FRGS in TCGA cohort was 0.777 at 1-year, 0.721 at 3-year and 0.725 at 5-year, significantly superior to the AUC of TNM stage (1-year: 0.701, 3-year: 0.691, 5-year: 0.686). A similar phenomenon was observed in GEO cohort 1 and 2. Multivariate Cox regression analysis indicted TNM stage and risk score were independent prognostic predictors. Finally, we built a nomogram with TNM stage and FRGS, the AUCs of which markedly higher than that of FRGS or TNM stage alone. Conclusion We constructed a prognostic FRGS with five ferroptosis-related genes and a nomogram for predicting the 1-, 3- and 5-year survival rate of LUAD patients, which may provide a new understanding of the prognostic value of ferroptosis progress in LUAD and will benefit prognosis assessment of LUAD patients.

Download Full-text

Identification of the Signature Associated with m6A RNA Methylation Regulators and m6A-Related Genes and Construction of the Risk Score for Prognostication in Early-Stage Lung Adenocarcinoma

10.21203/rs.3.rs-131013/v1 ◽

2020 ◽

Author(s):

Bingzhou Guo ◽

Hongliang Zhang ◽

Jinliang Wang ◽

Rilige Wu ◽

Junyan Zhang ◽

...

Keyword(s):

Lung Adenocarcinoma ◽

Cox Regression ◽

Risk Group ◽

Early Stage ◽

Gene Signature ◽

Cox Regression Analysis ◽

Rna Methylation ◽

M6a Rna Methylation ◽

Cox Analysis ◽

Random Survival Forest

Abstract Background: N6-methyladenosine (m6A) RNA modification play critical roles in multiple levels including the regulation of degradation, subcellular localization, splicing and local conformation changes of RNA transcripts. In this study, we aim to conduct comprehensive investigation on m6A RNA methylation regulators and m6A-related genes and their association with prognosis in early-stage Lung adenocarcinoma (LUAD). Methods: The relevant datasets which were used to analyze 21 m6A RNA methylation regulators and 887 m6A-related genes in m6Avar were downloaded from Gene Expression Omnibus database (GEO) and The Cancer Genome Atlas (TCGA) databases. Univariate cox regression analysis, random survival forest analysis, Kaplan-Meier anylysis, STRING and multivariate cox analysis were carried out, and a risk prognostic model based on five feature genes was constructed. Respectively, we treated GSE31210 (n=226) as training set, GSE50081 (n=128) and TCGA data (n=461) as test set. Results: By performing univariable cox regression and random survival forest algorithm in the training group, five prognosis-related genes including DENND1A, KBTBD6, KIF4A, BMPER, and YTHDC2 were screened out, which could divide LUAD patients into low-risk group and high-risk group (log rank p <0.001). The predictive efficacy of these genes was confirmed in the test group GSE50081 (log rank p <0.01) and TCGA datasets (log rank p <0.001). Cox analysis showed that this five-gene signature was an independent risk factor in LUAD. Further, genes in the signature were also external validated using online database. YTHDC2 played vital role of readers in m6A methylation. Conclusions: The findings of this study suggested that associated with m6A-related genes and m6A RNA methylation regulators, five-gene signature was reliable prognostic indicator for LUAD patients, indicating a clinical application prospect to serve as a potential therapeutic target.

Download Full-text

Development of an immune-related gene pairs signature for predicting clinical outcome in lung adenocarcinoma

Scientific Reports ◽

10.1038/s41598-021-83120-4 ◽

2021 ◽

Vol 11 (1) ◽

Author(s):

Chunlei Wu ◽

Quanteng Hu ◽

Dehua Ma

Keyword(s):

Overall Survival ◽

Regression Analysis ◽

Lung Adenocarcinoma ◽

Risk Score ◽

Cox Regression ◽

Related Gene ◽

Cox Regression Analysis ◽

Gene Pairs ◽

Immune Related Gene ◽

Pathological Subtype

AbstractLung adenocarcinoma (LUAD) is the main pathological subtype of Non-small cell lung cancer. We downloaded the gene expression profile and immune-related gene set from the TCGA and ImmPort database, respectively, to establish immune-related gene pairs (IRGPs). Then, IRGPs were subjected to univariate Cox regression analysis, LASSO regression analysis, and multivariable Cox regression analysis to screen and develop an IRGPs signature. The receiver operating characteristic curve (ROC) was applied for evaluating the predicting accuracy of this signature by calculating the area under ROC (AUC) and data from the GEO set was used to validate this signature. The relationship of 22 tumor-infiltrating immune cells (TIICs) to the immune risk score was also investigated. An IRGPs signature with 8 IRGPs was constructed. The AUC for 1- and 3-year overall survival in the TCGA set was 0.867 and 0.870, respectively. Similar results were observed in the AUCs of GEO set 1, 2 and 3 (GEO set 1 [1-year: 0.819; 3-year: 0.803]; GEO set 2 [1-year: 0.834; 3-year: 0.870]; GEO set 3 [1-year: 0.955; 3-year: 0.827]). Survival analysis demonstrated high-risk LUAD patients exhibited poorer prognosis. The multivariable Cox regression indicated that the risk score was an independent prognostic factor. The immune risk score was highly associated with several TIICs (Plasma cells, memory B cells, resting memory CD4 T cells, and activated NK cells). We developed a novel IRGPs signature for predicting 1- and 3- year overall survival in LUAD, which would be helpful for prognosis assessment of LUAD.

Download Full-text

A Six CT83-Related Genes Based Prognostic Signature for Lung Adenocarcinoma

Combinatorial Chemistry & High Throughput Screening ◽

10.2174/1871520621666210713112630 ◽

2021 ◽

Vol 24 ◽

Author(s):

Yongmei Wang ◽

Guimin Zhang ◽

Ruixian Wang

Keyword(s):

Regression Analysis ◽

Lung Adenocarcinoma ◽

Differential Expression ◽

Risk Score ◽

Prognostic Model ◽

Cox Regression ◽

Differential Expression Analysis ◽

Functional Enrichment ◽

Risk Scores ◽

Cox Regression Analysis

Background: This study aims to explore the prognostic values of CT83 and CT83-related genes in lung adenocarcinoma (LUAD). Methods: We downloaded the mRNA profiles of 513 LUAD patients (RNA sequencing data) and 246 NSCLC patients (Affymetrix Human Genome U133 Plus 2.0 Array) from TCGA and GEO databases. According to the median expression of CT83, the TCGA samples were divided into high and low expression groups, and differential expression analysis between them was performed. Functional enrichment analysis of differential expression genes (DEGs) was conducted. Univariate Cox regression analysis and LASSO Cox regression analysis were performed to screen the optimal prognostic DEGs. Then we established the prognostic model. A Nomogram model was constructed to predict the overall survival (OS) probability of LUAD patients. Results: CT83 expression was significantly correlated to the prognosis of LUAD patients. A total of 59 DEGs were identified, and a predictive model was constructed based on six optimal CT83-related DEGs, including CPS1, RHOV, TNNT1, FAM83A, IGF2BP1, and GRIN2A, could effectively predict the prognosis of LUAD patients. The nomogram could reliably predict the OS of LUAD patients. Moreover, the six important immune checkpoints (CTLA4, PD1, IDO1, TDO2, LAG3, and TIGIT) were closely correlated with the Risk Score, which was also differentially expressed between the LUAD samples with high and low-Risk Scores, suggesting that the poor prognosis of LUAD patients with high-Risk Score might be due to the immunosuppressive microenvironments. Conclusion: A prognostic model based on six optimal CT83 related genes could effectively predict the prognosis of LUAD patients.

Download Full-text

Prognostic Nomogram Based on Circular RNA-Associated Competing Endogenous RNA Network for Patients with Lung Adenocarcinoma

Oxidative Medicine and Cellular Longevity ◽

10.1155/2021/9978206 ◽

2021 ◽

Vol 2021 ◽

pp. 1-13

Author(s):

Yang Li ◽

Rongrong Sun ◽

Rui Li ◽

Yonggang Chen ◽

He Du

Keyword(s):

Overall Survival ◽

Regression Analysis ◽

Lung Adenocarcinoma ◽

Regulatory Mechanism ◽

Cox Regression ◽

Gene Signature ◽

Differentially Expressed ◽

Prognostic Signature ◽

Cox Regression Analysis ◽

Cerna Network

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.

Download Full-text

A Novel Ferroptosis-Related Genes Model for Prognosis Prediction of Lung Adenocarcinoma

10.21203/rs.3.rs-354347/v1 ◽

2021 ◽

Author(s):

Fei Li ◽

Dongcen Ge ◽

Shu-lan Sun

Keyword(s):

Regression Analysis ◽

High Risk ◽

Prediction Model ◽

Lung Adenocarcinoma ◽

Risk Score ◽

Cox Regression ◽

Risk Group ◽

Low Risk ◽

Cox Regression Analysis ◽

Validation Set

Abstract Background. Ferroptosis is a newly discovered form of cell death characterized by iron-dependent lipid peroxidation. The aim of this study is to investigate the relationship between ferroptosis and the prognosis of lung adenocarcinoma (LUAD).Methods. RNA-seq data was collected from the LUAD dataset of The Cancer Genome Altas (TCGA) database. We used ferroptosis-related genes as the basis, and identify the differential expression genes (DEGs) between cancer and paracancer. The univariate Cox regression analysis were used to screen the prognostic-related genes. We divided the patients into training and validation sets. Then, we screened out key genes and built a 5 genes prognostic prediction model by the applications of the least absolute shrinkage and selection operator (LASSO) 10-fold cross-validation and the multi-variate Cox regression analysis. We divided the cases by the median value of risk score and validated this model in the validation set. Meanwhile, we analyzed the somatic mutations, and estimated the score of immune infiltration in the high- and low-risk groups, as well as performed functional enrichment analysis of DEGs.Results. The result revealed that the high-risk score triggered the worse prognosis. The maximum area under curve (AUC) of the training set and the validation set of in this study was 0.7 and 0.69. Moreover, we integrated the age, gender, and tumor stage to construct the composite nomogram. The charts indicated that the AUC of cases with survival time of 1, 3 and 5 years are 0.698, 0.71 and 0.73. In addition, the mutation frequency of patients in the high-risk group was higher than that in the low-risk group. Simultaneously, DEGs were mainly enriched in ferroptosis-related pathways by analyzing the functional results.Conclusion. This study constructed a novel LUAD prognosis prediction model base on 5 ferroptosis-related genes, which can provide a prognostic evaluation tool for the clinical therapeutic decision.

Download Full-text

Prognostic Value of an Immune-Related Gene Signature in Oral Squamous Cell Carcinoma

Frontiers in Oncology ◽

10.3389/fonc.2021.776979 ◽

2021 ◽

Vol 11 ◽

Author(s):

Chao Zhu ◽

Liqun Gu ◽

Mianfeng Yao ◽

Jiang Li ◽

Changyun Fang

Keyword(s):

Overall Survival ◽

Squamous Cell Carcinoma ◽

Regression Analysis ◽

Oral Squamous Cell Carcinoma ◽

Risk Score ◽

Cox Regression ◽

Gene Signature ◽

Related Gene ◽

Cox Regression Analysis ◽

Immune Related Gene

The prognosis and immunotherapy response rates are unfavorable in patients with oral squamous cell carcinoma (OSCC). The tumor microenvironment is associated with tumor prognosis and progression, and the underlying mechanisms remain unclear. We obtained differentially expressed immune-related genes from OSCC mRNA data in The Cancer Genome Atlas (TCGA) database. Overall survival-related risk signature was constructed by univariate Cox regression analysis and LASSO Cox regression analysis. The prognostic performance was validated with receiver operating characteristic (ROC) analysis and Kaplan–Meier survival curves in the TCGA and Gene Expression Omnibus (GEO) datasets. The risk score was confirmed to be an independent prognostic factor and a nomogram was built to quantify the risk of outcome for each patient. Furthermore, a negative correlation was observed between the risk score and the infiltration rate of immune cells, as well as the expression of immunostimulatory and immunosuppressive molecules. Functional enrichment analysis between different risk score subtypes detected multiple immune-related biological processes, metabolic pathways, and cancer-related pathways. Thus, the immune-related gene signature can predict overall survival and contribute to the personalized management of OSCC patients.

Download Full-text

A methylation-based nomogram for predicting survival in patients with lung adenocarcinoma

BMC Cancer ◽

10.1186/s12885-021-08539-4 ◽

2021 ◽

Vol 21 (1) ◽

Author(s):

Xuelong Wang ◽

Bin Zhou ◽

Yuxin Xia ◽

Jianxin Zuo ◽

Yanchao Liu ◽

...

Keyword(s):

Risk Factors ◽

Dna Methylation ◽

Overall Survival ◽

Regression Analysis ◽

Lung Adenocarcinoma ◽

Risk Score ◽

Cox Regression ◽

Risk Group ◽

Clinical Risk Factors ◽

Cox Regression Analysis

Abstract Background DNA methylation alteration is frequently observed in Lung adenocarcinoma (LUAD) and may play important roles in carcinogenesis, diagnosis, and prognosis. Thus, this study aimed to construct a reliable methylation-based nomogram, guiding prognostic classification screening and personalized medicine for LUAD patients. Method The DNA methylation data, gene expression data and corresponding clinical information of lung adenocarcinoma samples were extracted from The Cancer Genome Atlas (TCGA) database. Differentially methylated sites (DMSs) and differentially expressed genes (DEGs) were obtained and then calculated correlation by pearson correlation coefficient. Functional enrichment analysis and Protein-protein interaction network were used to explore the biological roles of aberrant methylation genes. A prognostic risk score model was constructed using univariate Cox and LASSO analysis and was assessed in an independent cohort. A methylation-based nomogram that included the risk score and the clinical risk factors was developed, which was evaluated by concordance index and calibration curves. Result We identified a total of 1362 DMSs corresponding to 471 DEGs with significant negative correlation, including 752 hypermethylation sites and 610 hypomethylation sites. Univariate cox regression analysis showed that 59 DMSs were significantly associated with overall survival. Using LASSO method, we constructed a three-DMSs signature that was independent predictive of prognosis in the training cohort. Patients in high-risk group had a significant shorter overall survival than patients in low-risk group classified by three-DMSs signature (log-rank p = 1.9E-04). Multivariate cox regression analysis proved that the three-DMSs signature was an independent prognostic factor for LUAD in TCGA-LUAD cohort (HR = 2.29, 95%CI: 1.47–3.57, P = 2.36E-04) and GSE56044 cohort (HR = 2.16, 95%CI: 1.19–3.91, P = 0.011). Furthermore, a nomogram, combining the risk score with clinical risk factors, was developed with C-indexes of 0.71 and 0.70 in TCGA-LUAD and GSE56044 respectively. Conclusions The present study established a robust three-DMSs signature for the prediction of overall survival and further developed a nomogram that could be a clinically available guide for personalized treatment of LUAD patients.

Download Full-text

Expression Status and Prognostic Value of m6A RNA Methylation Regulators in Lung Adenocarcinoma

Life ◽

10.3390/life11070619 ◽

2021 ◽

Vol 11 (7) ◽

pp. 619

Author(s):

Xiuhong Li ◽

Zian Feng ◽

Rui Wang ◽

Jie Hu ◽

Xiaodong He ◽

...

Keyword(s):

Regression Analysis ◽

Survival Analysis ◽

Lung Adenocarcinoma ◽

Roc Curve ◽

Cox Regression ◽

Cox Regression Analysis ◽

Rna Methylation ◽

Kaplan Meier ◽

M6a Rna Methylation ◽

Prediction Efficiency

N6-methyladenosine (m6A) RNA modification is the most abundant modification method in mRNA, and it plays an important role in the occurrence and development of many cancers. This paper mainly discusses the role of m6A RNA methylation regulators in lung adenocarcinoma (LUAD) to identify novel prognostic biomarkers. The gene expression data of 19 m6A methylation regulators in LUAD patients and its relevant clinical parameters were extracted from The Cancer Genome Atlas (TCGA) database. We selected three significantly differentially expressed m6A regulators in LUAD to construct the risk signature, and evaluated its prognostic prediction efficiency using the receiver operating characteristic (ROC) curve. Kaplan–Meier survival analysis and Cox regression analysis were used to identify the independent prognostic significance of the risk signature. The ROC curve indicated that the area under the curve (AUC) was 0.659, which means that the risk signature had a good prediction efficiency. The results of the Kaplan–Meier survival analysis and Cox regression analysis showed that the risk score can be used as an independent prognostic factor for LUAD. In addition, we explored the differential signaling pathways and cellular processes related to m6A methylation regulators in LUAD.

Download Full-text

Identification of a circadian gene signature that predicts overall survival in lung adenocarcinoma

PeerJ ◽

10.7717/peerj.11733 ◽

2021 ◽

Vol 9 ◽

pp. e11733

Author(s):

Xinliang Gao ◽

Mingbo Tang ◽

Suyan Tian ◽

Jialin Li ◽

Wei Liu

Keyword(s):

Overall Survival ◽

Regression Analysis ◽

Lung Adenocarcinoma ◽

Circadian Clock ◽

Cox Regression ◽

Gene Signature ◽

Differentially Expressed ◽

Circadian Gene ◽

Cox Regression Analysis ◽

Circadian Genes

Background Lung adenocarcinoma (LUAD) is one of the most common subtypes of lung cancer which is the leading cause of death in cancer patients. Circadian clock disruption has been listed as a likely carcinogen. However, whether the expression of circadian genes affects overall survival (OS) in LUAD patients remains unknown. In this article, we identified a circadian gene signature to predict overall survival in LUAD. Methods RNA sequencing (HTSeq-FPKM) data and clinical characteristics were obtained for a cohort of LUAD patients from The Cancer Genome Atlas (TCGA). A multigene signature based on differentially expressed circadian clock-related genes was generated for the prediction of OS using Least Absolute Shrinkage and Selection Operator (LASSO)-penalized Cox regression analysis, and externally validated using the GSE72094 dataset from the GEO database. Results Five differentially expressed genes (DEGs) were identified to be significantly associated with OS using univariate Cox proportional regression analysis (P < 0.05). Patients classified as high risk based on these five DEGs had significantly lower OS than those classified as low risk in both the TGCA cohort and GSE72094 dataset (P < 0.001). Multivariate Cox regression analysis revealed that the five-gene-signature based risk score was an independent predictor of OS (hazard ratio > 1, P < 0.001). Receiver operating characteristic (ROC) curves confirmed its prognostic value. Gene set enrichment analysis (GSEA) showed that Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways related to cell proliferation, gene damage repair, proteasomes, and immune and autoimmune diseases were significantly enriched. Conclusion A novel circadian gene signature for OS in LUAD was found to be predictive in both the derivation and validation cohorts. Targeting circadian genes is a potential therapeutic option in LUAD.

Download Full-text