Skip to main content

Liver-specific lncRNA FAM99A may be a tumor suppressor and promising prognostic biomarker in hepatocellular carcinoma

Abstract

Background

Increasing evidence shows that liver-specific long non-coding RNAs (lncRNAs) play important roles in the development of hepatocellular carcinoma (HCC). We identified a novel liver-specific lncRNA, FAM99A, and examined its clinical significance and biological functions in HCC.

Methods

The expression level and clinical value of FAM99A in HCC were examined using The Cancer Genome Atlas (TCGA), International Cancer Genome Consortium (ICGC), and Gene Expression Omnibus (GEO) databases, and were further verified using quantitative real-time polymerase chain reaction (qRT–PCR) in our HCC cohort. Univariate and multivariate Cox proportional hazards regression models were also applied to identify independent prognostic indicators for HCC patients. Cell counting kit-8, colony formation, and Transwell assays were performed to evaluate the effects of FAM99A on the proliferation, migration, and invasion abilities of HCC cells in vitro. A subcutaneous xenograft tumor model was implemented to determine the effect of FAM99A on the tumor growth of HCC cells in vivo. RNA pull-down and mass spectrometry assays were performed to reveal the potential molecular mechanisms of FAM99A in HCC.

Results

The three public online databases and qRT–PCR data showed that FAM99A was frequently downregulated in HCC tissues and inversely correlated with microvascular invasion and advanced histological grade of HCC patients. Kaplan–Meier survival analysis indicated that decreased FAM99A was significantly associated with poor overall survival of HCC patients based on TCGA database (P = 0.040), ICGC data portal (P < 0.001), and our HCC cohort (P = 0.010). A multivariate Cox proportional hazards regression model based on our HCC cohort suggested that FAM99A was an independent prognostic factor of overall survival for HCC patients (hazard ratio: 0.425, P = 0.039). Upregulation of FAM99A suppressed the proliferation, colony formation, migration, and invasion capacities of HCC cells in vitro, and knockdown of FAM99A had the opposite effects. A subcutaneous xenograft tumor model demonstrated that overexpression of FAM99A significantly inhibited the tumor growth of HCC cells in vivo. Seven tumor-related proteins (PCBP1, SRSF5, SRSF6, YBX1, IGF2BP2, HNRNPK, and HNRNPL) were recognized as possible FAM99A-binding proteins by the RNA pull-down and mass spectrometry assays.

Conclusion

Our results suggest that FAM99A exerts cancer-inhibiting effects on HCC progression, and it may be a promising prognostic indicator for HCC patients.

Peer Review reports

Introduction

Liver cancer is a common solid malignancy. It was the sixth most common tumor and the third leading cause of cancer-related mortality worldwide according to the World Health Organization’s statistics in 2020 [1]. China has a high incidence of liver cancer, which account for more than half of all cases worldwide [2]. For Chinese males, liver cancer ranks second only to lung cancer in cancer-related deaths [3]. Of all types of liver cancer, hepatocellular carcinoma (HCC) is the most common type, and it is responsible for approximately 75-80% of all liver cancers. Hepatic resection and liver transplantation remain the first choice for patients with early- and mid-stage HCC. However, the 5-year recurrence rate of HCC is approximately 80% after radical treatment [4]. Most patients are diagnosed with late-stage HCC, who have limited therapeutic options and generally have expected median survival times of 6-8 months [5]. Therefore, more research is urgently needed to comprehensively understand the molecular mechanisms of HCC and identify new biomarkers and therapeutics for HCC patients.

Long non-coding RNAs (lncRNAs) are a class of endogenous non-coding RNAs that are longer than 200 nucleotides and exert pivotal effects in the pathogenesis and progression of cancers [6,7,8]. Accumulating studies have shown that lncRNAs play important roles in regulating various biological processes, including chromatin and genome dynamics, gene expression, development, and cell differentiation [9,10,11]. Many dysregulated lncRNAs are involved in the tumorigenesis and development of HCC. For example, lncRNA PCNAP1 accelerated hepatitis B virus (HBV) replication and hepatocarcinogenesis by modulating the miR-154/PCNA/HBV cccDNA signaling pathway [12]. LINC00662 was upregulated in HCC and promoted HCC progression by activating the Wnt/β-catenin signaling pathway and M2 macrophage polarization [13]. Another lncRNA, p53-stabilizing and activating RNA (PSTAR), was downregulated in HCC and suppressed HCC cell proliferation and tumorigenicity by inducing p53-mediated cell cycle arrest [14]. However, the mechanisms of HCC are not clear, and many dysregulated lncRNAs must be explored.

Compared to protein-coding genes, lncRNAs have more tissue-specific expression characteristics [6]. Many lncRNAs with liver-specific expression patterns were identified as significantly related to the occurrence and progression of HCC. For example, lncRNA HULC (highly upregulated in liver cancer) accelerates HCC progression and attenuates the sensitivity of HCC cells to chemotherapeutic agents [15,16,17,18]. LINC01093 is downregulated in HCC tissues and suppresses HCC growth and metastasis [19]. Our previous studies also identified three liver-specific lncRNAs, FAM99B [20], LINC02499 [21], and LINC01146 [22], that were all downregulated in HCC and exerted similar inhibitory effects on the proliferation, migration, and invasion of HCC cells.

The current study identified a novel liver-specific lncRNA, FAM99A (family with sequence similarity 99 member A), which was specifically expressed in normal liver tissues based on the RNA-sequencing data from the Genotype-Tissue Expression (GTEx) project (https://gtexportal.org/home/) (Additional file 1: Figure S1). Given the great importance of liver-specific lncRNAs in the development of HCC, we comprehensively examined the expression level of FAM99A in HCC tissues based on public online databases and our HCC cohort. The clinical significance and prognostic value of FAM99A in HCC patients were also investigated. The biological function of FAM99A on HCC cell proliferation, migration, invasion, and tumor growth was evaluated in vitro and in vivo. RNA pull-down assay and mass spectrometry analysis were also performed to investigate the potential molecular mechanisms of FAM99A in impeding the progression of HCC.

Methods

Human tissue specimens

The Ethics Committee of Guangxi Medical University approved the study. HCC tissues and corresponding paracancerous tissues were collected from 62 HCC patients who underwent radical surgical resection at the Affiliated Cancer Hospital of Guangxi Medical University between February 2016 and December 2019. The diagnosis of HCC was confirmed by pathological examination. The tissue samples were snap frozen in liquid nitrogen and stored in liquid nitrogen until use. Patients with a history of preoperative chemotherapy or radiotherapy were excluded. Informed consent was obtained from all included patients, and patients were followed up until December 2021.

Public online databases

The RNA sequencing (RNA-seq) data (371 HCC tissues and 50 adjacent normal tissues) of HCC patients (level 3) were extracted from The Cancer Genome Atlas (TCGA) database (https://cancergenome.nih.gov/) (up to January 14, 2019). Transcript expression data were calculated as transcripts per million (TPM) and normalized by converting to log2 (TPM + 1). The clinical parameters and follow-up information were also downloaded to assess the clinical significance of FAM99A. For overall survival (OS) analysis, 370 cases were included after eliminating one patient without follow-up OS data.

The RNA-seq and clinical data of another HCC cohort were obtained from the International Cancer Genome Consortium (ICGC) data portal (https://dcc.icgc.org/) (Data Release 28) on June 2, 2020. Due to a lack of paracancerous tissue expression data, the Liver Cancer, France [LICA-FR] cohort was excluded from our study. For the Liver Cancer, RIKEN, Japan [LIRI-JP] cohort, a total of 221 HCC tissues and 200 adjacent normal tissues were enrolled in the research after eliminating two metastatic tumor cases. The normalized read counts in this cohort were used for analysis, and OS analysis was performed to evaluate the prognostic significance of FAM99A.

To draw a comprehensive result, we also extracted microarray datasets containing FAM99A expression data from the Gene Expression Omnibus (GEO) database (http://www.ncbi.nlm.nih.gov/geo/) (up to October 4, 2020). The retrieved keywords are listed below: (long non-coding RNAs OR lncRNAs OR non-coding RNAs) AND (hepatocellular carcinoma OR HCC OR liver) AND (cancer OR tumor OR carcinoma OR neoplasm* OR malignant*). The “Top Organisms” was restricted to “Homo sapiens”. Studies that simultaneously met the following two inclusion criteria were included: 1) HCC tissues and peritumoral liver specimens were included in the study (≥ 5 samples each group); and 2) the RNA profiling included the expression data of FAM99A. The basic information and FAM99A expression were carefully extracted from eligible datasets. A meta-analysis based on GEO datasets was subsequently performed.

Cell culture and cell transfection

The human HCC cell lines (Huh-7, Hep 3B, HepG2, HCCLM3, MHCC97L and MHCC97H) were maintained in our laboratory, and their background information has been previously described in detail [20]. Cells were cultured in Dulbecco’s modified Eagle’s medium (DMEM, Gibco, USA) or Minimum Essential Medium (MEM, Gibco, USA) supplemented with 10% fetal bovine serum (FBS, Gibco, Australia) at 37 °C in a humidified incubator containing 5% CO2.

To obtain FAM99A overexpression cell lines, full-length FAM99A (1430 bp) was ligated into the LV5 (EF-1aF/GFP&Puro) vector. A nonsense oligonucleotide was used as a negative control (Lv-NC). The lentiviruses were synthesized by GenePharma Co., Ltd. (Shanghai, China). For FAM99A knockdown, short hairpin RNA (shRNA) targeting FAM99A was inserted into the GV493 (hU6-MCS-CBh-gcGFP-IRES-puromycin) vector. A non-silencing shRNA was used as a negative control (sh-NC). The lentiviruses were synthesized by Genechem Co., Ltd. (Shanghai, China). The following shRNA sequences were used: sh-NC, 5’-TTCTCCGAACGTGTCACGT-3’; sh-FAM99A, 5’-AATAAAAGTCACAGGACAA-3’. Cells (8 × 104 cells/well) were seeded into 6-well plates, incubated for 24 hours, then infected with lentiviruses according to the manufacturer’s instructions. Seventy-two hours after infection, the cells were exposed to puromycin (3.5 μg/ml) for two weeks to obtain stably transfected cell lines. The overexpression and knockdown efficiency of FAM99A in HCC cell lines were verified using qRT-PCR.

RNA extraction, reverse transcription, and qRT-PCR

Total RNA of liver cancer cells and HCC samples was extracted using TRIzol reagent (Invitrogen, USA), and 900 ng of total RNA was reverse transcribed to complementary DNA using the PrimeScript™ RT reagent Kit with gDNA Eraser (Takara, Japan). The relative RNA expression level was evaluated using the TB GreenTM Premix Ex TaqTM II Kit (Takara, Japan) in a real-time PCR system (Applied Biosystems StepOnePlus, USA). GAPDH was used as an internal control. The transcript level was calculated using the 2Ct method, and the primers for GAPDH and FAM99A are listed below:

GAPDH-forward: 5’-AGCCACATCGCTCAGACAC-3’,

GAPDH-reverse: 5’-GCCCAATACGACCAAATCC-3’;

FAM99A-forward: 5’-CTCTTGTCCAGGTCAGCATCTC-3’,

FAM99A-reverse: 5’-ACGCATCACAAAACAGCCAC-3’.

Subcellular fractionation analysis

Isolation of cytoplasmic and nuclear RNAs in Hep 3B cells was performed using the PARIS Kit (Invitrogen, USA) according to the manufacturer’s instructions. The relative expression of FAM99A in cytoplasmic and nuclear fractions was determined using qRT-PCR. GAPDH and U6 were used as cytoplasmic and nuclear controls, respectively. The following primers were used for U6: U6-forward: 5’-CGCTTCGGCAGCACATATA-3’; U6-reverse: 5’-TTCACGAATTTGCGTGTCAT-3’.

Cell viability and colony formation assays

For the Cell Counting Kit-8 (CCK-8, Dojindo, Japan) cell proliferation viability assay, stably transfected Huh-7 (3000 cells/well) and Hep 3B (2000 cells/well) cells were seeded into 96-well plates. A total of 10 μl of CCK-8 reagent was added to each well after 1, 2, 3, 4, and 5 days of incubation, and the absorbance was measured at 450 nm after incubation at 37 °C for another 2 hours. For the plate clone formation assay, stably transfected Huh-7 (1500 cells/well) and Hep 3B (500 cells/well) cells were seeded into 6-well plates and maintained in medium supplementing with 10% FBS for 14 days (Huh-7) or 10 days (Hep 3B). Colonies were fixed with methanol for 15 minutes, and stained with 0.1% crystal violet for another 15 minutes.

Cell migration and invasion assays

Transwell inserts (Costar, Corning, USA) with 8-μm polycarbonate membranes were used to assess the migration ability of cancer cells. The upper chambers coated with Matrigel (Corning, USA) were used for the invasion assay in vitro. Briefly, 200 μl serum-free medium containing 1×105 cells was added to the upper chambers. The bottom chambers received 600 μl culture medium containing 10% FBS. After incubating for 1 day (Huh-7) or 2 days (Hep 3B), the cells remaining in the upper chambers were gently wiped with wet cotton swabs. The migrated or invaded cells were fixed with methanol and stained with 0.1% crystal violet for 15 minutes. Five fields were randomly selected for imaging using the microscope EVOS FL Auto Cell Imaging System (EVOS FL, Thermo Fisher Scientific, USA).

In vivo tumorigenesis experiments

BALB/c nude mice (4 weeks, male) were provided by the Experimental Animal Center of Guangxi Medical University and raised under specific pathogen free (SPF) conditions. Nude mice were randomly divided into 2 groups (n = 8 per group). A total of 5×106 Huh-7 cells stably transfected with FAM99A overexpression (Lv-FAM99A) or negative control (Lv-NC) were resuspended in a mixture of 50 μl PBS and 50 μl Matrigel (Corning, USA) and subcutaneously injected into the right flank of mice. The tumor size was measured every 3 days using a Vernier Caliper, and the tumor volume was calculated by the formula: tumor volume (mm3) = 0.5×L×W2 (L, longest diameter; W: shortest diameter). Four weeks later, the nude mice were sacrificed, and their subcutaneous tumors were isolated and weighed. The tumors were divided into two parts. One part was frozen in liquid nitrogen for RNA extraction, and the other part was fixed using 4% paraformaldehyde for hematoxylin-eosin (HE) and immunohistochemical (IHC) staining. All animal experiments were performed in accordance with the Guiding Principles for Care and Use of Experimental Animals, and the Animal Care & Welfare Committee of Guangxi Medical University approved this study (approval number: 202012016).

HE & IHC

Paraformaldehyde-fixed tissues were trimmed, dehydrated in gradient alcohol and embedded in paraffin. After deparaffinization with xylene and rehydration in an ethanol gradient, paraffin sections were stained with hematoxylin for 5~10 minutes followed by several dips in 1% hydrochloric acid alcohol. After rinsing with distilled water, the sections were stained with a 1% eosin aqueous solution for 3 minutes, dehydrated in gradient alcohol and cleared with xylene.

For IHC staining, citrate buffer was used for antigen retrieval, and 3% hydrogen peroxide was used to block the activity of endogenous peroxidase. The paraffin sections were washed with PBS and incubated with the diluted primary antibody anti-Ki67 (1:200, Abcam, ab16667) for 1 hour at 37 °C. After incubation with horseradish peroxidase (HRP)-labeled secondary antibody for 30 minutes, sections were stained with diaminobenzidine (DAB) chromogen and counterstained with hematoxylin. The sections were scanned using the TissueFAXS PLUS system (TissueGnostic, Austria), and the sum of integrated optical density (IOD) was calculated using Image-Pro Plus 6.0 software.

RNA fluorescence in situ hybridization (RNA-FISH)

Cy3-labeled oligonucleotide probes specifically targeting FAM99A were designed and synthesized by Genechem (Shanghai, China). Human 18S FISH probe mix and U6 FISH probe mix were synthesized by RiboBio (Guangzhou, China). RNA-FISH was accomplished using a RiboTM Fluorescent in Situ Hybridization Kit (RiboBio, Guangzhou, China) according to the manufacturer’s protocol. Briefly, Hep 3B and Huh-7 cells were cultured on cell climbing slices, washed with PBS, and fixed with 4% paraformaldehyde for 10 minutes. After washing with PBS three times, the cells were permeabilized with 0.5% Triton X-100 for 5 minutes at 4 °C. The permeabilized cells were incubated with pre-hybridization buffer for 30 minutes and incubated with hybridization buffer containing probes overnight at 37 °C. Nuclei were counterstained with 4′-6-diamidino-2-phenylindole (DAPI) for 10 minutes, and cells were washed with PBS three times. Cells were visualized, and images were captured using a confocal microscope (LSM 800, Zeiss, Germany).

RNA pull-down assay and mass spectrometry analysis

RNA pull-down was performed using the Pierce Magnetic RNA-Protein Pull-Down Kit (Thermo Fisher Scientific, USA). The sense and antisense chains of FAM99A were transcribed in vitro using T7 RNA Polymerase (Roche, USA), and labeled with Biotin RNA Labeling Mix (Roche, USA) according to the manufacturer’s instructions. The biotin-labeled RNAs were incubated with streptavidin magnetic beads (Invitrogen, USA) at 4 °C overnight to obtain the bead-RNA complex. Cell lysates (approximately 1 mg protein) were added to the bead-RNA complexes and incubated at room temperature for 1 hour. The precipitated complexes were washed with washing buffer three times and boiled in SDS buffer. The retrieved RNA-binding proteins were separated using electrophoresis and visualized by silver staining. The eluted proteins were also collected to performed mass spectrometry analysis using a Q Exactive System (Thermo Fisher Scientific, USA).

Statistical analysis

Numerical data are presented as the means ± standard deviation (SD). Paired-samples and independent-samples t tests were implemented to compare differences between two groups. The expression level of FAM99A was divided into high expression and low expression groups according to their median values. The chi-squared test was performed to explore the clinicopathological significance of FAM99A. The OS between the high and low FAM99A expression groups was compared using Kaplan–Meier curves with the log-rank test (two sides). Univariate and multivariate Cox proportional hazards regression models were also applied to determine the independent prognostic factors of HCC patients. All of the above analyses were performed using SPSS 20.0 software, and the results with P < 0.05 were considered statistically significant. The graphs were drawn using GraphPad Prism 8.0 software.

To integrate the microarray results of GEO datasets, meta-analysis was implemented using STATA 15.0 software. The chi-squared and I-squared tests were used to examine the heterogeneity between the included datasets. When I2 > 50% or P < 0.05, significant heterogeneity existed among datasets, and the random-effect model was selected. Otherwise, the fixed-effect model was used for meta-analysis. A forest plot was drawn to obtain the standardized mean difference (SMD) and 95% confidence interval (CI). When P < 0.05, a pooled SMD < 0 indicated that the expression level of FAM99A was significantly downregulated in HCC tissues compared to adjacent normal tissues. In contrast, FAM99A was considered upregulated in HCC tissues when P < 0.005 and the combined SMD > 0.

Results

FAM99A is downregulated in HCC and negatively correlated with poor prognosis of HCC patients based on TCGA and ICGC databases

To investigate the role of FAM99A in HCC development, we extracted the RNA-seq data of FAM99A from TCGA database. The results indicated that the expression level of FAM99A in HCC tissues was significantly lower than noncancerous tissues (2.39 ± 2.14 vs. 4.95 ± 0.73; P < 0.001; Fig. 1A). The patients were divided into high (> 1.619) and low expression (≤ 1.619) groups according to the median value of FAM99A. The clinical significance of FAM99A in HCC patients was also examined. The results revealed that low FAM99A expression significantly correlated with vascular invasion (P < 0.001) and advanced histological grade (P = 0.004; Table 1). The Kaplan–Meier survival curve showed that HCC patients with lower FAM99A expression tended to have poor OS (χ2 = 4.199, P = 0.040; Fig. 1B).

Fig. 1
figure 1

FAM99A is frequently downregulated in HCC tissues based on public online databases. A The expression level of FAM99A in 371 HCC tissues and 50 noncancerous tissues based on TCGA database. B Kaplan-Meier analysis showed that low expression of FAM99A was related to poor overall survival of HCC patients in TCGA database (n = 370). C The expression level of FAM99A in 221 HCC tissues and 200 adjacent normal tissues based on the ICGC data portal. D Kaplan-Meier analysis revealed that downregulated expression of FAM99A remarkably correlated with shorter overall survival of HCC patients in the ICGC data portal (n = 221)

Table 1 The relationship between FAM99A expression and clinical parameters in TCGA dataset

We also examined FAM99A expression in the LIRI-JP cohort based on the ICGC database. The results showed that FAM99A was downregulated in HCC tissues versus adjacent normal tissues (8.93 ± 16.21 vs. 19.83 ± 13.65; P < 0.001; Fig. 1C). To examine the prognostic significance of FAM99A, 221 patients were divided into high and low expression groups based on the median (2.00) of FAM99A. Kaplan–Meier survival analysis demonstrated that patients in the low FAM99A expression group had a shorter OS time than patients in the high FAM99A expression group (χ2 = 13.495, P < 0.001; Fig. 1D).

Meta-analysis of FAM99A expression in HCC based on the GEO database

To support a comprehensive conclusion, we also retrieved microarray chip data containing FAM99A expression from the GEO database. Based on our inclusion criteria, 18 eligible GEO datasets were enrolled, and detailed information of these GEO datasets is listed in Table 2. As shown in Fig. 2, there was significant heterogeneity among these 18 GEO datasets (I2 = 87.7%; P < 0.001). Therefore, a random-effect model was applied for the meta-analysis. The pooled SMD suggested that the expression of FAM99A was decreased in HCC tissues compared to noncancerous tissues (SMD = -1.162, 95% CI (-1.541, -0.783); P < 0.001; Fig. 2).

Table 2 The basic information and FAM99A expression level of included GEO datasets
Fig. 2
figure 2

A forest plot of 18 GEO datasets for the FAM99A expression level in HCC tissues (using a random-effect model)

The decreased expression and clinical significance of FAM99A in HCC were verified based on our HCC cohort

To further confirm the expression level of FAM99A in HCC, qRT-PCR was performed to compare FAM99A expression between 62 pairs of HCC and corresponding peritumoral liver specimens. We found that FAM99A was remarkably downregulated in 98.39% (61/62) of HCC samples (1.84 ± 3.62 vs. 22.27 ± 14.27; P < 0.001; Fig. 3A, B). To investigate the clinical significance of FAM99A in HCC, patients were also divided into a high expression group (> 0.417) and a low expression group (≤ 0.417) on the basis of the median expression of FAM99A. Consistent with the TCGA results, low expression of FAM99A significantly correlated with microvascular invasion (P = 0.041; Table 3).

Fig. 3
figure 3

The decreased expression of FAM99A in HCC was verified based on our HCC cohort. A, B Relative expression level of FAM99A measured using qRT-PCR in 62 pairs of HCC tissues and corresponding adjacent normal tissues. C Kaplan–Meier analysis demonstrated that decreased FAM99A expression was remarkably associated with reduced overall survival of HCC patients in our cohort (n = 62)

Table 3 Correlation of FAM99A expression with clinicopathologic features of patients in our HCC cohort.

Kaplan–Meier survival curve was also performed, and the results showed that decreased FAM99A expression was remarkably associated with reduced OS of patients in our cohort (χ2 = 6.658, P = 0.010; Fig. 3C). The results of univariate Cox regression analysis found that tumor size, BCLC stage, Child-Pugh classification grade, microvascular invasion, Edmondson-Steiner grade, and FAM99A expression level correlated with the OS of HCC patients (P < 0.1; Table 4). Further multivariate Cox regression analysis demonstrated that the expression level of FAM99A was an independent prognostic factor for HCC patients (hazard ratio: 0.425, 95% CI: 0.189–0.958, P = 0.039; Table 4).

Table 4 Univariate and multivariate Cox proportional hazards regression models of overall survival in 62 HCC patients

FAM99A suppressed HCC cell viability and clonogenicity in vitro

To determine the function of FAM99A in HCC cell lines, qRT-PCR was performed to examine the expression level of FAM99A in six HCC cell lines, and the results are shown in Fig. 4A. Hep 3B and Huh-7 cells with relatively high expression of FAM99A were selected for gain-of-function and loss-of-function analyses. The qRT-PCR results showed that FAM99A was successfully overexpressed in the Lv-FAM99A group in Hep 3B (P < 0.001) and Huh-7 cells (P < 0.01; Fig. 4B) but was downregulated in the sh-FAM99A group in Hep 3B (P < 0.001) and Huh-7 cells (P < 0.01; Fig. 4C) compared to the corresponding negative control group.

Fig. 4
figure 4

FAM99A suppresses the proliferation ability of HCC cells in vitro. A Relative expression level of FAM99A in five HCC cell lines. B The overexpression efficiency of Lv-FAM99A was examined using qRT-PCR in Hep 3B and Huh-7 cells. C The knockdown efficiency of sh-FAM99A was examined using qRT-PCR in Hep 3B and Huh-7 cells. D CCK-8 assay showed that upregulation of FAM99A inhibited the proliferation of Hep 3B and Huh-7 cells. E CCK-8 assay indicated that downregulated FAM99A promoted the proliferation of Hep 3B and Huh-7 cells. F Colony formation assays showed that overexpressing FAM99A reduced the number of colonies. G Colony formation assays revealed that knocking down of FAM99A increased the number of colonies. All experiments were repeated three times. *P < 0.05; **P < 0.01; ***P < 0.001

The CCK-8 assay showed that the overexpression of FAM99A significantly inhibited the viability of Hep 3B and Huh-7 cells (both Ps < 0.05; Fig. 4D). In contrast, FAM99A knockdown evidently facilitated the viability of Hep 3B and Huh-7 cells (both Ps < 0.01; Fig. 4E). The colony formation assay also showed that overexpression of FAM99A remarkably impeded the clonogenicity of Hep 3B (P < 0.05) and Huh-7 cells (P < 0.001; Fig. 4F), but the colonies were significantly increased when FAM99A was knocked down in Hep 3B (P < 0.01) and Huh-7 cells (P < 0.05; Fig. 4G). These results indicated that FAM99A suppressed the proliferation of HCC cells.

FAM99A impeded the migration and invasion capacities of HCC cells in vitro

The effects of FAM99A on the migration and invasion of HCC cells were further investigated using Transwell migration and invasion assays. Upregulation of FAM99A significantly restrained the cell migration and invasion abilities of Hep 3B and Huh-7 cells (all Ps < 0.001; Fig. 5A, B). In contrast, silencing FAM99A remarkably enhanced the cell migration and invasion activities of Hep 3B and Huh-7 cells (all Ps < 0.001; Fig. 5C, D). These findings suggested that FAM99A inhibited the migration and invasion behavior of HCC cells.

Fig. 5
figure 5

FAM99A impedes the migration and invasion activities of HCC cells in vitro. A, B Overexpression of FAM99A suppressed the migration and invasion abilities of Hep-3B and Huh-7 cells. Scale bar: 400 μm. C, D FAM99A knockdown strengthened the migration and invasion abilities of Hep-3B and Huh-7 cells. Scale bar: 400 μm. ***P < 0.001.

FAM99A inhibited HCC cell tumor growth in vivo

To test whether FAM99A suppressed tumor growth in vivo, a subcutaneous xenograft tumor model was constructed in nude mice (Fig. 6A, B). As shown in Fig. 6C and D, the tumor volume and tumor weight of the Lv-FAM99A group were lower than the control group, which demonstrated that FAM99A markedly impeded the tumor growth of Huh-7 cells in vivo (both Ps < 0.05). The FAM99A expression level in solid tumors was confirmed using qRT-PCR, and the results revealed that FAM99A expression was significantly higher in the Lv-FAM99A group than the control group (P < 0.001; Fig. 6E). HE staining suggested that upregulation of FAM99A decreased cell necrosis and infiltration of tumor tissues compared to the negative control group (Fig. 6F). The staining intensity of the proliferation marker Ki-67 in the Lv-FAM99A group was evidently weakened compared to the control group (P < 0.001; Fig. 6G). These results confirmed that FAM99A suppressed HCC cell growth in vivo.

Fig. 6
figure 6

FAM99A inhibits the tumor growth of HCC cells in vivo. A, B Construction of the subcutaneous xenograft tumor model in nude mice. C The growth curves of tumor volume in nude mice with negative control and FAM99A overexpression. D The weights of subcutaneous tumors in nude mice with negative control and FAM99A overexpression. E Relative expression levels of FAM99A in subcutaneous tumors were verified using qRT-PCR. F Representative images of hematoxylin-eosin (HE) staining and Ki-67 immunohistochemical (IHC) staining from the tumor sections. G Sum of integrated optical density (IOD) of Ki-67 staining in the tumor sections. *P < 0.05; ***P < 0.001.

FAM99A may interact with seven important proteins to regulate HCC progression

The subcellular localization of lncRNAs determines the dominant mechanism of its molecular functions. We examined the subcellular localization of FAM99A using the lncATLAS website, which is a comprehensive resource of lncRNA localization in human cells based on RNA-seq datasets [23] (https://lncatlas.crg.eu/). The results suggested that FAM99A was primarily located in the nuclei of HepG2 cells and showed a characteristic liver-specific expression pattern (Fig. 7A). Subcellular fractionation and FISH assays were performed to confirm the results from the online database. Subcellular fractionation analysis indicated that FAM99A was mostly located in the nucleus of Hep 3B cells (93.16% in the nucleus vs. 6.84% in the cytoplasm; Fig. 7B). As shown in Fig. 7C, the FISH assay also suggested that FAM99A was primarily localized in the nuclei of Hep 3B and Huh-7 cells.

Fig. 7
figure 7

The subcellular localization of FAM99A. A The subcellular localization of FAM99A in human cells based on the lncATLAS website (https://lncatlas.crg.eu/). B Relative distribution proportion of FAM99A in nuclear and cytoplasmic fractions in Hep 3B cells was determined using qRT-PCR. C RNA-FISH assay showed the enrichment of FAM99A in the nuclei of Hep 3B and Huh-7 cells. Scale bar: 50 μm.

An RNA pull-down assay was performed to identify proteins that may combine with FAM99A, and mass spectrometry analysis identified 266 differential proteins of FAM99A compared with its antisense strand (Additional file 2: Table S1; Fig. 8A, B). These 266 differential proteins were imported into DAVID 6.8 [24] (https://david.ncifcrf.gov/), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment [25,26,27] was performed. The results suggested that these proteins were primarily enriched in the “Biosynthesis of antibiotics”, “Carbon metabolism”, “Ribosome”, “Complement and coagulation cascades”, “Biosynthesis of amino acids”, “Pentose phosphate pathway”, “Glycolysis/Gluconeogenesis”, and “Focal adhesion” pathways (Benjamin adjust P < 0.05; Table 5). We also predicted the RNA binding proteins of FAM99A using the RBPmap website [28] (http://rbpmap.technion.ac.il/). We found seven proteins that were recognized in the RNA pull-down and RBPmap results, including PCBP1 (Poly (rC) binding protein 1), SRSF5 (serine/arginine-rich splicing factor 5), SRSF6 (serine/arginine-rich splicing factor 6), YBX1 (Y-box-binding protein 1), IGF2BP2 (insulin-like growth factor 2 mRNA-binding protein 2), HNRNPK (heterogeneous nuclear ribonucleoprotein K), and HNRNPL (heterogeneous nuclear ribonucleoprotein L) (Fig. 8C).

Fig. 8
figure 8

The RNA pull-down assay and mass spectrometry analysis of FAM99A. A Silver staining of the indicated electrophoretic bands of proteins pulled down by biotin-labeled FAM99A in Hep 3B cells. B Venn diagram of the proteins pulled down by FAM99A (FAM99A+) and antisense FAM99A (FAM99A-). There were 266 differential proteins only identified by FAM99A+. C Venn diagram of the potential binding proteins of FAM99A identified using the RNA pull-down assay and the predicted result of the RBPmap website (http://rbpmap.technion.ac.il/)

Table 5 Significant KEGG pathways of 266 RNA-pull down proteins of FAM99A

Discussion

Liver-specific lncRNAs are emerging as pivotal regulators in the tumorigenesis and progression of HCC. We identified for the first time that FAM99A exhibited a highly liver-specific expression pattern, which suggests that FAM99A plays important roles in HCC. FAM99A was steadily downregulated in HCC tissues and negatively correlated with vascular invasion and advanced histological grade of HCC patients. The Kaplan-Meier curve analysis also revealed that patients with low expression of FAM99A tended to have poor OS based on the data from the TCGA database, ICGC database, and our HCC cohort. A multivariate Cox regression model demonstrated that FAM99A was an independent factor for the OS of HCC patients. In vitro experiments indicated that overexpression of FAM99A inhibited the cell proliferation, migration, and invasion abilities of HCC cell lines, and knocking down FAM99A produced the opposite effects. The subcutaneous tumor formation model suggested that FAM99A suppressed tumor growth of HCC cells in vivo. Our findings suggested that FAM99A exerted a cancer-inhibiting effect in HCC progression.

Patients with malignant clinical features tend to have a poor prognosis in HCC. Our study found that patients with reduced expression of FAM99A were more likely to develop microvascular invasion and advanced histological grade. Another study also revealed that the downregulation of FAM99A was significantly associated with incomplete tumor capsule, tumor differentiation, recurrence, and poor prognosis of HCC patients [29]. These two studies revealed that FAM99A was an independent prognostic indicator for the OS of HCC patients. These results suggested that decreased expression of FAM99A contributed to the progression of HCC and may be a promising prognostic indicator for HCC patients.

To the best of our knowledge, we are the first to identify FAM99A as a liver-specific lncRNA based on the GTEx project. A series of gain- and loss-of-function studies suggested that FAM99A suppressed the proliferation, migration, and invasion of HCC cells. Recent studies also reported the antitumorigenic functions of FAM99A in HCC. BX Zhao et al. found that FAM99A inhibited HCC metastasis and epithelial-mesenchymal transition by sponging miR-92a [29]. Another study demonstrated that FAM99A suppressed HCC cell viability and GLUT1-mediated glycolysis, which inhibited HCC progression [30]. Many abnormally expressed tissue-specific lncRNAs are involved in the tumorigenesis and progression of cancers. LINC00993 is a breast-specific lncRNA that is downregulated in triple-negative breast cancer (TNBC), and it suppressed the tumor growth of TBNC [31]. Testis developmental related 1 (TDRG1), also known as LINC00532, is expressed exclusively in the testis. TDRG1 is upregulated in testicular seminoma tissues and promotes the development, migration, and chemotherapy resistance of seminoma cells [32,33,34]. Prostate enriched lncRNA (PSLNR) is a prostate-specific lncRNA that inhibits prostate cancer progression via the p53-dependent pathway [35]. In contrast, another two prostate-specific lncRNAs (PCGEM1 and PCA3) are overexpressed in prostate cancer and promote the cell proliferation ability and inhibit the cell apoptosis ability of prostate cancer cells [36,37,38,39]. Notably, PCA3 is promising as a more efficiency diagnostic biomarker for prostate cancer than the currently used prostate-specific antigen [40, 41]. Given the important roles of tissue-specific lncRNAs in cancer, more studies are needed to further explore FAM99A as a diagnostic biomarker or therapeutic target for HCC patients.

The molecular mechanisms of lncRNAs primarily depend on its localization in cells. We revealed for the first time that FAM99A was primarily located in the nucleus using subcellular fractionation and RNA-FISH assays, which indicated that FAM99A may play key regulatory roles in pivotal nuclear processes, such as chromatin organization, transcriptional and post-transcriptional programs, subcellular structures, and nuclear structure organization. These nuclear processes required the interaction of lncRNAs with RNA binding proteins (RBPs) almost universally [42, 43]. Therefore, we performed an RNA pull-down assay and identified that FAM99A pulled down 266 proteins using mass spectrometry analysis. To reveal the underlying mechanisms of FAM99A in HCC, these 266 proteins were used in KEGG pathway enrichment analysis. The results demonstrated that these proteins were primarily involved in several important cancer-related pathways, including complement and coagulation cascades [44], the pentose phosphate pathway [45,46,47], glycolysis/gluconeogenesis [48], and focal adhesion [49, 50]. We also identified seven proteins (PCBP1, SRSF5, SRSF6, YBX1, IGF2BP2, HNRNPK, and HNRNPL) that may interact with FAM99A. These seven proteins play crucial roles in tumor occurrence and development. PCBP1 is a multifunctional RBP that regulates the alternative splicing, translation, and RNA stability of many cancer-related genes to exert its cancer-inhibiting effect [51,52,53]. SRSF5 and SRSF6 belong to the serine/arginine-rich (SR) protein family, which is an important class of splicing regulators. SRSF5 and SRSF6 play important roles in the development and progression of cancers [54,55,56,57]. YBX1, also known as YB-1, binds DNA and RNA, and it is closely related to various malignant phenotypes of cancer cells, including tumor cell proliferation, metastasis, angiogenesis, and drug resistance [58, 59]. IGF2BP2 is a member of the conserved oncofetal RNA-binding protein family that acts as a N6-methyladenosine (m6A) reader, and it is involved in the development and progression of various cancer types [60,61,62]. HNRNPK and HNRNPL belong to the heterogeneous nuclear ribonucleoprotein family and interact with tumor-associated lncRNAs to regulate the tumorigenesis and progression of various cancers, including HCC [63,64,65]. Given the great importance of these seven proteins, further validation and functional experiments are needed to reveal the relationships between FAM99A and these proteins.

To the best of our knowledge, the current study is the first study to fully examine the expression level and clinical and prognostic significance of FAM99A in HCC based on three public online databases and our own HCC cohort. The current study used meta-analysis to pool FAM99A expression in HCC based on 18 GEO datasets for the first time. The ICGC database was used for the first time to examine the expression level and survival significance of FAM99A in HCC. We also determined the subcellular localization of FAM99A in HCC cells for the first time and performed an RNA pull-down assay to enrich the mechanistic research of FAM99A in HCC development. Although we identified seven key RBPs that may interact with FAM99A, we were unable to verify their binding modes and mechanisms due to limitations of the experimental conditions and experimental levels. Therefore, further in-depth studies are needed to elucidate the comprehensive mechanisms of FAM99A in HCC development.

Conclusion

FAM99A is a liver-specific lncRNA that is downregulated in HCC and negatively associated with poor prognosis in HCC patients. It may exert its tumor-suppressing function via binding with some critical RBPs. Our study enriches the understanding of the effects of FAM99A in HCC and suggests that it may serve as a promising prognostic biomarker for HCC patients.

Availability of data and materials

Access to public data is described in the manuscript. Experimental data of our study is available from the corresponding author on reasonable request.

Abbreviations

lncRNAs:

Long non-coding RNAs

HCC:

Hepatocellular carcinoma

TCGA:

The Cancer Genome Atlas

ICGC:

International Cancer Genome Consortium

GEO:

Gene Expression Omnibus

qRT–PCR:

Quantitative real-time polymerase chain reaction

PSTAR:

p53-stabilizing and activating RNA

HULC:

Highly upregulated in liver cancer

GTEx:

Genotype-Tissue Expression

RNA-seq:

RNA sequencing

TPM:

Transcripts per million

OS:

overall survival

DMEM:

Dulbecco’s modified Eagle’s medium

MEM:

Minimum Essential Medium

FBS:

Fetal bovine serum

CCK-8:

Cell Counting Kit-8

SPF:

Specific pathogen free

HE:

Hematoxylin-eosin

IHC:

Immunohistochemical

FISH:

Fluorescence in situ hybridization

DAPI:

4′-6-diamidino-2-phenylindole

SD:

Standard deviation

SMD:

Standardized mean difference

CI:

Confidence interval

KEGG:

Kyoto Encyclopedia of Genes and Genomes

PCBP1:

Poly(rC) binding protein 1

SRSF5:

Serine/arginine-rich splicing factor 5

SRSF6:

Serine/arginine-rich splicing factor 6

YBX1:

Y-box-binding protein 1

IGF2BP2:

Insulin-like growth factor 2 mRNA-binding protein 2

HNRNPK:

Heterogeneous nuclear ribonucleoprotein K

HNRNPL:

Heterogeneous nuclear ribonucleoprotein L

TNBC:

Triple-negative breast cancer

RBPs:

RNA binding proteins

m6A:

N6-methyladenosine

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.

    Article  PubMed  Google Scholar 

  2. McGlynn KA, Petrick JL, London WT. Global epidemiology of hepatocellular carcinoma: an emphasis on demographic and regional variability. Clin Liver Dis. 2015;19:223–38.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Feng RM, Zong YN, Cao SM, Xu RH. Current cancer situation in China: good or bad news from the 2018 Global Cancer Statistics? Cancer Commun (Lond). 2019;39:22.

    Article  Google Scholar 

  4. Ishikawa T. Strategy for improving survival and reducing recurrence of HCV-related hepatocellular carcinoma. World J Gastroenterol. 2013;19:6127–30.

    Article  PubMed  PubMed Central  Google Scholar 

  5. European Association for the Study of the Liver. Electronic address, e. e. e., European Association for the Study of the, L. EASL Clinical Practice Guidelines: Management of hepatocellular carcinoma. J Hepatol. 2018;69:182–236.

  6. Bhan A, Soleimani M, Mandal SS. Long Noncoding RNA and Cancer: A New Paradigm. Cancer Res. 2017;77:3965–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Chi Y, Wang D, Wang J, Yu W, Yang J. Long Non-Coding RNA in the Pathogenesis of Cancers. Cells. 2019:8.

  8. Huang Z, Zhou JK, Peng Y, He W, Huang C. The role of long noncoding RNAs in hepatocellular carcinoma. Mol Cancer. 2020;19:77.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Constanty F, Shkumatava A. lncRNAs in development and differentiation: from sequence motifs to functional characterization. Development. 2021:148.

  10. Dykes IM, Emanueli C. Transcriptional and Post-transcriptional Gene Regulation by Long Non-coding RNA. Genomics Proteomics Bioinformatics. 2017;15:177–86.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Nair L, Chung H, Basu U. Regulation of long non-coding RNAs and genome dynamics by the RNA surveillance machinery. Nat Rev Mol Cell Biol. 2020;21:123–36.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Feng J, Yang G, Liu Y, Gao Y, Zhao M, Bu Y, et al. LncRNA PCNAP1 modulates hepatitis B virus replication and enhances tumor growth of liver cancer. Theranostics. 2019;9:5227–45.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Tian X, Wu Y, Yang Y, Wang J, Niu M, Gao S, et al. Long noncoding RNA LINC00662 promotes M2 macrophage polarization and hepatocellular carcinoma progression via activating Wnt/beta-catenin signaling. Mol Oncol. 2020;14:462–83.

    Article  CAS  PubMed  Google Scholar 

  14. Qin G, Tu X, Li H, Cao P, Chen X, Song J, et al. Long Noncoding RNA p53-Stabilizing and Activating RNA Promotes p53 Signaling by Inhibiting Heterogeneous Nuclear Ribonucleoprotein K deSUMOylation and Suppresses Hepatocellular Carcinoma. Hepatology. 2020;71:112–29.

    Article  CAS  PubMed  Google Scholar 

  15. Liu Y, Feng J, Sun M, Yang G, Yuan H, Wang Y, et al. Long non-coding RNA HULC activates HBV by modulating HBx/STAT3/miR-539/APOBEC3B signaling in HBV-related hepatocellular carcinoma. Cancer Lett. 2019;454:158–70.

    Article  CAS  PubMed  Google Scholar 

  16. Xin X, Wu M, Meng Q, Wang C, Lu Y, Yang Y, et al. Long noncoding RNA HULC accelerates liver cancer by inhibiting PTEN via autophagy cooperation to miR15a. Mol Cancer. 2018;17:94.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Xiong H, Ni Z, He J, Jiang S, Li X, He J, et al. LncRNA HULC triggers autophagy via stabilizing Sirt1 and attenuates the chemosensitivity of HCC cells. Oncogene. 2017;36:3528–40.

    Article  CAS  PubMed  Google Scholar 

  18. Zhang H, Liao Z, Liu F, Su C, Zhu H, Li Y, et al. Long noncoding RNA HULC promotes hepatocellular carcinoma progression. Aging (Albany NY). 2019;11:9111–27.

    Article  CAS  Google Scholar 

  19. He J, Zuo Q, Hu B, Jin H, Wang C, Cheng Z, et al. A novel, liver-specific long noncoding RNA LINC01093 suppresses HCC progression by interaction with IGF2BP1 to facilitate decay of GLI1 mRNA. Cancer Lett. 2019;450:98–109.

    Article  CAS  PubMed  Google Scholar 

  20. Mo M, Liu S, Ma X, Tan C, Wei L, Sheng Y, et al. A liver-specific lncRNA, FAM99B, suppresses hepatocellular carcinoma progression through inhibition of cell proliferation, migration, and invasion. J Cancer Res Clin Oncol. 2019;145:2027–38.

    Article  CAS  PubMed  Google Scholar 

  21. Ma X, Mo M, Tan HJJ, Tan C, Zeng X, Zhang G, et al. LINC02499, a novel liver-specific long non-coding RNA with potential diagnostic and prognostic value, inhibits hepatocellular carcinoma cell proliferation, migration, and invasion. Hepatol Res. 2020;50:726–40.

    Article  CAS  PubMed  Google Scholar 

  22. Ma X, Mo M, Tan C, Tan JHJ, Huang H, Liu B, et al. Liver-specific LINC01146, a promising prognostic indicator, inhibits the malignant phenotype of hepatocellular carcinoma cells both in vitro and in vivo. J Transl Med. 2022;20:57.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Mas-Ponte D, Carlevaro-Fita J, Palumbo E, Hermoso Pulido T, Guigo R, Johnson R. LncATLAS database for subcellular localization of long noncoding RNAs. RNA. 2017;23:1080–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Huang da W, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37:1–13.

    Article  PubMed  Google Scholar 

  25. Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 2019;28:1947–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Kanehisa M, Furumichi M, Sato Y, Ishiguro-Watanabe M, Tanabe M. KEGG: integrating viruses and cellular organisms. Nucleic Acids Res. 2021;49:D545–51.

    Article  CAS  PubMed  Google Scholar 

  27. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Zhao B, Ke K, Wang Y, Wang F, Shi Y, Zheng X, et al. HIF-1alpha and HDAC1 mediated regulation of FAM99A-miR92a signaling contributes to hypoxia induced HCC metastasis. Signal Transduct Target Ther. 2020;5:118.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Zheng X, Gou Y, Jiang Z, Yang A, Yang Z, Qin S. Icaritin-Induced FAM99A Affects GLUT1-Mediated Glycolysis via Regulating the JAK2/STAT3 Pathway in Hepatocellular Carcinoma. Front Oncol. 2021;11:740557.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Guo S, Jian L, Tao K, Chen C, Yu H, Liu S. Novel Breast-Specific Long Non-coding RNA LINC00993 Acts as a Tumor Suppressor in Triple-Negative Breast Cancer. Front Oncol. 2019;9:1325.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Gan Y, Wang Y, Tan Z, Zhou J, Kitazawa R, Jiang X, et al. TDRG1 regulates chemosensitivity of seminoma TCam-2 cells to cisplatin via PI3K/Akt/mTOR signaling pathway and mitochondria-mediated apoptotic pathway. Cancer Biol Ther. 2016;17:741–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Peng D, Wei J, Gan Y, Yang J, Jiang X, Kitazawa R, et al. Testis developmental related gene 1 regulates the chemosensitivity of seminoma TCam-2 cells to cisplatin via autophagy. J Cell Mol Med. 2019;23:7773–84.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Wang Y, Gan Y, Tan Z, Zhou J, Kitazawa R, Jiang X, et al. TDRG1 functions in testicular seminoma are dependent on the PI3K/Akt/mTOR signaling pathway. Onco Targets Ther. 2016;9:409–20.

    CAS  PubMed  PubMed Central  Google Scholar 

  35. Wang D, Wan X, Zhang Y, Kong Z, Lu Y, Sun X, et al. A novel androgen-reduced prostate-specific lncRNA, PSLNR, inhibits prostate-cancer progression in part by regulating the p53-dependent pathway. Prostate. 2019;79:1362–77.

    Article  CAS  PubMed  Google Scholar 

  36. Fu X, Ravindranath L, Tran N, Petrovics G, Srivastava S. Regulation of apoptosis by a prostate-specific and prostate cancer-associated noncoding gene, PCGEM1. DNA Cell Biol. 2006;25:135–41.

    Article  PubMed  Google Scholar 

  37. Petrovics G, Zhang W, Makarem M, Street JP, Connelly R, Sun L, et al. Elevated expression of PCGEM1, a prostate-specific gene with cell growth-promoting function, is associated with high-risk prostate cancer patients. Oncogene. 2004;23:605–11.

    Article  CAS  PubMed  Google Scholar 

  38. Srikantan V, Zou Z, Petrovics G, Xu L, Augustus M, Davis L, et al. PCGEM1, a prostate-specific gene, is overexpressed in prostate cancer. Proc Natl Acad Sci U S A. 2000;97:12216–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Zhang G, He X, Ren C, Lin J, Wang Q. Long noncoding RNA PCA3 regulates prostate cancer through sponging miR-218-5p and modulating high mobility group box 1. J Cell Physiol. 2019;234:13097–109.

    Article  CAS  PubMed  Google Scholar 

  40. Lee GL, Dobi A, Srivastava S. Prostate cancer: diagnostic performance of the PCA3 urine test. Nat Rev Urol. 2011;8:123–4.

    Article  PubMed  Google Scholar 

  41. Soares JC, Soares AC, Rodrigues VC, Melendez ME, Santos AC, Faria EF, et al. Detection of the Prostate Cancer Biomarker PCA3 with Electrochemical and Impedance-Based Biosensors. ACS Appl Mater Interfaces. 2019;11:46645–50.

    Article  CAS  PubMed  Google Scholar 

  42. Noh JH, Kim KM, McClusky WG, Abdelmohsen K, Gorospe M. Cytoplasmic functions of long noncoding RNAs. Wiley Interdiscip Rev RNA. 2018;9:e1471.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Sun Q, Hao Q, Prasanth KV. Nuclear Long Noncoding RNAs: Key Regulators of Gene Expression. Trends Genet. 2018;34:142–57.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Zhang J, Chen M, Zhao Y, Xiong H, Sneh T, Fan Y, et al. Complement and coagulation cascades pathway correlates with chemosensitivity and overall survival in patients with soft tissue sarcoma. Eur J Pharmacol. 2020;879:173121.

    Article  CAS  PubMed  Google Scholar 

  45. Ge T, Yang J, Zhou S, Wang Y, Li Y, Tong X. The Role of the Pentose Phosphate Pathway in Diabetes and Cancer. Front Endocrinol (Lausanne). 2020;11:365.

    Article  Google Scholar 

  46. Jiang P, Du W, Wu M. Regulation of the pentose phosphate pathway in cancer. Protein Cell. 2014;5:592–602.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Patra KC, Hay N. The pentose phosphate pathway and cancer. Trends Biochem Sci. 2014;39:347–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Grasmann G, Smolle E, Olschewski H, Leithner K. Gluconeogenesis in cancer cells - Repurposing of a starvation-induced metabolic pathway? Biochim Biophys Acta Rev Cancer. 2019;1872:24–36.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Eke I, Cordes N. Focal adhesion signaling and therapy resistance in cancer. Semin Cancer Biol. 2015;31:65–75.

    Article  CAS  PubMed  Google Scholar 

  50. Yam JW, Tse EY, Ng IO. Role and significance of focal adhesion proteins in hepatocellular carcinoma. J Gastroenterol Hepatol. 2009;24:520–30.

    Article  CAS  PubMed  Google Scholar 

  51. Huang S, Luo K, Jiang L, Zhang XD, Lv YH, Li RF. PCBP1 regulates the transcription and alternative splicing of metastasisrelated genes and pathways in hepatocellular carcinoma. Sci Rep. 2021;11:23356.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Wang X, Guo J, Che X, Jia R. PCBP1 inhibits the expression of oncogenic STAT3 isoform by targeting alternative splicing of STAT3 exon 23. Int J Biol Sci. 2019;15:1177–86.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Zhang T, Huang XH, Dong L, Hu D, Ge C, Zhan YQ, et al. PCBP-1 regulates alternative splicing of the CD44 gene and inhibits invasion in human hepatoma cell line HepG2 cells. Mol Cancer. 2010;9:72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Chen Y, Huang Q, Liu W, Zhu Q, Cui CP, Xu L, et al. Mutually exclusive acetylation and ubiquitylation of the splicing factor SRSF5 control tumor growth. Nat Commun. 2018;9:2464.

    Article  PubMed  PubMed Central  Google Scholar 

  55. She W, Shao J, Jia R. Targeting Splicing Factor SRSF6 for Cancer Therapy. Front Cell Dev Biol. 2021;9:780023.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Wan L, Yu W, Shen E, Sun W, Liu Y, Kong J, et al. SRSF6-regulated alternative splicing that promotes tumour progression offers a therapy target for colorectal cancer. Gut. 2019;68:118–29.

    Article  CAS  PubMed  Google Scholar 

  57. Yang S, Jia R, Bian Z. SRSF5 functions as a novel oncogenic splicing factor and is upregulated by oncogene SRSF3 in oral squamous cell carcinoma. Biochim Biophys Acta Mol Cell Res. 2018;1865:1161–72.

    Article  CAS  PubMed  Google Scholar 

  58. Alkrekshi A, Wang W, Rana PS, Markovic V, Sossey-Alaoui K. A comprehensive review of the functions of YB-1 in cancer stemness, metastasis and drug resistance. Cell Signal. 2021;85:110073.

    Article  PubMed  Google Scholar 

  59. Lasham A, Print CG, Woolley AG, Dunn SE, Braithwaite AW. YB-1: oncoprotein, prognostic marker and therapeutic target? Biochem J. 2013;449:11–23.

    Article  CAS  PubMed  Google Scholar 

  60. Bell JL, Wachter K, Muhleck B, Pazaitis N, Kohn M, Lederer M, et al. Insulin-like growth factor 2 mRNA-binding proteins (IGF2BPs): post-transcriptional drivers of cancer progression? Cell Mol Life Sci. 2013;70:2657–75.

    Article  CAS  PubMed  Google Scholar 

  61. Cao J, Mu Q, Huang H. The Roles of Insulin-Like Growth Factor 2 mRNA-Binding Protein 2 in Cancer and Cancer Stem Cells. Stem Cells Int. 2018;2018:4217259.

    Article  PubMed  PubMed Central  Google Scholar 

  62. Wang J, Chen L, Qiang P. The role of IGF2BP2, an m6A reader gene, in human metabolic diseases and cancers. Cancer Cell Int. 2021;21:99.

    Article  PubMed  PubMed Central  Google Scholar 

  63. Gu J, Chen Z, Chen X, Wang Z. Heterogeneous nuclear ribonucleoprotein (hnRNPL) in cancer. Clin Chim Acta. 2020;507:286–94.

    Article  CAS  PubMed  Google Scholar 

  64. Wang H, Ma P, Liu P, Guo D, Liu Z, Zhang Z. lncRNA SNHG6 promotes hepatocellular carcinoma progression by interacting with HNRNPL/PTBP1 to facilitate SETD7/LZTFL1 mRNA destabilization. Cancer Lett. 2021;520:121–31.

    Article  CAS  PubMed  Google Scholar 

  65. Wang Z, Qiu H, He J, Liu L, Xue W, Fox A, et al. The emerging roles of hnRNPK. J Cell Physiol. 2020;235:1995–2008.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We acknowledge the Guangxi Colleges and Universities Key Laboratory of Prevention and Control of Highly Prevalent Diseases, and the Guangxi Medical University Laboratory Animal Center for providing a good experimental environment. We also sincerely acknowledge the contributions from the TCGA, ICGC, and GEO databases.

Funding

This study was supported by the National Natural Science Foundation of China (81960613).

Author information

Authors and Affiliations

Authors

Contributions

MM, XM, YL, CT, BL, PT, QL contributed to data collection, performed the statistical analysis and experiments. SL, HY, DH assisted with the collection of clinical samples and data. MM wrote the first draft of the manuscript. XM and YL revised the manuscript. XQ and XZ designed the research and provided final approval of the manuscript. All the authors read and approved the final manuscript.

Author’s information

No applicable.

Corresponding authors

Correspondence to Xiaoyun Zeng or Xiaoqiang Qiu.

Ethics declarations

Ethics approval and consent to participate

Human tissue test: The informed consent was acquired from each patient and the study was approved by the Ethics Committee of Guangxi Medical University. All methods were carried out in accordance with the Declaration of Helsinki.

Animal models: The animal study was approved by the Animal Care &Welfare Committee of Guangxi Medical University. All methods were performed in accordance with relevant guidelines and regulations. This study was carried out in accordance with the ARRIVE guidelines.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Mo, M., Ma, X., Luo, Y. et al. Liver-specific lncRNA FAM99A may be a tumor suppressor and promising prognostic biomarker in hepatocellular carcinoma. BMC Cancer 22, 1098 (2022). https://doi.org/10.1186/s12885-022-10186-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12885-022-10186-2

Keywords