Investigation of the Clinical Significance and Prognostic Value of the lncRNA <i>ACVR2B-As1</i> in Liver Cancer

Nie, Yuanyuan; Jiao, Yan; Li, Yanqing; Li, Wei

doi:https://doi.org/10.1155/2019/4602371

BioMed Research International

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusions Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2019 | Article ID 4602371 | https://doi.org/10.1155/2019/4602371

Investigation of the Clinical Significance and Prognostic Value of the lncRNA ACVR2B-As1 in Liver Cancer

Yuanyuan Nie,¹Yan Jiao,²Yanqing Li,³and Wei Li¹

Academic Editor: Fengjie Sun

Received19 Jun 2019

Revised30 Aug 2019

Accepted07 Sept 2019

Published30 Oct 2019

Abstract

A refined liver cancer staging system and effective prognostic prediction can help clinicians make optimized treatment decisions, which is essential in our fight against cancer and for improving the unsatisfying survival rate of liver cancer globally. The prognosis of liver cancer is not only related to tumor status, it is also affected by the patients’ liver functions and the chosen treatment. Currently, several staging systems are being tested. Herein, we analyzed RNA-seq data from the TCGA database and identified a newly annotated lncRNA, ACVR2B-AS1, whose expression is upregulated in liver cancer. Higher ACVR2B-AS1 expression is an independent adverse prognostic factor for overall survival (OS) and relapse-free survival (RFS) in liver cancer patients. Our work suggests that the lncRNA ACVR2B-AS1 could be a candidate biomarker for liver cancer prognosis. Furthermore, ACVR2B-AS1 might serve as a potential therapeutic target, which is a possibility that is worthy of further study.

1. Introduction

Liver cancer is one of the leading causes of cancer-related deaths worldwide [1]. Although numerous resources have been invested into basic research and clinical trials, the survival rate of liver cancer patients remains dismal. For example, in the US, the latest five-year relative survival rate for liver cancer patients is only 18% [2]. This low rate is largely due to the aggressive nature of liver cancer and the fact that most patients diagnosed are already at advanced stages [3]. While the current staging system utilizes both patient assessment and proposed treatment responses, it remains important to refine the risk classification system, which is crucial for optimizing decision making and reducing cancer mortality [1].

Long noncoding RNAs (lncRNAs) represent a class of genes that are not translated into proteins and whose functions have only begun to emerge in the past few years [4, 5]. They are evolutionarily less conserved than coding genes, and it is estimated that the number of lncRNA genes is far larger than that of protein coding genes [5], although their exact abundance is not clear. Increasing evidence has demonstrated that lncRNAs have important biological functions and that their dysregulation is associated with diseases, including cancers [6]. However, the majority of annotated lncRNAs have not been functionally characterized, and their biological significance remains elusive [6].

ACVR2B-AS1 (ACVR2B-antisense RNA1) is a recently annotated lncRNA [7]. It is located on 3p22.2 and is transcribed from the opposite strand of ACVR2B. The function of ACVR2B-AS1 has not been studied. Here, using data extracted from the TCGA database, we identified ACVR2B-AS1 as an independent prognostic factor for liver cancer patient clinical outcomes. This may help us to refine liver cancer prognostic prediction and to further studies of ACVR2B-AS1.

2. Materials and Methods

2.1. Data Mining

RNA-seq data and patients’ clinical characteristics were downloaded from the TCGA database (https://portal.gdc.cancer.gov/). Transcript abundances were quantified using the RNA-Seq data via the Expectation-Maximization (RSEM) software [8].

2.2. Statistical Analysis

All statistical tests were performed in R [9]. Nonparametric Wilcoxon and Kruskal–Wallis tests were used for differential expression analysis among different subgroups. The diagnostic value of ACVR2B-AS1 was estimated using a receiver operating characteristic (ROC) curve drawn with the pROC package [10], and the area under curve (AUC) was calculated as previously described [10]. A cut-off value was determined by utilizing the ROC curve, and the overall population was then divided into two groups (ACVR2B-AS1 High and ACVR2B-AS1 Low) for the subsequent analysis. The associations between ACVR2B-AS1 and patient clinical features were analyzed via Fisher’s exact or Chi-squared tests. Overall survival and relapse-free survival were determined via Kaplan–Meier curves using the survival package in R [11]. The statistical significance of the differences noticed above was calculated using the log-rank test. The prognostic capabilities of ACVR2B-AS1 were assessed via univariate and multivariate Cox regression analysis using the Cox proportional hazard model. The ggplot2 package in R was used for data visualization [12].

3. Results

3.1. Patient Characteristics

RNA-seq data from a total of 371 patients diagnosed with liver cancer were extracted from the TCGA database for analysis. The clinical characteristics of the patients, including age, sex, tumor histological type, grade, stage, and vital status, are listed in Table 1.

3.2. ACVR2B-AS1 Expression Is Upregulated in Liver Cancer

ACVR2B-AS1 levels were significantly elevated in liver cancer specimens compared with their normal counterparts () (Figure 1). A subgroup analysis revealed that ACVR2B-AS1 was differentially expressed in tumors of different histological grades (), and grade IV tumors showed the highest ACVR2B-AS1 levels. Furthermore, by categorizing the patients based on vital status, we found that deceased patients showed higher ACVR2B-AS1 expression levels than did patients who were still alive at the time of sampling (). No significant differential expression was detected in the other subgroups analyzed.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

(k)

(l)

To further examine these findings, we generated an ROC curve to test the diagnostic abilities of ACVR2B-AS1 level for liver cancer diagnosis and histological grade classification. The results of the AUC analysis did not show any significance (Figure 2), indicating that ACVR2B-AS1 alone may not be an accurate diagnostic parameter.

(a)

(b)

(c)

(d)

(e)

3.3. ACVR2B-AS1 Levels Correlate with Patient Vital Status

Consistent with the previous results, correlation analysis confirmed that the vital status of the patients is correlated with ACVR2B-AS1 levels () (Table 2). The correlation between ACVR2B-AS1 levels and histological grade hardly reached statistical significance ().

3.4. ACVR2B-AS1 Is an Independent Adverse Prognostic Factor for OS in Liver Cancer

We next sought to verify the prognostic usefulness of ACVR2B-AS1 in liver cancer. A Kaplan–Meier curve was plotted against different ACVR2B-AS1 levels within the different subgroups (previously defined via the ROC curve against vital status) (Figure 3). The log-rank test revealed that patients with lower ACVR2B-AS1 expression levels had a significantly longer overall survival (OS) than those with higher ACVR2B-AS1 expression levels (). A subgroup analysis showed that ACVR2B-AS1-low patients displayed better OS among all of the subgroups analyzed.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

To better verify these results, we also performed a Cox regression analysis (Table 3). A univariate Cox regression analysis showed that patients with high ACVR2B-AS1 expression levels had significantly poorer outcomes (HR = 2.03, 95% CI (1.4–2.94)) compared with those with low ACVR2B-AS1 expression levels. Other common clinical parameters were also analyzed, and tumor clinical stage, together with T classification, and residual tumor were shown to be unfavorable factors for patient survival. To rule out possible interference among the variables, we performed a multivariate Cox regression analysis, which confirmed that the adverse effects of high ACVR2B-AS1 expression, T classification, and residual tumor remained significant. Together, these data demonstrated that high ACVR2B-AS1 expression (HR = 1.90, 95% CI (1.31–2.76)) is an independent prognostic factor for OS in liver cancer patients.

3.5. High ACVR2B-AS1 Expression Predicts Shorter RFS in Liver Cancer

We next evaluated the prognostic value of ACVR2B-AS1 expression level for relapse-free survival (RFS) (Figure 4). The log-rank test showed that patients with high ACVR2B-AS1 expression levels had significantly shorter RFS (). When we analyzed this parameter within the different categories, we found that although high ACVR2B-AS1 expression levels predicted poorer RFS in patients at histological stage I (), this effect failed to reach statistical significance in stage III patients (). There were no significant differences in RFS between high and low ACVR2B-AS1 patients in either the stage I or stage III groups, respectively. It was interesting to see that when taking sex into consideration, female patients with higher ACVR2B-AS1 expression levels had much shorter RFS than those females with lower expression levels (), although no difference was found in males. In both the young and old subgroups, patients with higher ACVR2B-AS1 expression levels had shorter times before relapse.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

Consistent with the aforementioned findings, a univariate Cox regression analysis showed that high ACVR2B-AS1 expression is an unfavorable prognostic factor for RFS (HR = 1.71, 95% CI [1.21–2.41]) (Table 4). An additional multivariate Cox analysis confirmed that ACVR2B-AS1 is an independent adverse prognostic factor for RFS, and the adjusted hazard ratio was 1.66 ().

4. Discussion

Liver cancer is one of the most life threatening cancers, and most people are diagnosed at later stages when surgical resection is not an option [3]. A refined staging system and more effective prognostic assessment is not only necessary for the selection of the best treatment for individual patients, but also important for designing clinical trials and coordinating the exchange of information between researchers with comparable criteria [1]. The old TNM staging system has gone out of fashion because establishing prognosis in liver cancer patients is highly sensitive to the patients’ liver functions, the chosen treatment, and other factors. Although there are several staging systems currently being tested, most of which include tumor status, liver function, patient status, and treatment responses, a more comprehensive system is urgently needed. Our previous work revealed several RNA biomarkers that are associated with cancer prognosis and can help us refine cancer staging systems [13–18]. Here, using data extracted from the TCGA database, we identified a novel lncRNA, ACVR2B-AS1, whose upregulation is common in liver cancer, and found that higher ACVR2B-AS1 expression predicted poorer OS and RFS in liver cancer patients.

ACVR2B-AS1 is a newly annotated lncRNA. So far, only two studies using TCGA data have mentioned potential roles of ACVR2B-AS1 in cancers. One of these studies that focused on endometrial cancer (UCEC) found that ACVR2B-AS1 was upregulated in cancer samples. A combined lncRNA-focus expression signature (LFES) that included 10 additional lncRNAs functioned as a superior unfavorable prognostic factor (AUC = 0.887) [19]. The other study, however, identified ACVR2B-AS1 as an independent favorable prognostic factor in breast cancer [20]. Our study revealed that ACVR2B-AS1 is an adverse factor in liver cancer. It is possible that ACVR2B-AS1 may have different roles in distinct biological contexts.

Although controversial results have been reported from genome-wide studies indicating possible regulatory roles of ACVR2B-AS1, its functions have not been experimentally dissected [19, 20]. Recent accumulating evidence has indicated that neighboring/overlapping genes tend to show correlated expression in eukaryotic genomes, indicating possible local regulatory mechanisms [21, 22]. Thus, as ACVR2B-AS1 is an intragenic antisense long noncoding RNA that overlaps with the promoter region of the ACVR2B protein coding gene, it is possible that ACVR2B-AS1 might regulate its neighbor gene ACVR2B in a cis manner. For example, ACVR2B-AS1 transcription might disturb the activation of ACVR2B transcription via “transcription interference” [23]. Alternatively, ACVR2B-AS1 activation might promote the transition of the local chromatin into a relatively open state to enhance ACVR2B transcription [24]. ACVR2B encodes a protein called activin receptor type IIB, which is a major receptor of the transforming growth factor beta (TGF-β) signaling pathway with many functions. For example, binding of ACVR2B to its ligand activin can activate the transcription of genes that inhibit muscle growth and is probably responsible for cancer-related cachexia [25, 26]. ACVR2B has also been reported to be a member of the MALAT1/miR-194-5p/ACVR2B signaling axis, which promotes clear cell kidney carcinoma (KIRC) progression [27]. In our study, ACVR2B-AS1 was associated with poorer prognosis, although it remains unknown if and how ACVR2B-AS1 might regulate tumorigenesis and progression via the mechanisms we mentioned above. Further experiments are required to answer these remaining questions.

It was worth noting that liver functions, another prognostic factor, are not included into the multivariate COX regression analysis, since such data are not available in the TCGA database. Nonetheless, our work here has demonstrated the prognostic value of ACVR2B-AS1 in multivariate COX regression analysis and managed to raise the potential possibility of incorporating ACVR2B-AS1 into liver cancer staging. To better testify its independent value as a prognosticator, further studies are encouraged, for example, in another patients cohort where more detailed information is available.

5. Conclusions

Overall, our data reported here identified the lncRNA ACVR2B-AS1 as an independent adverse prognostic factor in liver cancer. Attempts to decipher and identify the clinical significance of the enormous information hidden in RNA-seq data and to characterize the roles of ACVR2B-AS1 are encouraged, which might help refine the staging system in the future, especially at a time when RNA-seq are becoming more prevalent and more data will become available. In addition, more attention should be paid to basic research on ACVR2B-AS1, which might be a potential therapeutic target for liver cancer.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported in part by the National Natural Science Foundation (81670143, 81071920, 30370439, and 81372835) in China.

References

P. R. Galle, A. Forner, and J. M. Llovet, “EASL clinical practice guidelines: management of hepatocellular carcinoma,” Journal of Hepatology, vol. 69, no. 1, pp. 182–236, 2018.
View at: Publisher Site | Google Scholar
R. L. Siegel, K. D. Miller, and A. Jemal, “Cancer statistics, 2019,” CA: A Cancer Journal for Clinicians, vol. 69, no. 1, pp. 7–34, 2019.
View at: Publisher Site | Google Scholar
C.-Y. Liu, K.-F. Chen, and P.-J. Chen, “Treatment of liver cancer,” Cold Spring Harbor Perspectives in Medicine, vol. 5, no. 9, Article ID a021535, 2015.
View at: Publisher Site | Google Scholar
F. Kopp and J. T. Mendell, “Functional classification and experimental dissection of long noncoding RNAs,” Cell, vol. 172, no. 3, pp. 393–407, 2018.
View at: Publisher Site | Google Scholar
S. Djebali, C. A. Davis, A. Merkel et al., “Landscape of transcription in human cells,” Nature, vol. 489, no. 7414, pp. 101–108, 2012.
View at: Publisher Site | Google Scholar
M. Huarte, “The emerging role of lncRNAs in cancer,” Nature Medicine, vol. 21, no. 11, pp. 1253–1261, 2015.
View at: Publisher Site | Google Scholar
National Center for Biotechnology Information, https://www.ncbi.nlm.nih.gov/.
B. Li and C. N. Dewey, “RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome,” BMC Bioinformatics, vol. 12, no. 1, p. 323, 2011.
View at: Publisher Site | Google Scholar
R Foundation for Statistical Computing, Team RDCJC. R.: A Language and Environment for Statistical Computing, vol. 14, R Foundation for Statistical Computing, Vienna, Austria, 2009.
X. Robin, N. Turck, A. Hainard et al., “pROC: an open-source package for R and S+ to analyze and compare ROC curves,” BMC Bioinformatics, vol. 12, p. 77, 2011.
View at: Publisher Site | Google Scholar
T. M. Therneau and P. M. Grambsch, Modeling Survival Data: Extending the Cox Model, vol. 97, Springer, New York, NY, USA, 2000.
H. Wickham, “Ggplot2: elegant graphics for data analysis,” Journal of the Royal Statistical Society, vol. 174, no. 1, pp. 245-246, 2011.
View at: Google Scholar
Y. Jiao, Z. Fu, Y. Li, L. Meng, and Y. Liu, “High EIF2B5 mRNA expression and its prognostic significance in liver cancer: a study based on the TCGA and GEO database,” Cancer Management and Research, vol. 10, pp. 6003–6014, 2018.
View at: Publisher Site | Google Scholar
Y. Jiao, Z. Fu, Y. Li, W. Zhang, and Y. Liu, “Aberrant FAM64A mRNA expression is an independent predictor of poor survival in pancreatic cancer,” PLoS One, vol. 14, no. 1, Article ID e0211291, 2019.
View at: Publisher Site | Google Scholar
Y. Jiao, Y. Li, Z. Lu, and Y. Liu, “High trophinin-associated protein expression is an independent predictor of poor survival in liver cancer,” Digestive Diseases and Sciences, vol. 64, no. 1, pp. 137–143, 2019.
View at: Publisher Site | Google Scholar
Y. Jiao, Y. Li, P. Jiang, W. Han, and Y. Liu, “PGM5: a novel diagnostic and prognostic biomarker for liver cancer,” PeerJ, vol. 7, Article ID e7070, 2019.
View at: Publisher Site | Google Scholar
Y. Jiao, Y. Li, S. Liu, Q. Chen, and Y. Liu, “ITGA3 serves as a diagnostic and prognostic biomarker for pancreatic cancer,” OncoTargets and Therapy, vol. 12, pp. 4141–4152, 2019.
View at: Publisher Site | Google Scholar
Y. Li, Y. Jiao, Z. Fu, Z. Luo, J. Su, and Y. Li, “High miR-454-3p expression predicts poor prognosis in hepatocellular carcinoma,” Cancer Management and Research, vol. 11, pp. 2795–2802, 2019.
View at: Publisher Site | Google Scholar
M. Zhou, Z. Zhang, H. Zhao, S. Bao, and J. Sun, “A novel lncRNA-focus expression signature for survival prediction in endometrial carcinoma,” BMC Cancer, vol. 18, no. 1, p. 39, 2018.
View at: Publisher Site | Google Scholar
B. Pang, Q. Wang, S. Ning et al., “Landscape of tumor suppressor long noncoding RNAs in breast cancer,” Journal of Experimental & Clinical Cancer Research, vol. 38, no. 1, p. 79, 2019.
View at: Publisher Site | Google Scholar
A. T. Ghanbarian and L. D. Hurst, “Neighboring genes show correlated evolution in gene expression,” Molecular Biology and Evolution, vol. 32, no. 7, pp. 1748–1766, 2015.
View at: Publisher Site | Google Scholar
Q. Ning, Y. Li, Z. Wang, S. Zhou, H. Sun, and G. Yu, “The evolution and expression pattern of human overlapping lncRNA and protein-coding gene pairs,” Scientific Reports, vol. 7, no. 1, p. 42775, 2017.
View at: Publisher Site | Google Scholar
P. A. Latos, F. M. Pauler, M. V. Koerner et al., “Airn transcriptional overlap, but not its lncRNA products, induces imprinted Igf2r silencing,” Science, vol. 338, no. 6113, pp. 1469–1472, 2012.
View at: Publisher Site | Google Scholar
K. C. Wang, Y. W. Yang, B. Liu et al., “A long noncoding RNA maintains active chromatin to coordinate homeotic gene expression,” Nature, vol. 472, no. 7341, pp. 120–124, 2011.
View at: Publisher Site | Google Scholar
T. A. Nissinen, J. Hentila, F. Penna et al., “Treating cachexia using soluble ACVR2B improves survival, alters mTOR localization, and attenuates liver and spleen responses,” Journal of Cachexia, Sarcopenia and Muscle, vol. 9, no. 3, pp. 514–529, 2018.
View at: Publisher Site | Google Scholar
X. Gao, P. Zhao, J. Hu et al., “MicroRNA-194 protects against chronic hepatitis B-related liver damage by promoting hepatocyte growth via ACVR2B,” Journal of Cellular and Molecular Medicine, vol. 22, no. 9, pp. 4534–4544, 2018.
View at: Publisher Site | Google Scholar
Y. Ye, F. Zhang, Q. Chen, Z. Huang, and M. Li, “LncRNA MALAT1 modified progression of clear cell kidney carcinoma (KIRC) by regulation of miR-194-5p/ACVR2B signaling,” Molecular Carcinogenesis, vol. 58, no. 2, pp. 279–292, 2019.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2019 Yuanyuan Nie et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

725

Downloads

1273

Citations