- Research
- Open access
- Published:
Combining a machine-learning derived 4-lncRNA signature with AFP and TNM stages in predicting early recurrence of hepatocellular carcinoma
BMC Genomics volume 24, Article number: 89 (2023)
Abstract
Background
Near 70% of hepatocellular carcinoma (HCC) recurrence is early recurrence within 2-year post surgery. Long non-coding RNAs (lncRNAs) are intensively involved in HCC progression and serve as biomarkers for HCC prognosis. The aim of this study is to construct a lncRNA-based signature for predicting HCC early recurrence.
Methods
Data of RNA expression and associated clinical information were accessed from The Cancer Genome Atlas Liver Hepatocellular Carcinoma (TCGA-LIHC) database. Recurrence associated differentially expressed lncRNAs (DELncs) were determined by three DEG methods and two survival analyses methods. DELncs involved in the signature were selected by three machine learning methods and multivariate Cox analysis. Additionally, the signature was validated in a cohort of HCC patients from an external source. In order to gain insight into the biological functions of this signature, gene sets enrichment analyses, immune infiltration analyses, as well as immune and drug therapy prediction analyses were conducted.
Results
A 4-lncRNA signature consisting of AC108463.1, AF131217.1, CMB9-22P13.1, TMCC1-AS1 was constructed. Patients in the high-risk group showed significantly higher early recurrence rate compared to those in the low-risk group. Combination of the signature, AFP and TNM further improved the early HCC recurrence predictive performance. Several molecular pathways and gene sets associated with HCC pathogenesis are enriched in the high-risk group. Antitumor immune cells, such as activated B cell, type 1 T helper cell, natural killer cell and effective memory CD8 T cell are enriched in patients with low-risk HCCs. HCC patients in the low- and high-risk group had differential sensitivities to various antitumor drugs. Finally, predictive performance of this signature was validated in an external cohort of patients with HCC.
Conclusion
Combined with TNM and AFP, the 4-lncRNA signature presents excellent predictability of HCC early recurrence.
Introduction
The recent global cancer statistics study indicated that the new cases and deaths of liver cancer were 905,677 and 830,187 respectively and rank sixth in terms of incidence and third in terms of mortality [1]. Approximately 75–85% of primary liver cancer cases are caused by hepatocellular carcinoma (HCC) [1]. Although the main risk factors of HCC show regional differences, chronic hepatitis B or C infection remains the major causes of HCC [2, 3]. Benefit from vaccination against HBV, the incidence of HCC in high-risk countries of Eastern Asia has been dramatically reduced [3]. However, incidence rates of HCC in regions like Europe, Northern and South America, Australia/New Zealand, which were low-risk regions display the opposite trend or remain at a high-level plateau [4]. Thus, the overall global burden of liver cancer is increasing over time.
Long non-coding RNAs (lncRNA), a class of non-coding transcripts that exceed 200 nucleotides in length, have been identified as important regulators in the development of various cancers [5, 6]. Cancer-related lncRNAs are involved in genomic instability, sustained proliferation, activation of invasion and metastasis, and cell death resistance in cancer cells by means of diverse mechanisms [7] through binding with RNA, DNA, protein or encoding small peptides [8,9,10]. For example, lncRNA GHET1 promoted HCC cell tumorigenesis by activating H3K27 acetylation and regulating ATF1 [11]. LINC01234, a potential prognostic or therapeutic HCC marker, could modulate aspartate metabolic reprogramming and promote HCC progression [12]. Cancer-associated fibroblasts secreted exosomal lncRNA TUG1 facilitated HCC cell glycolysis, migration, and invasion via miR-524-5p/SIX1 axis [13]. LncRNA DANCR promoted HCC stemness by regulating mRNA stabilization [14]. LncRNA PVT1 promotes HCC cell proliferation and stemness by stabilizing NOP2 [15]. Besides experimental evidences, more lncRNA-disease associations were elucidated by powerful bioinformatics tools and models, which could help to underlie disease mechanisms at the level of lncRNA and facilitate the detection of biomarkers for diagnosis and prognosis, as well as disease prevention and treatment [16, 17]. Therefore, several lncRNAs have been reported to serve as diagnosis and prognosis biomarkers for HCC. Lnc-PCDH9-13:1 was upregulated in HCC tissues, serum and saliva of the patients and could serve as a biomarker for detecting early HCC [18]. Circulating exosomal lncRNA-ATB was identified as an independent predictor HCC overall survival and disease progression [19]. A panel of serum circulating lncRNA LINC00153, UCA1 and AFP was reported to have satisfactory sensitivity and specificity for HCC diagnosis [20]. A signature consisting of 50 lncRNA pairs could serve as an independent powerful prognostic indicator for HCC overall survival prediction [21]. Moreover, lncRNA signatures associated with genome instability, macrophages [22, 23], pyroptosis, ferroptosis, tumor microenvironment, m6A regulator, autophagy, hypoxia, glycolysis and EMT [24,25,26,27,28,29,30,31] were established for predicting overall survival in HCC.
Nearly 70% HCC patients had postsurgical recurrence in 5 years. Postsurgical recurrence is the primary limitation for the improvement of HCC prognosis [32]. Clinically, a recurrence within two years of surgery is referred to an early recurrence while a recurrence after two years is called a late recurrence [33]. Earlier studies has indicated that near 70% of recurrence was early recurrence, and HCC patients with early recurrence had a significantly lower 5-year overall survival rate compared to HCC patients with late recurrence [34]. Therefore, construction of signature predicting HCC early recurrence would enable an improved surveillance strategy and prognosis.
In this study, we collected data from TCGA-LIHC databases to construct a 4-lncRNA prognostic signature for HCC early recurrence. Multivariate Cox regression, Kaplan-Meier, nomogram and ROC analyses were performed to evaluate the predictive potential of this signature. KEGG, GO, GSEA were performed to explore the underlying mechanism of HCC early recurrence. Intratumor immune infiltration status and drug response prediction analyses were used to evaluate the potential of this signature in predicting therapeutic responses. In addition, the significance of this signature was further validated in external HCC cohorts (Figure S1).
Materials and methods
Data mining for candidate recurrence related dysregulated lncRNAs in hepatocellular carcinoma
Liver cancer RNA expression data and associated clinical features were downloaded from The Cancer Genome Atlas Liver Hepatocellular Carcinoma (TCGA-LIHC) database (https://portal.gdc.cancer.gov/projects/TCGA-LIHC). The expression profile and clinical features of 314 HCC patients with complete overall survival (OS) and disease free survival (DFS) record were reserved after carefully screening. These 314 patients were then randomly divided into the training group (N = 157) and the validation group (N = 157) by R package “caret” [35]. Clinical features of these patients including OS, DFS, TNM stages, cirrhotic status, vascular invasion, AFP, race, gender and age have been summarized (Table S1). Next, the differentially expressed lncRNAs (DElncs) were analyzed between the training group (N = 157) and the non-tumor group (N = 50). Differentially expressed gene (DEG) analyses were conducted by R packages “DESeq2” [36], “edgeR” [37, 38] and “limma” [39] with the cut-off value of |log2FC| > 1 and FDR < 0.05. Venn plots were drawn to identify the common dysregulated lncRNAs in all three DEG methods by R package “VennDiagram” [40]. Batch DFS survival analyses of the above dysregulated lncRNAs in the training group were then performed by R package “survival” [41] with both log-rank [42] and cox [43] methods with a cut-off value of P < 0.05. Finally, 81 candidate recurrence related dysregulated lncRNAs were obtained from the intersection of two survival analyses methods.
Construction and validation for lncRNA-based prognostic signature for HCC recurrence
Dimensionality reduction of candidate lncRNAs used for signature construction were further conducted by three different machine learning methods including least absolute shrinkage and selection operator (LASSO) [44], random forest [45] and Support Vector Machine Recursive Feature Elimination (SVM-RFE) [46] in the training group with 81 candidate recurrence related dysregulated lncRNAs and DFS. LASSO, random forest and SVM-RFE were conducted by using R package “glmnet” [47], “randomForest” [45], and “e1071” [48] separately. More specifically, the cut-off values were lambda.min for LASSO, 5-fold cross validation with min(error) and max(accuracy) for SVM-RFE and top 30 DFS related lncRNAs for random forest. Venn plot was then drawn in three machine learning methods to identify 11 lncRNAs. Next, multivariate cox analysis of these 11 lncRNAs was performed in the training group (N = 157) with DFS by R package “survival” and 4 lncRNAs (AC108463.1, AF131217.1, CMB9-22P13.1, TMCC1-AS1) with a cut-off value of P < 0.05 were finally obtained for signature construction. The coefficients of these 4 lncRNAs for signature construction were calculated by multivariate cox analysis in the training group with DFS, and the risk score (RS) for each patient was calculated by the formula \(risk score= \sum coefficient \times expression \left(gene\right)\). Receiver Operating Characteristic (ROC) analysis was performed by R package “pROC” [49] to evaluate the performances of the 4-lncRNA signature and three clinical features including AFP, TNM, as well as vascular invasion for HCC recurrence in the training group. All HCC patients were then further divided into the low- and high-risk group by median risk score from the training group.
Tissue samples and clinical information
A total number of 44 patients who had liver surgery and were diagnosed with HCC between October 2018 and December 2019 at Jinling Hospital were included. HCC and paracancerous tissues were collected and treated following the protocols (81YY-KYLL-19-05) approved by the Ethics Committee of Jinling Hospital and Shanghai University of Medicine and Health Sciences. In addition, clinical characteristics including age, gender, serum AFP, numbers of tumor, cirrhotic status, vascular invasion and T stage were collected or analyzed. The follow-up data of 24 patients with available overall survival and disease free survival record until March 2022 were enrolled for survival analyses (Table S5). All patients enrolled in this study signed the written informed consent.
Validation of the lncRNA-based prognostic signatures for HCC early recurrence
Since HCC early recurrence was defined as recurrence within 2 years in previous studies, those patients with follow-up time less than 24 months and no recurrence were excluded for the following studies of the 4-lncRNA signature with HCC early relapse. The HCC patients were 112 in the training group, 111 in the validation group and 223 in the entire group from the TCGA-LIHC cohort. To validate the prediction performance of the 4-lncRNA signature for HCC early recurrence, early recurrence cumulative event, cumulative hazard and probability were compared between the low- and high-risk group by R package “survminer” [50] in the training, validation and entire group. Multivariate cox analyses were also conducted to study the independent roles of the 4-lncRNA signature along with clinical characteristics like AFP, TNM stages, vascular invasion, gender and age. ROC analyses further identified that a combination model of lncRNA-based signature with AFP/TNM had the best prognostic prediction for HCC early recurrence. Nomogram was generated with the 4-lncRNA signature risk score, AFP and TNM stages and their corresponding coefficients from multivariate cox analyses, and the calibration curves were drawn by R package “regplot” [51]. Since no public external HCC dataset was available, we did further validation in clinical HCC samples by detecting the expression of 4 lncRNAs in clinical collected tumor and paired paracancerous tissues. Total RNA from 44 paired HCC and paracancerous tissue samples were collected by Jinling Hospital. cDNA was synthesized from total RNA by using the reverse transcription kit ReverTra Ace® qPCR RT Master Mix with gDNA Remover (TOYOBO). Real-time PCR reaction was conducted in 20 µL solution with Takara CYBR Premix Ex TaqII (Takara) in ABI BiosystemsTM 7500 Real-Time qPCR System (Applied Biosystems) by following the manufacturer’s protocol. Primers for AC108463.1, AF131217.1, CMB9-22P13.1, TMCC1-AS1 and 18s (Table S2) were purchased from GENEWIZ. To calculate the relative expression of lncRNAs, qRT-PCR results were interpreted by 2−ΔΔCT method with 18s as the housekeeping gene. Furthermore, 24 HCC patients from Jinling cohort were enrolled for early recurrence analysis following the same criteria of TCGA-LIHC cohort.
Comprehensive functional analyses of the 4-lnRNA prognostic signature
Total 223 HCC patients were divided into the low- and high-risk group by setting the 4-lncRNA prognostic signature median risk score as the cut-off. DEG analysis was performed between mRNA expression from the low- and high risk group, and all mRNAs were then arranged by descending |log2FC| for functional analyses. Kyoto Encyclopedia of Genes and Genomes (KEGG) [52], Gene Ontology (GO) and Gene Set Enrichment Analysis (GSEA) analyses were conducted by R package “clusterProfiler” [53]. The enriched KEGG pathways were determined by a cut-off value of |NES| > 1 and P < 0.05. The enriched GO terms included biological pathways (BP), cellular components (CC) and molecular functions (MF) were analyzed based on mRNAs with |logFC| > 1 and determined by a cut-off value of P < 0.05. GSEA analysis was performed with MSigDb C2: curated gene sets [54, 55] and enriched GSEA gene sets were determined by a cut-off value of |NES| > 1.5 and P < 0.05.
Immune infiltration and clinical treatment response analyses
Single sample Gene Set Enrichment Analysis (ssGSEA) was chosen for studying immune infiltration with R package “GSVA” [56] and normalized enrichment score (NES) was calculated for 28 immune cell types in each 223 HCC samples. The NES of 28 immune cell types were compared between the low- and the high-risk group, and the correlation between the 4-lncRNA signature risk score and cells NES score was conducted by R function “cor.test”. Immune therapy response prediction was conducted with Tumor Immune Dysfunction and Exclusion (TIDE) algorithm [57, 58] and SubMap modules of GenePattern [59] by mapping with a public dataset of immunotherapy responses of 47 melanoma patients [60]. Drug response prediction was performed by R package “pRRophetic” [61] with Genomics of Drug Sensitivity in Cancer (GDSC) pharmacogenomics database [60, 62]. Ridge regression was used for estimating the half maximal inhibitory concentration (IC50) and 10-fold cross validation was used for predicting the accuracy.
Statistical analyses
The sensitivity and specificity between two ROC curves were compared by DeLong’s test. The differences between Kaplan-Meier curves, cumulative hazard and events curves of survival analyses between two groups were compared by log-rank test. Univariate and multivariate analyses were conducted with the cox proportional hazards regression model. The comparisons of immune cells NES, drug IC50 and between the two groups, as well as the expression of lncRNAs between tumor and paracancerous tissue samples were analyzed by Wilcoxon test. The cut-off value of P < 0.05 was used for statistical significance.
Results
Identification of recurrence related lncRNAs
To construct a lncRNA signature to predict postsurgical recurrence of HCC, we started to identify dysregulated lncRNAs in the TCGA training group. According to the sequences of 16,193 annotated human lncRNAs in GENCODE V30, we collected 10,795 lncRNAs for DEG analysis after excluding those with extremely low expression. Three methods, DESeq2, edgeR and limma-voom, were employed to identify differentially expressed lncRNAs (DElncs) between the HCC samples of the training group (N = 157) and the liver tissues of the normal control group (N = 50) with cut-off value of |log2FC| > 1 and FDR < 0.05. Compared with normal controls, 2581 (2013 upregulated and 568 downregulated), 3430 (2913 upregulated and 517 downregulated) and 1631 (824 upregulated and 807 downregulated) DElncs were determined respectively (Figure S2). Venn diagram analysis revealed 1164 (801 upregulated and 363 downregulated) common DElncs (Fig. 1A). The PCA plot and heatmap generated by the 1164 common DElncs could clearly distinguish HCCs from normal controls (Fig. 1B and C), suggesting that the 1164 common DElncs might closely associate with the onset and development of HCC. To identify recurrence associated DElncs, the log-rank test and Cox regression analysis were performed in the training group to evaluate disease free survival (DFS). After combining the candidates from the log-rank test (149 DElncs) and Cox regression analysis (136 DElncs), 81 common DElncs were identified as recurrence related lncRNAs (Fig. 1D).
Construction of a 4-lncRNA prognostic signature for HCC recurrence
Based on the 81 recurrence associated lncRNAs, we then employed three classic machine learning methods, LASSO, Random Forest and SVM-RFE, to select important DElncs for predicting DFS in the training group. The LASSO, SVM-REF and Random Forest analyses selected 26, 66 and 30 candidates respectively (Fig. 2A to C). Venn diagram analysis collected 11 common Delncs for further analysis (Fig. 2D). Multivariate cox analysis of the 11 DElncs in the training group showed that 4 DElncs, AC108463.1, AF131217.1, CMB9-22P13.1, and TMCC1-AS1, are independent risk factors of DFS (Fig. 2E). We then constructed a prognostic signature based on the 4 DElncs, and calculated the risk score of individual HCC patients according to the linear combination of the regression coefficients and expression values of each DElncs [63]. Risk Score = (-0.0918*exp[AC108463.1]) + (-0.1112*exp[AF131217.1]) + (0.1484*exp[CMB9-22P13.1]) + (0.3737*exp[TMCC1-AS1]) (Table 1).
The 4-lncRNA prognostic signature predicts HCC early recurrence
Since HCC patients’ poor survival is largely attributed to early recurrence within two years after surgery, we intended to investigate whether the 4-lncRNA signature could provide a prognostic indication of HCC patients’ early recurrence. After excluding the patients whose follow-up data were collected less than 2-year and without recurrence records, 112 and 111 HCC patients were reserved in the training and the validation groups respectively. By setting the median risk score as a cut-off, HCC patients were categorized into the low- and high-risk groups. Kaplan-Meier survival analyses demonstrated that HCC patients from the high-risk group had shorter 2-year DFS in the training group (Fig. 3A, P < 0.0001), the validation group (Fig. 3B, P = 0.033) and the entire group (Fig. 3C, P < 0.0001). The HCC early recurrence rates in the high-risk group were up to 79% from all three groups (Fig. 3). In addition, cumulative event and cumulative hazard of the low-group patients were also compared to those of high-risk group patients. HCC patients in the high-risk group showed the higher cumulative events and cumulative hazard from all three groups (Figure S3). Thus, the findings by survival analyses indicate that the 4-lncRNA prognostic signature could predict HCC early recurrence within 2 years.
Combination of the 4-lncRNA signature risk score with AFP and TNM improves the prognostic performance for HCC early recurrence
To further evaluate the prognostic value of the 4-lncRNA signature, multivariate cox analyses of the risk score together with selected clinical features, including age, gender, AFP level, TNM stage and vascular invasion, were conducted in all 223 HCC patients with 2-year DFS. As shown in Fig. 4A, multivariate cox analyses suggest that the risk score (HR = 1.5, P = 0.015), AFP (HR = 1.74, P = 0.012) and TNM (HR = 2.01, P = 0.01 for stage III + IV) were independent risk indicators of HCC early recurrence. ROC analyses were then used for determining whether the combination of the independent risk indicators could improve prognostic performance. As shown in Fig. 4B, the combination of risk score with AFP and TNM showed the largest AUC (72.02%) for HCC early recurrence compared to risk score alone (AUC: 64.89%), risk score + AFP (AUC: 66.85%), and risk score + TNM (AUC: 70.80%). Meanwhile, the signature risk score, AFP and TMN stages were selected to establish a nomogram (Fig. 4C). The C-index for 1-year and 2-year DFS of this nomogram were 0.643 and 0.647, respectively. Moreover, a calibration curve revealed that the nomogram was good at predicting 1-year and 2-year DFS (Fig. 4D).
Enriched KEGG pathways, GO terms and GSEA gene sets in the low- and high-risk groups
To further understand the functional differences between the low- and high-risk groups, KEGG, GO and GSEA analyses were conducted based on the differences of mRNA expression between the two groups. The representative pathways activated in the high-risk group included “IL-17 signaling pathway”, “Pentose phosphate pathway”, “Pentose and glucuronate interconversions”, “Viral protein interaction with cytokine and cytokine receptor”, “NOD-like receptor signaling pathway” and “Transcriptional misregulation in cancer”, while the pathways suppressed in the high-risk group included “Aldosterone synthesis and secretion”, “Alanine, aspartate and glutamate metabolism”, “Vascular smooth muscle contraction” and “Glycosaminoglycan biosynthesis - heparan sulfate/ heparin” (Fig. 5A). Several GO terms from biological process (BP), cellular component (CC) and molecular function (MF) were also significantly activated or suppressed in the high-risk group. For example, transporter complex was activated in the high-risk group, while apical plasma membrane, presynaptic membrane and carbohydrate binding were suppressed in the high-risk group (Fig. 5B). Besides, GSEA analysis with C2 gene sets revealed activated and suppressed gene sets in the high-risk group (Fig. 5C). The most significantly activated gene sets in the high-risk group were “BOSCO_EPITHELIAL_DIFFERENTIATION_MODULE”, “CROMER_TUMORIGENESIS_DN” and “ANDERSEN_CHOLANGIOCARCINOMA_CLASS2”, while the most significantly suppressed gene sets in the high-risk group were “REACTOME_METALLOTHIONEINS_BIND_METALS”, “REACTOME_RESPONSE_TO_METAL_IONS”, “BOYAULT_LIVER_CANCER_SUBCLASS_G6_UP”, “CHIANG_LIVER_CANCER_SUBCLASS_CTNNB1_UP”, “BOYAULT_LIVER_CANCER_SUBCLASS_G123_DN” and “DESERT_PERIVENOUS_HEPATOCELLULAR_CARCINOMA_SUBCLASS_UP”.
Therapeutic responses prediction by the 4-lncRNA signature
Immunotherapy, as a novel treatment approach, has shown advantages in improving OS and DFS in HCC patients [64]. To investigate the potential of the 4-lncRNA signature in immunotherapy response prediction, the immune infiltration were compared between the low- and high-risk groups by calculating the NES of 28 immune cell types with ssGSEA. As shown in Fig. 6A, the intratumor infiltration of 10 immune cell types, including activated B cells, effector memory CD8 T cells, eosinophils, immature B cells, macrophages, mast cells, myeloid-derived suppressor cells, natural killer cells, regulatory T cells and type 1 T help cells, showed significantly higher NES in the low-risk group, whereas type 2 T help cells had significantly higher NES in the high-risk group. Correlation analyses illustrated that the NES of the above 10 immune cell types were negatively associated with risk scores (P < 0.05, Fig. 6B). Type 1 T helper cell, activated B cell and natural killer cell ranked the top3 negatively associated tumor infiltration lymphocytes (TILs) among them (|NES| > 0.25). However, although HCCs in the low-risk group showed greater immune cell infiltration, both TIDE predication and Submap analyses failed to show significant response advantages to anti-CTLA4 and anti-PD1 immunotherapy in the low-risk group (Fig. 6C). In addition, drug response prediction analysis indicated that the low-risk group HCC patients and high-risk group HCC patients showed significant differential responses to 27 drugs (Table S4). The low-risk group HCC patients may be more sensitive to AICAR, gefitinib and metforminin treatment, whereas the high-risk group HCC patients may better response to bexarotene, bleomycin, bortezomib, cisplatin, mitomycin C, paclitaxel, sorafenib, tipifarnib and vinorelbine (Fig. 6D). Thus, these findings indicate that the 4-lncRNA signature might act as a potential drug predictor.
Validation of the 4-lncRNAs prognostic signature in clinical HCC samples
Additionally, the expressions of these 4 lncRNAs were measured in a cohort of 44 paired HCC and their adjacent tissue samples. In line with our findings from TCGA-LIHC, AC108463.1, CMB9-22P13.1 and TMCC1-AS1 are highly expressed, and AF131217.1 is suppressed in HCC tissues compared to matched normal tissues (Fig. 7A-D). Furthermore, multivariate Cox analysis was conducted in 24 HCC patients whose 2-year DFS information was recorded. Results revealed that the 4-lncRNA signature risk score, T stage and AFP were independent risk factors for HCC early recurrence (Fig. 7E). Moreover, the 24 HCC patients were classified into low- and high-risk groups based on the combined model of 4-lncRNA signature risk score with AFP and T stage, Kaplan-Meier analysis illustrated that patients in the high-risk group had significantly poorer DFS within 2 years compared to those in the low-risk group (Fig. 7F, P < 0.0001).
Discussion
In the current study, we developed a novel 4-lncRNA prognostic signature for early recurrence prediction in HCC by combining multiple DEG analysis methods, survival analysis methods and machine learning methods. This 4-lncRNA signature could fairly predict the early recurrence of HCC in TCGA-LIHC cohort, and the prediction performance could be further improved by the combination of the 4-lncRNA signature with TNM stages and AFP. According to the risk scores derived from the signature, HCC patients could be categorized into low- and high-risk groups. Functional analyses including KEGG, GO and GSEA were conducted to reveal the underlying mechanisms for HCC early recurrence. Moreover, immune infiltration analysis was employed to find out the immune microenvironment differences between the two groups, and prediction analyses of immune therapy and drug response provided useful information for differential clinical treatment. Finally, the prognostic performance of this 4-lncRNA signature was evaluated in an external HCC cohort.
DESeq2, edgeR and limma-voom are three popularly adopted approaches for DEG analysis [65]. DESeq2 uses shrinkage estimators for dispersion and fold change [36], edgeR adopts a Poisson super dispersion model for account for both biological and technical variability [37], and limma-voom is based on the linear model [66]. In the current study, we employed all those three statistical methods with the same fold for detecting differential expressed lncRNAs (DElncs) between the HCC samples (N = 157) and normal samples (N = 50). The DEG analyses results showed that edgeR found the most DElncs (3430), while limma-voom found the least DElncs (1631). The purpose of time-to-event survival analysis was to find out the DElncs associated with disease free survival. We imported two common survival analyses methods, log-rank and cox, and determined 81 DFS related DElncs [67]. Log-rank test is a non-parametric test for comparing the differences in survival between groups of patients [42], while cox proportional-hazards model is a semiparametric regression model for investigating the impact of variables on patient’s survival [43]. Further dimensionality reduction of DFS related DElncs was conducted by three different machine learning methods including LASSO, random forest and SVM-RFE. LASSO is a method for estimation in liner models with favorable properties of both subset selection and ridge regression [44]. Random forest constructs regression trees in the way of using the best among a subset of predictors randomly chosen at each node to be split [45]. SVM is a popular tool for nonlinear classification, regression and outlier detection [68], and SVM-RFE uses the weight magnitude as ranking criterion [46]. All these three machine learning methods have been universally used in gene selection with their advantages [69,70,71,72,73,74]. For instance, Xiao et al. utilized random forest and SVM for screening prognostic gene in malignant pleural mesothelioma [71], Shen et al. applied SVM for evaluating selective mutant genes and constructing a model for predicting HCC DFS [70], Xiao et al. used SVM-RFE and LASSO for identify candidate hub genes related to colorectal cancer [72]. Our group had also developed a 25-lncRNA signature for predicting the early recurrence of HCC patients by LASSO, but 25 lncRNAs are too many for further validation and clinical application [75]. Considering the individual advantage of LASSO, random forest and SVM-RFE, we have chosen these three popular machine learning methods for feature gene selection in this study.
Kaplan-Meier plot confirmed that the 4-lncRNA signature could successfully classified HCC patients into the low- and high-risk groups and predict early recurrence. Several previous developed lncRNA-based signatures were reported for HCC survival prediction and showed better performance than clinicopathological factors [76,77,78,79,80]. For example, a 3-lncRNA signature by Gu et al. could well predict both recurrence free survival and overall survival in small HCC patients [76], a 15-lncRNA classifier by Zhang et al. effectively identified HCC recurrence [77], a 7-lncRNA classifier by Lv et al. predict early recurrence within 2 years after surgical resection for HCC [78]. The prognostic performance of this 4-lncRNA signature is comparable to other developed lncRNA-based signatures [76,77,78,79,80], while the differences of selected lncRNAs in each lncRNA-based signature might be due to the specific feature gene selection strategy and different training set. Moreover, the combined model of the 4-lncRNA signature risk score, AFP and TNM stages further improve the 2-year DFS prediction with an AUC of 72.02%. Functional analyses were performed to explore the differences between the low- and high-risk groups. Some KEGG pathways activated in the high-risk group are favorable in HCC pathogenesis. For example, IL-17 was reported to promote hepatocellular carcinoma progression [81, 82], the pentose phosphate pathway (PPP) is one of the essential components of cellular metabolism and plays a key role in HCC [83, 84]. Immune infiltration analyses were performed with ssGSEA to elucidate the heterogeneous immune environment in the low- and high-risk group HCC patients. More immune cells had higher NES in the low-risk group and were negatively associated with the 4-lncRNA signature risk score. In addition, the top3 TILs which are negatively related to risk score, Type 1 T helper cell, activated B cell and natural killer cell, are well-known antitumor immune cells participated in cancer immune therapy process [85,86,87], further indicating that this 4-lncRNA signature could potentially predict HCC patients’ prognosis after surgery. Although this signature failed in predicting the response to cancer immunotherapy of HCC patients, drug response prediction suggested that the low-risk group patients are more sensitive to AICAR, gefitinib and metformin treatments, whereas the high-risk group patients are more sensitive to bexarotene, bleomycin, bortezomib, cisplatin, mitomycin C, paclitaxel, sorafenib, tipifarnib and vinorelbine.
In this study, the expression regulation of the 4 lncRNAs has been validated in an external HCC cohort containing 44 paired tumor and matched normal samples. Multivariate cox analysis demonstrated that the risk score of this signature, T stage and AFP are three independent risk indicators for HCC early recurrence in this external cohort. An integrated model by combining this signature, T stage and AFP showed great prognostic potential in predicting HCC early recurrence in this external cohort. Additionally, the 4 lncRNAs involved in this signature have also been previously studied in HCC or other diseases. For instance, AF131217.1 was reported as a fluid shear force-sensitive RNA, which plays a protective role in atherosclerosis process [88]. AC108463.1 is related to gastric cancer progression [89]. CMB9-22P13.1 participates in the development of various cancer types including lung squamous cell carcinoma, breast cancer and hepatocellular carcinoma [90,91,92]. A recent study indicated that CMB9-22P13.1 could upregulate HOTTIP and activate HIF-1α/VEGF signaling, leading to enhanced hepatocellular carcinoma progression and angiogenesis [93]. The increase of TMCC1-AS1 facilitates proliferation, migration, invasion and EMT of HCC cells, resulting in poor outcome of liver cancer patients [94]. Given that the 4 lncRNAs have been selected to construct a prognostic signature for predicting HCC early recurrence, their roles in HCC progression should be intensively investigated. A variety of lncRNA database and developed computational models could be applied for lncRNA function and lncRNA-disease association prediction [9, 10, 16, 17], which may be further validated by experiments.
In conclusion, we developed a 4-lncRNA signature for predicting early recurrence in HCC. The integrated model of the 4-lncRNA signature risk score with TNM and AFP presents great prognostic performance for predicting HCC early recurrence. This signature might provide novel prognostic and therapeutic biomarkers for HCC and act to be a potential drug predictor.
Data availability
The dataset supporting the conclusions of this article is available in the TCGA-LIHC repository, http://cancergenome.nih.gov/.
Abbreviations
- HCC:
-
Hepatocellular carcinoma
- lncRNA:
-
Long non-coding RNA
- TCGA-LIHC:
-
The Cancer Genome Atlas Liver Hepatocellular Carcinoma
- OS:
-
overall survival
- DFS:
-
disease free survival
- DElncs:
-
differentially expressed lncRNAs
- DEG:
-
Differentially expressed gene
- LASSO:
-
least absolute shrinkage and selection operator
- SVM-RFE:
-
Support Vector Machine Recursive Feature Elimination
- RS:
-
risk score
- ROC:
-
Receiver Operating Characteristic
- KEGG:
-
Kyoto Encyclopedia of Genes and Genomes
- GO:
-
Gene Ontology
- GSEA:
-
Gene Set Enrichment Analysis
- BP:
-
biological pathways
- CC:
-
cellular components
- MF:
-
molecular functions
- ssGSEA:
-
Single sample Gene Set Enrichment Analysis
- NES:
-
normalized enrichment score
- TIDE:
-
Tumor Immune Dysfunction and Exclusion
- GDSC:
-
Genomics of Drug Sensitivity in Cancer
References
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(3):209–49.
Yang JD, Hainaut P, Gores GJ, Amadou A, Plymoth A, Roberts LR. A global view of hepatocellular carcinoma: trends, risk, prevention and management. Nat Rev Gastroenterol Hepatol. 2019;16(10):589–604.
de Martel C, Maucort-Boulch D, Plummer M, Franceschi S. World-wide relative contribution of hepatitis B and C viruses in hepatocellular carcinoma. Hepatology. 2015;62(4):1190–200.
Petrick JL, Florio AA, Znaor A, Ruggieri D, Laversanne M, Alvarez CS, et al. International trends in hepatocellular carcinoma incidence, 1978–2012. Int J Cancer. 2020;147(2):317–30.
Adams BD, Parsons C, Walker L, Zhang WC, Slack FJ. Targeting noncoding RNAs in disease. J Clin Invest. 2017;127(3):761–71.
Schmitt AM, Chang HY. Long noncoding RNAs in Cancer Pathways. Cancer Cell. 2016;29(4):452–63.
Huarte M. The emerging role of lncRNAs in cancer. Nat Med. 2015;21(11):1253–61.
Huang Z, Zhou JK, Peng Y, He W, Huang C. The role of long noncoding RNAs in hepatocellular carcinoma. Mol Cancer. 2020;19(1):77.
Zhou L, Wang Z, Tian X, Peng L. LPI-deepGBDT: a multiple-layer deep framework based on gradient boosting decision trees for lncRNA-protein interaction identification. BMC Bioinformatics. 2021;22(1):479.
Peng L, Tan J, Tian X, Zhou L, EnANNDeep. An ensemble-based lncRNA-protein Interaction Prediction Framework with adaptive k-Nearest neighbor classifier and deep models. Interdiscip Sci. 2022;14(1):209–32.
Ding G, Li W, Liu J, Zeng Y, Mao C, Kang Y, et al. LncRNA GHET1 activated by H3K27 acetylation promotes cell tumorigenesis through regulating ATF1 in hepatocellular carcinoma. Biomed Pharmacother. 2017;94:326–31.
Chen M, Zhang C, Liu W, Du X, Liu X, Xing B. Long noncoding RNA LINC01234 promotes hepatocellular carcinoma progression through orchestrating aspartate metabolic reprogramming.Mol Ther. 2022.
Lu L, Huang J, Mo J, Da X, Li Q, Fan M, et al. Exosomal lncRNA TUG1 from cancer-associated fibroblasts promotes liver cancer cell migration, invasion, and glycolysis by regulating the miR-524-5p/SIX1 axis. Cell Mol Biol Lett. 2022;27(1):17.
Yuan SX, Wang J, Yang F, Tao QF, Zhang J, Wang LL, et al. Long noncoding RNA DANCR increases stemness features of hepatocellular carcinoma by derepression of CTNNB1. Hepatology. 2016;63(2):499–511.
Wang F, Yuan JH, Wang SB, Yang F, Yuan SX, Ye C, et al. Oncofetal long noncoding RNA PVT1 promotes proliferation and stem cell-like property of hepatocellular carcinoma cells by stabilizing NOP2. Hepatology. 2014;60(4):1278–90.
Chen X, Yan CC, Zhang X, You ZH. Long non-coding RNAs and complex diseases: from experimental results to computational models. Brief Bioinform. 2017;18(4):558–76.
Chen X, Sun YZ, Guan NN, Qu J, Huang ZA, Zhu ZX, et al. Computational models for lncRNA function prediction and functional similarity calculation. Brief Funct Genomics. 2019;18(1):58–82.
Xie Z, Zhou F, Yang Y, Li L, Lei Y, Lin X, et al. Lnc-PCDH9-13:1 is a hypersensitive and specific biomarker for early hepatocellular carcinoma. EBioMedicine. 2018;33:57–67.
Lee YR, Kim G, Tak WY, Jang SY, Kweon YO, Park JG, et al. Circulating exosomal noncoding RNAs as prognostic biomarkers in human hepatocellular carcinoma. Int J Cancer. 2019;144(6):1444–52.
Huang J, Zheng Y, Xiao X, Liu C, Lin J, Zheng S, et al. A circulating long noncoding RNA panel serves as a diagnostic marker for Hepatocellular Carcinoma. Dis Markers. 2020;2020:5417598.
Bu X, Ma L, Liu S, Wen D, Kan A, Xu Y, et al. A novel qualitative signature based on lncRNA pairs for prognosis prediction in hepatocellular carcinoma. Cancer Cell Int. 2022;22(1):95.
Yan Y, Ren L, Liu Y, Liu L. Development and Validation of Genome Instability-Associated lncRNAs to predict prognosis and immunotherapy of patients with Hepatocellular Carcinoma. Front Genet. 2021;12:763281.
Chen GY, Wang D. Prognostic Value of Macrophage-Associated Long non-coding RNA expression for Hepatocellular Carcinoma. Cancer Manag Res. 2022;14:215–24.
Wu ZH, Li ZW, Yang DL, Liu J. Development and validation of a pyroptosis-related long non-coding RNA signature for Hepatocellular Carcinoma. Front Cell Dev Biol. 2021;9:713925.
Tao H, Zhang Y, Yuan T, Li J, Liu J, Xiong Y, et al. Identification of an EMT-related lncRNA signature and LINC01116 as an immune-related oncogene in hepatocellular carcinoma. Aging. 2022;14(3):1473–91.
Huang S, Zhang J, Lai X, Zhuang L, Wu J. Identification of Novel Tumor Microenvironment-Related long noncoding RNAs to determine the prognosis and response to Immunotherapy of Hepatocellular Carcinoma Patients. Front Mol Biosci. 2021;8:781307.
Wang L, Ge X, Zhang Z, Ye Y, Zhou Z, Li M, et al. Identification of a ferroptosis-related long noncoding RNA prognostic signature and its predictive ability to Immunotherapy in Hepatocellular Carcinoma. Front Genet. 2021;12:682082.
Jin C, Li R, Deng T, Li J, Yang Y, Li H, et al. Identification and validation of a Prognostic Prediction Model of m6A Regulator-Related LncRNAs in Hepatocellular Carcinoma. Front Mol Biosci. 2021;8:784553.
Deng Y, Zhang F, Sun ZG, Wang S. Development and validation of a prognostic signature Associated with Tumor Microenvironment based on autophagy-related lncRNA analysis in Hepatocellular Carcinoma. Front Med (Lausanne). 2021;8:762570.
Tang P, Qu W, Wang T, Liu M, Wu D, Tan L, et al. Identifying a hypoxia-related long non-coding RNAs signature to improve the prediction of prognosis and immunotherapy response in Hepatocellular Carcinoma. Front Genet. 2021;12:785185.
Bai Y, Lin H, Chen J, Wu Y, Yu S. Identification of Prognostic Glycolysis-Related lncRNA signature in Tumor Immune Microenvironment of Hepatocellular Carcinoma. Front Mol Biosci. 2021;8:645084.
Tabrizian P, Jibara G, Shrager B, Schwartz M, Roayaie S. Recurrence of hepatocellular cancer after resection: patterns, treatments, and prognosis. Ann Surg. 2015;261(5):947–55.
Wang MD, Li C, Liang L, Xing H, Sun LY, Quan B, et al. Early and late recurrence of Hepatitis B Virus-Associated Hepatocellular Carcinoma. Oncologist. 2020;25(10):e1541–e51.
Cheng Z, Yang P, Qu S, Zhou J, Yang J, Yang X, et al. Risk factors and management for early and late intrahepatic recurrence of solitary hepatocellular carcinoma after curative resection. HPB (Oxford). 2015;17(5):422–7.
Kuhn M. Caret: classification and regression training.Astrophysics Source Code Library. 2015:ascl: 1505.003.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.
McCarthy DJ, Chen Y, Smyth GK. Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation. Nucleic Acids Res. 2012;40(10):4288–97.
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47.
Chen H, Boutros MP. Package ‘VennDiagram’. Generate High-Resolution Venn and Euler Plots, Version. 2018;1:20.
Lin H, Zelterman D. Modeling survival data: extending the Cox model. Taylor & Francis; 2002.
Kaplan EL, Meier P. Nonparametric estimation from incomplete observations. J Am Stat Assoc. 1958;53(282):457–81.
Cox DR. Regression models and life-tables. J Roy Stat Soc: Ser B (Methodol). 1972;34(2):187–202.
Tibshirani R. Regression shrinkage and selection via the lasso. J Roy Stat Soc: Ser B (Methodol). 1996;58(1):267–88.
Liaw A, Wiener M. Classification and regression by randomForest. R news. 2002;2(3):18–22.
Guyon I, Weston J, Barnhill S, Vapnik V. Gene selection for cancer classification using support vector machines. Mach Learn. 2002;46(1):389–422.
Friedman J, Hastie T, Tibshirani R. Regularization Paths for generalized Linear Models via Coordinate Descent. J Stat Softw. 2010;33(1):1–22.
Meyer D, Dimitriadou E, Hornik K, Weingessel A, Leisch F, Chang C et al. Misc Functions of the Department of Statistics, Probability Theory Group (Formerly: E1071), TU Wien.(2020). Avalable at: https://cranr-project.org/web/packages/e1071.
Robin X, Turck N, Hainard A, Tiberti N, Lisacek F, Sanchez JC, et al. pROC: an open-source package for R and S + to analyze and compare ROC curves. BMC Bioinformatics. 2011;12:77.
Kassambara A, Kosinski M, Biecek P, Fabian S, survminer. Drawing Survival Curves using’ggplot2’. R package version 0.4. 3. Google Scholar. 2018.
Marshall R. regplot: Enhanced regression nomogram plot. R package version 10. 2020.
Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.
Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102(43):15545–50.
Liberzon A, Birger C, Thorvaldsdottir H, Ghandi M, Mesirov JP, Tamayo P. The Molecular Signatures database (MSigDB) hallmark gene set collection. Cell Syst. 2015;1(6):417–25.
Hanzelmann S, Castelo R, Guinney J. GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinformatics. 2013;14:7.
Jiang P, Gu S, Pan D, Fu J, Sahu A, Hu X, et al. Signatures of T cell dysfunction and exclusion predict cancer immunotherapy response. Nat Med. 2018;24(10):1550–8.
Fu J, Li K, Zhang W, Wan C, Zhang J, Jiang P, et al. Large-scale public data reuse to model immunotherapy response and resistance. Genome Med. 2020;12(1):21.
Reich M, Liefeld T, Gould J, Lerner J, Tamayo P, Mesirov JP. GenePattern 2.0. Nat Genet. 2006;38(5):500–1.
Lu X, Jiang L, Zhang L, Zhu Y, Hu W, Wang J, et al. Immune signature-based subtypes of cervical squamous cell Carcinoma tightly Associated with Human Papillomavirus Type 16 expression, molecular features, and clinical outcome. Neoplasia. 2019;21(6):591–601.
Geeleher P, Cox N, Huang RS. pRRophetic: an R package for prediction of clinical chemotherapeutic response from tumor gene expression levels. PLoS ONE. 2014;9(9):e107468.
Yang W, Soares J, Greninger P, Edelman EJ, Lightfoot H, Forbes S, et al. Genomics of Drug Sensitivity in Cancer (GDSC): a resource for therapeutic biomarker discovery in cancer cells. Nucleic Acids Res. 2013;41(Database issue):D955–61.
Therneau TM. A Package for Survival Analysis in R. R package version 3.1–12 ed2020.
Kole C, Charalampakis N, Tsakatikas S, Vailas M, Moris D, Gkotsis E et al. Immunotherapy for Hepatocellular Carcinoma: A 2021 Update.Cancers (Basel). 2020;12(10).
Stark R, Grzelak M, Hadfield J. RNA sequencing: the teenage years. Nat Rev Genet. 2019;20(11):631–56.
Law CW, Chen Y, Shi W, Smyth GK. Voom: Precision weights unlock linear model analysis tools for RNA-seq read counts. Genome Biol. 2014;15(2):R29.
Clark TG, Bradburn MJ, Love SB, Altman DG. Survival analysis part I: basic concepts and first analyses. Br J Cancer. 2003;89(2):232–8.
Cortes C, Vapnik V. Support-vector networks. Mach Learn. 1995;20(3):273–97.
Zhang R, Ye J, Huang H, Du X. Mining featured biomarkers associated with vascular invasion in HCC by bioinformatics analysis with TCGA RNA sequencing data. Biomed Pharmacother. 2019;118:109274.
Shen J, Qi L, Zou Z, Du J, Kong W, Zhao L, et al. Identification of a novel gene signature for the prediction of recurrence in HCC patients by machine learning of genome-wide databases. Sci Rep. 2020;10(1):4435.
Xiao Y, Huang W, Zhang L, Wang H. Identification of glycolysis genes signature for predicting prognosis in malignant pleural mesothelioma by bioinformatics and machine learning. Front Endocrinol (Lausanne). 2022;13:1056152.
Xiao Y, Zhang G, Wang L, Liang M. Exploration and validation of a combined immune and metabolism gene signature for prognosis prediction of colorectal cancer. Front Endocrinol (Lausanne). 2022;13:1069528.
Mahendran N, Durai Raj Vincent PM, Srinivasan K, Chang CY. Machine learning based computational gene selection models: a Survey, performance evaluation, Open Issues, and future research directions. Front Genet. 2020;11:603808.
Saeys Y, Inza I, Larranaga P. A review of feature selection techniques in bioinformatics. Bioinformatics. 2007;23(19):2507–17.
Fu Y, Wei X, Han Q, Le J, Ma Y, Lin X, et al. Identification and characterization of a 25-lncRNA prognostic signature for early recurrence in hepatocellular carcinoma. BMC Cancer. 2021;21(1):1165.
Gu J, Zhang X, Miao R, Ma X, Xiang X, Fu Y, et al. A three-long non-coding RNA-expression-based risk score system can better predict both overall and recurrence-free survival in patients with small hepatocellular carcinoma. Aging. 2018;10(7):1627–39.
Zhang Q, Ning G, Jiang H, Huang Y, Piao J, Chen Z, et al. 15-lncRNA-Based classifier-clinicopathologic Nomogram improves the prediction of recurrence in patients with Hepatocellular Carcinoma. Dis Markers. 2020;2020:9180732.
Lv Y, Wei W, Huang Z, Chen Z, Fang Y, Pan L, et al. Long non-coding RNA expression profile can predict early recurrence in hepatocellular carcinoma after curative resection. Hepatol Res. 2018;48(13):1140–8.
Gu JX, Zhang X, Miao RC, Xiang XH, Fu YN, Zhang JY, et al. Six-long non-coding RNA signature predicts recurrence-free survival in hepatocellular carcinoma. World J Gastroenterol. 2019;25(2):220–32.
Wang XX, Wu LH, Ai L, Pan W, Ren JY, Zhang Q, et al. Construction of an HCC recurrence model basedon the investigation of immune-relatedlncRNAs and related mechanisms. Mol Ther Nucleic Acids. 2021;26:1387–400.
Ma HY, Yamamoto G, Xu J, Liu X, Karin D, Kim JY, et al. IL-17 signaling in steatotic hepatocytes and macrophages promotes hepatocellular carcinoma in alcohol-related liver disease. J Hepatol. 2020;72(5):946–59.
Gu FM, Li QL, Gao Q, Jiang JH, Zhu K, Huang XY, et al. IL-17 induces AKT-dependent IL-6/JAK2/STAT3 activation and tumor progression in hepatocellular carcinoma. Mol Cancer. 2011;10:150.
Li M, Zhang X, Lu Y, Meng S, Quan H, Hou P, et al. The nuclear translocation of transketolase inhibits the farnesoid receptor expression by promoting the binding of HDAC3 to FXR promoter in hepatocellular carcinoma cell lines. Cell Death Dis. 2020;11(1):31.
Kowalik MA, Columbano A, Perra A. Emerging role of the Pentose phosphate pathway in Hepatocellular Carcinoma. Front Oncol. 2017;7:87.
Pages F, Galon J, Dieu-Nosjean MC, Tartour E, Sautes-Fridman C, Fridman WH. Immune infiltration in human tumors: a prognostic factor that should not be ignored. Oncogene. 2010;29(8):1093–102.
Zheng L, Zhang Z. Decoding the genetic basis of anti-tumor immunity. Immunity. 2021;54(2):199–201.
Galon J, Bruni D. Approaches to treat immune hot, altered and cold tumours with combination immunotherapies. Nat Rev Drug Discov. 2019;18(3):197–218.
Lu Q, Meng Q, Qi M, Li F, Liu B. Shear-Sensitive lncRNA AF131217.1 inhibits inflammation in HUVECs via Regulation of KLF4. Hypertension. 2019;73(5):e25–e34.
Piao HY, Guo S, Jin H, Wang Y, Zhang J. LINC00184 involved in the regulatory network of ANGPT2 via ceRNA mediated miR-145 inhibition in gastric cancer. J Cancer. 2021;12(8):2336–50.
Li G, Wang X, Luo Q, Gan C. Identification of key genes and long noncoding RNAs in celecoxibtreated lung squamous cell carcinoma cell line by RNAsequencing. Mol Med Rep. 2018;17(5):6456–64.
Ye IC, Fertig EJ, DiGiacomo JW, Considine M, Godet I, Gilkes DM. Molecular portrait of hypoxia in breast Cancer: a prognostic signature and novel HIF-Regulated genes. Mol Cancer Res. 2018;16(12):1889–901.
Zhang J, Lou W. A key mRNA-miRNA-lncRNA competing endogenous RNA triple sub-network linked to diagnosis and prognosis of Hepatocellular Carcinoma. Front Oncol. 2020;10:340.
Wei H, Xu Z, Chen L, Wei Q, Huang Z, Liu G, et al. Long non-coding RNA PAARH promotes hepatocellular carcinoma progression and angiogenesis via upregulating HOTTIP and activating HIF-1alpha/VEGF signaling. Cell Death Dis. 2022;13(2):102.
Chen C, Su N, Li G, Shen Y, Duan X. Long non-coding RNA TMCC1-AS1 predicts poor prognosis and accelerates epithelial-mesenchymal transition in liver cancer. Oncol Lett. 2021;22(5):773.
Cui H, Zhang Y, Zhang Q, Chen W, Zhao H, Liang J. A comprehensive genome-wide analysis of long noncoding RNA expression profile in hepatocellular carcinoma. Cancer Med. 2017;6(12):2932–41.
Zhao QJ, Zhang J, Xu L, Liu FF. Identification of a five-long non-coding RNA signature to improve the prognosis prediction for patients with hepatocellular carcinoma. World J Gastroenterol. 2018;24(30):3426–39.
Deng X, Bi Q, Chen S, Chen X, Li S, Zhong Z, et al. Identification of a five-autophagy-Related-lncRNA signature as a Novel Prognostic Biomarker for Hepatocellular Carcinoma. Front Mol Biosci. 2020;7:611626.
Acknowledgements
Not applicable.
Funding
This work is supported by grants from the National Natural Science Foundation of China (31870905 to H.W.), the Scientific Program of Shanghai Municipal Health Commission (201940352 to H.W.), and the Science and Technology Commission of Shanghai Municipality (22ZR1428100 to H.W.), 2020 “Shanghai University Young Teacher Training Funding Program” (A3-2601-20-201001-20, to Y.F.), the Hundred Teacher Talent Program (B3-0200-20-311008-23, to Y.F.) and the University-level Scientific Fund (E3-0200-21-201011-42, to Y.F.) of Shanghai University of Medicine and Health Sciences, the Zhejiang Province Major Science and Technology Project for Medicine and Health (No.WKJ-ZJ-2329, to S.L), State Key Project on Infectious Diseases of China (2018ZX10732202-003), the National Key Research and Development Program of China (2020YFA0909000), The National Natural Science Foundation of China (82127807) and Shanghai Key Laboratory of Molecular Imaging (18DZ2260400).
Author information
Authors and Affiliations
Author notes
Yi Fu, Anfeng Si and Xindong Wei equally contributed.
Contributions
Y.F. and H.W. developed the initial concept. Y.F., S.L., Y.S and H.W. designed the experiment. Y.F. performed bioinformatics analyses. A.S., X.L., X.W., and H.Q. performed the experiment and analyzed data. Y.F., A.S., X.L., X.W., Y.M., H.Q., Z.G., Y.P., Y.Z., X.K., S.L., Y.S and H.W. performed extensive literature search and discussion. Y.F. wrote the manuscript. S.L., Y.S and H.W. gave critical thoughts. H.W. revised the manuscript. All authors discussed the results and approved the final manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
The experimental protocol (81YY-KYLL-19-05) was established, according to the ethical guidelines of the Helsinki Declaration and was approved by the Human Ethics Committee of Jinling Hospital and Shanghai University of Medicine & Health Sciences. Written informed consent was obtained from individual or guardian participants.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Additional file 1. Figure S1.
Flowchart of our analysis strategy. Figure S2. Analyses of differentially expressed lncRNAs between the training group (N = 157) and normal (N = 50). Figure S3. HCC early recurrence analyses.
Additional file 2. Table S1.
Clinical characteristic of 314 TCGA-LIHC patients. Table S2. Differentially expressed lncRNAs associated with DFS. Table S3. Primers used in RT qPCR. Table S4. Drugs respond differently in the low- and high-risk group HCC patients. Table S5. Clinical characteristic of 24 HCC patients from Jinling cohort.
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.
About this article
Cite this article
Fu, Y., Si, A., Wei, X. et al. Combining a machine-learning derived 4-lncRNA signature with AFP and TNM stages in predicting early recurrence of hepatocellular carcinoma. BMC Genomics 24, 89 (2023). https://doi.org/10.1186/s12864-023-09194-8
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12864-023-09194-8