Jian Pu1†Ya Zhang2†Anmin Wang2Zebang Qin2Chenyi Zhuo1Wenchuan Li1Zuoming Xu1Qianli Tang1Jianchu Wang1
Huamei Wei3*- 1Department of Hepatobiliary Surgery, Affiliated Hospital of Youjiang Medical University for Nationalities, Baise, China
- 2Graduate College of Youjiang Medical University for Nationalities, Baise, China
- 3Department of Pathology, Affiliated Hospital of Youjiang Medical University for Nationalities, Baise, China
Background: Hepatocellular carcinoma (HCC) is one of the most aggressive malignancies. Increasing evidence revealed that long noncoding RNAs (lncRNAs) were frequently involved in various malignancies. Here, we explored the clinical significances, roles, and mechanisms of lncRNA ADORA2A antisense RNA 1 (ADORA2A-AS1) in HCC.
Methods: The clinical significances of ADORA2A-AS1 in HCC were analyzed using RNA sequencing (RNA-seq) data from The Cancer Genome Atlas (TCGA) project. The expressions of ADORA2A-AS1, Fascin Actin-Bundling Protein 1 (FSCN1), Matrix Metallopeptidase 2 (MMP2), and Baculoviral IAP Repeat Containing 7 (BIRC7) in HCC tissues and cells were measured by qRT-PCR. Cell Counting Kit-8 (CCK-8), 5-ethynyl-2’-deoxyuridine (EdU), caspase-3 activity assay, transwell migration and invasion assays, and xenograft growth and metastasis experiments were performed to evaluate the roles of ADORA2A-AS1 in HCC. RNA pull-down, RNA immunoprecipitation, qRT-PCR, Western blot, and RNA stability assay were performed to elucidate the mechanisms of ADORA2A-AS1 in HCC.
Results: ADORA2A-AS1 was identified as an HCC-related lncRNA, whose low expression was correlated with advanced stage and poor outcome in HCC. Gain- and loss-of functional experiments demonstrated that ADORA2A-AS1 inhibited HCC cell proliferation, induced cell apoptosis, repressed cell migration and invasion, and repressed xenograft growth and metastasis in vivo. Mechanistically, ADORA2A-AS1 competitively bound HuR (Hu Antigen R), repressed the binding of HuR to FSCN1 transcript, decreased FSCN1 transcript stability, and downregulated FSCN1 expression. The expression of FSCN1 was negatively correlated with ADORA2A-AS1 in HCC tissues. Through downregulating FSCN1, ADORA2A-AS1 repressed AKT pathway activation. Functional rescue assays showed that blocking of FSCN1/AKT axis abrogated the roles of ADORA2A-AS1 in HCC.
Conclusion: Low-expression ADORA2A-AS1 is correlated with poor survival of HCC patients. ADORA2A-AS1 exerts tumor-suppressive roles in HCC via binding HuR and repressing FSCN1/AKT axis.
Introduction
Liver cancer is one of the most aggressive malignancies, whose mortality nearly matches the incidence, with estimated 905,677 new cases and 830,180 deaths in 2020 globally (1). Hepatocellular carcinoma (HCC) is the major histological subtype of liver cancer (2). Although some advances have been achieved in the therapy of HCC, post-surgical recurrence or metastasis is still common for HCC patients, which leads to poor outcome of HCC (3, 4). Accordingly, it is urgent to elucidate the molecular mechanisms driving HCC development and establish new molecular targets for successful intervention.
Many genetic and epigenetic aberrations are involved in HCC initiation and progression (5, 6). Previous studies about HCC mainly focused on protein targets (7). However, protein-coding sequences only account for less than 2% of the human genome. About 70%~80% of human genome sequences transcribe into transcripts (8). Thus, noncoding transcripts are much more numerous than protein-coding mRNAs (9, 10). Among noncoding transcripts, long noncoding RNA (lncRNA) is the most numerous, which has over 200 nucleotides in length (11, 12). Accumulating evidence has documented that many lncRNAs were aberrantly expressed in various diseases, particular malignancies (13–16). Furthermore, many lncRNAs were documented to play critical roles in various pathophysiological processes including cancers (17–23). These cancer-related lncRNAs exert multiple oncogenic or tumor-suppressive functions, such as cellular proliferation, apoptosis, migration, invasion, drug resistance, and metastasis (24–26).
As a class of crucial regulatory RNAs, lncRNAs regulate the expression of target genes at transcriptional, posttranscriptional, translational, or posttranslational levels (27, 28). For example, lncRNAs bind and recruit epigenetic modification enzymes to target genes to regulate transcription of target genes (29, 30). LncRNAs directly bind mRNAs and regulate the stability and/or translation of the interacted mRNAs (31). Some lncRNAs were found to directly bind microRNAs and abrogate the repressive roles of microRNAs on mRNA stability and/or translation (32, 33). Other lncRNAs were also reported to bind proteins and change posttranslational modification of the interacted proteins, such as ubiquitination and phosphorylation (34, 35). These various roles and mechanisms of lncRNAs and the numerous number of lncRNAs suggested that further mining cancer-related lncRNAs may provide potential candidates for HCC targeting.
Here, via analyzing the RNA-sequencing (RNA-seq) data of HCC tissues from The Cancer Genome Atlas (TCGA) project, we identified a novel HCC-related lncRNA ADORA2A-AS1. We further investigated clinical significances, biological roles, and mechanisms of action of ADORA2A-AS1 in HCC in detail.
Materials and Methods
Public Datasets and Bioinformatics Analyses
RNA-seq data of HCC tissues and clinicopathological characteristics of these HCCs were retrieved from TCGA project using R2 (https://hgserver1.amc.nl/cgi-bin/r2/main.cgi) and OncoLnc (http://www.oncolnc.org/). Two in silico tools Coding Potential Assessment Tool (CPAT) (http://lilab.research.bcm.edu/) and TestCode (https://www.bioinformatics.org/sms/testcode.html) were used to calculate coding potential of ADORA2A-AS1 (36). The proteins bound to ADORA2A-AS1 were predicted using three in silico tools, RBPmap (mapping binding sites of RNA binding proteins) (http://rbpmap.technion.ac.il/), ATtRACT (a database of RNA binding proteins and associated motifs) (https://attract.cnic.es/#), and RNA-Protein Interaction Prediction (RPISeq) (http://pridb.gdcb.iastate.edu/RPISeq/index.html) (37, 38).
Clinical Subjects
Surgically resected HCC tissues (n = 76) were obtained from the Affiliated Hospital of Youjiang Medical University for Nationalities (Baise, China). Informed consents were acquired from all HCC patients. Clinicopathological characteristics of these 76 HCC cases were presented in Table 1. The Ethics Committee of Affiliated Hospital of Youjiang Medical University for Nationalities approved this study.
Table 1 Correlation between ADORA2A-AS1 expression and clinicopathological characteristics in 76 HCC cases.
Cell Culture and Treatment
Human HCC cell lines SNU-398, SK-HEP-1, and Hep3B were obtained from American Type Culture Collection (ATCC, Manassas, VA, USA). SNU-398 cells were cultured in RPMI 1640 medium (Invitrogen, Carlsbad, CA, USA) added with 10% fetal bovine serum (FBS, Invitrogen). SK-HEP-1 and Hep3B cells were cultured in Eagle’s minimum essential medium (Invitrogen) added with 10% FBS. All cells were authenticated by short-tandem repeat analyses. All cells were routinely detected as mycoplasma-free. Where indicated, cells were treated with 50 µM α-amanitin (Sigma-Aldrich, St. Louis, MO, USA) or 10 µM LY294002 (Selleck, Houston, TX, USA) for indicated time.
Determination of Subcellular Localization of ADORA2A-AS1
RNA fluorescence in situ hybridization (FISH) was performed to detect subcellular localization of ADORA2A-AS1 in SNU-398 cells as previously described (31). Briefly, the ADORA2A-AS1 probes were designed and synthesized by Advanced Cell Diagnostics (Newark, NJ, USA). RNA FISH was undertaken using RNAscope Fluorescent Multiplex Detection Kit (Advanced Cell Diagnostics) following the provided instruction. FISH results were detected using confocal laser scanning microscopy (Leica, Wetzlar, Germany). Furthermore, subcellular localization of ADORA2A-AS1 was evaluated using cytoplasmic and nuclear RNA isolation, followed by quantitative real-time polymerase chain reaction (qRT-PCR). Cytoplasmic and nuclear RNA isolation was undertaken using the PARIS Kit (Invitrogen). qRT-PCR was undertaken as described below.
RNA Extraction and qRT-PCR
Total RNA from HCC tissues and cells was extracted using TRIzol Reagent (Invitrogen). First-strand complementary DNA (cDNA) was generated using the M-MLV Reverse Transcriptase (Invitrogen) and random primers. qRT-PCR was undertaken using the TB Green Premix Ex Taq II (TaKaRa, Tokyo, Japan) on StepOnePlus Real-Time PCR System (Applied Biosystems, Foster City, CA, USA). The sequences of gene-specific primers were as follows: 5′-GGGCTTTCCTACATCCATC-3′ (sense) and 5′-CACAGCAAGTCCATTCCTCT-3′ (anti-sense) for ADORA2A-AS1, 5′-ATCAACCGCCCCATCATC-3′ (sense) and 5′-TGCCCACCGTCCAGTATTT-3′ (anti-sense) for FSCN1, 5′-CCGCAGTGACGGAAAGATG-3′ (sense) and 5′-CCTGGGACAGACGGAAGTT-3′ (anti-sense) for MMP2, 5′-TGGGACCCGTGGGAAGAA-3′ (sense) and 5′-GGCACTTTCAGACTGGACCTC-3′ (anti-sense) for BIRC7, 5′-TGGAGAAATAGTAGATGGC-3′ (sense) and 5′-GGTGAGGAAGTAAAAACAG-3′ (anti-sense) for Metastasis Associated Lung Adenocarcinoma Transcript 1 (MALAT1), 5′-ACACGGACAGGATTGACAGA-3′ (sense) and 5′-GGACATCTAAGGGCATCACA-3′ (anti-sense) for 18S rRNA, 5′-GTCGGAGTCAACGGATTTG-3′ (sense) and 5′-TGGGTGGAATCATATTGGAA-3′ (anti-sense) for glyceraldehyde 3-phosphate dehydrogenase (GAPDH). GAPDH was employed as an internal control. Relative expression levels of transcripts were calculated using the comparative Ct method.
Vectors and Lentivirus Productions
ADORA2A-AS1 full-length sequences were PCR-amplified with the primers 5′-GGGGTACCTTCGCCTCTCGGTTAGC-3′ (sense) and 5′-CGGGATCCAAAGTTTTAAAACATTCCTGTTAAC-3′ (anti-sense). Next, the PCR products were cloned into the Kpn I and BamH I sites of pcDNA3.1(+) vector (Invitrogen) to construct ADORA2A-AS1 overexpression vector pcDNA3.1-ADORA2A-AS1. The PCR products were also cloned into the Kpn I and BamH I sites of pSPT19 vector (Roche, Basel, Switzerland) to construct ADORA2A-AS1 in vitro transcription vector pSPT19-ADORA2A-AS1. Two pairs of cDNA oligonucleotides targeting ADORA2A-AS1 were synthesized and cloned into the shRNA lentivirus expressing plasmid pLV3/H1/GFP&Puro (GenePharma, Shanghai, China). ADORA2A-AS1-specific shRNA lentivirus was produced by GenePharma. Scrambled non-targeting shRNA lentivirus was used as negative control (NC). The sequences of shRNA oligonucleotide were as follows: 5′-GATCCGGATAGAGATCTCAAAGAAGTTTCAAGAGAACTTCTTTGAGATCTCTATCCTTTTTTG-3′ (sense) and 5′-AATTCAAAAAAGGATAGAGATCTCAAAGAAGTTCTCTTGAAACTTCTTTGAGATCTCTATCCG-3′ (anti-sense) for shRNA-ADORA2A-AS1-1, 5′-GATCCGCCCATTGCTTCCAAGCAACATTCAAGAGATGTTGCTTGGAAGCAATGGGCTTTTTTG-3′ (sense) and 5′-AATTCAAAAAAGCCCATTGCTTCCAAGCAACATCTCTTGAATGTTGCTTGGAAGCAATGGGCG-3′ (anti-sense) for shRNA-ADORA2A-AS1-2, 5′-GATCCGTTCTCCGAACGTGTCACGTTTCAAGAGAACGTGACACGTTCGGAGAACTTTTTTG-3′ (sense) and 5′-AATTCAAAAAAGTTCTCCGAACGTGTCACGTTCTCTTGAAACGTGACACGTTCGGAGAACG-3′ (anti-sense) for shRNA-NC.
Stable Cell Line Constructions
To generate HCC cells with ADORA2A-AS1 stable overexpression, ADORA2A-AS1 overexpression vector pcDNA3.1-ADORA2A-AS1was transfected into SNU-398 and SK-HEP-1 cells using Lipofectamine 3000 (Invitrogen). Empty vector pcDNA3.1 was concurrently transfected and used as negative control. Forty-eight hours after transfection, 800 µg/ml neomycin was added. Individual colonies were picked and expanded. The expression of ADORA2A-AS1 in these colonies was measured using qRT-PCR. We chose one colony with successful ADORA2A-AS1 overexpression for each cell line. To generate HCC cells with ADORA2A-AS1 stable silencing, SUN-398 and Hep3B cells were infected with ADORA2A-AS1-specific shRNA lentivirus. Ninety-six hours after infection, 2 µg/ml puromycin was added. Individual colonies were picked and expanded. The expression of ADORA2A-AS1 in these colonies was measured using qRT-PCR. We chose one colony with successful ADORA2A-AS1 downregulation for each shRNA lentivirus and each cell line. To generate HCC cells with ADORA2A-AS1 and FSCN1 overexpression, ADORA2A-AS1-overexpressing SUN-398 cells were infected with FSCN1 overexpression lentivirus (GenePharma). Ninety-six hours after infection, 2 µg/ml puromycin was added. Individual colonies were picked and expanded. The expression of FSCN1 in these colonies was measured using Western blot.
Examination of Cell Proliferation, Apoptosis, Migration, and Invasion
Cell proliferation was measured using Cell Counting Kit-8 (CCK-8) and 5-ethynyl-2’-deoxyuridine (EdU) incorporation assays as previously described (39). Briefly, for CCK-8 assay, HCC cells (2,000 per well) were plated into 96-well plate. After culture for 4 days, cell proliferation was measured using CCK-8 reagents (Dojindo, Kumamoto, Japan) following the provided protocol. EdU incorporation assays were performed using the Cell-Light EdU Apollo567 In Vitro Kit (RiboBio, Guangzhou, China) following the provided protocol. Cell apoptosis was measured using the Caspase-3 Activity Assay Kit (Cell Signaling Technology, Danvers, MA, USA) following the provided protocol. Cell migration was measured using transwell migration assay as we previously described (40). Cell invasion was measured using transwell invasion assay as we previously described (40).
Mouse Xenograft Models
Five-week-old male BALB/c nude mice were purchased from Shanghai SLAC Laboratory Animal Co. and maintained in specific pathogen-free conditions. SNU-398 cells with ADORA2A-AS1 overexpression or control were subcutaneously injected into the back flank of nude mice (5 × 106 cells per site). Subcutaneous xenograft volumes were measured every week following the formula volume = 0.5 × length × width2. At the 28th day after injection, subcutaneous xenografts were resected, weighed, and subjected to immunohistochemistry (IHC) staining using primary antibodies against Proliferating Cell Nuclear Antigen (PCNA) (#13110, 1:4,000; Cell Signaling Technology), Ki67 (#9027, 1:400; Cell Signaling Technology), or cleaved caspase-3 (#9664, 1:1,000; Cell Signaling Technology). The subcutaneous xenografts were also subjected to terminal deoxynucleotidyl transferase (TdT)-mediated dUTP nick end labeling (TUNEL) assay using the TUNEL Cell Apoptosis Detection Kit (Beyotime, Shanghai, China). SNU-398 cells with ADORA2A-AS1 overexpression or control were also intrasplenically injected into nude mice (2 × 106 cells per mouse). At the 35th day after injection, the livers were resected and subjected to hematoxylin and eosin (H&E) staining to count the diameter and number of metastatic nodules. All animal experiments were approved by the Animal Ethics Committee of Affiliated Hospital of Youjiang Medical University for Nationalities.
RNA Pull-Down
ADORA2A-AS1 was in vitro transcribed from pSPT19-ADORA2A-AS1 using the MEGAscript® Kit (Invitrogen) and Sp6 RNA polymerase (Invitrogen). ADORA2A-AS1 transcript was labeled using the Pierce™ RNA 3′ End Desthiobiotinylation Kit (Thermo Scientific, Waltham, MA, USA). RNA pull-down assays were performed in SNU-398 cells using ADORA2A-AS1 end-labeled with desthiobiotin and the Pierce™ Magnetic RNA-Protein Pull-Down Kit (Thermo Scientific) following the provided protocol. The protein present in the pull-down material was detected by Western blot.
Western Blot
Proteins extracted from indicated HCC cells or RNA pull-down material were subjected to sodium dodecyl sulfate–polyacrylamide gel electrophoresis and transferred to polyvinylidene fluoride membrane. After blocking, the membranes were incubated with primary antibodies against HuR (Hu Antigen R) (#03-102, 1:500; Millipore, Burlington, MA, USA), FSCN1 (#54545, 1:1,000; Cell Signaling Technology), phosphor-AKT (66444-1-Ig, 1:2,000; Proteintech, Chicago, IL, USA), AKT (10176-2-AP, 1:1,000; Proteintech), MMP2 (#40994, 1:1,000; Cell Signaling Technology), BIRC7 (#5471, 1:1,000; Cell Signaling Technology), or GAPDH (ab8245, 1:5,000; Abcam) overnight at 4°C. The secondary antibodies used were IRDye 680RD Goat anti-Mouse IgG and IRDye 800CW Goat anti-Rabbit IgG (Li-Cor, Lincoln, NE, USA).
RNA Immunoprecipitation
RNA immunoprecipitation (RIP) assays were undertaken in indicated HCC cells using the Magna RIP RNA-Binding Protein Immunoprecipitation Kit (Millipore) and a primary antibody against HuR (#03-102, 5 µg per reaction; Millipore) following the provided instruction. The precipitated RNAs were measured by qRT-PCR.
Statistical Analysis
All statistical analyses were performed using the GraphPad Prism 6.0 Software. Log-rank test, Pearson chi-square test, Student’s t-test, one-way ANOVA followed by Dunnett’s multiple comparisons test, and Pearson correlation analysis were used as indicated in the figure and table legends. p < 0.05 was considered statistically significant.
Results
Low Expression of ADORA2A-AS1 Was Associated With Advanced Stages and Poor Prognosis in Hepatocellular Carcinoma
To explore clinical significances of ADORA2A-AS1 in HCC, we analyzed the RNA-seq data of HCC tissues from TCGA project. The results revealed that low ADORA2A-AS1 expression was associated with worse overall survival of HCC patients (Figure 1A). Furthermore, correlation analyses between ADORA2A-AS1 expression and clinicopathological characteristics revealed that low ADORA2A-AS1 expression was associated with advanced American Joint Committee on Cancer (AJCC) pathologic t stage and pathologic tumor stage (Table 2). To further confirm the clinical significances of ADORA2A-AS1 in HCC, we randomly collected 76 HCC samples, measured ADORA2A-AS1 expression in these 76 HCC tissues, and analyzed the correlation between ADORA2A-AS1 expression and clinicopathological characteristics. Kaplan–Meier survival analyses revealed that low ADORA2A-AS1 expression was also associated with worse overall survival in our HCC cohort (Figure 1B). In this HCC cohort, low ADORA2A-AS1 expression was found to be correlated with bad encapsulation, microvascular invasion, and advanced Barcelona Clinic Liver Cancer (BCLC) stage (Table 1). These findings suggested ADORA2A-AS1 as a potential prognostic biomarker for HCC. Furthermore, the expressions of ADORA2A-AS1 in immortalized human liver cell THLE-2 and human HCC cells SNU-398, SK-HEP-1, and Hep3B were measured, and the results revealed that ADORA2A-AS1 was lowly expressed in HCC cells compared with that in immortalized liver cell (Figure 1C).
Figure 1 The correlation between ADORA2A-AS1 expression level and outcomes in hepatocellular carcinoma (HCC). (A) The correlation between ADORA2A-AS1 expression level and overall survival in HCC analyzed using the RNA sequencing (RNA-seq) data from The Cancer Genome Atlas (TCGA) project. n = 361, p = 0.0016, hazard ratio (HR) = 0.5710 by log-rank test. Median ADORA2A-AS1 expression level was used as cutoff. (B) Kaplan–Meier survival analysis of the correlation between ADORA2A-AS1 expression level and overall survival in our HCC cohort. n = 76, p = 0.0357, HR = 0.5178 by log-rank test. Median ADORA2A-AS1 expression level was used as cutoff. (C) ADORA2A-AS1 expression levels in immortalized human liver cell THLE-2 and human HCC cells SNU-398, SK-HEP-1, and Hep3B were measured by qRT-PCR. Results are presented as mean ± SD based on three independent experiments. **p < 0.01, ***p < 0.001 by one-way ANOVA followed by Dunnett’s multiple comparisons test.
Table 2 Correlation between ADORA2A-AS1 expression and clinicopathological characteristics in HCC according to TCGA dataset.
Two in silico tools CPAT (http://lilab.research.bcm.edu/) and TestCode (https://www.bioinformatics.org/sms/testcode.html) both indicated that ADORA2A-AS1 was a noncoding RNA (Figures S1A, B). RNA FISH results showed that ADORA2A-AS1 was mainly distributed in the cytoplasm (Figure S1C). Subcellular fractionation followed by qRT-PCR results also showed that ADORA2A-AS1 was mainly distributed in the cytoplasm (Figure S1D).
ADORA2A-AS1 Repressed Hepatocellular Carcinoma Cellular Proliferation and Promoted Hepatocellular Carcinoma Cellular Apoptosis
Due to the correlation between low ADORA2A-AS1 expression and bad encapsulation, microvascular invasion, advanced pathologic t stage, advanced pathologic tumor stage, advanced BCLC stage, and poor overall survival, we next detected the potential roles of ADORA2A-AS1 in HCC. We stably overexpressed ADORA2A-AS1 in SNU-398 and SK-HEP-1 cells via stable transfecting of ADORA2A-AS1 overexpression vector (Figure 2A). We also stably silenced ADORA2A-AS1 expression in SNU-398 and Hep3B cells via stable infecting of ADORA2A-AS1 shRNA lentiviruses (Figure 2B). CCK-8 assay results showed that ADORA2A-AS1 overexpression repressed cell proliferation of SNU-398 and SK-HEP-1 cells (Figure 2C). Conversely, CCK-8 assay results also showed that ADORA2A-AS1 silencing promoted cell proliferation of SNU-398 and Hep3B cells (Figure 2D). EdU incorporation assay results showed that SNU-398 and SK-HEP-1 cells with ADORA2A-AS1 overexpression had less EdU-positive cells compared with their control cells (Figure 2E), indicating that ADORA2A-AS1 overexpression repressed cell proliferation of SNU-398 and SK-HEP-1 cells. Conversely, SNU-398 and Hep3B cells with ADORA2A-AS1 silencing had more EdU-positive cells compared with their control cells (Figure 2F), indicating that ADORA2A-AS1 silencing promoted cell proliferation of SNU-398 and Hep3B cells. Caspase-3 activity assay results showed that ADORA2A-AS1 overexpression increased caspase-3 activity in SNU-398 and SK-HEP-1 cells (Figure 2G), indicating that ADORA2A-AS1 overexpression promoted cell apoptosis of SNU-398 and SK-HEP-1 cells. Conversely, caspase-3 activity assay results also showed that ADORA2A-AS1 silencing decreased caspase-3 activity in SNU-398 and Hep3B cells (Figure 2H), indicating that ADORA2A-AS1 silencing repressed cell apoptosis of SNU-398 and Hep3B cells. These findings demonstrated that ADORA2A-AS1 repressed HCC cellular proliferation and promoted HCC cellular apoptosis.
Figure 2 The roles of ADORA2A-AS1 in hepatocellular carcinoma (HCC) cellular proliferation and apoptosis. (A) ADORA2A-AS1 expression levels in SNU-398 and SK-HEP-1 cells with ADORA2A-AS1 stable overexpression or control were measured by qRT-PCR. (B) ADORA2A-AS1 expression levels in SNU-398 and Hep3B cells with ADORA2A-AS1 stable silencing or control were measured by qRT-PCR. (C) Cell proliferation of SNU-398 and SK-HEP-1 cells with ADORA2A-AS1 overexpression or control was evaluated by Cell Counting Kit-8 (CCK-8) assays. (D) Cell proliferation of SNU-398 and Hep3B cells with ADORA2A-AS1 stable silencing or control was evaluated by CCK-8 assays. (E) Cell proliferation of SNU-398 and SK-HEP-1 cells with ADORA2A-AS1 overexpression or control was evaluated by 5-ethynyl-2’-deoxyuridine (EdU) incorporation assays. Red color indicated EdU incorporated nucleus. Scale bar, 200 µm. (F) Cell proliferation of SNU-398 and Hep3B cells with ADORA2A-AS1 stable silencing or control was evaluated by EdU incorporation assays. Red color indicated EdU incorporated nucleus. Scale bar, 200 µm. (G) Cell apoptosis of SNU-398 and SK-HEP-1 cells with ADORA2A-AS1 overexpression or control was evaluated by caspase-3 activity assays. (H) Cell apoptosis of SNU-398 and Hep3B cells with ADORA2A-AS1 stable silencing or control was evaluated by caspase-3 activity assays. Results are presented as mean ± SD based on three independent experiments. *p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001 by Student’s t-test (A, C, E, G) or one-way ANOVA followed by Dunnett’s multiple comparisons test (B, D, F, H).
ADORA2A-AS1 Repressed Hepatocellular Carcinoma Cellular Migration and Invasion
Given that low ADORA2A-AS1 expression was associated with bad encapsulation and microvascular invasion, we next detected the potential roles of ADORA2A-AS1 in HCC cellular migration and invasion. Transwell migration assay results showed that SNU-398 and SK-HEP-1 cells with ADORA2A-AS1 overexpression had less migrated cells compared with their control cells (Figure 3A). Conversely, SNU-398 and Hep3B cells with ADORA2A-AS1 silencing had more migrated cells compared with their control cells (Figure 3B). Transwell invasion assay results showed that SNU-398 and SK-HEP-1 cells with ADORA2A-AS1 overexpression had less invasive cells compared with their control cells (Figure 3C). Conversely, SNU-398 and Hep3B cells with ADORA2A-AS1 silencing had more invasive cells compared with their control cells (Figure 3D). These findings demonstrated that ADORA2A-AS1 repressed HCC cellular migration and invasion.
Figure 3 The roles of ADORA2A-AS1 in hepatocellular carcinoma (HCC) cellular migration and invasion. (A) Cell migration of SNU-398 and SK-HEP-1 cells with ADORA2A-AS1 overexpression or control was evaluated by transwell migration assays. Scale bar, 100 µm. (B) Cell migration of SNU-398 and Hep3B cells with ADORA2A-AS1 stable silencing or control was evaluated by transwell migration assays. Scale bar, 100 µm. (C) Cell invasion of SNU-398 and SK-HEP-1 cells with ADORA2A-AS1 overexpression or control was evaluated by transwell invasion assays. Scale bar, 100 µm. (D) Cell invasion of SNU-398 and Hep3B cells with ADORA2A-AS1 stable silencing or control was evaluated by transwell invasion assays. Scale bar, 100 µm. Results are presented as mean ± SD based on three independent experiments. *p < 0.05, **p < 0.01 by Student’s t-test (A, C) or one-way ANOVA followed by Dunnett’s multiple comparisons test (B, D).
ADORA2A-AS1 Repressed Hepatocellular Carcinoma Xenograft Growth and Metastasis In Vivo
To investigate the roles of ADORA2A-AS1 in mouse xenograft model in vivo, SNU-398 cells with ADORA2A-AS1 stable overexpression or control were subcutaneously injected into nude mice. Subcutaneous xenograft growth curve showed that ADORA2A-AS1 overexpression slowed xenograft growth rate (Figure 4A). At the 28th day after injection, subcutaneous xenografts were resected and weighed. Obviously, the xenografts formed by ADORA2A-AS1-overexpressed SNU-398 cells were lighter and smaller than the xenografts formed by control SNU-398 cells (Figures 4B, C). The xenografts were further subjected to PCNA and Ki67 IHC staining to evaluate cell proliferation in vivo. Proliferation markers PCNA and Ki67 IHC staining both showed that ADORA2A-AS1 overexpression repressed cell proliferation in vivo (Figures 4D, E). Furthermore, cleaved caspase-3 IHC staining of the xenografts showed that ADORA2A-AS1 overexpression promoted cell apoptosis in vivo (Figure 4F). TUNEL assay results further confirmed that ADORA2A-AS1 overexpression promoted cell apoptosis in vivo (Figure 4G). To investigate the roles of ADORA2A-AS1 in HCC metastasis, SNU-398 cells with ADORA2A-AS1 stable overexpression or control were intrasplenically injected into nude mice. At the 35th day after injection, liver metastases were evaluated by H&E staining. As shown in Figures 4H–J, ADORA2A-AS1-overexpressed SNU-398 cells formed significantly smaller and less liver metastatic nodules. These findings demonstrated that ADORA2A-AS1 repressed HCC xenograft growth and metastasis in vivo.
Figure 4 The roles of ADORA2A-AS1 in hepatocellular carcinoma (HCC) xenograft growth and metastasis in vivo. (A) Subcutaneous xenografts volume over time after injection of SNU-398 cells with ADORA2A-AS1 overexpression or control. (B) Subcutaneous xenograft weights at the 28th day after injection of SNU-398 cells with ADORA2A-AS1 overexpression or control. (C) The photo of subcutaneous xenografts resected at the 28th day after injection of SNU-398 cells with ADORA2A-AS1 overexpression or control. (D) PCNA immunohistochemistry (IHC) staining of subcutaneous xenografts formed by SNU-398 cells with ADORA2A-AS1 overexpression or control. Scale bar, 50 µm. (E) Ki67 IHC staining of subcutaneous xenografts formed by SNU-398 cells with ADORA2A-AS1 overexpression or control. Scale bar, 50 µm. (F) Cleaved caspase-3 IHC staining of subcutaneous xenografts formed by SNU-398 cells with ADORA2A-AS1 overexpression or control. Scale bar, 50 µm. (G) In vivo cell apoptosis in subcutaneous xenografts formed by SNU-398 cells with ADORA2A-AS1 overexpression or control was evaluated by terminal deoxynucleotidyl transferase (TdT)-mediated dUTP nick end labeling (TUNEL) assay. Scale bar, 50 µm. (H–J) SNU-398 cells with ADORA2A-AS1 overexpression or control were intrasplenically injected into nude mice. At the 35th day after injection, maximal diameter and number of liver metastatic nodules were evaluated by H&E staining. Scale bar, 1,000 µm. Results are presented as mean ± SD based on n = 5 mice in each group. *p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001 by Student’s t-test.
ADORA2A-AS1 Downregulated FSCN1 via Binding HuR
The most common mechanism of action of lncRNAs is to bind proteins and regulate the function of the interacted proteins (20). To investigate whether ADORA2A-AS1 also binds protein and which proteins are bound to ADORA2A-AS1, two in silico tools RBPmap (mapping binding sites of RNA binding proteins) (http://rbpmap.technion.ac.il/) and ATtRACT (a database of RNA binding proteins and associated motifs) (https://attract.cnic.es/#) were used to predict the proteins bound to ADORA2A-AS1. Intriguingly, RNA-binding protein HuR was predicted to bind ADORA2A-AS1 by both tools (Tables S1, S2). The binding of ADORA2A-AS1 to HuR was further predicted by another in silico tool RPISeq (RNA-Protein Interaction Prediction) (http://pridb.gdcb.iastate.edu/RPISeq/index.html) with a score of 0.8. RNA pull-down assays using ADORA2A-AS1 end-labeled with desthiobiotin showed that HuR was specifically enriched in ADORA2A-AS1 group (Figure 5A). RIP assays using anti-HuR antibody showed that ADORA2A-AS1 was specifically enriched in anti-HuR antibody group (Figure 5B). These results demonstrated that ADORA2A-AS1 specifically bound to HuR. HuR is reported to bind and stabilize several mRNAs, including FSCN1, VEGFC, MAT2A, MDM2, CCND1, CASP3, USP7, FAS, PDGFA, CTSS, GATA3, and EGFR (41–46). The RNA-seq data of HCC tissues from TCGA project showed that among all these HuR targets, FSCN1 has the most significant correlation with ADORA2A-AS1 with r value of -0.4627 (Figure 5C and Figure S2). The significantly negative correlation between ADORA2A-AS1 and FSCN1 expression levels in HCC tissues was also found in our HCC cohort (Figure 5D). Therefore, we further explore the potential effects of ADORA2A-AS1 on FSCN1. RIP assays using anti-HuR antibody in SNU-398 cells with ADORA2A-AS1 stable overexpression or control showed that HuR bound to FSCN1 transcript and ADORA2A-AS1 overexpression repressed the binding of HuR to FSCN1 transcript (Figure 5E). Conversely, RIP assays using anti-HuR antibody in SNU-398 cells with ADORA2A-AS1 stable silencing or control showed that ADORA2A-AS1 silencing promoted the binding of HuR to FSCN1 transcript (Figure 5F). These findings suggested that ADORA2A-AS1 repressed the binding of HuR to FSCN1 transcript via competitively binding HuR. Given that HuR could bind and stabilize FSCN1 transcript, we next investigated the effects of ADORA2A-AS1 on FSCN1 transcript stability. SNU-398 cells with ADORA2A-AS1 stable overexpression or silencing were treated with α-amanitin to block new RNA transcription. FSCN1 transcript stability over time was measured, and the results showed that ADORA2A-AS1 overexpression shortened the half-life of FSCN1 transcript (Figure 5G). Conversely, ADORA2A-AS1 silencing elongated the half-life of FSCN1 transcript (Figure 5H). As expected, FSCN1 transcript level was decreased in SNU-398 cells with ADORA2A-AS1 overexpression and increased in SNU-398 cells with ADORA2A-AS1 silencing (Figures 5I, J). FSCN1 protein level was also decreased in SNU-398 cells with ADORA2A-AS1 overexpression and increased in SNU-398 cells with ADORA2A-AS1 silencing (Figures 5K, L). These findings suggested that ADORA2A-AS1 bound to HuR, repressed the binding of HuR to FSCN1 transcript, reduced FSCN1 transcript stability, and downregulated FSCN1 expression.
Figure 5 ADORA2A-AS1 upregulated FSCN1 expression via binding HuR. (A) Proteins were enriched by RNA pull-down assays using ADORA2A-AS1 end-labeled with desthiobiotin, followed by Western blot to detect the binding of HuR to ADORA2A-AS1. (B) RIP assays using anti-HuR antibody, followed by qRT-PCR to detect the binding of ADORA2A-AS1 to HuR. (C) The correlation between FSCN1 and ADORA2A-AS1 expression in hepatocellular carcinoma (HCC) tissues analyzed using the RNA-seq data from The Cancer Genome Atlas (TCGA) project. n = 360, p < 0.0001, r = -0.4627 by Pearson correlation analysis. (D) The correlation between FSCN1 and ADORA2A-AS1 expression in our HCC cohort. n = 76, p < 0.0001, r = -0.4975 by Pearson correlation analysis. (E) RIP assays in SNU-398 cells with ADORA2A-AS1 overexpression or control using anti-HuR antibody, followed by qRT-PCR to detect the binding of FSCN1 transcript to HuR. (F) RIP assays in SNU-398 cells with ADORA2A-AS1 silencing or control using anti-HuR antibody, followed by qRT-PCR to detect the binding of FSCN1 transcript to HuR. (G) FSCN1 transcript stability in SNU-398 cells with ADORA2A-AS1 overexpression or control over time was measured after blocking new RNA transcription using 50 µM α-amanitin. Here, 18S rRNA was used as endogenous control, whose transcription was not repressed by α-amanitin. (H) FSCN1 transcript stability in SNU-398 cells with ADORA2A-AS1 silencing or control over time was measured after blocking new RNA transcription using 50 µM α-amanitin. Here, 18S rRNA was used as endogenous control. (I) FSCN1 transcript levels in SNU-398 cells with ADORA2A-AS1 overexpression or control were measured by qRT-PCR. (J) FSCN1 transcript levels in SNU-398 cells with ADORA2A-AS1 silencing or control were measured by qRT-PCR. (K) FSCN1 protein levels in SNU-398 cells with ADORA2A-AS1 overexpression or control were measured by Western blot. (L) FSCN1 protein levels in SNU-398 cells with ADORA2A-AS1 silencing or control were measured by Western blot. Results are presented as mean ± SD based on three independent experiments. *p < 0.05, **p < 0.01, ***p < 0.001, ns, not significant, by Student’s t-test (B, E, G, I) or one-way ANOVA followed by Dunnett’s multiple comparisons test (F, H, J).
ADORA2A-AS1 Repressed PI3K/AKT Pathway
FSCN1, also named as Fascin, is an actin-bundling protein (47). High FSCN1 expression has been revealed to be associated with worse outcomes in HCC (48). Increasing evidence showed that FSCN1 is involved in cell proliferation, apoptosis, motility, chemotherapeutic resistance, tumor growth, and metastasis of several malignancies (49, 50). PI3K/AKT pathway was revealed by several reports to be a critical downstream target of FSCN1 (51, 52). Thus, we further investigate the effects of ADORA2A-AS1 on PI3K/AKT pathway. AKT phosphorylation level was decreased in SNU-398 cells with ADORA2A-AS1 overexpression compared with control SNU-398 cells (Figure 6A). Conversely, AKT phosphorylation level was increased in SNU-398 cells with ADORA2A-AS1 silencing (Figure 6B). MMP2 is a downstream target of AKT activation, which is frequently reported to be involved in cell motility (53). BIRC7, known as Livin, is another critical downstream target of AKT activation, which exerts antiapoptotic roles (51). MMP2 and BIRC7 expression levels were detected in SNU-398 cells with ADORA2A-AS1 overexpression or silencing. As shown in Figures 6C, D, mRNA and protein levels of MMP2 and BIRC7 were both decreased in SNU-398 cells with ADORA2A-AS1 overexpression compared with control SNU-398 cells. Conversely, mRNA and protein levels of MMP2 and BIRC7 were both increased in SNU-398 cells with ADORA2A-AS1 silencing (Figures 6E, F). The RNA-seq data of HCC tissues from TCGA project showed that MMP2 expression level was negatively correlated with ADORA2A-AS1 expression level (Figure 6G). The negative correlation between MMP2 and ADORA2A-AS1 expression levels was also confirmed in our HCC cohort (Figure 6H). BIRC7 expression level was also negatively correlated with ADORA2A-AS1 expression level in HCC tissues, analyzed using the RNA-seq data from TCGA project (Figure 6I). Our HCC cohort further supported the negative correlation between BIRC7 and ADORA2A-AS1 expression levels (Figure 6J). These findings suggested that ADORA2A-AS1 repressed PI3K/AKT pathway via downregulating FSCN1.
Figure 6 ADORA2A-AS1 repressed AKT pathway. (A) AKT phosphorylation levels in SNU-398 cells with ADORA2A-AS1 overexpression or control were measured by Western blot. (B) AKT phosphorylation levels in SNU-398 cells with ADORA2A-AS1 silencing or control were measured by Western blot. (C) MMP2 and BIRC7 transcript levels in SNU-398 cells with ADORA2A-AS1 overexpression or control were measured by qRT-PCR. (D) MMP2 and BIRC7 protein levels in SNU-398 cells with ADORA2A-AS1 overexpression or control were measured by Western blot. (E) MMP2 and BIRC7 transcript levels in SNU-398 cells with ADORA2A-AS1 silencing or control were measured by qRT-PCR. (F) MMP2 and BIRC7 protein levels in SNU-398 cells with ADORA2A-AS1 silencing or control were measured by Western blot. (G) The correlation between MMP2 and ADORA2A-AS1 expression in hepatocellular carcinoma (HCC) tissues analyzed using the RNA-seq data from The Cancer Genome Atlas (TCGA) project. n = 360, p < 0.0001, r = -0.3166 by Pearson correlation analysis. (H) The correlation between MMP2 and ADORA2A-AS1 expression in our HCC cohort. n = 76, p < 0.0001, r = -0.4625 by Pearson correlation analysis. (I) The correlation between BIRC7 and ADORA2A-AS1 expression in HCC tissues analyzed using the RNA-seq data from TCGA project. n = 360, p < 0.0001, r = -0.2428 by Pearson correlation analysis. (J) The correlation between BIRC7 and ADORA2A-AS1 expression in our HCC cohort. n = 76, p = 0.0004, r = -0.3941 by Pearson correlation analysis. For (C, E), results are presented as mean ± SD based on three independent experiments. **p < 0.01, ***p < 0.001 by Student’s t-test (C) or one-way ANOVA followed by Dunnett’s multiple comparisons test (E).
FSCN1 Overexpression Attenuated the Tumor-Suppressive Roles of ADORA2A-AS1 in Hepatocellular Carcinoma
To explore whether FSCN1 is the critical mediator of the roles of ADORA2A-AS1 in HCC, we stably overexpressed FSCN1 in SNU-398 cells with ADORA2A-AS1 stable overexpression (Figure 7A). CCK-8 assay results showed that FSCN1 overexpression reversed the decreased cell proliferation caused by ADORA2A-AS1 overexpression (Figure 7B). EdU incorporation assay results showed that FSCN1 overexpression reversed the decreased EdU-positive cell percentage caused by ADORA2A-AS1 overexpression (Figure 7C), indicating that FSCN1 overexpression rescued cell proliferation that was repressed by ADORA2A-AS1. Caspase-3 activity assay results showed that FSCN1 overexpression reversed the increased caspase-3 activity caused by ADORA2A-AS1 overexpression (Figure 7D), indicating that FSCN1 overexpression reversed cell apoptosis that was induced by ADORA2A-AS1. Transwell migration assay results showed that FSCN1 overexpression reversed the decreased migrated cell number caused by ADORA2A-AS1 overexpression (Figure 7E). Transwell invasion assay results showed that FSCN1 overexpression reversed the decreased invasive cell number caused by ADORA2A-AS1 overexpression (Figure 7F). These findings suggested that FSCN1 reversed the roles of ADORA2A-AS1 in repressing cell proliferation, inducing cell apoptosis, and repressing cell migration and invasion.
Figure 7 Functional rescue assays using FSCN1 overexpression. (A) FSCN1 protein levels in SNU-398 cells with ADORA2A-AS1 and FSCN1 overexpression were measured by Western blot. (B) Cell proliferation of SNU-398 cells with ADORA2A-AS1 and FSCN1 overexpression was evaluated by CCK-8 assays. (C) Cell proliferation of SNU-398 cells with ADORA2A-AS1 and FSCN1 overexpression was evaluated by 5-ethynyl-2’-deoxyuridine (EdU) incorporation assays. Red color indicated EdU-incorporated nucleus. Scale bar, 200 µm. (D) Cell apoptosis of SNU-398 cells with ADORA2A-AS1 and FSCN1 overexpression was evaluated by caspase-3 activity assays. (E) Cell migration of SNU-398 cells with ADORA2A-AS1 and FSCN1 overexpression was evaluated by transwell migration assays. Scale bar, 100 µm. (F) Cell invasion of SNU-398 cells with ADORA2A-AS1 and FSCN1 overexpression was evaluated by transwell invasion assays. Scale bar, 100 µm. Results are presented as mean ± SD based on three independent experiments. **p < 0.01, ns, not significant, by one-way ANOVA followed by Dunnett’s multiple comparisons test.
Inhibition of PI3K/AKT Pathway Attenuated the Oncogenic Roles of ADORA2A-AS1 Depletion in Hepatocellular Carcinoma
To further explore whether PI3K/AKT pathway is also the critical mediator of the roles of ADORA2A-AS1 in HCC, SNU-398 cells with ADORA2A-AS1 stable silencing were treated with PI3K/AKT pathway inhibitor LY294002. CCK-8 assay results showed that LY294002 reversed the increased cell proliferation caused by ADORA2A-AS1 silencing (Figure 8A). EdU incorporation assay results showed that LY294002 reversed the increased EdU-positive cell percentage caused by ADORA2A-AS1 silencing (Figure 8B), indicating that LY294002 reversed the increased cell proliferation induced by ADORA2A-AS1 silencing. Caspase-3 activity assay results showed that LY294002 reversed the reduced caspase-3 activity caused by ADORA2A-AS1 silencing (Figure 8C), indicating that LY294002 reversed the reduced cell apoptosis caused by ADORA2A-AS1 silencing. Transwell migration assay results showed that LY294002 reversed the increased migrated cell number caused by ADORA2A-AS1 silencing (Figure 8D). Transwell invasion assay results showed that LY294002 reversed the increased invasive cell number caused by ADORA2A-AS1 silencing (Figure 8E). These findings suggested that the oncogenic roles of ADORA2A-AS1 depletion in HCC were abolished by PI3K/AKT pathway inhibition.
Figure 8 PI3K/AKT pathway inhibitor LY294002 abolished the roles of ADORA2A-AS1 depletion in hepatocellular carcinoma (HCC). (A) Cell proliferation of ADORA2A-AS1-silenced SNU-398 cells treated with 10 µM LY294002 was evaluated by Cell Counting Kit-8 (CCK-8) assays. (B) Cell proliferation of ADORA2A-AS1-silenced SNU-398 cells treated with 10 µM LY294002 was evaluated by 5-ethynyl-2’-deoxyuridine (EdU) incorporation assays. Red color indicated EdU incorporated nucleus. Scale bar, 200 µm. (C) Cell apoptosis of ADORA2A-AS1-silenced SNU-398 cells treated with 10 µM LY294002 was evaluated by caspase-3 activity assays. (D) Cell migration of ADORA2A-AS1-silenced SNU-398 cells treated with 10 µM LY294002 was evaluated by transwell migration assays. Scale bar, 100 µm. (E) Cell invasion of ADORA2A-AS1-silenced SNU-398 cells treated with 10 µM LY294002 was evaluated by transwell invasion assays. Scale bar, 100 µm. Results are presented as mean ± SD based on three independent experiments. **p < 0.01, ns, not significant, by one-way ANOVA followed by Dunnett’s multiple comparisons test.
Discussion
In this study, we identified an HCC-related tumor-suppressive lncRNA ADORA2A-AS1. ADORA2A-AS1 locates at chromosome 22q11.23. ADORA2A-AS1 has limited protein-encoding ability. Low expression of ADORA2A-AS1 in HCC tissues was correlated with advanced stage, bad encapsulation, microvascular invasion, and poor survival of HCC patients. Our findings suggested ADORA2A-AS1 as a potential prognostic biomarker for HCC, which needs further multicentral investigation.
Gain and loss-of functional experiments demonstrated that ADORA2A-AS1 repressed HCC cellular proliferation, induced HCC cellular apoptosis, repressed HCC cellular migration and invasion, and restricted HCC xenograft growth and metastasis in vivo. Thus, our findings suggested ADORA2A-AS1 as a tumor-suppressive lncRNA in HCC. Enhancing ADORA2A-AS1 expression represented potential therapeutic strategy against HCC.
Mechanistic investigations identified HuR as a critical interaction partner of ADORA2A-AS1. HuR, also named ELAVL1 (embryonic lethal, abnormal vision Drosophila-like protein 1), is an RNA-binding protein. HuR posttranscriptionally regulates expression of target genes via directly binding target mRNAs and further stabilizing target mRNAs (54). Previous reports have shown that high expression of HuR was associated with poor survival in HCC (45). HuR promotes HCC malignant progression via regulating several mRNA targets (45).
Combining bioinformatics prediction, RNA pull-down, and RIP assays, we documented the specific physical binding between ADORA2A-AS1 and HuR. The binding of ADORA2A-AS1 to HuR competitively suppressed the binding of HuR to its mRNA targets. Here, we identified FSCN1 as one of the most significantly changed mRNA targets regulated by ADORA2A-AS1. ADORA2A-AS1 competitively suppressed the binding of HuR to FSCN1 transcript, leading to the reduction of FSCN1 transcript stability. Thus, ADORA2A-AS1 downregulated FSCN1 expression. The expression of FSCN1 was significantly negatively correlated with ADORA2A-AS1 in HCC tissues, which was verified in both TCGA project and our own HCC cohort. According to the RNA-seq data of TCGA project, the expression abundance of ADORA2A-AS1 was similar with that of FSCN1 transcript in HCC tissues, supporting the competition between ADORA2A-AS1 and FSCN1 transcript to bind HuR. Potential effects of ADORA2A-AS1 on other HuR targets need further investigation.
Conversely to ADORA2A-AS1, high expression of FSCN1 was correlated with poor outcome of HCC patients (48). FSCN1 was documented to exert oncogenic roles in several malignancies, including HCC (49, 50). For example, FSCN1 was reported to promote HCC cellular proliferation, migration, and invasion, inhibit HCC cellular apoptosis, and increase HCC cellular doxorubicin resistance (55–58). PI3K/AKT pathway is the classically oncogenic pathway, which has been documented to modulate cell proliferation, survival, and motility in various malignancies (59–62). PI3K/AKT pathway was also revealed to be the critical downstream target of FSCN1 (51). Consistently, in this study, we found that ADORA2A-AS1 repressed PI3K/AKT pathway through downregulating FSCN1. MMP2 and BIRC7 are two representative targets of FSCN1/AKT (51, 53). In this study, we further found that ADORA2A-AS1 downregulated MMP2 and BIRC7 via repressing FSCN1/AKT axis. The expressions of MMP2 and BIRC7 were both negatively correlated with ADORA2A-AS1 in HCC tissues, verified in both TCGA project and our own HCC cohort. Functional rescue experiments showed that FSCN1/AKT axis was the critical mediator of the roles of ADORA2A-AS1 in HCC. We acknowledge that this study also has some limitations, such as the potential effects of ADORA2A-AS1 on other HuR targets and FSCN1 targets. Nevertheless, FSCN1/AKT axis is one of the critical mediators of the functions of ADORA2A-AS1 in HCC.
In conclusion, ADORA2A-AS1 is identified as a novel HCC-related tumor-suppressive lncRNA, whose low expression is correlated with poor outcome of HCC patients. ADORA2A-AS1 exerts tumor-suppressive roles in HCC via competitively binding HuR, decreasing the binding of HuR to FSCN1 transcript, downregulating FSCN1 transcript stability, repressing FSCN1 expression, and repressing AKT pathway. ADORA2A-AS1/HuR/FSCN1/AKT represent potential therapeutic targets for HCC.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Ethics Statement
The studies involving human participants were reviewed and approved by the Ethics Committee of Affiliated Hospital of Youjiang Medical University for Nationalities. The patients/participants provided their written informed consent to participate in this study. The animal study was reviewed and approved by the Animal Ethics Committee of Affiliated Hospital of Youjiang Medical University for Nationalities.
Author Contributions
HW and JP designed the study. JP, YZ, AW, ZQ, CZ, WL, and ZX undertook the experiments. JP, HW, YZ, QT, and JW analyzed the data. JP and HW are the major contributors in writing the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the Guangxi Natural Science Foundation Project (2019GXNSFBA245023), Guangxi Science and Technology Project (AD17129025), and Baise Scientific Research and Technology Development Programme (BKZ2020-47-23).
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.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2021.754835/full#supplementary-material
References
1. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J Clin (2021) 71:209–49. doi: 10.3322/caac.21660
2. Kudo Y, Sugimoto M, Arias E, Kasashima H, Cordes T, Linares JF, et al. PKClambda/iota Loss Induces Autophagy, Oxidative Phosphorylation, and NRF2 to Promote Liver Cancer Progression. Cancer Cell (2020) 38:247–262 e11. doi: 10.1016/j.ccell.2020.05.018
3. Fountzilas C, Gupta M, Lee S, Krishnamurthi S, Estfan B, Wang K, et al. A Multicentre Phase 1b/2 Study of Tivozanib in Patients With Advanced Inoperable Hepatocellular Carcinoma. Br J Cancer (2020) 122:963–70. doi: 10.1038/s41416-020-0737-6
4. Song M, Pan Q, Yang J, He J, Zeng J, Cheng S, et al. Galectin-3 Favours Tumour Metastasis via the Activation of Beta-Catenin Signalling in Hepatocellular Carcinoma. Br J Cancer (2020) 123:1521–34. doi: 10.1038/s41416-020-1022-4
5. Chew SC, Choo SY, Chow PK. A New Perspective on the Immune Escape Mechanism in HCC: Onco-Foetal Reprogramming. Br J Cancer (2021) 124:1897–9. doi: 10.1038/s41416-021-01286-0
6. Johnson PJ, Dhanaraj S, Berhane S, Bonnett L, Ma YT. The Prognostic and Diagnostic Significance of the Neutrophil-to-Lymphocyte Ratio in Hepatocellular Carcinoma: A Prospective Controlled Study. Br J Cancer (2021) 125:714–6. doi: 10.1038/s41416-021-01445-3
7. Dolgormaa G, Harimoto N, Ishii N, Yamanaka T, Hagiwara K, Tsukagoshi M, et al. Mac-2-Binding Protein Glycan Isomer Enhances the Aggressiveness of Hepatocellular Carcinoma by Activating mTOR Signaling. Br J Cancer (2020) 123:1145–53. doi: 10.1038/s41416-020-0971-y
8. Iyer MK, Niknafs YS, Malik R, Singhal U, Sahu A, Hosono Y, et al. The Landscape of Long Noncoding RNAs in the Human Transcriptome. Nat Genet (2015) 47:199–208. doi: 10.1038/ng.3192
9. Pamedytyte D, Leipute E, Zilaitiene B, Sarauskas V, Dauksiene D, Dauksa A, et al. Different Stability of miRNAs and Endogenous Control Genes in Archival Specimens of Papillary Thyroid Carcinoma. Mol Med (2020) 26:100. doi: 10.1186/s10020-020-00218-7
10. Volovat SR, Volovat C, Hordila I, Hordila DA, Mirestean CC, Miron OT, et al. MiRNA and LncRNA as Potential Biomarkers in Triple-Negative Breast Cancer: A Review. Front Oncol (2020) 10:526850. doi: 10.3389/fonc.2020.526850
11. Luo H, Zhu G, Xu J, Lai Q, Yan B, Guo Y, et al. HOTTIP lncRNA Promotes Hematopoietic Stem Cell Self-Renewal Leading to AML-Like Disease in Mice. Cancer Cell (2019) 36:645–59.e8. doi: 10.1016/j.ccell.2019.10.011
12. Schmitt AM, Chang HY. Long Noncoding RNAs in Cancer Pathways. Cancer Cell (2016) 29:452–63. doi: 10.1016/j.ccell.2016.03.010
13. Acha-Sagredo A, Uko B, Pantazi P, Bediaga NG, Moschandrea C, Rainbow L, et al. Long Non-Coding RNA Dysregulation is a Frequent Event in Non-Small Cell Lung Carcinoma Pathogenesis. Br J Cancer (2020) 122:1050–8. doi: 10.1038/s41416-020-0742-9
14. Liu HT, Ma RR, Lv BB, Zhang H, Shi DB, Guo XY, et al. LncRNA-HNF1A-AS1 Functions as a Competing Endogenous RNA to Activate PI3K/AKT Signalling Pathway by Sponging miR-30b-3p in Gastric Cancer. Br J Cancer (2020) 122:1825–36. doi: 10.1038/s41416-020-0836-4
15. Mo S, Zhang L, Dai W, Han L, Wang R, Xiang W, et al. Antisense lncRNA LDLRAD4-AS1 Promotes Metastasis by Decreasing the Expression of LDLRAD4 and Predicts a Poor Prognosis in Colorectal Cancer. Cell Death Dis (2020) 11:155. doi: 10.1038/s41419-020-2338-y
16. Bispo S, Farias TDJ, de Araujo-Souza PS, Cintra R, Dos Santos HG, Jorge NAN, et al. Dysregulation of Transcription Factor Networks Unveils Different Pathways in Human Papillomavirus 16-Positive Squamous Cell Carcinoma and Adenocarcinoma of the Uterine Cervix. Front Oncol (2021) 11:626187. doi: 10.3389/fonc.2021.626187
17. Daneshvar K, Ardehali MB, Klein IA, Hsieh FK, Kratkiewicz AJ, Mahpour A, et al. lncRNA DIGIT and BRD3 Protein Form Phase-Separated Condensates to Regulate Endoderm Differentiation. Nat Cell Biol (2020) 22:1211–22. doi: 10.1038/s41556-020-0572-2
18. Li G, Kryczek I, Nam J, Li X, Li S, Li J, et al. LIMIT is an Immunogenic lncRNA in Cancer Immunity and Immunotherapy. Nat Cell Biol (2021) 23:526–37. doi: 10.1038/s41556-021-00672-3
19. Li Y, Yan G, Zhang J, Chen W, Ding T, Yin Y, et al. LncRNA HOXA11-AS Regulates Calcium Oxalate Crystal-Induced Renal Inflammation via miR-124-3p/MCP-1. J Cell Mol Med (2020) 24:238–49. doi: 10.1111/jcmm.14706
20. Liu SJ, Dang HX, Lim DA, Feng FY, Maher CA. Long Noncoding RNAs in Cancer Metastasis. Nat Rev Cancer (2021) 21:446–60. doi: 10.1038/s41568-021-00353-1
21. Zhang Y, Li Y, Hu Q, Xi Y, Xing Z, Zhang Z, et al. The lncRNA H19 Alleviates Muscular Dystrophy by Stabilizing Dystrophin. Nat Cell Biol (2020) 22:1332–45. doi: 10.1038/s41556-020-00595-5
22. Brex D, Barbagallo C, Mirabella F, Caponnetto A, Battaglia R, Barbagallo D, et al. LINC00483 Has a Potential Tumor-Suppressor Role in Colorectal Cancer Through Multiple Molecular Axes. Front Oncol (2020) 10:614455. doi: 10.3389/fonc.2020.614455
23. Jafarzadeh M, Soltani BM. Long Noncoding RNA LOC400043 (LINC02381) Inhibits Gastric Cancer Progression Through Regulating Wnt Signaling Pathway. Front Oncol (2020) 10:562253. doi: 10.3389/fonc.2020.562253
24. Li MM, Liu XH, Zhao YC, Ma XY, Zhou YC, Zhao YX, et al. Long Noncoding RNA KCNQ1OT1 Promotes Apoptosis in Neuroblastoma Cells by Regulating miR-296-5p/Bax Axis. FEBS J (2020) 287:561–77. doi: 10.1111/febs.15047
25. Sur S, Ray RB. Emerging Role of lncRNA ELDR in Development and Cancer. FEBS J (2021). doi: 10.1111/febs.15876
26. Zhu XT, Yuan JH, Zhu TT, Li YY, Cheng XY. Long Noncoding RNA Glypican 3 (GPC3) Antisense Transcript 1 Promotes Hepatocellular Carcinoma Progression via Epigenetically Activating GPC3. FEBS J (2016) 283:3739–54. doi: 10.1111/febs.13839
27. Yan X, Hu Z, Feng Y, Hu X, Yuan J, Zhao SD, et al. Comprehensive Genomic Characterization of Long Non-Coding RNAs Across Human Cancers. Cancer Cell (2015) 28:529–40. doi: 10.1016/j.ccell.2015.09.006
28. Farooqi AA, Nayyab S, Martinelli C, Berardi R, Katifelis H, Gazouli M, et al. Regulation of Hippo, TGFbeta/SMAD, Wnt/beta-Catenin, JAK/STAT, and NOTCH by Long Non-Coding RNAs in Pancreatic Cancer. Front Oncol (2021) 11:657965. doi: 10.3389/fonc.2021.657965
29. Li H, Yan R, Chen W, Ding X, Liu J, Chen G, et al. Long non Coding RNA SLC26A4-AS1 Exerts Antiangiogenic Effects in Human Glioma by Upregulating NPTX1 via NFKB1 Transcriptional Factor. FEBS J (2021) 288:212–28. doi: 10.1111/febs.15325
30. Li JK, Chen C, Liu JY, Shi JZ, Liu SP, Liu B, et al. Long Noncoding RNA MRCCAT1 Promotes Metastasis of Clear Cell Renal Cell Carcinoma via Inhibiting NPR3 and Activating P38-MAPK Signaling. Mol Cancer (2017) 16:111. doi: 10.1186/s12943-017-0681-0
31. Yuan JH, Liu XN, Wang TT, Pan W, Tao QF, Zhou WP, et al. The MBNL3 Splicing Factor Promotes Hepatocellular Carcinoma by Increasing PXN Expression Through the Alternative Splicing of lncRNA-PXN-As1. Nat Cell Biol (2017) 19:820–32. doi: 10.1038/ncb3538
32. Shen W, Huang B, He Y, Shi L, Yang J. Long Non-Coding RNA RP11-820 Promotes Extracellular Matrix Production via Regulating miR-3178/MYOD1 in Human Trabecular Meshwork Cells. FEBS J (2020) 287:978–90. doi: 10.1111/febs.15058
33. Yuan JH, Yang F, Wang F, Ma JZ, Guo YJ, Tao QF, et al. A Long Noncoding RNA Activated by TGF-Beta Promotes the Invasion-Metastasis Cascade in Hepatocellular Carcinoma. Cancer Cell (2014) 25:666–81. doi: 10.1016/j.ccr.2014.03.010
34. Fa J, Zhang X, Zhang X, Qi M, Zhang X, Fu Q, et al. Long Noncoding RNA lnc-TSSK2-8 Activates Canonical Wnt/beta-Catenin Signaling Through Small Heat Shock Proteins HSPA6 and CRYAB. Front Cell Dev Biol (2021) 9:660576. doi: 10.3389/fcell.2021.660576
35. Lin A, Li C, Xing Z, Hu Q, Liang K, Han L, et al. The LINK-A lncRNA Activates Normoxic HIF1alpha Signalling in Triple-Negative Breast Cancer. Nat Cell Biol (2016) 18:213–24. doi: 10.1038/ncb3295
36. Wang L, Park HJ, Dasari S, Wang S, Kocher JP, Li W. CPAT: Coding-Potential Assessment Tool Using an Alignment-Free Logistic Regression Model. Nucleic Acids Res (2013) 41:e74. doi: 10.1093/nar/gkt006
37. Muppirala UK, Honavar VG, Dobbs D. Predicting RNA-Protein Interactions Using Only Sequence Information. BMC Bioinformatics (2011) 12:489. doi: 10.1186/1471-2105-12-489
38. Paz I, Kosti I, Ares M Jr., Cline M, Mandel-Gutfreund Y. RBPmap: A Web Server for Mapping Binding Sites of RNA-Binding Proteins. Nucleic Acids Res (2014) 42:W361–7. doi: 10.1093/nar/gku406
39. Li J, Li MH, Wang TT, Liu XN, Zhu XT, Dai YZ, et al. SLC38A4 Functions as a Tumour Suppressor in Hepatocellular Carcinoma Through Modulating Wnt/beta-Catenin/MYC/HMGCS2 Axis. Br J Cancer (2021) 125:865–76. doi: 10.1038/s41416-021-01490-y
40. Pu J, Wang J, Wei H, Lu T, Wu X, Wu Y, et al. lncRNA MAGI2-AS3 Prevents the Development of HCC via Recruiting KDM1A and Promoting H3K4me2 Demethylation of the RACGAP1 Promoter. Mol Ther Nucleic Acids (2019) 18:351–62. doi: 10.1016/j.omtn.2019.08.020
41. Dufies M, Giuliano S, Ambrosetti D, Claren A, Ndiaye PD, Mastri M, et al. Sunitinib Stimulates Expression of VEGFC by Tumor Cells and Promotes Lymphangiogenesis in Clear Cell Renal Cell Carcinomas. Cancer Res (2017) 77:1212–26. doi: 10.1158/0008-5472.CAN-16-3088
42. Ghosh M, Aguila HL, Michaud J, Ai Y, Wu MT, Hemmes A, et al. Essential Role of the RNA-Binding Protein HuR in Progenitor Cell Survival in Mice. J Clin Invest (2009) 119:3530–43. doi: 10.1172/JCI38263
43. Ghosh U, Adhya S. Posttranscriptional Regulation of Cyclin D1 by ARE-Binding Proteins AUF1 and HuR in Cycling Myoblasts. J Biosci (2018) 43:685–91. doi: 10.1007/s12038-018-9788-8
44. Lin C, Zhang S, Wang Y, Wang Y, Nice E, Guo C, et al. Functional Role of a Novel Long Noncoding RNA TTN-AS1 in Esophageal Squamous Cell Carcinoma Progression and Metastasis. Clin Cancer Res (2018) 24:486–98. doi: 10.1158/1078-0432.CCR-17-1851
45. Papatheofani V, Levidou G, Sarantis P, Koustas E, Karamouzis MV, Pergaris A, et al. HuR Protein in Hepatocellular Carcinoma: Implications in Development, Prognosis and Treatment. Biomedicines (2021) 9:119. doi: 10.3390/biomedicines9020119
46. Vazquez-Chantada M, Fernandez-Ramos D, Embade N, Martinez-Lopez N, Varela-Rey M, Woodhoo A, et al. HuR/methyl-HuR and AUF1 Regulate the MAT Expressed During Liver Proliferation, Differentiation, and Carcinogenesis. Gastroenterology (2010) 138:1943–53. doi: 10.1053/j.gastro.2010.01.032
47. Lamb MC, Tootle TL. Fascin in Cell Migration: More Than an Actin Bundling Protein. Biology (Basel) (2020) 9:403. doi: 10.3390/biology9110403
48. Iguchi T, Aishima S, Umeda K, Sanefuji K, Fujita N, Sugimachi K, et al. Fascin Expression in Progression and Prognosis of Hepatocellular Carcinoma. J Surg Oncol (2009) 100:575–9. doi: 10.1002/jso.21377
49. Li S, Huang XT, Wang MY, Chen DP, Li MY, Zhu YY, et al. FSCN1 Promotes Radiation Resistance in Patients With PIK3CA Gene Alteration. Front Oncol (2021) 11:653005. doi: 10.3389/fonc.2021.653005
50. Lin S, Taylor MD, Singh PK, Yang S. How Does Fascin Promote Cancer Metastasis? FEBS J (2021) 288:1434–46. doi: 10.1111/febs.15484
51. Ghebeh H, Al-Khaldi S, Olabi S, Al-Dhfyan A, Al-Mohanna F, Barnawi R, et al. Fascin is Involved in the Chemotherapeutic Resistance of Breast Cancer Cells Predominantly via the PI3K/Akt Pathway. Br J Cancer (2014) 111:1552–61. doi: 10.1038/bjc.2014.453
52. Uyama N, Iimuro Y, Kawada N, Reynaert H, Suzumura K, Hirano T, et al. Fascin, a Novel Marker of Human Hepatic Stellate Cells, may Regulate Their Proliferation, Migration, and Collagen Gene Expression Through the FAK-PI3K-Akt Pathway. Lab Invest (2012) 92:57–71. doi: 10.1038/labinvest.2011.150
53. Guo B, Wu S, Zhu X, Zhang L, Deng J, Li F, et al. Micropeptide CIP2A-BP Encoded by LINC00665 Inhibits Triple-Negative Breast Cancer Progression. EMBO J (2020) 39:e102190. doi: 10.15252/embj.2019102190
54. Song X, Shi X, Li W, Zhang F, Cai Z. The RNA-Binding Protein HuR in Digestive System Tumors. BioMed Res Int (2020) 2020:9656051. doi: 10.1155/2020/9656051
55. Hayashi Y, Osanai M, Lee GH. Fascin-1 Expression Correlates With Repression of E-Cadherin Expression in Hepatocellular Carcinoma Cells and Augments Their Invasiveness in Combination With Matrix Metalloproteinases. Cancer Sci (2011) 102:1228–35. doi: 10.1111/j.1349-7006.2011.01910.x
56. Liu Y, Hong W, Zhou C, Jiang Z, Wang G, Wei G, et al. miR-539 Inhibits FSCN1 Expression and Suppresses Hepatocellular Carcinoma Migration and Invasion. Oncol Rep (2017) 37:2593–602. doi: 10.3892/or.2017.5549
57. Wang G, Zhu S, Gu Y, Chen Q, Liu X, Fu H. MicroRNA-145 and MicroRNA-133a Inhibited Proliferation, Migration, and Invasion, While Promoted Apoptosis in Hepatocellular Carcinoma Cells Via Targeting Fscn1. Dig Dis Sci (2015) 60:3044–52. doi: 10.1007/s10620-015-3706-9
58. Zhang Y, Lu Y, Zhang C, Huang D, Wu W, Zhang Y, et al. FSCN1 Increases Doxorubicin Resistance in Hepatocellular Carcinoma Through Promotion of Epithelial-Mesenchymal Transition. Int J Oncol (2018) 52:1455–64. doi: 10.3892/ijo.2018.4327
59. Marin-Bejar O, Rogiers A, Dewaele M, Femel J, Karras P, Pozniak J, et al. Evolutionary Predictability of Genetic Versus Nongenetic Resistance to Anticancer Drugs in Melanoma. Cancer Cell (2021) 39:1135–49. doi: 10.1016/j.ccell.2021.05.015
60. Pedroza DA, Rajamanickam V, Subramani R, Bencomo A, Galvez A, Lakshmanaswamy R. Progesterone Receptor Membrane Component 1 Promotes the Growth of Breast Cancers by Altering the Phosphoproteome and Augmenting EGFR/PI3K/AKT Signalling. Br J Cancer (2020) 123:1326–35. doi: 10.1038/s41416-020-0992-6
61. Perry JM, Tao F, Roy A, Lin T, He XC, Chen S, et al. Overcoming Wnt-Beta-Catenin Dependent Anticancer Therapy Resistance in Leukaemia Stem Cells. Nat Cell Biol (2020) 22:689–700. doi: 10.1038/s41556-020-0507-y
Keywords: hepatocellular carcinoma, long noncoding RNA, progression, HuR, Fascin, AKT pathway
Citation: Pu J, Zhang Y, Wang A, Qin Z, Zhuo C, Li W, Xu Z, Tang Q, Wang J and Wei H (2021) ADORA2A-AS1 Restricts Hepatocellular Carcinoma Progression via Binding HuR and Repressing FSCN1/AKT Axis. Front. Oncol. 11:754835. doi: 10.3389/fonc.2021.754835
Received: 07 August 2021; Accepted: 20 September 2021;
Published: 18 October 2021.
Edited by:
Rongxin Zhang, Guangdong Pharmaceutical University, ChinaReviewed by:
Cailin Xue, Nanjing Drum Tower Hospital, ChinaWafaa M. Rashed, Children’s Cancer Hospital, Egypt
Copyright © 2021 Pu, Zhang, Wang, Qin, Zhuo, Li, Xu, Tang, Wang and Wei. 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.
*Correspondence: Huamei Wei, weihuamei@yeah.net
†These authors have contributed equally to this work