frontiers in Endocrinology
Check for updates
Competing Endogenous RNA Network Analysis Reveals Pivotal ceRNAs in Adrenocortical Carcinoma
Weiwei Liang1* and Fangfang Sun2
1 Department of Endocrinology, The Second Affiliated Hospital, Zhejiang University School of Medicine, Hangzhou, China,
2 Cancer Institute (Key Laboratory of Cancer Prevention and Intervention, China National Ministry of Education; Key
Laboratory of Molecular Biology in Medical Sciences), The Second Affiliated Hospital, Zhejiang University School of Medicine, Hangzhou, China
Objective: To construct ceRNA network and identify pivotal competing endogenous RNAs (ceRNAs) in adrenocortical carcinoma (ACC) using ceRNA network analysis.
Methods: The RNA sequencing expression data of 77 ACCs in TCGA were obtained from GEPIA. Cancer specific ceRNAs, cancer specific microRNAs (miRNAs), and cancer specific messenger RNAs (mRNAs) were identified. The interaction of cancer specific miRNAs with cancer specific ceRNAs and cancer specific mRNAs were predicted. CeRNA network was constructed and visualized by Cytoscape 3.7.0 software. The genes in ceRNA network regulated GO terms and regulated pathways were performed by function analysis. Survival analysis of pivotal ceRNAs was performed for the pivotal lncRNAs.
Result: Twenty-eight cancer specific ceRNAs, 149 cancer specific miRNAs, and 104 mRNAs were identified. CeRNA network was constructed including 10 ceRNAs, 35 miRNAs, and 34 mRNAs. The genes in ceRNA network regulated GO terms and were classified into three groups: cellular component (CC), molecular function (MF), and biological process (BP). The genes in ceRNA network regulated the following pathways: leukocyte transendothelial migration, and proteoglycans in cancer. Survival analysis showed that CTB-63M22.1 and RP1-241P17.4 were significantly associated with ACC patient disease free survival and overall survival.
Conclusion: This study has constructed ceRNA networks in ACC. The study provides a set of pivotal ceRNAs for future investigation into the molecular mechanisms.
Keywords: adrenocortical carcinoma, non-coding RNA, ceRNA network analysis, lncRNA, ceRNA
INTRODUCTION
Cancers are often associated with aberrant transcriptomes. Recent analyses of the human transcriptome have revealed that 2% of the genome encode protein-coding transcripts even though over three-quarters of the genome are transcribed (1, 2). As genes do not function in isolation, they can be grouped into “networks” based on their interactions. A novel hypothesis, known as the ceRNA hypothesis presented by Salmena et al. (3), has emerged in recent years. CeRNAs, acting as miRNA sponges, can bind miRNAs competitively to indirectly regulate mRNAs using shared microRNA response elements. In theory, ceRNAs refer to all transcripts that may become
OPEN ACCESS
Edited by: Gabriella Castoria, Second University of Naples, Italy
Reviewed by: Vincenzo Pezzi, University of Calabria, Italy Dimitra Karagkouni, University of Thessaly, Greece
*Correspondence: Weiwei Liang helenliangww@zju.edu.cn
Specialty section:
This article was submitted to Cancer Endocrinology, a section of the journal Frontiers in Endocrinology
Received: 23 December 2018
Accepted: 26 April 2019
Published: 15 May 2019
Citation:
Liang W and Sun F (2019) Competing Endogenous RNA Network Analysis Reveals Pivotal ceRNAs in Adrenocortical Carcinoma. Front. Endocrinol. 10:301. doi: 10.3389/fendo.2019.00301
the targets of miRNA such as long non-coding RNA (lncRNA), pseudogene RNA, and circular RNA. The function of miRNAs is relatively well-understood. MiRNAs can interact with target mRNAs to regulate the expression of the gene (4). Studies have shown that miRNAs are involved in the initiation and progression of cancers (4). While studies on lncRNAs and transcribed pseudogenes are still in their infancy, some studies showed that lncRNAs and pseudogenes played key regulatory roles in cancer biology (5). Recently, investigations about the ceRNA networks provide a better understanding of tumorigenesis (3, 6). This hypothesis laid the foundations for an important discovery regarding cross talk between both coding and non-coding RNAs through miRNAs.
Though an increasing number of miRNAs (7, 8) and lncRNAs (9) have been discovered to play a role in endocrine cancer. Our knowledge about the ceRNA network in endocrine cancer is still limited. The regulation landscape of ceRNA crosstalk in cancers is still largely unknown.
In this study, we focused on adrenocortical carcinoma (ACC). ACC is a rare endocrine malignancy, often with an unfavorable prognosis (10, 11); 5-years overall survival is below 40% (11). According to a report based on the Surveillance, Epidemiology, and End Results (SEER) database, the incidence of ACC is ~0.72 per million cases per year leading to 0.2% of all cancer deaths in the United States (12). The recent identification of the genomic alterations had provided some molecular mechanisms (13), but ceRNA crosstalk in ACC is still unknown.
In this study, we have taken full advantage of the rich data from the TCGA consortium. The aberrant expression profiles of RNA in ACC were filtered. The aberrant ceRNA network was constructed. This paper represented a significant leap forward in understanding the biological functions of ceRNAs in endocrine cancers.
MATERIALS AND METHODS
Data Sources
The RNA sequencing expression data of 77 ACC patients in TCGA were obtained from GEPIA, a web server for cancer, normal gene expression profiling, and interactive analyses (http://gepia.cancer-pku.cn/index.html) (14). The dataset was strict to ACC, with [log2FC] > 3 and p-value < 0.01. ANOVA was used as differential methods. Both over-expressed and under- expressed expression data were downloaded.
Construction of ceRNA Network of ACC
The RNA expression data downloaded from GEPIA were annotated by R package RefSeq mRNA and RefSeq ncRNA. We then used R package transcript biotype to double check.
First, we identified cancer specific ceRNA, cancer specific miRNA, and cancer specific mRNA. The traditional view was that lncRNA functioned as a ceRNA. Recent articles suggested that pseudogenes also have ceRNA function (15). Therefore, in this paper lncRNA and pseudogenes were both identified as potential
A
B
NEAT1
MIR49THG
PYY2
H19
SLC11A1
GFBP3
8100A9
CYB561A3
ADHIB
is-min-505-55
Nu-miA-182-5p
THB81
UBE2C
hsa-miR-218-
hsa.mA-454.3p
amiR.93-30
BEX1
STEAP4
hsa-miR-1968-5phsa-let-7d-5p
sa-miR 3016
ba-niR-1962:50
hsa-min-17-4p
hsh-101-71-1-3p
hsa-miR-1068-5g
tam-miR-25a-3
sa-miR. 18-2
sa-miA-182.5p
ha@-miR-278-5p
968-le1-7d-6p
Fa-miRt-998-5
ha@-miA-508-3p
sam R-29b-2-50
CDH2
hos-miR-1308-3p
hsa-miR-156-60
haa-mil-106h-5p
taniR.19a.Jp
hsa-miR.156-3
a-miR-150-50
ba-miR. 424-5
hsa-let-71-1-3p
CYP11B2
bsa-miR-93-50
hsa-miA-29b-ap
ODOCTY
ILIALI
samiR.221-50
hoa-mIR-29c-
han-min-200c-
ha-mill-372-Jp
a-miF-15b-52
hsa-miR-222-5p
PHYHP
810049
PPS4Y1
sa-miR-29c-3p
CXCL12
UBE2C
hoa-miR-424-50
POGPRA
A
ADHIB
PYY2
MIND4P12
hsa-miR.296-3p
CCL2
LUG7L2
MMP2
SUGCT
EFIAY
hsa-miR-93-3p
RAPGEF4
hsa-miR-2Bb-2-5p
GR514
ABUIM1
RAB34
RAB34
RP11-649E7,
GTB-63M22
BHMT2
CDIQ
KGF2
PHYHIP
DOXJY
CCL2
KDM60
OFBP’S
13FBPS
THBS1
has-min-130a-ap
1
CYP1182
KLHDCBA
a-miR-2Ba-3p
SLC11A1
H19
EFIAY
INMT
CDKNIC
AVPRIA
BEX1
CYB561A3
hsa-miR-301a/p
RF11-316M1-12
sa-miR.28b-5p
GAB14
hau-miR-3016
MMP2
ILTRL1
NEAT!
KDMED
C
hoa-miR-454-3p
ARAP1-A81
risa-miR-218-5p
BUGCT
AP11-316M1.12
MTND4P12
CXOL12
RP1-241P17.4
AP1-241P17
RP11-649E7.5
POGFRA
CTB-63M22.1
hsa-miR.93-5p
ARAP1-AS1
MIR49THQ
hsa-miR-Job-3p
BHMT2
hsa-miA-99a-95
hsa-miR-29b-295p
bsa-miR-16-2-3p
hsa-let-7d.5p
sa-miR-150-5p hsa-let-71-1.3p
KLHDC8A
hsa-miR-372-3p
INMT
hsa-miR.608-5p
ADH1B
ABLIM1
hsa-miR-16b-Gp
hsa-miR-29b-3p
-miR-29c-3p
hsa-miR-18a-3p
Maa-miR-222-5p hsa-miR-93-3p
AVPR1A
hau-miR-20de.sp
STEAP4
IGF2
hsa-miR-221-5p
ÇOKNIC
hsa-miR-424-5p
hsa-miR-29a-3p
hea-min-15b-3p
haa-mi-150-5
hsa-miR-372.3p
hea-miR-221-5
haa-miR-508-5p
RAB34
65a-miA-222.5p
BEX1
hsa-mint-188-Jp
LUC7L2
hsa-miR-31-5p
haa-miR-93-5p
bsa-miR-27b-5p
hsa-miR-16-2-3g
nsa.miA-31-86ª miA-19a.3p
hsa-miR-106b.5p
hsa-miR-17-5p
has-miR-1968-5p
THBS1
IGFBPS
RAPGEF4
RPS4Y1
RPS4Y1
510049
CYB561A3
DDX3Y
13FBPS
MMP2
IGF2
ABLIM1
COKNIC
IGFBIP3
LUCTL2
CXCL12
CDH2
BHMT2
PHYHIP
RAPGEF4
PDGFRA
UBE2C
KLHDC8A
AVPAIA
FIGURE 1 | (A) ceRNA network in ACC. Diamond in yellow represented cancer specific ceRNA. Ellipse in gray represented cancer specific mRNA. Rectangle in blue represented cancer specific mRNA. (B) Down-regulated ceRNAs regulated ceRNA network. (C) Up-regulated ceRNAs regulated ceRNA network.
GO:0060326~cell chemotaxis -
☒
GO:0051897~positive regulation of protein kinase B signaling -
☒
☐
GO:0050729~positive regulation of inflammatory response -
☒
Count
☐
☐
3.0
GO:0044267~cellular protein metabolic process-
☒
☒ 3.5
☒
4.0
Biological Process
☐
☐
GO:0043410~positive regulation of MAPK cascade -
☒
☒
4.5
☐
☐
☐
☒
5.0
GO:0009612~response to mechanical stimulus -
☒
☐
☐
☐
-log10(PValue)
GO:0007568~aging -
☒
☐
☐
☐
☐
☐
☐
2.5
GO:0006955~immune response-
☒
☐ ☐
☐ ☐
☐
2.0
GO:0006954~inflammatory response -
☒
1.5
☐
GO:0006413~translational initiation -
☒
☐
GO:0001666~response to hypoxia-
☒
10
15
20
25
Fold Enrichment
FIGURE 2 | Gene ontology analysis for biological process of genes involved in the ceRNA network in ACC.
GO:0042567~insulin-like growth factor ternary complex -
☒
Count
2
☒ 3
☒
4
Cellular Component
GO:0016942~insulin-like growth factor binding protein complex -
☒ 5
☒
☒
6
☒ 7
☒ 8
-log10(PValue)
GO:0005615~extracellular space-
☒
2.2
2.0
1.8
GO:0005576~extracellular region -
☒
0
100
200
300
Fold Enrichment
FIGURE 3 | Gene ontology analysis for cellular component of genes involved in the ceRNA network in ACC.
cancer specific ceRNAs. Cancer specific miRNAs were retained from Oncomir (http://www.oncomir.org) (16). A list of miRNAs closely associated with patient survival or tumor development were retrieved with p < 0.05. mRNAs with absolute [log2FC] > 3 and p-value < 0.01 retained from GEPIA were identified as cancer specific mRNAs.
The ceRNA network was then constructed through the following steps: (1) Prediction of cancer specific ceRNA-cancer specific miRNA interaction: miRcode (http://www.mircode. org) and DIANA-LncBase v2 (www.microrna.gr/LncBase) (17) was used to predict cancer specific ncRNA-cancer specific miRNA interaction. (2) Cancer specific miRNA-cancer specific mRNA interaction: miRtarBase (http://mirtarbase. mbc.nctu.edu.tw/php/index.php) (18) was used to retrieve
cancer specific miRNA-cancer specific miRNA interaction. (3) CeRNA network construction and visualization; Cytoscape 3.7.0 software (19) was utilized to construct and visualize ceRNA network.
Functional Analysis
The biological function and pathway of the genes involved in the ceRNA network would demonstrate instructive information. Thus, we used Gene ontology analysis (GO) to identify characteristic biological attributes. Kyoto Encyclopedia of Genes and Genomes pathway (KEGG) enrichment analysis was performed to identify functional attributes. p < 0.05 was set as the cut-off criterion.
GO:0031995~insulin-like growth factor Il binding -
Count
2.00
2.25
2.50
Molecular Function
2.75
3.00
GO:0031994~insulin-like growth factor I binding -
-log10(PValue)
2.5
2.0
GO:0001968~fibronectin binding -
60
80
100
120
Fold Enrichment
FIGURE 4 | Gene ontology analysis for molecular function of genes involved in the ceRNA network in ACC.
| Term | Genes | P-value | Fold enrichment | FDR |
|---|---|---|---|---|
| GO:0009612~response to mechanical stimulus | CCL2, THBS1, CXCL12 | 5.29E-03 | 26.68 | 7.28 |
| GO:0060326~cell chemotaxis | CCL2, PDGFRA, CXCL12 | 6.38E-03 | 24.22 | 8.72 |
| GO:0050729~positive regulation of inflammatory response | CCL2, IL1RL1, S100A9 | 7.99E-03 | 21.57 | 10.80 |
| GO:0043410~positive regulation of MAPK cascade | IGF2, CDH2, IGFBP3 | 9.76E-03 | 19.44 | 13.05 |
| GO:0051897~positive regulation of protein kinase B signaling | IGF2, THBS1, IGFBP5 | 1.05E-02 | 18.74 | 13.93 |
| GO:0006413~translational initiation | EIF1AY, DDX3Y, RPS4Y1 | 2.63E-02 | 11.49 | 31.64 |
| GO:0007568~aging | CDKN1C, CCL2, IGFBP5 | 3.70E-02 | 9.54 | 41.61 |
| GO:0044267~cellular protein metabolic process | IGF2, IGFBP3, MMP2, IGFBP5 | 1.32E-03 | 17.79 | 1.86 |
| GO:0001666~response to hypoxia | CCL2, THBS1, CXCL12, MMP2 | 3.85E-03 | 12.20 | 5.34 |
| GO:0006954~inflammatory response | SLC11A1, CCL2, S100A9, THBS1, CXCL12 | 4.96E-03 | 6.92 | 6.85 |
| GO:0006955~immune response | SLC11A1, CCL2, IL1RL1, THBS1, CXCL12 | 7.17E-03 | 6.23 | 9.75 |
Survival Analysis of Pivotal lncRNAs
The ceRNAs in ACC ceRNA network were identified as pivotal ceRNAs. CeRNAs’ correlations with patient survival were featured in GEPIA. Available TCGA patient survival data were used for Kaplan-Meier survival analysis and to generate overall and disease-free survival plots. Group-cutoff was set as the median, the hazards ratio was calculated based on Cox PH Model and the 95% confidence interval was added as a dotted line. P < 0.05 was set as the cut-off criterion.
RESULTS
Identification of Cancer Specific ceRNAs, miRNAs, and mRNAs
Seventy-seven ACC patient data in TCGA were analyzed. Twenty-eight cancer specific ceRNAs (10 down regulated in ACC and 18 up regulated in ACC, Supplementary Table 1), 149 cancer specific miRNAs and 104 cancer specific mRNAs were identified (Supplementary Tables 2,3).
CeRNA Network Construction
To construct the ceRNA network, we next predicted the interaction of cancer specific miRNAs with cancer specific ceRNAs and cancer specific mRNAs. We overlapped the predicted targets of cancer specific ceRNAs with cancer specific miRNAs and overlapped the predicted targets of cancer specific miRNAs with cancer specific mRNAs. Finally, the ceRNA regulatory network was constructed (Figure 1A), including 10 ceRNAs, 35 miRNAs, and 34 mRNAs. In the down-regulated ceRNA regulatory network, there were 4 ceRNAs, 31 miRNAs, and 33 mRNAs involved (Figure 1B). And in the up-regulated ceRNA regulatory network, there were six ceRNAs, 24 miRNAs, and 25 mRNAs involved (Figure 1C).
Gene Ontology and Pathway Enrichment Analyses
We performed gene ontology analysis of the genes in the ceRNA network (Figures 2-4, Tables 1-3). The result showed that related to CC, the genes in the ceRNA network were mainly
enriched in extracellular region. Related to MF, the genes in the ceRNA network were mainly enriched in fibronectin binding. Related to BP, the genes in the ceRNA network were mainly enriched in inflammatory response and immune response.
We then performed pathway enrichment analysis of the genes in the ceRNA network and identified the ceRNAs regulated pathway. Functional analysis showed that the ceRNAs regulated pathways included leukocyte transendothelial migration, and proteoglycans in cancer.
Survival Analysis of Pivotal lncRNAs
We used GEPIA to perform survival analysis of pivotal ceRNAs (Figures 5, 6). p < 0.05 was set as the cut-off criterion. The results showed that CTB-63M22.1 and RP1-241P17.4 were significantly associated with ACC patient disease free survival and overall survival.
DISCUSSION
Salmena et al. (3) firstly proposed a hypothesis on how mRNAs and IncRNAs “talk” to each other using miRNA as letters of a new language. The last few years have nominated ceRNA hypothesis as a novel layer of gene regulation. This ceRNA regulatory network greatly expands the functional genetic information and will aid in deciphering the underlying pathogenesis of human polygenic disease at the post-transcriptional level. Previously, Xiao et al. (20) have analyzed the gene expression data of ACC in The Cancer Genome Atlas (TCGA) and identified the potential prognostic mRNA biomarkers. To our knowledge, this is the first study to investigate the ceRNA network in ACC.
In this study, we extracted RNA expression profile from cases of ACC from TCGA. Twenty-eight cancer specific ceRNAs, 149 cancer specific miRNAs, and 104 mRNAs were identified. A ceRNA network was then constructed, consisting of 10 ceRNAs, 35 miRNAs, and 34 miRNAs. The Kyoto Encyclopedia of Genes and Genomes (KEGG), and the Gene Ontology (GO) databases were used to identify biological and functional attributes of the genes in the ceRNA network. The ceRNA network-regulated GO
| Term | Genes | P-value | Fold enrichment | FDR |
|---|---|---|---|---|
| GO:0016942~insulin-like growth factor binding protein complex | IGFBP3, IGFBP5 | 5.09E-03 | 379.67 | 5.29 |
| GO:0042567~insulin-like growth factor ternary complex | IGFBP3, IGFBP5 | 6.79E-03 | 284.75 | 6.98 |
| GO:0005615~extracellular space | CCL2, S100A9, IGF2, THBS1, CXCL12, IGFBP3, MMP2 | 2.40E-02 | 2.96 | 22.77 |
| GO:0005576~extracellular region | CCL2, S100A9, IGF2, THBS1, CXCL12, IGFBP3, MMP2, IGFBP5 | 1.65E-02 | 2.83 | 16.20 |
| Term | Genes | P-value | Fold enrichment | FDR |
|---|---|---|---|---|
| GO:0031995~insulin-like growth factor Il binding | IGFBP3, IGFBP5 | 1.51E-02 | 127.89 | 16.02 |
| GO:0031994~insulin-like growth factor I binding | IGFBP3, IGFBP5 | 2.25E-02 | 85.26 | 23.04 |
| GO:0001968~fibronectin binding | THBS1, IGFBP3, IGFBP5 | 1.10E-03 | 59.02 | 1.26 |
A
Overall Survival
B
Overall Survival
C
Overall Survival
D
Overall Survival
E
Overall Survival
0
Low H19 TPM
2
High H19 TPM
Low MIR497HG TPM
1
High MIR497HG TPM
Low NEAT1 TPM
2
£
Low PYY2 TPM
9
High NEAT1 TPM
High PYY2 TPM
Low CTB-63M22.1 TPM
Logrank p=0.51
Logrank p=0.98
Logrank p=0.62
Logrank p=0.19
High CTB-63M22.1 TPM
3
HR(high)=1.3
3
HR(high)=0.99
HR(high)=1.2
Logrank p=0.0015
p(HR)=0.5
13
HR(high)=1.7
p(HR)=0.62
0
HR(high)=3.6
p(HR)=0.97
00
p(HR)=0.2
p(HR)=0.0028
Percent survival
n(high)=38
n(high)=38
n(high) 38
n(high)=38
0
n(low)=38
9
n(low)=38
n(low)=38
n(high)=38
3
n(low)=38
0
4
n(low)=38
5
3
5
à
:
2
S
3
S
8
8
00
:
8
8
0
50
100
150
0
50
100
150
0
50
100
150
0
50
100
150
0
50
100
150
Months
Months
Months
Months
Months
F
Overall Survival
G
Overall Survival
H
Overall Survival
I
Overall Survival
J Overall Survival
2
Low RP1-241P17.4 TPM
2
0
High RP1-241P17.4 TPM
Low TPM
High RP11-316M1.12 TPM
Low TPM
9
Low ARAP1-AS1 TPM
2
Low MIND4P12 TPM
Logrank p=0.014
HR(high)=2.6
Logrank p=0.83
High RP11-649E7.5 TPM
Logrank p=0:21
High ARAP1-AS1 TPM
Logrank p=0.81
High MTND4P12 TPM
Logrank p=0.87
3
p(HR)=0.017
8
HR(high)=0.92
0
HR(high)=1.7
p(HR)=0.83
p(HR)=0.21
0
·HR(high)=1.1
HR(high)=1.1
n(high)=38
p(HR)=0.81
2
p(HR)=0.87
Percent survival
n(high)=38
n(high)=38
n(high)=38
48
n(low)=38
n(high)=38
0
n(low)=38
0
n(low)=38
n(low)=38
48
n(low)=38
8
0.4
5
à
3
8
3
y
CH
2
2
0
8
:
g
18
0
50
100
150
0
50
100
150
0
50
100
150
0
50
100
150
0
50
100
150
Months
Months
Months
Months
Months
A
Disease Free Survival
B
Disease Free Survival
C
Disease Free Survival
D Disease Free Survival
E
Disease Free Survival
0
Low H19 TPM
0
High H19 TPM
Low MIR497HG TPM
High TPM
9
Low TPM
0
High NEAT1 TPM
Low PYY2 TPM
9
High
Low CTB-63M22.1 TPM
High
Logrank p=0.95
HR(high)=0.98
Logrank p=0.63
Logrank p=0.89
Logrank p=0.086
Logrank-p=7:3e-05
0.8
p(HR)=0.96
0
HR(high)=0.85
p(HR)=0.64
0
HR(high)=1
HR(high)=1.8
HR(high)=4.1
p(HR)=0.9
0
p(HR)=0.089
13
p(HR)=0.00023
Percent survival
n(high) 38
n(low)=38
n(high)=38
n(high)=38
n(high)=38
n(high)=38
0.6
n(low)=38
3
n[low)=38
n(low)=38
8
n[low)=38
8
CD
0
5
5
3
0
0
0
%
S
0
0
8
8
8
0
50
100
150
0
50
100
150
0
50
100
150
0
50
100
150
0
50
100
150
Months
Months
Months
Months
Months
F
Disease Free Survival
G
Disease Free Survival
H
Disease Free Survival
I Disease Free Survival
J Disease Free Survival
0
Low TPM
0
Low TPM
0
Low RP11-649E7.5 TPM
9
Low
High ARAP1-AS1 TPM
0
High RP1-241P17.4 TPM
High RP11-316M1.12 TPM
Logrank p=0.22
High TPM
Low MTND4P12 TPM
Logrank p=0.013
Logrank p=0.13
Logrank p=0.17
High
0.8
HR(high)=2.4
8
HR(high)=0.66
Logrank p=0.16
p(HR)-0.015
P(HR)=0.22
8
HR(high)=1.7
1
HR(high)=0.63
p(HR)=0.14
p(HR)=0.17
8
HR(high)=1.6
n(high) 38
n(high)=38
n(high) 38
n(high)=38
p(HR)=0.17
Percent survival
n(low)=38
n(low)=38
n(low)=38
n(low)=38
n(high) 38
0.6
0
8
0
8
n(low)=38
0
à
5
5
à
0
8
y
3
S
8
8
0
0
0
0
50
100
150
0
50
100
150
0
50
100
150
0
50
100
150
0
50
100
150
Months
Months
Months
Months
Months
FIGURE 6 | (A-J): Disease free survival plot of 10 pivotal ceRNAs using 77 ACC patient data from TCGA database.
terms were classified into three groups. The ceRNA network was mainly enriched in the following pathways: leukocyte transendothelial migration and proteoglycans in cancer. Survival analysis showed two potential prognostic associated ceRNAs: CTB-63M22.1 and RP1-241P17.4.
Of the genes that were involved in the ceRNA network, GO functional enrichment analysis showed that these genes were mainly enriched in the extracellular region, fibronectin binding, inflammatory response, and immune response. Perge’s (21) previous study showed that miRNAs in exosome (an extracellular
substance) might have potential diagnostic value in ACC. MiR-31 in serum exosomes possibly correlated with an aggressive profile in endocrine cancer metastasis (8). Interestingly, miR-31 was a miRNA identified in our ceRNA network. Griggs et al. (22) reported that fibronectin, a protein in the extracellular matrix, served as a growth factor delivery system in epithelial- mesenchymal transition of cancers. Host immune response status played an important role in the clinical outcome of patients with cancer. Systemic inflammatory response was related with cancer progression (23, 24).
Of the genes that were involved in ceRNA network, KEGG analysis showed that these genes were mainly enriched in leukocyte transendothelial migration and proteoglycans in cancer. Leukocyte transendothelial migration is generally activated in cancer progression (25). Proteoglycans (PGs), comprising structurally diverse constituents of the extracellular matrix and cell surfaces, have emerged as novel biomarkers and molecular players both within tumor cells and within their microenvironment. Recent studies have demonstrated that the expression of proteoglycans is regulated by miRNAs (26, 27). The ceRNA network identified in our study provided useful clues of further study.
By analyzing TCGA data, 35 cancer specific miRNA were identified in the ceRNA network. Several miRNAs in our list have been investigated in ACC: Doghman et al. (28) identified that miR-99a was differentially expressed in childhood adrenocortical tumors. Mir-99a might involve in tumor cell growth via regulating mTOR signaling. Bimpak et al. (29) found miR-200 was differentially expressed in massive macronodular adrenocortical disease and miR-200b directly targeted matrin 3 expression in an adrenocortical cancer cell line. The understanding of other miRNAs identified in this study is limited, therefore, further molecular studies are require. These ACC-specific miRNAs might become potential biomarkers with specificity in the diagnosis and progression of ACC in future.
In our study, we identified 10 pivotal ceRNAs. Liu’s study showed that H19 expression is maintained at low levels in hormonally active adrenocortical carcinomas, which is consistent with our study (30). There is some evidence indicating that the ceRNAs identified in our study are associated with other cancer types. Xiong et al. (31) reported that NEAT1, a ceRNA, accelerated lung adenocarcinoma deterioration. Teng et al.’s (32) recent study showed that lncRNA ARAP1-AS1, which functioned as a ceRNA, promoted the progression of bladder cancer by regulating the miR-4735-3p/NOTCH2 axis. Studies of other
ceRNAs are limited, but they provide a good direction for future research. Of the 10 pivotal ceRNAs, CTB-63M22.1, and RP1-241P17.4 were associated with tumor prognosis. There are currently no studies on these two ceRNAs, so they are worthy of future research.
Although the findings of our study have important clinical implications, the limitations must also be noted. In addition to lncRNA and pseudogenes, cirRNA and tRNA may also function as ceRN, and their role in ceRNA network needs further research. Further studies with much larger sample sizes and longer follow-up are required to verify our results. The biological roles of ceRNAs in ACC need to be further investigated. In the future, molecular biology methods including qPCR, luciferase reporter systems and co-immunoprecipitation assays are helpful to validate our findings, thus unraveling the molecular mechanisms of ceRNA networks. CeRNA activity is associated with a multitude of parameters such as the number of miRNA binding types residing on coding and non-coding transcripts, and these parameters need further investigation.
In conclusion, this study has constructed ceRNA network in ACC. The study provides a set of pivotal ceRNAs for future investigation into the molecular mechanisms. This study provided novel insights to explore the underlying mechanism of endocrine cancers.
AUTHOR CONTRIBUTIONS
WL drew the figures and wrote the paper. FS analyzed the data.
SUPPLEMENTARY MATERIAL
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo. 2019.00301/full#supplementary-material
REFERENCES
1. Djebali S, Davis CA, Merkel A, Dobin A, Lassmann T, Mortazavi A, et al. Landscape of transcription in human cells. Nature. (2012) 489:101- 8. doi: 10.1038/nature11233
2. Consortium EP. An integrated encyclopedia of DNA elements in the human genome. Nature. (2012) 489:57-74. doi: 10.1038/naturel1247
3. Salmena L, Poliseno L, Tay Y, Kats L, Pandolfi PP. A ceRNA hypothesis: the rosetta stone of a hidden RNA language? Cell. (2011) 146:353- 8. doi: 10.1016/j.cell.2011.07.014
4. Hayes J, Peruzzi PP, Lawler S. MicroRNAs in cancer: biomarkers, functions and therapy. Trends Mol Med. (2014) 20:460- 9. doi: 10.1016/j.molmed.2014.06.005
5. de Oliveira JC, Oliveira LC, Mathias C, Pedroso GA, Lemos DS, Salviano-Silva A, et al. LncRNAs in cancer: another layer of complexity. J Gene Med. (2018) 2018:e3065. doi: 10.1002/jgm.3065
6. Wee LM, Flores-Jasso CF, Salomon WE, Zamore PD. Argonaute divides its RNA guide into domains with distinct functions and RNA-binding properties. Cell. (2012) 151:1055-67. doi: 10.1016/j.cell.2012.10.036
7. Caglar O, Cayir A. Total circulating cell-free miRNA in plasma as a predictive biomarker of the thyroid diseases. J Cell Biochem. (2018) 120:9016- 22. doi: 10.1002/jcb.28173
8. Lima CR, Gomes CC, Santos MF. Role of microRNAs in endocrine cancer metastasis. Mol Cell Endocrinol. (2017) 456, 62-75. doi: 10.1016/j.mce.2017.03.015
9. Qin L, Luo JZ, Tang XL, Han CG. Identification of long noncoding RNA MIR22HG as a novel biomarker in thyroid cancer. Pathol Oncol Res. (2018) 25:703-10. doi: 10.1007/s12253-018-0521-6
10. Erickson LA, Rivera M, Zhang J. Adrenocortical carcinoma: review and update. Adv Anat Pathol. (2014) 21:151- 9. doi: 10.1097/PAP.0000000000000019
11. Jouinot A, Bertherat J. Management of endocrine disease: adrenocortical carcinoma: differentiating the good from the poor prognosis tumors. Eur J Endocrinol. (2018) 178:R215-30. doi: 10.1530/EJE-18-0027
12. Kebebew E, Reiff E, Duh QY, Clark OH, McMillan A. Extent of disease at presentation and outcome for adrenocortical carcinoma: have we made progress? World J Surg. (2006) 30:872-8. doi: 10.1007/s00268-0 05-0329-x
13. Pereira SS, Monteiro MP, Bourdeau I, Lacroix A, Pignatelli D. Mechanisms of endocrinology: cell cycle regulation in adrenocortical carcinoma. Eur J Endocrinol. (2018) 179:R95-110. doi: 10.1530/EJE-17-0976
14. 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-102. doi: 10.1093/nar/gkx247
15. Glenfield C, McLysaght A. Pseudogenes provide evolutionary evidence for the competitive endogenous RNA hypothesis. Mol Biol Evol. (2018) 35:2886- 99. doi: 10.1093/molbev/msy183
16. Wong NW, Chen Y, Chen S, Wang X. OncomiR: an online resource for exploring pan-cancer microRNA dysregulation. Bioinformatics. (2018) 34:713-5. doi: 10.1093/bioinformatics/btx627
17. Paraskevopoulou MD, Vlachos IS, Karagkouni D, Georgakilas G, Kanellos I, Vergoulis T, et al. DIANA-LncBase v2, indexing microRNA targets on non-coding transcripts. Nucleic Acids Res. (2016) 44:D231-8. doi: 10.1093/nar/gkv1270
18. Chou CH, Shrestha S, Yang CD, Chang NW, Lin YL, Liao KW, et al. miRTarBase update 2018, a resource for experimentally validated microRNA-target interactions. Nucleic Acids Res. (2018) 46:D296-302. doi: 10.1093/nar/gkx1067
19. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. (2003) 13:2498-504. doi: 10.1101/gr.1 239303
20. Xiao H, Xu D, Chen P, Zeng G, Wang X, Zhang X. Identification of five genes as a potential biomarker for predicting progress and prognosis in adrenocortical carcinoma. J Cancer. (2018) 9:4484-95. doi: 10.7150/jca.26698
21. Perge P, Butz H, Pezzani R, Bancos I, Nagy Z, Paloczi K, et al. Evaluation and diagnostic potential of circulating extracellular vesicle- associated microRNAs in adrenocortical tumors. Sci Rep. (2017) 7:5474. doi: 10.1038/s41598-017-05777-0
22. Griggs LA, Hassan NT, Malik RS, Griffin BP, Martinez BA, Elmore LW. Fibronectin fibrils regulate TGF-beta1-induced Epithelial-Mesenchymal Transition. Matrix Biol. (2017) 60-61:157- 75. doi: 10.1016/j.matbio.2017.01.001
23. Soria-Valles C, Lopez-Soto A, Osorio FG, Lopez-Otin C. Immune and inflammatory responses to DNA damage in cancer and aging. Mech Ageing Dev. (2017) 165(Pt A):10-16. doi: 10.1016/j.mad.2016.10.004
24. Dolan RD, McSorley ST, Horgan PG, Laird B, McMillan DC. The role of the systemic inflammatory response in predicting outcomes in patients with advanced inoperable cancer: systematic review and meta-analysis. Crit Rev Oncol Hematol. (2017) 116, 134-46. doi: 10.1016/j.critrevonc.20 17.06.002
25. Enarsson K, Lundin BS, Johnsson E, Brezicka T, Quiding-Jarbrink M. CD4+ CD25high regulatory T cells reduce T cell transendothelial migration in cancer patients. Eur J Immunol. (2007) 37:282-91. doi: 10.1002/eji.200636183
26. Bartel DP. MicroRNAs: target recognition and regulatory functions. Cell. (2009) 136:215-33. doi: 10.1016/j.cell.2009.01.002
27. Lujambio A, Lowe SW. The microcosmos of cancer. Nature. (2012) 482:347- 55. doi: 10.1038/nature10888
28. Duregon E, Rapa I, Votta A, Giorcelli J, Daffara F, Terzolo M. MicroRNA expression patterns in adrenocortical carcinoma variants and clinical pathologic correlations. Hum Pathol. (2014) 45:1555-62. doi: 10.1016/j.humpath.2014.04.005
29. Bimpaki EI, Iliopoulos D, Moraitis A, Stratakis CA. MicroRNA signature in massive macronodular adrenocortical disease and implications for adrenocortical tumourigenesis. Clin Endocrinol. (2010) 72:744-51. doi: 10.1111/j.1365-2265.2009.03725.x
30. Liu J, Kahri AI, Heikkila P, Ilvesmaki V, Voutilainen R. H19 and insulin-like growth factor-II gene expression in adrenal tumors and cultured adrenal cells. J Clin Endocrinol Metab. (1995) 80:492-6. doi: 10.1210/jcem.80.2.7531713
31. Xiong DD, Li ZY, Liang L, He RQ, Ma FC, Luo DZ. The LncRNA NEAT1 accelerates lung adenocarcinoma deterioration and binds to Mir-193a-3p as a competitive endogenous RNA. Cell Physiol Biochem. (2018) 48:905- 18. doi: 10.1159/000491958
32. Teng J, Ai X, Jia Z, Wang K, Guan Y, Guo Y. Long non-coding RNA ARAP1-AS1 promotes the progression of bladder cancer by regulating miR-4735-3p/NOTCH2 axis. Cancer Biol Ther. (2019) 20:552-61. doi: 10.1080/15384047.2018.1538613
Conflict of Interest Statement: 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.
Copyright @ 2019 Liang 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.