Genomics Score Based on Genome-Wide Network Analysis for Prediction of Survival in Gastric Cancer: A Novel Prognostic Signature

Purpose: Gastric cancer (GC) is a product of multiple genetic abnormalities, including genetic and epigenetic modifications. This study aimed to integrate various biomolecules, such as miRNAs, mRNA, and DNA methylation, into a genome-wide network and develop a nomogram for predicting the overall survival (OS) of GC.

Materials and Methods: A total of 329 GC cases, as a training cohort with a random of 150 examples included as a validation cohort, were screened from The Cancer Genome Atlas database. A genome-wide network was constructed based on a combination of univariate Cox regression and least absolute shrinkage and selection operator analyses, and a nomogram was established to predict 1-, 3-, and 5-year OS in the training cohort. The nomogram was then assessed in terms of calibration, discrimination, and clinical usefulness in the validation cohort. Afterward, in order to confirm the superiority of the whole gene network model and further reduce the biomarkers for the improvement of clinical usefulness, we also constructed eight other models according to the different combinations of miRNAs, mRNA, and DNA methylation sites and made corresponding comparisons. Finally, Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were also performed to describe the function of this genome-wide network.

Results: A multivariate analysis revealed a novel prognostic factor, a genomics score (GS) comprising seven miRNAs, eight mRNA, and 19 DNA methylation sites. In the validation cohort, comparing to patients with low GS, high-GS patients (HR, 12.886; P < 0.001) were significantly associated with increased all-cause mortality. Furthermore, after stratification of the TNM stage (I, II, III, and IV), there were significant differences revealed in the survival rates between the high-GS and low-GS groups as well (P < 0.001). The 1-, 3-, and 5-year C-index of whole genomics-based nomogram were 0.868, 0.895, and 0.928, respectively. The other models have comparable or relatively poor comprehensive performance, while they had fewer biomarkers. Besides that, DAVID 6.8 further revealed multiple molecules and pathways related to the genome-wide network, such as cytomembranes, cell cycle, and adipocytokine signaling.

Conclusion: We successfully developed a GS based on genome-wide network, which may represent a novel prognostic factor for GC. A combination of GS and TNM staging provides additional precision in stratifying patients with different OS prognoses, constituting a more comprehensive sub-typing system. This could potentially play an important role in future clinical practice.

Background

Gastric cancer (GC) is one of the most common malignant human tumors and the third leading cause of cancer-related mortalities worldwide. Reports estimate that nearly one million new cases and 800,000 deaths occur each year across the world (Torre et al., 2015). Despite the rapid research advancement, GC-related impacts on human life remain high around the globe. According to the global cancer burden data, hundreds of billions of dollars in economic losses are incurred each year due to GC. At the same time, stomach cancer has been reported to cause 19.1 million disability-adjusted life years, with 98% of these resulting from years of life lost and 2% from years lived with disability (Global Burden of Disease Cancer Collaboration et al., 2019).

Despite major breakthroughs in GC prevention, diagnosis, and treatment therapies reported over the past decade, prognosis remains a challenge at different TNM stages (Jiang et al., 2018a; Sun et al., 2019a,b). Notably, patients with similar clinical features and at the same tumor stage who receive uniform treatment have exhibited varying clinical outcomes (Bang et al., 2012; Jiang et al., 2018b). Such evidence indicates the existing challenges to traditional TNM staging (Serra et al., 2019), possibly due to a lack of molecular tools for effectively predicting the prognosis and the therapeutic effect of GC patients. Therefore, more rigorous and reliable systems that accurately reflect the heterogeneity of different patients and guide the development of treatment approaches are urgently needed (Duarte et al., 2018; Serra et al., 2019).

Tumors are a product of multiple genetic mutations, including genetic (gene expression) and epigenetic (DNA methylation and histone modification) modifications, as well as deregulations of tumor-suppressor genes and proto-oncogenes (Anna et al., 2018; Choi et al., 2019). In addition, changes in a set of genetic materials have been closely associated with cancer outcomes (Anna et al., 2018; Choi et al., 2019). Therefore, to effectively predict the prognosis of tumors, such as GC, a single biomarker is insufficient, necessitating the need for a gene network.

A variety of mRNAs have been associated with GC prognosis (Camargo et al., 2019), with microRNAs (miRNAs) also implicated in tumor prediction in the recent years (Li et al., 2010; Ueda et al., 2010; Camargo et al., 2019). These small, non-coding RNAs, comprising 22 nucleotides, primarily function to regulate protein translation by inhibiting the expression of target messenger RNAs (mRNAs). Apart from genetics, epigenetics is currently receiving considerable research attention. DNA methylation is the most common epigenetic event associated with cancer development and progression (Anna et al., 2018). Consequently, numerous studies have implicated DNA methylation in the diagnosis and the prognosis of various tumors, including GC (Camargo et al., 2019; Choi et al., 2019). Although these studies have revealed several biomarkers that have proved to be prognostic predictors in GC, only a handful have been adopted in clinical therapies or are used to build predictive models for the disease (Anna et al., 2018; Duarte et al., 2018; Camargo et al., 2019; Choi et al., 2019; Serra et al., 2019).

Previous studies have identified and recommended numerous biomarkers for GC. However, since malignant tumors often involve multiple layers and different levels of genetic changes, including the genome, transcriptome, and proteome, or even epigenetic content, selecting reasonable candidate factors from tens of thousands of biomarkers and comprehensively analyzing them as an independent feature is imperative to effectively develop a suitable prognostic target. Therefore, genetic networks containing a panel of abnormal factors from different regulatory levels represent the best chance for achieving prognostic value.

The whole genome-wide network analysis is reported in several other cancers, such as colorectal cancer, breast cancer, and lung cancer (Hou et al., 2018; Zhang et al., 2018), and it shows great value in differentiating the prognosis of these patients. Therefore, it is feasible and advantageous to apply genome-wide network analysis to GC.

In the current study, we performed a series of sophisticated statistical analyses and identified 33 genetic molecules that were highly correlated with the prognosis of GC. Specifically, we screened The Cancer Genome Atlas (TCGA), a genome project with 33 types of cancer, including gene expression, and DNA methylation as well as other biological information. Furthermore, we extended these independent prognostic factors to the “omics” concept. Then, a genome-wide network was constructed. Interestingly, the genomics score (GS) obtained herein could supplement TNM staging and enhance the prognostic value of different patients. Moreover, we developed multiple prognostic models, then validated, and compared them to ascertain their strengths and weaknesses. Finally, we performed pathway enrichment and gene oncology annotation analyses to elucidate the function of this gene network.

Materials and Methods

Data Acquisition and Preprocessing

Level 3 data were downloaded from the TCGA database using TCGA-Assembler Module A, in January 2019, which was then pretreated with Module B. The dataset comprised of clinical variables from 443 patients, including age, sex, stage, primary site, grade, treatment, and survival, as well as associated genome-wide data. In addition, the expression levels of 1,871 miRNAs, 20,531 mRNA, and 485,577 DNA methylation sites (Illumina methylation 450) were obtained from 384, 377, and 394 patients, respectively. Afterward, an intersection with a total of 332 samples was eventually retained. Furthermore, patients with missing active follow-up data were excluded from the analysis, leaving 329 patients in the final cohort (Figure 1). Moreover, genome-wide level 3 data whose expression levels for miRNAs, mRNA, and DNA methylation sites were missing in more than 50% of all samples were excluded from the final analysis. Finally, 329 GC patients with 566 miRNAs, 17,963 mRNA, and 396,081 DNA methylation sites were chosen for further analysis.

FIGURE 1

Figure 1. A Venn diagram displays the patients’ screening process.

Genome-Wide Network Analysis

Gene expression and DNA methylation data were normalized using R package before subsequent processing. Univariate and least absolute shrinkage and selection operator (LASSO) Cox regression models were combined and used to identify the most useful prognostic factors in miRNAs, mRNA, and DNA methylation sites associated with survival. Firstly, univariate Cox regression was performed on each candidate miRNA, mRNA, and DNA methylation site to elucidate its role in patient survival, then signatures with P value less than 0.05 were retained for subsequent analysis. Thereafter, the LASSO Cox regression model was applied to select and shrink the data (Supplementary Figure S3; Tibshirani, 1997). Finally, a GS, based on a genome-wide network comprising seven miRNAs, eight mRNAs, and 19 DNA methylation sites, was constructed for predicting survival. A summary of the whole screening process is displayed in Supplementary Figure S1.

Development and Comparison of Individualized Prediction Models

The TCGA cohort with 329 cases was used as the training set, with a random 150 cases from the total cohort included as a validation group. The random number is 1,356. Firstly, we developed the original GS based on 34 biomarkers (seven miRNAs, eight mRNAs, and 19 DNA methylation sites). Then, considering the complexity of the original GS and difficult clinical application, in order to obtain a more concise and effective GS, we also constructed eight other models according to the different combinations of miRNAs, mRNA, and DNA methylation and made corresponding comparisons. Finally, a total of nine GS models based on the genome-wide network from LASSO were adopted to screen for the most appropriate markers. These included the following models: genomics (seven miRNAs, eight mRNA, and 19 DNA methylation sites), miRNAs (seven miRNAs), mRNA (eight mRNA), methylation (19 DNA methylation sites), miRNAs + methylation (seven miRNAs and 19 DNA methylation sites), miRNAs + mRNA (seven miRNAs and eight mRNA), mRNA + methylation (eight mRNA and 19 DNA methylation sites), Cox-model 1 (two miRNAs, six mRNA, and nine DNA methylation sites), and Cox-model 2 (one miRNA, one mRNA, and seven DNA methylation sites). Among them, markers from Cox-model 1 were separately detected from miRNAs, mRNA, or DNA methylation sites using multivariate Cox regression analysis after LASSO (Supplementary Tables S2–S4). On the other hand, markers from Cox-model 2 depended on signatures from a multivariate Cox regression analysis combining the genome-wide network and the clinical characteristics (Supplementary Table S5). Thereafter, we constructed several nomograms by incorporating significant (P < 0.05) GS variables and other clinical features following multivariate Cox regression (Iasonos et al., 2008), and a clinical nomogram was built as a blank control. The equations used for calculating the GS of these models are listed in Supplementary Table S6.

To calculate the discrimination and the stability of different Cox regression models, we applied C-statistics and calibration. Additionally, we performed an analysis of time-dependent receiver operator characteristics (ROC), based on the 1-, 3-, and 5-year survival endpoints, to assess the prognostic accuracy of the different nomograms. Furthermore, we evaluated the potential net benefit of different predictive models using decision curve analysis (DCA). DCA compares the clinical usefulness of different indicators by calculating the potential net benefit of each decision strategy at each threshold probability. Thus, DCA was a significant novel approach for comparing the old and the new models (Vickers and Elkin, 2006).

Screening for Potential miRNA Target Genes

We predicted the potential target genes of the seven miRNAs, from LASSO, by screening the miRTarBase, miRDB, and TargetScan databases. Common genes from each miRNA across the three databases were then used for subsequent studies. More than 90% of the miRNAs showed negative regulation to target genes. Consequently, the expression data from TCGA were used to perform a batch of correlation analysis of each miRNA, with corresponding target genes, and the three genes with the largest absolute negative correlation were retained as the most likely targets. Additionally, at least three potential target genes from miRTarBase, which is co-expressed with miRNAs, were considered as equally important and were subjected to Cytoscape (version 3.7.2) for identification of miRNA–target genes co-expression network analysis (Supplementary Figure S2).

Functional Enrichment Analysis

The potential target genes that were negatively correlated with miRNAs in TCGA, as well as the coding sequences for mRNA and DNA methylation sites, were used for functional enrichment analysis using the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway and Gene Ontology (GO) using DAVID 6.8 (Supplementary Figure S2). Functional enrichment analysis indicates why the gene network produces images on the survival of GC from a molecular mechanism. Visualization was then done using the “ggplot2” package implemented in R.

Statistical Analysis

The patients were divided into low-risk and high-risk groups by the median GS as the cutoff point. Survival estimates were obtained according to the Kaplan–Meier method and compared using the log-rank test. Variables that reached significance, with P < 0.05, were entered into the multivariable analyses using the Cox proportional hazards model, with an entry stepwise approach to identify covariates associated with increased all-cause mortality, and then hazard ratio with 95% confidence intervals (CIs) of each variable was achieved. All the statistical significance values were set as two-sided (P < 0.05). LASSO Cox regression was performed through the “glmnet” package. Time-dependent ROC analysis at different follow-up times was implemented using the “timeROC” package of R project in order to further expound the performance of different GS models, and DCA was used to compare their clinical use by “rmda” package. Finally, nomogram based on the Cox regression model was constructed using the “rms” package. C-index and calibration to calculate the discrimination and the stability of these models were performed using c-statistics and Bootstrap sample. Harrell’s concordance index (C-index) indicated a better prognostic model if its value was closer to 1, and the calibration diagram showed that the better the prediction if the closer the correction line was to the diagonal. All statistical methods are applied to both the training group and the validation group. Statistical analyses were performed using SPSS statistical software (version 18.0) and R software (version 3.5.3).

Results

Patient Characteristics

Among the 329 GC patients analyzed in this study, 212 (64.4%) were male, whereas 117 (35.6%) were female. The average age of the entire study population was 65.0 ± 10.6 years. In terms of pathological stage, 38 (11.6%) cases were identified as stage I, 117 (35.6%) were stage II, 155 (47.1%) were stage III, and 19 (5.8%) were at stage IV. With regards to treatment, 303 (92.1%) patients received surgery (280 cases of R0 surgery, 14 R1, and nine R2), whereas 146 (44.4%) were subjected to fluorouracil-based chemotherapy. The genomic nomogram classified 165 samples into low GS (GS ≤ -0.137) and 164 into high GS (GS > -0.137) groups based on the median cutoff (Figure 2). A detailed description of tumor location, pathology grade, and Lauren classification is outlined in Supplementary Table S1, while a heat map of the genomic scores layered by clinicopathological

FIGURE 2

Figure 2. Distribution of patient cases and density based on genomics score (GS) in the total The Cancer Genome Atlas cohort (A,B). Scatter plots of genomics scores regarding the classification of low and high GS (C,D), and the bold line represents the median.

factors is illustrated in Supplementary Figure S4. The median (mean; 95% CI) survival time for OS was 1,043 (670.2–1,415.8) days in the total cohort, 466 (370.6–561.4) days in the high-GS group, and 2,613 (mean; 2209.4–2017.5) days in the low-GS group (Supplementary Figure S5). Toward the last follow-up, a total of 129 deaths and 200 censoring were recorded. The estimated cumulative 1-, 3-, and 5-year OS in the total cohort were 78.9, 48.4, and 36.7%, respectively, although these rates increase to 95.1, 74.2, and 68.5%, respectively, in the low-GS group. Conversely, the 1-, 3-, and 5-year OS decreased to 63.3, 23.6, and 15.4%, respectively, in the high-GS group. The baseline information of the validation cohort is also listed in Supplementary Table S1 and Supplementary Figure S6.

Survival Analysis

We identified a basic genome-wide network comprising seven miRNAs, eight mRNAs, and 19 DNA methylation sites as the prognostic factor for OS, from hundreds of thousands of univariate Cox regression and LASSO analyses. This network was then classified as other models in the training and the validation groups. Among the 34 features identified, poor prognosis was significantly associated with a high expression of seven miRNAs (hsa-mir-100, hsa-mir-1304, hsa-mir-136, hsa-mir-193b, hsa-mir-22, hsa-mir-653, and hsa-mir-6808), six mRNAs (NRP1|8829, RNF144A|9781, ZNF22|7570, DUSP1|1843, CPNE8|144402, MAGED1|9500, and LOC91450|91450), and seven DNA methylation sites (cg07020967, cg08859156, cg12485556, cg15861578, cg15861578, cg25161386, and cg22740006). Conversely, poor prognosis was strongly associated with a low expression of SOX14|8403 and 12 DNA methylation sites, including cg02223323, cg00481239, cg14791193, cg15486740, cg20100408, cg20350671, cg22395807, cg24361571, cg25361506, cg22813794, cg26014401, and cg26856948 (Table 2). Univariate analysis performed on clinical characteristics revealed a significant association between age, pathological stage, TNM, and surgery with OS (Table 1). On the other hand, results from multivariable Cox regression showed that age, pathological stage, and GS were significantly associated with all-cause mortality in GC (Table 1 and Figure 3). Furthermore, stratification of the pathological stage (I, II, III, and IV) revealed significant differences in survival rates between the high-GS and the low-GS groups (Figure 3). A similar result was found when the data were stratified by demographic variables (sex and age), clinical characteristics (primary site, grade, and Lauren classification) as well as treatments (surgery and chemotherapy; Supplementary Figures S7, S8). On the other hand, categorizing GS into high or low groups, using the median value across different models, indicated that the genomics nomogram had the highest HR value. Interestingly, HR was almost equal to miRNAs + methylation, mRNA + methylation, Cox-model 1, and Cox-model 2 nomograms, which contained fewer gene features. Moreover, the HR value showed a marked decrease in miRNAs, mRNA, methylation, and miRNAs + mRNA nomograms, which included the least characteristics (Table 3).

TABLE 1

Table 1. Univariable and multivariable analyses of the genomics score and the clinicopathological characteristics for overall survival in the training group and the validation group.

FIGURE 3

Figure 3. Kaplan–Meier curve of overall survival in all patients, then stratified by genomics score (GS), pathological stage, and age. Survival analysis in the low- and high-GS groups was further divided based on stages (stages I–IV).

TABLE 2

Table 2. miRNAs, mRNA, and DNA methylation whose expression levels showed a significant association with overall survival in least absolute shrinkage and selection operator.

TABLE 3

Table 3. Comparison of different genomics score models (based on the median value) for overall survival in the training group and the validation group.

Nomograms Based on Genome-Wide Network

A genomics nomogram was first constructed based on the genome-wide network, comprising 34 gene features (Figure 4). To obtain a more concise and effective nomogram, we also built a Cox-model 1 (17 gene features) and Cox-model 2 (nine gene features) nomograms (Supplementary Figures S9, S10). Next, a clinical nomogram, based on stage and age, was built as a control (Supplementary Figure S11). Thereafter, we performed internal and external validation to evaluate the feasibility of all nomograms using a three-grouped random bootstrap sampling (Figure 5 and Supplementary Figures S9–S11). We observed good predictive performance in the first three nomograms, but not in the simple clinical model.

FIGURE 4

Figure 4. Genomics nomogram to predict the probability of 1-, 3-, and 5-year overall survival (OS) in the training cohort (A) and the validation cohort (B): to determine how many points for each variable to the probability of OS, locate the variable on its axis, draw a line straight upward to the point axis, repeat this process for each variable, sum up the points achieved for each of the risk factors, locate the final sum on the total point axis, and draw a line straight down to find the patient’s probability of OS.

FIGURE 5

Figure 5. Time-dependent receiver operating characteristic curves on 1, 3, and 5 years of genomics nomogram. Calibration plot showing 1-, 3-, and 5-year overall survival (OS) probability: the nomogram-estimated OS is plotted on the x-axis, and the actual OS is plotted on the y-axis. The diagonal dotted line is a perfect estimation by an ideal model, in which the estimated outcome perfectly corresponds to the actual outcome. The solid line is the performance of the nomogram: A closer alignment with the diagonal dotted line represents a better estimation. Decision curve analysis for genomics nomogram and clinical nomogram: the y-axis measures the net benefit. The gray line or the horizontal black line represents a follow-up of all patients or no patients. The model makes more benefit with a higher position in curve.

Validation of the Nomograms Using ROC and DCA

To ensure a good comparison across different GS nomograms, we performed a time-dependent ROC (at 1, 3, and 5 years of follow-up) as well as DCA. In the validation group, genomics nomogram revealed the best comprehensive performance, with 1-, 3-, and 5-year area under the curve (AUC) values of 0.868, 0.895, and 0.928, respectively (Table 4), and Cox-model 1, miRNAs + methylation, and mRNA + methylation nomograms had a comparable performance, with 1-, 3-, and 5-year AUC values of 0.856–0.873, 0.884–0.905, and 0.907–0.919, respectively, but it had fewer biomarkers (Table 4). Although the Cox-model 2 nomograms had the least biomarkers, including miRNA, mRNA, and DNA methylation sites, it had a relatively poor performance with 1-, 3-, and 5-year AUC values of 0.835, 0.859, and 0.785, respectively. Besides that, the miRNA, mRNA, methylation, miRNAs + methylation, and miRNAs + mRNA nomograms recorded 1-, 3-, and 5-year AUC values of 0.729–0.877, 0.656–0.805, and 0.721–0.894, respectively. Finally, we found that, compared to miRNA (0.641, 0.729, and 0.736) and mRNA nomogram (0.806, 0.785, and 0.843), methylation nomogram had higher 1-, 3-, and 5-year AUC values of 0.866, 0.877, and 0.894. Nevertheless, all of them showed better performance than the clinical nomogram, which recorded 1-, 3-, and 5-year AUC values of 0.638, 0.598, and 0.721, respectively (Figures 6A,B and Supplementary Figure S12). The C-index based on different nomograms exhibited a similar effect (Supplementary Table S8). Additionally, DCA showed that the genomics, Cox-model 1, mRNA + methylation, and methylation nomograms had a significant net benefit compared to other GS models and the clinical nomogram (Figures 6C,D).

TABLE 4

Table 4. The area under the curve (AUC) values of different genomics score models in the training group and the validation group.

FIGURE 6

Figure 6. Time-dependent receiver operating characteristic curves on 1, 3, and 5 years of each nomogram and decision curve analysis for each nomogram.

Potential miRNA Target Genes

A total of 72 hsa-mir-22, 39 hsa-mir-100, 56 hsa-mir-136, 58 hsa-mir-193b, 23 hsa-mir-653, 96 hsa-mir-1304, and 285 hsa-mir-6808 potential target genes were identified from the miRTarBase, miRDB, and TargetScan databases (Supplementary Figure S14). We then performed a correlation analysis between each target gene and miRNAs and finally generated a miRNA–potential target gene plot (Supplementary Figure S15A) as well as a miRNA–target gene co-expression network (Supplementary Figure S15B) using Cytoscape.

Functional Analysis of Genome-Wide Network

We imported the 301 potential target genes, mRNA, and DNA methylation site-coding sequences, identified above, into DAVID for KEGG and GO analyses and identified biological processes, molecular functions as well as cellular components (Figures 7A–C). Their corresponding KEGG pathways were also plotted (Figure 7D).

FIGURE 7

Figure 7. Gene enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways for the genome-wide network. The biological process (A), molecular function (B), cellular component (C), and KEGG pathways (D).

Discussion

GC can be divided into two types or four main categories, according to the Lauren and World Health Organization (WHO) classifications (Lauren, 1965; Nagtegaal et al., 2019), although neither of these classifications is based on molecular markers. In the last decade, however, three novel molecular-based classification systems have been suggested for GC. The Singapore-Duke Group was the first to describe a classification with two intrinsic genomic subtypes, G-INT, and G-DIF, which had different gene expression (Tan et al., 2011; Serra et al., 2019). The subtypes have different levels of resistance to various chemotherapy drugs and show limited prognostic value. Later, TCGA used molecular evaluation to propose a new classification with four subtypes: EBV, MSI, GS, and CIN. The identification of these subtypes has provided a roadmap for patient stratification as well as targeted therapeutic trials (Cancer Genome Atlas Research, 2014). However, initial data on disease outcomes from this cohort did not show differences in survival among the four groups. A series of positive studies on prognosis based on TCGA classification was also reported (Sohn et al., 2017). In addition, the Asian Cancer Research group divided GC into four subtypes, MSI, EMT, MSS/TP53+, and MSS/TP53-, based on gene expression data and found significantly different survival outcomes across them (Cristescu et al., 2015; Serra et al., 2019). Despite the significant milestones of these studies, they are all mainly based on the analysis of gene expression (mRNA). Besides that, a 2019 study proposed a five-miRNA model, while it had a C-index of 0.72 only (Zhang et al., 2019). In the current study, we included methylation data and performed functional enrichment analysis, making our work stronger. The aforementioned classifications are also complicated and need further optimization to increase clinical applicability. Furthermore, they focused on typing and finding new targets, whereas our study reports on prognostic analysis.

Some of the biomarkers we identified herein, including hsa-mir-22, hsa-mir-100, hsa-mir-136, hsa-mir-193b, hsa-mir-1304, NRP1, DUSP1, and MAP7D2 (cg02223323), have previously been reported in GC (Grandclement and Borg, 2011; Chen et al., 2014; Mu et al., 2014; Zuo et al., 2015; Zheng et al., 2017; Chen et al., 2018; Kurata and Lin, 2018; Liu K.T. et al., 2018; Song et al., 2018; Teng et al., 2018; Liu et al., 2019; Pan et al., 2019; Wang et al., 2019). Others, such as CPNE8, MAGED1, RNF144A, SOX14, ACOT13 (cg15486740), EID1 (cg00481239), RPS4X (cg08859156), and TTRAP (cg15486740), have been identified in various tumors other than GC (Kamio et al., 2010; Zeng et al., 2012; Zhou et al., 2013; Kuang et al., 2017; Stanisavljevic et al., 2017; Liu X. et al., 2018; Tosic et al., 2018; Nagasawa et al., 2019; Yang et al., 2019). The remaining biomarkers, including hsa-mir-653, hsa-mir-6808, LOC91450, ZNF22, C1orf144 (cg14791193), GOLGA3 (cg26856948), HLA-DPB1 (cg20100408), LRFN4 (cg22740006), MDH2 (cg22813794), MIR365-2 (cg24361571), NUFIP2 (cg25161386), PREP (cg12485556), STYXL1 (cg22813794), TMEM117 (cg07020967), ZC3H10 (ZC3H10), IL1RAPL1 (cg20350671), PC (cg22740006), SHC4 (cg00481239), and ATXN10 (cg22395807), have not been previously reported.

Currently, focus has been directed on identifying prognostic miRNAs for GC. Particularly, one miRNA can regulate multiple targets, while multiple miRNAs can regulate a single mRNA. Therefore, a single miRNA may play an opposite role in cancer progression by regulating different target genes. For example, Mir-22 and Mir-100 were found to be tumor suppressors in various cancers, including GC (Chen et al., 2014; Zuo et al., 2015). Similarly, a high expression of Mir-136 was found to promote proliferation and invasion in GC cell lines by inhibiting PTEN expression (Chen et al., 2018), while a contrasting result was reported when HOXC10 was targeted (Zheng et al., 2017). Similarly, Mir-193b reportedly induced GC proliferation or apoptosis by mediating different mRNA expressions (Mu et al., 2014; Song et al., 2018), whereas a high Mir-1304 expression in GC was reported as a negative predictor for prognosis of lung and thyroid cancers (Kurata and Lin, 2018; Liu et al., 2019; Pan et al., 2019). However, the function of Mir-653 and Mir-6808 has not been previously reported. In the current study, we found an association between a high expression of all miRNAs and poor survival. Different outcomes may be observed in our study, relative to previous reports, owing to the huge number of corresponding miRNA target genes herein and lack of evidence on their role in GC development.

Messenger RNAs have been reported to play an essential role in GC cancer. For example, high NRP1 expression and hypermethylation were associated with poor GC prognosis (Wang et al., 2019), whereas another study indicated that it could be an anti-tumor target (Grandclement and Borg, 2011). In addition, high DUSP1 expression levels were found to promote progression, drug resistance, and poor prognosis of GC (Teng et al., 2018). On the other hand, SOX14 showed opposite prognostic values in cervical cancer and leukemia, with anti-tumor and carcinogenic effects, respectively (Stanisavljevic et al., 2017; Tosic et al., 2018). Studies have also implicated CPNE8, MAGED1, and RNF144A in ovarian and breast cancers (Zeng et al., 2012; Nagasawa et al., 2019; Yang et al., 2019). However, LOC91450 and ZNF22 have not been reported in cancer.

Accumulating evidence indicates that DNA methylation plays a significant role in cancer progression. However, only a handful of studies have described the relationship between levels of single-site methylation and GC prognosis. Particularly, high expressions of MAP7D2, ACOT13, EID1, RPS4X, and TTRAP have been associated with poor prognosis in gastric, lung, and pancreatic cancers as well as hepatic carcinoma, respectively, while a high TTRAP expression reportedly inhibits the growth of osteosarcoma (Kamio et al., 2010; Zhou et al., 2013; Kuang et al., 2017; Liu K.T. et al., 2018; Liu X. et al., 2018). Notably, the relationship between methylation levels and corresponding gene expression profiles is unknown, necessitating further research. Furthermore, the remaining DNA methylation sites and their corresponding genes have not been reported. Lastly, no study has described the prognostic significance using a genome-wide network.

Last but not least, in general, no study concerning their prognostic significance as a genome-wide network has been reported yet.

Tumorigenesis involves multiple interacting biological processes. In addition, an integrated genetic network is better at reflecting intra-tumor heterogeneity compared to a single biomarker. In the current study, we identified a novel, prognostic, signature genome-wide network, consisting of seven miRNAs, eight mRNA, and 19 DNA methylation sites after screening the entire TCGA cohort using training and random cohorts. This network was further divided into several other models.

Our results revealed that the integrative signature was an independent prognostic factor for survival in GC patients and performed better than any single biomarker or clinical characteristic. Moreover, stratification by other clinicopathological features, such as stage, age, sex, primary site, pathology grade, Lauren classification, and treatments, revealed significantly different prognosis values based on different GSs. In addition, staging was still an effective prognostic factor after dividing into low- and high-genomics-score groups, suggesting that GS and traditional staging can complement each other, and the genetic network could add prognostic value to traditional staging. Exclusion of patients with I staging showed that chemotherapy is a significant prognostic factor because I staging does not always need additional chemotherapy for effective prognosis.

We also developed and validated nomograms based on the GS. Particularly, results from ROC and DCA indicated that all of them had significantly better predictive performances than the traditional clinical nomogram. Comprehensive property (similar C-index) was not significantly different in genomics nomogram and Cox-model 1 nomogram, and compared to the genomics nomogram, Cox-model 1 nomogram had fewer biomarkers. In addition, Cox-model 1 nomogram performed well, with a higher positive net reclassification improvement (NRI). Therefore, Cox-model 1 nomogram might be more suitable for clinical application, which deserved further study. Besides that, the Cox-model 2 nomogram had the least feature (nine biomarkers) including miRNAs, mRNA, and DNA methylation sites for constructing a genome-wide network, while it had a lower C-index and a negative NRI. The other six models showed a relatively poor performance in ROC or DCA, with limited application value. What is more, it is possible that DNA methylation was the highest contributor to the survival prediction of this gene network. We suspect that this may be related to the larger number of DNA methylation sites compared to miRNA and mRNA.

We adopted GO and KEGG analyses to assess the influence of genome-wide network in the prognosis of GC. Generally, biological processes mainly involve various biological functions, such as methylation, phosphorylation, and endocrine regulation. Methylation pathway was related to the occurrence and the development of GC, which was consistent with our results. Besides that, functional enrichment analysis revealed that phosphorylation pathway was significantly enriched as well, which got more and more attention these years. On the other hand, the main components of participation included organelles, cytomembranes, extrinsic to membranes, nuclear and synapses, whereas molecular functions comprise nucleoside, ATP, RNA, and transcription factor binding as well as activity of various enzymes. Abnormal cell composition is closely related to the development of tumor. The abnormal protein may act on the nucleus, membrane, or cell matrix, thereby leading to the progression of cancer, such as NRP1 protein (Wang et al., 2019). In the current study, KEGG analysis indicated that the gene network function was a relevant pathway in cancer, cell cycle, and adipocytokine signaling, while the other pathways had been reported in small cell lung and bladder cancers. Further experiments to reveal the biological function of this gene network are needed.

We also employed a series of complex statistical analyses to construct and validate a genome-wide network based on different biomarkers and then divided it into different models. We recommend the resulting GS despite it not being an absolute representative of tumor heterogeneity. This network could complement the deficiency of traditional staging and generate a more accurate prediction of survival rates in GC patients. Additionally, it effectively distinguishes patients who could benefit from chemotherapy, thereby reducing unnecessary treatments. It is also possible that the network could be used to identify novel therapeutic targets for GC, although this requires further investigation.

Limitation

This study had several limitations. Firstly, information relating to patient co-morbidities and performance status was not available in the TCGA database. Secondly, the systemic chemotherapy regimens were not uniform, and most of them were based on fluoropyrimidines. Thirdly, the gene network contains too many biomarkers, increasing the difficulty of clinical use. Lastly, this was a retrospective study, without any independent external patient datasets as test. Despite some limitations, it was the first, to the best of our knowledge, to integrate miRNAs, mRNA, and DNA methylation sites as a genome-wide network to predict the OS of patients with GC, and we would try to design a validation in our hospital.

Conclusion

In summary, we used a TCGA cohort to develop and validate a novel genome-wide network comprising seven miRNAs, eight mRNAs, and 19 DNA methylation sites for the prognosis of GC. A combination of GS and TNM staging enhances its prognostic value, proposing a more comprehensive sub-typing system. The developed network is expected to aid in predicting GC patients who may benefit from chemotherapy to some degree.

Data Availability Statement

The datasets generated and analyzed during the current study are available in the TCGA database [https://portal.gdc.cancer.gov/].

Ethics Statement

Since TCGA was a public-use database, no additional permission was required from the Ethics Committee. In addition, this study was deemed exempt from institutional review board approval by Nanfang Hospital of Southern Medical University (Guangzhou, China).

Author Contributions

All the authors listed had made a substantial contribution to this work. YJ and GL put forward the conception and designed the study. WH, TL, and MZ collected and collated the data. ZS, HC, and ZH analyzed the data and wrote the manuscript together. HL, JY, and YH made contribution to proofread the article. Finally, all the authors took responsibility of the final manuscript and approved it for publication.

Funding

This work was supported by grants from the National Natural Science Foundation of China (81872013 and 81672446) and The Outstanding Youths Development Scheme of Nanfang Hospital, Southern Medical University (2018J007).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We would like to thank the staff members of TCGA program.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2020.00835/full#supplementary-material

FIGURE S1 | Screening process for the genome-wide network.

FIGURE S2 | Screening process for the Gene Ontology analysis and Kyoto Encyclopedia of Genes and Genomes pathways.

FIGURE S3 | Least absolute shrinkage and selection operator Cox regression performed for miRNAs, mRNA, and DNA methylation sites.

FIGURE S4 | A heat map based on genomic scores layered by clinicopathological factors.

FIGURE S5 | Median survival for patients in the low-genomics score (GS) and the high-GS groups.

FIGURE S6 | Distribution of patient cases and density based on genomics score in the validation group (A,B). Scatter plots of genomics score (GS) regarding the classification of low- and high-GS (C,D), and the bold line represents the median.

FIGURE S7 | Kaplan–Meier curve of overall survival for the low- and the high-genomics-score groups divided by age (young and old), sex (male and female), primary site (cardia, fundus/body, and antrum), pathology grade (I–II and III–IV), Lauren classification (intestinal type and diffused type), and treatment (R0 surgery, R1 surgery, chemotherapy, and no chemotherapy) in the training group.

FIGURE S8 | Kaplan–Meier curve of overall survival for the low- and the high-genomics-score groups divided by age (young and old), sex (male and female), primary site (cardia, fundus/body, and antrum), pathology grade (I-II and III-IV), Lauren classification (intestinal type and diffused type), and treatment (R0 surgery, R1 surgery, chemotherapy, and no chemotherapy) in the validation group.

FIGURE S9 | Cox-model 1 nomogram and its calibration plot in the training group and the validation group.

FIGURE S10 | Cox-model 2 nomogram and its calibration plot in the training group and the validation group.

FIGURE S11 | Clinical nomogram and its calibration plot in the training group and the validation group.

FIGURE S12 | Time-dependent receiver operating characteristic curves on 1, 3, and 5 years for each nomogram in the training group.

FIGURE S13 | Time-dependent receiver operating characteristic curves on 1, 3, and 5 years for each nomogram in the validation group.

FIGURE S14 | Venn diagrams for the target genes of each miRNA.

FIGURE S15 | Three most potential target genes for each miRNA (A) and the miRNA–target genes co-expression network (B).

FIGURE S16 | Time-dependent receiver operating characteristic curves on 1, 3, and 5 years for each signature in the genome-wide network. Each color represents a different signature: red 2 (genomics nomogram), orange (clinical nomogram), yellow (hsa-mir-100), blue (hsa-mir-1304), black (hsa-mir-193b), brown (hsa-mir-22), gray (hsa-mir-136), pale violet red 1 (hsa-mir-653), purple 2 (hsa-mir-6808), green 2 (NRP1|8829), dark red (SOX14|8403), dark orchid 3 (CPNE8|144402), dark orange 2 (MAGED1|9500), cyan (RNF144A|9781), corn silk 4 (ZNF22|7570), hot pink 1 (DUSP1|1843), chartreuse (LOC91450|91450), aquamarine 2 (cg02223323), dark green (cg00481239), maroon 1 (cg07020967), pale green (cg08859156) gold (cg12485556), purple 1 (cg14791193), antique white (cg15861578), khaki 2 (cg15486740) sea green 2 (cg20100408), light sky blue (cg20350671), maroon (cg22395807), Indian red 2 (cg24361571), tan 1 (cg25361506), peach puff (cg25622155), turquoise 1 (cg25161386), orchid 2 (cg22740006), dark golden rod 2 (cg22813794), dark golden rod 4 (cg26014401), and violet red 1 (cg26856948).

FIGURE S17 | Kaplan–Meier curve of overall survival in the low-genomics-score group stratified by clinical features (e.g., stage).

FIGURE S18 | Kaplan-Meier curve of overall survival in the high-genomics-score group stratified by clinical features (e.g., stage).

FIGURE S19 | Kaplan–Meier curve of overall survival in stages II and III patients.

TABLE S1 | Descriptive statistics of gastric cancer patients.

TABLE S2 | Multivariable analysis for miRNAs, mRNA, and DNA methylation sites and clinical characteristics.

TABLE S3 | Multivariable analysis for miRNAs.

TABLE S4 | Multivariable analysis for mRNA.

TABLE S5 | Multivariable analysis for DNA methylation sites.

TABLE S6 | The equation of different models.

TABLE S7 | The C-index value of different models.

TABLE S8 | The area under the curve values for each individual of the genome-wide network.

References

Anna, C. D. S., Junior, A. G. F., Soares, P., Tuji, F., Paschoal, E., Chaves, L. C., et al. (2018). Molecular biology as a tool for the treatment of cancer. Clin. Exper. Med. 18, 457–464.

Google Scholar

Bang, Y. J., Kim, Y. W., Yang, H. K., Chung, H. C., Park, Y. K., Lee, K. H., et al. (2012). Adjuvant capecitabine and oxaliplatin for gastric cancer after D2 gastrectomy (CLASSIC): a phase 3 open-label, randomised controlled trial. Lancet 379, 315–321. doi: 10.1016/s0140-6736(11)61873-4

CrossRef Full Text | Google Scholar

Camargo, M. C., Figueiredo, C., and Machado, J. C. (2019). Review: gastric malignancies: basic aspects. Helicobacter 24(Suppl. 1):e12642.

Google Scholar

Cancer Genome Atlas Research (2014). Comprehensive molecular characterization of gastric adenocarcinoma. Nature 513, 202–209. doi: 10.1038/nature13480

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, J., Zheng, B., Wang, C., Chen, Y., Du, C., Zhao, G., et al. (2014). Prognostic role of microRNA-100 in various carcinomas: evidence from six studies. Tumor Biol. 35, 3067–3071. doi: 10.1007/s13277-013-1398-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, X., Huang, Z., and Chen, R. (2018). Microrna-136 promotes proliferation and invasion ingastric cancer cells through Pten/Akt/P-Akt signaling pathway. Oncol. Lett. 15, 4683–4689.

Google Scholar

Choi, R. S., DiNardo, J. A., and Brown, M. L. (2019). Current and future molecular diagnostics of gastric cancer. Expert. Rev. Mol. Diagn. 19, 863–874.

Google Scholar

Cristescu, R., Lee, J., Nebozhyn, M., Kim, K. M., Ting, J. C., Wong, S. S., et al. (2015). Molecular analysis of gastric cancer identifies subtypes associated with distinct clinical outcomes. Nat. Med. 21, 449–456.

Google Scholar

Duarte, H. O., Gomes, J., Machado, J. C., and Reis, C. A. (2018). Gastric cancer: basic aspects. Helicobacter 23(Suppl. 1):e12523. doi: 10.1111/hel.12523

PubMed Abstract | CrossRef Full Text | Google Scholar

Global Burden of Disease Cancer Collaboration, Fitzmaurice, C., Akinyemiju, T. F., Al Lami, F. H., Alam, T., Alizadeh-Navaei, R., et al. (2019). Global, regional, and national cancer incidence, mortality, years of life lost, years lived with disability, and disability-adjusted life-years for 29 cancer groups, 1990 to 2017: a systematic analysis for the global burden of disease study. JAMA Oncol. 5, 1749–1768.

Google Scholar

Grandclement, C., and Borg, C. (2011). Neuropilins: a new target for cancer therapy. Cancers 3, 1899–1928. doi: 10.3390/cancers3021899

PubMed Abstract | CrossRef Full Text | Google Scholar

Hou, X., He, X., Wang, K., Hou, N., Fu, J., Jia, G., et al. (2018). Genome-wide network-based analysis of colorectal cancer identifies novel prognostic factors and an integrative prognostic index. Cell Physiol. Biochem. 49, 1703–1716. doi: 10.1159/000493614

PubMed Abstract | CrossRef Full Text | Google Scholar

Iasonos, A., Schrag, D., Raj, G. V., and Panageas, K. S. (2008). How to build and interpret a nomogram for cancer prognosis. J. Clin. Oncol. 26, 1364–1370. doi: 10.1200/jco.2007.12.9791

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, Y., Chen, C., Xie, J., Wang, W., Zha, X., Lv, W., et al. (2018a). Radiomics signature of computed tomography imaging for prediction of survival and chemotherapeutic benefits in gastric cancer. EBio Med. 36, 171–182. doi: 10.1016/j.ebiom.2018.09.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, Y., Zhang, Q., Hu, Y., Li, T., Yu, J., Zhao, L., et al. (2018b). ImmunoScore signature: a prognostic and predictive tool in gastric cancer. Ann. Surg. 267, 504–513. doi: 10.1097/sla.0000000000002116

PubMed Abstract | CrossRef Full Text | Google Scholar

Kamio, Y., Maeda, K., Moriya, T., Takasu, N., Takeshita, A., Hirai, I., et al. (2010). Clinicopathological significance of cell cycle regulatory factors and differentiation-related factors in pancreatic neoplasms. Pancreas 39, 345–352. doi: 10.1097/mpa.0b013e3181bb9204

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuang, J., Li, Q. Y., Fan, F., Shen, N. J., Zhan, Y. J., Tang, Z. H., et al. (2017). Overexpression of the X-linked ribosomal protein S4 predicts poor prognosis in patients with intrahepatic cholangiocarcinoma. Oncol. Lett. 14, 41–46. doi: 10.3892/ol.2017.6137

PubMed Abstract | CrossRef Full Text | Google Scholar

Kurata, J. S., and Lin, R. J. (2018). MicroRNA-focused CRISPR-Cas9 library screen reveals fitness-associated miRNAs. RNA 24, 966–981. doi: 10.1261/rna.066282.118

PubMed Abstract | CrossRef Full Text | Google Scholar

Lauren, P. (1965). The two histological main types of gastric carcinoma: diffuse and so-called intestinal-type carcinoma. An attempt at a histo-clinical classification. Acta Pathol. Microbiol. Scand. 64, 31–49. doi: 10.1111/apm.1965.64.1.31

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, X., Zhang, Y., Zhang, Y., Ding, J., Wu, K., Fan, D., et al. (2010). Survival prediction of gastric cancer by a seven-microRNA signature. Gut 59, 579–585. doi: 10.1136/gut.2008.175497

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, G. H., Shi, H., Deng, L., Zheng, H., Kong, W., Wen, X., et al. (2019). Circular RNA circ-FOXM1 facilitates cell progression as ceRNA to target PPDPF and MACC1 by sponging miR-1304-5p in non-small cell lung cancer. Biochem. Biophys. Res. Commun. 513, 207–212. doi: 10.1016/j.bbrc.2019.03.213

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, K. T., Yeh, I. J., Chou, S. K., Yen, M. C., and Kuo, P. L. (2018). Regulatory mechanism of fatty acid-CoA metabolic enzymes under endoplasmic reticulum stress in lung cancer. Oncol. Rep. 40, 2674–2682.

Google Scholar

Liu, X., Wu, J., Zhang, D., Bing, Z., Tian, J., Ni, M., et al. (2018). Identification of potential key GENES associated with the pathogenesis and prognosis of gastric cancer based on integrated bioinformatics analysis. Front. Genet. 9:265. doi: 10.3389/fonc.2019.00265

PubMed Abstract | CrossRef Full Text | Google Scholar

Mu, Y. P., Tang, S., Sun, W. J., Gao, W. M., Wang, M., Su, X. L., et al. (2014). Association of miR-193b down-regulation and miR-196a up-regulation with clinicopathological features and prognosis in gastric cancer. Asian Pac. J. Cancer Prev. 15, 8893–8900. doi: 10.7314/apjcp.2014.15.20.8893

PubMed Abstract | CrossRef Full Text | Google Scholar

Nagasawa, S., Ikeda, K., Horie-Inoue, K., Sato, S., Itakura, A., Takeda, S., et al. (2019). Systematic identification of characteristic genes of ovarian clear cell carcinoma compared with high-grade serous carcinoma based on RNA-sequencing. Int. J. Mol. Sci. 20:4330. doi: 10.3390/ijms20184330

PubMed Abstract | CrossRef Full Text | Google Scholar

Nagtegaal, I. D., Odze, R. D., Klimstra, D., Paradis, V., Rugge, M., and Schirmacher, P. (2019). The 2019 WHO classification of tumours of the digestive system. Histopathology 76, 182–188. doi: 10.1111/his.13975

PubMed Abstract | CrossRef Full Text | Google Scholar

Pan, Y., Xu, T., Liu, Y., Li, W., and Zhang, W. (2019). Upregulated circular RNA circ_0025033 promotes papillary thyroid cancer cell proliferation and invasion via sponging miR-1231 and miR-1304. Biochem. Biophys. Res. Commun. 510, 334–338. doi: 10.1016/j.bbrc.2019.01.108

PubMed Abstract | CrossRef Full Text | Google Scholar

Serra, O., Galán, M., Ginesta, M. M., Calvo, M., Sala, N., Salazar, R., et al. (2019). Comparison and applicability of molecular classifications for gastric cancer. Cancer Treat. Rev. 77, 29–34. doi: 10.1016/j.ctrv.2019.05.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Sohn, B. H., Hwang, J. E., Jang, H. J., Lee, H. S., Oh, S. C., Shim, J. J., et al. (2017). Clinical significance of four molecular subtypes of gastric cancer identified by the cancer genome atlas project. Clin. Cancer Res. 23, 4441–4449. doi: 10.1158/1078-0432.ccr-16-2211

PubMed Abstract | CrossRef Full Text | Google Scholar

Song, B., Du, J., Song, D. F., Ren, J. C., and Feng, Y. (2018). Dysregulation of NCAPG, KNL1, miR-148a-3p, miR-193b-3p, and miR-1179 may contribute to the progression of gastric cancer. Biol. Res. 51:44.

Google Scholar

Stanisavljevic, D., Petrovic, I., Vukovic, V., Schwirtlich, M., Gredic, M., Stevanovic, M., et al. (2017). SOX14 activates the p53 signaling pathway and induces apoptosis in a cervical carcinoma cell line. PLoS One 12:e0184686. doi: 10.1371/journal.pone.0184686

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, Z., Liu, H., Yu, J., Huang, W., Han, Z., Lin, T., et al. (2019a). Frequency and prognosis of pulmonary metastases in newly diagnosed gastric cancer. Front. Oncol. 9:671. doi: 10.3389/fonc.2019.00671

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, Z., Zheng, H., Yu, J., Huang, W., Li, T., Chen, H., et al. (2019b). Liver metastases in newly diagnosed gastric cancer: a population-based study from SEER. J. Cancer 10, 2991–3005. doi: 10.7150/jca.30821

PubMed Abstract | CrossRef Full Text | Google Scholar

Tan, I. B., Ivanova, T., Lim, K. H., Ong, C. W., Deng, N., Lee, J., et al. (2011). Intrinsic subtypes of gastric cancer, based on gene expression pattern, predict survival and respond differently to chemotherapy. Gastroenterology 141, 476–485.

Google Scholar

Teng, F., Xu, Z., Chen, J., Zheng, G., Zheng, G., Lv, H., et al. (2018). DUSP1 induces apatinib resistance by activating the MAPK pathway in gastric cancer. Oncol. Rep. 40, 1203–1222.

Google Scholar

Tibshirani, R. (1997). The lasso method for variable selection in the Cox model. Stat. Med. 16, 385–395. doi: 10.1002/(sici)1097-0258(19970228)16:4<385::aid-sim380>3.0.co;2-3

CrossRef Full Text | Google Scholar

Torre, L. A., Bray, F., Siegel, R. L., Ferlay, J., Lortet-Tieulent, J., Jemal, A., et al. (2015). Global cancer statistics, 2012. CA Cancer J. Clin. 65, 87–108. doi: 10.3322/caac.21262

PubMed Abstract | CrossRef Full Text | Google Scholar

Tosic, N., Petrovic, I., Grujicic, N. K., Davidovic, S., Virijevic, M., and Vukovic, N. S. (2018). Prognostic significance of SOX2, SOX3, SOX11, SOX14 and SOX18 gene expression in adult de novo acute myeloid leukemia. Leuk Res. 67, 32–38. doi: 10.1016/j.leukres.2018.02.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Ueda, T., Volinia, S., Okumura, H., Shimizu, M., Taccioli, C., Rossi, S., et al. (2010). Relation between microRNA expression and progression and prognosis of gastric cancer: a microRNA expression analysis. Lancet Oncol. 11, 136–146. doi: 10.1016/s1470-2045(09)70343-2

CrossRef Full Text | Google Scholar

Vickers, A. J., and Elkin, E. B. (2006). Decision curve analysis: a novel method for evaluating prediction models. Med. Decis. Mak. 26, 565–574. doi: 10.1177/0272989x06295361

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, G. H., Shi, B., Fu, Y., Zhao, S., Qu, K., Guo, Q., et al. (2019). Hypomethylated gene NRP1 is co-expressed with PDGFRB and associated with poor overall survival in gastric cancer patients. Biomed. Pharmacother. 111, 1334–1341. doi: 10.1016/j.biopha.2019.01.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, Y. L., Zhang, Y., Li, D. D., Zhang, F. L., Liu, H. Y., Liao, X. H., et al. (2019). RNF144A functions as a tumor suppressor in breast cancer through ubiquitin ligase activity-dependent regulation of stability and oncogenic functions of HSPA2. Cell Death Differ. 27, 1105–1118. doi: 10.1038/s41418-019-0400-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Zeng, Z. L., Wu, W. J., Yang, J., Tang, Z. J., Chen, D. L., Qiu, M. Z., et al. (2012). Prognostic relevance of melanoma antigen D1 expression in colorectal carcinoma. J. Transl. Med. 10:181. doi: 10.1186/1479-5876-10-181

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Yang, W., Li, D., Yang, J. Y., Guan, R., Yang, M. Q., et al. (2018). Toward the precision breast cancer survival prediction utilizing combined whole genome-wide expression and somatic mutation analysis. BMC Med. Genom. 11(Suppl. 5):104. doi: 10.1186/s12920-018-0419-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Z., Dong, Y., Hua, J., Xue, H., Hu, J., Jiang, T., et al. (2019). A five-miRNA signature predicts survival in gastric cancer using bioinformatics analysis. Gene 699, 125–134. doi: 10.1016/j.gene.2019.02.058

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, J., Ge, P., Liu, X., Wei, J., Wu, G., Li, X., et al. (2017). MiR-136 inhibits gastric cancer-specific peritoneal metastasis by targeting HOXC10. Tumour Biol. 39:1010428317706207.

Google Scholar

Zhou, C., Shen, Q., Xue, J., Ji, C., and Chen, J. (2013). Overexpression of TTRAP inhibits cell growth and induces apoptosis in osteosarcoma cells. BMB Rep. 46, 113–118. doi: 10.5483/bmbrep.2013.46.2.150

PubMed Abstract | CrossRef Full Text | Google Scholar

Zuo, Q. F., Cao, L. Y., Yu, T., Gong, L., Wang, L. N., Zhao, Y. L., et al. (2015). MicroRNA-22 inhibits tumor growth and metastasis in gastric cancer by directly targeting MMP14 and Snail. Cell Death Dis. 6:e2000. doi: 10.1038/cddis.2015.297

PubMed Abstract | CrossRef Full Text | Google Scholar