Background. Long noncoding RNA-based prognostic biomarkers have demonstrated great potential in the diagnosis and prognosis of cancer patients. However, systematic assessment of a multiple lncRNA-composed prognostic risk model is lacking in stomach adenocarcinoma (STAD). This study is aimed at constructing a lncRNA-based prognostic risk model for STAD patients. Methods. RNA sequencing data and clinical information of STAD patients were retrieved from The Cancer Genome Atlas (TCGA) database. Differentially expressed lncRNAs (DElncRNAs) were identified using the R software. Univariate and multivariate Cox regression analyses were performed to construct a prognostic risk model. The survival analysis, C-index, and receiver operating characteristic (ROC) curve were employed to assess the sensitivity and specificity of the model. The results were verified using the GEPIA online tool and our clinical samples. Pearson correlation coefficient analysis, Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment were performed to indicate the potential biological functions of the selected lncRNA. Results. A total of 1917 DElncRNAs were identified from 343 cases of STAD tissues and 30 cases of noncancerous tissues. According to univariate and multivariable Cox regression analyses, four DElncRNAs (AC129507.1, LINC02407, AL022316.1, and AP000695.2) were selected to establish a prognostic risk model. There was a significant difference in the overall survival between high-risk patients and low-risk patients based on this risk model. The C-index of the model was 0.652. The area under the curve (AUC) for the ROC curve was 0.769. GEPIA results confirmed the expression and prognostic significance of AP000695.2 in STAD. Our clinical data confirmed that upregulated expression of AP000695.2 was correlated with the T stage, distant metastasis, and TNM stage in STAD. GO and KEGG analyses demonstrated that AP000695.2 was closely related to the tumorigenesis process. Conclusions. In this study, we constructed a lncRNA-based prognostic risk model for STAD patients. Our study will provide novel insight into the diagnosis and prognosis of STAD patients.
Guangzhou General Science and Technology Project of Health and Family Planning20201A011007Fundamental Research Funds for the Central UniversitiesX2yx-D2175050Guangzhou Planned Project of Science and Technology201707010262201904010427Natural Science Foundation of Guangdong Province2016A030310109National Natural Science Foundation of China81502563816729001. Introduction
Stomach adenocarcinoma (STAD), as the most common type of gastric cancer (GC), is characterized by rapid growth and strong invasiveness [1]. It is among one of the leading causes of cancer-related death worldwide [2, 3]. Approximately 70% of stomach adenocarcinoma was diagnosed in developing countries [3]. To date, due to the late stage of diagnosis and lack of effective treatment strategies, the prognosis of stomach adenocarcinoma is unsatisfactory [4, 5]. Moreover, the development of targeted therapy and immune checkpoint inhibitors only benefits a very small population of gastric cancer patients [6–8]. The main reason for these phenomena is the lack of effective diagnosis and prognostic evaluation measures for STAD at an early stage [9, 10]. Therefore, identifying novel biomarkers which could diagnose and predict the survival of STAD patients is critical.
Long noncoding RNA (lncRNA) belongs to noncoding RNA molecules, which are longer than 200 nucleotides (nt) at length [11]. lncRNAs were initially considered to have no physiological function but a byproduct of RNA polymerase II transcription [12]. However, due to the development of next-generation sequencing technologies, lncRNAs emerge as crucial regulators in tumorigenesis [13–15]. Recent studies revealed that lncRNAs participate in various biological processes including transcriptional regulation, RNA editing, and posttranscriptional regulation of many genes [14, 16, 17]. There were close correlations between lncRNAs and various cellular events, including cell proliferation, migration, invasion, cell cycle, and apoptosis [18–20]. With regard to gastric cancer, lncRNAs have been detected and could function as oncogenes or tumor suppressors [21, 22]. Dysregulated expression of lncRNAs demonstrated important roles in predicting the relapse, metastasis, chemoresistance, and overall survival of gastric cancer [20, 23]. Upregulated expression of lncRNA AC093818.1 demonstrated high sensitivity in predicting metastatic gastric cancer and could accelerate gastric cancer metastasis [20]. Cisplatin-resistant gastric cancer patients demonstrated higher lncRNA HOXD-AS1 expression [23]. lncRNA LOC100130476 was related to tumor suppression and aberrant methylation in gastric cardia adenocarcinoma [24]. lncRNA N-BLR could promote the migration and invasion of gastric adenocarcinoma [25]. In addition, lncRNAs RP11-169F17.1 and RP11-669N7.2 were regarded as prognostic biomarkers of stomach adenocarcinoma [26]. lncRNA RP1-163G9.1 was associated with poor prognosis of gastric adenocarcinoma patients [12]. All these studies indicate that lncRNAs could serve as prognostic biomarkers for gastric cancer patients. However, the specificity and sensitivity of single lncRNA as a biomarker are limited. An integrated lncRNA prognostic risk model would play more vital roles in the diagnosis and prognosis of STAD patients [12, 24–26].
In this study, we established a prognostic risk model by analyzing RNA sequencing data from TCGA dataset. The results were verified using the GEPIA online tool and our clinical data. This study will provide novel insight into the prognostic values of the lncRNA-based model, which could be used in tumor diagnosis and survival prognosis for STAD patients.
2. Material and Methods2.1. Ethical Statement
A total of 78 STAD patients who underwent surgery between January 2010 and December 2017 were included in this study. None of them had received any antitumor therapies prior to the collection of the tissue samples. There were 25 females and 53 males, with a median age of 52 years (range, 21-75 years). All tumor tissues and adjacent normal tissues were obtained with informed consent, and the study was approved by the Research Ethics Committee of Guangzhou First People’s Hospital, South China University of Technology.
2.2. Data Retrieval and Identification of Differently Expressed lncRNAs (DElncRNAs)
The RNA sequence data and corresponding clinical information from 343 STAD patients were downloaded from TCGA dataset (up to December 2019). There were 343 cases of tumor tissues and 30 cases of adjacent nontumoral tissues included in this study. The data was retrieved using the Perl software. The DElncRNAs were compared between tumor tissues and adjacent nontumoral tissues using the “edgeR” package based on the R language. DElncRNAs were identified using the following criteria: ∣logfoldchange∣>1 and adjusted P value < 0.05. The heatmap was drawn using the “pheatmap” package in the R software.
2.3. Survival Analysis Using Kaplan-Meier Survival Analysis
The correlations between DElncRNAs and overall survival of STAD patients were analyzed using Kaplan-Meier survival analysis with the log-rank test in the R software, where P<0.05 was considered statistically significant.
2.4. Cox Regression Analysis
Univariate and multivariate Cox regression analyses were performed to establish a prognostic risk model which could independently predict the survival of STAD patients. After being subjected to univariate Cox regression analysis, the DElncRNAs with P value < 0.001 were further subjected to a multivariate Cox proportional hazards model in the R software.
2.5. Assessment of the Prognostic Risk Model
Kaplan-Meier survival analysis was used to assess the survival curve in the high-risk and low-risk groups using the “survival” package in the R software. The C-index value of the prognostic risk model was calculated using the “survcomp” package. Receiver operating characteristic (ROC) curve analysis was performed using the “survivalROC” package in the R software.
2.6. Validation of the Expression and Prognostic Values of the Independent Prognostic lncRNAs
The expression and prognostic significance of these independent prognostic lncRNAs were further verified using the online tool Gene Expression Profiling Interactive Analysis (GEPIA, http://http://gepia.cancer-pku.cn/index.html).
2.7. Real-Time Quantitative PCR (qRT-PCR)
Total RNAs were collected from 78 paired STAD tumor tissues and adjacent normal tissues and extracted by the TRIzol reagent (Invitrogen, USA). The synthesis of cDNA was conducted by using HiScript II 1st Strand cDNA Synthesis Kit (Vazyme Biotech, Nanjing, China) from isolated RNA [27, 28]. The cDNA was subjected to the real-time quantitative PCR using the SYBR Green Master Mix. GAPDH and 18S were used as internal controls. The sequences of the sense and antisense primers were as follows: 5′-GGACACTCTGAAGGAACTC-3′ (F) and 5′-GATGACCATTAGCCAACAAG-3′ (R) for AP000695.2; 5′-CTCCTCCTGTTCGACAGTCAGC-3′ (F) and 5′-CCCAATACGACCAAATCCGTT-3′ (R) for GAPDH; and 5′-CGGCGACGACCCATTCGAAC-3′ (F) and 5′-GAATCGAACCCTGATTCCCCGTC-3′ (R) for 18S. The comparative threshold cycle (2−ΔΔCT) method was used to calculate the relative mRNA values. The relative expression levels of AP000695.2 were normalized to the value of GAPDH.
2.8. Correlation between lncRNA and Coexpressed Protein Coding Genes
The correlations between lncRNAs and their coexpressed core genes were assessed using the “limma” package in the R software. The results were considered statistically significant when the Pearson correlation coefficient was greater than 0.4 and P value < 0.05.
2.9. Gene Oncology (GO) Annotation and KEGG Pathway Enrichment Analysis
GO enrichment analysis and KEGG pathway enrichment were performed using the “clusterProfiler” package in the R software. A P value < 0.05 was considered statistically significant.
2.10. Statistical Analysis
Data are presented as mean±SD. The statistical analyses were assessed using the GraphPad Prism software (version 7.0). Two-tailed Student’s t-test was used to assess the statistical differences between the two groups. The chi-square test or Fisher’s exact test was used to compare the differences between AP000695.2 expression and clinical parameters. A two-tailed P value < 0.05 was considered statistically significant.
3. Results3.1. Identification of Differentially Expressed lncRNAs (DElncRNAs) in STAD Patients
The overview of the screening strategy is shown in Figure 1(a). The identification of DElncRNAs was compared between 343 cases of tumor tissues and 30 cases of adjacent normal tissues based on TCGA database. Volcano plots demonstrated the distribution of DElncRNAs in STAD patients (Figure 1(b)). A total of 1917 DElncRNAs were identified, including 438 downregulated and 1497 upregulated lncRNAs. Figure 1(c) shows the heatmap of the DElncRNAs compared between STAD tumor tissues and adjacent nontumoral tissues.
Identification of differentially expressed lncRNAs (DElncRNAs) in stomach adenocarcinoma patients. (a) The main flowchart of the current study. The data was retrieved from TCGA database. Differentially expressed lncRNAs were compared between 343 stomach adenocarcinoma tissues and 30 adjacent nontumoral tissues. (b) Volcano map of DElncRNAs. Red plots demonstrate upregulated lncRNAs, while green plots demonstrate downregulated lncRNAs. (c) Heatmap of the DElncRNAs between tumor tissues and adjacent nontumoral tissues. Red, black, and green colors represent the relatively high, medium, and low expression of lncRNAs, respectively. STAD: stomach adenocarcinoma; TCGA: The Cancer Genome Atlas; ROC: receiver operating characteristic curve; GO: Gene Ontology; KEGG: Kyoto Encyclopedia of Genes and Genomes; GEPIA: Gene Expression Profiling Interactive Analysis.
3.2. Construction of a Prognostic Risk Model in STAD Patients
To study the potential values of DElncRNAs in predicting the overall survival of STAD patients, we constructed a prognostic risk model based on the univariate and multivariate COX analysis. The univariate analysis was performed to identify possible diagnostic and prognostic biomarkers for STAD patients. Table 1 shows the identified DElncRNAs with P value < 0.001. With a standard of P value < 0.001, 7 DElncRNAs (AL022316.1, AC129507.1, AP000695.1, FLJ16779, LINC01537, LINC02407, and AP000695.2) were subjected to further multivariate Cox regression analysis. We constructed a prognostic risk model containing 4 DElncRNAs (AC129507.1, LINC02407, AL022316.1, and AP000695.2). The forest plots demonstrated that AC129507.1 (P=0.003), AL022316.1 (P<0.001), and AP000695.2 (P=0.012) could be used as independent biomarkers for STAD patients (Figure 2(a)).
Univariate analysis results of the statistically significant lncRNAs (P<0.001).
Gene
HR
Z
P value
AL022316.1
0.79131398
-3.91180144
9.16E−05
AC129507.1
1.31492429
3.70563892
0.000210859
AP000695.1
1.28943052
3.51995878
0.000431614
FLJ16779
1.15788429
3.50127829
0.000463032
LINC01537
1.24779825
3.46898743
0.000522424
LINC02407
1.27355052
3.30854576
0.000937819
AP000695.2
1.23639943
3.30519073
0.000949118
Construction and performance of a prognostic risk model. (a) Multivariate Cox regression analysis and forest plots demonstrate that key DElncRNAs could be viewed as independent prognostic biomarkers. (b) Kaplan-Meier curves demonstrate the difference between the high-risk group and low-risk group in OS based on the prognostic risk model. (c) ROC analysis of the prognostic risk model. (d) Scatter diagram shows the distributions of risk scores of STAD patients. (e) Scatter diagram shows the survival status of the patients based on this prognostic risk model. (f) Heatmap of the 4 prognostic lncRNA expression profiles.
We next assessed the sensitivity and specificity of this risk model. Based on the prognostic risk model, patients could be divided into the high-risk group and low-risk group. As shown in Figure 2(b), the patients classified into the high-risk group demonstrated a poorer survival rate than patients in the low-risk group (P=4.418e−09). The C-index of this risk model is 0.652 (95% CI, 0.606-0.698, Table 2). The ROC was calculated, and the value of the AUC is 0.769 (Figure 2(c)). All these data indicated that this prognostic risk model demonstrated good sensitivity and specificity. Figure 2(d) shows the risk score distribution of the STAD patients. Figure 2(e) shows the survival status of the STAD patients based on the risk model. The heatmap shows the lncRNA expression profiles of the prognostic risk model based on the risk score (Figure 2(f)).
C-index analysis of the prognostic risk model.
C-index
SE
Lower
Upper
P value
0.6518902
0.02357405
0.6056859
0.6980945
1.1705E−10
3.3. Expression and Prognostic Significance of the DElncRNAs in the Prognostic Risk Model Based on TCGA Database
The different expression levels of the independent prognostic lncRNAs are summarized in Table 3. Survival analysis for the independent prognostic lncRNAs was calculated using Kaplan-Meier survival analysis. As shown in Figure 3, upregulated expression of AC129507.1 (P=5.925e−09, Figure 3(a)), LINC02407 (P=2.101e−03, Figure 3(b)), and AP000695.2 (P=4.447e−02, Figure 3(c)) was positively correlated with poor overall survival, while high expression of AL022316.1 (P=8.944e−04, Figure 3(d)) was positively correlated with better OS.
Expression of DElncRNAs in the prognostic risk model based on TCGA database.
DElncRNAs
logFC
P value
FDR
AC129507.1
-1.66040229
1.40E−18
4.35E-17
LINC02407
2.69475885
7.73E−15
1.40E−13
AL022316.1
2.04469210
2.04E−08
1.34E−07
AP000695.2
1.59480729
1.45E−07
8.08E−07
The correlations between prognostic lncRNAs and overall survival of STAD patients. Kaplan-Meier curves with the log-rank test for AC129507.1 (a), LINC02407 (b), AP000695.2 (c), and AL022316.1 (d).
3.4. Expression and Prognostic Values of the DElncRNAs in the Prognostic Risk Model Using GEPIA
The online tool GEPIA was conducted to verify the expression and prognostic values of the prognostic DElncRNAs. A total of 408 tumor samples and 211 normal samples were included in GEPIA. As shown in Figure 4(a), upregulated expression of AP000695.2 was detected in the STAD tumor tissues. These data were in accordance with the data retrieved from TCGA database (Table 3). Overall survival analyses were also in accordance with the results of TCGA database. As shown in Figure 4(b), the high AP000695.2 expression (P=0.024, Figure 4(b)) was positively correlated with poorer OS.
The expression and prognostic value of AP000695.2 were analyzed by GEPIA. There were 408 tumor samples and 211 normal samples in the GEPIA database. (a) The expression of AP000695.2 in STAD tumor tissues and normal tissues. (b) Kaplan-Meier survival curve for AP000695.2.
3.5. Validation of the Expression and Prognostic Value of AP000695.2 in STAD Patients
The expression of AP000695.2 in 78 paired STAD tumor tissues and adjacent normal tissues was compared using real-time quantitative PCR (qRT-PCR). As shown in Figure 5(a), GAPDH was used as an endogenous control and the results demonstrated that the relative mRNA level of AP000695.2 was much higher in tumor tissues than in normal tissues (P<0.001). Application of 18S as another internal control also confirmed the upregulated expression of AP000695.2 in STAD tumor tissues (P<0.001, Supplementary Figure 1). The mean relative expression level of AP000695.2 was 2.141 (range, 0.19-5.04). The patients were divided into two groups based on the mean mRNA value of AP000695.2: low-expression group (N=45) and high-expression group (N=33). Moreover, the correlations between AP000695.2 expression and clinicopathological parameters were assessed (Figures 5(b)–5(f)). Upregulated expression of AP000695.2 was positively correlated with the T stage (P=0.025, Figure 5(c)), distant metastasis (P=0.001, Figure 5(e)), and TNM stage (P=0.007, Figure 5(f)). However, there was no statistical significance between AP000695.2 expression and tumor size (P=0.099, Figure 5(b)) and lymph node metastasis (P=0.285, Figure 5(d)). Taken together, TCGA database, GEPIA online database, and our clinical samples indicated that lncRNA AP000695.2 might work as an oncogene in STAD.
Validation of AP000695.2 expression and prognostic value in STAD patients. (a) The mRNA expression level of AP000695.2 in 78 paired STAD tumor tissues and adjacent normal tissues using real-time quantitative PCR (qRT-PCR). GAPDH was used as an endogenous control. The correlations between AP000695.2 expression and tumor size (b), T stage (c), lymph node metastasis (d), metastasis (e), and TNM stage (f) were analyzed.
3.6. Coexpressed Genes of AP000695.2
We further studied the coexpressed protein coding genes of AP000695.2 to clarify its potential biological functions. According to the cut-off criteria (Pearson correlation coefficient>0.4 or <-0.4 and P value < 0.001), AP000695.2 coexpressed with 59 mRNAs. Figure 6 shows the top 5 correlating mRNAs of AP000695.2. AP000695.2 expression was positively correlated with CLDN14 (Figure 6(a)), PLAUR (Figure 6(b)), TMP1 (Figure 6(c)), TNFRSF12A (Figure 6(d)), and SPOCD1 (Figure 6(e)).
The correlations between AP000695.2 and its coexpressed genes. The correlation plots of top 5 correlated mRNAs of AP000695. Pearson correlation coefficient>0.4 and P<0.01.
3.7. GO and KEGG Pathway Enrichment Analysis of the Independent Prognostic DElncRNAs
To further clarify the potential biological functions of AP000695.2 in STAD, we analyzed GO and KEGG pathway enrichment analysis. As shown in Figure 7(a), the GO enrichment results demonstrated obvious correlation with “extracellular matrix structural constituent,” “endopeptidase activity,” “metalloendopeptidase activity,” “metallopeptidase activity,” and “cell adhesion molecule binding” for AP000695.2. KEGG pathway enrichment analysis demonstrated that these genes were correlated with “complement and coagulation cascades,” “proteoglycans in cancer,” “ECM-receptor interaction,” “protein digestion and adsorption,” and “AGE-RAGE signaling pathway in diabetic complications” pathways (Figure 7(b)). All these data suggested that AP000695.2 was involved in multiple cellular processes, including the processes of tumor development.
Functional enrichment analysis for AP000695.2 coexpressed protein coding genes: (a) bubble chart of GO enrichment and (b) bubble chart of KEGG pathway enrichment.
4. Discussion
Most of the stomach adenocarcinoma patients were diagnosed at an advanced stage. Thus, the treatment strategies for most of these patients were limited [29]. The advancement in targeted therapies and immunotherapies failed to demonstrate promising results in STAD patients [30, 31]. And they are still lacking molecular drugs that can specifically target STAD. With the development of the next-generation sequence, studies about noncoding RNAs have greatly accelerated [32, 33]. The development of RNA sequencing technology has facilitated the identification of novel lncRNAs [34]. A bunch of lncRNAs have been reported to play important roles in gastric cancer cell proliferation, invasion, and metastasis [35, 36]. Dysregulation of lncRNAs has been reported to be correlated with tumorigenesis and progression of stomach adenocarcinoma. Moreover, the expression and biological functions of lncRNAs in STAD have been widely explored [37].
Nowadays, lncRNA-based biomarkers have been widely used in the diagnosis and prognosis of gastric cancer patients [38, 39]. lncRNA SNHG15 demonstrated a prospective biomarker for the diagnosis and therapeutics for gastric cancer patients [39]. lncRNA MAGI2-AS3 was identified as epithelial-mesenchymal transition- (EMT-) related lncRNA, which could predict the progression of gastric cancer patients [38]. Moreover, some lncRNAs could be detected in plasma and serve as noninvasive biomarkers, which makes lncRNA an ideal biomarker for cancer patients [40]. Zheng et al. demonstrated that serum lncRNAs FAM49B-AS, GUSBP11, and CTDHUT could be viewed as biomarkers for gastric cancer patients [40]. In this study, we identified 1917 differentially expressed lncRNAs between STAD tumor tissues and adjacent nontumoral tissues. Then, we constructed a four lncRNA- (AC129507.1, LINC02407, AL022316.1, and AP000695.2) based prognostic risk model based on these DElncRNAs. This prognostic risk model demonstrated good sensitivity and specificity in predicting the risk and survival of STAD patients. According to this model, our clinical data also confirmed that upregulated expression of independent prognostic lncRNA AP000695.2 was positively correlated with the T stage, M stage, and TNM stage. All these data indicated that AP000695.2 might work as an oncogenic gene in STAD patients. However, there were no articles concerning the expression and the biological functions of AP000695.2 in cancer patients.
lncRNAs work in multiple ways to regulate the expression and functions of target genes, including signal, decoy, guide, and scaffold [41, 42]. According to the results of GO enrichment and KEGG pathway enrichment, AP000695.2 coexpressed genes were enriched in biological functions and pathways including extracellular matrix structural constituent, ECM-receptor interaction, and proteoglycans in cancer. All these data indicated that lncRNA AP000695.2 might be involved in the cellular process of tumorigenesis. Thus, further studies are needed to elucidate the biological functions and molecular mechanisms of AP000695.2 in STAD.
5. Conclusion
In conclusion, DElncRNAs significantly associated with oncogenesis and prognosis of STAD were identified from TCGA database. The independent prognostic lncRNAs demonstrated good sensitivity and specificity in predicting the survival of STAD patients. Moreover, our preliminary studies revealed that upregulated expression of AP000695.2 was correlated with poor overall survival and the biological functions of AP000695.2 were predicted. However, there are no researchers once clarifying the functions and mechanisms of these DElncRNAs. Our study provides a novel lncRNA-based prognostic risk model and a new direction for future studies, where we will further explore the biological functions and mechanisms of these lncRNAs in STAD. Meanwhile, we will detect the expression of these DElncRNAs using our own clinical samples to further confirm the sensitivity and specificity of these DElncRNAs in predicting the survival of STAD patients.
AbbreviationsSTAD:
Stomach adenocarcinoma
TCGA:
The Cancer Genome Atlas
DElncRNAs:
Differentially expressed lncRNAs
ROC:
Receiver operating characteristic curve
GO:
Gene Ontology
KEGG:
Kyoto Encyclopedia of Genes and Genomes
lncRNA:
Long noncoding RNA
GC:
Gastric cancer
CLDN14:
Claudin 14
PLAUR:
Plasminogen activator urokinase receptor
TMP1:
Thymidylate synthase gene
TNFRSF12A:
Tumor necrosis factor receptor superfamily member 12A
SPOCD1:
Spen paralogue and orthologue C-terminal domain containing 1.
Data Availability
All data in our study are available from the corresponding authors upon reasonable request.
Conflicts of Interest
All the authors declare that there are no conflicts of interest.
Acknowledgments
This research was funded by the National Natural Science Foundation of China (Nos. 81672900 and 81502563), the Natural Science Foundation of Guangdong Province of China (No. 2016A030310109), the Guangzhou Planned Project of Science and Technology (Nos. 201904010427 and 201707010262), the Fundamental Research Funds for the Central Universities (No. X2yx-D2175050), and the Guangzhou General Science and Technology Project of Health and Family Planning (grant number 20201A011007).
Supplementary Materials
Supplementary Figure 1: the mRNA expression level of AP000695.2 in 78 paired STAD tumor tissues and adjacent normal tissues using real-time quantitative PCR (qRT-PCR). 18S was used as an endogenous control.
GuZ.WangL.YaoX.LongQ.LeeK.LiJ.YueD.YangS.LiuY.LiN.LiY.ClC-3/SGK1 regulatory axis enhances the olaparib-induced antitumor effect in human stomach adenocarcinoma2020111089810.1038/s41419-020-03107-333093458HaradaK.LopezA.ShanbhagN.BadgwellB.BabaH.AjaniJ.Recent advances in the management of gastric adenocarcinoma patients2018710.12688/f1000research.15133.12-s2.0-85052972400FerlayJ.SoerjomataramI.DikshitR.EserS.MathersC.RebeloM.ParkinD. M.FormanD.BrayF.Cancer incidence and mortality worldwide: sources, methods and major patterns in GLOBOCAN 201220151365E359E38610.1002/ijc.292102-s2.0-8501733085725220842BanksM.GrahamD.JansenM.GotodaT.CodaS.di PietroM.UedoN.BhandariP.PritchardD. M.KuipersE. J.Rodriguez-JustoM.NovelliM. R.RagunathK.ShepherdN.Dinis-RibeiroM.British Society of Gastroenterology guidelines on the diagnosis and management of patients at risk of gastric adenocarcinoma20196891545157510.1136/gutjnl-2018-3181262-s2.0-8506866422831278206HaradaK.Mizrak KayaD.ShimodairaY.AjaniJ. A.Global chemotherapy development for gastric cancer201720S19210110.1007/s10120-016-0655-82-s2.0-84990929612WangB.-C.ZhangZ.-J.FuC.WangC.Efficacy and safety of anti-PD-1/PD-L1 agents vs chemotherapy in patients with gastric or gastroesophageal junction cancer: a systematic review and meta-analysis20199847, article e1805410.1097/MD.000000000001805431764833BrarG.ShahM. A.The role of pembrolizumab in the treatment of PD-L1 expressing gastric and gastroesophageal junction adenocarcinoma201912, article 175628481986976710.1177/17562848198697672-s2.0-85072065982IlsonD. H.Advances in the treatment of gastric cancer201935655155410.1097/MOG.00000000000005772-s2.0-8507150846431436556Uchôa GuimarãesC.Ferreira MartinsN.Cristina da Silva OliveiraK.AlmeidaC.PinheiroT.GigekC.Roberto de Araújo CavalleroS.AssumpçãoP. P.Cardoso SmithM. A.BurbanoR. R.CalcagnoD. Q.Liquid biopsy provides new insights into gastric cancer2018919151441515610.18632/oncotarget.245402-s2.0-8504359268329599934YoshidaT.YamaguchiT.MaekawaS.TakanoS.KunoT.TanakaK.IwamotoF.TsukuiY.KobayashiS.AsakawaY.ShindoH.FukasawaM.NakayamaY.InoueT.UetakeT.OhtakaM.SatoT.MochizukiK.EnomotoN.Identification of early genetic changes in well-differentiated intramucosal gastric carcinoma by target deep sequencing201922474275010.1007/s10120-019-00926-y2-s2.0-85061495507SchmittA. M.ChangH. Y.Long noncoding RNAs in cancer pathways201629445246310.1016/j.ccell.2016.03.0102-s2.0-8496273933627070700SunB.DangY.ZhangF.LiK.OuyangX.WangK.HuangQ.Long non-coding RNA RP1-163G9.1 is downregulated in gastric adenocarcinoma and is associated with a poor prognosis20194163575358510.3892/or.2019.71272-s2.0-8506445782431002370HuX.SoodA. K.DangC. V.ZhangL.The role of long noncoding RNAs in cancer: the dark matter matters20184881510.1016/j.gde.2017.10.0042-s2.0-8503150061529054012TangY.CheungB. B.AtmadibrataB.MarshallG. M.DingerM. E.LiuP. Y.LiuT.The regulatory role of long noncoding RNAs in cancer2017391121910.1016/j.canlet.2017.01.0102-s2.0-8501090368028111137HuarteM.The emerging role of lncRNAs in cancer201521111253126110.1038/nm.39812-s2.0-8494619688026540387SunQ.HaoQ.PrasanthK. V.Nuclear long noncoding RNAs: key regulators of gene expression201834214215710.1016/j.tig.2017.11.0052-s2.0-8503944159729249332ShihJ. W.KungH. J.Long non-coding RNA and tumor hypoxia: new players ushered toward an old arena20172415310.1186/s12929-017-0358-42-s2.0-8502691525328789687WangY.LuJ. H.WuQ. N.JinY.WangD. S.ChenY. X.LiuJ.LuoX. J.MengQ.PuH. Y.WangY. N.HuP. S.LiuZ. X.ZengZ. L.ZhaoQ.DengR.ZhuX. F.JuH. Q.XuR. H.LncRNA LINRIS stabilizes IGF2BP2 and promotes the aerobic glycolysis in colorectal cancer201918117410.1186/s12943-019-1105-031791342HanY.WuN.JiangM.ChuY.WangZ.LiuH.CaoJ.LiuH.XuB.XieX.Long non-coding RNA MYOSLID functions as a competing endogenous RNA to regulate MCL-1 expression by sponging miR-29c-3p in gastric cancer2019526, article e1267810.1111/cpr.126782-s2.0-8507190689431497917BaM. C.BaZ.LongH.CuiS. Z.GongY. F.YanZ. F.LinK. P.WuY. B.TuY. N.LncRNA AC093818.1 accelerates gastric cancer metastasis by epigenetically promoting PDK1 expression20201116410.1038/s41419-020-2245-231988283ZhengJ.ZhangH.MaR.LiuH.GaoP.Long non-coding RNA KRT19P3 suppresses proliferation and metastasis through COPS7A-mediated NF-κB pathway in gastric cancer201938457073708810.1038/s41388-019-0934-z2-s2.0-8507081939631409899FengL.LiH.LiF.BeiS.ZhangX.LncRNA KCNQ1OT1 regulates microRNA-9-LMX1A expression and inhibits gastric cancer cell progression202012170771710.18632/aging.10265131915311YeY.YangS.HanY.SunJ.XvL.WuL.MingL.HOXD-AS1 confers cisplatin resistance in gastric cancer through epigenetically silencing PDCD4 via recruiting EZH2201999, article 19006810.1098/rsob.19006831551012GuoW.DongZ.ShiY.LiuS.LiangJ.GuoY.GuoX.ShenS.WangG.Methylation-mediated downregulation of long noncoding RNA LOC100130476 in gastric cardia adenocarcinoma201633549750810.1007/s10585-016-9794-x2-s2.0-8496921147827189370YounY. H.ByunH. J.YoonJ. H.ParkC. H.LeeS. K.Long noncoding RNA N-BLR upregulates the migration and invasion of gastric adenocarcinoma201913442142910.5009/gnl184082-s2.0-8506929460930970439YangJ.SongH.Identification of long noncoding RNA RP11-169F17.1 and RP11-669N7.2 as novel prognostic biomarkers of stomach adenocarcinoma based on integrated bioinformatics analysis201911111307132110.2217/epi-2019-01152-s2.0-8507223402131368349LuL.duH.HuangH.WangC.WangP.ZhaZ.WuY.LiuX.WengC.FangX.LiB.MaoH.WangL.GuanM.LiuG.CCR9 promotes migration and invasion of lung adenocarcinoma cancer stem cells202017791292010.7150/ijms.4086432308544WangW.LvL.PanK.ZhangY.ZhaoJ. J.ChenJ. G.ChenY. B.LiY. Q.WangQ. J.HeJ.ChenS. P.ZhouZ. W.XiaJ. C.Reduced expression of transcription factor AP-2α is associated with gastric adenocarcinoma prognosis201169e2489710.1371/journal.pone.00248972-s2.0-8005313565621966377WagnerA. D.SynN. L. X.MoehlerM.GrotheW.YongW. P.TaiB. C.HoJ.UnverzagtS.Chemotherapy for advanced gastric cancer20178, article CD00406410.1002/14651858.CD004064.pub42-s2.0-84941982813Abdel-RahmanO.Immune checkpoints aberrations and gastric cancer; assessment of prognostic value and evaluation of therapeutic potentials201697657110.1016/j.critrevonc.2015.08.0152-s2.0-84952639687StricklerJ. H.EGFR amplification as a target in gastroesophageal adenocarcinoma: do anti-EGFR therapies deserve a second chance?20188667968110.1158/2159-8290.CD-18-01912-s2.0-8504827120029858225OoiW. F.NargundA. M.LimK. J.ZhangS.XingM.MandoliA.LimJ. Q.HoS. W. T.GuoY.YaoX.LinS. J.NandiT.XuC.OngX.LeeM.TanA. L. K.LamY. N.TeoJ. X.KanedaA.WhiteK. P.LimW. K.RozenS. G.TehB. T.LiS.SkanderupA. J.TanP.Integrated paired-end enhancer profiling and whole-genome sequencing reveals recurrentCCNE1andIGF2enhancer hijacking in primary gastric adenocarcinoma20206961039105210.1136/gutjnl-2018-3176122-s2.0-8507261837531542774SalemM. E.PucciniA.XiuJ.RaghavanD.LenzH. J.KornW. M.ShieldsA. F.PhilipP. A.MarshallJ. L.GoldbergR. M.Comparative molecular analyses of esophageal squamous cell carcinoma, esophageal adenocarcinoma, and gastric adenocarcinoma201823111319132710.1634/theoncologist.2018-01432-s2.0-8505628538329866946ZhuoW.LiuY.LiS.GuoD.SunQ.JinJ.RaoX.LiM.SunM.JiangM.XuY.TengL.JinY.SiJ.LiuW.KangY.ZhouT.Long Noncoding RNA _GMAN_ , Up-regulated in Gastric Cancer Tissues, Is Associated With Metastasis in Patients and Promotes Translation of Ephrin A1 by Competitively Binding _GMAN-AS_20191563676691.e1110.1053/j.gastro.2018.10.0542-s2.0-8506101401830445010CaoY.XiongJ. B.ZhangG. Y.LiuY.JieZ. G.LiZ. R.Long noncoding RNA UCA1 regulates PRL-3 expression by sponging microRNA-495 to promote the progression of gastric cancer20201985386410.1016/j.omtn.2019.10.02031982772LiZ. T.ZhangX.WangD. W.XuJ.KouK. J.WangZ. W.YongG.LiangD. S.SunX. Y.Overexpressed lncRNA GATA6-AS1 inhibits LNM and EMT via FZD4 through the Wnt/β-catenin signaling pathway in GC20201982784010.1016/j.omtn.2019.09.03431981860GuoW.LvP.LiuS.XuF.GuoY.ShenS.LiangJ.KuangG.DongZ.Aberrant methylation-mediated downregulation of long noncoding RNA C5orf66-AS1 promotes the development of gastric cardia adenocarcinoma201857785486510.1002/mc.228062-s2.0-8504484998629566283LiD.WangJ.ZhangM.HuX.SheJ.QiuX.ZhangX.XuL.LiuY.QinS.LncRNA MAGI2-AS3 is regulated by BRD4 and promotes gastric cancer progression via maintaining ZEB1 overexpression by sponging miR-141/200a20201910912310.1016/j.omtn.2019.11.00331837602ShuaiY.MaZ.LuJ.FengJ.LncRNA SNHG15: a new budding star in human cancers2020531, article e1271610.1111/cpr.1271631774607ZhengR.LiangJ.LuJ.LiS.ZhangG.WangX.LiuM.WangW.ChuH.TaoG.ZhaoQ.WangM.duM.QiangF.ZhangZ.Genome-wide long non-coding RNAs identified a panel of novel plasma biomarkers for gastric cancer diagnosis201922473174110.1007/s10120-018-00915-72-s2.0-8505959750130603909UlitskyI.BartelD. P.lincRNAs: genomics, evolution, and mechanisms20131541264610.1016/j.cell.2013.06.0202-s2.0-8487998778923827673ZhangC.PengG.Non-coding RNAs: an emerging player in DNA damage response201576320221110.1016/j.mrrev.2014.11.0032-s2.0-84925220515