- 1Department of Radiation Oncology, Inner Mongolia Cancer Hospital and The Affiliated People’s Hospital of Inner Mongolia Medical University, Hohhot, China
- 2Institute for Endemic Fluorosis Control, Center for Endemic Disease Control, Chinese Center for Disease Control and Prevention, Harbin Medical University, Harbin, China
- 3Department of Obstetrics and Gynecology, The First Affiliated Hospital of Harbin Medical University, Harbin, China
- 4Department of Pathology, Harbin Medical University Cancer Hospital, Harbin, China
The deregulation of long non-coding RNAs (lncRNAs) by epigenetic alterations has been implicated in cancer initiation and progression. However, the epigenetically regulated lncRNAs and their association with clinical outcome and therapeutic response in ovarian cancer (OV) remain poorly investigated. This study performed an integrative analysis of DNA methylation data and transcriptome data and identified 419 lncRNAs as potential epigenetically regulated lncRNAs. Using machine-learning and multivariate Cox regression analysis methods, we identified and developed an epigenetically regulated lncRNA expression signature (EpiLncRNASig) consisting of five lncRNAs from the list of 17 epigenetically regulated lncRNAs significantly associated with outcome. The EpiLncRNASig could stratify patients into high-risk groups and low-risk groups with significantly different survival and chemotherapy response in different patient cohorts. Multivariate Cox regression analyses, after adjusted by other clinical features and treatment response, demonstrated the independence of the DEpiLncSig in predicting survival. Functional analysis for relevant protein-coding genes of the DEpiLncSig indicated enrichment of known immune-related or cancer-related biological pathways. Taken together, our study not only provides a promising prognostic biomarker for predicting outcome and chemotherapy response but also will improve our understanding of lncRNA epigenetic regulation mechanisms in OV.
Introduction
Ovarian cancer (OV) is one of the most lethal gynecologic cancers and is the eighth leading cause of cancer-related deaths in women (Coburn et al., 2017; Momenimovahed et al., 2019). According to the American Cancer Society, approximately 21,750 new cases and 13,940 estimated deaths occur in 2020 in the United States (Siegel et al., 2020). Surgery followed by chemo-based treatments (platinum-based or taxane-based) is the first-line treatment, and approximately 80% of OV patients initially respond to treatment. However, most patients with advanced disease still experienced recurrent disease after a period of chemotherapy, which has led to poor outcomes with five survival rates of less than 50%, even if significant advances have been made in surgical and chemo-based treatments for OV patients. Therefore, it is critical to identify reliable and useful biomarkers for improving outcomes of OV patients.
Increasing efforts in studying molecular omics have improved our understanding of the molecular mechanisms of OV carcinogenesis and progression and contributed to the identification and development of novel molecular biomarkers and specific therapies (Cancer Genome Atlas Research Network, 2011; Lu et al., 2014; Clifford et al., 2018; Li et al., 2020). Molecular profiles have been extensively investigated and characterized during the past years, leading to the identification of a number of dysregulated molecules associated with development, progression, recurrence, metastasis, and therapeutic response of OV (Adib et al., 2004; Barrett et al., 2015; Sallinen et al., 2019; Wang et al., 2019; Zhao H. et al., 2020). Long non-coding RNAs (lncRNAs) are a significant class of non-coding RNAs (ncRNAs) that have been discovered in the last 10 years. The accumulated evidence has shown that lncRNAs are essential regulators in gene expression networks and performed a variety of functional roles in almost all biological processes through diverse mechanisms at transcriptional, post-transcriptional, and epigenetic levels (Mercer et al., 2009; Marchese et al., 2017; Yao et al., 2019). Dysregulated lncRNAs have been observed in various cancers and have also been implicated as novel biomarkers in cancer diagnosis, prognosis, and therapeutic response (Zhou et al., 2016a, 2020; Hosseini et al., 2017; Huo et al., 2017; Tian et al., 2018; Bao et al., 2020; Sun et al., 2020). Previous studies have reported that the expression of some lncRNAs can be regulated by epigenetic modification similar to protein-coding genes (Zhao et al., 2015; Wang et al., 2018; Li et al., 2019). However, the epigenetically regulated lncRNAs and their association with clinical outcome and therapeutic response in OV remain poorly investigated.
In this study, we performed an integrative analysis of DNA methylation data and transcriptome data to identify epigenetically regulated lncRNAs associated with OV carcinogenesis. Together with clinical data, we further investigated the clinical value of these epigenetically regulated lncRNAs in predicting outcome and chemotherapy response. Finally, we used the relevant genes of epigenetically regulated lncRNAs to infer their potential functions.
Materials and Methods
Acquisition and Analysis of DNA Methylation Data for OV Patients
DNA methylation data profiled by Illumina 27k methylation array for 583 OV tumor tissues and 12 non-cancer tissues were obtained from the UCSC Xena Browser1. Then probes were filtered and prepossessed as follows: (i) SNP-enriched probes were removed; (ii) those probes with missing values in more than 10% samples were removed; and (iii) imputation was performed to replace missing values by calculating the median methylation level for each sample across all CpG sites. Finally, 24793 CpG probes were kept for further analysis. Differential methylation analysis on the site level between OV tumor tissues and non-cancer tissues was conducted using the limma package. Differentially methylated CpG sites were determined based on FDR adjusted p-value < 0.05 and absolute mean methylation difference >0.3.
Acquisition and Analysis of lncRNA Data of OV Patient
RNA-seq data of OV tumor tissues and non-cancer tissues were obtained from the UCSC Xena Browser (see text footnote 1). A total of 14,614 were obtained from the RNA-seq data based on GENCODE annotations2. The epigenetically regulated lncRNAs were identified by measuring the association between lncRNA expression and CpG levels through the Pearson correlation coefficient. Those lncRNAs correlated with CpG sites with | r| > 0.4 and p < 0.001 were considered as epigenetically regulated lncRNAs (DEpiLncRNAs). Hierarchical clustering analysis was conducted using the R package “pheatmap” with “ward.D2” method.
Development of an Epigenetically Regulated Five-lncRNA Signature (EpiLncRNASig)
Univariate Cox regression analysis of DEpiLncRNAs with OS was performed to identify candidate DEpiLncRNAs with prognostic roles. Using a bidirectional elimination strategy, the stepwise regression was performed to select optimal lncRNA biomarkers from the list of candidate prognostic DEpiLncRNAs by an automatic procedure. Specifically, at each step, we added variables whose inclusion gave the most significant improvement and removed variables whose exclusion gave the most insignificant deterioration to the quality of the prediction model, which is assessed by the Akaike information criterion (AIC). Then, we repeated this process until none of them improved the model to a statistically significant extent (Zhou et al., 2020). Then a prognostic lncRNA-focused score model (EpiLncRNASig) was constructed based on the linear combination of the expression of optimal DEpiLncRNAs biomarkers, weighted by their coefficients from the multivariate regression analysis as in previous studies (Bao et al., 2020, 2021; Sun et al., 2020).
Statistical Analysis
Univariate and multivariate Cox regression analyses were performed on the individual clinical variables with and without the EpiLncRNASig. Hazard ratios (HR) and 95% confidence intervals (CI) were calculated. Kaplan-Meier survival plots and log-rank tests were applied to compare the differences in survival time between different patient groups using the R package “survival.” All statistical analyses were performed using R Statistical Software (version 3.6.3).
Function Enrichment Analysis
The association between expression levels of lncRNAs and mRNAs was measured by calculating the Pearson correlation coefficient, and the top 100 mRNAs were considered as lncRNA-related mRNAs. Then function enrichment analysis of Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) was performed to infer possible functional roles of lncRNA biomarkers using the R package “clusterProfiler” (Yu et al., 2012). GO terms or KEGG pathways with adjusted p-value < 0.05 were considered to be significantly enriched.
Results
Identification of Dysregulated Epigenetically Regulated lncRNAs (DEpiLncRNAs) Associated With OV Development
To identify potential GpG sites associated with OV development, we performed compared analysis for GpG levels between OV tumor tissues and non-cancer tissues and identified 605 differentially methylated sites (FDR adjusted p-value < 0.05 and absolute mean methylation difference >0.3) using the R package “RnBeads.” Hierarchical clustering analysis showed that these 605 differentially methylated sites could discriminate between OV patients with and healthy controls, as shown in Figure 1A. To identify DEpiLncRNAs in OV, we calculated the Pearson correlation coefficient to evaluate the association between lncRNA expression and CpG levels and identified 1,497 lncRNA-CpG pairs, including 419 lncRNAs, which were defined as DEpiLncRNAs. The network between lncRNA and CpG sites were visualized using the Cytoscape software (Figure 1B).
Figure 1. Identification of dysregulated epigenetically regulated lncRNAs (DEpiLncRNAs) associated with OV development. (A) Heatmap of differentially methylated sites could discriminate between OV patients with and healthy controls. (B) Interaction network of dysregulated epigenetically regulated lncRNAs and CpG.
Development of an Epigenetically Regulated lncRNA Signature (DEpiLncSig) for Predicting the Outcome
All TCGA-OV patients with survival information were split equally into the discovery cohort (n = 187) and validation cohort (n = 186). We first conducted the univariate analysis for 419 lncRNAs with OS and identified 17 lncRNAs that were significantly associated with OS (p < 0.01). Then we used a stepwise regression model and identified five optimal lncRNA biomarkers from the list of 17 lncRNAs significantly associated with OS (Table 1). Of them, three lncRNAs (LINC00189, CACNA1G-AS1, and AC105384.1) seem to be risk factors, and the remaining two lncRNAs (AL133467.1 and CHRM3-AS2) are protective factors. The five optimal lncRNA biomarkers and their CpGs were listed in Supplementary Table 1. Finally, we performed multivariate Cox regression analysis for five optimal lncRNA biomarkers and constructed a lncRNA-based risk score model based on the linear combination of the expression of five optimal lncRNA biomarkers, weighted by their coefficients from the multivariate regression analysis as follows: EpiLncRNASig = (0.197)∗expression (LINC00189) + (0.167)∗ expression (CACNA1G-AS1) + (−0.432)∗ expression (AL133467.1) + (−0.374)∗ expression (CHRM3-AS2) + (0.395)∗ expression (AC105384.1). Each patient in the discovery cohort was assigned a risk score and was subsequently defined as high-risk with higher EpiLncRNASig or low-risk with lower EpiLncRNASig according to the median value of the risk score (0.0538). Survival analysis showed that patients in the low-risk group have significantly improved OS (median time 1,742 days) compared to those in the high-risk group (median time 1,039 days) (log-rank test p < 0.0001) (Figure 2A). The 3- and 5-year survival rates of patients in the low-risk group are 75 and 48%, respectively, whereas the corresponding rates are 45 and 11%, respectively, in the high-risk group. Furthermore, there were significant differences in disease-free survival (median time 678 vs. 500 days), disease-specific survival (median time 1,933 vs. 1,039 days), and progression-free survival (median time 690 vs. 427 days) (Figure 2A). As shown in Figure 2B, three lncRNAs (LINC00189, CACNA1G-AS1, and AC105384.1) revealed significantly higher expression in the high-risk group compared to those in the low-risk group, and two lncRNAs (AL133467.1 and CHRM3-AS2) are significantly upregulated in the low-risk group and down-regulated in the high-risk group (Figure 2B).
Figure 2. Prognostic performance of epigenetically regulated lncRNA signature (DEpiLncSig) in the discovery cohort. (A) Kaplan-Meier survival curves of survival times between the high-risk group and low-risk group. (B) Boxplots of expression levels of prognostic lncRNAs in the high-risk group and low-risk group, respectively.
Further Validation of the DEpiLncSig in the Independent Patient Cohort
To test the robustness of the DEpiLncSig, the DEpiLncSig was applied to 186 OV patients in the validation cohort. The DEpiLncSig classified 186 OV patients of the validation cohort into the high-risk group (n = 92) and low-risk group (n = 94) according to the risk cutoff derived from the discovery cohort. As observed in the training group, the OS time of the low-risk group patients (median 1,736 days) was significantly better than that of high-risk group patients (median 1,264 days) (log-rank test p = 0.021) (Figure 3A). The 3- and 5-year survival rates of patients in the low-risk group were 70 and 41%, respectively, whereas the corresponding rates are 62 and 28% in the high-risk group, respectively. Furthermore, the high-risk patients have significantly shorter disease-free survival (median time 554 vs. 1,042 days, log-rank test p = 0.0087), disease-specific survival (median time 1,319 vs. 1,767 days, log-rank test p = 0.032), and progression-free survival (median time 454 vs. 553 days, log-rank test p = 0.042) compared to those in the low-risk group (Figures 3B–D). When the DEpiLncSig was further tested in the entire TCGA-OV patients, all TCGA-OV patients were divided into the high-risk group and low-risk group with significantly different overall survival (median time 1,163 vs. 1,736 days, log-rank test p < 0.0001), disease-free survival (median time 535 vs. 818 days, log-rank test p < 0.0001), disease-specific survival (median time 1,199 vs. 1,784 days, log-rank test p < 0.0001) and progression-free survival (median time 447 vs. 628 days, log-rank test p = 0.00011) (Figures 3E–H). These results demonstrated the stable and reliable prognostic performance of the DEpiLncSig in predicting the outcome of OV patients.
Figure 3. Independent validation of prognostic performance of epigenetically regulated lncRNA signature. Kaplan-Meier survival curves between the high-risk group and low-risk group for overall survival (A), progression-free survival (B), disease-specific survival (C), and disease-free survival (D) in the independent patient cohort. Kaplan-Meier survival curves between the high-risk group and low-risk group for overall survival (E), disease-free survival (F), disease-specific survival (G), and progression-free survival (H) in the entire TCGA-OV patient cohort.
The Prognostic Performance of the DEpiLncSig Is Independent of Other Clinical Features
To test whether the prognostic performance of the DEpiLncSig is independent of other clinical features, we conducted univariate and multivariate Cox regression analyses for the individual clinical variables (age, stage, grade, and treatment response) with and without the EpiLncRNASig. As shown in Table 2, the EpiLncRNASig and treatment response were all significantly associated with overall survival in the univariate and multivariate analysis in the discovery cohort. In the validation cohort, although the EpiLncRNASig (HR = 1.56, 95% CI = 1.07–2.28, p = 0.022), age (HR = 1.59, 95% CI = 1.09–2.31, p = 0.015) and treatment response (HR = 0.2, 95% CI = 0.13–0.33, p = 4.50E-11) were all significantly associated with overall survival in the univariate analysis, only the EpiLncRNASig (HR = 1.71, 95% CI = 1.09–2.67, p = 0.019) and treatment response (HR = 0.19, 95% CI = 0.12–0.31, p = 3.80E-11) were significant in the multivariate analysis. In the entire TCGA-OV cohort, univariate analysis showed that the EpiLncRNASig (HR = 2.09, 95% CI = 1.6–2.74, p = 6.80E-08), age (HR = 1.35, 95% CI = 1.04–1.75, p = 0.022), and treatment response (HR = 0.23, 95% CI = 0.17–0.32, p = 1.70E-18) were significantly associated with overall survival in the univariate analysis. However, the EpiLncRNASig (HR = 2.21, 95% CI = 1.61–3.05, p = 1.20E-06) still maintained a significant association with overall survival when after adjusted by other clinical features (Table 2). Results from multivariate Cox regression analyses suggested that the prognostic performance of the DEpiLncSig is independent of other clinical features and treatment response.
Table 2. Univariate and Multivariate Cox regression analysis with overall survival in different patient cohorts.
Association of the EpiLncRNASig With Chemotherapy Response
To further examine the association of the EpiLncRNASig with a chemotherapy response, we first compared the distribution of the risk score of the EpiLncRNASig for patients with complete response (CR) and non-CR, and we found that patients with CR have significantly lower risk scores than those with non-CR (Wilcoxon rank-sum test p = 0.0011) (Figure 4A). We then assessed the relationship between the EpiLncRNASig and the likelihood of CR by plotting the percentage of OV patients achieving CR as a function of the risk score and found that there was a significant association between risk score and the likelihood of CR (Pearson correlation coefficient r = −0.77, p = 0.0085) (Figure 4B). As shown in Figure 4C, patients with low EpiLncRNASig seem to have a higher likelihood of CR than those with high EpiLncRNASig. To further examine whether the EpiLncRNASig is predictive for both patients achieving CR and non-CR, we performed a stratification analysis. The results of the stratification analysis showed that the EpiLncRNASig could further subdivide patients achieving CR into the high-risk group and low-risk group with significantly different survival times (log-rank p < 0.0001) (Figure 4D). Similar results were observed for patients achieving non-CR (log-rank p = 0.011) (Figure 4E).
Figure 4. Association of the EpiLncRNASig with chemotherapy response. (A) Boxplots of EpiLncRNASig scores for patients with CR and non-CR. (B) Correlation of the EpiLncRNASig with complete response. (C) The percentage of patients achieving CR in the high-risk group and low-risk group. Kaplan-Meier survival curves of overall survival times between the high-risk group and low-risk group for patients with CR (D) and non-CR (E).
In silico Analysis for lncRNA Biomarker Function
To infer potential biological roles of these epigenetically regulated lncRNA biomarkers, we calculated the Pearson correlation coefficient between expression levels of lncRNAs and mRNAs, and selected top 100 mRNAs were considered as lncRNA-related mRNAs (Supplementary Table 2). Then we performed GO and KEGG function enrichment analysis for 471 lncRNA-related mRNAs. Results of GO analysis showed that lncRNA-related mRNAs were enriched in the extracellular matrix organization, connective tissue development, lymphocyte differentiation, and T cell differentiation and activation (Figure 5A). Results of KEGG analysis showed that lncRNA-related mRNAs were enriched in immune-related or cancer-related biological pathways, including primary immunodeficiency, cell adhesion molecules, proteoglycans in cancer, cytoline-cytokine receptor interaction and the T cell receptor signaling pathway, Th17 cell differentiation, Th1, and Th2 cell differentiation (Figure 5B).
Figure 5. Functional prediction of the EpiLncRNASig. (A) Functional enrichment analysis for GO terms. (B) Functional enrichment analysis for KEGG pathways.
Discussion
Accumulating evidence has shown that altered molecular profiles are an indispensable laboratory tool to improve cancer diagnosis, prognosis and therapeutic response (Malone et al., 2020), and overcome the limitations of typical clinical and imaging characteristics owing to the genetic and molecular heterogeneity of cancers. During the past years, lncRNAs have become a rising star in the field of biomarker research and were extensively studied and characterized in various cancers. Aberrant lncRNA expressions have also been observed in OV development, progression, recurrence, metastasis, and therapeutic response, indicating the potential roles as biomarkers for precision medicine of OV. Several lncRNA-focused expression signatures have been proposed for the prediction and monitoring of disease status, prognosis, and chemotherapeutic sensitivity. For example, Pan et al. (2020) proposed two lncRNA signatures to predict the prognosis and efficiency of chemotherapy. Zhou et al. (2016b) constructed progression-related lncRNA-associated ceRNA networks based on the “ceRNA hypothesis” and subsequently identified a 10-lncRNA signature associated with the outcome of OV patients. However, further investigation is needed to identify epigenetically regulated lncRNAs and explore their potential function and clinical application in OV.
For the above purpose, we performed an integrative analysis of DNA methylation data and transcriptome data and identified 419 lncRNAs as potential epigenetically regulated lncRNAs. These predicted epigenetically regulated lncRNAs in this study would provide an available resource for studying the interplay between lncRNAs and epigenetic regulation. A recent study performed by Li et al. (2019) has characterized the crosstalk between DNA methylation and lncRNA expression in thyroid Cancer (THCA) and methylation-driven 5-lncRNA-based signature (EpiLncPM) with potential clinical application in predicting the prognosis of THCA. Therefore, to further investigate whether epigenetically regulated lncRNAs identified here also have potential clinical value, we performed univariate Cox regression analysis with OS and found that only 17 of 419 epigenetically regulated lncRNAs identified here were significantly associated with patients’ outcome. By using machine-learning and multivariate Cox regression analysis methods, we identified and developed an epigenetically regulated lncRNA expression signature (EpiLncRNASig) consisting of five lncRNAs from the list of 17 epigenetically regulated lncRNAs associated with outcome. Furthermore, we tested and validated these in different patient cohorts, demonstrating a similar effective predictive performance in predicting OS. Moreover, multivariate Cox regression analyses after adjusted by other clinical features and treatment response demonstrated the independence of the DEpiLncSig in predicting OS. Further examination of the association of the EpiLncRNASig with chemotherapy response provided evidence supporting that the DEpiLncSig is not only a prognostic factor but may also be an indicator for chemotherapy response.
Compared to the huge number of lncRNAs identified and recorded in the public database, the number of functionally well-characterized lncRNAs is relatively small. Among five lncRNAs in the EpiLncRNASig, several lncRNAs have been well studied in some cancers. For example, LINC00189 has been reported to be associated with several cancers, including ovarian cancer, cervical cancer (Zhang et al., 2020), clear cell renal cell carcinoma (Xu et al., 2020), and urinary bladder cancer (Zhang et al., 2016). By analyzing CACNA1G-AS1 expression levels in 122 pairs of non-small cell lung cancer (NSCLC) and normal tissue samples as well as in NSCLC cell lines, Yu et al. (2018) identified CACNA1G-AS1 as an oncogene to promote cell migration, invasion, and epithelial-mesenchymal transition via regulating HNRNPA2B1. A recent study performed by Wei revealed higher expression of CACNA1G-AS1 in colorectal cancer (CRC) tissues compared to adjacent normal tissues and found that CACNA1G-AS1 was able to enhance proliferative and invasive abilities of CRC cells by downregulating p53 via forming a carcinogenic complex with EZH2 (Wei et al., 2020). Another recent study has reported that lncRNA CACNA1G-AS1 can act as competing endogenous RNAs (ceRNAs) to regulate miR-205 expression, and it promotes proliferation and invasion in human keloid fibroblasts (Zhao X. et al., 2020). To further gain novel insights into the functional roles of the EpiLncRNASig in OV, we performed functional enrichment analysis for relevant protein-coding genes of epigenetically regulated lncRNAs by considering their co-expression relation relationships and found that the EpiLncRNASig was involved in known immune-related or cancer-related biological pathways.
In conclusion, we reported and provided a knowledge base of novel epigenetically regulated lncRNAs in OV, which will improve our understanding of lncRNA epigenetic regulation mechanisms in OV. Finally, we proposed an epigenetically regulated five-lncRNA signature as promising prognostic biomarkers for predicting outcome and chemotherapy response with further validation in prospective cohort studies.
Data Availability Statement
Publicly available datasets were analyzed in this study. This data can be found here: DNA methylation and RNA-seq data of OV patients were obtained from the UCSC Xena Browser (https://xena.ucsc.edu/) under the cohort accession: GDC TCGA Ovarian Cancer (OV).
Author Contributions
XL and DW conceived and designed the experiments. HY, LG, MZ, NN, and YW performed the experiments and analyzed the data. DW and XL wrote the manuscript. All authors read and approved the final manuscript.
Funding
This study was supported by the National Natural Science Foundation of China (81860534 and 81760552), the Inner Mongolia Autonomous Region Science and Technology Planning Project (2019GG039 and 2019GG086), the Haiyan Fund of The Affiliated Tumor Hospital of Harbin Medical University (JJZD2020-14), and the Scientific Research Innovation Fund of The First Affiliated Hospital of Harbin Medical University (2020L05).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2021.644940/full#supplementary-material
Footnotes
References
Adib, T. R., Henderson, S., Perrett, C., Hewitt, D., Bourmpoulia, D., Ledermann, J., et al. (2004). Predicting biomarkers for ovarian cancer using gene-expression microarrays. Br. J. Cancer 90, 686–692. doi: 10.1038/sj.bjc.6601603
Bao, S., Hu, T., Liu, J., Su, J., Sun, J., Ming, Y., et al. (2021). Genomic instability-derived plasma extracellular vesicle-microRNA signature as a minimally invasive predictor of risk and unfavorable prognosis in breast cancer. J. Nanobiotechnol. 19:22. doi: 10.1186/s12951-020-00767-3
Bao, S., Zhao, H., Yuan, J., Fan, D., Zhang, Z., Su, J., et al. (2020). Computational identification of mutator-derived lncRNA signatures of genome instability for improving the clinical outcome of cancers: a case study in breast cancer. Brief. Bioinform. 21, 1742–1755. doi: 10.1093/bib/bbz118
Barrett, C. L., DeBoever, C., Jepsen, K., Saenz, C. C., Carson, D. A., and Frazer, K. A. (2015). Systematic transcriptome analysis reveals tumor-specific isoforms for ovarian cancer diagnosis and therapy. Proc. Natl. Acad. Sci. U.S.A. 112, E3050–E3057. doi: 10.1073/pnas.1508057112
Cancer Genome Atlas Research Network (2011). Integrated genomic analyses of ovarian carcinoma. Nature 474, 609–615. doi: 10.1038/nature10166
Clifford, C., Vitkin, N., Nersesian, S., Reid-Schachter, G., Francis, J. A., and Koti, M. (2018). Multi-omics in high-grade serous ovarian cancer: biomarkers from genome to the immunome. Am. J. Reprod. Immunol. 80:e12975. doi: 10.1111/aji.12975
Coburn, S. B., Bray, F., Sherman, M. E., and Trabert, B. (2017). International patterns and trends in ovarian cancer incidence, overall and by histologic subtype. Int. J. Cancer 140, 2451–2460. doi: 10.1002/ijc.30676
Hosseini, E. S., Meryet-Figuiere, M., Sabzalipoor, H., Kashani, H. H., Nikzad, H., and Asemi, Z. (2017). Dysregulated expression of long noncoding RNAs in gynecologic cancers. Mol. Cancer 16:107. doi: 10.1186/s12943-017-0671-2
Huo, X., Han, S., Wu, G., Latchoumanin, O., Zhou, G., Hebbard, L., et al. (2017). Dysregulated long noncoding RNAs (lncRNAs) in hepatocellular carcinoma: implications for tumorigenesis, disease progression, and liver cancer stem cells. Mol. Cancer 16:165. doi: 10.1186/s12943-017-0734-4
Li, J. Y., Li, C. J., Lin, L. T., and Tsui, K. H. (2020). Multi-omics analysis identifying key biomarkers in ovarian cancer. Cancer Control 27:1073274820976671. doi: 10.1177/1073274820976671
Li, Q., Wang, P., Sun, C., Wang, C., and Sun, Y. (2019). Integrative analysis of methylation and transcriptome identified epigenetically regulated lncRNAs with prognostic relevance for thyroid cancer. Front. Bioeng. Biotechnol. 7:439. doi: 10.3389/fbioe.2019.00439
Lu, J., Wu, D., Li, C., Zhou, M., and Hao, D. (2014). Correlation between gene expression and mutator phenotype predicts homologous recombination deficiency and outcome in ovarian cancer. J. Mol. Med. 92, 1159–1168. doi: 10.1007/s00109-014-1191-9
Malone, E. R., Oliva, M., Sabatini, P. J. B., Stockley, T. L., and Siu, L. L. (2020). Molecular profiling for precision cancer therapies. Genome Med. 12:8. doi: 10.1186/s13073-019-0703-1
Marchese, F. P., Raimondi, I., and Huarte, M. (2017). The multidimensional mechanisms of long noncoding RNA function. Genome Biol. 18:206. doi: 10.1186/s13059-017-1348-2
Mercer, T. R., Dinger, M. E., and Mattick, J. S. (2009). Long non-coding RNAs: insights into functions. Nat. Rev. Genet. 10, 155–159. doi: 10.1038/nrg2521
Momenimovahed, Z., Tiznobaik, A., Taheri, S., and Salehiniya, H. (2019). Ovarian cancer in the world: epidemiology and risk factors. Int. J. Womens Health 11, 287–299. doi: 10.2147/IJWH.S197604
Pan, Y., Jia, L. P., Liu, Y., Han, Y., Li, Q., Zou, Q., et al. (2020). A novel signature of two long non-coding RNAs in BRCA mutant ovarian cancer to predict prognosis and efficiency of chemotherapy. J. Ovarian Res. 13:112. doi: 10.1186/s13048-020-00712-w
Sallinen, H., Janhonen, S., Polonen, P., Niskanen, H., Liu, O. H., Kivela, A., et al. (2019). Comparative transcriptome analysis of matched primary and distant metastatic ovarian carcinoma. BMC Cancer 19:1121. doi: 10.1186/s12885-019-6339-0
Siegel, R. L., Miller, K. D., and Jemal, A. (2020). Cancer statistics, 2020. CA Cancer J. Clin. 70, 7–30. doi: 10.3322/caac.21590
Sun, J., Zhang, Z., Bao, S., Yan, C., Hou, P., Wu, N., et al. (2020). Identification of tumor immune infiltration-associated lncRNAs for improving prognosis and immunotherapy response of patients with non-small cell lung cancer. J Immunother Cancer 8:e000110. doi: 10.1136/jitc-2019-000110
Tian, T., Wang, M., Lin, S., Guo, Y., Dai, Z., Liu, K., et al. (2018). The impact of lncRNA dysregulation on clinicopathology and survival of breast cancer: a systematic review and meta-analysis. Mol. Ther. Nucleic Acids 12, 359–369. doi: 10.1016/j.omtn.2018.05.018
Wang, J., Dean, D. C., Hornicek, F. J., Shi, H., and Duan, Z. (2019). RNA sequencing (RNA-Seq) and its application in ovarian cancer. Gynecol. Oncol. 152, 194–201. doi: 10.1016/j.ygyno.2018.10.002
Wang, Z., Yang, B., Zhang, M., Guo, W., Wu, Z., Wang, Y., et al. (2018). lncRNA epigenetic landscape analysis identifies EPIC1 as an oncogenic lncRNA that interacts with MYC and promotes cell-cycle progression in cancer. Cancer Cell 33, 706.e709–720.e709. doi: 10.1016/j.ccell.2018.03.006
Wei, L. J., Bai, D. M., Wang, Z. Y., and Liu, B. C. (2020). Upregulated lncRNA CACNA1G-AS1 aggravates the progression of colorectal cancer by downregulating p53. Eur. Rev. Med. Pharmacol. Sci. 24, 130–136. doi: 10.26355/eurrev_202001_19902
Xu, D., Dang, W., Wang, S., Hu, B., Yin, L., and Guan, B. (2020). An optimal prognostic model based on gene expression for clear cell renal cell carcinoma. Oncol. Lett. 20, 2420–2434. doi: 10.3892/ol.2020.11780
Yao, R. W., Wang, Y., and Chen, L. L. (2019). Cellular functions of long noncoding RNAs. Nat. Cell Biol. 21, 542–551. doi: 10.1038/s41556-019-0311-8
Yu, G., Wang, L. G., Han, Y., and He, Q. Y. (2012). clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS 16, 284–287. doi: 10.1089/omi.2011.0118
Yu, P. F., Kang, A. R., Jing, L. J., and Wang, Y. M. (2018). Long non-coding RNA CACNA1G-AS1 promotes cell migration, invasion and epithelial-mesenchymal transition by HNRNPA2B1 in non-small cell lung cancer. Eur. Rev. Med. Pharmacol. Sci. 22, 993–1002. doi: 10.26355/eurrev_201802_14381
Zhang, X., Zhang, M., Hou, Y., Xu, L., Li, W., Zou, Z., et al. (2016). Single-cell analyses of transcriptional heterogeneity in squamous cell carcinoma of urinary bladder. Oncotarget 7, 66069–66076. doi: 10.18632/oncotarget.11803
Zhang, Y., Zhang, X., Zhu, H., Liu, Y., Cao, J., Li, D., et al. (2020). Identification of potential prognostic long non-coding RNA biomarkers for predicting recurrence in patients with cervical cancer. Cancer Manag. Res. 12, 719–730. doi: 10.2147/CMAR.S231796
Zhao, H., Gu, S., Bao, S., Yan, C., Zhang, Z., Hou, P., et al. (2020). Mechanistically derived patient-level framework for precision medicine identifies a personalized immune prognostic signature in high-grade serous ovarian cancer. Brief Bioinform. doi: 10.1093/bib/bbaa069 Online ahead of print
Zhao, X., Jie, X., Gao, Y. K., Nie, B., and Jiang, H. (2020). Long non-coding RNA CACNA1G-AS1 promotes proliferation and invasion and inhibits apoptosis by regulating expression of miR-205 in human keloid fibroblasts. Biosci. Rep. 40:BSR20192839. doi: 10.1042/BSR20192839
Zhao, T., Xu, J., Liu, L., Bai, J., Wang, L., Xiao, Y., et al. (2015). Computational identification of epigenetically regulated lncRNAs and their associated genes based on integrating genomic data. FEBS Lett. 589, 521–531. doi: 10.1016/j.febslet.2015.01.013
Zhou, M., Sun, Y., Sun, Y., Xu, W., Zhang, Z., Zhao, H., et al. (2016a). Comprehensive analysis of lncRNA expression profiles reveals a novel lncRNA signature to discriminate nonequivalent outcomes in patients with ovarian cancer. Oncotarget 7, 32433–32448. doi: 10.18632/oncotarget.8653
Zhou, M., Wang, X., Shi, H., Cheng, L., Wang, Z., Zhao, H., et al. (2016b). Characterization of long non-coding RNA-associated ceRNA network to reveal potential prognostic lncRNA biomarkers in human ovarian cancer. Oncotarget 7, 12598–12611. doi: 10.18632/oncotarget.7181
Keywords: ovarian cancer, long non-coding RNAs, epigenetic modification, chemotherapy response, signature
Citation: Yang H, Gao L, Zhang M, Ning N, Wang Y, Wu D and Li X (2021) Identification and Analysis of An Epigenetically Regulated Five-lncRNA Signature Associated With Outcome and Chemotherapy Response in Ovarian Cancer. Front. Cell Dev. Biol. 9:644940. doi: 10.3389/fcell.2021.644940
Received: 22 December 2020; Accepted: 03 February 2021;
Published: 23 February 2021.
Edited by:
Jianzhong Su, Wenzhou Medical University, ChinaReviewed by:
Hao Cai, First Affiliated Hospital of Gannan Medical University, ChinaShengli Li, Shanghai Jiao Tong University, China
Copyright © 2021 Yang, Gao, Zhang, Ning, Wang, Wu and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Di Wu, bWFzdGVyd3U3MjhAMTYzLmNvbQ==; Xiaomei Li, ZmFubGl3ZW5xaUAxNjMuY29t
†These authors have contributed equally to this work and share first authorship