Research Paper Volume 12, Issue 17 pp 17418—17435
Identification of subtype-specific genes signature by WGCNA for prognostic prediction in diffuse type gastric cancer
- 1 Department of General Surgery, The Second Affiliated Hospital of Nanchang University, Nanchang, Jiangxi Province, China
Received: April 6, 2020 Accepted: July 7, 2020 Published: September 11, 2020https://doi.org/10.18632/aging.103743
How to Cite
Copyright: © 2020 Zhou et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Background: Gastric cancer is a common malignancy and had poor response to treatment due to its strong heterogeneity. This study aimed to identify essential genes associated with diffuse type gastric cancer and construct a powerful prognostic model.
Results: We conducted a weighted gene co-expression network analysis (WGCN) using transcripts per million (TPM) expression data from The Cancer Genome Atlas (TCGA) to find out the module related with diffuse type gastric cancer. Combining Least Absolute Shrinkage and Selection Operator (LASSO) with multi-cox regression, the 10 specific genes risk score model of diffuse type gastric cancer was established. The concordance index (0.97), the area under the respective ROC curves (AUCs) (1-years: 0.98; 3-years: 1; 5-years: 1) and survival difference of high- and low risk groups (p=2.84e-10) of this model in TCGA dataset were obtained. The moderate predicting performance was observed in the independent cohort of GSE15459 and GSE62254. The results of the gene set enrichment analysis (GSEA) using high-and low risk group as phenotype indicated differential expression of tumor-related pathways.
Conclusion: Thus, we constructed a reliable prognostic model for diffuse type gastric cancer, which should be beneficial for clinical therapeutic decision-making.
Gastric cancer (GC) is the sixth most common cancer and the second leading cause of cancer-related deaths worldwide. In 2018, 1,033,701 people were diagnosed and 782,685 people died from GC . Patients with GC usually have an unfavorable prognosis, as the majority reach the advanced stages of disease prior to diagnosis . The lack of precision treatment and evaluation strategies have prompted researchers to investigate carcinogenic abnormalities of GC to assess survival rates and guide medical decisions. Identifying therapeutic targets and prognostic biomarkers for early GC and developing appropriate therapeutic methods are a prospective method to identify the subtypes of GC and improve the prognosis of patients with advanced gastric cancer. However, the underlying heterogeneity and complexity of GC make it difficult to identify reliable factors for effective clinical treatment [3, 4].
GC has many classification systems, such as: the Lauren classification system, and the World Health Organization (WHO) classification systems [5, 6]. Lauren classification mainly includes intestinal type and diffuse type. Because of its strong perceiving of histology and biology of gastric cancer, it has been widely used in clinical practice . Intestinal type GC cells are tubular or glandular, more densely arranged and more cohesive, whereas diffuse type GC cells are usually diffuse and have poor adhesion, resulting in less glandular formation and easier diffusion. The two histological types differ in their clinical and molecular features to the point of representing distinct entities . Diffuse type GC usually had characteristic mutations in genes that participate in adhesion, chromatin integrity, or cell motility . Intestinal type GC exhibited aneuploidy or other genetic features more frequently.
In the past, it was not accurate to analyze GC as a whole. In this study, we divided the samples from The Cancer Genome Atlas (TCGA; https://cancergenome.nih.gov/) database to diffuse- and intestinal- type GC. we identified genes significantly associated with diffuse type based on weighted gene co-expression network analysis (WGCNA). Ten genes were obtained to construct the predicting system, which was proved an effective prognostic system for diffuse type GC.
Detection of gene co-expression modules correlated with diffuse type GC cohort
The data was processed and analyzed by following the workflow in Figure 1. Top 50% most variable genes (9752 genes) were used for WGCNA. An obvious outlier was removed (Supplementary Figure 1A) and a soft threshold = 4 was selected to construct a scale-free network (Supplementary Figure 1B, 1C). A total of 30 gene modules were identified after setting the minimum cluster size as 30 (Figure 2A). The grey module contained genes not attributed to any modules. The identified genes associated with a clinical trait were of great value in the exploration of the molecular characteristics of that trait. In the present study, the clinical parameters of diffuse type GC patients, including age at diagnosis, gender, race, tumor differentiation, pathological stage, pathological T category, pathological N category, M category, lymphatic metastasis, cancer status, helicobacter pylori infection, status and overall survival (OS) time was extracted for analysis. we found that blue module (R=0.26, p=3e-07) and green module (R=0.27, p=7e-08) were significantly associated with diffuse type GC (Figure 2B–2D) and the blue (R=0.22, p=3e-05) and green (R=0.26, p=4e-07) module were also significantly correlated with tumor differentiation.
Figure 2. Identification of modules associated with the diffuse type GC. (A) Dendrogram of 9752 genes clustered based on a dissimilarity measure (1-TOM). (B) Correlation of module eigengenes with all traits. Each unit contains the corresponding correlation coefficient and P-value. The table is color-coded by correlation according to the color legend. (C, D) Scatter diagrams between blue and green modules and diffuse type GC.
Hub genes in blue and green module
Before conducting the univariate cox analysis, we first adopted the merge function in R to integrate the expression profiles of the 1794 module genes with corresponding 71 diffuse type GC patients’ survival time and status information. 63 overall survival (OS)-related hub genes (p<0.1) were picked out by univariate cox regression analysis. According to the characteristics of variable selection and regularization, while fitting the generalized linear model, LASSO regression was performed to select hub genes for predicting the prognosis of high-performance patients (Figure 3A, 3B). This approach is popular in machine learning and is implemented through the “glmnet” package. 16 hub genes ("RALA", "DDX3Y", "ERP29", "SRSF5", "SLC9A3R1", "FBXO9", "GMDS", "CCNI", "LEF1", "RFX5", "CAST", "ELMO1", "FRZB", "TMEM92", "SELP" and "NMB") were identified.
Figure 3. Identification of hub genes using LASSO regression. (A) The trajectory of each prognosis-related candidate gene’s coefficient in diffuse type GC was observed in the LASSO coefficient profiles with the changing of the lambda in LASSO algorithm. (B) After the 10-fold cross-validation, a confidence interval was got for partial likelihood deviance as the lambda changed.
Construction of risk score model based on multivariate Cox regression
First, followed by multivariate cox regression, the optimal 10 prognostic signatures in diffuse type GC samples (Table 1), including: “RALA”, “DDX3Y”, “SRSF5”, “SLC9A3R1”, “GMDS”, “LEF1”, “RFX5”, “CAST”, “FRZB” and “SELP”, were picked out. Then, the risk score (RS) model for OS was identified. All 71 diffuse type GC samples were endowed with a RS by RS model calculation and divided into high- and low- RS groups using the median value as the cut-off point. The concordance index (C-index) of this model was 0.97, indicating that this model had quite high reliability. Figure 4A showed that patients in the low-risk group had longer OS (p < 2.8e-10) than those of the high-risk group. The relation of high- and low expression of 10 genes to OS were viewed in Supplementary Figures 2–11. To determine the predictive accuracy of this prognostic model, we performed a receiver operating characteristic (ROC) curve analysis, which demonstrated that the area under the curve (AUC) was 0.98 for 1-year survival, 1 for 3-year survival and 1 for 5-year (Figure 4B). And the calibration curve supported the predicting model of 10 genes (Figure 4C–4E). The survival, 10 genes expression and risk score for TCGA samples could be viewed in Supplementary Table 1.
Figure 4. The prognostic performance of the 10 genes model in the TCGA-STAD. (A) Survival analysis of the high-risk group and the low-risk group using Kaplan–Meier curves. (B) The prognostic efficiency of the 10 genes model for survival time. ROC curves of the 10 genes signature for predicting 1-, 3 -and 5- year survival were analyzed. (C–E) The comparison between predicted and actual outcome for 1-, 3-, and 5-year survival probabilities was showed in the calibration plots.
Table 1. The univariate and multivariate Cox regression analysis between 10 markers and OS in diffuse type GC.
|Notes: Hazard ratio (HR), 95% confidence intervals (95%CI)|
Validation of 10-genes RS model in external independent cohort
75 diffuse type GC patients from GSE15459 were remaining after removing the non-diffuse type and survival time as zero. Consistent with the results in the TCGA cohort, the low RS group had good performance in OS (P=0.00633473) than in the high RS group (Figure 5A). Moreover, the AUCs for 1-year, 3-year and 5-year survival in the validation cohort were 0.717, 0.727 and 0.7 (Figure 5B), respectively. There were 135 diffuse type GC patients in GSE62254 by the above same processing. The low RS group had good performance in OS (P=2E-8) than the high RS group as same as above (Figure 5C). the AUCs for 1-year, 3-year and 5-year survival in the validation cohort were 0.661, 0.752 and 0.758 (Figure 5D). Thus, the RS model had again been proved to be reliable. The survival, 10 genes expression and risk score for GSE15459 and GSE62254 samples could be viewed in Supplementary Tables 2 and 3.
Figure 5. The prognostic performance of the 10 genes model in the GSE15459 and GSE62254. (A) Survival analysis of the high-risk group and the low-risk group using Kaplan-Meier curves in the GSE15459. (B) The prognostic efficiency of the 10 genes model for survival time. ROC curves of the 10 genes signature for predicting 1-, 3 -and 5- year survival were analyzed in the GSE15459. (C) Survival analysis of the high-risk group and the low-risk group using Kaplan-Meier curves in the GSE62254. (D) The prognostic efficiency of the 10 genes model for survival time. ROC curves of the 10 genes signature for predicting 1-, 3 -and 5- year survival were analyzed in the GSE62254.
In order to explore the difference of functions and pathways of high- and low- RS groups, the gene set enrichment analysis (GSEA) was performed using the risk score as the reference phenotype. The GSEA analysis revealed that the complement and coagulation cascades, neuroactive ligand-receptor interaction, hypertrophic cardiomyopathy (HCM), steroid hormone biosynthesis and dilated cardiomyopathy were upregulated in the high-risk group (Figure 6A). Then, the spliceosome, one carbon pool by folate, nucleotide excision repair, cell cycle, RNA degradation, mismatch repair, DNA replication, ubiquitin mediated proteolysis, homologous recombination, p53 signaling pathway, basal transcription factors and base excision repair were upregulated in the low-risk group (Figure 6B). To investigate the enrichment pathways of 10 prognostic genes, we divided the samples into two groups, a group samples with the expression of upper quantile for one of 10 prognostic genes and another group samples with the expression of lower quantile for one of 10 prognostic genes. All 10 genes related pathway enrichment could be viewed in Supplementary Table 4.
Gastric cancer is a common malignancy and had poor response to treatment for its strong heterogeneity. Previous studies investigated gastric cancer samples as a whole  , which made researchers miss some important information and even obtaining imprecise conclusion. As stated above, for Lauren extensively used in clinical practice for the ability of perceiving the histology and biology of GC, diffuse type GC patients usually progressed faster after the diagnosis or were diagnosed in the late stage. Since the significant heterogeneity of two type, it was urgent to find a predicting model for OS of diffuse type GC facilitating clinical decision-making.
WGCNA, the algorithm aiming to investigate the relationships between genes and phenotype of samples, can be applied to identify complex biological mechanisms responsible for the target phenotypes. The unsupervised hierarchical clustering method selected by WGCNA avoided potential biases and subjective decisions attributed to the selection of the candidate genes previously reported associated with diffuse type GC. We applied a system biology approach, namely WGCNA, to analyze TPM expression dataset to identify the networks and genes associated with diffuse type GC. The lasso regression algorithm as the precision and efficiency of variable selection reduced the dimension of model and cox regression algorithm was performed to identify 10 prognostic genes model. The C-index, AUCs and survival difference of high- and low- risk groups well demonstrated perfect performance of predicting survival of model in diffuse type GC. And the similar phenomenon was observed in the independent cohort of GSE15459 and GSE62254.
The green module and blue module were associated with measles, HTLV-I (human T-cell lymphotropic virus type I) infection, Epstein-Barr virus (EB) infection, measles, pathways in cancer, focal adhesion, cell adhesion molecules (CAMs), apoptosis and many immune relevant pathways (the detailed KEGG enrichment can be found in Supplementary Table 5), indicated that not only EB but also HTLV-1 and measles might be related with the occurrence and progression of diffuse type GC [11, 12]. However, the study about HTLV-1 and gastric cancer were few [13, 14]. Matsumoto S et.al  drew a conclusion that HTLV-1 infection likely reduced the risk of helicobacter pylori infection and proliferation and, thereby, the risk of gastric cancer. However, helicobacter pylori infection only played an important role in the intestinal type cancer . Reverse effect was possible to be present in the diffuse type GC. Certainly, the real results would need to get in the subsequent investigation. The measles was thought to be relevant to lung cancer, whereas poor evidence supported . Then, immune relevant pathways enrichment was consistent with the previous studies [17–19]. Interestingly, the two hub modules with diffuse type GC were highly correlated with poorly differentiated tumors, that mean diffuse type GC with poor differentiation. The two points proved that the WGCNA approach of identifying genes related with diffuse type GC was enough reliable.
Then, we analyzed the KEGG enrichment between high-and low- RS groups using GSEA. The nucleotide excision repair, cell cycle, RNA degradation, mismatch repair, DNA replication, p53 signaling pathway and base excision repair were upregulated in low RS patients. However, only complement and coagulation cascades, neuroactive ligand receptor interaction, steroid hormone biosynthesis, hypertrophic cardiomyopathy (HCM) and dilated cardiomyopathy (nominal p<0.05) were upregulated in high RS patients. The phenomenon was involved with not enough sample size. The significant enrichment differences may be a classification method for diffuse type GC patients or an idea for identifying the subtype of other malignancies. Thus, we should validate our findings in more sample size in the future. According to GSEA analysis of every hub genes, many pathways associated with tumor were up-regulated or down-regulated between the upper quantile expression group and the lower quantile expression group. But the implied mechanism under these genes needed a large sample size to investigate in the late studies.
Our study has several limitations. There was not enough sample size due to the small percentage of diffuse type GC in all GC. As a retrospective study, the patient cohort was heterogeneous, and the significance and robustness of the results and hub genes in the prognostic assessment should be validated in prospective cohorts.
In conclusion, our study is the first study to screen the characteristic hub genes of diffuse type GC using WGCNA and to construct a prognostic model based on hub genes. The prognostic predictive model of 10 genes was proved to be able to accurately investigate the prognosis of diffuse type GC. This model might be applied to identify the high-risk patients, and assess the prognosis, so as to facilitate the precise treatment in diffuse type GC.
Materials and Methods
Data collection and processing
Public gene-expression data and full clinical annotation were searched in the Gene-Expression Omnibus (GEO) and TCGA database. The procedure used for data set selection in the GEO database was as follows. The following search parameters were used: (cancer) OR tumor) OR carcinoma) OR adenocarcinoma) AND ((gastric) OR Stomach)) AND "Homo sapiens"[porgn: txid9606]. In the initial search, 755 items were recognized. The eligible criteria included that: 1) owning Lauren classification; 2) owning survival information and 3) at least 50 diffuse type gastric cancer. We removed the datasets that don’t meet the criteria by checking them one by one carefully. In total, we gathered three patient cohorts with gastric cancer for this study: TCGA-STAD, GSE62254 and GSE15459.
The TCGA-STAD RNA-seq data and clinical data of the 375 GC samples were downloaded by using the “TCGAbiolinks” package in R , which was used as a training set for prognostic prediction of the multi-gene signature. The RNA-seq data for 19505 genes measured as fragments per kilobase of transcript per million mapped reads (FPKM), which were converted to transcripts per million (TPM) after removing duplicated genes and zero expression genes . We obtained the raw expression and clinical data of GC samples from GSE622254 and GSE15459 via GEO database which was used as a validation set. The raw expression data was processed by log10-transformed. The demographics are listed in Table 2.
Table 2. Basic characteristics of the datasets.
|Age (Mean ± SD)||62.25±10.47||58.58±12.54||59.87±13.20|
|OS(Mean ± SD)||1.23±0.75||3.87±2.73||3.18±3.67|
Weighted Gene co-expression network construction
To find modules of highly correlated with diffuse type GC, WGCNA was performed using the WGCNA R package  and carried out on top 50% most variable genes (9752 genes). An unsupervised co-expression relationship was initially built on the basis of the adjacency matrix of connection strengths by using Pearson’s correlation coefficients for gene pairs. This matrix was increased to β = 4 based on the scale-free topology criterion. The adjacency matrix of gene expression data for GC patients was then clustered using topological overlap matrix analysis. Finally, the dynamic tree cut algorithm was applied to the dendrogram for module identification with the mini-size of module gene numbers set as 30 and a cut height of 0.95. The module eigengenes (MEs) as the first principal component was performed with the expression data for each co-expressed module in all GC samples. The module that had the strongest association with diffuse type GC was selected for further analysis. The WGCNA algorithm was described in detail by Zhang Bin et al. .
The construction of multi-gene signature risk score model
To identify the prognostic genes in the module correlation with diffuse type GC, the univariate Cox regression analysis was applied using “survival” package. The genes with P<0.1 were defined to be related to the over survival. Given the already detected prognostic genes, we further investigated the significant signature associated with survival across the diffuse type GC samples by the LASSO regression model using “glmnet” package. Then, multivariate Cox regression analysis was performed with the “survival” package to screen out independent prognostic factors from these robust markers, which were conducted the risk score model: prognostic score =∑(C×EXPmRNA), where EXP is the TPM value of the gene and C is the regression coefficient for the corresponding gene in multivariate Cox hazard model analysis. ROC plot with AUC and calibration curve were derived to assess predictive significance of the model using “rms”, “survival” and “timeROC” package in R. Additionally, survival time difference between high- and low-risk levels were estimated via Kaplan-Meier analysis. The risk score model was validated with an independent data set GSE15459 and GSE62254.
Functional pathway analysis
To investigate the biological function of different risk score groups, we further conducted GSEA using the risk score as the phenotype . With the GSEA 4.0.3 software via the Java platform, we derived the “c2.cp.kegg.v7.1.symbols.gmt gene sets” as the reference set. We divided the samples into two groups, a group samples with the expression of upper quantile for one of 10 prognostic genes and another group samples with the expression of lower quantile for one of 10 prognostic genes. The enriched signaling pathways with FDR < 0.25 or nominal p < 0.05 were defined as statistically significant.
GC: gastric cancer; WHO: World Health Organization; TCGA: The Cancer Genome Atlas; WGCNA: weighted gene co-expression network analysis; TPM: transcripts per million; FPKM: per million mapped reads; GEO: Gene Expression Omnibus; MEs: module eigengenes; LASSO: Least Absolute Shrinkage and Selection Operator; ROC: receiver operating characteristic; AUC: the area under the respective ROC curves; OS: overall survival; GSEA: the gene set enrichment analysis; HCM: hypertrophic cardiomyopathy; C-index: The concordance index; CAMs: cell adhesion molecules; htlv-1: human T-cell lymphotropic virus type I; EB: Epstein-Barr virus.
Qi Zhou, Li-Qiang Zhou and Shi-Hao Li designed experiments; the preliminary preparation work was carried out by Qi Zhou and Li-Qiang Zhou; Qi Zhou, Li-Qiang Zhou, Yi-Wu Yuan, Shi-Hao Li, Li Liu, Jin-Liang Wang, Deng-Zhong Wu and You Wu contributed to a lot of statistical analysis using R; Qi Zhou edited the manuscript; Shi-Hao Li and Qi Zhou helped perform the analysis with constructive discussions; all authors read and approved the manuscript.
Conflicts of Interest
The authors declare no potential conflicts of interest.
This work was supported by the National Natural Science Foundation of China (No. 81872480, No. 81760549, No. 81560492).
- 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:394–424. https://doi.org/10.3322/caac.21492 [PubMed]
- 2. Orditura M, Galizia G, Sforza V, Gambardella V, Fabozzi A, Laterza MM, Andreozzi F, Ventriglia J, Savastano B, Mabilia A, Lieto E, Ciardiello F, De Vita F. Treatment of gastric cancer. World J Gastroenterol. 2014; 20:1635–49. https://doi.org/10.3748/wjg.v20.i7.1635 [PubMed]
- 3. Tan P, Yeoh KG. Genetics and molecular pathogenesis of gastric adenocarcinoma. Gastroenterology. 2015; 149:1153–62.e3. https://doi.org/10.1053/j.gastro.2015.05.059 [PubMed]
- 4. Wong SS, Kim KM, Ting JC, Yu K, Fu J, Liu S, Cristescu R, Nebozhyn M, Gong L, Yue YG, Wang J, Ronghua C, Loboda A, et al. Genomic landscape and genetic heterogeneity in gastric adenocarcinoma revealed by whole-genome sequencing. Nat Commun. 2014; 5:5477. https://doi.org/10.1038/ncomms6477 [PubMed]
- 5. Fléjou JF. [WHO classification of digestive tumors: the fourth edition]. Ann Pathol. 2011; 31:S27–31. https://doi.org/10.1016/j.annpat.2011.08.001 [PubMed]
- 6. Lauren P. 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. 1965; 64:31–49. https://doi.org/10.1111/apm.1918.104.22.168 [PubMed]
- 7. Van Cutsem E, Sagaert X, Topal B, Haustermans K, Prenen H. Gastric cancer. Lancet. 2016; 388:2654–64. https://doi.org/10.1016/S0140-6736(16)30354-3 [PubMed]
- 8. Shah MA, Khanin R, Tang L, Janjigian YY, Klimstra DS, Gerdes H, Kelsen DP. Molecular classification of gastric cancer: a new paradigm. Clin Cancer Res. 2011; 17:2693–701. https://doi.org/10.1158/1078-0432.CCR-10-2203 [PubMed]
- 9. Liu D, Zhou B, Liu R. A transcriptional co-expression network-based approach to identify prognostic biomarkers in gastric carcinoma. PeerJ. 2020; 8:e8504. https://doi.org/10.7717/peerj.8504 [PubMed]
- 10. Zhang L, Kang W, Lu X, Ma S, Dong L, Zou B. Weighted gene co-expression network analysis and connectivity map identifies lovastatin as a treatment option of gastric cancer by inhibiting HDAC2. Gene. 2019; 681:15–25. https://doi.org/10.1016/j.gene.2018.09.040 [PubMed]
- 11. Shinozaki-Ushiku A, Kunita A, Fukayama M. Update on epstein-barr virus and gastric cancer (review). Int J Oncol. 2015; 46:1421–34. https://doi.org/10.3892/ijo.2015.2856 [PubMed]
- 12. Chen XZ, Chen H, Castro FA, Hu JK, Brenner H. Epstein-barr virus infection and gastric cancer: a systematic review. Medicine (Baltimore). 2015; 94:e792. https://doi.org/10.1097/MD.0000000000000792 [PubMed]
- 13. Hirata T, Nakamoto M, Nakamura M, Kinjo N, Hokama A, Kinjo F, Fujita J. Low prevalence of human T cell lymphotropic virus type 1 infection in patients with gastric cancer. J Gastroenterol Hepatol. 2007; 22:2238–41. https://doi.org/10.1111/j.1440-1746.2006.04740.x [PubMed]
- 14. Matsumoto S, Yamasaki K, Tsuji K, Shirahama S. Human T lymphotropic virus type 1 infection and gastric cancer development in Japan. J Infect Dis. 2008; 198:10–15. https://doi.org/10.1086/588733 [PubMed]
- 15. Liu KS, Wong IO, Leung WK. Helicobacter pylori associated gastric intestinal metaplasia: treatment and surveillance. World J Gastroenterol. 2016; 22:1311–20. https://doi.org/10.3748/wjg.v22.i3.1311 [PubMed]
- 16. Benharroch D, Ariad S, Tadmor N, Nalbandyan K, Lazarev I. Relevance of the measles virus expression in cancer - an update. Pathol Oncol Res. 2016; 22:661–66. https://doi.org/10.1007/s12253-016-0080-7 [PubMed]
- 17. Milne AN, Carneiro F, O’Morain C, Offerhaus GJ. Nature meets nurture: molecular genetics of gastric cancer. Hum Genet. 2009; 126:615–28. https://doi.org/10.1007/s00439-009-0722-x [PubMed]
- 18. Werner M, Becker KF, Keller G, Höfler H. Gastric adenocarcinoma: pathomorphology and molecular pathology. J Cancer Res Clin Oncol. 2001; 127:207–16. https://doi.org/10.1007/s004320000195 [PubMed]
- 19. Ge S, Xia X, Ding C, Zhen B, Zhou Q, Feng J, Yuan J, Chen R, Li Y, Ge Z, Ji J, Zhang L, Wang J, et al. A proteomic landscape of diffuse-type gastric cancer. Nat Commun. 2018; 9:1012. https://doi.org/10.1038/s41467-018-03121-2 [PubMed]
- 20. Colaprico A, Silva TC, Olsen C, Garofano L, Cava C, Garolini D, Sabedot TS, Malta TM, Pagnotta SM, Castiglioni I, Ceccarelli M, Bontempi G, Noushmehr H. TCGAbiolinks: an r/bioconductor package for integrative analysis of TCGA data. Nucleic Acids Res. 2016; 44:e71. https://doi.org/10.1093/nar/gkv1507 [PubMed]
- 21. Wagner GP, Kin K, Lynch VJ. Measurement of mRNA abundance using RNA-seq data: RPKM measure is inconsistent among samples. Theory Biosci. 2012; 131:281–85. https://doi.org/10.1007/s12064-012-0162-3 [PubMed]
- 22. Langfelder P, Horvath S. Fast R functions for robust correlations and hierarchical clustering. J Stat Softw. 2012; 46:i11. https://doi.org/10.18637/jss.v046.i11 [PubMed]
- 23. Zhang B, Horvath S. A general framework for weighted gene co-expression network analysis. Stat Appl Genet Mol Biol. 2005; 4:Article17. https://doi.org/10.2202/1544-6115.1128 [PubMed]
- 24. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005; 102:15545–50. https://doi.org/10.1073/pnas.0506580102 [PubMed]