Identification and Validation of a Seven m6A-related lncRNAs Signature Predicting Prognosis of Ovarian Cancer

Author(s):  
Yan Li ◽  
Xiaoying Wang ◽  
Yue Han ◽  
Xun Li

Abstract Background: Long non-coding RNAs (lncRNAs) play an important role in angiogenesis, immune response, inflammatory response and tumor development and metastasis. m6 A (N6 - methyladenosine) is one of the most common RNA modifications in eukaryotes. The aim of our research was to investigate the potential prognostic value of m6A-related lncRNAs in ovarian cancer (OC).Methods: The data we need for our research was downloaded from the Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) database. Pearson correlation analysis between 21 m6A regulators and lncRNAs was performed to identify m6A-related lncRNAs. Univariate Cox regression analysis was implemented to screen for lncRNAs with prognostic value. A least absolute shrinkage and selection operator (LASSO) Cox regression and multivariate Cox regression analyses was used to further reduct the lncRNAs with prognostic value and construct a m6A-related lncRNAs signature for predicting the prognosis of OC patients. Results: 275 m6A-related lncRNAs were obtained using pearson correlation analysis. 29 m6A-related lncRNAs with prognostic value was selected through univariate Cox regression analysis. Then, a seven m6A-related lncRNAs signature was identified by LASSO Cox regression. Each patient obtained a riskscore through multivariate Cox regression analyses and the patients were classified into high-and low-risk group using the median riskscore as a cutoff. Kaplan-Meier curve revealed that the patients in high-risk group have poor outcome. The receiver operating characteristic curve revealed that the predictive potential of the m6A-related lncRNAs signature for OC was powerful. The predictive potential of the m6A-related lncRNAs signature was successfully validated in the GSE9891, GSE26193 datasets and our clinical specimens. Multivariate analyses suggested that the m6A-related lncRNAs signature was an independent prognostic factor for OC patients. Moreover, a nomogram based on the expression level of the seven m6A-related lncRNAs was established to predict survival rate of patients with OC. Finally, a competing endogenous RNA (ceRNA) network associated with the seven m6A-related lncRNAs was constructed to understand the possible mechanisms of the m6A-related lncRNAs involed in the progression of OC.Conclusions: In conclusion, our research revealed that the m6A-related lncRNAs may affect the prognosis of OC patients and identified a seven m6A-related lncRNAs signature to predict the prognosis of OC patients.

2021 ◽  
Vol 20 ◽  
pp. 153303382110414
Author(s):  
Xiaoyong Li ◽  
Jiaqong Lin ◽  
Yuguo pan ◽  
Peng Cui ◽  
Jintang Xia

Background: Liver progenitor cells (LPCs) play significant roles in the development and progression of hepatocellular carcinoma (HCC). However, no studies on the value of LPC-related genes for evaluating HCC prognosis exist. We developed a gene signature of LPC-related genes for prognostication in HCC. Methods: To identify LPC-related genes, we analyzed mRNA expression arrays from a dataset (GSE57812 & GSE 37071) containing LPCs, mature hepatocytes, and embryonic stem cell samples. HCC RNA-Seq data from The Cancer Genome Atlas (TCGA) were used to explore the differentially expressed genes (DEGs) related to prognosis through DEG analysis and univariate Cox regression analysis. Lasso and multivariate Cox regression analyses were performed to construct the LPC-related gene prognostic model in the TCGA training dataset. This model was validated in the TCGA testing set and an external dataset (International Cancer Genome Consortium [ICGC] dataset). Finally, we investigated the relationship between this prognostic model with tumor-node-metastasis stage, tumor grade, and vascular invasion of HCC. Results: Overall, 1770 genes were identified as LPC-related genes, of which 92 genes were identified as DEGs in HCC tissues compared with normal tissues. Furthermore, we randomly assigned patients from the TCGA dataset to the training and testing cohorts. Twenty-six DEGs correlated with overall survival (OS) in the univariate Cox regression analysis. Lasso and multivariate Cox regression analyses were performed in the TCGA training set, and a 3-gene signature was constructed to stratify patients into 2 risk groups: high-risk and low-risk. Patients in the high-risk group had significantly lower OS than those in the low-risk group. Receiver operating characteristic curve analysis confirmed the signature's predictive capacity. Moreover, the risk score was confirmed to be an independent predictor for patients with HCC. Conclusion: We demonstrated that the LPC-related gene signature can be used for prognostication in HCC. Thus, targeting LPCs may serve as a therapeutic alternative for HCC.


2020 ◽  
Author(s):  
Zhihao Wang ◽  
Kidane Siele Embaye ◽  
Qing Yang ◽  
Lingzhi Qin ◽  
Chao Zhang ◽  
...  

Abstract Background: Given that metabolic reprogramming has been recognized as an essential hallmark of cancer cells, this study sought to investigate the potential prognostic values of metabolism-related genes(MRGs) for hepatocellular carcinoma (HCC) diagnosis and treatment. Methods: The metabolism-related genes sequencing data of HCC samples with clinical information were obtained from the International Cancer Genome Consortium(ICGC) and The Cancer Genome Atlas (TCGA). The differentially expressed MRGs were identified by Wilcoxon rank sum test. Then, univariate Cox regression analysis were performed to identify metabolism-related DEGs that related to overall survival(OS). A novel metabolism-related prognostic signature was developed using the least absolute shrinkage and selection operator (Lasso) and multivariate Cox regression analyses . Furthermore, the signature was validated in the TCGA dataset. Finally, cox regression analysis was applied to identify the prognostic value and clinical relationship of the signature in HCC. Results: A total of 178 differentially expressed MRGs were detected between the ICGA dataset and the TCGA dataset. We found that 17 MRGs were most significantly associated with OS by using the univariate Cox proportional hazards regression analysis in HCC. Then, the Lasso and multivariate Cox regression analyses were applied to construct the novel metabolism-relevant prognostic signature, which consisted of six MRGs. The prognostic value of this prognostic model was further successfully validated in the TCGA dataset. Further analysis indicated that this signature could be an independent prognostic indicator after adjusting to other clinical factors. Six MRGs (FLVCR1, MOGAT2, SLC5A11, RRM2, COX7B2, and SCN4A) showed high prognostic performance in predicting HCC outcomes, and were further associated with tumor TNM stage, gender, age, and pathological stage. Finally, the signature was found to be associated with various clinicopathological features. Conclusions: In summary, our data provided evidence that the metabolism-based signature could serve as a reliable prognostic and predictive tool for overall survival in patients with HCC.


2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Jianfeng Zheng ◽  
Jialu Guo ◽  
Benben Cao ◽  
Ying Zhou ◽  
Jinyi Tong

Abstract Background Both N6-methyladenosine (m6A) modification and lncRNAs play an important role in the carcinogenesis and cancer inhibition of ovarian cancer (OC). However, lncRNAs involved in m6A regulation (LI-m6As) have never been reported in OC. Herein, we aimed to identify and validate a signature based on LI-m6A for OC. Methods RNA sequencing profiles with corresponding clinical information associated with OC and 23 m6A regulators were extracted from TCGA. The Pearson correlation coefficient (PCC) between lncRNAs and 23 m6A regulators (|PCC|> 0.4 and p < 0.01) was calculated to identify LI-m6As. The LI-m6As with significant prognostic value were screened based on univariate Cox regression analysis to construct a risk model by LASSO Cox regression. Gene Set Enrichment Analysis (GSEA) was implemented to survey the biological functions of the risk groups. Several clinicopathological characteristics were utilized to evaluate their ability to predict prognosis, and a nomogram was constructed to evaluate the accuracy of survival prediction. Besides, immune microenvironment, checkpoint, and drug sensitivity in the two risk groups were compared using comprehensive algorithms. Finally, real-time qPCR analysis and cell counting kit-8 assays were performed on an alternative lncRNA, CACNA1G-AS1. Results The training cohort involving 258 OC patients and the validation cohort involving 111 OC patients were downloaded from TCGA. According to the PCC between the m6A regulators and lncRNAs, 129 LI-m6As were obtained to perform univariate Cox regression analysis and then 10 significant prognostic LI-m6As were identified. A prognostic signature containing four LI-m6As (AC010894.3, ACAP2-IT1, CACNA1G-AS1, and UBA6-AS1) was constructed according to the LASSO Cox regression analysis of the 10 LI-m6As. The prognostic signature was validated to show completely opposite prognostic value in the two risk groups and adverse overall survival (OS) in several clinicopathological characteristics. GSEA indicated that differentially expressed genes in disparate risk groups were enriched in several tumor-related pathways. At the same time, we found significant differences in some immune cells and chemotherapeutic agents between the two groups. An alternative lncRNA, CACNA1G-AS1, was proven to be upregulated in 30 OC specimens and 3 OC cell lines relative to control. Furthermore, knockdown of CACNA1G‐AS1 was proven to restrain the multiplication capacity of OC cells. Conclusions Based on the four LI-m6As (AC010894.3, ACAP2-IT1, CACNA1G-AS1, and UBA6-AS1), the risk model we identified can independently predict the OS and therapeutic value of OC. CACNA1G‐AS1 was preliminarily proved to be a malignant lncRNA.


2022 ◽  
Vol 11 ◽  
Author(s):  
Yue Wang ◽  
Bao Xuan Li ◽  
Xiang Li

Ovarian cancer (OC) is a highly heterogeneous disease with different cellular origins reported; thus, precise prognostic strategies and effective new therapies are urgently needed for patients with OC. A growing number of studies have shown that most malignancies have intensive angiogenesis and rapid growth. Therefore, angiogenesis plays an important role in the development of tumor metastasis. However, the prognostic value of angiogenesis-related genes (ARGs) in OC remains to be further elucidated. In this study, the expression data and corresponding clinical data from patients with OC and normal control samples were downloaded with UCSC XENA. A total of 1,960 differentially expressed ARGs were screened and functionally annotated through Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Univariate Cox regression analysis was performed to identify ARGs associated with prognosis. New ARGs signatures (including ESM1, CXCL13, TPCN2, PTPRD, FOXO1, and ELK3) were constructed for the prediction of overall survival (OS) in OC based on the least absolute shrinkage and selection operator (LASSO) and multivariate Cox regression analysis. Patients were divided based on their median risk score. In the The Cancer Genome Atlas (TCGA) training dataset, the survival analysis showed that overall survival was lower in the high-risk group than that in the low-risk group (p &lt; 0.0001). The International Cancer Genome Consortium (ICGC) database was used for validation, and the receiver operating characteristic (ROC) curves showed good performance. Univariate and multivariate Cox analyses were conducted to identify independent predictors of OS. The nomogram, including the risk score, age, stage, grade, and position, can not only show good predictive ability but also can explore the correlation analysis based on ARGs for immunogenicity, immune components, and immune phenotypes with risk score. Risk scores were correlated strongly with the type of immune infiltration. Furthermore, homologous recombination defect (HRD), NtAIscore, LOH score, LSTm score, stemness index (mRNAsi), and stromal cells were significantly correlated with risk score. The present study suggests that the novel signature constructed from six ARGs may serve as effective prognostic biomarkers for OC and contribute to clinical decision making and personalized prognostic monitoring of OC.


2021 ◽  
Author(s):  
Yan Li ◽  
Yue Han ◽  
Xiaoyin Wang

Abstract Background: The mortality rate of ovarian cancer (OC) ranks the first in gynecological tumors, which seriously threatens women's health and life. In recent years, alternative splicing (AS) has gradually been considered to play a key role in immune infiltrates of tumor. However, the prognostic significance of AS events related to immune infiltrates in OC remains unknow. The aim of our research was to investigate the potential prognostic value of AS events associated with immune infiltrates in OC.Methods: The RNA sequences (RNA-seq) and clinical data were downloaded from the Cancer Genome Atlas (TCGA) database. The AS event data was obtained from TCGA SpliceSeq database. Single sample gene set enrichment analysis (ssGSEA) was performed to calculate the abundance of 28 immune cell types in samples from TCGA-OV dataset. A consensus clustering algorithm was used to group the OC patients. Differential expression analysis was used to identify differentially expressed AS events (DEASs) between groups. Univariate Cox regression analysis was implemented to screen for AS events with prognostic value. A least absolute shrinkage and selection operator (LASSO) Cox regression and multivariate Cox regression analyses was used to further narrow the AS events with prognostic value and construct an alternative splicing prognostic signature for predicting the prognosis of OC patients. Results: The OC patients from TCGA database were classified into two groups (cluster.A and cluster.B) based on the 28 types of TIIC using a consensus clustering algorithm. The patients in the cluster.A group had increased immune infiltrates compared with the cluster.B group. 3616 DEASs were acquired and 171 DEASs have prognostic value (p<0.05). 28 DEASs with prognostic value (p<0.001) were fitted into LASSO Cox regression and multivariate Cox regression analyses. A prognostic signature with 18 DEASs was constructed to predict the prognosis of OC patients. Each patient obtained a riskscore and the patients were classified into high-and low-risk group using the median riskscore as a cutoff. Kaplan-Meier curve revealed that the patients in high-risk group have poor outcome. Conclusions: Collectively, our research identified an alternative splicing prognostic signature associated with immune infiltrates of OC, which may provide new directions for the immunotherapy of OC patients.


2022 ◽  
Vol 2022 ◽  
pp. 1-10
Author(s):  
Xiaoting Zhang ◽  
Yue Su ◽  
Xian Fu ◽  
Jing Xiao ◽  
Guicheng Qin ◽  
...  

Lung squamous cell carcinoma (LUSC) is the most common type of lung cancer accounting for 40% to 51%. Long noncoding RNAs (lncRNAs) have been reported to play a significant role in the invasion, migration, and proliferation of lung cancer tissue cells. However, systematic identification of lncRNA signatures and evaluation of the prognostic value for LUSC are still an urgent problem. In this work, LUSC RNA-seq data were collected from TCGA database, and the limma R package was used to screen differentially expressed lncRNAs (DElncRNAs). In total, 216 DElncRNAs were identified between the LUSC and normal samples. lncRNAs associated with prognosis were calculated using univariate Cox regression analysis. The overall survival (OS) prognostic model containing 10 lncRNAs and the disease-free survival (DFS) prognostic model consisting of 11 lncRNAs were constructed using a machine learning-based algorithm, systematic LASSO-Cox regression analysis. We found that the survival rate of samples in the high-risk group was lower than that in the low-risk group. Results of ROC curves showed that both the OS and DFS risk score had better prognostic effects than the clinical characteristics, including age, stage, gender, and TNM. Two lncRNAs (LINC00519 and FAM83A-AS1) that were commonly identified as prognostic factors in both models could be further investigated for their clinical significance and therapeutic value. In conclusion, we constructed lncRNA prognostic models with considerable prognostic effect for both OS and DFS of LUSC.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Yan Ouyang ◽  
Kaide Xia ◽  
Xue Yang ◽  
Shichao Zhang ◽  
Li Wang ◽  
...  

AbstractAlternative splicing (AS) events associated with oncogenic processes present anomalous perturbations in many cancers, including ovarian carcinoma. There are no reliable features to predict survival outcomes for ovarian cancer patients. In this study, comprehensive profiling of AS events was conducted by integrating AS data and clinical information of ovarian serous cystadenocarcinoma (OV). Survival-related AS events were identified by Univariate Cox regression analysis. Then, least absolute shrinkage and selection operator (LASSO) and multivariate Cox regression analysis were used to construct the prognostic signatures within each AS type. Furthermore, we established a splicing-related network to reveal the potential regulatory mechanisms between splicing factors and candidate AS events. A total of 730 AS events were identified as survival-associated splicing events, and the final prognostic signature based on all seven types of AS events could serve as an independent prognostic indicator and had powerful efficiency in distinguishing patient outcomes. In addition, survival-related AS events might be involved in tumor-related pathways including base excision repair and pyrimidine metabolism pathways, and some splicing factors might be correlated with prognosis-related AS events, including SPEN, SF3B5, RNPC3, LUC7L3, SRSF11 and PRPF38B. Our study constructs an independent prognostic signature for predicting ovarian cancer patients’ survival outcome and contributes to elucidating the underlying mechanism of AS in tumor development.


2021 ◽  
Author(s):  
Shaopei Ye ◽  
Wenbin Tang ◽  
Ke Huang

Abstract Background: Autophagy is a biological process to eliminate dysfunctional organelles, aggregates or even long-lived proteins. . Nevertheless, the potential function and prognostic values of autophagy in Wilms Tumor (WT) are complex and remain to be clarifed. Therefore, we proposed to systematically examine the roles of autophagy-associated genes (ARGs) in WT.Methods: Here, we obtained differentially expressed autophagy-related genes (ARGs) between healthy and Wilms tumor from Therapeutically Applicable Research To Generate Effective Treatments(TARGET) and The Cancer Genome Atlas (TCGA) database. The functionalities of the differentially expressed ARGs were analyzed using Gene Ontology. Then univariate COX regression analysis and multivariate COX regression analysis were performed to acquire nine autophagy genes related to WT patients’ survival. According to the risk score, the patients were divided into high-risk and low-risk groups. The Kaplan-Meier curve demonstrated that patients with a high-risk score tend to have a poor prognosis.Results: Eighteen DEARGs were identifed, and nine ARGs were fnally utilized to establish the FAGs based signature in the TCGA cohort. we found that patients in the high-risk group were associated with mutations in TP53. We further conducted CIBERSORT analysis, and found that the infiltration of Macrophage M1 was increased in the high-risk group. Finally, the expression levels of crucial ARGs were verifed by the experiment, which were consistent with our bioinformatics analysis.Conclusions: we emphasized the clinical significance of autophagy in WT, established a prediction system based on autophagy, and identified a promising therapeutic target of autophagy for WT.


2021 ◽  
Author(s):  
Sijia Li ◽  
Hongyang Zhang ◽  
Wei Li

Abstract Background: The purpose of our study is establishing a model based on ferroptosis-related genes predicting the prognosis of patients with head and neck squamous cell carcinoma (HNSCC).Methods: In our study, transcriptome and clinical data of HNSCC patients were from The Cancer Genome Atlas, ferroptosis-related genes and pathways were from Ferroptosis Signatures Database. Differentially expressed genes (DEGs) were screened by comparing tumor and adjacent normal tissues. Functional enrichment analysis of DEGs, protein-protein interaction network and gene mutation examination were applied. Univariate Cox regression analysis and least absolute shrinkage and selection operator (LASSO) regression were used to identified DEGs. The model was constructed by multivariate Cox regression analysis and verified by Kaplan-Meier analysis. The relationship between risk scores and other clinical features was also analyzed. Univariate and multivariate Cox analysis was used to verified the independence of our model. The model was evaluated by receiver operating characteristic analysis and calculation of the area under the curve (AUC). A nomogram model based on risk score, age, gender and TNM stages was constructed.Results: We analyzed data including 500 tumor tissues and 44 adjacent normal tissues and 259 ferroptosis-related genes, then obtained 73 DEGs. Univariate Cox regression analysis screened out 16 genes related to overall survival, and LASSO analysis fingered out 12 of them with prognostic value. A risk score model based on these 12 genes was constructed by multivariate Cox regression analysis. According to the median risk score, patients were divided into high-risk group and low-risk group. The survival rate of high-risk group was significantly lower than that of low-risk group in Kaplan-Meier curve. Risk scores were related to T and grade. Univariate and multivariate Cox analysis showed our model was an independent prognostic factor. The AUC was 0.669. The nomogram showed high accuracy predicting the prognosis of HNSCC patients.Conclusion: Our model based on 12 ferroptosis-related genes performed excellently in predicting the prognosis of HNSCC patients. Ferroptosis-related genes may be promising biomarkers for HNSCC treatment and prognosis.


2021 ◽  
Vol 39 (15_suppl) ◽  
pp. e17543-e17543
Author(s):  
Xiaoxiang Chen ◽  
Jing Ni ◽  
Xia Xu ◽  
Wenwen Guo ◽  
Xianzhong Cheng ◽  
...  

e17543 Background: Homologous recombination deficiency (HRD) is the first phenotypically defined predictive biomarker for Poly (ADP-ribose) polymerase inhibitors (PARPi) in ovarian cancer. However, the proportion of HRD positive in real world and the relationship of HRD status with PARPi in Chinese ovarian cancer patients remains unknown. Methods: A total of sixty-four ovarian cancer patients underwent PARPi, both Olaparib and Niraparib, were enrolled from August 2018 to January 2021 in Jiangsu Institute of Cancer Hospital. HRD score which was the sum of loss of heterozygosity (LOH), telomeric allelic imbalance (TAI) and large-scale state transitions (LST) events were calculated using tumor DNA-based next generation sequencing (NGS) assays. HRD-positive was defined by either BRCA1/2 pathogenic or likely pathogenic mutation or HRD score ≥42. Progression-free survival (PFS) was analyzed with a log-rank test using HRD status and summarized using Kaplan-Meier methodology. Univariate and multiple cox-regression analysis were conducted to investigate all possible clinical factors. Results: 71.9% (46/64) patients were HRD positive and the rest 28.1% (18/64) were HRD negative, which was higher than the HRD positive proportion reported in Western countries. The PFS among HRD positive patients was significantly longer than those HRD negative patients (medium PFS 8.9 m vs 3.6 m, hazard ratio [HR]: 0.22, p < 0.001). Among them, 23 patients who were BRCA wild type but HRD positive had longer PFS than those with BRCA wild type and HRD negative (medium PFS 9.2 m vs 3.6 m, HR: 0.20, p < 0.001). Univariate cox-regression analysis found that HRD status, previous treatment lines, secondary cytoreductive surgery (SCS) were significantly associated with PFS after PARPi treatment. After multiple regression correction, HRD status (HR: 0.39, 95% CI: [0.20-0.76], p = 0.006), ECOG score (HR: 2.53, 95% CI: [1.24-5.17], p = 0.011) and SCS (HR: 2.21, 95% CI: [1.09-4.48], p = 0.028) were the independent factors. Subgroup analysis in ECOG = 0 subgroup (N = 36), HRD positive patients had significant longer PFS than HRD negative patients (medium PFS 10.3 m vs 5.8 m, HR: 0.14, p < 0.001). Also in the subgroup of patients without SCS, PFS in patients with HRD was longer than patients without HRD (medium PFS 10.2 m vs 5.7 m, HR: 0.29, p = 0.003). Conclusions: This is the first real-world data of HRD status in ovarian cancer patients from China and demonstrate that HRD is a valid biomarker for PARP inhibitors in Chinese ovarian cancer patients.


Sign in / Sign up

Export Citation Format

Share Document