frontiers in Oncology
Check for updates
SPRY4-AS1, A Novel Enhancer RNA, Is a Potential Novel Prognostic Biomarker and Therapeutic Target for Hepatocellular Carcinoma
Mu Ye 1,2, Sheng Wang1,3, Jing-bo Qie 1,3* and Pei-long Sun 1,2*
1 Center for Tumor Diagnosis and Therapy, Jinshan Hospital, Fudan University, Shanghai, China, 2 Department of General Surgery, Jinshan Hospital, Fudan University, Shanghai, China, 3 Institutes of Biomedical Sciences, Fudan University, Shanghai, China
A growing number of evidence have demonstrated the involvement of enhancer RNAs (eRNAs) in tumor progression. However, the possible functions of eRNAs in hepatocellular carcinoma (HCC) remain largely unclear. Our present research aimed to screen critical eRNAs and to further delve into the clinical significance of eRNAs in HCC patients. In this study, we identified 124 prognosis-related eRNAs by analyzing The Cancer Genome Atlas (TCGA) datasets. Among them, SPRY4 antisense RNA 1 (SPRY4-AS1) may be a key eRNA involved in HCC progression. SPRY4 was a regulatory target of SPRY4-AS1. High SPRY4-AS1 expression was associated with poor prognosis of HCC patients. Kyoto Encyclopedia of Genes and Genomes (KEGG) assays revealed that the mainly enriched biological process included Human papillomavirus infection, Hippo signaling pathway, and Proteoglycans in cancer. Besides, RT-PCR and immunohistochemical staining confirmed SPRY4-AS1 as an overexpressed eRNA in HCC specimens. The pan-cancer assays revealed that SPRY4-AS1 was associated with glioblastoma multiforme (GBM), adrenocortical carcinoma (ACC), brain lower grade glioma (LGG) and mesothelioma (MESO). Positive associations were observed between SPRY4-AS1 and SPRY4 (its target gene) in 16 tumor types. Collectively, our findings reveal a novel eRNA SPRY4-AS1 for HCC progression and suggest that SPRY4-AS1 may be a potential biomarker and therapeutic target for HCC.
Keywords: enhancer RNA, hepatocellular carcinoma, biomarker, therapeutic target, SPRY4-AS1
INTRODUCTION
Hepatocellular carcinoma (HCC) is a high-prevalence cancer worldwide, especially in China (1, 2). Its incidence as well as mortality rate will still be on the increase over the forthcoming decades (3). Although more and more optimizations in the treatment strategies have resulted in distinct improvement of the survival of HCC patients, HCC remains a highly deadly neoplasm, with offensive malignancy and late diagnosis (4, 5). More optimized treatment strategies could be
OPEN ACCESS
Edited by: Fan Feng, The 302nd Hospital of PLA, China
Reviewed by:
Bo-an Li, The Fifth Medical Center of Chinese PLA General Hospital, China Ming-bo Cao, First Affiliated Hospital of Zhengzhou University, China Qiyu Jiang, Fifth Medical Center of the PLA General Hospital, China Xu Che, Chinese Academy of Medical Sciences and Peking Union Medical College, China
*Correspondence: Jing-bo Qie
Pei-long Sun
Specialty section:
This article was submitted to Molecular and Cellular Oncology, a section of the journal Frontiers in Oncology
Received: 27 August 2021
Accepted: 09 September 2021 Published: 04 October 2021
Citation:
Ye M, Wang S, Qie J-b and Sun P-l (2021) SPRY4-AS1, A Novel Enhancer RNA, Is a Potential Novel Prognostic Biomarker and Therapeutic Target for Hepatocellular Carcinoma. Front. Oncol. 11:765484. doi: 10.3389/fonc.2021.765484
provided by prognosis prediction (6). In addition, these prognostic biomarkers may be potential targets for therapeutic intervention, thereby improving clinical outcomes (7). At present, as a commonly used biomarker, a-fetoprotein (AFP) displayed a relatively low sensitivity and specificity (8). Therefore, other valuable biomarkers must be investigated.
Long noncoding RNA (lncRNA) is in the molecule form with no possibility of being translated into a protein, and its nucleotides are over 200 (9). An example is enhancer RNA (eRNA) transcribed from putative enhancer regions (10). Its functions are still unclear, except that enhancer transcription is considered a noisy by-product of transcription machinery (11). Over the past years, substantial research has evidenced the crucial role of eRNAs in gene regulation for embryonic development and diseases (12, 13). Among thousands of eRNAs found in the cells of human beings, many display a mediating role in activating target genes (14). Specifically, eRNAs have been found to participate in the progression of cancers (15, 16). In recent years, the tumor-related functions of eRNAs have been reported in HCC. For instance, lncRNA SNHG17 was found to promote the proliferation and metastasis of HCC cells via regulating miR-3180-3p/RFX1 (17). LncRNA TTN-AS1 was observed to intensify sorafenib resistance in HCC through the modulation of miR-16-5p/cyclin E1 axis. However, what potential functions eRNA have in HCC remained unclear even though it is key to gene transcription control.
This research intended to identify prognostic eRNAs involved in HCC progression. Our findings suggest a strategy for targeting eRNA as a potential biomarker for the treatment of HCC patients.
MATERIALS AND METHODS
Identification of Functional eRNAs in Hepatocellular Carcinoma Using Bioinformatics Technology
The Cancer Genome Atlas (TCGA) database (https://portal.gdc. cancer.gov/) was used to collect the data of 33 cancer types, including clinical and survival information, as well as RNA expression profiles. Here, 33 cancer types were included: pancreatic adenocarcinoma (PAAD), prostate adenocarcinoma (PRAD), rectum
Abbreviations: ACC, Adrenocortical Carcinoma; BLCA, Bladder Urothelial Carcinoma; BRCA, Breast Invasive Carcinoma; CESC, Cervical Squamous Cell Carcinoma and Endocervical Adenocarcinoma; CHOL, Cholangiocarcinoma; COAD, Colon Adenocarcinoma; DLBC, Lymphoid Neoplasm Diffuse Large B- cell Lymphoma; ESCA, Oesophageal Carcinoma; GBM, Glioblastoma Multiforme; HNSC, Head and Neck Squamous Cell Carcinoma; KICH, Kidney Chromophobe; KIRC, Kidney Renal Clear Cell Carcinoma; KIRP, Kidney Renal Papillary Cell Carcinoma; LAML, Acute Myeloid Leukaemia; LGG, Brain Lower Grade Glioma; LIHC, Liver Hepatocellular Carcinoma; LUAD, Lung Adenocarcinoma; LUSC, Lung Squamous Cell Carcinoma; MESO, Mesothelioma; OV, Ovarian Serous Cystadenocarcinoma; PAAD, Pancreatic Adenocarcinoma; PCPG, Pheochromocytoma and Paraganglioma; PRAD, Prostate Adenocarcinoma; READ, Rectum Adenocarcinoma; SARC, Sarcoma; SKCM, Skin Cutaneous Melanoma; STAD, Stomach Adenocarcinoma; TGCT, Testicular Germ Cell Tumours; THCA, Thyroid Carcinoma; THYM, Thymoma; UCEC, Uterine Corpus Endometrial Carcinoma; UCS, Uterine Carcinosarcoma; UVM, Uveal Melanoma.
adenocarcinoma (READ), skin cutaneous melanoma (SKCM), stomach adenocarcinoma (STAD), colon adenocarcinoma (COAD), lymphoid neoplasm diffuse large B-cell lymphoma (DLBC), oesophageal carcinoma (ESCA), glioblastoma multiforme (GBM), bladder urothelial carcinoma (BLCA), head and neck squamous cell carcinoma (HNSC), kidney chromophobe (KICH), kidney renal clear cell carcinoma (KIRC), kidney renal papillary cell carcinoma (KIRP), brain lower grade glioma (LGG), liver hepatocellular carcinoma (LIHC), and uterine carcinosarcoma (UCS). The relevant eRNA information was obtained from the eRNA source literatures. By the use of the limma R software package, the eRNA expressing matrix and HCC patients’ survival information were combined. Kaplan-Meier assays were applied to determine the survival-associated eRNAs. Based on the median expression of each eRNA, patients were divided into high or low groups, low expression and high expression. By the use of Spearman’s correlation assays, candidate key eRNAs that were related to survival and target genes in HCC were obtained.
Gene Enrichment Analysis
For exploring mechanisms whereby eRNAs may influence the clinical outcome of HCC patients, we used the “cluterprofile” R package to perform functional enrichment analyses of Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) (18, 19). The “ggplot2” R package was applied to identify the significant pathways. To avert a high false discovery rate (FDR) because of complex comparisons, we examined q-values. FDR- adjusted q-value <0.25 and p-value <0.05 were considered statistically significant.
Tumor Samples
Ten paired HCC and matched normal non-tumor specimens were collected from 10 HCC patients (ages between 27 and 78 years old) at the Jinshan Hospital, Fudan University for RT-qPCR assays between June 2020 and May 2021. None of the patients received preoperative chemotherapy or radiotherapy. All cases were histopathologically confirmed as HCC by two independent pathologists. All tissues were preserved at -80℃ until use. The Clinical Research & Ethics Committee of our hospital provided its approval of the present research. All patients offered informed consent in a written form.
Quantitative Real-Time PCR Analysis
From all tumor and normal specimens, TRIzol reagent (Invitrogen, Suzhou, Jiangsu, China) was used to extract total
| Name | Bidirectional primer sequence (5'-3') |
|---|---|
| SPRY4: Forward primer | TCTGACCAACGGCTCTTAGAC |
| SPRY4: Reverse primer | GTGCCATAGTTGACCAGAGTC |
| SPRY4-AS1: Forward primer | GACCTGCTCGACCTGACCCTC |
| SPRY4-AS1: Reverse primer | CCACCTCGAACCACAATTCA |
| GAPDH: Forward primer | GGAGCGAGATCCCTCCAAAAT |
| GAPDH: Reverse primer | GGCTGTTGTCATACTTCTCATGG |
GAPDH, glyceraldehyde 3-phosphate dehydrogenase; SPRY4-AS1, SPRY4 antisense RNA 1.
RNA under the manufacturer’s instructions. RNA was reverse transcribed to cDNA from 2 µg of total RNA by Reverse Transcription Kits (Takara, Hangzhou, Zhejiang, China). qRT-PCR assays were carried out by a protocol from Power SYBR Green (Takara, Hangzhou, Zhejiang, China). The calculation and normalization of the relative expressions of SPRY4 antisense RNA 1 (SPRY4-AS1) and SPRY4 were achieved with 2-44Ct methods relative to glyceraldehyde 3-phosphate dehydrogenase (GAPDH). Specific primer sequences are shown in Table 1.
Immunohistochemical Staining
Immunohistochemical staining was performed according to a previous study (20). Briefly, the cancer samples were fixed, paraffin-embedded, and sectioned. To obtain antigen retrieval, the sections were then incubated for 30 min. Then, for blocking endogenous peroxidase activities, we cultured the specimens 3% H2O2 for 20 min. Subsequently, primary antibodies targeting SPRY4 were applied to cultivate the sections of tumor specimens. Finally, after phosphate buffered saline (PBS) was applied to wash the sections, they are incubated with horseradish peroxidase (HRP)-conjugated secondary antibody followed by a light microscope for the photograph. Primary antibodies targeting SPRY4 were provided by Abcam (Cambridge, UK).
SPRY4-AS1 Expression and Its Prognostic Value in Pan-Cancer
Firstly, the R limma package was used to obtain the expression data of SPRY4-AS1 and the corresponding target gene SPRY4 in pan-cancer, followed by combination with the survival data of pan-cancer. Samples were divided, based on the median value of SPRY4-AS1 expression, into two groups, high expression and low expression, followed by Kaplan-Meier methods determining the survival difference. The associations between SPRY4-AS1 and SPRY in pan-cancer were examined by the use of Spearman’s coefficient.
Statistical Analysis
Through the SPSS statistical software package (standard version 18.0, SPSS Inc., Chicago, IL, USA) or RStudio (v.3.6.1, RStudio Team, 2016; RStudio Integrated Development for R; RStudio, Inc., Boston, USA), we conducted statistical analyses. Student’s t- test or chi-square test was applied to the analysis of group differences using SPSS. Combining Kaplan-Meier curves with log-rank tests with the “survival” R package, we analyzed survival results. Pearson correlation was used to analyze genes co- expressed with eRNAs in various tumors using RStudio. A p < 0.05 denoted statistical significance.
RESULTS
Screening the Key eRNA in Hepatocellular Carcinoma
To screen prognostic eRNAs, we analyzed TCGA datasets including 371 HCC cases (Table 2) and performed Kaplan-
| Characteristic | n | Proportion (%) | |
|---|---|---|---|
| Total | 371 | 100 | |
| Median follow-up (days) | 557 (1-3,675) | 371 | 100 |
| Age | 59.5 ±13.0 | 371 | 100 |
| Sex | Male | 251 | 67.7 |
| Female | 120 | 32.3 | |
| Tumor grade | I | 54 | 14.6 |
| II | 178 | 48.0 | |
| III | 122 | 32.9 | |
| IV | 12 | 3.2 | |
| Unknown | 5 | 1.3 | |
| Stage | I | 174 | 46.9 |
| II | 86 | 23.2 | |
| III | 84 | 22.6 | |
| IV | 5 | 1.3 | |
| Unknown | 22 | 5.9 | |
| T stage | I | 183 | 49.3 |
| II | 94 | 25.3 | |
| III | 80 | 21.6 | |
| IV | 12 | 3.2 | |
| Unknown | 2 | 0.5 | |
| N stage | Without metastasis | 254 | 68.5 |
| With metastasis | 4 | 1.1 | |
| Unknown | 113 | 30.5 | |
| M stage | Without metastasis | 268 | 72.2 |
| With metastasis | 4 | 1.1 | |
| Unknown | 99 | 26.7 |
HCC, hepatocellular carcinoma.
Meier methods, identifying 124 prognosis-related eRNAs with p < 0.05 (Table S1). To further identify 124 eRNAs with a distinct association with their target genes related to HCC, we performed Spearman’s correlation. SPRY4-AS1 exhibited the highest cor value. The results of Kaplan-Meier methods for SPRY4-AS1 suggested that the high-expression group had a shorter overall survival compared with the other group (p = 0.0002; Figure 1A). Besides, SPRY4 was found in Figure 1B to positively correlate with SPRY4-AS1 (R = 0.56, p < 2.2e10-15). Our group further examined the associations between HCC patients’ clinical features and the SPRY4-AS1 expression. We observed SPRY4-AS1 expression was not associated with age (p = 0.079; Figure 2A), However, we observed that the expression of SPRY4-AS1 in HCC was distinctly linked to gender (p = 0.00016; Figure 2B), grade (Figure 2C), stage (Figure 2D), and cancer status (p = 0.0039; Figure 2E).
Gene Enrichment Analysis
A distinct association was found between 672 transcripts and SPRY4-AS1 (p < 0.05), including SPRY4. For the 1,407 target genes, KEGG pathway analysis and GO enrichment analysis revealed possible mechanisms responsible for the SPRY4-AS1 function. In biological process (BP), the association of terms with actin filament organization and the regulation of intracellular transport, cellular protein localization, cell cycle process, and intracellular transport was found; in Cellular Component (CC), they were related to focal adhesion, cell-substrate junction, cell leading edge, intrinsic component of organelle membrane, and
exocytic vesicle. In Molecular Function (MF), term enrichment mainly involved actin binding, cadherin binding, GTPase regulator activity, and actin filament binding (Figure 3A). Figure 3B presents the 30 most important pathways. As suggested by KEGG assays, the BP mainly included Human papillomavirus infection, Hippo signaling pathway, and Proteoglycans in cancer (adjusted p < 0.001).
The Distinct Upregulation of SPRY4-AS1 and SPRY4 in Hepatocellular Carcinoma Specimens
To determine the expression of SPRY4-AS1 and SPRY4 in HCC specimens, we applied 10 HCC specimens and 10 non-tumor specimens to RT-PCR. As shown in Figures 4A, B, we observed that compared to the non-tumor specimens, the expression of
A
B
SPRY4-AS1 level
6
high
low
R= 0.56, p < 2.2e-16
1.00-
Survival probability
0.75
4
0.50
SPRY4
0.25
p=0.002
0.00
2
..
0
1
2
3
4
5
6
7
8
9
10
SPRY4-AS1 level
Time(years)
high
184
113
66
39 54
30
23
16
5
3
1
0
low
184
151
77
36
20
13
4
3
3
1
0
1
2
3
4
5
6
7
8
9
10
0.0
0.5
1.0
Time(years)
SPRY4-AS1
1.5
A
B
C
Age
8 ⇐ 86 E
>65
Gelder
Ę
female
male
Grade EI G1 E 02 E G3 G4
0.079
0.00016
0.46
1.5
1.5
0.32
0.32
2
0.081
SPRY4-AS1 expression
SPRY4-AS1 expression
SPRY4-AS1 expression
0.0018
1.0
1.0
0.0089
1
0.5
0.5
0.0
0.0
0
D
⇐ 65
>65
Age
female
Gender
male
G1
G2
G3
G4
Grade
D
E
Stage E stage l E stage il stage Ill E stage IV
Cancer_status = TUMOR FREE WITH TUMOR
0.7
0.0039
0.82
0.61
0.55
0.75
SPRY4-AS1 expression
2
0.0037
SPRY4-AS1 expression
0.016
0.50
1.
0.25
0
0.00
stage I
stage II
stage III
stage IV
TUMOR FREE
WITH TUMOR
Stage
Cancer_status
FIGURE 2 | The associations between SPRY4 antisense RNA 1 (SPRY4-AS1) and clinical features. (A) Age. (B) Gender. (C) Grade. (D) Stage. (E) Cancer status.
A
B
actin filament organization
regulation of intracellular transport-
Human papillomavirus infection
positive regulation of cellular protein localization-
Hippo signaling pathway
positive regulation of cell cycle process-
Proteoglycans in cancer
positive regulation of intracellular transport
positive regulation of nervous system development
€
Hepatocellular carcinoma
regulation of protein localization to membrane-
Fc gamma R-mediated phagocytosis
Cytoskeleton-dependent
Adrenergic signaling in cardiomyocytes
Count
cytokinesis
substrate adhesion-dependent cell spreading
Cushing syndrome
Count
5.0
Cell cycle
focal adhesion-
7.5
10
cell-substrate junction-
O
20
Oocyte meiosis
10.0
cell leading edge
intrinsic component of organelle membrane-
30
Glycerophospholipid metabolism
12.5
Glucagon signaling pathway
15.0
exocylic vesicle
lamelipodium
8
qvalue
ECM-receptor interaction
melanosome
Hypertrophic cardiomyopathy
pigment granule
0.05
pvalue
0.10
GnRH signaling pathway
lateral plasma membrane-
actomyosin
0.15
Dilated cardiomyopathy
0.01
0.02
actin binding
0.20
Prostate cancer
003
cadherin binding
Melanogenesis
0.04
GTPase regulator activity
Ether lipid metabolism
actin filament binding
small GTPase binding-
Central carbon metabolism in cancer
calcium ion transmembrane transporter activity-
编
SNARE interactions in vesicular transport
calcium channel activity
Sphingolipid metabolism
cytoskeletal anchor activity-
ligand-gated calcium channel activity
Pyrimidine metabolism
chloride channel inhibitor activity
Basal cell carcinoma
0.01
0.02
0.03
0.04
0.05
0.02
0.03
0.04
0.05
0.06
GeneRatio
GeneRatio
A
B
C
Relative SPRY4 expression
Relative SPRY4-AS1 expression
10
p<0.01
8
p<0.01
Relative SPRY4 expression
8
r=0.8009
8-
6
6
p=0.0054
6
4
4
4
2
2.
2
0
0
0
Normal (n=10)
Tumor (n=10)
Normal (n=10)
Tumor (n=10)
0
2
4
6
8
10
Relative SPRY4-AS1 expression
D
Tumor specimens
Normal specimens
SPRY4-AS1 and SPRY4 was distinctly increased in HCC specimens. Correlation assays revealed that SPRY4-AS1 expressions were positively associated with SPRY4 in 10 HCC tissues (Figure 4C). Furthermore, immunohistochemistry showed that the expressions of SPRY4 at the protein level were lower in non-tumor specimens than those in tumor specimens (Figure 4D).
Pan-Cancer Verification
The purpose of survival and correlation analyses was to demonstrate SPRY4-AS1’s prognostic value in pan-cancer. SPRY4-AS1 was found to associate with survival in ACC (Figure 5A), GBM (Figure 5B), LGG (Figure 5C), and MESO (Figure 5D). Besides, we observed that SPRY4-AS1 and SPRY4 were associated with 16 types of tumors (Figure 6).
A
Cancer: ACC
B
Cancer: GBM
SPRY4-AS1 level
+
high
+
low
SPRY4-AS1 level
+
high
+ low
1.00-
1.00
Survival probability
Survival probability
0.75
0.75
0.50
0.50
0.25
p=0.003
0.25
p=0.006
SPRY4-AS1 level
0.00
0.00
0
1
2
3
4
5
6
7
8
9
10
11
12
SPRY4-AS1 level
0
1
2
3
4
5
6
7
8
Time(years)
Time(years)
high
39
35
24
8
4
ON
high
83
34
9
2
00
0
0
0
0
low
40
40
34
18
12
18
00 00
3 0 - chw
3
3
JON
26
16
12
4
1
low
84
47
16
6
2
1
1
O
0
1
2
3
4
5
6
7
8
9
10
11
12
0
1
2
3
4
5
6
7
8
Time(years)
Time(years)
C
Cancer: LGG
D Cancer: MESO
SPRY4-AS1 level
+
high
+
low
SPRY4-AS1 level
+
high
+
low
1.00
1.00-
Survival probability
Survival probability
0.75
0.75
0.50
0.50
0.25
p=0.004
0.25
p=0.031
0.00
0.00
0
2
4
6
8
10
12
14
16
18
20
SPRY4-AS1 level
0
1
2
3
4
5
6
7
8
SPRY4-AS1 level
Time(years)
Time(yea’s)
high - low
262
113
262
141
46
29
13
12
5
27
Nw
1
49
7
4
0
00
0
0
high
low
42
42
28
11
3
28
19
11
1
7
0
5
0
3
0
1
0
14
0
2
4
6
8
10
12
14
16
18
20
0
1
2
3
4
5
6
7
8
Time(years)
Time(years)
FIGURE 5 | Survival curves for SPRY4 antisense RNA 1 (SPRY4-AS1) in pan-cancer, including ACC (A), GBM (B), LGG (C), MESO (D).
DISCUSSION
HCC is a widely known malignant tumor in the world, and its incidence in many countries shows an upward trend (21). In the early stages, many HCC patients could achieve a favorable long- term survival after operative treatments and adjuvant therapies (22, 23). However, it is hard for patients with metastasis to receive operation, and their clinical outcomes remained very poor (24). Over the past decades, as the understanding of molecular mechanisms involved in HCC progression is deepened, various targeted therapies for HCC patients were developed (25, 26). Sensitive biomarkers are necessary to guide the application of targeted therapies. However, the related biomarkers were limited in clinical practice. In addition, the early detection also contributed to a favorable prognosis for HCC patients. Over the past 10 years, more and more studies focused on the huge potential of ncRNAs used as novel biomarkers for HCC patients (27, 28).
A growing number of studies have reported the involvement of eRNAs in several tumors. For instance, FOXP4-AS1, a prognosis-related eRNA, was reported to be lowly expressed in ovarian cancer specimens and predicted a shorter 5-year overall
survival of ovarian cancer patients (29). The eRNA SMAD7e was found to promote the proliferation and metastasis of bladder cancer (30). Located in the tissue-specific enhancer of a tumor suppressor gene EMX2, eRNA EMX2OS was lowly expressed in kidney renal clear cell carcinoma specimens and associated with poor prognosis of tumor patients (31). However, in HCC, the effects of eRNAs were rarely reported.
In this study, we analyzed TCGA datasets and identified a novel HCC-related prognostic eRNA SPRY4-AS1. To date, the effects of SPRY4-AS1 in tumors have not been investigated. Moreover, SPRY4-AS1 was found to be the most critical eRNA candidate sequence in HCC, which regulated SPRY4. The expression of SPRY4-AS1 was found in clinical correlation analysis to vary with cancer status, gender, stage, and grade. Interestingly, we observed that the expression of SPRY4-AS1 in female patients was distinctly higher than that in male patients. The reason was unclear. More studies are needed to explore its reason. KEGG pathway enrichment results concluded the possible role of SPRY4-AS1 in the survival outcome of HCC patients. This could be achieved by Proteoglycans in cancer and Hippo signaling pathway. Previously, dysregulation of the Hippo pathway has been recognized in a variety of human cancers, including HCC (32, 33).
Cancer: ACC
Cancer: BLCA
Cancer: BRCA
Cancer: CESC
R=0.56, p <2.2e-16
=
6
R=0.41, p <2.2e=16
6
6
R=0.67.p=9.7e-12
.
6
R=0.67. p<2.2e-16
..
4
4
SPRY4
SPRY4
4
SPRY4
%
SPRY4
4
2
2
2
2
·
.
0
.
0.0
0.4
08
0
1.O
1
SPRY4-AS1
1.2
0.0
0.5
SPRY4-AS1
1.5
20
0
2
3
SPRY4-AS1
0.0
0.5
SPRY4-AS1
1.0
Cancer: ESCA
Cancer: GBM
Cancer: KIRC
Cancer: KIRP
5
R=0.44. p= 1e-08.
6
R=0.5. p <2.2e-16
R=0.45, p <2.2€
6
6
R = 0165, p <2.2e-16
4
4-
SPRY4
SPRY4
4
SPRY4
4
SPRY4
3
2
2
:
2
2
1
·
0.00
0.25
0.50
SPRY4-AS1
0.75
00
0.3
0.6
SPRY4-AS1
0.9
1.2
00
0.5
SPRY4-AS1
1.0
0.0
0.5
SPRY4-AS1
10
1.5
Cancer: LGG
Cancer: LIHC
Cancer: LUAD
Cancer: LUSC
R=0.8, p < 2.2e-16
6
R=0.44, p < 2.2e-16
R=0.47. p<2.2e-16
R=0.56, p <2.2e-16
8
5
6
4
4
SPRY4
4
SPRY4
SPRY4
SPRY4
A
:
2
2
2
2
N
1
A
0
0.0
0.5
Yo
SPRY4-AS1
1.5
0.0
0.5
SPRY4-AS1
1.0
1.5
00
0.5
1.0
SPRY4-AS1
1.5
2.0
0.0
0.5
SPRY4-AS1
1.0
Cancer: OV
Cancer: PCPG
Cancer: PRAD
Cancer: SARC
8
7
10.0
R=0.64, p< 2.2e-16
6
R = 0.42, p <2.2e=16
R=0.67. p<2.2e-16
6
6
R=0.73, p<2.2e=16
7.5
5
:
4
SPRY4
SPRY4
SPRY4
SPRY4
4
5.0
:
4
2
25
2
3
0.0
.
2
0.0
05
1.0
SPRY4-AS1
1.5
2.0
00
0.3
0.6
1.0
SPRY4-AS1
09
00
0.5
SPRY4-AS1
1.0
1.5
0.0
0.5
SPRY4-AS1
1.5
20
FIGURE 6 | The association between SPRY4 antisense RNA 1 (SPRY4-AS1) and SPRY4 in pan-cancer.
The activity of Hippo signaling pathways was confirmed to modulate the proliferation and metastasis of HCC cells (34, 35). Furthermore, SPRY4-AS1 expressions in surgical specimens were verified by RT-PCR and immunohistochemistry, revealing an association with survival in ACC, GBM, LGG, and MESO. The above findings suggested that SPRY4-AS1 may be a novel prognostic biomarker for HCC patients.
Importantly, because this predictive model was based on the data from TCGA set, it remained unknown whether SPRY4-AS1 had similar predictive power beyond molecular subtypes in HCC
patients of other Chinese hospitals. In addition, we did not unveil the molecular pathway through which SPRY4-AS1 exerted its function; this needs a further in-depth investigation.
CONCLUSION
To conclude, a novel prognostic eRNA signature for HCC was examined in this paper. The research results may promote the
prognostic evaluation of HCC and guide further study on targeted therapies in HCC.
DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/ Supplementary Material.
ETHICS STATEMENT
The studies involving human participants were reviewed and approved by Jinshan Hospital, Fudan University. The patients/ participants provided their written informed consent to participate in this study.
AUTHOR CONTRIBUTIONS
MY, J-BQ, and P-LS contributed to the conception. MY and SW contributed to the design and revision of the article. MY and J-BQ contributed to the analysis and interpretation of data. All authors contributed to the article and approved the submitted version.
FUNDING
This work was supported in part by the Shanghai Municipal Health Commission (201740041).
SUPPLEMENTARY MATERIAL
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2021. 765484/full#supplementary-material
REFERENCES
1. Siegel RL, Miller KD, Jemal A. Cancer Statistics, 2018. CA: Cancer J Clin (2018) 68:7-30. doi: 10.3322/caac.21442
2. Xu XF, Xing H, Han J, Li ZL, Lau WY, Zhou YH, et al. Risk Factors, Patterns, and Outcomes of Late Recurrence After Liver Resection for Hepatocellular Carcinoma: A Multicenter Study From China. JAMA Surg (2019) 154:209-17. doi: 10.1001/jamasurg.2018.4334
3. Grandhi MS, Kim AK, Ronnekleiv-Kelly SM, Kamel IR, Ghasebeh MA, Pawlik TM. Hepatocellular Carcinoma: From Diagnosis to Treatment. Surg Oncol (2016) 25:74-85. doi: 10.1016/j.suronc.2016.03.002
4. Sim HW, Knox J. Hepatocellular Carcinoma in the Era of Immunotherapy. Curr Problems Cancer (2018) 42:40-8. doi: 10.1016/j.currproblcancer.2017.10.007
5. Lu J, Qian J, Yin S, Zhou L, Zheng S, Zhang W. Mechanisms of RNA N(6)- Methyladenosine in Hepatocellular Carcinoma: From the Perspectives of Etiology. Front Oncol (2020) 10:1105. doi: 10.3389/fonc.2020.01105
6. Piñero F, Dirchwolf M, Pessoa MG. Biomarkers in Hepatocellular Carcinoma: Diagnosis, Prognosis and Treatment Response Assessment. Cells (2020) 9 (6):1370. doi: 10.3390/cells9061370
7. Rebouissou S, Nault JC. Advances in Molecular Classification and Precision Oncology in Hepatocellular Carcinoma. J Hepatol (2020) 72:215-29. doi: 10.1016/j.jhep.2019.08.017
8. Sauzay C, Petit A, Bourgeois AM, Barbare JC, Chauffert B, Galmiche A, et al. Alpha- Foetoprotein (AFP): A Multi-Purpose Marker in Hepatocellular Carcinoma. Clinica chimica acta; Int J Clin Chem (2016) 463:39-44. doi: 10.1016/j.cca.2016.10.006
9. Kopp F, Mendell JT. Functional Classification and Experimental Dissection of Long Noncoding RNAs. Cell (2018) 172:393-407. doi: 10.1016/j.cell.2018.01.011
10. Sartorelli V, Lauberth SM. Enhancer RNAs are an Important Regulatory Layer of the Epigenome. Nat Struct Mol Biol (2020) 27:521-8. doi: 10.1038/s41594- 020-0446-0
11. Tippens ND, Vihervaara A, Lis JT. Enhancer Transcription: What, Where, When, and Why? Genes Dev (2018) 32:1-3. doi: 10.1101/gad.311605.118
12. Bose DA, Berger SL. eRNA Binding Produces Tailored CBP Activity Profiles to Regulate Gene Expression. RNA Biol (2017) 14:1655-9. doi: 10.1080/ 15476286.2017.1353862
13. Ye R, Cao C, Xue Y. Enhancer RNA: Biogenesis, Function, and Regulation. Essays Biochem (2020) 64:883-94. doi: 10.1042/EBC20200014
14. Jiao AL, Slack FJ. RNA-Mediated Gene Activation. Epigenetics (2014) 9:27-36. doi: 10.4161/epi.26942
15. Wu M, Shen J. From Super-Enhancer Non-Coding RNA to Immune Checkpoint: Frameworks to Functions. Front Oncol (2019) 9:1307. doi: 10.3389/fonc.2019.01307
16. Zhang Z, Lee JH, Ruan H, Ye Y, Krakowiak J, Hu Q, et al. Transcriptional Landscape and Clinical Utility of Enhancer RNAs for eRNA-Targeted Therapy in Cancer. Nat Commun (2019) 10:4562. doi: 10.1038/s41467-019-12543-5
17. Ma T, Zhou X, Wei H, Yan S, Hui Y, Liu Y, et al. Long Non-Coding RNA SNHG17 Upregulates RFX1 by Sponging miR-3180-3p and Promotes Cellular Function in Hepatocellular Carcinoma. Front Genet (2020) 11:607636. doi: 10.3389/fgene.2020.607636
18. Kanehisa M, Furumichi M, Tanabe M, Sato Y, Morishima K. KEGG: New Perspectives on Genomes, Pathways, Diseases and Drugs. Nucleic Acids Res (2017) 45:D353-d61. doi: 10.1093/nar/gkw1092
19. Kanehisa M. Toward Understanding the Origin and Evolution of Cellular Organisms. Protein Sci Publ Protein Soc (2019) 28:1947-51. doi: 10.1002/pro.3715
20. Sun F, Wang J, Sun Q, Li F, Gao H, Xu L, et al. Interleukin-8 Promotes Integrin ß3 Upregulation and Cell Invasion Through PI3K/Akt Pathway in Hepatocellular Carcinoma. J Exp Clin Cancer Res (2019) 38:449. doi: 10.1186/ s13046-019-1455-x
21. Wang W, Wei C. Advances in the Early Diagnosis of Hepatocellular Carcinoma. Genes Dis (2020) 7:308-19. doi: 10.1016/j.gendis.2020.01.014
22. Wang T, Zhang KH. New Blood Biomarkers for the Diagnosis of AFP- Negative Hepatocellular Carcinoma. Front Oncol (2020) 10:1316. doi: 10.3389/fonc.2020.01316
23. Shiani A, Narayanan S, Pena L, Friedman M. The Role of Diagnosis and Treatment of Underlying Liver Disease for the Prognosis of Primary Liver Cancer. Cancer control J Moffitt Cancer Center (2017) 24:1073274817729240. doi: 10.1177/1073274817729240
24. Amorim J, França M, Perez-Girbes A, Torregrosa A, Marti-Bonmatí L. Critical Review of HCC Imaging in the Multidisciplinary Setting: Treatment Allocation and Evaluation of Response. Abdominal Radiol (New York) (2020) 45:3119-28. doi: 10.1007/s00261-020-02470-1
25. Shariati M, Meric-Bernstam F. Targeting AKT for Cancer Therapy. Expert Opin Investig Drugs (2019) 28:977-88. doi: 10.1080/13543784.2019.1676726
26. Li J, Chen F, Cona MM, Feng Y, Himmelreich U, Oyen R, et al. A Review on Various Targeted Anticancer Therapies. Targeted Oncol (2012) 7:69-85. doi: 10.1007/s11523-012-0212-2
27. Chen X, Yang T, Wang W, Xi W, Zhang T, Li Q, et al. Circular RNAs in Immune Responses and Immune Diseases. Theranostics (2019) 9:588-607. doi: 10.7150/thno.29678
28. Beermann J, Piccoli MT, Viereck J, Thum T. Non-Coding RNAs in Development and Disease: Background, Mechanisms, and Therapeutic Approaches. Physiol Rev (2016) 96:1297-325. doi: 10.1152/physrev. 00041.2015
29. Hua T, Tian YJ, Wang RM, Zhao CF, Kong YH, Tian RQ, et al. FOXP4-AS1 is a Favorable Prognostic-Related Enhancer RNA in Ovarian Cancer. Bioscience Rep (2021) 41(4):BSR20204008. doi: 10.1042/BSR20204008
30. Che W, Ye S, Cai A, Cui X, Sun Y. CRISPR-Cas13a Targeting the Enhancer RNA- SMAD7e Inhibits Bladder Cancer Development Both In Vitro and In Vivo. Front Mol Biosci (2020) 7:607740. doi: 10.3389/fmolb.2020.607740
31. Jiang H, Chen H, Wan P, Song S, Chen N. Downregulation of Enhancer RNA EMX2OS is Associated With Poor Prognosis in Kidney Renal Clear Cell Carcinoma. Aging (Albany NY) (2020) 12:25865-77. doi: 10.18632/aging.202151
32. Kim W, Khan SK, Gvozdenovic-Jeremic J, Kim Y, Dahlman J, Kim H, et al. Hippo Signaling Interactions With Wnt/B-Catenin and Notch Signaling Repress Liver Tumorigenesis. J Clin Invest (2017) 127:137-52. doi: 10.1172/ JCI88486
33. Taha Z, Janse van Rensburg HJ, Yang X. The Hippo Pathway: Immunity and Cancer. Cancers (2018) 10(4):94. doi: 10.3390/cancers10040094
34. Hu H, Xu L, Luo SJ, Xiang T, Chen Y, Cao ZR, et al. Retinal Dehydrogenase 5 (RHD5) Attenuates Metastasis via Regulating HIPPO/YAP Signaling Pathway in Hepatocellular Carcinoma. Int J Med Sci (2020) 17:1897-908. doi: 10.7150/ ijms.46091
35. Huang Z, Zhou JK, Wang K, Chen H, Qin S, Liu J, et al. PDLIM1 Inhibits Tumor Metastasis Through Activating Hippo Signaling in Hepatocellular
Carcinoma. Hepatol (Baltimore Md) (2020) 71:1643-59. doi: 10.1002/ hep.30930
Conflict of Interest: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Copyright @ 2021 Ye, Wang, Qie and Sun. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.