- Open Access
Identification and validation of an eight-lncRNA signature that predicts prognosis in patients with esophageal squamous cell carcinoma
Cellular & Molecular Biology Letters volume 27, Article number: 39 (2022)
Esophageal squamous cell carcinoma (ESCC) is correlated with worse clinical prognosis and lacks available targeted therapy. Thus, identification of reliable biomarkers is required for the diagnosis and treatment of ESCC.
We downloaded the GSE53625 dataset as a training dataset to screen differentially expressed RNAs (DERs) with the criterion of false discovery rate (FDR) < 0.05 and |log2fold change (FC)| > 1. A support vector machine classifier was used to find the optimal feature gene set that could conclusively distinguish different samples. An eight-lncRNA signature was identified by random survival forest algorithm and multivariate Cox regression analysis. The RNA sequencing data from The Cancer Genome Atlas (TCGA) database were used for external validation. The predictive value of the signature was assessed using Kaplan–Meier test, time-dependent receiver operating characteristic (ROC) curves, and dynamic area under the curve (AUC). Furthermore, a nomogram to predict patients’ 3-year and 5-year prognosis was constructed. CCK-8 assay, flow cytometry, and transwell assay were conducted in ESCC cells.
A total of 1136 DERs, including 689 downregulated mRNAs, 318 upregulated mRNAs, 74 downregulated lncRNAs and 55 upregulated lncRNAs, were obtained in the GES53625 dataset. From the training dataset, we identified an eight-lncRNA signature, (ADAMTS9-AS1, DLX6-AS1, LINC00470, LINC00520, LINC01497, LINC01749, MAMDC2-AS1, and SSTR5-AS1). A nomogram based on the eight-lncRNA signature, age, and pathologic stage was developed and showed good accuracy for predicting 3-year and 5-year survival probability of patients with ESCC. Functionally, knockdown of LINC00470 significantly suppressed cell proliferation, G1/S transition, and migration in two ESCC cell lines (EC9706 and TE-9). Moreover, knockdown of LINC00470 downregulated the protein levels of PCNA, CDK4, and N-cadherin, while upregulating E-cadherin protein level in EC9706 and TE-9 cells.
Our eight-lncRNA signature and nomogram can provide theoretical guidance for further research on the molecular mechanism of ESCC and the screening of molecular markers.
Esophageal cancer (EC) is the seventh most common type of malignancy , which is histologically divided into two subtypes: esophageal squamous cell carcinoma (ESCC) and esophageal adenocarcinoma (EAC) . Accounting for > 90% of EC cancers, ESCC is the main EC histologic type, particularly in high-incidence areas of Asia and Africa [2, 3]. Recently, major progress has been made in diagnostic and medical management, especially surgical techniques, chemotherapy, and radiotherapy. Unfortunately, most patients with ESCC have suffered extremely poor outcome mainly due to being diagnosed at advanced stage[4, 5]. Hence, there is an urgent need for identification of reliable biomarkers and targets associated with the prognosis of ESCC.
Nowadays, long noncoding RNAs (lncRNAs) are defined as a class of non-protein-coding RNA transcripts larger than 200 nucleotides in length , which have important regulatory roles in multiple biological processes, including cell differentiation, proliferation, glucose metabolism, and immune response [7, 8]. Aberrantly expressed lncRNAs have contributed to the progression of ESCC pathogenesis from the view of prognosis and cellular functions. For example, upregulation of LINC01296 was associated with poor prognosis and promoted cell proliferation and migration in ESCC . Gao et al.  highlighted the pivotal role of lncRNA CASC9 as a novel diagnostic, prognostic biomarker, and a potential therapeutic target of ESCC. Similarly, LOC100133669 was upregulated in ESCC tissues, and high LOC100133669 expression was associated with poor prognosis of patients with ESCC . Nevertheless, our knowledge on the prognostic role of lncRNAs in ESCC is far from sufficient. Currently, the advancement of high-throughput microarray platforms has helped us perform comprehensive and systemic analysis of lncRNA profiling analysis in ESCC prognosis.
Two major online databases have provided comprehensive cancer genomic datasets: Gene Expression Omnibus (GEO; http://www.ncbi.nlm.nih.gov/geo/) database, a comprehensive library of gene expression in the National Center of Biotechnology Information (NCBI) , and The Cancer Genome Atlas (TCGA, https://gdc-portal.nci.nih.gov/), launched in 2006 by the National Cancer Institute (NCI) and the National Human Genome Research Institute (NHGRI), which contains RNA sequencing (RNA-seq) data and is the database with the most large-scale sequencing results . The methods of mining these two databases mainly focus on the screening of differentially expressed RNAs (DERs) and the analysis of gene regulation networks.
Considering the updated gene expression data and related prognostic information in GEO and TCGA databases, we downloaded lncRNA data, screened DERs, constructed support vector machine (SVM) classifier, and established and validated a risk prediction model for survival prognosis. In addition, we validated the roles of the target gene in vitro.
Materials and methods
The gene expression profile GSE53625 , including 179 ESCC tumor samples and matched controls, was downloaded from Gene Expression Omnibus (GEO: http://www.ncbi.nlm.nih.gov/geo/) database  under the GPL18109 platform (Agilent human lncRNA + mRNA array V.2.0). These 179 samples from GEO were used as a training set. Meanwhile, the data of RNA-seq expression, including 161 tumor tissue samples (80 squamous carcinoma and 81 adenocarcinoma) and 11 controls (platform: Illumina HiSeq 2000 RNA Sequencing), were obtained from the TCGA database. We kept 80 squamous carcinoma sample as the validation set. Statistical clinical information of patients in the training set and validation set is summarized in Table 1.
Identification of significantly DERs
Differential expression analyses were performed for the identification of differentially expressed RNAs (DERs), including lncRNAs and mRNAs (hereafter referred to as “DElncRNAs” and “DEmRNAs,” respectively) between 179 tumor samples and 179 control samples using Limma package version 3.34.7 in R3.4.1 language . The same cutoff value (FDR < 0.05 and |log2FC|) was taken as the inclusion criteria for selection of DElncRNAs and DEmRNAs. According to the value of DERs in training set, pheatmap version 1.0.8 in R3.4.1 language  based on centered Pearson correlation algorithm  was utilized to perform bidirectional hierarchical clustering for describing the gene expression differences between tumor samples and control samples.
Construction and evaluation of SVM classifier
Combined with survival information in training set, we performed univariate Cox regression analysis from survival package version 2.41–1 in R3.4.1 language  to screen significantly prognostic-related DERs (PDERs, including PDElncRNAs and PDEmRNAs) with log-rank p-value < 0.05 as the cutoff criterion. The screened PDElncRNAs were used to conduct recursive feature elimination (RFE) analysis in caret package in R3.4.1 language [20, 21] to extract the optimal feature genes with the minimum root mean square error (RMSE) obtained by the 100-fold cross-validation. Subsequently, these optimal feature genes were applied to construct Sigmoid kernel support vector machine (SVM) model using the R3.4.1 e1071 package (https://cran.r-project.org/web/packages/e1071) . We then evaluated the model’s performance in GSE53625 training set and TCGA validation set using area under the curve (AUC) in receiver operating characteristic (ROC) curve. Meanwhile, we calculated each index value of the ROC curve, including sensitivity, specificity, positive prediction value (PPV), and negative prediction value (NPV).
Identification of signature lncRNAs and RS calculation
On the basis of the optimal feature genes, signature lncRNAs correlated with independent prognosis were identified using a multivariable Cox proportional hazards model implemented with the R3.4.1 survival package version 2.41–1  with log-rank p-value < 0.05 as the cutoff criterion. Then, we calculated risk score (RS) following the risk formula: ∑βlncRNA × ExplncRNA, where βlncRNA indicates the coefficient and ExplncRNA indicates the expression level of signature lncRNA. Afterwards, all patients in training set and validation set were divided into high-risk and low-risk groups according to their median risk score. We used the Kaplan–Meier method in R3.4.1 survival package version 2.41–1  to analyze the overall survival of the two groups and verified the prediction value of the model by plotting ROC curves for the training set and validation set.
Independent prognosis analysis and nomogram construction
The prognostic value of clinical variables and the RS calculated based on lncRNA signature in training set was initially assessed in univariate Cox proportional hazards regression analyses. Subsequently, each significantly different variable was further evaluated in a multivariate Cox proportional hazards regression analysis. The log-rank p-value < 0.05 was served as the cutoff criterion. Furthermore, a nomogram to predict patients’ 3-year and 5-year prognosis was constructed using R3.4.1 rms package version 5.1–2 (https://cran.r-project.org/web/packages/rms/index.html) [23, 24].
Prediction analysis of signature lncRNA-related genes and functional enrichment
To evaluate the function of signature lncRNAs, we first identified mRNAs significantly related to the signature lncRNAs via calculating the Pearson correlation coefficient (PCC) between 8 signature lncRNAs and 92 PDEmRNAs in the data from the training set using the cor.test function in R3.4.1 language . After screening the connection pairs with RCC > 0.6, signature lncRNA and PDEmRNAs co-expression network was constructed and visualized using Cytoscape version 3.6.1 . Subsequently, these PDEmRNAs in co-expression network were inputted into David website (https://david.ncifcrf.gov) to perform GO biological process and KEGG pathway enrichment analysis, with p < 0.05 as the cutoff value.
Clinical samples and cell lines
The tissue samples used were collected from the Harbin Medical University Cancer Hospital between September 2018 and October 2019, including 15 ESCC tissues and 15 adjacent tissues, all from surgically removed specimens. The study was approved by the ethics committee of the Harbin Medical University Cancer Hospital, and each patient signed a written informed consent form.
Two ESCC cell lines (EC9706 and TE-9) were purchased from the Cell Bank of Type Culture Collection of Chinese Academy of Sciences (Shanghai, China), which were cultured in DMEM with 10% FBS (Gibco, USA) at 37 °C containing 5% CO2.
For gene knockdown, EC9706 and TE-9 cells were seeded into six-well plates at a density of 3 × 105 cells per well to 80% confluence and transfected with small interfering RNA targeting LINC00470 (si-LINC00470) or negative control (si-NC) generated by GenePharma (Shanghai, China) in accordance with the instructions of Lipofectamine 3000 Reagents (Invitrogen, USA). After 48 h, cells were harvested for further analysis.
Quantitative real-time PCR analysis
Total RNA was extracted from tissues and cells using TRIzol reagent (TakaRa, Dalian, China), and reverse transcription was performed with PrimeScript RT Reagent Kit with gDNA Eraser (TakaRa, Dalian, China). Quantitative real-time PCR analysis was conducted on LightCycler 480 II Real-Time PCR System (Roche, Basel, Switzerland) using SYBR Premix Ex Taq II (TakaRa). The primers used in our study were as follows: LINC00470, forward 5′-CGTAAGGTGACGAGGAGCTG-3′ and reverse 5′-GGGGAATGGCTTTTGGGTCA-3′; GAPDH forward 5′- GTCAACGGATTTGGTCTGTATT-3′ and reverse 5′- AGTCTTCTGGGTGGCAGTGAT-3′. The relative expression level LINC00470 was calculated using 2−ΔΔCT method and normalized to GAPDH.
Cell proliferation assay
CCK-8 assay was performed to evaluate the cell proliferation ability in ESCC cells. In brief, transfected cells were inoculated into 96-well plates at a density of 3000 cells per well. At the indicated timepoint (0, 24, 48, and 72 h, respectively), 10 µl of CCK-8 solution (Sigma-Aldrich, USA) was added to each well. After 2 h incubation, the absorbance in each well was measured at 450 nm under a microplate reader.
The cell cycle distribution was analyzed using flow cytometry. Briefly, transfected cells (1 × 106) were harvested, washed with PBS, and fixed by ice-cold ethanol (70%) overnight at 4 °C. Afterwards, cells were washed with PBS twice and stained with propidium iodide (PI) for 30 min at 37 °C. The DNA content of stained cells was determined using BD FACSCalibur flow cytometer (BD Biosciences, Franklin Lakes, NJ, USA) and analyzed with ModFitLT.
Cell migration assay
Cell migration was measured using transwell 24-well chambers (Corning Inc, Corning, NY, USA). In brief, transfected cells (5 × 105) were harvested and resuspended in serum-free medium. Then, the cell suspensions were added to the upper chamber, and 600 µl medium containing 15% FBS was added to the lower chamber. After 12 h culture, the migratory cells in the lower chamber were fixed with 4% paraformaldehyde for 10 min and stained in 0.5% crystal violet (Sigma-Aldrich, USA) for 30 min. Finally, migratory cells were photographed and counted from five random fields under a light microscope.
Western blot analysis
Total protein sample was extracted from cell lines with RIPA lysis buffer (Beyotime Institute of Biotechnology, Shanghai, China). Proteins of equal amounts (30 μg) were separated by 10% SDS-PAGE and transferred to PVDF membranes (Millipore). After blocking with 5% nonfat milk, the membranes were incubated with primary antibodies against PCNA (1:1000, ab18197, Abcam), CDK4 (1:1000, ab226474, Abcam), E-cadherin (1:1000, ab219332, Abcam), N-cadherin (1:1000, ab76059, Abcam), and GAPDH (1:5,000; ab8245; Abcam) overnight at 4 °C. After an incubation with horseradish-peroxidase-conjugated secondary antibody (1:5000, SC-2005, Santa Cruz, Inc.) for 2 h at room temperature, the protein bands were visualized with the enhanced chemiluminescence (ECL) Plus kit (Beyotime Institute of Biotechnology).
All quantitative data were analyzed using GraphPad Prism 5 (La Jolla, CA, USA) and expressed as mean ± standard deviation (SD). Differences between si-NC and si-LINC00470 groups were assessed using Student’s t-test. A p-value of < 0.05 was considered statistically significant.
Identification of significantly DERs
Significant DERs were first identified among 179 tumor samples compared with 179 control samples in the training set. A total of 129 DElncRNAs (74 downregulated and 55 upregulated) and 1007 DEmRNAs (689 downregulated and 318 upregulated) were identified and are listed in Additional file 1: Table S1. These data were used to build the volcano plot of DElncRNAs and DEmRNAs (Fig. 1A) and the bidirectional hierarchical clustering heatmap (Fig. 1B), indicating the samples tend to cluster in two distinct directions.
Optimal feature gene selection
A total of 114 PDERs, including 22 PDElncRNAs and 92 PDEmRNAs, were obtained after univariate Cox regression analysis and are listed in Additional file 2: Table S2. Based on the screened 22 PDElncRNAs, the lncRNA combination with the lowest RMSE was selected as the optimal feature genes in the RFE recursive algorithm screening. As shown in Fig. 2, when the number of lncRNAs was 13, the optimal parameter (minimum RMSE = 0.1352) was obtained, and corresponding 13 optimal feature genes are summarized in Additional file 3: Table S3. A classification model was constructed in training set, whose performance was assessed in the GSE53625 training set and TCGA validation set. The classification results of samples based on the classifier are shown in the scatter diagram in Fig. 3 (left), in which the points with two different colors and shapes are clearly distinguished. The area under the ROC curve is shown in Fig. 3 (right), and corresponding index values of the ROC curve are presented in Table 2. ROC curve analysis revealed an AUC of 0.997 in the training set and 0.901 in the validation set. These results indicate that these optimal feature genes could be used as effective and accurate ESCC diagnostic biomarkers.
Identification and validation of an eight-signature lncRNAs
Multivariate Cox regression analysis was used to develop signature lncRNAs that are independent predictors of the optimal feature genes in the SVM model. An eight-lncRNA signature was identified, including ADAMTS9-AS1, DLX6-AS1, LINC00470, LINC00520, LINC01497, LINC01749, MAMDC2-AS1, and SSTR5-AS1. The risk coefficients suggested that ADAMTS9-AS1, LINC01497, and MAMDC2-AS1 were risk factors for ESCC (coef > 0), whereas DLX6-AS1, LINC00470, LINC00520, LINC01749, and SSTR5-AS1 appeared to be protective factors (coef < 0) (Table 3). The RS of each patient in the training set and validation set was calculated with the following formula: RS = (0.147172) × ExpADAMTS9-AS1 + (−0.063991) × ExpDLX6-AS1 + (−0.112843) × ExpLINC00470 + (−0.065239) × ExpLINC00520 + (0.184709) × ExpLINC01497 + (−0.166036) × ExpLINC01749 + (0.104274) × ExpMAMDC2-AS1 + (−0.163769) × ExpSSTR5-AS1. The higher the risk score, the worse the clinical prognosis. Accordingly, patients were divided into high- and low-risk groups depending on their median risk score to assess the score’s ability to accurately predict survival in a Cox regression model (Additional file 4: Table S4). Kaplan–Meier analysis showed that patients in the low-risk group had better prognosis than those in the high-risk group in the training set (Fig. 4A) and validation set (Fig. 4B). The AUC of the ROC curve was 0.989 in the training set and 0.865 in the validation set (Fig. 4C). These results confirmed that the risk score could be an independent predictor of overall survival.
The eight-lncRNA signature was an independent predictor of ESCC prognosis
To investigate whether the eight-lncRNA signature was an independent predictor of prognosis among patients with ESCC in the training set, we performed univariate and multivariate Cox regression analyses. As illustrated in Table 4, the age, pathologic N, pathologic stage, adjuvant therapy, and RS model status were significantly correlated with patients’ overall survival in the univariate Cox regression. Moreover, the age, pathologic stage, and RS model status based on the eight-lncRNA signature remained three independent predictors. In addition, the results from Kaplan–Meier analysis showed that age (Fig. 5A) and pathologic stage (Fig. 5B) had a significant impact on the prognosis of patients with ESCC (with a log-rank test p-value less than 0.0001). Furthermore, a nomogram was constructed that integrated age, pathologic stage, and RS model status to analyze the relationship between these three predictors and survival prognosis (Fig. 6A), which indicated that a higher total number of points on the nomogram presented a worse prognosis. Further analysis suggested that the predicted 3-year and 5-year survival rates by the survival model in the histogram were consistent with the actual 3-year and 5-year survival rates (Fig. 6B).
Functional characteristics of signature lncRNA-related genes
We first calculated the PCC between expression levels of 92 PDEmRNAs and eight-lncRNA signature and obtained 279 connection pairs with PCC > 0.6 (Additional file 1: Table S5). A total of 82 nodes, including 8 signature lncRNAs and 74 PDEmRNAs, were obtained in the constructed co-expression network (Fig. 7). Then we performed GO and KEGG functional enrichment analysis for these 74 PDEmRNAs. As shown in Fig. 8 and Table 5, these mRNAs were mainly enriched in the differentiation and development of epidermal and epithelial cells in GO biological process analysis, as well as the secretion of digestive juices in KEGG enrichment analysis.
Validation of the expression levels of eight-lncRNA signature in ESCC tissues
Quantitative real-time PCR analysis was performed to determine the expression levels of eight-lncRNA signature in 15 pairs of tumor tissues and matched adjacent tissues derived from patients with ESCC. As shown in Fig. 9, the expression levels of DLX6-AS1 and LINC00470 were significantly upregulated, while LINC01479, LINC01749, and SSTR5-AS1 were markedly downregulated in ESCC tissues compared with adjacent tissues. However, there was no significant differences in expression levels of ADAMTS9-AS1, LINC00520, or MAMDC2-AS1 between two groups. According to the higher fold change, we selected LINC00470 for subsequent functional assays.
Knockdown of LINC00470 suppresses ESCC cell proliferation, G1/S transition, and migration
To investigate the function of LINC00470 in ESCC in vitro, LINC00470 expression was first knocked down in EC9706 and TE-9 cells by using si-LINC00470 transfection, which was demonstrated by quantitative real-time PCR analysis (Fig. 10A). CCK-8 assay showed that knockdown of LINC00470 resulted in growth retardation of EC9706 and TE-9 cells (Fig. 10B). Moreover, the percentage of cells at G0/G1 phase was significantly increased, in accordance with S and G2/M phase being decreased in si-LINC00470 group compared with si-NC group in both EC9706 (Fig. 10C) and TE-9 (Fig. 10D) cells. In addition, transwell assay indicated that knockdown of LINC00470 markedly inhibited the cell migration ability in EC9706 and TE-9 cells (Fig. 10E). At the molecular level, knockdown of LINC00470 downregulated the protein levels of PCNA, CDK4, and N-cadherin, while upregulating E-cadherin protein level in EC9706 and TE-9 cells (Fig. 10F). The above results demonstrate that knockdown of LINC00470 can inhibit the proliferation and migration of ESCC cells.
To the best of our best knowledge, the tumor–node–metastasis (TNM) staging system acts as the main transitional algorithm to direct the treatment strategies and also serves as a prognostic predictor, but fails to consider the genetic alterations in most types of cancers, including ESCC [27, 28]. In recent years, identification of lincRNA-based signatures has received great attention for its potential to aid in the prognosis of cancers, including hepatocellular carcinoma , bladder cancer , and pancreatic cancer .
In the present study, we first identified 1136 significantly DEGs between tumor tissues and normal tissues in GEO data and confirmed 114 DEGs correlated with prognosis. Finally, eight-lncRNA signature (DLX6-AS1, LINC00470, LINC01479, LINC01749, SSTR5-AS1, ADAMTS9-AS1, LINC00520, and MAMDC2-AS1) was constructed for ESCC. Importantly, a robust nomogram consisting of age, pathologic stage, and RS model status based on the eight-lncRNAs signature was constructed for prediction of prognosis for patients with ESCC. Further analysis suggested the predicted 3-year and 5-year survival rates by the survival model in the histogram were consistent with the actual 3- and 5-year survival rates. By integrating diverse prognostic variables based on clinical characteristics, nomogram has been a widely used tool in oncology that could determine individual probability . Here, our data suggest that our constructed nomogram had better predictive accuracy than each factor alone. Similar to our data, Khalil et al.  established a three-lncRNA signature and demonstrated that it could precisely predict overall survival and disease-free survival for ESCC. Three-lncRNA signature (RP11-366H4.1.1, LINC00460, and AC093850.2) was constructed by random forest algorithm and support vector machine algorithm and identified to be potential predictor of overall survival for patients with ESCC . In addition, Mao et al.  identified a robust seven-lncRNA signature associated with overall survival that was independent of classical prognostic factors and molecular subtypes in ESCC. The different lncRNA signatures identified in ESCC might be mainly ascribed to different sample resources, sample sizes, and analysis methods. Subsequently, our data showed that 74 PDEmRNAs in co-expression network were mainly enriched in the differentiation and development of epidermal and epithelial cells, as well as the secretion of digestive juices. Consistently, ESCC progression was closely associated with epidermal and epithelial cell differentiation and growth [35, 36].
Subsequently, we confirmed that the expression levels of DLX6-AS1 and LINC00470 were significantly upregulated, while LINC01479, LINC01749, and SSTR5-AS1 were markedly downregulated in ESCC tissues compared with adjacent tissues. By searching published articles, we found that no review had explored the intriguing mechanisms of these five lncRNAs in ESCC, except DLX6-AS1. Several studies have demonstrated that DLX6-AS1 is associated with malignant progression and promotes cell growth and metastasis in ESCC cells [37,38,39]. Considering the relatively higher increased fold change in expression level, we selected LINC00470 for further functional experiments. As expected, knockdown of LINC00470 significantly suppressed cell proliferation, G1/S transition, and migration in two ESCC cell lines (EC9706 and TE-9). In fact, LINC00470 has been reported to be an oncogene in other malignant tumors. For instance, Wu et al.  reported that LINC00470 promoted glioma cell proliferation and invasion and attenuated chemosensitivity. Yan et al.  performed overexpression and knockdown experiments to demonstrate the oncogenic functions of LINC00470 on gastric cancer cell proliferation, migration, and invasion. The findings by Huang et al.  indicated that knockdown of LINC00470 expression inhibited cell proliferation and cell cycle progression, while overexpression of LINC00470 showed the opposite effects in hepatocellular carcinoma. In addition, LINC00470 promoted invasiveness, migration, and angiogenesis of endometrial cancer cells . Knockdown of LINC00470 could significantly inhibit the melanoma cell proliferation and migration, and suppress the growth of tumor in vivo . On the basis of this evidence, we speculate that high LINC00470 expression appears to be related to poor prognosis in ESCC. It must be mentioned that there are several limitations to this study, including lack of further in vitro experimental study and in vivo data to validate the prognostic performance of our proposed lncRNA signature.
In summary, our findings identified and validated an eight-lincRNA signature and nomogram as reliable prognostic tools for ESCC. These eight hub genes (ADAMTS9-AS1, DLX6-AS1, LINC00470, LINC00520, LINC01497, LINC01749, MAMDC2-AS1, and SSTR5-AS1) may offer novel therapeutic strategies for patients with ESCC.
Availability of data and materials
All datasets generated for this study are included in the manuscript.
Esophageal squamous cell carcinoma
National Center of Biotechnology Information
National Cancer Institute
Differentially expressed genes
Support vector machine
Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424.
Arnold M, Soerjomataram I, Ferlay J, Forman D. Global incidence of oesophageal cancer by histological subtype in 2012. Gut. 2015;64(3):381–7.
Herszenyi L, Tulassay Z. Epidemiology of gastrointestinal and liver tumors. Eur Rev Med Pharmacol Sci. 2010;14(4):249–58.
Wang WL, Chang WL, Yang HB, Wang YC, Chang IW, Lee CT, et al. Low disabled-2 expression promotes tumor progression and determines poor survival and high recurrence of esophageal squamous cell carcinoma. Oncotarget. 2016;7(44):71169–81.
Aquino JL, Said MM, Pereira DA, Cecchino GN, Leandro-Merhi VA. Complications of the rescue esophagectomy in advanced esophageal cancer. Arq Bras Cir Dig. 2013;26(3):173–8.
Quinn JJ, Chang HY. Unique features of long non-coding RNA biogenesis and function. Nat Rev Genet. 2016;17(1):47–62.
Rinn JL, Chang HY. Genome regulation by long noncoding RNAs. Annu Rev Biochem. 2012;81:145–66.
Nie L, Wu HJ, Hsu JM, Chang SS, Labaff AM, Li CW, et al. Long non-coding RNAs: versatile master regulators of gene expression and crucial players in cancer. Am J Transl Res. 2012;4(2):127–50.
Wang B, Liang T, Li J. Long noncoding RNA LINC01296 is associated with poor prognosis in ESCC and promotes ESCC cell proliferation, migration and invasion. Eur Rev Med Pharmacol Sci. 2018;22(14):4524–31.
Gao GD, Liu XY, Lin Y, Liu HF, Zhang GJ. LncRNA CASC9 promotes tumorigenesis by affecting EMT and predicts poor prognosis in esophageal squamous cell cancer. Eur Rev Med Pharmacol Sci. 2018;22(2):422–9.
Guan Z, Wang Y, Wang Y, Liu X, Wang Y, Zhang W, et al. Long non-coding RNA LOC100133669 promotes cell proliferation in oesophageal squamous cell carcinoma. Cell Prolif. 2020;53(4):e12750.
Edgar R, Domrachev M, Lash AE. Gene expression omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002;30(1):207–10.
Weinstein JN, Collisson EA, Mills GB, Shaw KR, Ozenberger BA, Cancer Genome Atlas Research N, et al. The cancer genome atlas pan-cancer analysis project. Nat Genet. 2013;45(10):1113–20.
Li J, Chen Z, Tian L, Zhou C, He MY, Gao Y, et al. LncRNA profile study reveals a three-lncRNA signature associated with the survival of patients with oesophageal squamous cell carcinoma. Gut. 2014;63(11):1700–10.
Barrett T, Suzek TO, Troup DB, Wilhite SE, Ngau WC, Ledoux P, et al. NCBI GEO: mining millions of expression profiles—database and tools. Nucleic Acids Res. 2005;33:D562–6.
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.
Wang L, Cao C, Ma Q, Zeng Q, Wang H, Cheng Z, et al. RNA-seq analyses of multiple meristems of soybean: novel and alternative transcripts, evolutionary and functional implications. BMC Plant Biol. 2014;17(14):169.
Eisen MB, Spellman PT, Brown PO, Botstein D. Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci USA. 1998;95(25):14863–8.
Wang P, Wang Y, Hang B, Zou X, Mao JH. A novel gene expression-based prognostic scoring system to predict survival in gastric cancer. Oncotarget. 2016;7(34):55343–51.
Lu X, Yang Y, Wu F, Gao M, Xu Y, Zhang Y, et al. Discriminative analysis of schizophrenia using support vector machine and recursive feature elimination on structural MRI images. Medicine. 2016;95(30):e3973.
Deist TM, Dankers F, Valdes G, Wijsman R, Hsu IC, Oberije C, et al. Machine learning algorithms for outcome prediction in (chemo)radiotherapy: an empirical comparison of classifiers. Med Phys. 2018;45(7):3449–59.
Wang Q, Liu X. Screening of feature genes in distinguishing different types of breast cancer using support vector machine. Onco Targets Ther. 2015;8:2311–7.
Anderson WI, Schlafer DH, Vesely KR. Thyroid follicular carcinoma with pulmonary metastases in a beaver (Castor canadensis). J Wildl Dis. 1989;25(4):599–600.
Eng KH, Schiller E, Morrell K. On representing the prognostic value of continuous gene expression biomarkers with the restricted mean survival curve. Oncotarget. 2015;6(34):36308–18.
Zou KH, Tuncali K, Silverman SG. Correlation and simple linear regression. Radiology. 2003;227(3):617–22.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.
Gerlinger M, Rowan AJ, Horswell S, Math M, Larkin J, Endesfelder D, et al. Intratumor heterogeneity and branched evolution revealed by multiregion sequencing. N Engl J Med. 2012;366(10):883–92.
Navin N, Kendall J, Troge J, Andrews P, Rodgers L, McIndoo J, et al. Tumour evolution inferred by single-cell sequencing. Nature. 2011;472(7341):90–4.
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.
He A, He S, Peng D, Zhan Y, Li Y, Chen Z, et al. Prognostic value of long non-coding RNA signatures in bladder cancer. Aging. 2019;11(16):6237–51.
Wu B, Wang K, Fei J, Bao Y, Wang X, Song Z, et al. Novel three-lncRNA signature predicts survival in patients with pancreatic cancer. Oncol Rep. 2018;40(6):3427–37.
Mao Y, Fu Z, Zhang Y, Dong L, Zhang Y, Zhang Q, et al. A seven-lncRNA signature predicts overall survival in esophageal squamous cell carcinoma. Sci Rep. 2018;8(1):8823.
Khalil AM, Guttman M, Huarte M, Garber M, Raj A, Rivea Morales D, et al. Many human large intergenic noncoding RNAs associate with chromatin-modifying complexes and affect gene expression. Proc Natl Acad Sci USA. 2009;106(28):11667–72.
Huang GW, Xue YJ, Wu ZY, Xu XE, Wu JY, Cao HH, et al. A three-lncRNA signature predicts overall survival and disease-free survival in patients with esophageal squamous cell carcinoma. BMC Cancer. 2018;18(1):147.
Cui L, Pan XM, Ma CF, Shang-Guan J, Yu HB, Chen GX, et al. Association between epidermal growth factor polymorphism and esophageal squamous cell carcinoma susceptibility. Dig Dis Sci. 2010;55(1):40–5.
Yoshioka M, Ohashi S, Ida T, Nakai Y, Kikuchi O, Amanuma Y, et al. Distinct effects of EGFR inhibitors on epithelial- and mesenchymal-like esophageal squamous cell carcinoma cells. J Exp Clin Cancer Res. 2017;36(1):101.
Wang M, Li Y, Yang Y, Liu X, Zang M, Li Y, et al. Long non-coding RNA DLX6-AS1 is associated with malignant progression and promotes proliferation and invasion in esophageal squamous cell carcinoma. Mol Med Rep. 2019;19(3):1942–50.
Wu SB, Wang HQ. Upregulation of long noncoding RNA DLX6-AS1 promotes cell growth and metastasis in esophageal squamous cell carcinoma via targeting miR-577. Eur Rev Med Pharmacol Sci. 2020;24(14):7557.
Wu SB, Wang HQ. Upregulation of long noncoding RNA DLX6-AS1 promotes cell growth and metastasis in esophageal squamous cell carcinoma via targeting miR-577. Eur Rev Med Pharmacol Sci. 2020;24(3):1195–201.
Wu C, Su J, Long W, Qin C, Wang X, Xiao K, et al. LINC00470 promotes tumour proliferation and invasion, and attenuates chemosensitivity through the LINC00470/miR-134/Myc/ABCC1 axis in glioma. J Cell Mol Med. 2020;24(20):12094–106.
Yan J, Huang X, Zhang X, Chen Z, Ye C, Xiang W, et al. LncRNA LINC00470 promotes the degradation of PTEN mRNA to facilitate malignant behavior in gastric cancer cells. Biochem Biophys Res Commun. 2020;521(4):887–93.
Huang W, Liu J, Yan J, Huang Z, Zhang X, Mao Y, et al. LncRNA LINC00470 promotes proliferation through association with NF45/NF90 complex in hepatocellular carcinoma. Hum Cell. 2020;33(1):131–9.
Yi T, Song Y, Zuo L, Wang S, Miao J. LINC00470 stimulates methylation of PTEN to facilitate the progression of endometrial cancer by recruiting DNMT3a Through MYC. Front Oncol. 2021;11:646217.
Huang T, Wang YJ, Huang MT, Guo Y, Yang LC, Liu XJ, et al. LINC00470 accelerates the proliferation and metastasis of melanoma through promoting APEX1 expression. Cell Death Dis. 2021;12(5):410.
We would like to thank all participants enrolled in the present study.
Ethics approval and consent to participate
This study was conducted in accordance with the Declaration of Helsinki (1975) and approved by the ethics committee of the Harbin Medical University Cancer Hospital (approval no. HMUC-M54G, 2018.8.23, Heilongjiang Province, China). Each patient signed a written informed consent form.
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Identification of DElncRNAs and DEmRNAs.
Table S2. List of total PDERs after univariate cox regression analysis.
Table S3. List of optial feature genes.
Table S4. Summary of patients in high- and low-risk groups.
Table S5. List of lncRNA signature and corresponding connection pairs.
About this article
Cite this article
Zhang, J., Ling, X., Fang, C. et al. Identification and validation of an eight-lncRNA signature that predicts prognosis in patients with esophageal squamous cell carcinoma. Cell Mol Biol Lett 27, 39 (2022). https://doi.org/10.1186/s11658-022-00331-x
- Esophageal squamous cell carcinoma
- Long noncoding RNA