- 1First People’s Hospital of Hangzhou Lin’an District, Affiliated Lin’an People’s Hospital, Hangzhou Medical College, Hangzhou, China
- 2Department of Research and Development, Zhejiang Zhongwei Medical Research Center, Hangzhou, China
- 3Department of Obstetrics and Gynecology, The Second Affiliated Hospital of Wenzhou Medical University, Wenzhou, China
Background: Hepatocellular carcinoma (HCC) is a high-burden cancer. The molecular mechanism of HCC has not been fully elucidated. Notably, current research has revealed a significant function for long non-coding RNAs (lncRNAs) in the prognosis of patients with HCC. Here, this study aims to construct a regulated lncRNA-mediated ceRNA network and find biological targets for the treatment of HCC.
Methods: Based on the RNA expression patterns from the TCGA, we did an analysis to determine which genes were expressed differently between liver tumor tissues and noncancerous tissues. Then, using bioinformatic tools, we built a lncRNA-miRNA-mRNA ceRNA network and used GO and KEGG functional analyses on the DEmRNAs connected to ceRNA networks. The main lncRNAs in the subnetwork were chosen, and we next looked at the relationships between these lncRNAs and the clinical characteristics of patients with HCC. The prognosis-related genes and immune cells were identified using Kaplan-Meier and Cox proportional hazard analyses, and CIBERSORT was utilized to separate the 22 immune cell types. CCK8 assay was performed to measure cell viability in HCC cells after lncRNA HOTTIP modulation.
Results: Differentially expressed mRNA and lncRNAs in HCC and paracancerous tissues were identified. There are 245 lncRNAs, 126 miRNAs, and 1980 mRNAs that are expressed differently in liver tumour tissues than in noncancerous cells. Function analysis showed that mRNAs in ceRNA network were significantly enriched in G1/S transition of mototiv cell cycle, positive regulation of cell cycle process, hepatocellular carcinoma, and cancer related pathways. CD8 T cells and T follicular helper cells had a favourable link with a 0.65 correlation coefficient. Additionally, there was a strong correlation between Eosinophils, activated NK cells, and B memory cells. Strikingly, depletion of lncRNA HOTTIP inhibited viability of HCC cells. In addition, miR-205 upregulation suppressed viability of HCC cells, while miR-205 downregulation repressed viability of HCC cells. Notably, miR-205 depletion rescued HOTTIP depletion-mediated suppression of cell viability in HCC.
Conclusion: A ceRNA network was created by examining the lncRNA, miRNA, and mRNA expression profiles of liver tumours from the TCGA database. LncRNA HOTTIP promoted cell viability via inhibition of miR-205 in HCC cells.
Introduction
Hepatocellular carcinoma (HCC) is one of the most common tumors in the world (1). The epidemiology of HCC shows two main patterns: one in North America and Western Europe, and the other in non-Western countries, such as sub-Saharan Africa, Central and South-east Asia, and the Amazon Basin (2). Men are generally more affected than women, and most commonly between the ages of 30 and 50, and HCC causes 662,000 deaths worldwide each year, about half of them in China (3). By 2022, there will be 41,260 new cases and 30,520 new deaths according to cancer statistics (4).
The treatment of HCC can be divided into surgical treatment and non-surgical treatment (5). Surgical treatment includes resection, cryoablation and liver transplantation. Nonsurgical treatments include either liver-specific therapy (i.e., percutaneous ethanol injection, radiofrequency/microwave ablation, trans-arterial embolization, external radiation therapy) or systemic therapy (chemotherapy, molecular targeted therapy, immune checkpoint inhibitor therapy) (6). Surgical treatment is the main treatment for HCC, but the 5-year survival rate after surgical resection is only about 30% (7). Currently, the molecular mechanism of HCC occurrence has not been fully elucidated (8), so it is of great significance to explore the molecular mechanism of HCC occurrence and development and screen biological markers for the early diagnosis of HCC for the treatment of this deadly disease.
Long non-coding RNA (lncRNA) is a class of non-coding RNA with a length greater than 200bp, which has no protein-coding function, poor conserved between species, and strong tissue specificity and spatio-temporal specificity (9). Recently, more and more lncRNAs have been found to act as oncogenic or tumor suppressor factors, and many oncogenic genomic loci are mainly transcribed from those lncRNAs that play an important role in cancer induction (10–13). With the rapid development of cancer research, a large number of lncRNAs related to tumorigenesis have been identified (14–16). Studies have shown that lncRNA MEG3 can promote the proliferation of HCC cells (17, 18). At the same time, some studies have proved that lncRNAs HOTTIP and UCA1 can promote the deterioration of HCC. Epigenetic related molecules (such as SNRPC), noncoding RNAs (such as HSA-Mir-221) and immune-related molecules (such as DCK) can be used as potential biomarkers for diagnosis, treatment and prognosis of HCC (19–22).
Studies have shown that lncRNAs have a wide range of biological functions, such as participating in RNA generation and processing, transcription regulation and chromatin remodeling (23). However, the mechanism of action is complex, such as participating in transcriptional regulation by binding characteristic proteins or taking part in post-transcriptional regulation as ceRNA (24). According to the ceRNA hypothesis proposed in 2011, when circRNA and mRNA share the same miRNA response element (MRE), they will competitively bind MRE to regulate the expression level of related genes. The most common mechanism of action between lncRNAs and miRNAs is that lncRNAs directly complement and pair with miRNA seed sequences, adsorb miRNAs to form complexes, and reduce the number of miRNAs involved in downstream gene regulation in cells (25). Studies have shown that the abnormal expression of lncRNAs affects the expression of mRNAs through sponge miRNAs, leading to tumorigenesis and cancer progression (26, 27). Therefore, in order to better study the functions of lncRNAs, there is an urgent need to better understand the regulatory networks of lncRNAs. This article aimed to perform a comprehensive analysis of TCGA transcriptome data to identify differentially expressed lncRNAs, miRNAs, and mRNAs to generate a regulated ceRNA network.
Methods:
Database
Data were downloaded from the TCGA database (https://portal.gdc.cancer.gov/) (28). The data were used to study the differential expression of lncRNA, mRNA, miRNA among HCC patients. In addition, edgeR(3.28.0) (29) of R language was mainly used to identify differentially expressed RNA. Screening criteria is: |logFC| > 2, and P value < 0.01 for mRNA and miRNA, and |logFC| > 4, and P value < 0.01 for lncRNA.
Construction of ceRNA network and prognosis model
Using the mircode database, the correlation analysis of differentially expressed lncRNAs and mRNAs was carried out. Then, a ceRNA network was built using the miRDB, miRTarBase, and TargetScan databases. By adding edge and nodal gene information to the interaction network with Cytoscape 3.8.0, the HCC ceRNA network was displayed. The Kyoto Encyclopedia of Genes and Genomes and Gene Ontology (GO) (30) were then used to annotate the differentially expressed mRNAs in the ceRNA network. The R language’s “clusterProfiler,” (31) “org.Hs.eg.db,” “enrichplot,” and “ggplot2” (32) packages were used to carry out the aforementioned analysis. According to the median model risk score, all patients were separated into high-risk and low-risk groups, and Kaplan-Meier survival curves were created to analyze the variations in survival between the two groups (33). It can be regarded as statistically significant when P <0.05. Using the area under the curve (AUC), the model’s predictive effectiveness was assessed (visualized by the “timeROC” package in R language). After that, the “rms” package of R was used to build the nomogram survival prediction map, and the calibration curve was used to assess the accuracy of the map in predicting the survival rate of HCC patients.
Assessing immune cells
Using the “CIBERSORT” software tool in the R language, the abundance of 22 immune cell subsets in 424 tumour specimens were estimated (34). The threshold used was P <0.05. The “Pheatmap” programme was used in the study to show immune cell infiltration in HCC. To examine the immune cells that are connected to survival, we employed univariate Cox regression (35). Lasso regression and multivariate Cox regression analysis were used together to create the prognosis related model (36), and the model’s accuracy was assessed using the AUC value. We estimated the prognostic model’s quantification using the multi-index AUC value, and we determined the nomogram’s predictive power using the calibration curve. To investigate the connection between important genes in the ceRNA network and immune cells in the model, we employed Pearson correlation analysis (37).
Construction of nomogram model
For the derived independent prognostic indicators, a nomogram model was created and shown. The concordance index (C-index) of the independent prognostic factors and the Nomogram model, comprised of risk group fitting with the Coxph model, were determined to validate the nomogram’s prediction abilities. The statistical test was then run with resampling, and the significant P < value was determined. The degree to which the independent prognostic factors and compound factors fitted the Coxph model was evaluated, and the factor with the lowest P value was chosen to draw the calibration curve. The model created by this component performed better when the calibration curve was closer to 45°.
Transfection
The HOTTIP siRNA and negative control siRNA (NC), miR-205 mimics and miR-205 inhibitors were obtained from GenePharm (Shanghai, China). The HOTTIP and miR-205 were transfected into HCC cells by lipofectamine 3000 (Invitrogen) based on the manufacturer’s protocol.
Reverse transcription-quantitative PCR analysis
Total RNA was isolated from HCC cell lines using TRIzol reagent based on the manufacturer’s protocol. RT-PCR was performed as described before. The primers are as follows: HOTTIP: forward primer: 5’-CCT AAA GCC ACG CTT CTT TG-3’; reverse primer: 5’-TGC AGG CTG GAG ATC CTA CT-3’; GAPDH: forward primer: 5’-GTC AAC GGA TTT GGT CTG TATT-3’; reverse primer: 5’-AGT CTT CTG GGT GGC AGT GAT-3’.
Cell counting assay Kit-8 assay
Huh-7 and HepG2 HCC cell lines were cultured in Dulbecco’s minimal essential medium (DMEM) supplemented with 5% fetal bovine serum, 100 Units/ml penicillin and 100 μg/ml streptomycin, at 37°C under a 5% CO2 atmosphere. CCK-8 assay was performed to analyze cell proliferation based on the manufacture’s protocol in HCC cell lines after 48 and 72 h transfection. The OD450nm was measured by a microplate reader.
Statistical analysis
Statistical analyses for CCK-8 assay were performed by GraphPad Prism 5.0 (Graph Pad Software, La Jolla, CA). The results were presented as mean ± standard deviation (SD). Multiple groups were analyzed by one-way analysis of variance (ANOVA) test followed Turkey’s posttest. P < 0.05 was indicated statistically significant.
Results:
Expression of DERNAs
In order to better understand the expression of RNAs in HCC, we conducted differential analysis on mRNAs, lncRNAs and miRNAs, and the results are shown in Figures 1A–F. The Volcano diagram and Heatmap demonstrated that 1980 DEmRNAs had 1770 up-regulated and 210 down-regulated mRNAs (Figures 1A, B). There were 245 DElncRNAs, including 242 up-regulated and 3 down-regulated lncRNAs (Figures 1C, D). Moreover, 126 DEmiRNAs had 123 up-regulated and 3 down-regulated miRNAs (Figures 1E, F).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1042928/fonc-12-1042928-HTML-r1/image_m/fonc-12-1042928-g001.jpg)
Figure 1 Volcano plots and heatmaps of differentially expressed mRNA, miRNA, and lncRNA. (A) Volcano plot of differentially expressed mRNA. (B) Heatmap of differentially expressed mRNA. (C) Volcano plot of differentially expressed lncRNA. (D) Heatmap of differentially expressed lncRNA. (E) Volcano plot of differentially expressed miRNA. (F) Heatmap of differentially expressed miRNA.
Go and KEGG pathway analysis
The biofunctional analysis of the ceRNA co-modules is undertaken to further investigate the functional role of mRNAs/miRNAs/lncRNAs in carcinogenesis and if the ceRNA co-modules identified in the study are connected with HCC. Function analysis showed that mRNAs in ceRNA network were significantly enriched in G1/S transition of mototiv cell cycle, positive regulation of cell cycle process, hepatocellular carcinoma, and cancer related pathways (Figures 2A, B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1042928/fonc-12-1042928-HTML-r1/image_m/fonc-12-1042928-g002.jpg)
Figure 2 lncRNA-miRNA-mRNA ceRNA network. (A) GO analysis (B) KEGG Pathway Analysis (C) lncRNA-miRNA-mRNA ceRNA network generated by Cytoscape (version 3.6.1, https://cytoscape.org/). Red: upregulated RNAs. Green: downregulated RNAs. Triangle: miRNA. Diamond: lncRNA. Circular: mRNA.
Construction of a lncRNA-miRNA-mRNA ceRNA network
Based on the several databases, including miRcode database, a ceRNA network was constructed by differentially expressed lncRNAs, miRNAs and mRNAs. 8 DElncRNAs, 11 DEmiRNAs, 98 DEmRNAs are shown as nodes, and 156 interactions are shown as edges. Visualization is achieved by Cytoscape (Figure 2C).
Prognosis model of hepatocellular carcinoma
Our data imply that 24 mRNAs may be protective factors (Figure 3). We screened 24 genes associated with HCC prognosis (Figure 3). The prognosis of 31 genes was determined by univariate Cox regression analysis. Then, the results of Lasso regression showed that ERVMER61-1, SLC7A11, CEP55, CBX2, EZH2, PBK, AP002478.1 and miR-137 were associated with the HCC prognosis (P <0.05) (Figure 4A). The median risk score was used to categorise patients into high-risk and low-risk groups. The high-risk samples had a lower survival duration than the low-risk samples (Figures 4B–D). The AUC at 1 years, 3 years, 5 years were 0.764,0.61 and 0.802 (Figure 4E). The low-risk group had a worse chance of survival than the high-risk group (Figure 4F).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1042928/fonc-12-1042928-HTML-r1/image_m/fonc-12-1042928-g004.jpg)
Figure 4 LASSO and multivariate regression models screening 4 genes. (A) different characteristics of genes and their corresponding coefficients. (B) cross-validation to build the model. (C) risk score with survival status cases. (D) risk score with survival time. (E) ROC curve plotted under the prediction model. (F) KM curve of this prediction model.
The independent prognostic indicators were then used to create the Nomogram, which forecasts survival rates for the first, three to five years. There is a nomogram model made up of CBX2, PBK and AP002478.1 (Figures 5A, B). The nomogram model’s calibration curves demonstrated good agreement between the anticipated 1/3/5-year overall survival rate and the actual survival rate, indicating that the model was accurate in its predictions (Figure 5C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1042928/fonc-12-1042928-HTML-r1/image_m/fonc-12-1042928-g005.jpg)
Figure 5 The association of CBX2, PBK, AP002478.1 and survival. (A) Heatmap of differentially expressed genes. (B) Nomogram based on multiple Cox regression. (C) Calibration curve for 3-year survival. The closer to the diagonal, the better the prediction effect. (D) Co-expression patterns among fractions of immune cells.
Immune infiltration in HCC
CD8 T cells were positively correlated with T follicular helper cells, and the correlation coefficient was 0.65. Additionally, there was a positive link between Eosinophils and NK activated cells and B memory cells with a correlation coefficient of 0.75 and 0.68, respectively (Figure 5D). The expression of dendritic resting cells was higher in G3/4 and lower in G1/2 (Figure 6A). Similarly, T cell expression was higher in the G3/4 than in the G1/2 (Figure 6A). Neutrophils was highly expressed in T3/4 (Figure 6A). Macrophage M2 was highly expressed in T3/4 and stage III/IV (Figure 6B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1042928/fonc-12-1042928-HTML-r1/image_m/fonc-12-1042928-g006.jpg)
Figure 6 Immune cell infiltration in HCC. (A, B) Correlation between immune cells and clinical shape. (C, D): different characteristics of immune cells and their corresponding coefficients; (E) ROC curve plotted under the prediction model). (F) KM curve of this prediction model.
To prevent overfitting and resolve severe collinearity, we performed lasso regression (Figures 6C, D). The AUC curve demonstrates that the nomogram survival prediction model had a high accuracy (Figure 6E). The AUC at 1 year, 3 years, and 5 years were 0.764, 0.61, and 0.812, respectively (Figure 6E). The high-risk group had worse survival than the low-risk group (Figure 6F). Figure 7A depicted that dendritic activated cells, and eosinophil were highly expressed in tumor cells. The Pearson correlation test was used to investigate the relationship between major genes in the ceRNA network and prognosis-related immune cells (Figure 7B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1042928/fonc-12-1042928-HTML-r1/image_m/fonc-12-1042928-g007.jpg)
Figure 7 Immune cells and CBX2, PBK, AP002478.1 in HCC. (A) Heatmap of the three immune cells in the Cox regression model. (B) correlation between immune cells and genes.
LncRNA HOTTIP promotes cell viability via targeting miR-205
Our ceRNA network showed the role of lncRNA HOTTIP and miR-205 in HCC progression. To validate the function of lncRNA HOTTIP on cell viability of HCC, HOTTIP was downregulated by siRNA transfection in HCC cells (Figure 8A). CCK-8 data showed that downregulation of HOTTIP led to inhibition of cell viability in HepG2 and Huh7 cells (Figure 8B). Moreover, miR-205 expression level was increased in HCC cells after miR-205 mimic transfection, while miR-205 was downregulated in HCC cells after miR-205 inhibitor treatment (Figure 8C). Furthermore, inhibition of miR-205 increased viability of HCC cells, whereas upregulation of miR-205 decreased viability of HCC cells (Figures 8D, E). Notably, inhibition of miR-205 rescued HOTTIP siRNA-mediated suppression of cell viability in HCC (Figure 8F). Altogether, TOTTIP enhanced cell viability via targeting miR-205 in HCC cells.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1042928/fonc-12-1042928-HTML-r1/image_m/fonc-12-1042928-g008.jpg)
Figure 8 LncRNA HOTTIP promoted cell viability via miR-205 in HCC. (A) RT-PCR was done to measure the expression of HOTTIP in HCC cells after HOTTIP siRNA transfection. (B) CCK8 assay was done to measure cell viability in HCC cells after HOTTIP siRNA transfection. (C): RT-PCR was done to measure the expression of miR-205 in HCC cells after miR-205 inhibitor or mimics transfection. (D, E) CCK8 assay was done to measure cell viability in HCC cells after miR-205 inhibitor or mimics transfection. (F) CCK8 assay was done to measure cell viability in HCC cells after HOTTIP siRNA transfection or miR-205 inhibitor or combination treatments. NC, negative control; inh, inhibitor; CTR, control. *P < 0.05; **P < 0.01 compared with NC or CTR; &P<0.05 compared with siRNA alone or inhibitor alone.
Discussion
HCC is one of the most common malignant tumors in the gastrointestinal system (38). At present, the molecular pathogenesis of HCC is still unclear (39). Therefore, it is very important to understand the molecular mechanism and process of HCC and identify new therapeutic targets to improve the clinical prognosis of patients. LncRNAs have been considered as an underappreciated novel therapeutic target (40). LncRNAs play important roles in gene regulation, including regulation of gene activation and silencing, X chromosome inactivation, alternative splicing, and post-translational regulation (41). Besides, lncRNAs have also been shown to be associated with tumor progression (12). Therefore, in this study, we studied the relationship between lncRNAs and liver cancer, and constructed a ceRNA regulatory network based on lncRNAs.
The current work used lncRNA-associated ceRNA to identify key biomarkers related to the prognosis of HCC. The lncRNAs LINC00491, DSCR4, AC061975.6, LINC00221, TDRG1 and HOTTIP served as hub nodes in the ceRNA network, which targeted other miRNAs and mRNAs. The ceRNA modulatory network of CBX2, PBK and AP002478.1 may impact HCC progression, according to KM survival and correlation study. Consistent with our study, one recent study suggested that LINC00491 promotes tumor growth and lung metastasis in mouse xenografts. LINC00491 is highly expressed in HCC patients and is associated with poor prognosis. More deeply, LINC00491 promotes HCC progression by sponging miR-324-5p/ROCK1 and may be a potential therapeutic target for HCC (42). In addition, overexpression of lncRNA TDRG1 promotes the viability, invasion and migration of endometrial cancer cells, inhibits cell apoptosis, and upregulates the expression of VEGF-A, PI3K, Bcl-2, MMP2 and surviving (43). Similarly, TDRG1 overexpression also promoted GC growth and metastasis in vitro and in vivo, which was regulated by the miR-873-5p/HDGF pathway (44).
LncRNA HOTTIP has been reported to regulate HCC development and progression (45). In our study, we found that HOTTIP can act as a node of ceRNA and is closely related to the development of HCC. Meanwhile, a study also demonstrated that HOTTIP acts as a molecular sponge for miR-148a-3p to increase WNT1 expression, thereby modulating the CSC-like properties of breast cancer, suggesting that HOTTIP is a new target for breast cancer treatment (46). LncRNA HOTTIP expression was correlated with tumor recurrence in HCC patients following liver transplantation (47). Similarly, lncRNA HOTTIP/HOXA13 expression was linked to tumor progression and therapeutic outcome in HCC patients (19). Additionally, lncRNA HOTTIP was elevated in HCC and was regulated by miR-125b (48). Ge et al. found that miR-192 and miR-204 reduced the expression of lncRNA HOTTIP and impaired GLS1-involved glutaminolysis in HCC (49). Tang et al. reported that solamargine inhibited the expression of lncRNA HOTTIP and TUG1, increased miR-4726-5p expression and inhibited MUC1 expression in HCC cells, leading to inhibition of cell growth and promotion of sorafenib efficacy (50). Our in vitro data further identified that HOTTIP promoted viability of HCC cells via regulation of miR-205. Dong et al. discovered that lncRNA HOTTIP promoted cisplatin resistance via targeting miR-205 and modulating ZEB2 in ovarian cancer cells (51). It is pivotal to determine whether lncRNA HOTTIP enhanced viability of HCC cells via regulation of miR-205/ZEB2 axis.
In our study, we found that dendritic activated cells and eosinophil are closely related to liver cancer. Dendritic cells are the most powerful professional antigen-presenting cells, widely distributed in peripheral lymphoid tissues, and play an important role in the immune regulation of the body. Liver DCs play an important role in the regulation of hepatic immune responses, and their precursor cells can induce antigen-specific immune tolerance in T cells that have not been stimulated by specific antigens, while for T memory cells, they are strong immune stimulators (52–54). In vitro studies have shown that once hepatic DC precursor cells are induced by antigen, they can also stimulate the proliferation of non-antigen-sensitized allogeneic T cells (55). There are some limitations in our study. First of all, the main disadvantage of this study is the lack of in vivo validation. Secondly, this study only included the TCGA database, and more databases should be included for research in the future.
Conclusion
The TCGA database’s analysis of the lncRNA, miRNA, and mRNA expression profiles of HCC resulted in the creation of a ceRNA network. It was discovered that lncRNAs may play a role in carcinogenesis and may help predict the overall survival of patients with HCC. To ascertain their viability as therapeutic or diagnostic biomarkers, it is important to investigate the probable processes of these lncRNAs in HCC. However, further research on the clinical practice of this nomogram combined model and the exact molecular mechanism of the two lncRNAs identified in this study are required.
Data availability statement
The original contributions presented in the study are included in the article/supplementary material. Further inquiries can be directed to the corresponding authors.
Author contributions
WC, FC, MG, ZJ and LS performed the experiments, analyzed the data. WC and JW wrote the manuscript. Z-WW edited the manuscript. JW supervised this study. All authors contributed to the article and approved the submitted version.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Abbreviations
ceRNA, competing endogenous RNA; circRNA, circular RNAs; CSC, cancer stem cells; DSCR4, Down Syndrome critical region 4; EZH2, Enhancer of Zeste 2; GO, gene ontology; HCC, hepatocellular carcinoma; HOTTIP, HOXA distal transcript antisense RNA; HOXA13, homeobox A13; KEGG, The Kyoto Encyclopedia of Genes and Genomes; lncRNA, long noncoding RNA; MEG3, maternally expressed 3; miRNA, microRNA; mRNA, messenger RNA; SLC7A11, solute carrier family 7 member 11; TCGA, The Cancer Genome Atlas; TDRG1, testis development related 1; UCA1, urothelial cancer associated 1.
References
1. Kim YS, Shin SW. Hepatocellular carcinoma. New Engl J Med (2019) 381(1):e2. doi: 10.1056/NEJMc1906565
2. Gan L, Liu Z and Sun C. Obesity linking to hepatocellular carcinoma: A global view. Biochim Biophys Acta (BBA) - Rev Cancer (2018) 1869(2):97–102. doi: 10.1016/j.bbcan.2017.12.006
3. Raoul JL, Edeline J. Systemic treatment of hepatocellular carcinoma: standard of care in China and elsewhere. Lancet Oncol (2020) 21(4):479–81. doi: 10.1016/s1470-2045(20)30082-6
4. Siegel RL, Miller KD, Fuchs HE and Jemal A. Cancer statistics, 2022. CA: A Cancer J Clin (2022) 72(1):7–33. doi: 10.3322/caac.21708
5. Allaire M, Nault J-C. Molecular targets for HCC and future treatments. J Hepatol (2017) 66(1):234–5. doi: 10.1016/j.jhep.2016.07.034
6. Yang JD, Heimbach JK. New advances in the diagnosis and management of hepatocellular carcinoma. BMJ (2020) 371:m3544. doi: 10.1136/bmj.m3544
7. Torimura T, Iwamoto H. Treatment and the prognosis of hepatocellular carcinoma in Asia. Liver Int (2022) 42(9):2042–54. doi: 10.1111/liv.15130
8. Caruso S, O’Brien DR, Cleary SP, Roberts LR, Zucman-Rossi J. Genetics of hepatocellular carcinoma: Approaches to explore molecular diversity. Hepatology (2021) 73(S1):14–26. doi: 10.1002/hep.31394
9. Yang M, Lu H, Liu J, Wu S, Kim P and Zhou X. lncRNAfunc: A knowledgebase of lncRNA function in human cancer. Nucleic Acids Res (2021) 50(D1):D1295–306. doi: 10.1093/nar/gkab1035
10. Nojima T, Proudfoot NJ. Mechanisms of lncRNA biogenesis as revealed by nascent transcriptomics. Nat Rev Mol Cell Biol (2022) 23(6):389–406. doi: 10.1038/s41580-021-00447-6
11. Chen X, Liu Y, Liu H, Wang ZW and Zhu X. Unraveling diverse roles of noncoding RNAs in various human papillomavirus negative cancers. Pharmacol Ther (2022) 238:108188. doi: 10.1016/j.pharmthera.2022.108188
12. Jiang W, Pan S, Chen X, Wang ZW and Zhu X. The role of lncRNAs and circRNAs in the PD-1/PD-L1 pathway in cancer immunotherapy. Mol Cancer (2021) 20(1):116. doi: 10.1186/s12943-021-01406-7
13. Xie W, Chu M, Song G, Zuo Z, Han Z, Chen C, et al. Emerging roles of long noncoding RNAs in chemoresistance of pancreatic cancer. Semin Cancer Biol (2022) 83:303–18. doi: 10.1016/j.semcancer.2020.11.004
14. Liu J, Shang G. The roles of noncoding RNAs in the development of osteosarcoma stem cells and potential therapeutic targets. Front Cell Dev Biol (2022) 10:773038. doi: 10.3389/fcell.2022.773038
15. Chen T, Liu J, Zhang H, Li J and Shang G. Long intergenic noncoding RNA00265 enhances cell viability and metastasis via targeting miR-485-5p/USP22 axis in osteosarcoma. Front Oncol (2022) 12:907472. doi: 10.3389/fonc.2022.907472
16. Jiang W, Xia J, Xie S, Zou R, Pan S, Wang ZW, et al. Long non-coding RNAs as a determinant of cancer drug resistance: Towards the overcoming of chemoresistance via modulation of lncRNAs. Drug Resist Update (2020) 50:100683. doi: 10.1016/j.drup.2020.100683
17. Ren Y, Li R-Q, Cai Y-R, Xia T, Yang M and Xu F-J. Effective codelivery of lncRNA and pDNA by pullulan-based nanovectors for promising therapy of hepatocellular carcinoma. Advanced Funct Materials (2016) 26(40):7314–25. doi: 10.1002/adfm.201603041
18. Zhang L, Yang Z, Trottier J, Barbier O and Wang L. Long noncoding RNA MEG3 induces cholestatic liver injury by interaction with PTBP1 to facilitate shp mRNA decay. Hepatology (2017) 65(2):604–15. doi: 10.1002/hep.28882
19. Quagliata L, Matter MS, Piscuoglio S, Arabi L, Ruiz C, Procino A, et al. Long noncoding RNA HOTTIP/HOXA13 expression is associated with disease progression and predicts outcome in hepatocellular carcinoma patients. Hepatology (2014) 59(3):911–23. doi: 10.1002/hep.26740
20. Kamel MM, Matboli M, Sallam M, Montasser IF, Saad AS and El-Tawdi AHF. Investigation of long noncoding RNAs expression profile as potential serum biomarkers in patients with hepatocellular carcinoma. Trans Res (2016) 168:134–45. doi: 10.1016/j.trsl.2015.10.002
21. Cai J, Zhou M and Xu J. N6-methyladenosine (m6A) RNA methylation regulator SNRPC is a prognostic biomarker and is correlated with immunotherapy in hepatocellular carcinoma. World J Surg Oncol (2021) 19(1):241. doi: 10.1186/s12957-021-02354-8
22. Zhang D-Y, Zou X-J, Cao C-H, Zhang T, Lei L, Qi X-L, et al. Identification and functional characterization of long non-coding RNA MIR22HG as a tumor suppressor for hepatocellular carcinoma. Theranostics (2018) 8(14):3751–65. doi: 10.7150/thno.22493
23. Yin Z, Dong C, Jiang K, Xu Z, Li R, Guo K, et al. Heterogeneity of cancer-associated fibroblasts and roles in the progression, prognosis, and therapy of hepatocellular carcinoma. J Hematol Oncol (2019) 12(1):101. doi: 10.1186/s13045-019-0782-x
24. Zhong Y, Du Y, Yang X, Mo Y, Fan C, Xiong F, et al. Circular RNAs function as ceRNAs to regulate and control human cancer progression. Mol Cancer (2018) 17(1):79. doi: 10.1186/s12943-018-0827-8
25. Salmena L, Poliseno L, Tay Y, Kats L, Pandolfi Pier P. A ceRNA hypothesis: The Rosetta stone of a hidden RNA language? Cell (2011) 146(3):353–8. doi: 10.1016/j.cell.2011.07.014
26. Wang Z, Yang B, Zhang M, Guo W, Wu Z, Wang Y, et al. lncRNA epigenetic landscape analysis identifies EPIC1 as an oncogenic lncRNA that interacts with MYC and promotes cell-cycle progression in cancer. Cancer Cell (2018) 33(4):706–720.e709. doi: 10.1016/j.ccell.2018.03.006
27. Yu H, Yao X, Meng X. P-132 - a novel LncRNA (LOC105371049) regulates colorectal cancer proliferation, metastasis and metabolism. Ann Oncol (2019) 30:: iv35–iv36. doi: 10.1093/annonc/mdz155.131
28. Tomczak K, Czerwińska P and Wiznerowicz M. Review the cancer genome atlas (TCGA): An immeasurable source of knowledge. Contemp Oncology/Współczesna Onkologia (2015) 2015(1):68–77. doi: 10.5114/wo.2014.47136
29. Robinson MD, McCarthy DJ and Smyth GK. edgeR: a bioconductor package for differential expression analysis of digital gene expression data. bioinformatics (2010) 26(1):139–40. doi: 10.1093/bioinformatics/btp616
30. Consortium GO. The gene ontology (GO) database and informatics resource. Nucleic Acids Res (2004) 32(suppl_1):D258–61. doi: 10.1093/nar/gkh036
31. Yu G, Wang L-G, Han Y and He Q-Y. clusterProfiler: an r package for comparing biological themes among gene clusters. Omics: J Integr Biol (2012) 16(5):284–7. doi: 10.1089/omi.2011.0118
32. Wickham H, Chang W and Wickham MH. Package ‘ggplot2’. create elegant data visualisations using the grammar of graphics. Version (2016) 2(1):1–189.
33. Goel MK, Khanna P and Kishore J. Understanding survival analysis: Kaplan-Meier estimate. Int J Ayurveda Res (2010) 1(4):274–8. doi: 10.4103/0974-7788.76794
34. Chen B, Khodadoust MS, Liu CL, Newman AM and Alizadeh AA. Profiling tumor infiltrating immune cells with CIBERSORT. Methods Mol Biol (2018) 1711:243–59. doi: 10.1007/978-1-4939-7493-1_12
35. Lunn M, McNeil D. Applying cox regression to competing risks. Biometrics (1995) 51(2):524–32. doi: 10.2307/2532940
37. Cleophas TJ, Zwinderman AH. Bayesian Pearson Correlation analysis. In: Modern Bayesian statistics in clinical research. Switzerland: Springer (2018). p. 111–8.
38. von Felden J, Garcia-Lezana T, Schulze K, Losic B and Villanueva A. Liquid biopsy in the clinical management of hepatocellular carcinoma. Gut (2020) 69(11):2025–34. doi: 10.1136/gutjnl-2019-320282
39. Llovet JM, Pinyol R, Kelley RK, El-Khoueiry A, Reeves HL, Wang XW, et al. Molecular pathogenesis and systemic therapies for hepatocellular carcinoma. Nat Cancer (2022) 3(4):386–401. doi: 10.1038/s43018-022-00357-2
41. Liang Y, Song X, Li Y, Chen B, Zhao W, Wang L, et al. RETRACTED ARTICLE: LncRNA BCRT1 promotes breast cancer progression by targeting miR-1303/PTBP3 axis. Mol Cancer (2020) 19(1):85. doi: 10.1186/s12943-020-01206-5
42. Wang W, Yang T, Li D, Huang Y, Bai G and Li Q. LINC00491 promotes cell growth and metastasis through miR-324-5p/ROCK1 in liver cancer. J Transl Med (2021) 19(1):504. doi: 10.1186/s12967-021-03139-z
43. Ma Y, Xu XL, Huang HG, Li YF and Li ZG. LncRNA TDRG1 promotes the aggressiveness of gastric carcinoma through regulating miR-873-5p/HDGF axis. Biomedicine Pharmacotherapy (2020) 121:109425. doi: 10.1016/j.biopha.2019.109425
44. Chen S, L-l W, Sun K-x, Liu Y, Guan X, Zong Z-h, et al. LncRNA TDRG1 enhances tumorigenicity in endometrial carcinoma by binding and targeting VEGF-a protein. Biochim Biophys Acta (BBA) - Mol Basis Dis (2018) 1864:3013–21. doi: 10.1016/j.bbadis.2018.06.013
45. Zhang Y, Huang JC, Cai KT, Yu XB, Chen YR, Pan WY, et al. Long noncoding RNA HOTTIP promotes hepatocellular carcinoma tumorigenesis and development: A comprehensive investigation based on bioinformatics, qRTPCR and metaanalysis of 393 cases. Int J Oncol (2017) 51(6):1705–21. doi: 10.3892/ijo.2017.4164
46. Han L, Yan Y, Zhao L, Liu Y, Lv X, Zhang L, et al. LncRNA HOTTIP facilitates the stemness of breast cancer via regulation of miR-148a-3p/WNT1 pathway. J Cell Mol Med (2020) 24(11):6242–52. doi: 10.1111/jcmm.15261
47. Wu L, Yang Z, Zhang J, Xie H, Zhou L and Zheng S. Long noncoding RNA HOTTIP expression predicts tumor recurrence in hepatocellular carcinoma patients following liver transplantation. Hepatobiliary Surg Nutr (2018) 7(6):429–39. doi: 10.21037/hbsn.2018.10.07
48. Tsang FH, Au SL, Wei L, Fan DN, Lee JM, Wong CC, et al. Long non-coding RNA HOTTIP is frequently up-regulated in hepatocellular carcinoma and is targeted by tumour suppressive miR-125b. Liver Int (2015) 35(5):1597–606. doi: 10.1111/liv.12746
49. Ge Y, Yan X, Jin Y, Yang X, Yu X, Zhou L, et al. MiRNA-192 and miRNA-204 directly suppress lncRNA HOTTIP and interrupt GLS1-mediated glutaminolysis in hepatocellular carcinoma. PloS Genet (2015) 11(12):e1005726. doi: 10.1371/journal.pgen.1005726
50. Tang Q, Li X, Chen Y, Long S, Yu Y, Sheng H, et al. Solamargine inhibits the growth of hepatocellular carcinoma and enhances the anticancer effect of sorafenib by regulating HOTTIP-TUG1/miR-4726-5p/MUC1 pathway. Mol Carcinog (2022) 61(4):417–32. doi: 10.1002/mc.23389
51. Dong YJ, Feng W and Li Y. HOTTIP-miR-205-ZEB2 axis confers cisplatin resistance to ovarian cancer cells. Front Cell Dev Biol (2021) 9:707424. doi: 10.3389/fcell.2021.707424
52. Thomson AW, Drakes ML, Zahorchak AF, O’Connell PJ, Steptoe RJ, Qian S, et al. Hepatic dendritic cells: Immunobiology and role in liver transplantation. J Leukoc Biol (1999) 66(2):322–30. doi: 10.1002/jlb.66.2.322
53. Doherty DG, O’Farrelly C. Dendritic cells: regulators of hepatic immunity or tolerance? J Hepatol (2001) 34(1):156–60. doi: 10.1016/s0168-8278(00)00020-9
54. Khanna A, Morelli AE, Zhong C, Takayama T, Lu L and Thomson AW. Effects of liver-derived dendritic cell progenitors on Th1- and Th2-like cytokine responses in vitro and in vivo. J Immunol (2000) 164(3):1346–54. doi: 10.4049/jimmunol.164.3.1346
55. Abe M, Akbar SM, Horiike N and Onji M. Induction of cytokine production and proliferation of memory lymphocytes by murine liver dendritic cell progenitors: Role of these progenitors as immunogenic resident antigen-presenting cells in the liver. J Hepatol (2001) 34(1):61–7. doi: 10.1016/s0168-8278(00)00004-0
Keywords: HCC, lncRNA, HOTTIP, MiR-205, prognosis, immune microenvironment
Citation: Chen W, Chen F, Gong M, Jin Z, Shu L, Wang Z-w and Wang J (2022) Comprehensive analysis of lncRNA-mediated ceRNA networkfor hepatocellular carcinoma. Front. Oncol. 12:1042928. doi: 10.3389/fonc.2022.1042928
Received: 13 September 2022; Accepted: 30 September 2022;
Published: 21 October 2022.
Edited by:
Aamir Ahmad, University of Alabama at Birmingham, United StatesReviewed by:
Mohd Farhan, King Faisal University, Saudi ArabiaQiang You, Guangzhou Medical University, China
Copyright © 2022 Chen, Chen, Gong, Jin, Shu, Wang and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jianjiang Wang, 13615713000@163.com; Zhi-wei Wang, wangzhiwei@wmu.edu.cn