Research Paper Volume 12, Issue 12 pp 11843—11863
Serum exosomal long noncoding RNAs lnc-FAM72D-3 and lnc-EPC1-4 as diagnostic biomarkers for hepatocellular carcinoma
- 1 Department of General Surgery, The Third Affiliated Hospital of Sun Yat-sen University, Guangzhou 510000, China
- 2 Department of Cell-gene Therapy Translational Medicine Research Center, The Third Affiliated Hospital of Sun Yat-sen University, Guangzhou 510000, China
- 3 Department of Liver Disease Lab, The Third Affiliated Hospital of Sun Yat-sen University, Guangzhou 510000, China
- 4 Department of Hepatobilliary Surgery, The Third Affiliated Hospital of Sun Yat-sen University, Guangzhou 510000, China
Received: October 22, 2019 Accepted: May 1, 2020 Published: June 18, 2020
https://doi.org/10.18632/aging.103355How to Cite
Copyright © 2020 Yao et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Abstract
Long noncoding RNAs (lncRNAs), such as LINC00462, HOTAIR, and MALAT1, are significantly upregulated in hepatocellular carcinoma (HCC) tissues. However, lncRNA expression in the serum of HCC patients is still unclear. To identify candidate lncRNAs for HCC diagnosis, we purified exosomal total RNA from the serum of healthy volunteers (controls) and hepatitis, cirrhosis, and HCC patients. To assess the function of lncRNAs, small interfering RNAs and overexpression vectors were designed and cell viability and cell apoptosis assays conducted. The exosomes of the control group had a larger number of lncRNAs with a high amount of alternative splicing compared to hepatic disease patients. qPCR assays showed that lnc-FAM72D-3, lnc-GPR89B-15, lncZEB2-19, and lnc-EPC1-4 are differentially expressed in HCC. Furthermore, the expression level of lnc-EPC1-4 correlated with age. While the expression levels of lnc-GPR89B-15 and lnc-EPC1-4 correlated with serum alpha-fetoprotein level. lnc-FAM72D-3 knockdown decreased cell viability and promoted cell apoptosis, indicating that lnc-FAM72D-3 functions as an oncogene in HCC. In contrast, lnc-EPC1-4 overexpression inhibited cell proliferation and induced cell apoptosis, indicating that it functions as a tumor suppressor gene. Collectively, these findings show that lnc-FAM72D-3 and lnc-EPC1-4 play a novel role that might contribute to hepatocarcinogenesis and identify potential candidate biomarkers for HCC diagnosis.
Introduction
Liver cancer, also called hepatocellular carcinoma (HCC), is the fifth-most common cancer worldwide, with a mortality rate that increased by almost 3% per year from 2010 to 2014 [1, 2]. Surgical treatment of HCC might be effective in the initial stage of the disease. However, most HCC patients present with intra- and extrahepatic metastasis upon diagnosis, which limits their treatment options to radiotherapy and chemotherapy. Unfortunately, the treatment efficacy of radiotherapy and chemotherapy is limited [3]. Therefore, early diagnosis of HCC is of utmost importance. Currently, the most commonly used methods of diagnosing HCC involve imaging or detection of tumor biomarkers, such as alpha-fetoprotein (AFP), from the patient’s serum. However, these techniques are less sensitive for early diagnosis of HCC [4]. Therefore, it is important to identify novel diagnostic and prognostic molecular markers to improve early diagnosis of HCC.
Long noncoding RNAs (lncRNAs) are non-protein-coding transcripts with a length of >200 nucleotides. lncRNAs participate in various biological processes, such as tumorigenesis, metastasis, and proliferation [5–8], so they are considered diagnostic biomarkers for multiple cancers, such as gastric, bladder, colorectum, prostate, and renal cancers [9–12]. lncRNAs, such as LINC00462, HOTAIR, MALAT1, CCAT1, CCAT2, LINC00161, and SPRY4-IT1, are significantly upregulated in HCC tissues [13–19]. However, lncRNA expression in the serum of HCC patients is still unclear.
Circulating exosomes are small 30–150 nm membrane vesicles that are released into extracellular environments, such as serum or plasma [20]. They usually contain proteins and RNAs that perform various functions, such as cell–cell communication and signaling. Researchers have focused on identifying biomarkers in blood exosomes as a minimally invasive sample source that reflects the signatures of pathological conditions such as cancer and atherosclerosis without the need for surgical sample collection [21]. Although various microRNAs (miRNAs) and proteins are found inside exosomes, indicating that they could potentially be used as diagnostic biomarkers for HCC, little is known about the serum exosomal lncRNAs in HCC [22–24].
This study determined whether exosomal lncRNAs have functions in hepatic diseases and could potentially be used as diagnostic biomarkers. We purified exosomal lncRNAs from 45 healthy volunteers, 45 hepatitis patients, 45 cirrhosis patients, and 45 HCC patients. We conducted RNA sequencing (RNA-seq) to obtain expression profiles and performed correlative and functional analyses of candidate lncRNAs.
Results
Exosome particle and protein concentrations modified in hepatic diseases
To investigate exosomal lncRNAs involved in hepatic diseases, we collected serum from healthy volunteers and hepatitis, cirrhosis, and HCC patients. Exosomes were purified by ultracentrifugation, and their sizes were determined by NTA. As shown in Supplementary Figure 1A, particle diameters of exosomes ranged from ~50 to 150 nm for the four groups, and there was no significant difference in the size distribution. CD63, TSG101, and CD9, which are surface antigens commonly expressed in exosomes, were detected in all four groups (Supplementary Figure 1B). In contrast, calnexin protein expression was not detected in any group (Supplementary Figure 1B). These findings indicated successful purification of exosomes and also that they were of uniform size across the four groups.
To further evaluate potential differences between exosomes of healthy volunteers and hepatic disease patients, we measured the particle and protein concentrations of exosomes. The control and hepatitis groups had similar particle concentrations, but the particle concentration was 50% lower for the cirrhosis group and 20% higher for the HCC group (Table 1). The hepatitis, cirrhosis, and HCC groups also had higher protein concentrations compared to the control group, indicating that hepatic diseases induce changes in both particle concentration and protein concentration of exosomes.
Table 1. Comparison of serum exosome particle concentrations and protein concentrations in healthy volunteers and patients with hepatic diseases.
| Exosome grouping | Particle concentration (particles/ml) | Exosome protein concentration (μg/μl) | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| Normal | 3.40E+11 | 9.127 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| Hepatitis | 3.81E+11 | 11.345 * | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| Cirrhosis | 1.97E+11 * | 11.399 * | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| Hepatocellular carcinoma | 4.54E+11 * | 11.299 * | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| Compared with the control group, * p < 0.05 was considered to be significantly different. | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Bioinformatics study of exosomal lncRNAs
To further qualify differences between the exosomes of healthy volunteers and hepatic disease patients, we purified exosomal total RNA and conducted RNA-seq. RNAs were annotated, and lncRNAs were extracted from the data for further analysis. For all four groups, lncRNAs were mostly intronic. In the control group, we determined 946 relevant genes to be involved in lncRNA formation, with 1955 associated alternative splicing events (Figure 1A). In contrast, exosomes of the hepatitis, cirrhosis, and HCC groups had obviously fewer lncRNAs with fewer alternative splice events (ranging from 4- to 40-fold reduction), indicating that liver disease downregulates lncRNA production and alternative splicing of lncRNAs. However, exosomes of all four groups had similar length distributions for lncRNAs, the majority of which were ~500 nt (Figure 1B). The heatmap showed that compared to the control group, 34,175 lncRNAs were upregulated and 15,553 downregulated in the hepatitis group, 16,673 lncRNAs were upregulated and 22,033 downregulated in the cirrhosis group, and 28,475 lncRNAs were upregulated and 16,453 downregulated in HCC group (Figure 1C and Supplementary Tables 2–4). We also examined chromosomal distributions of lncRNAs. In the hepatitis, cirrhosis, and HCC groups, upregulated lncRNAs were most abundant on chromosomes 1 and 2, being about twice as much as downregulated lncRNAs. Interestingly, in the cirrhosis group, upregulated lncRNAs were also enriched on chromosomes 5 and 6 and were six times more than downregulated lncRNAs (Figure 1D), indicating that lncRNAs that play a role in liver disease are more distributed on chromosomes 1 and 2.
Figure 1. Bioinformatic analyses of lncRNAs. (A) Number of genes (red) and alternate splicing events (green) for lncRNAs from control, hepatitis, cirrhosis, and HCC groups. (B) Length distribution of lncRNAs in control, hepatitis, cirrhosis, and HCC groups. (C) Cluster heatmaps of lncRNA expression profiles in the control group compared to hepatitis, cirrhosis, and HCC groups. Gradient green represents downregulated lncRNAs, while gradient red represents upregulated lncRNAs. (D) Chromosomal distributions of lncRNAs in the control group compared to hepatitis, cirrhosis, and HCC groups. lncRNA, long noncoding RNA; HCC, hepatocellular carcinoma.
To determine the potential relevance of differentially expressed lncRNAs from the hepatitis, cirrhosis, and HCC groups, we conducted KEGG and GO analysis. In KEGG pathway annotation, the pathways were mainly RNA transport and protein processing in the endoplasmic reticulum (Figure 2A): Ras, Rap1, and phosphatidylinositol-3 kinase/protein kinase B (PI3K/Akt) signaling pathways, which were enriched in the control group compared to hepatitis and HCC groups (Figure 2C, 2E). In addition, GO analysis showed that significantly expressed lncRNAs are more focused on the “cell” and “cell part” of cellular components, the “binding” of molecular function, and the “cellular process” and “biology regulation” of biological process in the control group compared to hepatitis and HCC groups (Figure 2B, 2D, 2F), indicating that pathways in KEGG analysis and biological processes in GO analysis play important roles in the occurrence and development of liver diseases. Regulatory networks of common significantly differentially expressed lncRNAs revealed many clusters with low interconnectedness between the genes (Figure 3), indicating that the underlying mechanism is complicated and that each lncRNA might perform distinct functions in regulating specific genes.
Figure 2. Functional analyses of significantly expressed lncRNAs. (A, C, E) KEGG pathway enrichment of the differentially expressed lncRNAs in control, hepatitis, and HCC groups. (B, D, F) GO analysis of differentially expressed lncRNAs in control, hepatitis, and HCC groups, respectively. lncRNA, long noncoding RNA; HCC, hepatocellular carcinoma.
Figure 3. Regulation network of differentially expressed lncRNAs. lncRNA, long noncoding RNA.
Validation of the expression of differentially expressed lncRNAs
To further investigate the expression differences in exosomal lncRNAs of the four groups, we performed qPCR and examined the expression levels of nine lncRNAs. The expression levels of lnc-GPR89B-15 and lnc-ZEB2-19 verified by qPCR were consistent with RNA-seq results in hepatitis, cirrhosis, and HCC groups (Figure 4). lnc-FAM72D-3 expression in hepatitis, cirrhosis, and HCC groups was significantly high compared to the control group, while lnc-EPC1-4 expression in hepatitis, cirrhosis, and HCC groups was significantly low compared to the control group (Figure 4). Four lncRNAs (lnc-GPR89B-15, lnc-FAM72D-3, lncEPC1-4, and lnc-ZEB2-19) with differential expression in the HCC groups were selected for further evaluation.
Figure 4. Differentially expressed lncRNAs in exosomes of control, hepatitis, cirrhosis, and HCC groups validated by real-time qPCR. lncRNA, long noncoding RNA; HCC, hepatocellular carcinoma; qPCR, quantitative polymerase chain reaction.
Next, to determine the relationship between exosomal expression and diseases, we compared the expression levels of these lncRNAs between the control group and hepatitis, cirrhosis, and HCC groups (Table 2). As indicated by the area under the curve (AUC) of the receiver operating characteristic (ROC) curves, each of the four lncRNAs showed a close correlation with hepatic diseases Figure 5). Of them, lnc-ZEB2-19 had the highest AUC in all three comparisons, indicating a significant correlation between lnc-ZEB2-19 exosomal expression and hepatic diseases. Since most lncRNA expression levels were correlated with hepatitis, cirrhosis, and HCC to some extent, we combined the data from the patients and calculated the overall correlation between different patient characteristics and lncRNA expression (Table 3). The results suggested that the expression of lnc-EPC1-4 was significantly correlated with age, while the expression of lnc-GPR89B-15 and lnc-EPC1-4 was significantly correlated with AFP concentration (P < 0.05).
Figure 5. Clinical analysis of AUC of (A) lnc-FAM72D-3, (B) lnc-EPC1-4, (C) lnc-ZEB2-19, and (D) lnc-GPR89B-15 in the control group compared to hepatitis, cirrhosis, and HCC groups. N VS H, normal versus hepatitis; N VS C, normal versus cirrhosis; N VS L, normal versus HCC; AUC, area under the curve; ROC, receiver operating characteristic; HCC, hepatocellular carcinoma.
Table 2. Characteristics of the patients.
| Group | Variable | No. | Group | Variable | No. |
| Healthy (n=45) | Age, years | Cirrhosis (n = 45) | Age, years | ||
| mean | 29.608 | mean | 51.496 | ||
| SD | 11.331 | SD | 10.339 | ||
| Male | 26 | Male | 34 | ||
| female | 19 | female | 11 | ||
| Serum AFP | Serum AFP | ||||
| <20 ng/ml | 45 | <20 ng/ml | 33 | ||
| ≥20 ng/ml | 0 | =20–400 ng/ml | 10 | ||
| >400 ng/ml | 2 | ||||
| Hepatitis (n=45) | Age, years | liver cancer (n = 45) | Age, years | ||
| mean | 41.143 | mean | 52.333 | ||
| SD | 13.155 | SD | 11.909 | ||
| Male | 38 | Male | 36 | ||
| female | 7 | female | 9 | ||
| Serum AFP | Serum AFP | ||||
| <20 ng/ml | 20 | <20 ng/ml | 16 | ||
| =20–400 ng/ml | 19 | =20–400 ng/ml | 12 | ||
| >400 ng/ml | 6 | >400 ng/ml | 17 |
Table 3. Correlation between expression level of the differentially expressed lncRNAs and clinical characteristics.
| Variable | No. of case | lnc-FAM72D-3 | lnc-EPC1-4 | lnc-ZEB2-19 | lnc-GPR89B-15 | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| Age (year) | MEAN | R | P | MEAN | R | P | MEAN | R | P | MEAN | R | P | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| <35 | 57 | 2.847±0.185 | 0.468 | 0 | 21.842±3.410 | 0.165 | 0.02 | 1.68±0.217 | −0.263 | 0 | 3.720±0.374 | −0.266 | 0 | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| ≥35 | 123 | 12.075±1.076 | 42.430±5.82 | 0.77±0.057 | 2.638±0.148 | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| Gender | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| Male | 134 | 19.114±1.898 | −0.493 | 0 | 40.992±6.441 | −0.05 | 0.649 | 0.970±0.131 | 0.044 | 0.414 | 2.752±0.046 | 0.046 | 0.388 | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| Female | 46 | 4.099±0.307 | 63.036±0.20 | 1.108±0.101 | 32.801±5.274 | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| AFP | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| <20 ng/ml | 114 | −0.069 | 0.197 | 0.16 | 0.003 | −0.037 | 0.493 | −0.118 | 0.027 | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| =20–400 ng/ml | 42 | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| >400 ng/ml | 24 | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| Abbreviations: AFP, alpha fetoprotein; lnc-FAM72D-3, long non-coding RNA FAM72D-3; lnc-EPC1-4, long non-coding RNA; lnc-ZEB2-19, long non-coding RNA ZEB2-19; lnc-GPR89B-15 long non-coding RNA GPR89B-15. | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Function analysis of differentially expressed LncRNAs
The expression levels of lnc-GPR89B-15, lnc-FAM72D-3, lncEPC1-4, and lnc-ZEB2-19 were detected in different cell lines using qPCR (Supplementary Figure 2). lnc-FAM72D-3 and lnc-EPC1-4 expression was significantly different in different cell lines (Supplementary Figure 2). Therefore, we generated siRNAs and overexpression vectors to inhibit and promote lnc-FAM72D-3 and lnc-EPC1-4 expression. On the basis of their high expression levels, we selected HepG2 and SNU-423 cells for evaluation by siRNA-mediated knockdown, while on the basis of their low expression levels, we selected MHCC-LM3 and MHCC97-H cells for vector-mediated overexpression. The lnc-FAM72D-3 siRNA-1 significantly downregulated lnc-FAM72D-3 expression in HepG2 cells with a small but significant decrease in SNU-423 cells, while all three lnc-EPC1-4 siRNAs downregulated lnc-EPC1-4 expression (Supplementary Figure 3, left). Therefore, we selected lnc-FAM72D-3 siRNA-1 and lnc-EPC1-4 siRNA-1 for subsequent analysis (Supplementary Figure 3, right). qPCR results also confirmed that MHCC-LM3 and MHCC97-H overexpression was successful (>150-fold increase). Next, we identified candidate miRNAs predicted bioinformatically to be targeted by these lncRNAs. qPCR results confirmed that hsa-miR-5787 is targeted by lnc-FAM72D (Figure 6A). hsa-miR-29b-1-5p, hsa-miR-19b-1-5p, hsa-miR-511-5p, and SLCO1B1 levels were each significantly downregulated by lnc-EPC1-4 siRNA-1, while CCND1 levels were slightly upregulated (Figure 6B–6D). These results indicated that hsa-miR-5787 might be a target of lnc-FAM72D-3, while lnc-EPC1-4 might target several genes, including hsa-miR-29b-1-5p, hsa-miR-511-5p, and SLCO1B1. To further assess the functional effects of lnc-FAM72D-3 expression modulation, we performed cell proliferation and apoptosis assays. HepG2 cells transfected with lnc-FAM72D-3 siRNA-1 or lnc-EPC1-4 siRNA-1 and MHCC-LM3 cells transfected with pcDNA3.1-FAM72D-3 or pcDNA3.1-EPC1-4 were evaluated by MTS and flow cytometry assays. When HepG2 cells were transfected with lnc-FAM72D-3 siRNA-1 for 48 h, ~50% of HepG2 cells were viable, which was low compared to the siRNA negative control (NC) group, indicating that lnc-FAM72D-3 siRNA-1 significantly inhibits HepG2 cells proliferation (Figure 7A). In addition, pcDNA3.1-FAM72D-3 enhanced HepG2 cell viability. Also, the early apoptotic percentage of HepG2 cells transfected with lnc-FAM72D-3 siRNA-1 was 9.55%, which was low compared to the siRNA NC group, indicating that lnc-FAM72D-3 siRNA-1 induces HepG2 cell apoptosis. In addition, pcDNA3.1-FAM72D-3 induced slight inhibition of cell apoptosis compared to the pcDNA3.1 group (Figure 7B). Lnc-EPC1-4 siRNA and pcDNA3.1-EPC1-4 had the opposite effect: lncEPC1-4 siRNA increased cell viability and decreased apoptosis, while pcDNA3.1-EPC1-4 decreased cell viability and increased apoptosis (Figure 7C, 7D). These results indicated that lnc-FAM72D-3 functions as an oncogene that is upregulated in HCC, while lnc-EPC1-4 functions as a tumor suppressor gene that is downregulated in HCC.
Figure 6. Effect of lnc-FAM72D-3 and lnc-EPC1-4 suppression and overexpression on candidate targets. (A) has-miR-5787 and EIF5 expression when lnc-FAM72D-3 was suppressed or overexpressed in HepG2 and MHCC-LM3 cells. (B) has-miR-29b-1-5p and STAT3 expression when lnc-EPC1-4 was suppressed or overexpressed in HepG2 and MHCC-LM3 cells. (C) has-miR-511-5p and SLCO1B1 expression when lnc-EPC1-4 was suppressed or overexpressed in HepG2 and MHCC-LM3 cells. (D) has-miR-19b-1-5p and CCND1 expression when lnc-EPC1-4 was suppressed or overexpressed in HepG2 and MHCC-LM3 cells. STAT3, signal transducer and activator of transcription 3.
Figure 7. (A) Cell viability and (B) cell apoptosis in cells with lnc-FAM72D-3 suppression or overexpression. (C) Cell viability and (D) cell apoptosis in cells with lnc-EPC1-4 suppression or overexpression.
Discussion
In this study, we examined exosomal lncRNAs that are differentially expressed in hepatic disease patients. Increased lnc-FAM72D-3 and decreased lnc-EPC1-4 levels are clearly correlated with HCC; these findings were verified by qPCR for a group of 180 healthy volunteers and hepatic disease patients. Functional validation with MTS and cell apoptosis assays showed that lnc-FAM72D-3 knockdown induces significant inhibition of cell proliferation and promotes cell apoptosis, while upregulation of lnc-EPC1-4 expression has the same effect. There are currently no reports on the exact underlying mechanism; however, we have elucidated new roles for them in HCC and other hepatic diseases. On the basis of our bioinformatic prediction of downstream targets of lnc-FAM72D-3 and lnc-EPC1-4, several candidate target genes and miRNAs have been shown to regulate cell proliferation and cell apoptosis. These included hsa-miR-511-5p, which we had demonstrated increases by pcDNA3.1-lnc-EPC1-4 transfection. Consistent with our findings, miR-511 acts as a tumor suppressor gene that attenuates tumor cell growth and proliferation by inhibiting tribbles pseudokinase 2 (TRIB2) [25]. TRIB2 is involved in the etiology of various cancers, including leukemia, melanoma, lung cancer, and HCC. In addition, miR-511 regulates the AKT signaling module, which plays a pivotal function in cell proliferation, differentiation, metabolism, and apoptosis in HCC [26–29]. MiR-29b, which we had shown is induced by pcDNA3.1-lnc-EPC1-4 and inhibited by lnc-EPC1-4 siRNA, has also been suggested to participate in the inhibition of angiogenesis and tumorigenesis of cancer through the Akt3 signaling pathway [30]. Therefore, decreased lnc-EPC1-4 siRNA levels in HCC might enhance cell growth by activating the AKT pathway.
We also demonstrated that SLCO1B1 is targeted by lnc-EPC1-4. The sequential administration of sorafenib (front-line treatment) and regorafenib (second-line treatment) efficiently improves therapeutic efficacy in advanced HCC. Influx transporters (SLCO1B1, SLCO1B3, and SLC22A1) are highly involved in the metabolism and clearance of these chemotherapy agents in the liver, while downregulation of these influx transporters might inhibit cellular uptake, thereby discounting their clinical benefits [31–33]. Hu et al. [31] reported that SLCO1B1 expression is significantly downregulated in HCC, which is consistent with our findings. In addition, CDCA7, as a target gene of hsa-miR-29b-1-5p, is overexpressed in several malignancies [34, 35]. lnc-EPC1-4 overexpression induces miR-29b-1-5p expression, which provides a potential mechanism for CDCA7 overexpression and might also explain, in part, the antitumor effects of lnc-EPC1-4.
Our results also demonstrated that lnc-FAM72D-3 siRNA-1 induces hsa-miR-5787 expression, while pcDNA3.1-lnc-FAM72D-3 inhibits hsa-miR-5787 expression. However, modulation of lnc-FAM72D-3 fails to affect eukaryotic translation initiation factor 5 (EIF5), which is a target of has-miR-5787 [36]. EIF5 not only plays an important role in the initiation of translation of some proteins but also regulates the occurrence of cancer and promotes cell proliferation, senescence, and apoptosis. Therefore, a mechanism that involves alternate hsa-miR-5787 targets might account for the oncogenic effects of lnc-FAM72D-3 in HCC cells.
In summary, lnc-FAM72D-3 and lnc-EPC1-4 might be used as diagnostic and prognostic biomarkers for HCC. Future studies that characterize the expression of additional candidate targets of these lncRNAs might help us better define their mechanisms, which might ultimately lead to the development of new diagnostic and therapeutic tools for HCC.
Materials and Methods
Exosome purification from serum
Whole blood was collected from 45 healthy volunteers (controls), 45 hepatitis patients, 45 cirrhosis patients, and 45 HCC patients. Serum was isolated from the whole blood and ultracentrifuged using a Beckman L8 ultracentrifuge (GMI, Ramsey, MN, USA) at 100,000×g for 20 h at 10°C. The supernatants were discarded, and pellets containing exosomes were resuspended in phosphate buffer solution (PBS). The size of the exosomes was validated by using nanoparticle tracking analysis (NTA; Particle Metrix, Germany), and the content of membrane proteins, such as CD63, CD9, and TSG101, was verified using western blotting.
Western blotting
Exosomes, HepG2 cells, and MHCC-LM3 cells were lysed in 1% sodium dodecyl sulfate (SDS) lysis buffer. Protein concentrations were evaluated using the bicinchoninic acid (BCA) protein assay kit (Beyotime, Shanghai, China). Proteins were subjected to 10% sodium dodecyl sulfate–polyacrylamide gel electrophoresis (SDS-PAGE) and then transferred onto nitrocellulose membranes. The nitrocellulose membranes were blocked with 5% nonfat milk in PBS for 1 h at room temperature and incubated overnight at 4°C with primary antibodies anti-CD63 (ab134045; Abcam, UK), anti-TSG101 (ab125011; Abcam), anti-CD9 (ab92726; Abcam), and anti-calnexin (ab22595; Abcam). After several washes, the nitrocellulose membranes were incubated in a blocking buffer with a secondary antibody coupled to horseradish peroxidase (HRP) for 2 h at room temperature. The complexes were detected by electrochemiluminescence (ECL) (Amersham Biosciences/GE Healthcare, Velizy, France).
Cell culture
HCC cell lines (HepG2, Hep3B, SNU-423, and PLC/PRF/5) were obtained from the American Type Culture Collection (ATCC, Manassas, VA, USA). The human normal liver L02 cell line and Huh-7, MHCC-LM3, and MHCC97-H cells were obtained from CELLCOOK (Guangzhou, China) and cultured with 5% CO2 at 37°C. The cells were maintained in high-glucose Dulbecco’s Modified Eagle Medium supplemented with 10% fetal bovine serum (FBS; Gibco, Thermo Fisher Scientific, Waltham, MA, USA) and passaged every 3 days.
LncRNA suppression and overexpression
Full-length lncRNAs were amplified by polymerase chain reaction (PCR) using Thermo Scientific Phusion Flash High-Fidelity PCR Master Mix (Thermo Fisher Scientific), and corresponding complementary DNAs (cDNAs) were subcloned into the NheI and KpnI sites of pcDNA3.1(+) (Genepharma, China) according to the manufacturer’s instructions. The vectors constructed were verified by Sanger sequencing. Small interfering RNAs (siRNAs) for candidate lncRNAs were designed and synthesized. The vectors and siRNAs were transfected into cell lines using lipofectamine 2000 (Thermo Fisher Scientific) according to the manufacturer’s instructions. The expression level of lncRNAs and their targeted miRNAs and messenger RNAs (mRNAs) were determined by real-time quantitative PCR (qPCR) (Supplementary Table 1).
MTS assay
Cells were seeded at a density of 1 × 106/mL into 96-well plates (100 μL/well), treated with overexpression vectors or siRNAs, and examined after 24, 48, and 72 h. A 3-(4,5-dimethylthiazol-2-yl)-5-(3-carboxymethoxyphenyl)-2-(4-sulfophenyl)-2H-tetrazolium (MTS) mixture (Promega Corporation, Madison, WI, USA) was added for 3 h, and the optical density of the cells was detected at a frequency of 490 nm.
Cell apoptosis assay
Cells were resuspended in 1× binding buffer, and 5 μL of fluorochrome-conjugated annexin V and 5 μL of 7-aminoactinomycin D (7-AAD) staining solution (Thermo Fisher Scientific) were added. The cells were assessed using a fluorescence-activated cell sorting Calibur flow cytometer, and the percentage of the apoptotic cells was measured.
RNA-seq
Total RNA was purified from exosomes from the isolated serum from 45 healthy volunteers (controls), 45 hepatitis patients, 45 cirrhosis patients, and 45 HCC patients using the RNeasy Mini Kit (QIAGEN, Germany). RNA integrity was evaluated on the basis of the RNA integrity number (RIN) using Agilent Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA, USA). RNA was further purified using the RNA Clean XP Kit (Beckman Coulter, Brea, CA, UA), and the DNA residue was removed using the RNase-free DNase Set (QIAGEN). RNA quality and concentration were determined using NanoDrop 2000 (Thermo Fisher Scientific). Ribosomal RNA (rRNA) was removed using the NEBNext rRNA Depletion kit (New England Biolabs, Ipswich, MA, USA). Then, 1 μg of total RNA was used for library preparation using the VAHTSTM mRNA-seq v2 library Prep Kit (Vazyme, Nanjing, China) according to the manufacturer’s instructions. RNA was fragmented and then double-stranded cDNA was synthesized. End-polishing was performed, and the cDNA fragments were ligated with adapters and then subjected to universal PCR amplification in order to obtain a sufficient library for sequencing. The library quality was assessed using Agilent Bioanalyzer 2100. RNA-seq was performed using Illumina Hiseq 4000 (Illumina, San Diego, CA, USA). The data were assembled and annotated with corresponding symbols for transcripts. Differentially expressed lncRNAs were screened using R software (q < 0.05 and fold change > 2.0).
Gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis
The Database for Annotation, Visualization and Integrated Discovery was used to annotate the potential functions of various signaling pathways for differentially expressed lncRNAs. The functional annotations of parental genes were predicted by GO functional annotation. Scatter plots were used to visualize GO analysis results. KEGG pathway annotation was used to identify relevant pathways of differentially expressed lncRNAs.
Construction of an lncRNA-miRNA-mRNA network
We predicted miRNA binding to lncRNA using miRanda software. By inputting miRNA and lncRNA sequences, we predicted the miRNA that lncRNA might combine with and filtered out the miRNA that was differentially expressed in RNA-seq. miRNA target genes were downloaded in miRTarbase, while an lncRNA-miRNA-mRNA regulatory network was constructed using Cytoscape ver. 3.5.1.
qPCR
Total RNA was extracted from exosomes and gene expression verified by PCR. cDNA was synthesized from total RNA using M-MLV Reverse Transcriptase (Promega Corporation) according to the manufacturer’s instructions. PCR was performed using the GoTaq qPCR Master Mix (Promega Corporation). PCR amplification was performed using the ABI 7500 system (Applied Biosystems, Foster City, CA, USA). Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) was used as an internal control.
Statistical analysis
Statistical analysis was performed using SPSS Statistics 19 (IBM Corporation. Armonk, NY, USA) and GraphPad Prism (GraphPad Software, CA, USA). Comparisons were performed using independent t-tests, one-way analysis of variance (ANOVA), and two-way ANOVA when data corresponded to a normal distribution. P < 0.05 was considered statistically significant.
Ethics statement
This study was approved by the ethics committee of Third Affiliated Hospital of Sun Yat-sen University. The patients provided informed written consent.
Supplementary Materials
Supplementary Figures
Supplementary Table 1
Supplementary Table 2
Supplementary Table 3
Supplementary Table 4
Author Contributions
Conceptualization and Methodology: Qi Zhang, Meihai Deng, and Zhicheng Yao; Validation: Changchang Jia, Yan Tai, and Hao Liang. Formal analysis: Qi Zhang, Zhicheng Yao, and Yan Tai. Investigation: Qi Zhang, Zhicheng Yao, Yan Tai, and Zhaozhong Zhong. Writing-original draft: Qi Zhang and Zhicheng Yao.
Conflicts of Interest
These authors declare no conflicts of interst.
Funding
The present study was supported by the National Natural Science Foundation of China (grant no. 81870128), the Science and Technology Planning Project of Guangdong (2014A020212185), the Natural Science Foundation of Guangdong Province (2018A030313641), Young teacher cultivation project of Sun Yat-sen university (19ykpy34).
References
- 1. GBD 2013 Mortality and Causes of Death Collaborators. Global, regional, and national age-sex specific all-cause and cause-specific mortality for 240 causes of death, 1990-2013: a systematic analysis for the global burden of disease study 2013. Lancet. 2015; 385:117–71. https://doi.org/10.1016/S0140-6736(14)61682-2 [PubMed]
- 2. Siegel RL, Miller KD, Jemal A. Cancer statistics, 2017. CA Cancer J Clin. 2017; 67:7–30. [PubMed]
- 3. Kumar N, Choudhary NS, Saraf N, Saigal S. Hepatocellular carcinoma presenting as dyspnea at rest due to innumerable lung metastasis. J Clin Exp Hepatol. 2016; 6:319–20. https://doi.org/10.1016/j.jceh.2016.05.005 [PubMed]
- 4. Zhou L, Liu J, Luo F. Serum tumor markers for detection of hepatocellular carcinoma. World J Gastroenterol. 2006; 12:1175–81. https://doi.org/10.3748/wjg.v12.i8.1175 [PubMed]
- 5. Yang L, Lin C, Jin C, Yang JC, Tanasa B, Li W, Merkurjev D, Ohgi KA, Meng D, Zhang J, Evans CP, Rosenfeld MG. lncRNA-dependent mechanisms of androgen-receptor-regulated gene activation programs. Nature. 2013; 500:598–602. https://doi.org/10.1038/nature12451 [PubMed]
- 6. Huarte M. The emerging role of lncRNAs in cancer. Nat Med. 2015; 21:1253–61. https://doi.org/10.1038/nm.3981 [PubMed]
- 7. Yan X, Zhang D, Wu W, Wu S, Qian J, Hao Y, Yan F, Zhu P, Wu J, Huang G, Huang Y, Luo J, Liu X, et al. Mesenchymal stem cells promote hepatocarcinogenesis via lncRNA-MUF interaction with ANXA2 and miR-34a. Cancer Res. 2017; 77:6704–16. https://doi.org/10.1158/0008-5472.CAN-17-1915 [PubMed]
- 8. Yin J, Luo W, Zeng X, Zeng L, Li Z, Deng X, Tan X, Hu W. UXT-AS1-induced alternative splicing of UXT is associated with tumor progression in colorectal cancer. Am J Cancer Res. 2017; 7:462–72. [PubMed]
- 9. Martens-Uzunova ES, Böttcher R, Croce CM, Jenster G, Visakorpi T, Calin GA. Long noncoding RNA in prostate, bladder, and kidney cancer. Eur Urol. 2014; 65:1140–51. https://doi.org/10.1016/j.eururo.2013.12.003 [PubMed]
- 10. Ellinger J, Alam J, Rothenburg J, Deng M, Schmidt D, Syring I, Miersch H, Perner S, Müller SC. The long non-coding RNA lnc-ZNF180-2 is a prognostic biomarker in patients with clear cell renal cell carcinoma. Am J Cancer Res. 2015; 5:2799–807. [PubMed]
- 11. Shi J, Li X, Zhang F, Zhang C, Guan Q, Cao X, Zhu W, Zhang X, Cheng Y, Ou K, Chen Q, Hu S. Circulating lncRNAs associated with occurrence of colorectal cancer progression. Am J Cancer Res. 2015; 5:2258–65. [PubMed]
- 12. Zhang K, Shi H, Xi H, Wu X, Cui J, Gao Y, Liang W, Hu C, Liu Y, Li J, Wang N, Wei B, Chen L. Genome-wide lncRNA microarray profiling identifies novel circulating lncRNAs for detection of gastric cancer. Theranostics. 2017; 7:213–27. https://doi.org/10.7150/thno.16044 [PubMed]
- 13. Deng L, Yang SB, Xu FF, Zhang JH. Long noncoding RNA CCAT1 promotes hepatocellular carcinoma progression by functioning as let-7 sponge. J Exp Clin Cancer Res. 2015; 34:18. https://doi.org/10.1186/s13046-015-0136-7 [PubMed]
- 14. Jing W, Gao S, Zhu M, Luo P, Jing X, Chai H, Tu J. Potential diagnostic value of lncRNA SPRY4-IT1 in hepatocellular carcinoma. Oncol Rep. 2016; 36:1085–92. https://doi.org/10.3892/or.2016.4859 [PubMed]
- 15. Zhang H, Xing Z, Mani SK, Bancel B, Durantel D, Zoulim F, Tran EJ, Merle P, Andrisani O. RNA helicase DEAD box protein 5 regulates polycomb repressive complex 2/hox transcript antisense intergenic RNA function in hepatitis B virus infection and hepatocarcinogenesis. Hepatology. 2016; 64:1033–48. [PubMed]
- 16. Chen F, Bai G, Li Y, Feng Y, Wang L. A positive feedback loop of long noncoding RNA CCAT2 and FOXM1 promotes hepatocellular carcinoma growth. Am J Cancer Res. 2017; 7:1423–34. [PubMed]
- 17. Gong J, Qi X, Zhang Y, Yu Y, Lin X, Li H, Hu Y. Long noncoding RNA linc00462 promotes hepatocellular carcinoma progression. Biomed Pharmacother. 2017; 93:40–47. https://doi.org/10.1016/j.biopha.2017.06.004 [PubMed]
- 18. Hou Z, Xu X, Fu X, Tao S, Zhou J, Liu S, Tan D. HBx-related long non-coding RNA MALAT1 promotes cell metastasis via up-regulating LTBP3 in hepatocellular carcinoma. Am J Cancer Res. 2017; 7:845–56. [PubMed]
- 19. Xu LC, Chen QN, Liu XQ, Wang XM, Chang QM, Pan Q, Wang L, Wang YL. Erratum: up-regulation of LINC00161 correlates with tumor migration and invasion and poor prognosis of patients with hepatocellular carcinoma. Oncotarget. 2019; 10:1474. https://doi.org/10.18632/oncotarget.26712 [PubMed]
- 20. van Niel G, Porto-Carreiro I, Simoes S, Raposo G. Exosomes: a common pathway for a specialized function. J Biochem. 2006; 140:13–21. https://doi.org/10.1093/jb/mvj128 [PubMed]
- 21. Mathivanan S, Ji H, Simpson RJ. Exosomes: extracellular organelles important in intercellular communication. J Proteomics. 2010; 73:1907–20. https://doi.org/10.1016/j.jprot.2010.06.006 [PubMed]
- 22. Sugimachi K, Matsumura T, Hirata H, Uchi R, Ueda M, Ueo H, Shinden Y, Iguchi T, Eguchi H, Shirabe K, Ochiya T, Maehara Y, Mimori K. Identification of a bona fide microRNA biomarker in serum exosomes that predicts hepatocellular carcinoma recurrence after liver transplantation. Br J Cancer. 2015; 112:532–38. https://doi.org/10.1038/bjc.2014.621 [PubMed]
- 23. Moris D, Beal EW, Chakedis J, Burkhart RA, Schmidt C, Dillhoff M, Zhang X, Theocharis S, Pawlik TM. Role of exosomes in treatment of hepatocellular carcinoma. Surg Oncol. 2017; 26:219–28. https://doi.org/10.1016/j.suronc.2017.04.005 [PubMed]
- 24. Zhang Z, Li X, Sun W, Yue S, Yang J, Li J, Ma B, Wang J, Yang X, Pu M, Ruan B, Zhao G, Huang Q, et al. Loss of exosomal miR-320a from cancer-associated fibroblasts contributes to HCC proliferation and metastasis. Cancer Lett. 2017; 397:33–42. https://doi.org/10.1016/j.canlet.2017.03.004 [PubMed]
- 25. Zhang C, Chi YL, Wang PY, Wang YQ, Zhang YX, Deng J, Lv CJ, Xie SY. miR-511 and miR-1297 inhibit human lung adenocarcinoma cell proliferation by targeting oncogene TRIB2. PLoS One. 2012; 7:e46090. https://doi.org/10.1371/journal.pone.0046090 [PubMed]
- 26. Leng J, Han C, Demetris AJ, Michalopoulos GK, Wu T. Cyclooxygenase-2 promotes hepatocellular carcinoma cell growth through akt activation: evidence for akt inhibition in celecoxib-induced apoptosis. Hepatology. 2003; 38:756–68. [PubMed]
- 27. Nakanishi K, Sakamoto M, Yamasaki S, Todo S, Hirohashi S. Akt phosphorylation is a risk factor for early disease recurrence and poor prognosis in hepatocellular carcinoma. Cancer. 2005; 103:307–12. [PubMed]
- 28. Schmitz KJ, Wohlschlaeger J, Lang H, Sotiropoulos GC, Malago M, Steveling K, Reis H, Cicinnati VR, Schmid KW, Baba HA. Activation of the ERK and AKT signalling pathway predicts poor prognosis in hepatocellular carcinoma and ERK activation in cancer tissue is associated with hepatitis C virus infection. J Hepatol. 2008; 48:83–90. https://doi.org/10.1016/j.jhep.2007.08.018 [PubMed]
- 29. Nakajima T, Yasui K, Zen K, Inagaki Y, Fujii H, Minami M, Tanaka S, Taniwaki M, Itoh Y, Arii S, Inazawa J, Okanoue T. Activation of b-myb by E2F1 in hepatocellular carcinoma. Hepatol Res. 2008; 38:886–95. [PubMed]
- 30. Teng Y, Zhang Y, Qu K, Yang X, Fu J, Chen W, Li X. MicroRNA-29B (mir-29b) regulates the warburg effect in ovarian cancer by targeting AKT2 and AKT3. Oncotarget. 2015; 6:40799–814. https://doi.org/10.18632/oncotarget.5695 [PubMed]
- 31. Hu DG, Marri S, McKinnon RA, Mackenzie PI, Meech R. Deregulation of the genes that are involved in drug absorption, distribution, metabolism, and excretion in hepatocellular carcinoma. J Pharmacol Exp Ther. 2019; 368:363–81. https://doi.org/10.1124/jpet.118.255018 [PubMed]
- 32. Sekine S, Ogawa R, Ojima H, Kanai Y. Expression of SLCO1B3 is associated with intratumoral cholestasis and CTNNB1 mutations in hepatocellular carcinoma. Cancer Sci. 2011; 102:1742–47. [PubMed]
- 33. Herraez E, Lozano E, Macias RI, Vaquero J, Bujanda L, Banales JM, Marin JJ, Briz O. Expression of SLC22A1 variants may affect the response of hepatocellular carcinoma and cholangiocarcinoma to sorafenib. Hepatology. 2013; 58:1065–73. [PubMed]
- 34. Martín-Cortázar C, Chiodo Y, Jiménez RP, Bernabé M, Cayuela ML, Iglesias T, Campanero MR. CDCA7 finely tunes cytoskeleton dynamics to promote lymphoma migration and invasion. Haematologica. 2020; 105:730–40. https://doi.org/10.3324/haematol.2018.215459 [PubMed]
- 35. Ye L, Li F, Song Y, Yu D, Xiong Z, Li Y, Shi T, Yuan Z, Lin C, Wu X, Ren L, Li X, Song L. Overexpression of CDCA7 predicts poor prognosis and induces EZH2-mediated progression of triple-negative breast cancer. Int J Cancer. 2018; 143:2602–13. https://doi.org/10.1002/ijc.31766 [PubMed]
- 36. Yoo H, Yoo JK, Lee J, Lee DR, Ko JJ, Oh SH, Choo YK, Kim JK. The hsa-miR-5787 represses cellular growth by targeting eukaryotic translation initiation factor 5 (eIF5) in fibroblasts. Biochem Biophys Res Commun. 2011; 415:567–72. https://doi.org/10.1016/j.bbrc.2011.10.103 [PubMed]