Next Article in Journal
Levels of Evidence Supporting United States Guidelines in Pancreatic Adenocarcinoma Treatment
Next Article in Special Issue
Single-Cell Profiling of the Immune Atlas of Tumor-Infiltrating Lymphocytes in Endometrial Carcinoma
Previous Article in Journal
Defective DNAM-1 Dependent Cytotoxicity in Hepatocellular Carcinoma-Infiltrating NK Cells
Previous Article in Special Issue
Identification of a Novel Cuproptosis-Related Gene Signature for Prognostic Implication in Head and Neck Squamous Carcinomas
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cancers
Volume 14
Issue 16
10.3390/cancers14164056
cancers-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The m6A-Related Long Noncoding RNA Signature Predicts Prognosis and Indicates Tumor Immune Infiltration in Ovarian Cancer

by
Rui Geng
1,†,
Tian Chen
2,†,
Zihang Zhong
1,
Senmiao Ni
1,
Jianling Bai
1,* and
Jinhui Liu
2,*
1
Department of Biostatistics, School of Public Heath, Nanjing Medical University, 101 Longmian Avenue, Jiangning District, Nanjing 211166, China
2
Department of Gynecology, The First Affiliated Hospital of Nanjing Medical University, Nanjing 210029, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Cancers 2022, 14(16), 4056; https://doi.org/10.3390/cancers14164056
Submission received: 5 July 2022 / Revised: 14 August 2022 / Accepted: 17 August 2022 / Published: 22 August 2022
(This article belongs to the Collection Application of Bioinformatics in Cancers)
Browse Figures
Review Reports Versions Notes

Abstract

:

Simple Summary

OV is the most lethal gynecological malignancy. M6A and lncRNAs have a great impact on OV development and patient immunotherapy response. This study provided an accurate prognostic signature for patients with OV and elucidated the potential mechanism of the mRLs in immune modulation and treatment response, giving new insights into identifying new therapeutic targets.

Abstract

Background: OV is the most lethal gynecological malignancy. M6A and lncRNAs have a great impact on OV development and patient immunotherapy response. In this paper, we decided to establish a reliable signature of mRLs. Method: The lncRNAs associated with m6A in OV were analyzed and obtained by co-expression analysis of the TCGA-OV database. Univariate, LASSO and multivariate Cox regression analyses were employed to establish the model of mRLs. K-M analysis, PCA, GSEA and nomogram based on the TCGA-OV and GEO database were conducted to prove the predictive value and independence of the model. The underlying relationship between the model and TME and cancer stemness properties were further investigated through immune feature comparison, consensus clustering analysis and pan-cancer analysis. Results: A prognostic signature comprising four mRLs, WAC-AS1, LINC00997, DNM3OS and FOXN3-AS1, was constructed and verified for OV according to the TCGA and GEO database. The expressions of the four mRLs were confirmed by qRT-PCR in clinical samples. Applying this signature, one can identify patients more effectively. The samples were divided into two clusters, and the clusters had different overall survival rates, clinical features and tumor microenvironments. Finally, pan-cancer analysis further demonstrated that the four mRLs were significantly related to immune infiltration, TME and cancer stemness properties in various cancer types. Conclusions: This study provided an accurate prognostic signature for patients with OV and elucidated the potential mechanism of the mRLs in immune modulation and treatment response, giving new insights into identifying new therapeutic targets.

1. Introduction

Ovarian cancer is a gynecological malignancy which is prevalent in women aged 40 to 79 [1]. Additionally, OV is a common type of this cancer, and the occurrence of OV is ninety percent in all patients with ovarian cancer [2]. Due to the lack of specific initial symptoms and sensitive biomarkers for early diagnosis, most OV patients cannot be diagnosed at an early stage and the progression is rapid [3,4]. Consequently, developing novel and reliable signatures to diagnose and prognose OV at an early stage is an urgent need.
The results of randomized trials show that transvaginal ultrasound and CA-125 testing are of little help to diagnose ovarian cancer and change the mortality rate [5].
Therefore, more attention has been given to discovering key regulatory genes involved in cancer pathogenesis and progression and constructing prognostic signatures based on them. Recently, as a reversible epigenetic modification of various RNAs, the vital role and molecular mechanisms of m6A RNA modification in cancer pathogenesis, drug response and developing new targets for therapy have been intensely studied and confirmed [6]. Furthermore, it has also been found that m6A RNA modification functions in immunity, which provides insights into cancer immunotherapy [7]. M6A is dynamically regulated by three types of proteases with distinct functions, namely “writers, “erasers” or “readers” [8,9]. Changing the structure of RNA can influence a variety of cell processes. As a result, the effect of mRLs may play an important role in the transfer of cancer cells.
LncRNA is a group of RNA molecules. Although lncRNA cannot be translated into protein, it has a vital effect on many biological processes such as tumorigenesis and development [10,11,12]. The relationship between m6A and lncRNAs is under in-depth study. Extensive studies have revealed that the biogenesis and functions of lncRNAs depend on m6A modification [13,14,15]. YTHDF1 and YTHDF2, m6A readers, modify the stability of oncogenic lncRNA THOR, thus promoting the biological behavior of cancer cells [16]. LncRNA RP11 dependent on m6A induction can trigger colorectal cancer metastasis by preventing the degradation of Zeb1 [17]. In addition, m6A modification is also adversely affected by lncRNAs [18], and they can function together to regulate protein translation [18,19]. For instance, LINRIS stabilize IGF2BP2 (an m6A reader) by preventing its K139 ubiquitination, thus promoting colorectal cancer cell line growth [20]. In OV, the stability of RHPN1-AS1 is improved by m6A modification, promoting the proliferation and metastasis of OV cell [21]. Hence, linking m6A and lncRNA, and further studying the function and molecular mechanism of mRLs in tumorigenesis and progression can affect the overall survival rate of OV. However, such research on OV is still lacking.
This study established and verified a mRL prognostic signature for OV using the TCGA and GEO database. The underlying correlation between the signature and TME, as well as therapeutic response, was also explored. Moreover, clustering subgroups and pan-cancer analysis were also used to verify the application value of the model in distinguishing patients in terms of prognosis and therapeutic response.

2. Materials and Methods

2.1. Collecting and Disposing Data

The flow chart of our study is performed in Figure 1a. We collected gene expression profiling and clinicopathological data of OV samples from the TCGA dataset (https://cancergenome.nih.gov/, (accessed on 5 February 2022)) and the GSE9891 database (http://www.ncbi.nlm.nih.gov/geo/, (accessed on 5 February 2022)). The single gene expression in OV samples with different stages was obtained from GEPIA (http://gepia.cancer-pku.cn/index.html, (accessed on 5 February 2022)) [22]. To identify the role of the four mRLs in more cancers, we obtained the pan-cancer data of 33 cancers from the UCSC Xena database (https://xenabrowser.net/datapages/, (accessed on 5 February 2022)) [23], including RNA-seq, immune subtypes, prognosis profiles and stemness scores.

2.2. Identification of mRLs Associated with OV Prognosis

Firstly, the expression matrix of 23 m6A regulators (METTL3, METTL14, METTL16, WTAP, VIRMA, ZC3H13, RBM15, RBM15B, YTHDC1, YTHDC2, YTHDF1, YTHDF2, YTHDF3, HNRNPC, FMR1, LRPPRC, HNRNPA2B1, IGFBP1, IGFBP2, IGFBP3, RBMX, FTO and ALKBH5) was achieved in previous studies [24]. Secondly, Pearson correlation analysis was implemented to select 419 mRLs according to the standard of p < 0.05 and absolute correlation coefficient >0.3, then 61 mRLs shared in both TCGA and GEO were obtained after the removal of batch effects through combat via the “sva” R package. Finally, we collected 374 samples with sufficient prognostic clinical data and lncRNAs expression information from TCGA. The 374 samples from the TCGA-OV dataset were assigned into TCGA-train (n = 187) and TCGA-test (n = 187) according to 1:1 randomly, and we employed the train set to build the signature. The test set (187 samples), the entire set (374 samples) and the GSE9891 database (278 samples) were applied for validation and further research. Univariate Cox regression analysis was applied to screen the relationships between mRLs related to the patients’ prognoses in the TCGA-train set using the survival package of R (https://www.bioconductor.org/packages/release/bioc/html/rbsurv.html, (accessed on 5 February 2022)).

2.3. RNA Isolation and qRT-PCR

The research was approval by the First Affiliated Hospital of Nanjing Medical University Ethics Committee and the participants who took part in the study all provided written informed consent. Total RNA of the OV and normal tissues was obtained through TRIzol reagent (Thermo Fisher Scientific, Waltham, MA, USA), whose integrity was estimated by the Agilent Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA, USA). cDNA synthesis was conducted using the high-capacity reverse transcription kits (TaKaRa, Shiga, Japan) and then qRT-PCR was based on SYBR Green PCR Kit (Thermo Fisher Scientific, Waltham, MA, USA) and the 2−ΔΔCt method on Light Cycler 480 (Roche, Switzerland). GAPDH was the endogenous control. Primer sequences for GAPDH and four mRLs were presented in Table S1.

2.4. Establish and Prove an mRL Prognostic Signature for OV

The model was established using LASSO in TCGA-train. Then, the multivariate regression analysis was conducted to determine mRLs with independence. The format to assess the risk score was: risk   score   = ∑ i = 1 n coef i ∗ lncRNA i expression . The patients were divided into two groups in terms of median risk score of train set. Kaplan–Meier analysis of the OS of the groups was conducted through the “survival” R package. Through the “ROC” R package, the ROC was depicted and the AUC was calculated to determine its specificity and sensitivity. PCA was used to prove the grouping ability of the model by grouping the visualization of high-dimensional data of the risk model [25]. We also conducted subgroup analysis to verify the prognostic ability of the mRL model. Univariate and multivariate Cox analysis were applied to assess whether the signature was an independent factor.

2.5. Establishment and Validation of a Nomogram

To improve the predictive ability of the mRL signature, we established a nomogram composed of the signature, patient age, grade and tumor stage. Calibration curves were applied to explore the accuracy and reliability of the nomogram.

2.6. Gene Set Enrichment Analysis

As carried out in our previous study [26], we performed function annotation through GSEA to reveal potential mechanisms in different risk groups. KEGG in GSEA was applied to select predefined gene sets; 5000 permutations were conducted for the gene set to calculate p-values. GSEA was analyzed using “clusterProfiler” package [27].

2.7. The Correlation between the Signature and Tumor Infiltration Immune Cells (TIICs)

CIBERSORT was employed to assess the abundance of TIIC profiles in all tumor tissues. We assessed the proportion of 22 types of TIICs in every patient according to the CIBERSORT score [28]. Furthermore, the correlation between the lncRNAs and TIICs score was evaluated to identify which immune cells had distinguished between cancers.

2.8. The Correlation to Other Immune Features

ESTIMATE was used to compute the proportion of immune and stromal components in TME for each sample [29]. The relationships between lncRNA expression and the risk score were calculated by Spearman correlation. Furthermore, we identified six immune types and used the analysis of variance (ANOVA) to evaluate the relationships between these subtypes and lncRNA expression. The data of cancer stem cell-like properties were employed to assess the stemness characteristics of the cancer cells. The association between stemness features and lncRNA was also evaluated by Spearman analysis.

2.9. Immunotherapy Response Prediction

Immunophenoscores (IPSs) are calculated based on four major categories of gene expression z-scores to evaluate and compare the potential reaction to ICI between the two groups, and the high scores represent high immunogenicity [30]. IPS data of each OV patient in TCGA-entire were achieved from The Cancer Immunome Atlas (TCIA) (https://tcia.at/home, (accessed on 5 February 2022)).

2.10. Assess Drug Sensitivity

The newest CellMiner version (https://discover.nci.nih.gov/cellminer/, (accessed on 5 February 2022)) could help researchers achieve the NCI-60 data flexibly for correlations between genomic, molecular and pharmacologic parameters [31]. We collected the expression of 4 m6A-related genes and z-score for cell sensitivity data (GI50) from the webpage and used Pearson’s correlation coefficient to assess the influence of the m6A-related gene expression on drug sensitivity.

2.11. Consensus Clustering

The expression of the lncRNAs was utilized to recognize the subtypes in OV patients via “ConsensusClusterPlus” R package (http://www.bioconductor.org/, (accessed on 5 February 2022)). We used the log-rank test and K–M curve to calculate the OS difference between groups.

2.12. Statistical Analysis

R 4.0.2 was used to carry out statistical analyses. The relationships between the risk score or the mRL expression level and stemness score, estimate score and drug sensitivity were calculated using Spearman or Pearson correlation analysis. Univariate or multivariate regression analyses were applied to calculate the relationships between the four mRLs with patient OS. Kaplan-Meier analysis could compare the survival state between the groups. Time-dependent ROC curve analysis could test the predict ability of the model. Subgroup analysis could evaluate the stability of the prognostic signature in subgroups stratified by clinical characteristics. Student’s t-test and ANOVA were applied to find the differences between the different groups. Linear mixed-effect models were used for pan-cancer analysis. The hazard ratio and 95% confidence interval were assessed to screen prognosis-related genes. p < 0.05 was regarded as statistically significant.

3. Results

3.1. Identification of mRLs in OV

Firstly, we abstracted the expression of lncRNAs from the TCGA-OV dataset and identified 419 mRLs (p < 0.05, absolute correlation coefficient > 0.3) through Pearson’s correlation analysis with m6A regulators. The relational network between the m6A regulators and mRLs was shown in Figure 1b. After collecting the expression of lncRNAs in the GSE9891 dataset, we found 61 mRLs shared in both datasets. When combined with clinical information, we screened five mRLs related to prognosis through the univariate regression analysis of the TCGA-train (Table S2).

3.2. Establish and Verify a Prognostic Signature Based on mRLs in OV Patients

Firstly, 374 OV patients were randomly classified into train and test sets at a 1:1 ratio, including 187 OV patients, respectively. Secondly, we performed the LASSO (Figure S1a,b) and multivariate regression analysis (Figure S1c) in the train set to build a reliable prognostic signature. The formula was as follows:
Risk core = (−0.068088177 * expr (WAC-AS1)) + (−0.276777737 * expr (LINC00997)) + (0.09155959 * expr (DNM3OS)) + (−0.137781856 * (FOXN3-AS1), expr means the expression value. WAC-AS1, LINC00997, DNM3OS and FOXN3-AS1 were finally identified as key mRLs. Both the outcomes of the expression analyses in the GEPIA dataset and the qRT-PCR results indicated that the expression of DNM3OS and FOXN3-AS1 was higher in normal tissue compared to the tumor tissue (p < 0.05) (Figure S2a–h). Afterward, we assessed the risk score of OV patients in the TCGA-train and separated them into two groups. The distributions of risk score, survival status of each sample (Figure 2a) and the heatmap of key gene expression patterns (Figure 2e) in the TCGA-train are shown in Figure 2. People in the high-risk group had poor survival prognosis compared to the low-risk groups in the TCGA-train (Figure 2i). Time-dependent receiver-operating characteristic curve (ROC) analysis based on the signature revealed that the AUC of 1-, 2- and 3-year survival was 0.624, 0.694 and 0.630, respectively (Figure 2m). To verify the mRL prognostic signature, we calculated the risk score of the patients in the TCGA-test (N = 187), TCGA-entire (N = 374) and GEO datasets (N = 278) and divided the patients into either group according to their median risk score of the TCGA-train (Figure 2b–d). The heatmap showed that, with the exception of DNM3OS, the expression of WAC-AS1, LINC0099 and FOXN3-AS1 was higher in the low-risk groups (Figure 2e–h). The Kaplan–Meier analyses of the TCGA-test, TCGA-entire and GEO datasets revealed the same trend, that the OV patients with a low risk had higher OS than the high-risk group (p < 0.01) (Figure 2j–l). Furthermore, the time-dependent ROC curve highlighted the capability of the prognostic signature based on mRLs. The AUC of 1-, 2- and 3-year survival in the three sets is shown in Figure 2n–p. In addition, the distributions of the groups are different in the PCA plots (Figure S3). These results all demonstrate the accuracy and reliability of the mRL prognostic signature in OV prediction.

3.3. The Relationships between the Model and Clinical Factors

The correlations between the risk core and clinical characteristics are performed in Figure S4a–f. The risk score in the age ≤ 60 group was lower than the age > 60 group (p < 0.05, Figure S4d), and the proportion of older patients with a high-risk score was 52%, while the people with a low-risk score was 40% (Figure S4a). However, the difference in the risk score in the groups assigned by tumor grade and tumor stage was not significant (Figure S4b–c,e–f). Furthermore, the risk score of the high immune score group was markedly higher than the lower score patients (p < 0.05, Figure S4g), and the significant difference could be found in the risk core of Cluster 1 and Cluster 2 subgroups (p < 0.05, Figure S4h). Stratification analysis grouped by patient age, tumor grade and tumor stage highlighted that the OS of the patients in the high-risk group was worse than the low-risk group (Figure S4i–n). When comparing the patients’ states between the groups, we also found a higher mortality of patients in the high-risk group (Figure S5a), and the risk scores of dead patients were statistically higher than the alive patients (p = 0.002, Figure S5b). K–M curves of OS for patients subjected to chemotherapy demonstrated that patients at a high risk had poorer prognoses than those at a lower risk (p < 0.001, Figure S5c). In patients with BRCA1, the low-risk group also had a favorable survival outcome (Figure S5d).

3.4. Clinical Application of the Signature

Univariate and multivariate regression analyses in the TCGA-train, TCGA-test, TCGA-entire and GEO datasets all indicated the independence of the mRL model (Table S3). Moreover, patient age and tumor stage were also crucial poor prognostic factors for OV. To predict performance of the signature, we built a nomogram comprising the risk score, age, grade and stage in both the TCGA (Figure 3a) and GEO datasets (Figure S6). Figure 3b–d revealed that the predicted rates of OS were highly consistent with the observed rates. In addition, the prognostic signature showed superior predictive ability compared to other clinical characteristics, and when combined with clinical characteristics, the model showed better predictive power than the signature when used alone (Figure 3e).

3.5. The Correlations between the mRL Prognostic Signature and TME

GSEA indicated that the pathway related to immunity was enriched in people with a high-risk score (Figure 4a,b). Firstly, we implemented the ESTIMATE algorithm which indicated that the stromal score, immune score and ESTIMATE score were higher in the high-risk group (p < 0.05), and the risk score had a positive association with them (Figure 4c–h). The same trends were found in the GEO dataset (Figure 4i–n). Then, by comparing the 22 TIICs between the groups, we identified that the resting dendritic cells and M2 macrophages were positively related to the risk score and (p < 0.05) and the infiltration level of the T follicular helper cells and regulatory T cells (Tregs) had a negative correlation with the risk score (p < 0.05) in the TCGA dataset (Figure 5a–i). In the GEO dataset, the risk score had a positive correlation with gamma delta T cells and resting memory CD4 T cells, while the risk score has a negative relationship with Tregs, resting NK cells and activated dendritic cells (Figure S7). The differences in immune score and abundance of TIICs between the groups showed that patients with a high risk had a repressive immune phenotype, which could partly explain the poorer OS of people with a high risk score.
To identify the immunological role of the four mRLs, we assessed their association with the TIICs using the CIBERSORT algorithm. We found that DNM3OS had a significant association with T follicular helper cells, M2 macrophages and memory B cells (p < 0.001); WAC-AS1 was correlated with T follicular helper cells and M2 macrophages (p < 0.001) (Figure S8a). Moreover, WAC-AS1 (R = −0.15, p = 0.0041) and LINC00997 (R = −0.23, p < 0.001) had negative relationships with the immune score, whereas DNM3OS (R = 0.13, p = 0.008) and FOXNS-AS1 (R = 0.19, p < 0.001) were positively related to the immune score (Figure S8b). Next, to clarify the potential mechanisms of the lncRNAs in tumor progression, we used GSEA analysis to establish the enriched pathways of the mRLs in OV (Figure S8c–f). DNM3OS was highly enriched in the calcium signaling pathway, the focal adhesion pathway, the hematopoietic cell lineage pathway, the neuroactive ligand–receptor interaction pathway and cancer pathways.
The expression level of PD-L2 was higher in tumor tissues than in normal tissues (p < 0.05, Figure 6a), and it was markedly higher in people with a high risk score in both the TCGA and GEO dataset (p < 0.05, Figure 6c–f), suggesting that people in the high-risk group will benefit less from immune checkpoint inhibitors (ICIs). In addition, LINC00997 showed a high correlation with PD-L2 (Figure 6b); thus, the function of LINC00997 in the tumor immune microenvironment deserves further research. Comparing the IPSs between the groups, low-risk groups had significantly higher IPSs, suggesting a more immunogenic phenotype (Figure 6g–j). These results all demonstrate that the signature can predict the efficacy of immunotherapy for OV patients.

3.6. Consensus Clustering for mRLs Related to OV Prognosis and TME

After the consensus clustering, OV patients in the TCGA dataset were classified into Cluster 1 and Cluster 2 subgroups following the expression of the significant prognostic mRLs (Figure 7a–d). The heatmap of the expression pattern between Cluster 1 and Cluster 2 subgroups was performed in Figure 7e. In addition, the stroma, immune and ESTIMATE scores were markedly higher and the tumor purity was substantially lower (p < 0.05) in Cluster 2 (Figure 7g–j). Moreover, the infiltration of naïve B cells (p < 0.05), T cell activated memory CD4 cells (p < 0.05), T follicular helper cells (p < 0.05), M1 macrophages (p < 0.05) and resting mast cells (p < 0.05) was higher in the Cluster 1 subgroup, whereas the abundance of M0 macrophages (p < 0.05), M2 macrophages (p < 0.05), activated mast cells (p < 0.05) and neutrophils (p < 0.05) was higher in the Cluster 2 subgroup (Figure 7f); thus, the Cluster 2 subgroup tended to behave more as immunosuppressive phenotype. Moreover, the results of Kaplan–Meier analyses highlighted that the OS of Cluster 2 was dramatically worse than that of Cluster 1 (Figure 7k), and the expression of PD-L2 was higher in Cluster 2 (Figure 7l); this indicates that patients in the Cluster 2 may respond more sensitively to immunotherapy. The above results all demonstrate that the expression of the mRLs influences the TME, leading to different prognoses.

3.7. The Expression and Immune Status of lncRNAs in Pan-Cancer

We used the pan-cancer data of 33 cancers from TCGA to further study the relationships between the four mRLs and immune features, stem-like properties and patient prognoses in pan-cancer. The expression pattern of lncRNAs in 33 cancers (Figure 8a–d) and the heatmap (Figure S9a) showed inter-tumor heterogeneity, but in most cancers, the expression level of WAC-AS1 was higher in tumor tissues, whereas the expression level of DNM3OS was the opposite, which was consistent with the expression trend in OV (Figure S2). The expression distribution of the four mRLs across all 33 cancer types showed the same trend as the above results (Figure 8e). Moreover, there is no significant correlation between the key mRLs (Figure S9b). The Kaplan–Meier analyses (Figure S9c) and univariate regression (Figure S10a) analyses showed that the role of the mRLs in the survival of patients with different cancers varied.
To identify the relationships between the key mRLs and immune features in pan-cancer, we investigated the immune subtype and the landscape of correlation with stromal score, immune score, stemness scores based on DNA-methylation (DNAss) and stemness scores based on mRNA (RNAss). Through an extensive immunogenomic analysis of pan-cancer, six immune subtypes, including wound healing (C1), IFN-γ dominant (C2), inflammatory (C3), lymphocyte depleted (C4), immunologically quiet (C5) and TGF-β dominant (C6), were identified according to differences in immune, genetic and clinical features, in which C3 had great OS, C2 and C1 had poor prognosis, while C4 and C6 had the least favorable outcomes [32]. We evaluated the correlations between the mRL expression level and the six immune subtypes and found that all of them were related to pan-cancer immune subtypes and elevated C1, C2 and C6 subtypes with upregulated expression of DNM3OS (Figure 9a). Additionally, DNM3OS had a positive correlation with stromal score and immune score, whereas WAC-AS1 and LINC00997 were negatively related to stromal and immune score (Figure 9b–d), which showed similar correlation in OV. All four mRLs were negatively related to RNAss (Figure 9e). Therefore, the four mRLs took part in the formation of specific TIME. Finally, we applied the CellMiner database to evaluate the correlation between the mRLs and drug sensitivity. A higher z-score means higher sensitivity to one specific drug. Significantly, the expression of FOXNS-AS1 influenced drug sensitivity in various cell lines (Figure S10b).

4. Discussion

Due to the rapid progression and challenge of early diagnosis, people with OV face a poor OS, with five-year survival rates below 45% [33]. More and more efforts have been carried out to discover robust and sensitive predictive models for early diagnosis and prognosis prediction, including models based on immune genes, autophagy-related genes, m6A regulators and so on [34,35,36,37,38]. We found emerging evidence proving that m6A modification had a crucial effect in the formation of TME landscape heterogeneity and complexity [39,40,41]. Zhang et al. constructed an m6A score and demonstrated that a low m6A score was correlated with a high mutation burden, an inflamed TME phenotype and improved survival [42]. As the main modification for RNAs, m6A modifications almost regulate all biological pathways and their changes in tumors reprogram the tumor immune microenvironment, helping tumor cells escape and metastasize [15,43]. This study focused on mRLs, established a signature for OV patients and validated its predictive capability using TCGA and GEO datasets.
The model was established based on four mRLs associated with patient prognoses, including DNM3OS, WAC-AS1, FOXNS-AS1 and LINC00997, among which DNM3OS was correlated with a high-risk score, while the others were correlated with a low-risk score. It was proved that the overexpression of DNM3OS played a role in the development of ovarian cancer [44]. WAC-AS1, a competing endogenous RNA, has a great influence on the regulation of tumor glycolysis [45]. Additionally, the prognoses of glioma patients who had low WAC-AS1 expression were better than the high-expression group [46]. In contrast to our results, it was reported that LINC00997 had a positive relationship with the metastasis and development of colorectal cancer [47]. Subsequently, OV patients were separated into two groups in terms of the median risk score, and the patients in higher risk group had a worse prognosis. Multivariate Cox regression analysis highlighted that the model had a great effect on OS. ROC analysis proved that the model was more accurate than single clinical characteristics to predict the outcomes of OV. We constructed a nomogram for clinical application. GSEA revealed that people with a high-risk score were significantly enriched in immune response and immune system process pathways. We believe that the mRLs can take part in the mechanisms of tumorigenesis and progression of OV, and the signature might provide an accurate prediction and a theoretical foundation for the treatment of OV.
Due to the resistance to chemotherapy, patients with OV ultimately die of recurring lesions after surgical removal. To improve the prognoses of OV patients, new treatment modalities have emerged, especially immunotherapy, which could boost the patient’s immune system to find and kill tumor cells. The immune landscape of an OV patient determines their specific response to various treatments and subsequent prognosis. The proportion of intratumoral CD8+TILs, CD4+ Treg, macrophages and MDSC was significantly correlated with patients’ outcomes [48,49,50]. This finding has led to several stratification strategies, such as based on the number of CD8+TILs [50], six immune subtypes [32] or cytotoxic immunophenotypes [51] to identify those who benefit from immunotherapy. After validating the mRL prognostic signature in TCGA and GEO datasets, we performed immune subtype, consensus clustering and pan-cancer analyses to further research their function when regulating TME and the underlying molecular mechanism of interaction between lncRNA and m6A regulators in OV.
TIME is associated with cancer occurrence and progression [48,52]. This study calculated the correlation between the signature and TME in OV. According to the above results, high immune scores means higher risk scores. Some published articles showed that OV higher immune score patients had a poorer prognosis [53,54], consistent with our findings. After that, positive associations were identified between the risk score and resting dendritic cells, M2 macrophages, gamma delta T cells as well as resting memory CD4 T cells, while negative associations were identified between the risk score and T follicular helper cell infiltration, Tregs and resting NK cells. TIICs in the TME have been found to be related to the prognoses of OV patients [55]. For example, tumor-infiltrating macrophages with an M2 phenotype showed immunosuppressive activity, negatively related to OS [56]. It was also reported that the function of resting dendritic cells could be inhibited through GDF-15 interacting with CD44, thus facilitating ovarian cancer immune escape [57]. Therefore, the high level of M2 macrophages is closely correlated with a higher risk score and poorer prognosis, while resting NK cell infiltration improved the survival of patients with ovarian cancer [58,59]. Therefore, bad clinical prognosis of high-risk patients may be associated with TIICs, and mRLs may act as significant regulators of TIIC infiltration in OV.
Cancer cells can create an immunosuppressive microenvironment by activating immune checkpoint pathways, thus achieving immune escape [60]. Therefore, blocking the role of immune checkpoints in cancer treatment has attracted attention. Some ICIs have been approved for clinical application in several human malignancies [61]. Although no study has found an improvement in the survival of patients in clinical trials, including using ICIs targeting PD1, PD-L1 and CTLA4 in OV [62], identifying treatment targets, predicting treatment response, screening potential drugs and providing new immunotherapy will provide new breakthroughs for combined and personal therapy of OV [63]. The relationships between mRLs and immune checkpoint expression indicated that the risk score had a positive relationship with the PD-L2 expression level. The high expression of an immune checkpoint ligand PD-L2 represents a suboptimal response to ICIs [64]. The low-risk groups exhibited higher IPSs, indicating they may benefit more from immunotherapy, which is consistent with previous research that found that the anti-tumor responses mediated by ICIs depend on basic immunophenotypes [65,66,67,68]. Based on the above results, we propose that our risk signature may help predict immunotherapy efficacy for OV. Additionally, our study highlights the biological mechanism of mRLs in the occurrence and progression of OV.
Furthermore, we conducted consensus clustering analysis to screen m6A-related genes of OV and achieved two clusters according to the expression of four mRLs. We found that the expression of WAC-AS1, LINC00997 and FOXN3-AS1 was higher, and the expression of DNM3OS was lower in Cluster 1; thus, patients in Cluster 1 had lower risk scores. We compared the biological characteristics between the two clusters to further clarify the relationships between immune features and the mRLs. Cluster 1 showed higher infiltration of niave B cells, memory-activated CD4 T cells, T follicular helper cells and M1 macrophages than Cluster 2, while Cluster 2 had a higher infiltration of M2 macrophages than Cluster 1. The estimate score was higher in Cluster 2, and the tumor purity was higher in Cluster 1, which was consistent with the above results. The Kaplan–Meier analysis suggested that Cluster 2 had a poor prognosis. Furthermore, the expression of PD-L2 was higher in Cluster 2, which suggested that Cluster 1 respond more sensitively to ICIs. How m6A modification influences the tumor microenvironment, especially the TIIC types, has been researched [42,69,70]; our study might provide more evidence for further research.
Finally, we used pan-cancer analysis to further study the role of the mRLs in different cancers. The expression pattern of the four mRLs in pan-cancer showed intertumor heterogeneity, but in most cancer types, WAC-AS1 had a higher expression in tumor tissues, whereas the expression level of DNM3OS was the opposite, which was consistent with the expression trend in OV. By calculating the correlation between the lncRNAs and immune features, we found that DNM3OS was closely related to M2 macrophages (p < 0.001) in OV, elevated C1, C2 and C6 subtypes and an immune score with upregulated expression of DNM3OS in pan-cancer. Das et al. proved that the overexpression of DNM3OS promoted macrophage inflammatory phenotype, immune response genes and phagocytosis [71]. It was also reported that the function of DNM3OS was cooperated by miR-214 through the proinflammatory TLR4/IFN-γ/STAT1 pathways regulated pericystic macrophage accumulation [72]. These findings all demonstrate that mRLs play a role in TME and significantly influence OS.
Our study still has several limitations. Firstly, the results were entirely based on open access databases, limiting the validation of our cohort and external multicenter cohorts. Secondly, the interactions between mRLs and m6A regulators and the mechanism of the function of lncRNAs in TME were not demonstrated by the experiments. Thirdly, the prognostic signature needs to be applied in a real-world setting in order to explore its sensitivity and efficiency in classifying patients for target therapies.

5. Conclusions

In conclusion, our study established an independent and reliable prediction signature of mRLs, and systematically evaluated its predictive accuracy and role in TIME. Consensus clustering analysis and pan-cancer analysis demonstrated the vital role of the specific mRLs in different cancer types, which may help the discovery of novel therapeutic targets for OV.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cancers14164056/s1, Figure S1: Construction of a risk signature based on m6A-related prognostic lncRNAs in the TCGA-train set. LASSO regression and multivariate Cox regression were performed, calculating the minimum criteria (a,b) and coefficients (c); Figure S2: The expression of the four prognostic m6A-related lncRNAs in TCGA dataset and our own samples. (a–d) The expression of the four prognostic m6A-related lncRNAs in tumor (N = 426) and normal tissue (N = 88) in TCGA dataset. (e–h) The expression levels of the four prognostic m6A-related lncRNAs mRNA in normal (N = 16) and tumor (N = 16) clinical tissue samples. * p < 0.05, ** p < 0.01, and *** p < 0.001; Figure S3: Principal component analysis between the high- and low-risk groups based on the m6A-related lncRNA prognostic signature in the TCGA-train (a, N = 187), TCGA-test (b, N = 187), TCGA-entire (c, N = 374) and GEO (d, N = 278) dataset.; Figure S4: Risk Scores Correlated with Clinicopathological Features in TCGA train set. (a–c) The proportion of patients with different clinicopathologic features in high- and low-risk groups. (d–h) Risk scores in groups stratified by patient age, tumor grade, tumor stage, immune score and consensus cluster. (i–n) The Kaplan–Meier curve for patients with different risk score in subgroups stratified by clinical features; Figure S5: Risk score correlated with other clinical features. (a) The proportion of patients’ state in high- and low-risk groups. (b) The risk score of patients with different state. (c) Kaplan-Meier curves of OS for patients with chemotherapy. (d) Kaplan-Meier curves of OS for patients with BRCA1; Figure S6: A nomogram for OV patients in GEO dataset. (a) A nomogram for predicting the 1-, 2-, and 3-year OS of OV patients. (b–d) Calibration curves for the prediction of 1-, 2- or 3-year overall survival of OV patients; Figure S7: Relationships between the risk score and immune cell infiltration in GEO dataset. (a) The proportion of 22 immune cells infiltration in high- and low-risk groups. (b–i) Correlation of expression of the m6A-related lncRNAs and infiltration of specific immune cell type; Figure S8: Association of expression of the m6A-related lncRNAs with immune features. (a) The heatmap of correlation between the m6A-related lncRNAs and immune infiltrated cells. (b) The correlation between the m6A-related lncRNAs and immune score. (c–f) GSEA analysis indicating the most enriched pathways of the m6A-related lncRNAs in OV. Spearman correlation analysis was used to calculate the association; Figure S9: Expression levels of 4 m6A-related lncRNAs in cancers and adjacent normal tissues. (a) Heatmap showing the expression level of the m6A-related lncRNAs in tumor tissue compared with normal tissue in 18 cancer types which were composed of more than 5 normal samples. (b) Correlation plot based on Spearman Correlation analysis showing the correlation of gene expression among the 4 m6A-related lncRNAs across all 33 cancer types. (c) The association of 4 m6A-related lncRNAs with patients’ prognoses in pan-cancer. The representative top 2 cancer types among the m6A-related lncRNAs are shown according to the p-value; Figure S10: Association of expression of the m6A-related lncRNAs with patients’ prognoses and drug sensitivity in pan-cancer. (a) The forest plots with 95% confidence intervals and hazard ratios for overall survival for different cancer types in pan-cancer. (b) The scatter plot of the correlation between the m6A-related lncRNAs expression and drug sensitivity (the Z score of the CellMiner interface) for the Pearson correlation test using NCI-60 cell line data. top 16 associations are shown, ordered by p-value; Table S1: Primers’ information; Table S2: A total of 5 m6A-related prognostic lncRNAs in the TCGA-train dataset; Table S3: The conjoint univariate and multivariate Cox regression analyses of the independent predictor of OV.

Author Contributions

All authors contributed to the study conception and design. Conceptualization, J.B. and J.L.; Writing—Original Draft Preparation, R.G. and T.C.; Writing—Review and Editing, Z.Z. and S.N. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the Jiangsu Province Nature Science Foundation (Grant No. BK20220729).

Institutional Review Board Statement

The research was approval by the First Affiliated Hospital of Nanjing Medical University Ethics Committee. Protocol code:2022-SRFA-029. Approval date: 28 February 2022.

Informed Consent Statement

Not applicable.

Data Availability Statement

The datasets supporting the conclusions of this article are included within the article.

Acknowledgments

We would like to thank the researchers and study participants for their contributions.

Conflicts of Interest

The authors have no relevant financial or non-financial interest to disclose.

Abbreviations

ANOVAanalysis of variance
AUCarea under the curve
CIconfidence interval
DNAssstemness scores based on DNA methylation
GEOGene Expression Omnibus
GSEAgene set enrichment analysis
HRhazard ratio
ICIsimmune checkpoint inhibitors
IPSsimmunophenoscores
LASSOleast absolute shrinkage and selection operator
lncRNAslong noncoding RNAs
m6AN6-methyladenosine
mRLsN6-methyladenosine-related long noncoding RNAs
OSoverall survival
OVovarian serous cystadenocarcinoma
PCAprincipal component analysis
qRT-PCRquantitative real-time PCR
RNAssstemness scores based on mRNA
ROCreceiver-operating characteristic curves
ssGSEASingle Sample Gene Set Enrichment Analysis
TCGAThe Cancer Genome Atlas
TCIAThe Cancer Immunome Atlas
TIICstumor-infiltration immune cells
TMEtumor microenvironment
Tregregulatory T cells

References

  1. Siegel, R.L.; Miller, K.D.; Jemal, A. Cancer statistics. CA Cancer J. Clin. 2019, 69, 7–34. [Google Scholar] [CrossRef] [PubMed]
  2. Torre, L.A.; Trabert, B.; DeSantis, C.E.; Miller, K.D.; Samimi, G.; Runowicz, C.D.; Gaudet, M.M.; Jemal, A.; Siegel, R.L. Ovarian cancer statistics. CA Cancer J. Clin. 2018, 68, 284–296. [Google Scholar] [CrossRef] [PubMed]
  3. Andersen, C.L.; Sikora, M.J.; Boisen, M.M.; Ma, T.; Christie, A.; Tseng, G.; Park, Y.; Luthra, S.; Chandran, U.; Haluska, P.; et al. Active Estrogen Receptor-alpha Signaling in Ovarian Cancer Models and Clinical Specimens. Clin. Cancer Res. 2017, 23, 3802–3812. [Google Scholar] [CrossRef] [PubMed]
  4. Meng, C.; Zhou, J.-Q.; Liao, Y.-S. Autophagy-related long non-coding RNA signature for ovarian cancer. J. Int. Med. Res. 2020, 48, 300060520970761. [Google Scholar] [CrossRef]
  5. Henderson, J.T.; Webber, E.M.; Sawaya, G.F. Screening for Ovarian Cancer: Updated Evidence Report and Systematic Review for the US Preventive Services Task Force. JAMA 2018, 319, 595–606. [Google Scholar] [CrossRef]
  6. Huang, H.; Weng, H.; Chen, J. m6A Modification in Coding and Non-coding RNAs: Roles and Therapeutic Implications in Cancer. Cancer Cell 2020, 37, 270–288. [Google Scholar] [CrossRef]
  7. Shulman, Z.; Stern-Ginossar, N. The RNA modification N6-methyladenosine as a novel regulator of the immune system. Nat. Immunol. 2020, 21, 501–512. [Google Scholar] [CrossRef]
  8. Roundtree, I.A.; Evans, M.E.; Pan, T.; He, C. Dynamic RNA Modifications in Gene Expression Regulation. Cell 2017, 169, 1187–1200. [Google Scholar] [CrossRef]
  9. Wang, T.; Kong, S.; Tao, M.; Ju, S. The potential role of RNA N6-methyladenosine in Cancer progression. Mol. Cancer 2020, 19, 88. [Google Scholar] [CrossRef]
  10. Bhan, A.; Soleimani, M.; Mandal, S.S. Long Noncoding RNA and Cancer: A New Paradigm. Cancer Res. 2017, 77, 3965–3981. [Google Scholar] [CrossRef]
  11. Botti, G.; Scognamiglio, G.; Aquino, G.; Liguori, G.; Cantile, M. LncRNA HOTAIR in Tumor Microenvironment: What Role? Int. J. Mol. Sci. 2019, 20, 2279. [Google Scholar] [CrossRef]
  12. Peng, W.X.; Koirala, P.; Mo, Y.Y. LncRNA-mediated regulation of cell signaling in cancer. Oncogene 2017, 36, 5661–5667. [Google Scholar] [CrossRef]
  13. Chen, Y.; Lin, Y.; Shu, Y.; He, J.; Gao, W. Interaction between N6-methyladenosine (m6A) modification and noncoding RNAs in cancer. Mol. Cancer 2020, 19, 94. [Google Scholar] [CrossRef]
  14. Coker, H.; Wei, G.; Brockdorff, N. m6A modification of non-coding RNA and the control of mammalian gene expression. Biochim. Biophys. Acta Gene Regul. Mech. 2019, 1862, 310–318. [Google Scholar] [CrossRef]
  15. Yi, Y.-C.; Chen, X.-Y.; Zhang, J.; Zhu, J.-S. Novel insights into the interplay between m6A modification and noncoding RNAs in cancer. Mol. Cancer 2020, 19, 121. [Google Scholar] [CrossRef]
  16. Liu, H.; Xu, Y.; Yao, B.; Sui, T.; Lai, L.; Li, Z. A novel N6-methyladenosine (m6A)-dependent fate decision for the lncRNA THOR. Cell Death Dis. 2020, 11, 613. [Google Scholar] [CrossRef]
  17. Wu, Y.; Yang, X.; Chen, Z.; Tian, L.; Jiang, G.; Chen, F.; Li, J.; An, P.; Lu, L.; Luo, N.; et al. m6A-induced lncRNA RP11 triggers the dissemination of colorectal cancer cells via upregulation of Zeb1. Mol. Cancer 2019, 18, 87. [Google Scholar] [CrossRef]
  18. Dai, F.; Wu, Y.; Lu, Y.; An, C.; Zheng, X.; Dai, L.; Guo, Y.; Zhang, L.; Li, H.; Xu, W.; et al. Crosstalk between RNA m6A Modification and Non-coding RNA Contributes to Cancer Growth and Progression. Mol. Ther.-Nucleic Acids 2020, 22, 62–71. [Google Scholar] [CrossRef]
  19. Song, P.; Yang, F.; Jin, H.; Wang, X. The regulation of protein translation and its implications for cancer. Signal Transduct. Target. Ther. 2021, 6, 68. [Google Scholar] [CrossRef]
  20. Wang, Y.; Lu, J.-H.; Wu, Q.-N.; Jin, Y.; Wang, D.-S.; Chen, Y.-X.; Liu, J.; Luo, X.-J.; Meng, Q.; Pu, H.-Y.; et al. LncRNA LINRIS stabilizes IGF2BP2 and promotes the aerobic glycolysis in colorectal cancer. Mol. Cancer 2019, 18, 174. [Google Scholar] [CrossRef]
  21. Wang, J.; Ding, W.; Xu, Y.; Tao, E.; Mo, M.; Xu, W.; Cai, X.; Chen, X.; Yuan, J.; Wu, X. Long non-coding RNA RHPN1-AS1 promotes tumorigenesis and metastasis of ovarian cancer by acting as a ceRNA against miR-596 and upregulating LETM1. Aging 2020, 12, 4558–4572. [Google Scholar] [CrossRef]
  22. Tang, Z.; Li, C.; Kang, B.; Gao, G.; Li, C.; Zhang, Z. GEPIA: A web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017, 45, W98–W102. [Google Scholar] [CrossRef]
  23. Goldman, M.J.; Zhang, J.; Fonseca, N.A.; Cortés-Ciriano, I.; Xiang, Q.; Craft, B.; Piñeiro-Yáñez, E.; O’Connor, B.D.; Bazant, W.; Barrera, E.; et al. A user guide for the online exploration and visualization of PCAWG data. Nat. Commun. 2020, 11, 3400. [Google Scholar] [CrossRef]
  24. Zaccara, S.; Ries, R.J.; Jaffrey, S.R. Reading, writing and erasing mRNA methylation. Nat. Rev. Mol. Cell Biol. 2019, 20, 608–624. [Google Scholar] [CrossRef]
  25. Li, X.; Li, Y.; Yu, X.; Jin, F. Identification and validation of stemness-related lncRNA prognostic signature for breast cancer. J. Transl. Med. 2020, 18, 331. [Google Scholar] [CrossRef]
  26. Liu, J.; Xu, W.; Li, S.; Sun, R.; Cheng, W. Multi-omics analysis of tumor mutational burden combined with prognostic assessment in epithelial ovarian cancer based on TCGA database. Int. J. Med. Sci. 2020, 17, 3200–3213. [Google Scholar] [CrossRef]
  27. Yu, G.; Wang, L.-G.; Han, Y.; He, Q.-Y. clusterProfiler: An R Package for Comparing Biological Themes among Gene Clusters. OMICS 2012, 16, 284–287. [Google Scholar] [CrossRef]
  28. Newman, A.M.; Liu, C.L.; Green, M.R.; Gentles, A.J.; Feng, W.; Xu, Y.; Hoang, C.D.; Diehn, M.; Alizadeh, A.A. Robust enumeration of cell subsets from tissue expression profiles. Nat. Methods 2015, 12, 453–457. [Google Scholar] [CrossRef]
  29. Yoshihara, K.; Shahmoradgoli, M.; Martínez, E.; Vegesna, R.; Kim, H.; Torres-Garcia, W.; Trevino, V.; Shen, H.; Laird, P.W.; Levine, D.A.; et al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nat. Commun. 2013, 4, 2612. [Google Scholar] [CrossRef]
  30. Charoentong, P.; Finotello, F.; Angelova, M.; Mayer, C.; Efremova, M.; Rieder, D.; Hackl, H.; Trajanoski, Z. Pan-cancer Immunogenomic Analyses Reveal Genotype-Immunophenotype Relationships and Predictors of Response to Checkpoint Blockade. Cell Rep. 2017, 18, 248–262. [Google Scholar] [CrossRef]
  31. Reinhold, W.C.; Sunshine, M.; Varma, S.; Doroshow, J.H.; Pommier, Y. Using CellMiner 1.6 for Systems Pharmacology and Genomic Analysis of the NCI-60. Clin. Cancer Res. 2015, 21, 3841–3852. [Google Scholar] [CrossRef] [PubMed]
  32. Thorsson, V.; Gibbs, D.L.; Brown, S.D.; Wolf, D.; Bortone, D.S.; Ou Yang, T.H.; Porta-Pardo, E.; Gao, G.F.; Plaisier, C.L.; Eddy, J.A.; et al. The Immune Landscape of Cancer. Immunity 2018, 48, 812–830.e14. [Google Scholar] [CrossRef] [PubMed]
  33. Webb, P.M.; Jordan, S.J. Epidemiology of epithelial ovarian cancer. Best Pract. Res. Clin. Obstet. Gynaecol. 2017, 41, 3–14. [Google Scholar] [CrossRef] [PubMed]
  34. Shen, S.; Wang, G.; Zhang, R.; Zhao, Y.; Yu, H.; Wei, Y.; Chen, F. Development and validation of an immune gene-set based Prognostic signature in ovarian cancer. eBioMedicine 2019, 40, 318–326. [Google Scholar] [CrossRef]
  35. Millstein, J.; Budden, T.; Goode, E.; Anglesio, M.; Talhouk, A.; Intermaggio, M.; Leong, H.; Chen, S.; Elatre, W.; Gilks, B.; et al. Prognostic gene expression signature for high-grade serous ovarian cancer. Ann. Oncol. 2020, 31, 1240–1250. [Google Scholar] [CrossRef]
  36. Liu, J.; Meng, H.; Nie, S.; Sun, Y.; Jiang, P.; Li, S.; Yang, J.; Sun, R.; Cheng, W. Identification of a prognostic signature of epithelial ovarian cancer based on tumor immune microenvironment exploration. Genomics 2020, 112, 4827–4841. [Google Scholar] [CrossRef]
  37. Zheng, M.; Hu, Y.; Gou, R.; Nie, X.; Li, X.; Liu, J.; Lin, B. Identification three LncRNA prognostic signature of ovarian cancer based on genome-wide copy number variation. Biomed. Pharmacother. 2020, 124, 109810. [Google Scholar] [CrossRef]
  38. Fan, L.; Lin, Y.; Lei, H.; Shu, G.; He, L.; Yan, Z.; Rihan, H.; Yin, G. A newly defined risk signature, consisting of three m(6)A RNA methylation regulators, predicts the prognosis of ovarian cancer. Aging 2020, 12, 18453–18475. [Google Scholar] [CrossRef]
  39. Han, D.; Liu, J.; Chen, C.; Dong, L.; Liu, Y.; Chang, R.; Huang, X.; Liu, Y.; Wang, J.; Dougherty, U.; et al. Anti-tumour immunity controlled through mRNA m6A methylation and YTHDF1 in dendritic cells. Nature 2019, 566, 270–274. [Google Scholar] [CrossRef]
  40. Wang, H.; Hu, X.; Huang, M.; Liu, J.; Gu, Y.; Ma, L.; Zhou, Q.; Cao, X. Mettl3-mediated mRNA m6A methylation promotes dendritic cell activation. Nat. Commun. 2019, 10, 1898. [Google Scholar] [CrossRef]
  41. Wang, Q.; Chen, C.; Ding, Q.; Zhao, Y.; Wang, Z.; Chen, J.; Jiang, Z.; Zhang, Y.; Xu, G.; Zhang, J.; et al. METTL3-mediated m6A modification of HDGF mRNA promotes gastric cancer progression and has prognostic significance. Gut 2020, 69, 1193–1205. [Google Scholar] [CrossRef]
  42. Zhang, B.; Wu, Q.; Li, B.; Wang, D.; Wang, L.; Zhou, Y.L. m6A regulator-mediated methylation modification patterns and tumor microenvironment infiltration characterization in gastric cancer. Mol. Cancer 2020, 19, 53. [Google Scholar] [CrossRef]
  43. Lan, Q.; Liu, P.Y.; Haase, J.; Bell, J.L.; Hüttelmaier, S.; Liu, T. The Critical Role of RNA m6A Methylation in Cancer. Cancer Res. 2019, 79, 1285–1292. [Google Scholar] [CrossRef]
  44. He, L.; He, G. DNM3OS Facilitates Ovarian Cancer Progression by Regulating miR-193a-3p/MAP3K3 Axis. Yonsei Med. J. 2021, 62, 535–544. [Google Scholar] [CrossRef]
  45. Xia, X.; Zhang, H.; Xia, P.; Zhu, Y.; Liu, J.; Xu, K.; Yuan, Y. Identification of Glycolysis-Related lncRNAs and the Novel lncRNA WAC-AS1 Promotes Glycolysis and Tumor Progression in Hepatocellular Carcinoma. Front. Oncol. 2021, 11, 733595. [Google Scholar] [CrossRef]
  46. Zheng, J.; Zhou, Z.; Qiu, Y.; Wang, M.; Yu, H.; Wu, Z.; Wang, X.; Jiang, X. A Prognostic Ferroptosis-Related lncRNAs Signature Associated With Immune Landscape and Radiotherapy Response in Glioma. Front. Cell Dev. Biol. 2021, 9, 675555. [Google Scholar] [CrossRef]
  47. Shi, Z.; Shen, C.; Yu, C.; Yang, X.; Shao, J.; Guo, J.; Zhu, X.; Zhou, G. Long non-coding RNA LINC00997 silencing inhibits the progression and metastasis of colorectal cancer by sponging miR-512-3p. Bioengineered 2021, 12, 627–639. [Google Scholar] [CrossRef]
  48. Odunsi, K. Immunotherapy in ovarian cancer. Ann. Oncol. 2017, 28 (Suppl. S8), viii1–viii7. [Google Scholar] [CrossRef]
  49. Hartnett, E.G.; Knight, J.; Radolec, M.; Buckanovich, R.J.; Edwards, R.P.; Vlad, A.M. Immunotherapy Advances for Epithelial Ovarian Cancer. Cancers 2020, 12, 3733. [Google Scholar] [CrossRef]
  50. Osorio Ovarian Tumor Tissue Analysis (OTTA) Consortium; Goode, E.L.; Block, M.S.; Kalli, K.R.; Vierkant, R.A.; Chen, W.; Fogarty, Z.C.; Gentry-Maharaj, A.; Toloczko, A.; Hein, A.; et al. Dose-Response Association of CD8+ Tumor-Infiltrating Lymphocytes and Survival Time in High-Grade Serous Ovarian Cancer. JAMA Oncol. 2017, 3, e173290. [Google Scholar] [CrossRef]
  51. Tamborero, D.; Rubio-Perez, C.; Muiños, F.; Sabarinathan, R.; Piulats, J.M.; Muntasell, A.; Dienstmann, R.; Lopez-Bigas, N.; Gonzalez-Perez, A. A Pan-cancer Landscape of Interactions between Solid Tumors and Infiltrating Immune Cell Populations. Clin. Cancer Res. 2018, 24, 3717–3728. [Google Scholar] [CrossRef]
  52. Kreuzinger, C.; Geroldinger, A.; Smeets, D.; Braicu, E.I.; Sehouli, J.; Koller, J.; Wolf, A.; Darb-Esfahani, S.; Joehrens, K.; Vergote, I.; et al. A Complex Network of Tumor Microenvironment in Human High-Grade Serous Ovarian Cancer. Clin. Cancer Res. 2017, 23, 7621–7632. [Google Scholar] [CrossRef]
  53. Zheng, J.; Guo, J.; Cao, B.; Zhou, Y.; Tong, J. Identification and validation of lncRNAs involved in m6A regulation for patients with ovarian cancer. Cancer Cell Int. 2021, 21, 363. [Google Scholar] [CrossRef]
  54. Wei, Q.; Yang, D.; Liu, X.; Zhao, H.; Yang, Y.; Xu, J.; Liu, T.; Yi, P. Exploration of the Role of m6 A RNA Methylation Regulators in Malignant Progression and Clinical Prognosis of Ovarian Cancer. Front. Genet. 2021, 12, 650554. [Google Scholar] [CrossRef]
  55. Olbrecht, S.; Busschaert, P.; Qian, J.; Vanderstichele, A.; Loverix, L.; Van Gorp, T.; Van Nieuwenhuysen, E.; Han, S.; Van den Broeck, A.; Coosemans, A.; et al. High-grade serous tubo-ovarian cancer refined with single-cell RNA sequencing: Specific cell subtypes influence survival and determine molecular subtype classification. Genome Med. 2021, 13, 111. [Google Scholar] [CrossRef]
  56. Tan, Q.; Liu, H.; Xu, J.; Mo, Y.; Dai, F. Integrated analysis of tumor-associated macrophage infiltration and prognosis in ovarian cancer. Aging 2021, 13, 23210–23232. [Google Scholar] [CrossRef]
  57. Gao, Y.; Xu, Y.; Zhao, S.; Qian, L.; Song, T.; Zheng, J.; Zhang, J.; Chen, B. Growth differentiation factor-15 promotes immune escape of ovarian cancer via targeting CD44 in dendritic cells. Exp. Cell Res. 2021, 402, 112522. [Google Scholar] [CrossRef] [PubMed]
  58. Zhu, H.; Blum, R.H.; Bjordahl, R.; Gaidarova, S.; Rogers, P.; Lee, T.T.; Abujarour, R.; Bonello, G.B.; Wu, J.; Tsai, P.-F.; et al. Pluripotent stem cell–derived NK cells with high-affinity noncleavable CD16a mediate improved antitumor activity. Blood 2020, 135, 399–410. [Google Scholar] [CrossRef]
  59. Klingemann, H.; Boissel, L.; Toneguzzo, F. Natural Killer Cells for Immunotherapy—Advantages of the NK-92 Cell Line over Blood NK Cells. Front. Immunol. 2016, 7, 91. [Google Scholar] [CrossRef]
  60. Motz, G.T.; Coukos, G. Deciphering and Reversing Tumor Immune Suppression. Immunity 2013, 39, 61–73. [Google Scholar] [CrossRef]
  61. Banchereau, R.; Leng, N.; Zill, O.; Sokol, E.; Liu, G.; Pavlick, D.; Maund, S.; Liu, L.F.; Kadel, E., 3rd; Baldwin, N.; et al. Molecular determinants of response to PD-L1 blockade across tumor types. Nat. Commun. 2021, 12, 3969. [Google Scholar] [CrossRef] [PubMed]
  62. Maiorano, B.A.; Maiorano, M.F.P.; Lorusso, D.; Maiello, E. Ovarian Cancer in the Era of Immune Checkpoint Inhibitors: State of the Art and Future Perspectives. Cancers 2021, 13, 4438. [Google Scholar] [CrossRef] [PubMed]
  63. Yang, C.; Xia, B.-R.; Zhang, Z.-C.; Zhang, Y.-J.; Lou, G.; Jin, W.-L. Immunotherapy for Ovarian Cancer: Adjuvant, Combination, and Neoadjuvant. Front. Immunol. 2020, 11, 577869. [Google Scholar] [CrossRef]
  64. Miao, Y.R.; Thakkar, K.N.; Qian, J.; Kariolis, M.S.; Huang, W.; Nandagopal, S.; Yang, T.T.C.; Diep, A.N.; Cherf, G.M.; Xu, Y.; et al. Neutralization of PD-L2 is Essential for Overcoming Immune Checkpoint Blockade Resistance in Ovarian Cancer. Clin. Cancer Res. 2021, 27, 4435–4448. [Google Scholar] [CrossRef]
  65. Liu, Y.-T.; Sun, Z.-J. Turning cold tumors into hot tumors by improving T-cell infiltration. Theranostics 2021, 11, 5365–5386. [Google Scholar] [CrossRef]
  66. Chen, D.S.; Mellman, I. Elements of cancer immunity and the cancer-immune set point. Nature 2017, 541, 321–330. [Google Scholar] [CrossRef]
  67. Bruni, D.; Angell, H.K.; Galon, J. The immune contexture and Immunoscore in cancer prognosis and therapeutic efficacy. Nat. Rev. Cancer 2020, 20, 662–680. [Google Scholar] [CrossRef]
  68. Hegde, P.S.; Karanikas, V.; Evers, S. The Where, the When, and the How of Immune Monitoring for Cancer Immunotherapies in the Era of Checkpoint Inhibition. Clin. Cancer Res. 2016, 22, 1865–1874. [Google Scholar] [CrossRef]
  69. Gu, J.; Bi, F. Significance of N6-Methyladenosine RNA Methylation Regulators in Immune Infiltrates of Ovarian Cancer. Front. Genet. 2021, 12, 671179. [Google Scholar] [CrossRef]
  70. Jiang, Y.; Wan, Y.; Gong, M.; Zhou, S.; Qiu, J.; Cheng, W. RNA demethylase ALKBH5 promotes ovarian carcinogenesis in a simulated tumour microenvironment through stimulating NF-kappaB pathway. J. Cell Mol. Med. 2020, 24, 6137–6148. [Google Scholar] [CrossRef]
  71. Das, S.; Reddy, M.A.; Senapati, P.; Stapleton, K.; Lanting, L.; Wang, M.; Amaram, V.; Ganguly, R.; Zhang, L.; Devaraj, S.; et al. Diabetes Mellitus–Induced Long Noncoding RNA Dnm3os Regulates Macrophage Functions and Inflammation via Nuclear Mechanisms. Arter. Thromb. Vasc. Biol. 2018, 38, 1806–1820. [Google Scholar] [CrossRef] [PubMed]
  72. Lakhia, R.; Yheskel, M.; Flaten, A.; Ramalingam, H.; Aboudehen, K.; Ferrè, S.; Biggers, L.; Mishra, A.; Chaney, C.; Wallace, D.P.; et al. Interstitial microRNA miR-214 attenuates inflammation and polycystic kidney disease progression. JCI Insight 2020, 5, e133785. [Google Scholar] [CrossRef] [PubMed]
Cancers 14 04056 g001 550
Figure 1. (a) Flow chart of the study and (b) co-expression network diagram of m6A regulators and lncRNAs in the TCGA dataset.
Figure 1. (a) Flow chart of the study and (b) co-expression network diagram of m6A regulators and lncRNAs in the TCGA dataset.
Cancers 14 04056 g001
Cancers 14 04056 g002 550
Figure 2. Construction and validation of the mRL prognostic signature. Distribution of risk score, OS status (a–d) and heatmap of the expression of the four prognostic mRLs (e–h) in the TCGA-train, test, entire and GEO dataset. (i–l) Kaplan–Meier curves of OS for OV patients with a high- or low-risk score in the TCGA-train, test, entire and GEO dataset. (m–p) Time-dependent ROC analysis of risk score in predicting the prognosis in the TCGA-train, test, entire and GEO dataset. ROC, receiver-operating characteristic. In the TCGA-train set, N = 187; in the TCGA-test set, N = 187; in the TCGA-entire set, N = 374; in the GEO set, N = 278.
Figure 2. Construction and validation of the mRL prognostic signature. Distribution of risk score, OS status (a–d) and heatmap of the expression of the four prognostic mRLs (e–h) in the TCGA-train, test, entire and GEO dataset. (i–l) Kaplan–Meier curves of OS for OV patients with a high- or low-risk score in the TCGA-train, test, entire and GEO dataset. (m–p) Time-dependent ROC analysis of risk score in predicting the prognosis in the TCGA-train, test, entire and GEO dataset. ROC, receiver-operating characteristic. In the TCGA-train set, N = 187; in the TCGA-test set, N = 187; in the TCGA-entire set, N = 374; in the GEO set, N = 278.
Cancers 14 04056 g002
Cancers 14 04056 g003 550
Figure 3. A nomogram for OV patients in TCGA dataset. (a) A nomogram for predicting the 1-, 2- and 3-year OS of OV patients. (b–d). Calibration curves for the prediction of 1-, 2- or 3-year overall survival of OV patients. (e) Time-dependent ROC curves for each parameter and the combination of independent prognostic risk factors in the TCGA dataset (for predicting 1, 2 and 3-year OS). In the TCGA-entire set, N = 374.
Figure 3. A nomogram for OV patients in TCGA dataset. (a) A nomogram for predicting the 1-, 2- and 3-year OS of OV patients. (b–d). Calibration curves for the prediction of 1-, 2- or 3-year overall survival of OV patients. (e) Time-dependent ROC curves for each parameter and the combination of independent prognostic risk factors in the TCGA dataset (for predicting 1, 2 and 3-year OS). In the TCGA-entire set, N = 374.
Cancers 14 04056 g003
Cancers 14 04056 g004 550
Figure 4. The correlation between risk score and immune features in the TCGA dataset. (a,b) GSEA showed that immune response and immune system process were enriched in the high-risk group. (c–e) Stroma, immune and ESTIMATE scores in the high- and low-risk groups in TCGA dataset. (f–h) The correlation between risk score and stroma, immune and ESTIMATE scores in TCGA dataset. (i–k) Stroma, immune and ESTIMATE scores in the high- and low-risk groups in GEO dataset. (l–n) The correlation between risk score and stroma, immune and ESTIMATE scores in GEO dataset. In TCGA-entire set, N = 374; in GEO set, N = 278.
Figure 4. The correlation between risk score and immune features in the TCGA dataset. (a,b) GSEA showed that immune response and immune system process were enriched in the high-risk group. (c–e) Stroma, immune and ESTIMATE scores in the high- and low-risk groups in TCGA dataset. (f–h) The correlation between risk score and stroma, immune and ESTIMATE scores in TCGA dataset. (i–k) Stroma, immune and ESTIMATE scores in the high- and low-risk groups in GEO dataset. (l–n) The correlation between risk score and stroma, immune and ESTIMATE scores in GEO dataset. In TCGA-entire set, N = 374; in GEO set, N = 278.
Cancers 14 04056 g004
Cancers 14 04056 g005 550
Figure 5. Relationships between the risk score and immune cell infiltration and the role of the mRLs in tumor microenvironment in TCGA dataset. (a) The proportion of 22 immune cells infiltration in high- and low-risk groups. (b–i) Correlation of expression of the mRLs and infiltration of specific immune cell type. In TCGA-entire set, N = 374.
Figure 5. Relationships between the risk score and immune cell infiltration and the role of the mRLs in tumor microenvironment in TCGA dataset. (a) The proportion of 22 immune cells infiltration in high- and low-risk groups. (b–i) Correlation of expression of the mRLs and infiltration of specific immune cell type. In TCGA-entire set, N = 374.
Cancers 14 04056 g005
Cancers 14 04056 g006 550
Figure 6. IPSs and immunotherapy gene expression analysis. (a) The expression of PD-L2 in tumor and normal tissues. (b) The correlation between the mRLs and PD-L2. (c–f) The expression of PD-L2 in low- and high-risk groups in TCGA and GEO datasets. (g–j) The IPSs of low- and high-risk groups. * p < 0.05, ** p < 0.01 and *** p < 0.001. In TCGA-entire set, N = 374; In GEO set, N = 278.
Figure 6. IPSs and immunotherapy gene expression analysis. (a) The expression of PD-L2 in tumor and normal tissues. (b) The correlation between the mRLs and PD-L2. (c–f) The expression of PD-L2 in low- and high-risk groups in TCGA and GEO datasets. (g–j) The IPSs of low- and high-risk groups. * p < 0.05, ** p < 0.01 and *** p < 0.001. In TCGA-entire set, N = 374; In GEO set, N = 278.
Cancers 14 04056 g006
Cancers 14 04056 g007 550
Figure 7. Prognosis and immune infiltrations in consensus clustering subgroups of OV. (a) Consensus clustering matrix for k = 2. (b) Consensus clustering cumulative distribution function (CDF) for k = 2 to 9. (c) Relative change in area under the CDF curve for k = 2 to 9. (d) Tracking plot for k = 2 to 9. (e) Expression pattern of the mRLs in Cluster 1 and Cluster 2 subgroups. (f) The abundance of 21 immune cells in Cluster 1 and Cluster 2 subgroups. (g–j) Stroma, immune and ESTIMATE scores and tumor purity in Cluster 1 and Cluster 2 subgroups. (k) Kaplan–Meier analysis of patients in Cluster 1 and Cluster 2 subgroups. (l) The expression of PD-L2 in Cluster 1 and Cluster 2 subgroups.
Figure 7. Prognosis and immune infiltrations in consensus clustering subgroups of OV. (a) Consensus clustering matrix for k = 2. (b) Consensus clustering cumulative distribution function (CDF) for k = 2 to 9. (c) Relative change in area under the CDF curve for k = 2 to 9. (d) Tracking plot for k = 2 to 9. (e) Expression pattern of the mRLs in Cluster 1 and Cluster 2 subgroups. (f) The abundance of 21 immune cells in Cluster 1 and Cluster 2 subgroups. (g–j) Stroma, immune and ESTIMATE scores and tumor purity in Cluster 1 and Cluster 2 subgroups. (k) Kaplan–Meier analysis of patients in Cluster 1 and Cluster 2 subgroups. (l) The expression of PD-L2 in Cluster 1 and Cluster 2 subgroups.
Cancers 14 04056 g007
Cancers 14 04056 g008 550
Figure 8. Expression of pan-cancer mRLs. (a–d) Boxplot showing the expression level of the mRLs in tumor tissue compared with normal tissue in 18 cancer types which were composed of at least 5 normal samples. (e) Boxplot showing the expression distribution of the mRLs across pan-cancer. * p < 0.05, ** p < 0.01 and *** p < 0.001.
Figure 8. Expression of pan-cancer mRLs. (a–d) Boxplot showing the expression level of the mRLs in tumor tissue compared with normal tissue in 18 cancer types which were composed of at least 5 normal samples. (e) Boxplot showing the expression distribution of the mRLs across pan-cancer. * p < 0.05, ** p < 0.01 and *** p < 0.001.
Cancers 14 04056 g008
Cancers 14 04056 g009 550
Figure 9. Association of expression of the mRLs with immune subtypes, tumor microenvironment and cancer stemness in pan-cancer. (a) Association of expression of the m6A-related lncRNA with immune infiltrate subtypes across all the cancer types tested with ANOVA. (b,c) Correlation matrix between tumor microenvironment stromal scores (b) and immune scores (c) and the mRL expression using the ESTIMATE algorithm. (d,e) Correlation matrix plots showing the association between the m6A-related gene expression and cancer stemness DNA score (d) and RNA score (e). * p < 0.05, ** p < 0.01 and *** p < 0.001.
Figure 9. Association of expression of the mRLs with immune subtypes, tumor microenvironment and cancer stemness in pan-cancer. (a) Association of expression of the m6A-related lncRNA with immune infiltrate subtypes across all the cancer types tested with ANOVA. (b,c) Correlation matrix between tumor microenvironment stromal scores (b) and immune scores (c) and the mRL expression using the ESTIMATE algorithm. (d,e) Correlation matrix plots showing the association between the m6A-related gene expression and cancer stemness DNA score (d) and RNA score (e). * p < 0.05, ** p < 0.01 and *** p < 0.001.
Cancers 14 04056 g009
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Geng, R.; Chen, T.; Zhong, Z.; Ni, S.; Bai, J.; Liu, J. The m6A-Related Long Noncoding RNA Signature Predicts Prognosis and Indicates Tumor Immune Infiltration in Ovarian Cancer. Cancers 2022, 14, 4056. https://doi.org/10.3390/cancers14164056

AMA Style

Geng R, Chen T, Zhong Z, Ni S, Bai J, Liu J. The m6A-Related Long Noncoding RNA Signature Predicts Prognosis and Indicates Tumor Immune Infiltration in Ovarian Cancer. Cancers. 2022; 14(16):4056. https://doi.org/10.3390/cancers14164056

Chicago/Turabian Style

Geng, Rui, Tian Chen, Zihang Zhong, Senmiao Ni, Jianling Bai, and Jinhui Liu. 2022. "The m6A-Related Long Noncoding RNA Signature Predicts Prognosis and Indicates Tumor Immune Infiltration in Ovarian Cancer" Cancers 14, no. 16: 4056. https://doi.org/10.3390/cancers14164056

APA Style

Geng, R., Chen, T., Zhong, Z., Ni, S., Bai, J., & Liu, J. (2022). The m6A-Related Long Noncoding RNA Signature Predicts Prognosis and Indicates Tumor Immune Infiltration in Ovarian Cancer. Cancers, 14(16), 4056. https://doi.org/10.3390/cancers14164056

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop