ncRNA-mediated low expression of P2RY14 correlates with poor prognosis and tumor immune infiltration in ovarian carcinoma
Original Article

ncRNA-mediated low expression of P2RY14 correlates with poor prognosis and tumor immune infiltration in ovarian carcinoma

Chenggan Shu1#, Lifen Liu2#, Xiaoping Chen1, Jinling Xue1, Jiahong Fei1, Jianqing Wang1, Xiaoyue Yang3, Qi Peng4, Huaqin Yuan5

1Department of Obstetrics and Gynecology, Yancheng First Hospital, Affiliated Hospital of Nanjing University Medical School & the First People’s Hospital of Yancheng, Yancheng, China; 2Department of Obstetrics and Gynecology, The Second Affiliated Hospital of Soochow University, Suzhou, China; 3Department of Obstetrics and Gynecology, Jiangsu University Affiliated People’s Hospital, Zhenjiang, China; 4Department of Gynecology and Obstetrics, Yangzhou Jiangdu People’s Hospital, Yangzhou, China; 5Department of Oncology, Nanjing Gaochun People’s Hospital Affiliated to Yangzhou University, Nanjing, China

Contributions: (I) Conception and design: C Shu, H Yuan; (II) Administrative support: L Liu; (III) Provision of study materials or patients: J Xue, J Fei; (IV) Collection and assembly of data: X Chen, J Wang; (V) Data analysis and interpretation: X Yang, Q Peng; (VI) Manuscript writing: All authors; (VII) Final approval of manuscript: All authors.

#These authors contributed equally to this work.

Correspondence to: Huaqin Yuan. Department of Oncology, Nanjing Gaochun People’s Hospital Affiliated to Yangzhou University, Nanjing 211316, China. Email: huaqin_yuan@163.com.

Background: Ovarian cancer (OV) has been puzzling clinicians because of its poor prognosis. More and more evidence show that the G protein coupled receptor P2RY14 plays a key role in the initiation and progression of various types of human cancer. The purpose of our study is to explore the correlation between P2RY14 and the prognosis of ovarian cancer patients and the relevant mechanism.

Methods: First, the differentially expressed gene P2RY14 was screened from The Cancer Genome Atlas (TCGA) database. Explored possible P2RY14 related miRNAs and lncRNAs through multiple public databases, predicted and analyzed the expression level of candidate miRNAs and candidate lncRNAs that can bind to candidate miRNAs in OV through StarBase database. The TIMER database was used to comprehensively analyze the expression of tumor infiltrating immune cells, and to analyze the correlation between the expression level of P2RY14 and the level of immune cell infiltration in OV or the expression level of immune checkpoints.

Results: Patients with P2RY14 overexpression had better overall survival (OS) and progression-free interval (PFI). In the Targetscan database, 22 upstream miRNAs that may bind to P2RY14 were predicted. According to the regulatory network constructed by the Cytoscape software, correlation analysis and the role of miRNAs in the prognosis of OV, we first determined that the candidate miRNAs were miR-34c-5p. Then, we predicted the upstream lncRNAs of miR-34c-5p in the StarBase database, the expression level of these lncRNAs in OV in the Gene Expression Profiling Interactive Analysis (GEPIA) database, and the role in prognosis. We determined that LINC00665 is the most potential lncRNA upstream of ovarian cancer miRNA (hsa-miR-34c-5p)-P2RY14. Then, we analyzed the results in the Timer database, suggesting that P2RY14 expression was positively correlated with CD8+T Cell, CD4+T Cell, Macrophage, Neutral and Dendritic cells, and negatively correlated with B cells. Meanwhile, P2RY14 was positively correlated with CD274 and PDCD1.

Conclusions: P2RY14 can be used as a new predictive biomarker of ovarian cancer. Intervention of P2RY14 can affect the prognosis of ovarian cancer by affecting LINC00665-miR-34c-5p-P2RY14 axis. These findings provide a potential target for the development of anti-cancer strategies for ovarian cancer.

Keywords: P2RY14; LINC00665; miR-34c-5p; ovarian cancer (OV); prognosis


Submitted Oct 28, 2022. Accepted for publication Dec 19, 2022. Published online Jan 09 2023.

doi: 10.21037/atm-22-6120


Highlight box

Key findings

• P2RY14 is down regulated in many types of human cancers (including OV), and OV patients with high P2RY14 expression have better prognosis. This paper studies and determines the upstream regulation mechanism of P2RY14 in OV, namely LINC00665-miR-34c-5p-P2RY14 axis.

What is known and what is new?

• Ovarian cancer is a highly invasive tumor. It is of great significance to reveal the potential molecular mechanism of OV carcinogenesis and progression.

• In OV, we also explored the ncRNA related regulation of P2RY14, including miRNA and lncRNA. This study suggests that P2RY14 expression is related to immune cell infiltration, immune cell biomarkers or OV immune checkpoint.

What is the implication, and what should change now?

• It needs to be further verified through basic experiments and clinical trials in the future.


Introduction

Ovarian cancer (OV) is a highly aggressive tumor and the most prevalent female reproductive system malignancy (1). Annually, there are about 220,000 new cases of OV and over 140,000 fatalities worldwide (2). For tumors that are confined to the ovary, the effective treatment rate can reach 90%. However, the 5-year survival rate for OV is about 30%, and only 25% of patients may receive a diagnosis before their condition deteriorates (3). Even after surgery and chemotherapy, the overall survival (OS) rates of patients do not improve considerably; they are more likely to experience recurrence after surgery and acquire treatment resistance, which contributes to the high death rate of OV (4). The lack of identifiable tumor markers has become a serious issue in the early detection and monitoring of OV therapy.

The first purinoceptors, referred to known as P1 (adenosine) and P2 (ATP) receptors, were identified in 1978 (5). P2 (ATP) receptor is divided into P2X and P2Y, namely P2X (ion channel receptor) and P2Y (G protein coupled receptor). A subclass of G protein-coupled nucleotide receptors, known as P2Y receptors, primarily functions through adenine nucleoside triphosphate (ATP) (6). Thus far, nine functional P2Y receptors have been successfully cloned and confirmed using human tissue cells: P2Y1, 2, 4, 6, 11, 12, 13, and 14. The vacant numbers represent the P2Y receptor subtypes only found in non-mammals, whose function is unclear. The P2Y14 receptor is the only member of the P2 receptor family that reacts completely to uridine diphosphate glucose (UDP)-sugars, such as UDP-galactose and UDP-glucose (UDPG) (7). The up-regulation of the GPCR P2RY14 is associated with the activation of extracellular regulated protein kinases (ERK) signaling. It plays a role in patient survival, reducing interleukin-6 production in microgliomas and inhibiting glioma cell proliferation (8). The activity of P2RY14 via its ligand UDPG results in the selective phosphorylation of ERK1/2 (9). One study showed that P2RY14 could significantly prolong the OS of LUAD patients, which could be an important indicator of the prognosis of LUAD patients. New LUAD predictive biomarkers. P2RY14 may be related to LUAD immune invasion and play an important role in inhibiting tumor cell immune escape within LUAD microenvironment (10). However, at present, comprehensive research on the expression, prognosis, and mechanisms of P2RY14 in OV is insufficient. Moreover, the function of P2RY14 in tumor immune infiltration in OV is currently unknown.

With the rapid development of high-throughput sequencing technology, about 2% of ribonucleic acid (RNA) in mammalian genomes can encode proteins, and RNA that does not encode proteins is called non-coding RNA (ncRNA). ncRNA can be divided into two main categories according to length: <200 nucleotide ncRNAs, such as microRNA (miRNA) and small nuclear RNA (snRNA); >200 nucleotide ncRNAs, such as long non-coding RNA (lncRNA) and circular RNA (circRNA). Ncrnas are not only involved in the normal biological functions of the body, but also play an important role in the pathophysiological regulation of cells at various levels, such as epigenetics, transcription and post-transcription, which are further involved in the progression of a variety of diseases, especially in the occurrence and development of malignant tumors (11).

In the present study, we examined the relationship between P2RY14 expression and survival in various human cancer types. Subsequently, we investigated the regulation of P2RY14 by ncRNA, especially long non-coding RNA (lncRNA) and microRNA (miRNA), in OV. A correlation between P2RY14 expression, immune cell infiltration, biomarkers, and the OV immune checkpoint was established. This suggests that the ncRNA-mediated down-regulation of P2RY14 in OV patients is associated with a poor prognosis and immune involvement of the tumor. We present the following article in accordance with the REMARK reporting checklist (available at https://atm.amegroups.com/article/view/10.21037/atm-22-6120/rc).


Methods

Downloading, processing, and analysis of The Cancer Genome Atlas (TCGA) data

From TCGA database (https://genome-cancer.ucsc.edu/), we acquired the messenger RNA (mRNA) expression data of different types of cancer. We subsequently standardized this data and used the R program (https://www.r-project.org/) to perform a differential expression analysis of P2RY14.

Gene Expression Profiling Interactive Analysis (GEPIA) database analysis

We used the GEPIA database (http://gepia.cancer-pku.cn/) to determine the expressions of P2RY14 and lncRNAs in various cancer types and assess the prognostic significance of potential lncRNAs in OV. This database was also used to evaluate the relationship between P2RY14 expression and immunological checkpoints in OV.

Candidate miRNA prediction

To understand whether upstream genes regulate P2RY14, we further explored related miRNAs that may be bound to the 3'UTR (three prime untranslated region) end of P2RY14 mRNA using the TargetScan database (https://www.TargetScan.org).

StarBase database analysis

Using the StarBase database (https://starbase.sysu.edu.cn), we conducted expression correlation analysis on the miRNA-P2RY14 and lncRNA-miR-34c-5p in OV. The miR-34c-5p expression levels in OV patients and healthy controls were also investigated. Additionally, potential lncRNAs that might bind to miR-34c-5p were predicted.

Tumor Immune Estimation Resource (TIMER) database analysis

Using the TIMER database (https://cistrome.shinyapps.io/timer/), we thoroughly examined the expression of immune cells that had infiltrated tumors and studied the association between immune checkpoint expression, immune cell infiltration, and P2RY14 expression levels in OV. The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013).

Statistical analysis

As mentioned above, online resources were used to automatically calculate the statistical analysis for this investigation. A P value or log-rank P value of <0.05 was considered statistically significant.


Results

Expression of P2RY14 in different tumors

We evaluated the expression of P2RY14 in 33 human malignancies from TCGA database to investigate its potential role in carcinogenesis. P2RY14 was considerably up-regulated compared to standard samples in liver hepatocellular carcinoma (LIHC) and kidney renal clear cell carcinoma (KIRC). At the same time, it was significantly down-regulated in bladder urothelial carcinoma (BLCA), breast invasive carcinoma (BRCA), colon adenocarcinoma (COAD), kidney chromophobe (KICH), kidney renal papillary cell carcinoma (KIRP), lung adenocarcinoma (LUAD), lung squamous cell carcinoma (LUSC), prostate adenocarcinoma (PRAD), rectum adenocarcinoma (READ), stomach adenocarcinoma (STAD), and uterine corpus endometrial carcinoma (UCEC) (Figure 1A). Then, using the GEPIA database, the level of P2RY14 in various tumors was also investigated. We discovered that P2RY14 was up-regulated in LAML (acute myeloid leukemia) but down-regulated in LUSC, LUAD, BLCA, COAD, and READ (Figure 1B).

Figure 1 Expression of P2RY14 in several tumor types. P2RY14 expression in various malignancies found in TCGA database (A); P2RY14 expression in various tumors found in the GEPIA database (B). *, P<0.05; **, P<0.01; ***, P<0.001; ns, not significant. TPM, transcripts per million; COAD, colon adenocarcinoma; READ, rectum adenocarcinoma; LUSC, lung squamous cell carcinoma; BLCA, bladder urothelial carcinoma; LUAD, lung adenocarcinoma; LAML, acute myeloid leukemia; TCGA, The Cancer Genome Atlas; GEPIA, Gene Expression Profiling Interactive Analysis.

The prognosis of P2RY14 in different tumors

Using the TCGA database, we also explored the relationship between P2RY14 and several tumor prognoses, including OS and progression-free survival [or progression-free interval (PFI)]. Individuals with high P2RY14 expression in the LUAD, KIRC, OV, cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC), sarcoma (SARC), and oral squamous cell carcinoma (OSCC) had a superior OS. However, those with elevated expression in the BLCA and STAD had a worse OS, as shown in Figure 2A. Next, we observed that patients with low expression of P2RY14 in testicular germ cell tumors (TGCT), STAD, and ESCC had a better PFI, while those with high P2RY14 expression in CESC, KIRC, and OV had a better PFI (Figure 2B).

Figure 2 P2RY14 prognosis in several tumor types from TCGA database. (A) The association between P2RY14 expression and OS in various malignancies; (B) the association between PFI and P2RY14 expression in various cancers. P<0.05 was considered statistically significant. BLCA, bladder urothelial carcinoma; LUAD, lung adenocarcinoma; KIRC, kidney renal clear cell carcinoma; STAD, stomach adenocarcinoma; OV, ovarian cancer; CESC, cervical squamous cell carcinoma and endocervical adenocarcinoma; SARC, sarcoma; OSCC, oral squamous cell carcinoma; TGCT, testicular germ cell tumors; ESCC, European Space Components Coordination; TCGA, The Cancer Genome Atlas; OS, overall survival; PFI, progression-free interval.

P2RY14 upstream miRNA identification and validation

It is widely believed that ncRNAs control how genes are expressed. We initially anticipated the upstream miRNAs that would bind to P2RY14 in the TargetScan database, and ultimately discovered 22 candidate miRNAs to determine whether certain ncRNAs regulate P2RY14. The details of these miRNAs are shown in Table 1. Cytoscape software was used to create a miRNA-P2RY14 regulatory network to improve visualization (Figure 3). There must be a negative relationship between miRNA and P2RY14 based on how miRNA regulates the expression of its gene products. As a result, we analyzed expression correlation in the StarBase database, as presented in Table 2. At the same time, we constructed a scatter plot to express this correlation more intuitively (Figure 4). We initially screened and obtained 13 negatively correlated candidate miRNAs. The expression of putative miRNAs with prognostic significance in the prediction of OV was examined using the StarBase database (Figure 5). We identified several potential miRNAs, including hsa-miR-4770 and hsa-miR-34c-5p. Finally, we chose hsa-miR-34c-5p as our final candidate miRNA based on previously published research results.

Table 1

Predicted miRNAs in the TargetScan database that bind to the 3'UTR end of P2RY14

Candidate miRNA Position in
the UTR
Seed match Context++ score Context++ score percentile Weighted context++ score Conserved branch length PCT Predicted relative KD
hsa-miR-140-5p 285–292 8mer −0.46 99 −0.46 0.233 <0.1 −3.823
hsa-miR-449a 106–113 8mer −0.48 98 −0.48 0.131 <0.1 −4.531
hsa-miR-449b-5p 106–113 8mer −0.49 98 −0.49 0.131 <0.1 −4.221
hsa-miR-34a-5p 106–113 8mer −0.5 98 −0.5 0.131 <0.1 −4.639
hsa-miR-34c-5p 106–113 8mer −0.52 98 −0.52 0.131 <0.1 −4.221
hsa-miR-129-2-3p 453–460 8mer −0.39 98 −0.39 1.654 <0.1 −4.019
hsa-miR-129-1-3p 453–460 8mer −0.39 98 −0.39 1.654 <0.1 −4.019
hsa-miR-143-3p 96–102 7mer-m8 −0.32 98 −0.32 1.554 0.23 −4.053
hsa-miR-4770 96–102 7mer-m8 −0.32 98 −0.32 1.554 0.23 −4.076
hsa-miR-6088 96–102 7mer-m8 −0.27 97 −0.27 1.554 0.23 −4.104
hsa-miR-219a-5p 529–536 8mer −0.34 96 −0.34 0.048 <0.1 −4.895
hsa-miR-6766-3p 529–536 8mer −0.34 96 −0.34 0.048 <0.1 −5.152
hsa-miR-4782-3p 529–536 8mer −0.34 96 −0.34 0.048 <0.1 −5.152
hsa-miR-106b-5p 630–636 7mer-m8 −0.16 92 −0.16 0.403 <0.1 −4.427
hsa-miR-20b-5p 630–636 7mer-m8 −0.13 91 −0.13 0.403 <0.1 −5.412
hsa-miR-106a-5p 630–636 7mer-m8 −0.15 91 −0.15 0.403 <0.1 −5.029
hsa-miR-20a-5p 630–636 7mer-m8 −0.15 91 −0.15 0.403 <0.1 −4.491
hsa-miR-17-5p 630–636 7mer-m8 −0.14 91 −0.14 0.403 <0.1 −5.537
hsa-miR-9-5p 15–21 7mer-1A −0.16 91 −0.16 0.066 <0.1 −2.103
hsa-miR-302e 629–635 7mer-m8 −0.16 90 −0.16 0.403 <0.1 −4.39
hsa-miR-372-3p 629–635 7mer-m8 −0.18 90 −0.18 0.403 <0.1 −4.965
hsa-miR-93-5p 630–636 7mer-m8 −0.12 90 −0.12 0.403 <0.1 −5.486

UTR, untranslated region; PCT, percentage; KD, k-dimensional.

Figure 3 Cytoscape-generated miRNA-P2RY14 regulatory network.

Table 2

Correlation between the candidate miRNAs in the StarBase database and P2RY14 in ovarian cancer

Candidate miRNA Gene R value P value
hsa-miR-34a-5p P2RY14 0.001 9.86e-1
hsa-miR-449b-5p P2RY14 −0.225 1.10e-5
hsa-miR-449a P2RY14 −0.199 1.00e-4
hsa-miR-140-5p P2RY14 0.073 1.58e-1
hsa-miR-34c-5p P2RY14 −0.046 3.72e-1
hsa-miR-129-2-3p P2RY14 0.028 5.84e-1
hsa-miR-129-1-3p P2RY14 −0.066 2.02e-1
hsa-miR-143-3p P2RY14 0.131 1.09e-2
hsa-miR-4770 P2RY14 −0.096 6.22e-2
hsa-miR-6088 P2RY14 0.000 1.00e+0
hsa-miR-4782-3p P2RY14 0.073 1.60e-1
hsa-miR-6766-3p P2RY14 0.058 2.60e-1
hsa-miR-219a-5p P2RY14 −0.106 3.95e-2
hsa-miR-106b-5p P2RY14 −0.163 1.52e-3
hsa-miR-9-5p P2RY14 −0.034 5.15e-1
hsa-miR-17-5p P2RY14 −0.102 4.73e-2
hsa-miR-20a-5p P2RY14 −0.108 3.59e-2
hsa-miR-372-3p P2RY14 −0.094 6.94e-2
hsa-miR-302e P2RY14 0.000 1.00e+0
hsa-miR-20b-5p P2RY14 −0.071 1.71e-1
hsa-miR-106a-5p P2RY14 0.002 9.63e-1
hsa-miR-93-5p P2RY14 −0.022 6.75e-1
Figure 4 Correlation between miRNA and P2RY14 in the StarBase database. P<0.05 was considered statistically significant. OV, ovarian cancer; ENROCI, The Encyclopedia of RNA Interactomes (http://starbase.sysu.edu.cn/index.php); FPKM, fragments per kilobase of exon model per million mapped fragments; RPM, reads per million mapped reads.
Figure 5 Expression of candidate miRNAs in the StarBase database in the prognosis of OV. P<0.05 was considered statistically significant. OV, ovarian cancer.

LncRNA prediction and analysis upstream of miR-34c-5p

Next, Cytoscape software was used to generate a lncRNA-miR-34c-5p regulatory network using the TargetScan database to determine the upstream lncRNA of miR-34c-5p (Figure 6). Based on the assumption that the expression differential of lncRNAs in OV and surrounding tissues was statistically significant according to the binding of base pairs, we preliminarily retrieved six potential lncRNAs at the same time. The specifics of these lncRNAs are displayed in Table 3. Using the GEPIA database, we then calculated the overexpression of these lncRNAs in OV. LINC00665, LINC01123, and LINC01132 were considerably up-regulated in OV tumor tissues compared to the normal controls, but LINC00667, LINC00921, and NEAT1 were significantly down-regulated (Figure 7). The predictive significance of these six lncRNAs in OV was then assessed in the StarBase database. Only LINC00665 and LINC01132 were quite different, and the prognosis for OV patients with elevated levels of both genes was favorable (Figure 8). In the competitive endogenous RNA (ceRNA) concept, lncRNAs can efficiently enhance mRNA expression by adhering to conserved miRNAs. Therefore, there must be a negative association between miRNA and lncRNA, as demonstrated in Table 4 and Figure 9.

Figure 6 Cytoscape software was used to design the lncRNA-miR-34c-5p regulatory network.

Table 3

Upstream lncRNAs predicted to regulate hsa-miR-34c-5p in ovarian cancer from the TargetScan database

Gene ID Gene name clipExpNum RBP merClass miRseq Align pancancerNum
ENSG00000232677 LINC00665 5 AGO1-4 7mer-m8 cguUAGUC--GAUUGAU-GUGACGGa | ||| |: ||: ||||||| 6
ENSG00000204588 LINC01123 4 AGO1-4, AGO2 7mer-m8 cguUAGUCGAUUGAUGUGACGGa |||||| |:| ||||||| 5
ENSG00000227630 LINC01132 1 AGO1-4 8mer cgUUAGUCGAUUGAUGUGACGGa || ::| ::|| ||||||| 3
ENSG00000245532 NEAT1 4 AGO1-4, AGO2 7mer-m8 cguuaGUCGAUUGAUGUGACGGa |||| ||:||||||| 3
ENSG00000263753 LINC00667 3 AGO2 7mer-m8 cgUUAGUCGAU-UG---AUGUGACGGa |||| ||| :| |:||||||| 2
ENSG00000281005 LINC00921 1 AGO2 7mer-m8 cguUAGUC-GAUU---GAUGUGACGGa ||||| |||| | :||||||| 1
Figure 7 Expression of candidate lncRNAs in OV from the GEPIA database. *, P<0.05. TPM, transcripts per million; OV, ovarian cancer; GEPIA, Gene Expression Profiling Interactive Analysis.
Figure 8 Correlation between candidate lncRNAs in the StarBase database and the prognosis of OV patients. P<0.05 was considered statistically significant. TPM, transcripts per million; HR, hazard ratio; OV, ovarian cancer.

Table 4

Results of a correlation study of predicted lncRNA and miRNA (hsa-miR-34c-5p)

lncRNA miRNA R value P value
LINC00665 hsa-miR-34c-5p −0.013 7.96e-1
LINC01132 hsa-miR-34c-5p 0.093 7.25e-2

P<0.05 was considered statistically significant.

Figure 9 Correlation between LINC00665 and LINC01132 and miR-34c-5p in the StarBase database. P<0.05 was considered statistically significant. ACC, adrenocortical carcinoma; OV, ovarian cancer; ENROCI, The Encyclopedia of RNA Interactomes (http://starbase.sysu.edu.cn/index.php); FPKM, fragments per kilobase of exon model per million mapped fragments; RPM, reads per million mapped reads.

Furthermore, the correlation between these two lncRNAs and miR-34c-5p in OV was determined using the StarBase database. The results of the expression, survival, and correlation analyses showed that LINC00665 is a promising lncRNA upstream of the OV miRNA (hsa-miR-34c-5p)-P2RY14.

Correlation analysis between P2RY14 in the TIMER database and immune cell infiltration in OV

A previous study revealed that P2RY14 also has a role in the immune system. The activity and mechanisms of P2RY14 can be better understood through correlation analysis. Therefore, we examined P2RY14 in OV and immune cell invasion in the TIMER database. P2RY14 expression was positively correlated with dendritic cells, neutrophils, macrophages, CD8+ T-cells, and CD4+ T-cells but negatively associated with B-cells (Figure 10).

Figure 10 Analysis of the correlation between P2RY14 and immune cell invasion in OV from the TIMER database. **, P<0.01; ***, P<0.001 were considered statistically significant. OV, ovarian cancer; TPM, transcripts per million.

Relationship between P2RY14 and immune cell biomarker expression in OV

To broadly comprehend the significance of P2RY14 in tumor immunity, we compared P2RY14 expression with immune cell markers in OV using the GEPIA database. As demonstrated in Table 5, dendritic cells (NRP1, CD1C, HLA-DPA1, HLA-DQB1, HLA-DPB1, and HLA-DRA), M1 biomarkers (IRF5, PTGS2, and NOS2), M2 biomarkers (VSIG4), B-cell biomarkers (CD79A), and M2 biomarkers (ITGAM) were all positively associated.

Table 5

GEPIA database investigation of the relationship between immune cell markers and P2RY14 in ovarian cancer

Immune cell Biomarker R value P value
B cell CD19 0.0033 0.95
CD79A 0.21 7.8e-06
CD8+ T cell CD8A 0.45 0
CD8B 0.043 0.38
CD4+ T cell CD4 0.43 0
M1 IRF5 0.13 0.007
PTGS2 0.15 0.0023
NOS2 0.12 0.012
M2 CD163 0.44 0
VSIG4 0.36 7.5e-15
MS4A4A 0.61 0
Neutrophils CEACAM8 0.06 0.22
ITGAM 0.38 4.4e-16
CCR7 0.47 0
Dendritic cell CD1C 0.32 1.9e-11
HLA-DPA1 0.27 8.8e-09
HLA-DRA 0.26 6.3e-08
HLA-DQB1 0.18 0.00022
HLA-DPB1 0.29 1.4e-09
NRP1 0.29 1.5e-09
ITGAX 0.39 0

GEPIA, Gene Expression Profiling Interactive Analysis.

Connection between immunological checkpoints and P2RY14 in OV

Important immunological checkpoints are responsible for tumor immune escape, including cytotoxic T lymphocyte-associated antigen-4 (CTLA-4) and programmed death receptor 1/programmed death receptor ligand-1 (PD1/PD-L1). So, we assessed the connection between P2RY14 and PD-L1, CTLA-4, or PD1. We discovered that CD274 and programmed cell death 1 (PDCD1) had strong positive correlations with P2RY14. We also identified a substantial positive association between P2RY14, CD274, and PDCD1 in OV, which was similar to the TIMER database analysis results (Figure 11). These findings imply that P2RY14-mediated OV carcinogenesis may involve tumor immune escape.

Figure 11 Relationship between P2RY14 expression and CD274, CTLA4, and PDCD1 expression in OV. P<0.05 was considered statistically significant. TPM, transcripts per million; OV, ovarian cancer.

Discussion

At present, OV is still associated with a poor prognosis. Understanding the molecular process of P2RY14 could offer important insights for designing effective treatment targets or identifying potential prognostic indicators. Increasing evidence suggests that P2RY14 may be essential for the onset and development of numerous cancers and illnesses. However, research on P2RY14 in OV is still lacking; thus, further investigation and exploration are required. In this study, TCGA database was utilized to analyze P2RY14 expression in a wide variety of cancers. Then, the GEPIA database was employed to validate our findings. At the same time, we discovered that the P2RY14 expression was up-regulated in paracancerous tissues, and individuals with high P2RY14 expression in OV had an improved OS and PFI. These findings reflect the significance of our study on target molecules.

Circular RNAs (circRNAs), miRNAs, lncRNAs, and other ncRNAs interact via the ceRNA mechanism and control gene expression, as has been well-established (12-16). After estimating the upstream miRNAs that could attach to P2RY14 in the TargetScan database, we discovered 22 potential miRNAs. These candidates had various expressions and functions in OV, including overexpression of the lncRNA H19. We found that the miR-140-5p gene was down-regulated and the phosphatidylinositide 3-kinases/protein kinase B (PI3K/AKT) signaling pathway was activated, encouraging OV cells to invade, migrate, and undergo epithelial-mesenchymal transition (17). MicroRNA-449a down-regulates NOTCH1 in OV cells, which decreases cell survival and improves the effects of cisplatin-induced cytotoxicity by platinum (18). Through the control of mitogen-activated protein kinase 3 (MAP3K3) by miR-129-3p, circSETDB1 knockdown prevents the malignant evolution of serous OV (19). Also, the transforming growth factor-β activating kinase 1 (TAK1) axis, which is activated by the lncRNA MCM3AP antisense RNA 1 (non-protein coding) (MCM3AP-AS1), induces OV progression (20). After correlation and prognostic analysis, we found that hsa-miR-34c-5p and hsa-miR-4770 were our potential candidate miRNAs. Referring to relevant research reports, we found that there are few studies on miR-4770 at present, and tumor related studies are only colorectal cancer. However, there are abundant researches on miR-34c-5p, for example, the upregulation of miR-34c-5p inhibits nasopharyngeal carcinoma cells through mediating NOTCH1 (21). MiR-34c-5p mediates the malignant behavior of oral squamous cell carcinoma through targeting and TRIM29 (22). LncRNA-ES3/miR-34c-5p/BMF axis participates in regulating VSMC calcification/senescence induced by high glucose (23). At the same time, a study reported that the characteristics of ovarian cancer extracellular vesicle miRNA (OCEM) developed based on miRNAs including miR-34c-5p can improve the diagnosis and treatment of ovarian cancer (24). At the same time, the malignant ascites derived extracellular vesicles can significantly promote the invasion of ovarian cancer cells and the growth of ascites derived organs (25). After analyzing the correlation and predictive data, we ultimately decided that miR-34c-5p was the best potential miRNA. Next, we projected the upstream lncRNAs and identified LINC00665 as the most likely lncRNA upstream of OV miRNA (hsa-miR-34c-5p)-P2RY14 through correlation, survival, and expression analyses. Additionally, substantial evidence suggests that the lncRNA LINC00665 targets miR-181a-5p/FH2 domain containing (FHDC) and stimulates OV cell proliferation while inhibiting apoptosis (26). Furthermore, a possible OV regulatory mechanism was found to involve the LINC00665-miR-34c-5p-P2RY14 axis.

Numerous studies have shown that tumor immune cell invasion affects the effectiveness of chemotherapy, immunotherapy, or radiation as well as the survival of cancer patients (27-29). Our findings indicate a strong association between P2RY14 and various immune system cells in OV, including B cells, CD8+ T cells, CD4+ T cells, macrophages, neutrophils, and dendritic cells. Additionally, P2RY14 was strongly linked to these invading immune cell indicators. These results imply that the function of P2RY14 in OV could be influenced by tumor immune participation.

Also, proper immunological checkpoint expression is necessary for immunotherapy effectiveness as well as the infiltration of immune cells into the tumor microenvironment (30). Consequently, the relationship between P2RY14 and immunological checkpoints was also investigated. We discovered that P2RY14 was significantly linked to CD274 and PDCD1, which further substantiated the idea that focusing on P2RY14 could improve the effectiveness of immunotherapy in OV.


Conclusions

In summary, our research shows that P2RY14 is down-regulated in various human malignancies, including OV, and that OV patients with high P2RY14 expression levels tend to have better prognoses. We determined the LINC00665-miR-34c-5p-P2RY14 axis as the upstream regulation mechanism of P2RY14 in OV. However, according to our most recent research, P2RY14 may modify the tumor immunological microenvironment, which may impact the prognosis of OV. Therefore, in the future, basic research and clinical studies will be needed to confirm the current evidence.


Acknowledgments

Funding: This study was supported by the Open Project of Soochow University (No. GZK1202106).


Footnote

Reporting Checklist: The authors have completed the REMARK reporting checklist. Available at https://atm.amegroups.com/article/view/10.21037/atm-22-6120/rc

Conflicts of Interest: All authors have completed the ICMJE uniform disclosure form (available at https://atm.amegroups.com/article/view/10.21037/atm-22-6120/coif). The authors have no conflicts of interest to declare.

Ethical Statement: The authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013).

Open Access Statement: This is an Open Access article distributed in accordance with the Creative Commons Attribution-NonCommercial-NoDerivs 4.0 International License (CC BY-NC-ND 4.0), which permits the non-commercial replication and distribution of the article with the strict proviso that no changes or edits are made and the original work is properly cited (including links to both the formal publication through the relevant DOI and the license). See: https://creativecommons.org/licenses/by-nc-nd/4.0/.


References

  1. Siegel RL, Miller KD, Jemal A. Cancer Statistics, 2017. CA Cancer J Clin 2017;67:7-30. [Crossref] [PubMed]
  2. Torre LA, Trabert B, DeSantis CE, et al. Ovarian cancer statistics, 2018. CA Cancer J Clin 2018;68:284-96. [Crossref] [PubMed]
  3. Shinmura H, Yoneyama K, Harigane E, et al. Use of tumor markers to distinguish endometriosis-related ovarian neoplasms from ovarian endometrioma. Int J Gynecol Cancer 2020;30:831-6. [Crossref] [PubMed]
  4. Li SS, Ma J, Wong AST. Chemoresistance in ovarian cancer: exploiting cancer stem cell metabolism. J Gynecol Oncol 2018;29:e32. [Crossref] [PubMed]
  5. Kolachala VL, Bajaj R, Chalasani M, et al. Purinergic receptors in gastrointestinal inflammation. Am J Physiol Gastrointest Liver Physiol 2008;294:G401-10. [Crossref] [PubMed]
  6. Grobben B, Claes P, Roymans D, et al. Ecto-nucleotide pyrophosphatase modulates purinoceptor-mediated-mediated signal transduction and is inhibited by purinoceptor antagonists. Br J Pharmacol 2000;130:139-45. [Crossref] [PubMed]
  7. Chambers JK, Macdonald LE, Sarau HM, et al. A G protein-coupled receptor for UDP-glucose. J Biol Chem 2000;275:10767-71. [Crossref] [PubMed]
  8. Curet MA, Watters JJ. P2Y14 receptor activation decreases interleukin-6 production and glioma GL261 cell proliferation in microglial transwell cultures. J Neurooncol 2018;137:23-31. [Crossref] [PubMed]
  9. Shah K, Moharram SA, Kazi JU. Acute leukaemia cells resistant to PI3K/mTOR inhibition display upregulation of P2RY14 expression. Clin Epigenetics 2018;10:83. [Crossref] [PubMed]
  10. Xu T, Xu S, Yao Y, et al. P2RY14 downregulation in lung adenocarcinoma: a potential therapeutic target associated with immune infiltration. J Thorac Dis 2022;14:515-35. [Crossref] [PubMed]
  11. Beg A, Parveen R, Fouad H, et al. Role of different non-coding RNAs as ovarian cancer biomarkers. J Ovarian Res 2022;15:72. [Crossref] [PubMed]
  12. Lou W, Ding B, Wang J, et al. The Involvement of the hsa_circ_0088494-miR-876-3p-CTNNB1/CCND1 Axis in Carcinogenesis and Progression of Papillary Thyroid Carcinoma. Front Cell Dev Biol 2020;8:605940. [Crossref] [PubMed]
  13. Gao S, Ding B, Lou W. microRNA-Dependent Modulation of Genes Contributes to ESR1's Effect on ERalpha Positive Breast Cancer. Front Oncol 2020;10:753. [Crossref] [PubMed]
  14. Ghafouri-Fard S, Shoorei H, Anamag FT, et al. The Role of Non-Coding RNAs in Controlling Cell Cycle Related Proteins in Cancer Cells. Front Oncol 2020;10:608975. [Crossref] [PubMed]
  15. Razavi ZS, Tajiknia V, Majidi S, et al. Gynecologic cancers and noncoding RNAs: Epigenetic regulators with emerging roles. Crit Rev Oncol Hematol 2021;157:103192. [Crossref] [PubMed]
  16. Fabrizio FP, Sparaneo A, Muscarella LA. NRF2 Regulation by Noncoding RNAs in Cancers: The Present Knowledge and the Way Forward. Cancers (Basel) 2020;12:3621. [Crossref] [PubMed]
  17. Xu H, Ding Y, Yang X. Overexpression of Long Noncoding RNA H19 Downregulates miR-140-5p and Activates PI3K/AKT Signaling Pathway to Promote Invasion, Migration and Epithelial-Mesenchymal Transition of Ovarian Cancer Cells. Biomed Res Int 2021;2021:6619730. [Crossref] [PubMed]
  18. Zhou Y, Chen Q, Qin R, et al. MicroRNA-449a reduces cell survival and enhances cisplatin-induced cytotoxicity via downregulation of NOTCH1 in ovarian cancer cells. Tumour Biol 2014;35:12369-78. [Crossref] [PubMed]
  19. Li B, Zhang L. CircSETDB1 knockdown inhibits the malignant progression of serous ovarian cancer through miR-129-3p-dependent regulation of MAP3K3. J Ovarian Res 2021;14:160. [Crossref] [PubMed]
  20. Wen J, Han S, Cui M, et al. Long noncoding RNA MCM3APAS1 drives ovarian cancer progression via the microRNA1433p/TAK1 axis. Oncol Rep 2020;44:1375-84. [PubMed]
  21. Xiao S, Guo J, Zhang W, et al. A Six-microRNA Signature Nomogram for Preoperative Prediction of Tumor Deposits in Colorectal Cancer. Int J Gen Med 2022;15:675-87. [Crossref] [PubMed]
  22. Xu X, Yan H, Zhang L, et al. Up-regulation of miR-34c-5p inhibits nasopharyngeal carcinoma cells by mediating NOTCH1. Biosci Rep 2020;40:BSR20200302. [Crossref] [PubMed]
  23. Shen Y, Sun C, Zhao B, et al. miR-34c-5p mediates the cellular malignant behaviors of oral squamous cell carcinoma through targeted binding of TRIM29. Ann Transl Med 2021;9:1537. [Crossref] [PubMed]
  24. Lin X, Zhan JK, Zhong JY, et al. lncRNA-ES3/miR-34c-5p/BMF axis is involved in regulating high-glucose-induced calcification/senescence of VSMCs. Aging (Albany NY) 2019;11:523-35. [Crossref] [PubMed]
  25. Wang W, Jo H, Park S, et al. Integrated analysis of ascites and plasma extracellular vesicles identifies a miRNA-based diagnostic signature in ovarian cancer. Cancer Lett 2022;542:215735. [Crossref] [PubMed]
  26. Wang S, Wang Y, Lu J, et al. LncRNA LINC00665 Promotes Ovarian Cancer Cell Proliferation and Inhibits Apoptosis via Targeting miR-181a-5p/FHDC. Appl Biochem Biotechnol 2022;194:3819-32. [Crossref] [PubMed]
  27. Waniczek D, Lorenc Z, Śnietura M, et al. Tumor-Associated Macrophages and Regulatory T Cells Infiltration and the Clinical Outcome in Colorectal Cancer. Arch Immunol Ther Exp (Warsz) 2017;65:445-54. [Crossref] [PubMed]
  28. Zhang H, Liu H, Shen Z, et al. Tumor-infiltrating Neutrophils are Prognostic and Predictive for Postoperative Adjuvant Chemotherapy Benefit in Patients With Gastric Cancer. Ann Surg 2018;267:311-8. [Crossref] [PubMed]
  29. Lyu L, Yao J, Wang M, et al. Overexpressed Pseudogene HLA-DPB2 Promotes Tumor Immune Infiltrates by Regulating HLA-DPB1 and Indicates a Better Prognosis in Breast Cancer. Front Oncol 2020;10:1245. [Crossref] [PubMed]
  30. Chae YK, Arya A, Iams W, et al. Current landscape and future of dual anti-CTLA4 and PD-1/PD-L1 blockade immunotherapy in cancer; lessons learned from clinical trials with melanoma and non-small cell lung cancer (NSCLC). J Immunother Cancer 2018;6:39. [Crossref] [PubMed]

(English Language Editor: A. Kassem)

Cite this article as: Shu C, Liu L, Chen X, Xue J, Fei J, Wang J, Yang X, Peng Q, Yuan H. ncRNA-mediated low expression of P2RY14 correlates with poor prognosis and tumor immune infiltration in ovarian carcinoma. Ann Transl Med 2023;11(1):10. doi: 10.21037/atm-22-6120

Download Citation