Survival-Related lncRNA Landscape Analysis Identifies LINC01614 as an Oncogenic lncRNA in Gastric Cancer

Background: Long non-coding RNAs (lncRNAs) reportedly play important roles in biomarker and tumorigenesis of gastric cancer (GC). This study aimed to determine the potential application of prognostic lncRNA signature and identified the role of LINC01614 in carcinogenesis in GC. Material and Methods: Data accessed from the Cancer Genome Atlas database was used to construct a lncRNA signature. Joint effect analysis of the signature and clinical parameters was performed to verify the clinical value of the signature. Co-expression analysis was conducted for prognostic lncRNAs and protein-coding genes. Moreover, the relative expression of LINC01614 was validated in GC tissues and cell lines. In vitro and in vivo experiments were conducted to analyze the biological functions of the newly identified gene in GC cells. Results: A seven-lncRNA (LINC01614, LINC01537, LINC01210, OVAAL, LINC01446, CYMP-AS1, and SCAT8) signature was identified as a promising prognostic signature in GC. Results indicated that the seven-lncRNA was involved in tumorigenesis and progression pathways. LINC01614 expression was identified and found to be upregulated in GC tissues and cells. The study findings revealed that LINC01614 promoted cell proliferation, migration, invasion, and epithelial-mesenchymal transition. Knockdown of LINC01614 arrested cell cycle distribution at the G2/M phase. Further, LINC01614 also promoted tumor growth in vivo. Conclusion: We developed an independent seven-lncRNA biomarker for prognostic prediction and identified LINC01614 as an oncogenic lncRNA in GC.


INTRODUCTION
Gastric cancer (GC) is prevalent worldwide (Siegel et al., 2020). Among the Chinese population, GC is one of the most commonly reported cancers and poses a considerable public health burden (Chen et al., 2018). Despite the remarkable developments in diagnosis and treatment of GC, such as chemotherapy, radiotherapy, and surgical techniques, the overall survival has remained unsatisfactory (Yang et al., 2018). Malignant transformation and tumor progression in GC patients is diverse owing to the heterologous nature of the tumor. Moreover, the lack of availability of reliable diagnostic methods in the early stage leads to poor overall survival (OS) in GC patients (Van Cutsem et al., 2016). Recently, utilization of precision medicine has necessitated disease characterization by molecular biology approaches; the same plays a vital in GC. Replacement of the suboptimal methods with more refined and advanced methods may benefit the patients. All these objective factors contribute to the emerging demand for reliable prognosis biomarkers for GC.
The transcriptome plays an important role in disease progression and is indispensable in the investigation of carcinogenesis and pathophysiology of GC. With the development of next-generation sequencing, long noncoding RNAs (lncRNAs) have received considerable attention in the recent years. After being initially misunderstood and overlooked as "junk" in the transcriptome, lncRNAs have recently been proven to be a significant contributor in tumorigenesis. The lncRNA transcription has been reported to partly provide the signals necessary for malignant transformation (Huarte, 2015). Furthermore, lncRNAs are capable of establishing interactions with the DNA, protein, and RNA to perform regulation in multiple cancer phenotypes (Schmitt and Chang, 2016).
Utilizing the well-established methodology of genome-wide screening analysis, numerous studies have shown that the aberrant expression of lncRNA is associated with the development of various cancer types, including GC, and lncRNAs may serve as promising biomarkers in GC. (Yang et al., 2016;Zhuo et al., 2019) In this regard, systematic research is warranted to identify the function of lncRNA and to construct a promising signature based on the differential expression level of lncRNA. This study aimed to identify a promising lncRNA signature that could be associated with the OS of patients and to obtain a key lncRNA in GC.

RNA Sequencing Data and Annotation
The Cancer Genome Atlas (TCGA) database was accessed to extract data on 407 GC RNA sequencing (RNA-Seq) profiles. Ensemble Genomes annotation was used to define the lncRNA.

Clinical Information
Clinical data on 380 patients corresponding to the 407 RNA-Seq profiles were obtained from TCGA. Among the data acquired, data on 368 GC patients were included in the survival analysis; data on 12 patients were excluded owing to a lack of complete survival information. Patients who met the following criteria were included in the study: 1) diagnosed GC, 2) complete RNA-Seq data, 3) complete prognostic information. The data are publicly accessible and open in TCGA, data acquisition and application in the research complied to TCGA responsible use policy and publication guidelines. Consequently, there is no need for additional agreement by any local ethics committees.

Expression Analysis of Differentially Expressed lncRNAs in Gastric Cancer
Based on the RNA-Seq data, data on lncRNA expression profiles were extracted according to annotation. To efficiently identify the relevant biomarkers in GC, lncRNA data with a mean expression value of less than one was first excluded from the subsequent steps of screening. The edgeR package in the R language was used to screen out the differentially expressed lncRNAs (DELs) (Robinson et al., 2010). The threshold of |log2 fold change (log2FC)|> 1 as well as false discovery rate (FDR) < 0.05 were considered to screen the significant DELs.

Construction and Analysis of the DEL Expression-Based Prognostic Signature
Association between DEL and OS of patients with GC was performed by using a univariate Cox regression model in the "Survival" package. Considering p-value < 0.01, DELs were regarded as statistically significant and considered to be potential candidates of the prognostic DELs. Subsequently, multivariate Cox regression and the "step" function which could choose a formula-based model by Akaike Information Criteria in a stepwise algorithm method in R language were used to assess the optimal combination and to identify the best model of prognostic DEL signature. The prognostic signature model was established based on a linear combination of the prognostic DEL expression levels and the multivariate Cox model. The weight was set according to the coefficient (β) that was used to multiply the regression model values calculated and output. The risk score formula was calculated as shown in the index: The Prognosis Risk Score Index (β 1 * expression level of DEL 1 ) + (β 2 * expression level of DEL 2 ) +. . .+ (β n * expression level of DEL n ). Risk score for each GC patient is calculated based on the index. Median value of the risk score for all patients was set as a cut-off point, which was then used to categorize the patients into high-and low-risk groups. Time-dependent receiver operating characteristic (ROC) curves were generated within one, three, and 5 years. The Kaplan-Meier survival curves for the cases predicted to be high-risk or low-risk were generated. Cox regression analysis and stratified analysis were further performed to demonstrate whether the signature was independent for other clinical features. Moreover, joint effect analysis was applied to explore the combined role of DEL expression-based signature and the clinical factors in prognosis.

Prognostic DELs for Prediction and Enriched Pathway Analysis
Based on the guilt-by-association theory, the Pearson correlation coefficient was used to examine the correlation between prognostic lncRNAs and protein-coding genes (PCGs) (Rinn and Chang, 2012). The PCGs with |Pearson correlation coefficient| > 0.4 and p < 0.05 were considered to be potential lncRNA-related PCGs and subsequently used for conducting enrichment analysis. Gene Ontology (GO) analysis comprising molecular function (MF), cellular component (CC), and Frontiers in Genetics | www.frontiersin.org October 2021 | Volume 12 | Article 698947 biological process (BP) was performed. The Kyoto Encyclopedia of Genes and Genomes (KEGG) is a database resource that provides information on the biological function of the genes and the associated pathways. Using the "clusterProfiler" package in R, PCGs exhibiting co-expression pattern with target DELs were analyzed based on GO and KEGG (Yu et al., 2012). The enriched terms with Benjamini-Hochberg-adjusted p < 0.05 were considered statistically significant.

Clinical Samples
Twenty-two paired GC tissues and adjacent normal tissues (distance to cancer >5 cm) were obtained from patients who underwent GC resection at Guangxi Medical University Cancer Hospital. All selected samples were diagnosed by pathologists and the patients were never subjected to chemoradiotherapy or other treatments before surgery. The Ethics Committee of the Guangxi Medical University Cancer Hospital approved this study. Both informed and written consent was obtained from the participants.

Cell Culture and Transfection
Human GC cell lines (SGC-7901, HGC-27, MGC-803, and AGS) and immortalized human gastric epithelial cell line (GES-1) were purchased from the Cell Bank of the Chinese Academy of Science, Shanghai, China. All cells were incubated in the Dulbecco's modified Eagle medium (DMEM) (GIBCO, Waltham, MA, United States) supplemented with 10% fetal bovine serum (FBS; GIBCO). Cells were incubated at 37°C with humidified 5% carbon dioxide. At 50% cell confluency, silencing small interfering RNAs (siRNAs; GenePharma, Shanghai, China) were transfected into cells by using Lipofectamine 3,000 (Invitrogen, CA, United States). The target sequences of siRNAs are listed in Supplementary  Table S1. LINC01614-NC and LINC01614-overexpressing plasmids were designed and synthesized by GenePharma Co., Ltd. and packaged into lentivirus vectors. The lentivirus vector was then transfected into SGC-7901 cells; the stably transfected SGC-7901-NC and SGC-7901-LINC01614 cells were selected using puromycin.

Quantitative Real-Time Polymerase Chain Reaction (qRT-PCR)
Based on the instructions of the RNAiso reagent (TaKaRa Biotechnology, Dalian, China), total RNA was isolated from the tissues and cells. All samples were quantified by using Nanodrop One (Thermo Fisher Scientific Inc., MA, United States). Primers were synthesized by Tsingke (Beijing, China). Total RNA (20 ng) was subjected to reverse transcription and then used for performing real-time PCR using a qRT-PCR Kit, according to the manufacturer's instructions (TaKaRa Biotechnology). And qTOWER 3G (Analytik jena, jena, German) was used to determine the expression level of LINC01614. The relative expression of the target gene was quantified using the 2 −ΔΔCt method; β-actin was used as the internal reference. The primer sequences are shown in Supplementary Table S1.

Cell Proliferation Assay
Cell proliferation assay was performed using the Cell Counting Kit-8 (CCK-8; Dojindo Laboratory, Tokyo, Japan). Cells were plated in 96-well plates at a density of 1 × 10 3 cells/well. CCK-8 was used to determine the cell proliferation ability, according to the manufacturers' instructions. The absorbance of each well was recorded at 450 nm using a microplate reader (BioTek, VT, United States).

Colony Formation Assay
Cells in the logarithmic phase were isolated by pancreatic enzyme digestion and then incubated in DMEM. Thereafter, the cells were added into 6-well plates at a density of 1 × 10 3 cells/well and cultured in an incubator until visible cell colony formation. Cells were then subjected to washing steps using phosphate-buffered saline (PBS), and fixed and stained with 0.1% crystal violet and methanol. Each experiment was repeated thrice.

Western Blotting
Cells in each group were subjected to washing steps thrice using cold PBS and underwent lysis on ice for 10 min using RIPA (Beyotime), protease inhibitors, and phosphatase inhibitors. Total protein was extracted according to manufacturer's instructions and the Enhanced BCA Protein Assay Kit (Beyotime) was used to determine the protein concentration, according to a selected protocol. Target proteins were isolated using 10% sodium dodecyl sulfate polyacrylamide gel electrophoresis. Subsequently, the proteins were transferred onto a polyvinylidene fluoride membrane (Sigma-Aldrich, MI, United States). The membrane was incubated with primary antibodies including Vimentin rabbit polyclonal antibody (1: 1,000; Cell Signaling Technology, MA, United States), E-cadherin rabbit polyclonal antibody (1:1,000; Cell Signaling Technology) and Glyceraldehyde 3-phosphate dehydrogenase (GADPH) rabbit polyclonal antibody (1:1,000; Cell Signaling Technology) at 4°C overnight. GADPH was used as the loading control. Finally, the secondary antibody, goat antirabbit IgG H&L (HRP) (1:1,000; Beyotime), was added to the membrane and subjected to incubation. Thereafter, the membrane was subjected to washing steps using PBS containing Tween-20 and the bands were visualized using a chemiluminescent detection system.

Tumor Xenograft Model
Four-to six-week-old male BALB/c nude mice were purchased from Guangxi Medical University Animal Center. All animal experiments were approved by the Animal Care & Welfare Committee of Guangxi Medical University, and animals were maintained at a specific pathogen-free environment and following the National Institutes of Health Guide for the Care and Use of Laboratory Animals. Ten mice were randomly assigned into two groups, namely SGC-7901-LINC01614 and SGC-7901-NC. For establishment of the tumor growth model, stably transfected SGC-7901 cells (1 × 10 7 ) were subcutaneously injected into the right lower groin of nude mice. The nude mice were sacrificed after 25 days. The tumor volume was calculated with the formula: volume length × width 2 /2.

Cell Cycle Assay
Propidium iodide (PI) stain was used to determine the cell cycle distribution. Cells were harvested at 48 h post-transfection. Thereafter, the cells were subjected to washing steps and were resuspended with PBS and centrifuged at 1,200 rpm for 5 min. After aspiration of the supernatant, cells were subjected to fixation using 1 ml pre-cooled (-20°C) 75% ethanol for 3 h at 4°C under dark conditions. Cells were then centrifuged at 12,000 rpm for 5 min and subjected to washing steps using PBS. Finally, the cells were incubated with 500 μL of the PI/ RNase Staining Buffer (BD Life Sciences, NJ, United States) for 15 min at room temperature under dark conditions. The Accuri C6 Plus flow cytometer (BD Life Sciences) was used for cell cycle observation.

Wound Healing Assay
At 48 h post-transfection, 70 μL cells were seeded into each well of ibidi Culture-Inserts. Remove the Culture-Insert to create the gap and record the image until they form an optically confluent monolayer. The cell debris was washed away using PBS. Images were recorded at 0, 36 and 48 h under ×100 magnification by Olympus I×73 inverted microscope (Olympus, Japan), respectively.

Transwell Assay
Cell invasion was assessed using 24-well transwell chambers (Corning Life Sciences, MA, United States), with the upper surface paved using 40 μL diluted Matrigel (BD Biosciences). GC cells (5 × 10 4 cells/well) were transferred to the upper chamber and 500 μL of DMEM containing 10% FBS was added to the lower chamber. Cells were then incubated 48 h. Thereafter, the solution in the upper chamber was removed and the cells were subjected to washing steps thrice and fixation using 4% paraformaldehyde for 30 min was performed. Cells were stained with 0.1% crystal violet and the cells that passed through the membrane were counted in five randomly chosen visual fields under ×100 magnification using an inverted microscope (Olympus, Japan).

Gene Set Enrichment Analysis (GSEA)
The data were categorized to indicate high-and low-LINC0614 (as control group) phenotypes based on the median value. Data were analyzed using "clusterProfiler" and KEGG enrichment was performed. The number of gene set permutations was set as 1,000. Nominal p value and normalized enrichment score were used to perform sorting of the enriched pathways.

Statistical Analyses
Multiple testing using the Benjamini-Hochberg procedure was performed to adjust the FDR. Log-rank test and univariate Cox regression were performed to examine the significance of the prognostic signature. The Chi-square test or the Fisher's exact test was used to compare the validation and current dataset. Clinical features with p < 0.1 were included in the multivariate Cox regression model and for conducting further stratified analysis. All bioinformatics statistical analyses were performed using the R software (version 3.5.3).
The statistical analyses of experiments were performed using GraphPad Prism 8.0. All data have been expressed as mean ± standard deviation or median. Differences were calculated using the Student's t-test or one-way analysis of variance followed by post hoc Dunnett's test. A p-value < 0.05 was considered statistically significant.

Differential Expression of lncRNAs in Gastric Cancer
Expression profiles of 23,151 genes obtained using RNA-Seq were downloaded from TCGA. Among these, 2,174 genes were annotated as lncRNAs and 858 lncRNAs with mean value > 1 were selected and identified as DELs. A total of 192 and 666 lncRNAs were upregulated and downregulated in GC, respectively. The heatmap of all 858 lncRNAs has been presented in Figure 1A.

Construction of the DELs-Based Signature
Among the 858 lncRNAs, 20 with p < 0.01 were considered potential candidates as prognostic DELs. Subsequently, using the "step" function, the following seven lncRNAs were selected to construct a prognostic signature: LINC01614, LINC01537, LINC01210, OVAAL, LINC01446, CYMP-AS1, and SCAT8. Moreover, multivariate Cox regression was used to determine the relative contribution of potential prognostic lncRNAs toward OS prediction; additionally, our study validated that all potential DELs were independent prognostic indicators except OVAAL ( Table 1). The risk score formula was imputed as per methods mentioned before. Kaplan-Meier curves showed that high expression levels of LINC01614, LINC01537, OVAAL, LINC01446, CYMP-AS1, and SCAT8 were associated with an increased risk of death; LINC01210 seemed to be the only lncRNA that followed an opposite trend and it seemed not to be a prognostic factor in GC based on log-rank test ( Figure 1B). After performing calculation of the risk score for each patient, median risk score was set as a cut-off point based on which the patients were categorized into high-and low-risk groups. Survival status analysis revealed that the high-risk group was subjected to an increased risk of death. Heatmap illustration showed that the protective lncRNA, LINC01210, exhibited a decreased expression in the low-risk group, while others exhibited an increasing trend ( Figure 1C). Patients in the high-risk group presented with a significantly decreasing median survival time of 560 days compared to those in the low-risk group (2,100 days; p < 0.001; Table 2 and Figure 1D). The ROC curve analysis showed an area under the curve (AUC) of 0.679, 0.756, and 0.786 for 1-, 3-, and 5-years OS prediction in GC, respectively ( Figure 1E).

Stratified Analysis and Analyses Combined With Risk Score and Clinical Parameters
To further confirm the combined prognostic value of clinical parameters and the prognostic signature, univariate Cox Frontiers in Genetics | www.frontiersin.org October 2021 | Volume 12 | Article 698947 Frontiers in Genetics | www.frontiersin.org October 2021 | Volume 12 | Article 698947 5 regression model was used. As shown in Table 2, risk score, age, tumor stage, and residual tumor of patients were associated with OS of GC patients. To consider all confounding factors, age, tumor stage, histological grade, and residual tumor with p < 0.1 were included in the multivariate Cox model and stratified analysis. After performing adjustments for the confounding factors, multivariate Cox model-based analysis showed that patients in the high-risk group were exposed to a significantly high risk of death (p < 0.0001, adjusted hazard ratio 2.338; 95% confidence interval 1.666-3.281 for OS). Stratified analysis conducted to explore the prognostic signature value under all clinical circumstances showed that the risk score could be efficiently used to predict the tumor stage, histological stage, and residual tumor in favorable and adverse strata, after subjection to adjustment for the confounding factors ( Figure 2A). The joint effect analysis based on the combination of the seven-lncRNA signature and clinical features including tumor stage, histological stage, and residual tumor, indicated that the signature showed appreciable performance in predicting GC patients' survival ( Figures 2B-D; Table 3).

Functional Analysis and Pathway Prediction
A total of 3,165 PCGs were identified as seven-lncRNAcorrelated genes (Figure 3). GO enrichment using "clusterProfiler" indicated that the seven lncRNAs were significantly associated with 1065 GO terms. The top ten associated GO terms of BP, CC, and MF, including extracellular matrix and cell adhering binding, are shown in Figure 2E. KEGG enrichment confirmed that a total of 321   pathways were enriched by the target genes. The top 20 enriched pathways included phosphatidylinositol 3-kinaseprotein kinase B (PI3K-Akt), focal adhesion, and cell cycle pathway ( Figure 2F).

LINC01614 Expression Is Upregulated in Gastric Cancer
Similar to the findings reported by a previous study, the seven-lncRNA signature might be considered a prognostic signature in GC; LINC01614 was the most significant lncRNA (p 5 × 10 -4 ).
To validate the findings, LINC01614 was selected in the following study. Through analyzing the dataset GSE95667 in GEO database, it is found that the expression of LINC01614 in gastric cancer tissues is also upregulated with the log2FC of 1.953 (p < 0.05). Consistent with TCGA analysis results, analysis based on qRT-PCR revealed the existence of an increased expression of LINC01614 in the GC tissues compared to the adjacent normal tissues ( Figure 4A). Consistently, LINC01614 expression also increased in SGC-7901, HGC-27, AGS, and MGC-803 compared to GES-1 ( Figure 4B, p < 0.05).

LINC01614 Promotes Proliferation in Vitro and in Vivo
AGS and MGC-803 with relatively high expression levels of LINC01614 were used in subsequent experiments. After  Vimentin expression increased; LINC01614 overexpression promoted the epithelial-mesenchymal transition (EMT) process ( Figure 4G). By performing subcutaneous implantation of the transfected SGC-7901, we found that tumor size in the SGC-7901-LINC01614 group was remarkably bigger than that observed in the SGC-7901-NC group; the observation was consistent with the outcomes of the in vitro experiment ( Figures 4H,I).  Figure 5A). Results indicated that the downregulated expression of LINC01614 might lead to the arrest of cells at the G2/M phase. Wound healing assay showed that the migration ability was weakened in MGC-803 and AGS cells transfected with si-LINC01614-1 or si-LINC01614-2 compared to the NC group ( Figure 5B). The number of invasive cells decreased in groups showing inhibited LINC01614 expression compared to the NC group in transwell assay ( Figure 5C). Thus, downregulated expression of LINC01614 inhibited the migration and invasion of GC cells.

GSEA Analysis of LINC01614
GSEA was used to further validate the involvement of LINC01614 in GC pathogenesis. Enriched pathways, including PI3K-Akt signaling pathway, mitogen-activated protein kinase signaling pathway, focal adhesion, Wnt signaling pathway, Toll-like receptor signaling pathway, tumor growth factor-beta signaling pathway, and ECM-receptor interaction revealed that LINC01614 might exert potential regulatory effect in GC through multiple mechanisms ( Figure 5D).

DISCUSSION
TCGA, a milestone-like program, established a comprehensive dataset of mutations that occur in 33 cancer types (Weinstein et al., 2013). Based on the data on lncRNA expression profiles extracted from TCGA, the present study showed that aberrant expression levels of LINC01614, LINC01537, LINC01210, OVAAL, LINC01446, CYMP-AS1, and SCAT8 may associated with OS in GC patients and these seven lncRNAs could be used to construct a prognostic signature in GC. The ROC analysis indicated that the signature showed appreciable performance in survival prediction within 5 years. Kaplan-Meier survival curves reiterated the same finding. And Our study showed an improvement in the survival prediction compared to the signature reported before. (Miao et al., 2017;. Aberrant expression of lncRNAs and their role in GC pathogenesis have attracted the attention of multiple researchers. Xu et al. have validated that PVT1 with upregulated expression in GC through activation of FOXM1 transcription acts as a oncogene. (Lian et al., 2020). Additionally, a study confirmed AK023391 as a novel oncogenic lncRNA in GC that exhibited functions by mediating the PI3K-Akt signaling pathway (Huang et al., 2017). All such studies have highlighted the convincing regulatory role of lncRNAs in tumorigenesis and progression, and also lays the foundation for the conduction of future studies dedicated to unravelling the functions of lncRNA in GC. Inspiringly, some of the GC survival-related lncRNAs show commonality with the previous studies conducted on cancer.
OVAAL, has been previously identified to play a significant role in cell survival, proliferation, and evasion from cellular senescence and contributes to cancer cell survival (Sang et al., 2018). LINC01446 functions as an oncogene and promotes GC and glioblastoma progression and indicates a poor prognosis, an effect similar to that predicted here. (Lian et al., 2020;Zhang et al., 2018).
LINC01614 is located on human chromosome 2q35, so far, many studies have shown that it promotes the progression of different cancers, such as lung cancer , breast cancer (Vishnubalaji et al., 2019), glioma (Wang Y et al., 2020) and osteosarcoma (Cai et al., 2021). Similarly, Chen et al.' work addresses part of the same issues and reports similar results as those presented here. And the two independent pieces of work arrive at essentially the same conclusions that LINC01614 acts as an oncogene in GC and gives added weight to the finding of both studies. Chen et al. reported that LINC01614 has the potential to segregate a subgroup of L1 patients which was defined by authors with worse survival, however, the experiment to illustrate the relative expression of LINC01614 in GC samples and in vivo study was not performed. (Chen et al., 2021). And in our study, in vitro and in vivo experiments revealed the oncogenetic character of LINC01614 which is the most significant candidate of the signature in GC. In the study, LINC01614 expression was observed to be upregulated in GC tissue obtained compared to the adjacent normal tissue in our cohort; upregulated expression of LINC01614 in GC tissues and cell lines were consistent with the bioinformatics analysis results as well. These findings suggest that LINC01614 may play a role in tumorigenesis. We then examined the function of LINC01614 in GC cells in vitro, the results revealed that LINC01614 stimulated the development of gastric cancer, it promotes cell proliferation and colony formation. The effects of knockdown of LINC01614 using siRNAs transfection were revealed as followed: the attenuation of cell proliferation, inhibition of cell migration, invasion and EMT, leading to G2/M cell cycle arrest.
Nevertheless, the limitations of this study should be acknowledged. Despite the significant upregulation of LINC01614 characterized here in GC along with previous reports in other cancers, the molecular mechanisms still are not clearly seen, and more explorations are also necessary to detect upstream or downstream effectors of LINC01614 in GC progression. For another, the prognostic value of the lncRNAs signature and the expression of the lncRNA after treatment should be validated in other cohorts with large sample in future. In conclusion, this study identified a seven-lncRNA signature comprising LINC01614 and six lncRNAs to predict OS in GC patients. Moreover, LINC01614 expression was reported to be markedly upregulated in GC, and it may promote gastric carcinogenesis.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.