Home A four-lncRNA signature for predicting prognosis of recurrence patients with gastric cancer
Article Open Access

A four-lncRNA signature for predicting prognosis of recurrence patients with gastric cancer

  • Qiang Chen , Zunqi Hu , Xin Zhang , Ziran Wei , Hongbing Fu , DeJun Yang EMAIL logo and Qingping Cai EMAIL logo
Published/Copyright: April 3, 2021

Abstract

Purpose

This study aimed to develop a multi-long noncoding RNA (lncRNA) signature for the prediction of gastric cancer (GC) based on differential gene expression between recurrence and nonrecurrence patients.

Methods

By repurposing microarray expression profiles of RNAs from The Cancer Genome Atlas (TCGA), we performed differential expression analysis between recurrence and nonrecurrence patients. A prognostic risk prediction model was constructed based on data from TCGA database, and its reliability was validated using data from Gene Expression Omnibus database. Furthermore, the lncRNA-associated competing endogenous RNA (ceRNA) network was constructed, namely, DIANA-LncBasev2 and starBase database.

Results

We identified 363 differentially expressed RNAs (317 mRNAs, 18 lncRNAs, and 28 microRNAs [miRNAs]). Principal component analysis showed that the seven-feature lncRNAs screened by support vector machine–recursive feature elimination algorithm was more informative for predicting recurrence of GC in comparison with the eight-feature lncRNAs screened by random forest–out-of-bag algorithm. Four of the seven-feature lncRNAs including LINC00843, SNHG3, C21orf62-AS1, and MIR99AHG were chosen to develop a four-lncRNA risk score model. This risk score model was able to distinguish patients with high and low risk of recurrence, and was tested in two independent validation sets. The ceRNA network of this four-lncRNA signature included 10 miRNAs and 178 mRNAs. The mRNAs significantly related to the Wnt-signaling pathway and relevant biological processes.

Conclusion

A useful four-lncRNA signature recurrence was established to distinguish GC patients with high and low risk of recurrence. Regulating the relevant miRNAs and Wnt pathway might partly affect GC metastasisby.

1 Introduction

Gastric cancer (GC) is the fifth common cancer characterized by high incidence and mortality [1]. Current management strategies for GC mainly include surgical resection guided by endoscopic detection and chemotherapy or chemoradiotherapy as adjuvant therapy [2]. Despite advancements in treatments, GC patients have unsatisfactory prognosis [3]. Conventional tumor node metastasis (TNM) stage system is inadequate for outcome prediction of GC patients [4], and novel prognostic biomarkers are complementary and necessary for identifying potential high-risk GC patients and contribute to better outcome in GC patients.

Long noncoding RNAs (lncRNAs) play important regulatory roles in cancer biogenesis. A large number of lncRNAs have been shown to be dysregulated in GC, participate in gastric tumorigenesis and progression through interacting with DNA, RNA, and proteins [5,6]. Moreover, lncRNAs are associated with the prognosis of GC patients, and several lncRNAs-based signatures have been reported for outcome prediction [7,8,9]. High incidence of recurrence following curative resection is a primary cause of undesirable prognosis in patients with advanced GC [10,11]. Therefore, we utilized gene profiling data of recurrence and nonrecurrence GC patients from The Cancer Genome Atlas (TCGA) to explore aberrantly expressed lncRNAs associated with recurrence and develop an lncRNA-based signature for prognosis stratification of GC patients. Two validation data sets from Gene Expression Omnibus (GEO) were used to confirm the prognostic ability of the signature. Furthermore, we unraveled recurrence-related differentially expressed microRNAs (miRNAs) and mRNAs as well as to study the relationships between signature lncRNAs, miRNAs, and potential targeted miRNAs by constructing competing endogenous RNA (ceRNA) network, thereby providing insights into the regulatory mechanisms of these signature lncRNAs in GC.

2 Methods

2.1 Data and preprocessing

We obtained RNA sequencing data (including mRNA and lncRNA) of 407 GC samples and miRNA sequencing data of 477 samples with the corresponding clinical information from the publicly accessible TCGA database (https://gdc-portal.nci.nih.gov/) based on Illumina HiSeq 2000 RNA Sequencing platform. A total of 287 samples with paired RNA and miRNA data as well as the corresponding overall survival (OS) information were selected as the training set. The detailed clinical information of samples in TCGA data sets is shown in Table S1.

Meanwhile, we launched a search in GEO database (https://www.ncbi.nlm.nih.gov/geo/) using “gastric cancer,” “stomach cancer,” and “homo sapiens” as the key words. Potential data sets were selected when the following criteria were met: the number of total samples ≥250 and the number of GC samples with corresponding clinical data ≥200. As a result, two data sets including GSE26253 and GSE62254 were chosen. GSE26253 [12] (GPL8432 Illumina HumanRef-8 WG-DASL v3.0 platform) contained gene expression data of 432 GC tissue samples with clinical information available (validation set 1), while GSE62254 [2] (Affymetrix Human Genome U133 Plus 2.0 Array platform) included gene profiling data of 300 GC tissue samples, among which 282 samples had the corresponding clinical data (validation set 2). The detailed clinical information of samples in GSE26253 and GSE62254 data sets is shown in Tables S2 and S3.

In order to supplement the information of targeted molecular therapy, we collected information about the targeted drugs from databases TCGA, xena.ucsc, and cbioporta. Unfortunately, we only got the information on whether the patients received “targeted treatment” or not; however, there were no specific information of targeted drug molecules. Meanwhile, the detailed clinical characteristics of the patients were analyzed and shown in Table S4. In order to exclude the effect of therapeutic schedule on the screening of the present lncRNAs, correlation analysis was used to eliminate the doubt in this aspect, and the results showed that four important lncRNAs were not related to radiotherapy, chemotherapy, and targeted therapy. Their detailed information is shown in Table S5.

Raw data were standardized as previously described by Chaudhary et al. [13]. HUGO Gene Nomenclature Committee (HGNC) [14] repository (http://www.genenames.org/) enrolled 4,313 lncRNAs, 19,197 protein-coding genes, and 1,914 miRNAs. Using the HGNC database, we annotated 13,105 mRNAs, 1,051 lncRNAs, and 413 miRNAs from the abovementioned data sets, according to RefSeq ID information.

2.2 Differential expression analysis between recurrence and nonrecurrence GC samples

We screened differentially expressed RNAs (DERs) including lncRNAs, mRNAs, and miRNAs between recurrence and nonrecurrence samples in the training set, using FDR <0.05 and |log2 FC| > 0.263 as a selection threshold. The identified significant lncRNAs, miRNAs, and mRNAs consequently underwent two-way hierarchical clustering analysis based on centered Pearson correlation algorithm [15].

2.3 Prognostic model building and validation

In order to identify recurrence-related feature lncRNAs from the preselected differentially expressed lncRNAs, we employed and compared the support vector machine–recursive feature elimination (SVM–RFE) and random forest–out-of-bag (RF–OOB) algorithms for performance. With regard to SVM–RFE [16] algorithm (100-fold cross validation), the lncRNAs subset with the best accuracy was chosen to be the signature lncRNAs. Using RF–OOB [17] algorithm, the subset of lncRNAs with the minimal value of OOB error was selected to be the optimal feature lncRNAs. Principal component analysis (PCA) method was applied to compare the performances of the two algorithms. The feature lncRNAs identified by the superior approach was applied in further analysis.

Using training set, we did univariable Cox regression analysis on the identified recurrence-related signature lncRNAs. The significant lncRNAs were further included in multivariable Cox regression analysis. log-rank p value <0.05 defined the significance. The lncRNAs that were independent prognostic indicators were selected to construct a prognostic score formula as following:

Risk score (RS) = β lncRNAn × Exp lncRNAn

where here β lncRNA suggests multivariable Cox regression coefficient of lncRNAn; ExplncRNAn indicates the expression value of lncRNAn.

We calculated the risk score for each patient in the training set based on the formula. Patients in the training set were categorized by the median risk score into a high-risk group and a low-risk group. Similarly, the lncRNAs-based risk score formula was applied to distinguish patients in two validation sets. Recurrence-free survival (RFS) time of two risk groups was compared using Kaplan–Meier method and log-rank test. Sensitivity and specificity of the risk score model were assessed using the ROC curve.

Using the training set data, we conducted uni- and multivariable Cox regression analysis to evaluate the association of clinical features and risk score model status with RFS time of patients. Integrating independent prognostic clinical features with risk score status, nomogram was built; and the calibration curves were plotted to ascertain its predictive performance.

2.4 Statistical analysis

A wide range of packages in R software (version 3.4.1) were utilized for bioinformatics and statistical analyses of our study: limma package (version 3.34.7) for differential expression analysis; pheatmap package (version 1.0.8) for two-way hierarchical clustering analysis; e1071 [18] (version1.7-1, https://cran.r-project.org/web/packages/e1071); caret [19] (version 6.0-76, https://cran.r-project.org/web/packages/caret) packages for SVM–RFE method; bootstrap algorithm of randomForest [20] package (https://cran.r-project.org/web/packages/randomForest/) for RF–OOB method; psych [21] package (version 1.8.12) for PCA; survival package (http://bioconductor.org/packages/survivalr/) for uni- and multivariable Cox regression analysis; Kaplan–Meier curves; pROC package (https://cran.r-project.org/web/packages/pROC/index.html) for ROC curve analysis; and rms package (version 5.1-2) for nomogram building.

2.5 Construction of lncRNA-associated ceRNA network

For dissecting the underlying molecular mechanisms of the prognostic signature lncRNAs in GC biology, we designed a three-phase study. Initially, the relations between the signature lncRNAs with the differentially expressed miRNAs were predicted based on DIANA-LncBasev2 [22] database. The negatively correlated lncRNA–miRNA pairs were selected to build an lncRNA–miRNA network.

Second, starBase database (version 2.0, http://starbase.sysu.edu.cn/) was used for the prediction of potential target mRNAs of the miRNAs included in the lncRNA–miRNA network. starBase program compiles prediction results from TargetScan, PicTar, RNA22, PITA, and miRanda [23]. The differentially expressed mRNAs were mapped to the predicted target mRNAs that were enrolled in at least three of the five programs. The miRNA–mRNA pairs with negative correlation were selected to develop an miRNA–mRNA network.

Third, an lncRNA–miRNA–mRNA network was built with predicted lncRNA–miRNA and miRNA–mRNA pairs with negative correlation. Networks were visualized with Cytoscape [24] software (version 3.6.1). All mRNAs in the networks were subject to gene ontology (GO) function and kyoto encyclopedia of genes and genomes (KEGG) pathway enrichment analysis by Enrichr [25] tool (http://amp.pharm.mssm.edu/Enrichr/). P value <0.05 was considered significant.

3 Results

3.1 Identification of differentially expressed mRNAs, lncRNAs, and miRNAs between recurrence and nonrecurrence GC patients

The training set (N = 287) consisted of 61 recurrence patients and 226 nonrecurrence patients. Totally 363 DERs between the recurrent and nonrecurrent patients were identified by performing differential gene expression analysis, which were consisted of 317 mRNAs (115 downregulated and 202 upregulated mRNAs; Table S6), 18 lncRNAs (8 downregulated and 10 upregulated lncRNAs; Table S7), and 28 miRNAs (9 downregulated and 19 upregulated miRNAs; Table S8; Figure 1a and b).

Figure 1 
                  Differentially expressed lncRNAs, miRNAs, and mRNAs of recurrence and nonrecurrence patients. (a) Left image is a scatterplot for differentially expressed RNAs. Lateral axis represents nonrecurrence samples; vertical axis represents recurrence samples. Upregulated, downregulated, and nonregulated RNAs are shown in red, green, and blue spots, respectively. Right image shows ratios of up- and downregulated RNAs. (b) Heatmap displays differentially expressed RNAs in recurrence and nonrecurrence samples.
Figure 1

Differentially expressed lncRNAs, miRNAs, and mRNAs of recurrence and nonrecurrence patients. (a) Left image is a scatterplot for differentially expressed RNAs. Lateral axis represents nonrecurrence samples; vertical axis represents recurrence samples. Upregulated, downregulated, and nonregulated RNAs are shown in red, green, and blue spots, respectively. Right image shows ratios of up- and downregulated RNAs. (b) Heatmap displays differentially expressed RNAs in recurrence and nonrecurrence samples.

3.2 Construction of a four-lncRNA prognostic score risk model

We used SVM–RFE and RF–OOB algorithms to select feature lncRNAs informative of recurrence based on the abovementioned differentially expressed lncRNAs. Using SVF-RFM algorithm, we obtained a combination of eight-feature lncRNAs when achieving the highest accuracy = 0.895 (Figure 2a, Table S9). A set of seven lncRNAs was determined by RF–OOB algorithm (minimal OOB error = 0.173; Figure 2b; Table S10). The two groups of feature lncRNAs were compared by PCA. Regarding the eight-feature lncRNAs identified by SVM–RFE, the top four principal components in PCA approximately account for 80% of the total observed variances (Figure 3a). In contrast, the top three principal components for the seven-feature lncRNAs selected by RF–OOB algorithm could describe almost 80% of the total observed variances (Figure 3b). These results suggested that the seven-feature lncRNAs were more indicative of recurrence than the eight-feature lncRNAs and were thus used in further analysis.

Figure 2 
                  Feature lncRNAs selected by SVM–RFE (a) and RF–OOB (b). The red dashed line shows the number of lncRNAs achieving the highest accurracy (a) or the smallest OOB error (b).
Figure 2

Feature lncRNAs selected by SVM–RFE (a) and RF–OOB (b). The red dashed line shows the number of lncRNAs achieving the highest accurracy (a) or the smallest OOB error (b).

Figure 3 
                  Comparative analysis of SVM–RFE (a) and RF–OOB (b) by PCA approach. Upper images illustrate 3D PCA plot of PC1, PC2, and PC3. Recurrence and nonrecurrence samples are labeled in red triangles and blue balls, respectively. Below images exhibit cumulative contribution of different PCs. The red dashed line indicates cumulative proportion of 80%.
Figure 3

Comparative analysis of SVM–RFE (a) and RF–OOB (b) by PCA approach. Upper images illustrate 3D PCA plot of PC1, PC2, and PC3. Recurrence and nonrecurrence samples are labeled in red triangles and blue balls, respectively. Below images exhibit cumulative contribution of different PCs. The red dashed line indicates cumulative proportion of 80%.

By subjecting the seven-feature lncRNAs to a univariable Cox regression analysis, we found six prognosis-related lncRNAs (p value <0.05), which further underwent multivariable Cox regression analysis. Four lncRNAs (LINC00843, SNHG3, C21orf62-AS1, and MIR99AHG), which were independent predictors of prognosis, were depicted in Table 1. A risk score formula was created based on the expression of the four signature lncRNAs for prognosis prediction as follows:

Risk score = ( 1.0637 ) × Exp LINC 00843 + ( 0.4163 ) × Exp SNHG 3 + ( 1.4731 ) × Exp C 21 orf 62 -AS 1 + ( 0.2341 ) × Exp MIR 99 AHG

Table 1

Characteristics of four independent prognostic lncRNAs

ID Coefficient P value HR 95% CI
LINC00843 −1.0637 6.71 × 10−3 0.345 0.111–0.878
SNHG3 −0.4163 8.97 × 10−3 0.659 0.208–0.967
C21orf62-AS1 1.4731 9.17 × 10−3 4.363 1.788–7.171
MIR99AHG 0.2341 4.96 × 10−2 1.264 1.045–2.477

We calculated the risk score for each patient in the training set and ranked all the patients based on four-lncRNA signature. With the median risk score as the cutoff, the training set was categorized into a high-risk group (N = 144) and a low-risk group (N = 143). Significantly longer RFS time was observed in the low-risk patients relative to the high-risk patients (p value = 1.448 × 10−3, HR = 2.329[1.363–3.979]; Figure 4a). We validated the four-lncRNA signature risk score in two validation sets. As shown in Figure 4b and c, either validation set was divided by the risk score into two risk groups with statistical significance in RFS time (validation set 1: p value = 3.733 × 10−2, HR = 1.368[1.017–1.839]; validation set 2: p value = 1.444 × 10−2, HR = 1.555[1.089–2.221]). The risk score of each sample in the data set TCGA, GSE26253, and GSE62254 is shown in Table S11(1–3), respectively. The AUC values of ROC curves were 0.936, 0.827, and 0.822 for training set, validation set 1, and validation set 2, separately (Figure 4d). These observations demonstrated robust predictive performance of the four-lncRNA signature risk score in GC.

Figure 4 
                  Kaplan–Meier estimates of recurrence-free survival and ROC curves for patients in training set and two validation sets by the four-lncRNA signature. (a–c) Kaplan–Meier curves for training set, validation set 1, and validation set 2. (d) ROC curves and AUC values of three data sets.
Figure 4

Kaplan–Meier estimates of recurrence-free survival and ROC curves for patients in training set and two validation sets by the four-lncRNA signature. (a–c) Kaplan–Meier curves for training set, validation set 1, and validation set 2. (d) ROC curves and AUC values of three data sets.

3.3 Building nomogram based on gender, histologic grade, and risk score model status

By performing univariable Cox regression analyses to assess the relationship of clinical variables and risk score model status with RFS time of patients in the training set, gender (HR = 1.761, 95% CI = 0.967–1.014, p value = 4.92 × 10−2), neoplasm histologic grade (HR = 2.181, 95% CI = 1.232–3.856, p value = 6.19 × 10−3), and the four-lncRNA risk score model status (HR = 2.329, 95% CI = 1.363–3.979, p value = 1.45 × 10−3) were statistically significant (Table 2). Furthermore, multivariable Cox regression analysis incorporating risk score model status with gender and neoplasm histologic grade was performed. Gender (HR = 1.735, 95% CI = 1.179–3.078, p value = 4.93 × 10−2), neoplasm histologic grade (HR = 2.063, 95% CI = 1.156–3.682, p value = 1.43 × 10−2), and risk score model status (HR = 2.059, 95% CI = 1.198–3.537, p value = 8.91 × 10−3) were found to be independent prognostic predictors (Table 2), indicating that prognostic value of the four-lncRNA risk score is independent of other clinical features.

Table 2

Uni- and multivariable Cox regression analysis of clinical features and risk score model status

Clinical characteristics TCGA (N = 287) Univariables Cox Multivariables Cox
HR 95% CI p HR 95% CI p
Age (years, mean ± sd) 65.06 ± 10.65 0.991 0.967–1.014 4.12 × 10−1
Gender (male/female) 181/106 1.761 0.994–3.117 4.92 × 10 −2 1.735 1.179–3.078 4.93 × 10 −2
Pathologic_M (M0/M1/–) 259/14/14 1.149 0.359–3.679 8.15 × 10−1
Pathologic_N (N0/N1/N2/N3/–) 89/77/58/53/10 1.151 0.918–1.440 2.24 × 10−1
Pathologic_T (T1/T2/T3/T4/–) 16/64/125/78/4 0.796 0.599–1.056 1.12 × 10−1
Pathologic_stage (I/II/III/IV/–) 41/95/121/19/11 1.029 0.763–1.387 8.53 × 10−1
Neoplasm histologic grade (G1/G2/G3/–) 5/101/172/9 2.181 1.232–3.856 6.19 × 10 −3 2.063 1.156–3.682 1.43 × 10 −2
Radiation therapy (Yes/No/–) 54/229/4 0.627 0.321–1.225 1.69 × 10−1
Targeted molecular therapy (Yes/No/–) 130/151/6 1.375 0.824–2.293 2.21 × 10−1
Helicobacter pylori infection (Yes/No/–) 15/107/165 0.314 0.119– 2.358 2.35 × 10−1
Residual tumor (R0/R1/R2/R3/–) 239/10/6/32 1.774 0.922–3.413 1.28 × 10−1
Chemotherapy (Yes/No/–) 136/151 1.313 0.772–2.231 3.13 × 10−1
RS model status (High/Low) 143/144 2.329 1.363–3.979 1.45 × 10 −3 2.059 1.198–3.537 8.91 × 10 −3
Recurrence (Yes/No) 61/226
Recurrence-free survival time (months, mean ± sd) 19.70 ± 18.48

Statistically significant data are indicated by bold for significant at p < 0.05.

A composite nomogram on the basis of gender, neoplasm histologic grade, and risk score status was established to predict the probability of 3-year and 5-year RFS (Figure 5a). Calibration plot for goodness of fit of the model exhibited good consistence between the predicted and actual RFS time (Figure 5b).

Figure 5 
                  Nomogram to predict 3-year and 5-year probability of recurrence-free survival using TCGA set. (a) nomogram based on gender, neoplasm grade, and risk score model status. RS model status = four-lncRNA risk score model status. (b) calibration curves of actual and predicted probabilities.
Figure 5

Nomogram to predict 3-year and 5-year probability of recurrence-free survival using TCGA set. (a) nomogram based on gender, neoplasm grade, and risk score model status. RS model status = four-lncRNA risk score model status. (b) calibration curves of actual and predicted probabilities.

3.4 Characterization of ceRNA network of four signature lncRNAs and functional annotation

Using DIANA-LncBasev2, we analyzed the correlations between the four signature lncRNAs and the differentially expressed miRNAs. As a result, we procured 10 negatively correlated lncRNA–miRNA pairs and constructed an lncRNA–miRNA network. There were 3 lncRNAs (upregulated MIR99AHG, downregulated LINC00843, and SNHG3) and 10 miRNAs (five upregulated and five downregulated miRNAs) in the network (Figure 6). As to the 10 miRNAs (miR-7, miR-552, miR-4676, miR-1304, miR-2110, miR-216a, miR-205, miR-487a, miR-551b, and miR-34b), we predicted potential target mRNAs by using starBase, which were further mapped by the 321 differentially expressed mRNAs. The resulting 351 miRNA–mRNA pairs with negative correlation were used to develop a miRNA–mRNA network. Figure 7 exhibited this network of 10 miRNAs and 178 mRNAs (54 downregulated and 124 upregulated mRNAs). Finally, the preselected negatively correlated lncRNA–miRNA pairs and miRNA–mRNA pairs were integrated into a ceRNA regulatory network. Three signature lncRNAs, 10 miRNAs, and 178 mRNAs were observed in the network (Figure 8). Their detailed information are shown in Tables S12 and S13.

Figure 6 
                  Signature lncRNA–miRNA network. Diamonds and triangles stand for lncRNAs and miRNAs, respectively.
Figure 6

Signature lncRNA–miRNA network. Diamonds and triangles stand for lncRNAs and miRNAs, respectively.

Figure 7 
                  The miRNA–mRNA network. Balls and triangles denote mRNAs and miRNAs, respectively.
Figure 7

The miRNA–mRNA network. Balls and triangles denote mRNAs and miRNAs, respectively.

Figure 8 
                  A ceRNA network of signature lncRNAs, miRNAs, and mRNAs. Diamonds, triangles, and balls represent signature lncRNAs, miRNAs, and mRNAs, separately. Red and black links suggest lncRNA–miRNA and miRNA–mRNA interactions, respectively.
Figure 8

A ceRNA network of signature lncRNAs, miRNAs, and mRNAs. Diamonds, triangles, and balls represent signature lncRNAs, miRNAs, and mRNAs, separately. Red and black links suggest lncRNA–miRNA and miRNA–mRNA interactions, respectively.

The GO function and KEGG pathway-enrichment analysis for the mRNAs in the ceRNA network was conducted. These mRNAs were significantly enriched in 19 GO terms predominately related to Wnt-signaling pathway, cell proliferation, migration, and other biological processes (Table 3). Five significant KEGG-signaling pathways were noted, including Wnt-signaling pathway, cell adhesion molecules, GC, regulation of lipolysis in adipocytes, and steroid hormone biosynthesis pathways (Table 3). These results revealed that four signature lncRNAs-related ceRNA regulation participated in a variety of biological functions and signaling pathways in GC.

Table 3

Results of GO function and KEGG pathway enrichment analysis

Type Term Count of significantly enriched genes P value
GO biology process Regulation of canonical Wnt-signaling pathway 10 2.27 × 10−5
Negative regulation of canonical Wnt-signaling pathway 7 3.82 × 10−4
Regulation of cell migration 10 5.67 × 10−4
Neuron projection morphogenesis 7 6.77 × 10−4
Regulation of cell proliferation 16 9.86 × 10−4
Negative regulation of Wnt-signaling pathway 7 9.91 × 10−4
Cell morphogenesis involved in neuron differentiation 5 1.84 × 10−3
Positive regulation of cell growth 5 1.92 × 10−3
Wnt-signaling pathway 5 3.04 × 10−3
Positive regulation of cell motility 6 5.55 × 10−3
Positive regulation of cellular process 11 7.01 × 10−3
Skeletal system development 5 1.02 × 10−2
Positive regulation of cell proliferation 9 1.41 × 10−2
Positive regulation of cell migration 6 1.47 × 10−2
Extracellular matrix organization 6 1.72 × 10−2
Positive regulation of multicellular organismal process 5 3.54 × 10−2
Negative regulation of cell proliferation 7 4.47 × 10−2
Proteolysis 6 4.69 × 10−2
KEGG pathway Wnt-signaling pathway 7 5.43 × 10−4
Cell adhesion molecules (CAMs) 4 4.08 × 10−2
Gastric cancer 4 4.43 × 10−2
Regulation of lipolysis in adipocytes 2 4.86 × 10−2
Steroid hormone biosynthesis 2 4.98 × 10−2

4 Discussion

Although a growing number of prognostic lncRNAs for GC have been uncovered [26], some limitations are present, such as small number of lncRNAs, limited sample size, and insufficient validation. GC patients often experience recurrence following surgical resection [27]. To develop a recurrence-related multi-lncRNA signature for the prediction of RFS in GC patients, we repurposed the existing microarray data downloaded from TCGA to profile lncRNAs, miRNAs, and mRNAs in GC patients. A total of 363 DERs between recurrence and nonrecurrence patients were obtained, comprising 18 lncRNAs, 317 mRNAs, and 28 miRNAs.

Another highlight of this study was that we applied and compared SVM–RFE and RF–OOB to identify feature lncRNAs’ most informative of recurrence from the differentially expressed lncRNAs. SVM–RFE is considered as an efficient method to select informative genes for cancer classification, in which all the features are listed based on some score function together with the removal of the features with the lowest scores [16]. One shortcoming of this method is that it only aims to identify the optimal combination for classification. RF is a highly data-adaptive classification tool based on decision trees that is especially suitable for high-dimension genomic data analysis with OBB error to assess the predictive performance of RF [28]. According to the results of PCA, the RF-OBB-based seven feature lncRNAs was more informative of recurrence in comparison with the SVM–RFE-based eight feature lncRNAs in the current study. From the seven feature lncRNAs, we identified a four-lncRNA signature that was significantly associated with patients’ RFS and had independent prognostic value. Moreover, a four-lncRNA risk score for outcome prediction was developed using TCGA set and validated using two GEO data sets. Our results indicated that this four-lncRNA risk score could successfully distinguish GC patients with high risk from GC patients with low risk. In addition, our study showed that predictive performance of the four lncRNAs-based risk score was independent of other clinical variables. Potential clinical application of this risk score would be beneficial to improving individualized treatment decision-making for GC patients.

The four identified prognostic lncRNAs were LINC00843, SNHG3, C21orf62-AS1, and MIR99AHG. lncRNA SNHG3 is reported to be implicated in development of various types of cancers, such as colorectal cancer and ovarian cancer [29,30]. SNHG3 is overexpressed in hepatocellular carcinoma (HCC), showing correlation with the survival of HCC patients [31]. SNHG3 is upregulated lncRNA in GC [32]. However, biological functions and prognostic value of SNHG3 in GC remain elusive. lncRNA C21orf62-AS1 is abnormally expressed in chromophobe renal cell carcinoma, correlating with OS of patients [33]. Additionally, C21orf62-AS is stimulated by interferon-beta in patients with multiple sclerosis [34]. lncRNA MIR99AHG is significantly downregulated in colorectal cancer [35]. There is evidence that MIR99AHG is positively relevant to OS of patients with lung squamous cell carcinoma [36]. There is little information concerning lncRNA LINC00843. As far as we know, this is the first time that these lncRNAs are found to be prognostic biomarkers for GC recurrence.

It has been demonstrated that ceRNAs act as key regulators among different RNA transcripts and lncRNAs sponges miRNAs, thereby regulating the expression of targeted mRNAs [37]. The present study established a ceRNA network comprising three signature lncRNAs (LINC00843, SNHG3, and MIR99AHG), 10 miRNAs, and 178 target mRNAs by bioinformatics prediction. Among the 10 miRNAs, miR-7 plays an antimetastatic role in GC through targeting insulin-like growth factor-1 receptor [38]. The miR-487a strengthens cell proliferation and suppresses cell apoptosis, driving GC progression by targeting T-cell intracellular antigen-1 [39]. Furthermore, three miRNAs including miR-216a, miR-205, and miR-551b are related to epithelial–mesenchymal transition (EMT) and metastasis of GC [4042]. These findings are suggestive of implication of the three signature lncRNAs in metastasis of GC through influencing relevant miRNAs and target mRNAs. Noticeably, results of GO and KEGG pathway enrichment analysis demonstrated that mRNAs in the ceRNA network were associated with several Wnt-signaling pathway-related biological processes and KEGG-signaling pathway. Wnt-signaling pathway plays a fundamental role in progression and metastasis of GC, participating in regulating GC cell growth and apoptosis [43]. Involvement of Wnt-signaling pathway in EMT is controlled by miRNAs [44]. Besides, Wnt-signaling pathway is reported to mediate the oncogenic role of miR-552 by Dachshund family transcription factor 1 in colorectal cancer [45]. It can be speculated that Wnt-signaling pathway may partly mediate the effect of these signature lncRNAs on progression and metastasis of GC through miRNAs. Our study only contains results based on gene mining approaches. Clinical experimental studies and large prospective studies are necessary to verify our findings.

In conclusion, we generated a recurrence-related four-lncRNA signature predictive of individual mortality risk of DFS in GC patients. Prognostic capability of the lncRNAs-based signature had been successfully validated using two independent data sets and showed independence of other clinical features. The four-lncRNA signature functionally involved several metastasis-related miRNAs, their targeted mRNAs, and Wnt-signaling pathway. This study suggested potential prognostic biomarkers and therapeutic targets for recurrence GC and provided novel insights into the underlying mechanisms of GC progression. Validation of our findings and investigation on functional mechanisms warrant future studies.

Abbreviations

GC

Gastric cancer

lncRNA

Long non-coding RNA

TCGA

The Cancer Genome Atlas

GEO

Gene Expression Omnibus

ceRNA

competing endogenous RNAs

miRNA

microRNA

OS

overall survival

HGNC

HUGO Gene Nomenclature Committee

DERs

differentially expressed RNAs

SVM

Support Vector Machine

RFE

Recursive Feature Elimination

RF

Random Forest

OOB

Out Of Bag

PCA

Principal Component Analysis

RFS

Recurrence-free survival

GO

Gene ontology

KEGG

kyoto encyclopedia of genes and genomes

CAMs

cell adhesion molecules

HCC

hepatocellular carcinoma

EMT

epithelial-mesenchymal transition

DACH

Dachshund family transcription factor


Co-first authors.

tel: + 86-13818889536

Research funding

This study was supported by Science and Technology Commission of Shanghai Municipality (Grant 18411952800), Shanghai Changzheng Hospital Pyramid Talent Project (Grant No. 201711015), and the Natural Science Foundation of Shanghai (Grant 17ZR1439300).

  1. Conflict of interest: Authors state no conflict of interest.

  2. Data availability statements: All data generated or analyzed during this study are included in this published article and its supplementary information files.

References

[1] 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.10.3322/caac.21492Search in Google Scholar

[2] Cristescu R, Lee J, Nebozhyn M, Kim K-M, Ting JC, Wong SS, et al. Molecular analysis of gastric cancer identifies subtypes associated with distinct clinical outcomes. Nat Med. 2015;21(5):449–56.10.1038/nm.3850Search in Google Scholar

[3] Chen W, Zheng R, Baade PD, Zhang S, Zeng H, Bray F, et al. Cancer statistics in China, 2015. CA Cancer J Clin. 2016;66(2):115–32.10.3322/caac.21338Search in Google Scholar

[4] Marano L, Boccardi V, Braccio B, Esposito G, Grassia M, Petrillo M, et al. Comparison of the 6th and 7th editions of the AJCC/UICC TNM staging system for gastric cancer focusing on the “N” parameter-related survival: the monoinstitutional NodUs Italian study. World J Surg Oncol. 2015;13(1):215.10.1186/s12957-015-0633-3Search in Google Scholar

[5] Li T, Mo X, Fu L, Xiao B, Guo J. Molecular mechanisms of long noncoding RNAs on gastric cancer. OncoTarget. 2016;7(8):8601.10.18632/oncotarget.6926Search in Google Scholar

[6] Sun M, Nie F, Wang Z, De W. Involvement of lncRNA dysregulation in gastric cancer. Histol Histopathol. 2016;31(1):33–9.Search in Google Scholar

[7] Zhu X, Tian X, Yu C, Shen C, Yan T, Hong J, et al. A long non-coding RNA signature to improve prognosis prediction of gastric cancer. Mol Cancer. 2016;15(1):60.10.1186/s12943-016-0544-0Search in Google Scholar

[8] Miao Y, Sui J, Xu S-Y, Liang G-Y, Pu Y-P, Yin L-H. Comprehensive analysis of a novel four-lncRNA signature as a prognostic biomarker for human gastric cancer. OncoTarget. 2017;8(43):75007.10.18632/oncotarget.20496Search in Google Scholar

[9] Fan ZY, Liu W, Chao Y, Zhu ZL, Wei X, Li JF, et al. Identification of a five-lncRNA signature for the diagnosis and prognosis of gastric cancer. Tumor Biol. 2016;37(10):13265–77.10.1007/s13277-016-5185-9Search in Google Scholar

[10] Maehara Y, Hasuda S, Koga T, Tokunaga E, Kakeji Y, Sugimachi K. Postoperative outcome and sites of recurrence in patients following curative resection of gastric cancer. Br J Surg. 2000;87(3):353–7.10.1046/j.1365-2168.2000.01358.xSearch in Google Scholar

[11] Van Cutsem E, Sagaert X, Topal B, Haustermans K, Prenen H. Gastric cancer. Lancet. 2016;388(10060):2654–64.10.1016/S0140-6736(16)30354-3Search in Google Scholar

[12] Lee J, Sohn I, Do IG, Kim KM, Park SH, Park JO, et al. Nanostring-based multigene assay to predict recurrence for gastric cancer patients after surgery. PLoS One. 2014;9(3):e90133.10.1371/journal.pone.0090133Search in Google Scholar PubMed PubMed Central

[13] Chaudhary K, Poirion OB, Lu L, Garmire LX. Deep learning-based multi-omics integration robustly predicts survival in liver cancer. Clin Cancer Res. 2018;24(6):1248–59.10.1158/1078-0432.CCR-17-0853Search in Google Scholar PubMed PubMed Central

[14] Wright MW. A short guide to long non-coding RNA gene nomenclature. Hum Genomics. 2014;8:7.10.1186/1479-7364-8-7Search in Google Scholar PubMed PubMed Central

[15] 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.10.1073/pnas.95.25.14863Search in Google Scholar PubMed PubMed Central

[16] Samb ML, Camara F, Ndiaye S, Slimani Y, Esseghir MA. A novel RFE-SVM-based feature selection approach for classification. Int J Adv Sci Technol. 2012;43(1):27–36.Search in Google Scholar

[17] Zapf A, Brunner E, Konietschke F. A wild bootstrap approach for the selection of biomarkers in early diagnostic trials. BMC Med Res Methodol. 2015;15:43.10.1186/s12874-015-0025-ySearch in Google Scholar PubMed PubMed Central

[18] Wang Q, Liu X. Screening of feature genes in distinguishing different types of breast cancer using support vector machine. OncoTargets Ther. 2015;8:2311–7.10.2147/OTT.S85271Search in Google Scholar PubMed PubMed Central

[19] 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.10.1002/mp.12967Search in Google Scholar PubMed PubMed Central

[20] Tolosi L, Lengauer T. Classification with correlated features: unreliability of feature ranking and solutions. Bioinformatics (Oxford, Engl). 2011;27(14):1986–94.10.1093/bioinformatics/btr300Search in Google Scholar PubMed

[21] Shan S, Chen W, Jia JD. Transcriptome analysis revealed a highly connected gene module associated with cirrhosis to hepatocellular carcinoma development. Front Genet. 2019;10:305.10.3389/fgene.2019.00305Search in Google Scholar PubMed PubMed Central

[22] Paraskevopoulou MD, Vlachos IS, Karagkouni D, Georgakilas G, Kanellos I, Vergoulis T, et al. DIANA-lncBase v2: indexing microRNA targets on non-coding transcripts. Nucleic Acids Res. 2016;44(D1):D231–8.10.1093/nar/gkv1270Search in Google Scholar PubMed PubMed Central

[23] Ou M, Zhang X, Dai Y, Gao J, Zhu M, Yang X, et al. Identification of potential microRNA – target pairs associated with osteopetrosis by deep sequencing, iTRAQ proteomics and bioinformatics. Eur J Hum Genet. 2014;22(5):625.10.1038/ejhg.2013.221Search in Google Scholar PubMed PubMed Central

[24] Kohl M, Wiese S, Warscheid B. Cytoscape: software for visualization and analysis of biological networks. Methods Mol Biol. 2011;696:291–303.10.1007/978-1-60761-987-1_18Search in Google Scholar PubMed

[25] Chen EY, Tan CM, Kou Y, Duan Q, Wang Z, Meirelles GV, et al. Enrichr: interactive and collaborative HTML5 gene list enrichment analysis tool. BMC Bioinform. 2013;14:128.10.1186/1471-2105-14-128Search in Google Scholar PubMed PubMed Central

[26] Tian X, Zhu X, Yan T, Yu C, Shen C, Hong J, et al. Differentially expressed lncRNAs in gastric cancer patients: a potential biomarker for gastric cancer prognosis. J Cancer. 2017;8(13):2575.10.7150/jca.19980Search in Google Scholar PubMed PubMed Central

[27] Smyth E, Verheij M, Allum W, Cunningham D, Cervantes A, Arnold D. Gastric cancer: ESMO Clinical Practice Guidelines for diagnosis, treatment and follow-up. Ann Oncol. 2016;27(suppl_5):v38–49.10.1093/annonc/mdw350Search in Google Scholar PubMed

[28] Chen X, Ishwaran H. Random forests for genomic data analysis. Genomics. 2012;99(6):323–9.10.1016/j.ygeno.2012.04.003Search in Google Scholar PubMed PubMed Central

[29] Huang W, Tian Y, Dong S, Cha Y, Li J, Guo X, et al. The long non-coding RNA SNHG3 functions as a competing endogenous RNA to promote malignant development of colorectal cancer. Oncol Rep. 2017;38(3):1402–10.10.3892/or.2017.5837Search in Google Scholar PubMed PubMed Central

[30] Li N, Zhan X, Zhan X. The lncRNA SNHG3 regulates energy metabolism of ovarian cancer by an analysis of mitochondrial proteomes. Gynecol Oncol. 2018;150(2):343–54.10.1016/j.ygyno.2018.06.013Search in Google Scholar PubMed

[31] Zhang T, Cao C, Wu D, Liu L. SNHG3 correlates with malignant status and poor prognosis in hepatocellular carcinoma. Tumor Biol. 2016;37(2):2379–85.10.1007/s13277-015-4052-4Search in Google Scholar PubMed

[32] Cao W-J, Wu H-L, He B-S, Zhang Y-S, Zhang Z-Y. Analysis of long non-coding RNA expression profiles in gastric cancer. World J Gastroenterol. 2013;19(23):3658.10.3748/wjg.v19.i23.3658Search in Google Scholar PubMed PubMed Central

[33] He H-T, Xu M, Kuang Y, Han X-Y, Wang M-Q, Yang Q. Biomarker and competing endogenous RNA potential of tumor-specific long noncoding RNA in chromophobe renal cell carcinoma. OncoTargets Ther. 2016;9:6399.10.2147/OTT.S116392Search in Google Scholar PubMed PubMed Central

[34] Satoh Ji, Takitani M, Miyoshi J, Kino Y. RNA-Seq data analysis identifies the comprehensive profile of in vivo interferon-β-stimulated genes in multiple sclerosis. Clin Exp Neuroimmunol. 2016;7(1):39–51.10.1111/cen3.12268Search in Google Scholar

[35] Huang F, Wen C, Zhuansun Y, Huang L, Chen W, Yang X, et al. A novel long noncoding RNA OECC promotes colorectal cancer development and is negatively regulated by miR-143-3p. Biochem Biophys Res Commun. 2018;503(4):2949–55.10.1016/j.bbrc.2018.08.075Search in Google Scholar PubMed

[36] Sui J, Xu S-Y, Han J, Yang S-R, Li C-Y, Yin L-H, et al. Integrated analysis of competing endogenous RNA network revealing lncRNAs as potential prognostic biomarkers in human lung squamous cell carcinoma. OncoTarget. 2017;8(39):65997.10.18632/oncotarget.19627Search in Google Scholar PubMed PubMed Central

[37] Wang Y, Hou J, He D, Sun M, Zhang P, Yu Y, et al. The emerging function and mechanism of ceRNAs in cancer. Trends Genet. 2016;32(4):211–24.10.1016/j.tig.2016.02.001Search in Google Scholar PubMed PubMed Central

[38] Zhao X, Dou W, He L, Liang S, Tie J, Liu C, et al. MicroRNA-7 functions as an anti-metastatic microRNA in gastric cancer by targeting insulin-like growth factor-1 receptor. Oncogene. 2013;32(11):1363.10.1038/onc.2012.156Search in Google Scholar PubMed

[39] Yang X, Wang M, Lin B, Yao D, Li J, Tang X, et al. miR-487a promotes progression of gastric cancer by targeting TIA1. Biochimie. 2018;154:119–26.10.1016/j.biochi.2018.08.006Search in Google Scholar PubMed

[40] Tao Y, Yang S, Wu Y, Fang X, Wang Y, Song Y, et al. MicroRNA-216a inhibits the metastasis of gastric cancer cells by targeting JAK2/STAT3-mediated EMT process. OncoTarget. 2017;8(51):88870.10.18632/oncotarget.21488Search in Google Scholar PubMed PubMed Central

[41] Xu C, Li Me, Zhang L, Bi Y, Wang P, Li J, et al. MicroRNA-205 suppresses the invasion and epithelial-mesenchymal transition of human gastric cancer cells. Mol Med Rep. 2016;13(6):4767–73.10.3892/mmr.2016.5118Search in Google Scholar PubMed

[42] Song G, Zhang H, Chen C, Gong L, Chen B, Zhao S, et al. miR-551b regulates epithelial-mesenchymal transition and metastasis of gastric cancer by inhibiting ERBB4 expression. OncoTarget. 2017;8(28):45725.10.18632/oncotarget.17392Search in Google Scholar PubMed PubMed Central

[43] Zheng R, Deng Q, Liu Y, Zhao P. Curcumin inhibits gastric carcinoma cell growth and induces apoptosis by suppressing the Wnt/β-catenin signaling pathway. Med Sci Monit Int Med J Exp Clin Res. 2017;23:163.10.12659/MSM.902711Search in Google Scholar

[44] Ghahhari NM, Babashah S. Interplay between microRNAs and Wnt/β-catenin signalling pathway regulates epithelial–mesenchymal transition in cancer. Eur J Cancer. 2015;51(12):1638–49.10.1016/j.ejca.2015.04.021Search in Google Scholar PubMed

[45] Cao J, Yan XR, Liu T, Han XB, Yu JJ, Liu SH, et al. MicroRNA-552 promotes tumor cell proliferation and migration by directly targeting DACH1 via the Wnt/β‑catenin signaling pathway in colorectal cancer. Oncol Lett. 2017;14(3):3795–802.10.3892/ol.2017.6600Search in Google Scholar PubMed PubMed Central

Received: 2020-08-22
Revised: 2021-01-12
Accepted: 2021-01-27
Published Online: 2021-04-03

© 2021 Qiang Chen et al., published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Articles in the same Issue

  1. Research Articles
  2. Identification of ZG16B as a prognostic biomarker in breast cancer
  3. Behçet’s disease with latent Mycobacterium tuberculosis infection
  4. Erratum
  5. Erratum to “Suffering from Cerebral Small Vessel Disease with and without Metabolic Syndrome”
  6. Research Articles
  7. GPR37 promotes the malignancy of lung adenocarcinoma via TGF-β/Smad pathway
  8. Expression and role of ABIN1 in sepsis: In vitro and in vivo studies
  9. Additional baricitinib loading dose improves clinical outcome in COVID-19
  10. The co-treatment of rosuvastatin with dapagliflozin synergistically inhibited apoptosis via activating the PI3K/AKt/mTOR signaling pathway in myocardial ischemia/reperfusion injury rats
  11. SLC12A8 plays a key role in bladder cancer progression and EMT
  12. LncRNA ATXN8OS enhances tamoxifen resistance in breast cancer
  13. Case Report
  14. Serratia marcescens as a cause of unfavorable outcome in the twin pregnancy
  15. Spleno-adrenal fusion mimicking an adrenal metastasis of a renal cell carcinoma: A case report and embryological background
  16. Research Articles
  17. TRIM25 contributes to the malignancy of acute myeloid leukemia and is negatively regulated by microRNA-137
  18. CircRNA circ_0004370 promotes cell proliferation, migration, and invasion and inhibits cell apoptosis of esophageal cancer via miR-1301-3p/COL1A1 axis
  19. LncRNA XIST regulates atherosclerosis progression in ox-LDL-induced HUVECs
  20. Potential role of IFN-γ and IL-5 in sepsis prediction of preterm neonates
  21. Rapid Communication
  22. COVID-19 vaccine: Call for employees in international transportation industries and international travelers as the first priority in global distribution
  23. Case Report
  24. Rare squamous cell carcinoma of the kidney with concurrent xanthogranulomatous pyelonephritis: A case report and review of the literature
  25. An infertile female delivered a baby after removal of primary renal carcinoid tumor
  26. Research Articles
  27. Hypertension, BMI, and cardiovascular and cerebrovascular diseases
  28. Case Report
  29. Coexistence of bilateral macular edema and pale optic disc in the patient with Cohen syndrome
  30. Research Articles
  31. Correlation between kinematic sagittal parameters of the cervical lordosis or head posture and disc degeneration in patients with posterior neck pain
  32. Review Articles
  33. Hepatoid adenocarcinoma of the lung: An analysis of the Surveillance, Epidemiology, and End Results (SEER) database
  34. Research Articles
  35. Thermography in the diagnosis of carpal tunnel syndrome
  36. Pemetrexed-based first-line chemotherapy had particularly prominent objective response rate for advanced NSCLC: A network meta-analysis
  37. Comparison of single and double autologous stem cell transplantation in multiple myeloma patients
  38. The influence of smoking in minimally invasive spinal fusion surgery
  39. Impact of body mass index on left atrial dimension in HOCM patients
  40. Expression and clinical significance of CMTM1 in hepatocellular carcinoma
  41. miR-142-5p promotes cervical cancer progression by targeting LMX1A through Wnt/β-catenin pathway
  42. Comparison of multiple flatfoot indicators in 5–8-year-old children
  43. Early MRI imaging and follow-up study in cerebral amyloid angiopathy
  44. Intestinal fatty acid-binding protein as a biomarker for the diagnosis of strangulated intestinal obstruction: A meta-analysis
  45. miR-128-3p inhibits apoptosis and inflammation in LPS-induced sepsis by targeting TGFBR2
  46. Dynamic perfusion CT – A promising tool to diagnose pancreatic ductal adenocarcinoma
  47. Biomechanical evaluation of self-cinching stitch techniques in rotator cuff repair: The single-loop and double-loop knot stitches
  48. Review Articles
  49. The ambiguous role of mannose-binding lectin (MBL) in human immunity
  50. Case Report
  51. Membranous nephropathy with pulmonary cryptococcosis with improved 1-year follow-up results: A case report
  52. Fertility problems in males carrying an inversion of chromosome 10
  53. Acute myeloid leukemia with leukemic pleural effusion and high levels of pleural adenosine deaminase: A case report and review of literature
  54. Metastatic renal Ewing’s sarcoma in adult woman: Case report and review of the literature
  55. Burkitt-like lymphoma with 11q aberration in a patient with AIDS and a patient without AIDS: Two cases reports and literature review
  56. Skull hemophilia pseudotumor: A case report
  57. Judicious use of low-dosage corticosteroids for non-severe COVID-19: A case report
  58. Adult-onset citrullinaemia type II with liver cirrhosis: A rare cause of hyperammonaemia
  59. Clinicopathologic features of Good’s syndrome: Two cases and literature review
  60. Fatal immune-related hepatitis with intrahepatic cholestasis and pneumonia associated with camrelizumab: A case report and literature review
  61. Research Articles
  62. Effects of hydroxyethyl starch and gelatin on the risk of acute kidney injury following orthotopic liver transplantation: A multicenter retrospective comparative clinical study
  63. Significance of nucleic acid positive anal swab in COVID-19 patients
  64. circAPLP2 promotes colorectal cancer progression by upregulating HELLS by targeting miR-335-5p
  65. Ratios between circulating myeloid cells and lymphocytes are associated with mortality in severe COVID-19 patients
  66. Risk factors of left atrial appendage thrombus in patients with non-valvular atrial fibrillation
  67. Clinical features of hypertensive patients with COVID-19 compared with a normotensive group: Single-center experience in China
  68. Surgical myocardial revascularization outcomes in Kawasaki disease: systematic review and meta-analysis
  69. Decreased chromobox homologue 7 expression is associated with epithelial–mesenchymal transition and poor prognosis in cervical cancer
  70. FGF16 regulated by miR-520b enhances the cell proliferation of lung cancer
  71. Platelet-rich fibrin: Basics of biological actions and protocol modifications
  72. Accurate diagnosis of prostate cancer using logistic regression
  73. miR-377 inhibition enhances the survival of trophoblast cells via upregulation of FNDC5 in gestational diabetes mellitus
  74. Prognostic significance of TRIM28 expression in patients with breast carcinoma
  75. Integrative bioinformatics analysis of KPNA2 in six major human cancers
  76. Exosomal-mediated transfer of OIP5-AS1 enhanced cell chemoresistance to trastuzumab in breast cancer via up-regulating HMGB3 by sponging miR-381-3p
  77. A four-lncRNA signature for predicting prognosis of recurrence patients with gastric cancer
  78. Knockdown of circ_0003204 alleviates oxidative low-density lipoprotein-induced human umbilical vein endothelial cells injury: Circulating RNAs could explain atherosclerosis disease progression
  79. Propofol postpones colorectal cancer development through circ_0026344/miR-645/Akt/mTOR signal pathway
  80. Knockdown of lncRNA TapSAKI alleviates LPS-induced injury in HK-2 cells through the miR-205/IRF3 pathway
  81. COVID-19 severity in relation to sociodemographics and vitamin D use
  82. Clinical analysis of 11 cases of nocardiosis
  83. Cis-regulatory elements in conserved non-coding sequences of nuclear receptor genes indicate for crosstalk between endocrine systems
  84. Four long noncoding RNAs act as biomarkers in lung adenocarcinoma
  85. Real-world evidence of cytomegalovirus reactivation in non-Hodgkin lymphomas treated with bendamustine-containing regimens
  86. Relation between IL-8 level and obstructive sleep apnea syndrome
  87. circAGFG1 sponges miR-28-5p to promote non-small-cell lung cancer progression through modulating HIF-1α level
  88. Nomogram prediction model for renal anaemia in IgA nephropathy patients
  89. Effect of antibiotic use on the efficacy of nivolumab in the treatment of advanced/metastatic non-small cell lung cancer: A meta-analysis
  90. NDRG2 inhibition facilitates angiogenesis of hepatocellular carcinoma
  91. A nomogram for predicting metabolic steatohepatitis: The combination of NAMPT, RALGDS, GADD45B, FOSL2, RTP3, and RASD1
  92. Clinical and prognostic features of MMP-2 and VEGF in AEG patients
  93. The value of miR-510 in the prognosis and development of colon cancer
  94. Functional implications of PABPC1 in the development of ovarian cancer
  95. Prognostic value of preoperative inflammation-based predictors in patients with bladder carcinoma after radical cystectomy
  96. Sublingual immunotherapy increases Treg/Th17 ratio in allergic rhinitis
  97. Prediction of improvement after anterior cruciate ligament reconstruction
  98. Effluent Osteopontin levels reflect the peritoneal solute transport rate
  99. circ_0038467 promotes PM2.5-induced bronchial epithelial cell dysfunction
  100. Significance of miR-141 and miR-340 in cervical squamous cell carcinoma
  101. Association between hair cortisol concentration and metabolic syndrome
  102. Microvessel density as a prognostic indicator of prostate cancer: A systematic review and meta-analysis
  103. Characteristics of BCR–ABL gene variants in patients of chronic myeloid leukemia
  104. Knee alterations in rheumatoid arthritis: Comparison of US and MRI
  105. Long non-coding RNA TUG1 aggravates cerebral ischemia and reperfusion injury by sponging miR-493-3p/miR-410-3p
  106. lncRNA MALAT1 regulated ATAD2 to facilitate retinoblastoma progression via miR-655-3p
  107. Development and validation of a nomogram for predicting severity in patients with hemorrhagic fever with renal syndrome: A retrospective study
  108. Analysis of COVID-19 outbreak origin in China in 2019 using differentiation method for unusual epidemiological events
  109. Laparoscopic versus open major liver resection for hepatocellular carcinoma: A case-matched analysis of short- and long-term outcomes
  110. Travelers’ vaccines and their adverse events in Nara, Japan
  111. Association between Tfh and PGA in children with Henoch–Schönlein purpura
  112. Can exchange transfusion be replaced by double-LED phototherapy?
  113. circ_0005962 functions as an oncogene to aggravate NSCLC progression
  114. Circular RNA VANGL1 knockdown suppressed viability, promoted apoptosis, and increased doxorubicin sensitivity through targeting miR-145-5p to regulate SOX4 in bladder cancer cells
  115. Serum intact fibroblast growth factor 23 in healthy paediatric population
  116. Algorithm of rational approach to reconstruction in Fournier’s disease
  117. A meta-analysis of exosome in the treatment of spinal cord injury
  118. Src-1 and SP2 promote the proliferation and epithelial–mesenchymal transition of nasopharyngeal carcinoma
  119. Dexmedetomidine may decrease the bupivacaine toxicity to heart
  120. Hypoxia stimulates the migration and invasion of osteosarcoma via up-regulating the NUSAP1 expression
  121. Long noncoding RNA XIST knockdown relieves the injury of microglia cells after spinal cord injury by sponging miR-219-5p
  122. External fixation via the anterior inferior iliac spine for proximal femoral fractures in young patients
  123. miR-128-3p reduced acute lung injury induced by sepsis via targeting PEL12
  124. HAGLR promotes neuron differentiation through the miR-130a-3p-MeCP2 axis
  125. Phosphoglycerate mutase 2 is elevated in serum of patients with heart failure and correlates with the disease severity and patient’s prognosis
  126. Cell population data in identifying active tuberculosis and community-acquired pneumonia
  127. Prognostic value of microRNA-4521 in non-small cell lung cancer and its regulatory effect on tumor progression
  128. Mean platelet volume and red blood cell distribution width is associated with prognosis in premature neonates with sepsis
  129. 3D-printed porous scaffold promotes osteogenic differentiation of hADMSCs
  130. Association of gene polymorphisms with women urinary incontinence
  131. Influence of COVID-19 pandemic on stress levels of urologic patients
  132. miR-496 inhibits proliferation via LYN and AKT pathway in gastric cancer
  133. miR-519d downregulates LEP expression to inhibit preeclampsia development
  134. Comparison of single- and triple-port VATS for lung cancer: A meta-analysis
  135. Fluorescent light energy modulates healing in skin grafted mouse model
  136. Silencing CDK6-AS1 inhibits LPS-induced inflammatory damage in HK-2 cells
  137. Predictive effect of DCE-MRI and DWI in brain metastases from NSCLC
  138. Severe postoperative hyperbilirubinemia in congenital heart disease
  139. Baicalin improves podocyte injury in rats with diabetic nephropathy by inhibiting PI3K/Akt/mTOR signaling pathway
  140. Clinical factors predicting ureteral stent failure in patients with external ureteral compression
  141. Novel H2S donor proglumide-ADT-OH protects HUVECs from ox-LDL-induced injury through NF-κB and JAK/SATA pathway
  142. Triple-Endobutton and clavicular hook: A propensity score matching analysis
  143. Long noncoding RNA MIAT inhibits the progression of diabetic nephropathy and the activation of NF-κB pathway in high glucose-treated renal tubular epithelial cells by the miR-182-5p/GPRC5A axis
  144. Serum exosomal miR-122-5p, GAS, and PGR in the non-invasive diagnosis of CAG
  145. miR-513b-5p inhibits the proliferation and promotes apoptosis of retinoblastoma cells by targeting TRIB1
  146. Fer exacerbates renal fibrosis and can be targeted by miR-29c-3p
  147. The diagnostic and prognostic value of miR-92a in gastric cancer: A systematic review and meta-analysis
  148. Prognostic value of α2δ1 in hypopharyngeal carcinoma: A retrospective study
  149. No significant benefit of moderate-dose vitamin C on severe COVID-19 cases
  150. circ_0000467 promotes the proliferation, metastasis, and angiogenesis in colorectal cancer cells through regulating KLF12 expression by sponging miR-4766-5p
  151. Downregulation of RAB7 and Caveolin-1 increases MMP-2 activity in renal tubular epithelial cells under hypoxic conditions
  152. Educational program for orthopedic surgeons’ influences for osteoporosis
  153. Expression and function analysis of CRABP2 and FABP5, and their ratio in esophageal squamous cell carcinoma
  154. GJA1 promotes hepatocellular carcinoma progression by mediating TGF-β-induced activation and the epithelial–mesenchymal transition of hepatic stellate cells
  155. lncRNA-ZFAS1 promotes the progression of endometrial carcinoma by targeting miR-34b to regulate VEGFA expression
  156. Anticoagulation is the answer in treating noncritical COVID-19 patients
  157. Effect of late-onset hemorrhagic cystitis on PFS after haplo-PBSCT
  158. Comparison of Dako HercepTest and Ventana PATHWAY anti-HER2 (4B5) tests and their correlation with silver in situ hybridization in lung adenocarcinoma
  159. VSTM1 regulates monocyte/macrophage function via the NF-κB signaling pathway
  160. Comparison of vaginal birth outcomes in midwifery-led versus physician-led setting: A propensity score-matched analysis
  161. Treatment of osteoporosis with teriparatide: The Slovenian experience
  162. New targets of morphine postconditioning protection of the myocardium in ischemia/reperfusion injury: Involvement of HSP90/Akt and C5a/NF-κB
  163. Superenhancer–transcription factor regulatory network in malignant tumors
  164. β-Cell function is associated with osteosarcopenia in middle-aged and older nonobese patients with type 2 diabetes: A cross-sectional study
  165. Clinical features of atypical tuberculosis mimicking bacterial pneumonia
  166. Proteoglycan-depleted regions of annular injury promote nerve ingrowth in a rabbit disc degeneration model
  167. Effect of electromagnetic field on abortion: A systematic review and meta-analysis
  168. miR-150-5p affects AS plaque with ASMC proliferation and migration by STAT1
  169. MALAT1 promotes malignant pleural mesothelioma by sponging miR-141-3p
  170. Effects of remifentanil and propofol on distant organ lung injury in an ischemia–reperfusion model
  171. miR-654-5p promotes gastric cancer progression via the GPRIN1/NF-κB pathway
  172. Identification of LIG1 and LIG3 as prognostic biomarkers in breast cancer
  173. MitoQ inhibits hepatic stellate cell activation and liver fibrosis by enhancing PINK1/parkin-mediated mitophagy
  174. Dissecting role of founder mutation p.V727M in GNE in Indian HIBM cohort
  175. circATP2A2 promotes osteosarcoma progression by upregulating MYH9
  176. Prognostic role of oxytocin receptor in colon adenocarcinoma
  177. Review Articles
  178. The function of non-coding RNAs in idiopathic pulmonary fibrosis
  179. Efficacy and safety of therapeutic plasma exchange in stiff person syndrome
  180. Role of cesarean section in the development of neonatal gut microbiota: A systematic review
  181. Small cell lung cancer transformation during antitumor therapies: A systematic review
  182. Research progress of gut microbiota and frailty syndrome
  183. Recommendations for outpatient activity in COVID-19 pandemic
  184. Rapid Communication
  185. Disparity in clinical characteristics between 2019 novel coronavirus pneumonia and leptospirosis
  186. Use of microspheres in embolization for unruptured renal angiomyolipomas
  187. COVID-19 cases with delayed absorption of lung lesion
  188. A triple combination of treatments on moderate COVID-19
  189. Social networks and eating disorders during the Covid-19 pandemic
  190. Letter
  191. COVID-19, WHO guidelines, pedagogy, and respite
  192. Inflammatory factors in alveolar lavage fluid from severe COVID-19 pneumonia: PCT and IL-6 in epithelial lining fluid
  193. COVID-19: Lessons from Norway tragedy must be considered in vaccine rollout planning in least developed/developing countries
  194. What is the role of plasma cell in the lamina propria of terminal ileum in Good’s syndrome patient?
  195. Case Report
  196. Rivaroxaban triggered multifocal intratumoral hemorrhage of the cabozantinib-treated diffuse brain metastases: A case report and review of literature
  197. CTU findings of duplex kidney in kidney: A rare duplicated renal malformation
  198. Synchronous primary malignancy of colon cancer and mantle cell lymphoma: A case report
  199. Sonazoid-enhanced ultrasonography and pathologic characters of CD68 positive cell in primary hepatic perivascular epithelioid cell tumors: A case report and literature review
  200. Persistent SARS-CoV-2-positive over 4 months in a COVID-19 patient with CHB
  201. Pulmonary parenchymal involvement caused by Tropheryma whipplei
  202. Mediastinal mixed germ cell tumor: A case report and literature review
  203. Ovarian female adnexal tumor of probable Wolffian origin – Case report
  204. Rare paratesticular aggressive angiomyxoma mimicking an epididymal tumor in an 82-year-old man: Case report
  205. Perimenopausal giant hydatidiform mole complicated with preeclampsia and hyperthyroidism: A case report and literature review
  206. Primary orbital ganglioneuroblastoma: A case report
  207. Primary aortic intimal sarcoma masquerading as intramural hematoma
  208. Sustained false-positive results for hepatitis A virus immunoglobulin M: A case report and literature review
  209. Peritoneal loose body presenting as a hepatic mass: A case report and review of the literature
  210. Chondroblastoma of mandibular condyle: Case report and literature review
  211. Trauma-induced complete pacemaker lead fracture 8 months prior to hospitalization: A case report
  212. Primary intradural extramedullary extraosseous Ewing’s sarcoma/peripheral primitive neuroectodermal tumor (PIEES/PNET) of the thoracolumbar spine: A case report and literature review
  213. Computer-assisted preoperative planning of reduction of and osteosynthesis of scapular fracture: A case report
  214. High quality of 58-month life in lung cancer patient with brain metastases sequentially treated with gefitinib and osimertinib
  215. Rapid response of locally advanced oral squamous cell carcinoma to apatinib: A case report
  216. Retrieval of intrarenal coiled and ruptured guidewire by retrograde intrarenal surgery: A case report and literature review
  217. Usage of intermingled skin allografts and autografts in a senior patient with major burn injury
  218. Retraction
  219. Retraction on “Dihydromyricetin attenuates inflammation through TLR4/NF-kappa B pathway”
  220. Special Issue Computational Intelligence Methodologies Meets Recurrent Cancers - Part I
  221. An artificial immune system with bootstrap sampling for the diagnosis of recurrent endometrial cancers
  222. Breast cancer recurrence prediction with ensemble methods and cost-sensitive learning
Downloaded on 8.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/med-2021-0241/html
Scroll to top button