Abstract
Gastric cancer (GC) is the fifth most common cancer and the third leading cause of cancer-related deaths worldwide. The high mortality might be attributed to delay in detection and is closely related to lymph node metastasis. Therefore, it is of great importance to explore the mechanism of lymph node metastasis and find strategies to block GC metastasis. Messenger RNA (mRNA), microRNA (miRNA) and long non-coding RNA (lncRNA) expression data and clinical data were downloaded from The Cancer Genome Atlas (TCGA) database. A total of 908 differentially expressed factors with variance >0.5 including 542 genes, 42 miRNA, and 324 lncRNA were screened using significant analysis microarray algorithm, and interaction networks were constructed using these differentially expressed factors. Furthermore, we conducted functional modules analysis in the network, and found that yellow and turquoise modules could separate samples efficiently. The groups classified in the yellow and turquoise modules had a significant difference in survival time, which was verified in another independent GC mRNA dataset (GSE62254). The results suggested that differentially expressed factors in the yellow and turquoise modules may participate in lymph node metastasis of GC and could be applied as potential biomarkers or therapeutic targets for GC.
Gastric cancer; Lymph node metastasis; Co-expression network; Modular analysis
Introduction
Gastric cancer (GC) is the fifth most common cancer and the third leading cause of cancer-related deaths worldwide (11. Torre LA, Bray F, Siegel RL, Ferlay J, Lortet-Tieulent J, Jemal A. Global cancer statistics. CA Cancer J Clin 2015; 65: 87–108, doi: 10.3322/caac.21262.
https://doi.org/10.3322/caac.21262...
). The high mortality of GC might be attributed to delay in detection and is closely related to metastasis and recurrence. Lymph node metastasis occurs in 70% of patients with advanced GC (22. Xue L, Lü N, Li A, Zou S, Lin D, He Z, et al. A clinicopathological analysis of gastric lymphoma. Chinese J Pathol 2005; 34: 332–336.,33. Roviello F, Rossi S, Marrelli D, Pedrazzani C, Corso G, Vindigni C, et al. Number of lymph node metastases and its prognostic significance in early gastric cancer: a multicenter Italian study. J Surg Oncol 2006; 94: 275–280, doi: 10.1002/jso.20566.
https://doi.org/10.1002/jso.20566...
). It is an early event in GC metastasis and an independent prognostic factor, and can significantly affect the prognosis of patients (44. Fatourou E, Ziogas D, Baltogiannis G. Moving from lymph node metastasis in gastric cancer to biological markers. World J Surg 2010; 34: 1140–1141, doi: 10.1007/s00268-009-0368-9.
https://doi.org/10.1007/s00268-009-0368-...
). Therefore, predicting, diagnosing and investigating lymph node metastasis in GC is very important for the prognosis and treatment of patients.
The molecular mechanism of lymph node metastasis has been preliminarily clarified, and mainly includes cell migration and degradation of extracellular matrix, tumor cell apoptosis and immune escape, formation of new lymphatic vessels, and other aspects (55. Chan A, Wong B, Lam S-K. Gastric cancer: past, present and future. Can J Gastroenterol 2001; 15: 469–474, doi: 10.1155/2001/850308.
https://doi.org/10.1155/2001/850308...
,66. Nishida T. Treatment of gastric cancer in Asia: the missing link. Lancet Oncol 2009; 10: 1027–1028, doi: 10.1016/S1470-2045(09)70322-5.
https://doi.org/10.1016/S1470-2045(09)70...
). A series of growth factors, cytokines, chemokines, miRNA and long non-coding RNA (lncRNA) associated with lymph node metastasis have been discovered, which were found to interact with each other to form a complex regulatory network and are involved in various processes of lymph node metastasis in GC (77. Ji R-C. Lymphatic endothelial cells, tumor lymphangiogenesis and metastasis: New insights into intratumoral and peritumoral lymphatics. Cancer Metast Rev 2006; 25: 677–694, doi: 10.1007/s10555-006-9026-y.
https://doi.org/10.1007/s10555-006-9026-...
–1313. Liu N, Zuo C, Wang X, Chen T, Yang D, Wang J, et al. miR-942 decreases TRAIL-induced apoptosis through ISG12a downregulation and is regulated by AKT. Oncotarget 2014; 5: 4959, doi: 10.18632/oncotarget.2067.
https://doi.org/10.18632/oncotarget.2067...
). For example, miRNA-375 is downregulated in gastric carcinomas and regulates cell survival by targeting PDK1 and 14-3-3ζ (1414. Tsukamoto Y, Nakada C, Noguchi T, Tanigawa M, Nguyen LT, Uchida T, et al. MicroRNA-375 is downregulated in gastric carcinomas and regulates cell survival by targeting PDK1 and 14-3-3ζ. Cancer Res 2010; 70: 2339–2349, doi: 10.1158/0008-5472.CAN-09-2777.
https://doi.org/10.1158/0008-5472.CAN-09...
), miRNA-7 functions as an anti-metastatic miRNA in GC by targeting insulin-like growth factor-1 receptor (1515. 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 :1363–1372, doi: 10.1038/onc.2012.156.
https://doi.org/10.1038/onc.2012.156...
), and miR-148a contributes to the maintenance of homeostasis in normal stomach tissue and plays an important role in GC invasion by regulating MMP7 expression. (1616. Sakamoto N, Naito Y, Oue N, Sentani K, Uraoka N, Zarni Oo H, et al. MicroRNA‐148a is downregulated in gastric cancer, targets MMP7, and indicates tumor invasiveness and poor prognosis. Cancer Sci 2014; 105: 236–243, doi: 10.1111/cas.12330.
https://doi.org/10.1111/cas.12330...
). As for lncRNAs, the HOTAIR functions as a competing endogenous RNA to regulate HER2 expression by sponging MiR-331-3p (1212. Liu X-h, Sun M, Nie F-q, Ge Y-b, Zhang E-b, Yin D-d, et al. Lnc RNA HOTAIR functions as a competing endogenous RNA to regulate HER2 expression by sponging miR-331-3p in gastric cancer. Mol Cancer 2014; 13: 1, doi: 10.1186/1476-4598-13-1.
https://doi.org/10.1186/1476-4598-13-1...
), HMlincRNA717 may play crucial roles during cancer occurrence and progression (88. Shao Y, Chen H, Jiang X, Chen S, Li P, Ye M, et al. Low expression of lncRNA-HMlincRNA717 in human gastric cancer and its clinical significances. Tumor Biol 2014; 35: 9591–9595, doi: 10.1007/s13277-014-2243-z.
https://doi.org/10.1007/s13277-014-2243-...
9. Wang Y, Zhang D, Wu K, Zhao Q, Nie Y, Fan D. Long noncoding RNA MRUL promotes ABCB1 expression in multidrug-resistant gastric cancer cell sublines. Mol Cell Bio 2014; 34: 3182–3193, doi: 10.1128/MCB.01580-13.
https://doi.org/10.1128/MCB.01580-13...
), and ATB plays an important role in epithelial-mesenchymal transition to promote invasion and metastasis through the TGFβ/miR-200s/ZEB axis, resulting in a poor prognosis in GC (1010. Saito T, Kurashige J, Nambara S, Komatsu H, Hirata H, Ueda M, et al. A long non-coding RNA activated by transforming growth factor-β is an independent prognostic marker of gastric cancer. Ann Surg Oncol 2015; 22: 915–922, doi: 10.1245/s10434-015-4554-8.
https://doi.org/10.1245/s10434-015-4554-...
11. Ueda T, Volinia S, Okumura H, Shimizu M, Taccioli C, Rossi S, et al. Relation between microRNA expression and progression and prognosis of gastric cancer: a microRNA expression analysis. Lancet Oncol 2010; 11: 136–146, doi: 10.1016/S1470-2045(09)70343-2.
https://doi.org/10.1016/S1470-2045(09)70...
). Although researchers have explored the factors that affect lymph node metastasis, most of the studies focus on one or several factors, and the molecular mechanism of lymph node metastasis is still unclear.
Recent developments in bioinformatics and statistical genomics provide biological systems approaches to better understand the organization of the transcriptome and transcriptional regulation. Among all of the systematic biology approaches, gene network analysis is a powerful approach that considers gene interactions. It has been widely applied in gene expression studies of humans and model organisms (1717. Lee I, Date SV, Adai AT, Marcotte EM. A probabilistic functional network of yeast genes. Science 2004; 306: 1555–1558, doi: 10.1126/science.1099511.
https://doi.org/10.1126/science.1099511...
18. Oldham MC, Konopka G, Iwamoto K, Langfelder P, Kato T, Horvath S, et al. Functional organization of the transcriptome in human brain. Nature Neurosci 2008; 11: 1271–1282, doi: 10.1038/nn.2207.
https://doi.org/10.1038/nn.2207...
–1919. Ponomarev I, Wang S, Zhang L, Harris RA, Mayfield RD. Gene coexpression networks in human brain identify epigenetic modifications in alcohol dependence. J Neurosc 2012; 32: 1884–1897, doi: 10.1523/JNEUROSCI.3136-11.2012.
https://doi.org/10.1523/JNEUROSCI.3136-1...
).
In the present study, we used large quantities of messenger RNA (mRNA)-seq and microRNA (miRNA)-seq data in GC patients in The Cancer Genome Atlas (TCGA) database to screen mRNA, miRNA, and lncRNA with differential expression between the samples with and without lymph node metastasis, and then construct the co-expression networks based on the differential expression of various factors. The network module and Cox regression model were combined to screen survival-related genes. Moreover, the correlation between different expression levels of these genes and prognosis of GC patients was verified in another independent dataset.
Material and Methods
Data sources
Samples of gastric cancer were selected from the TCGA database (http://cancergenome.nih.gov/), in which mRNAs data and miRNAs data were profiled from the Illumina platform. The lncRNAs data were annotated from mRNA transcriptomic database through matching to the HGNC database. In total, 396 samples with mRNA-miRNA-lncRNA paired samples were obtained, and 356 samples remained after removing 40 samples with unclear status in lymph node metastasis. A total of 210 lymph node metastasis samples and 146 non-lymph node metastasis samples were included.
Data preprocessing
After download from TCGA, the expression profiles in level 3 were merged and prepared. The mRNA, miRNAs and lncRNAs with ≥20% missing values were removed, while those with <20% missing values were replaced by mean values. The values of lncRNA and mRNA were assessed by PRKM, and miRNA values were assessed by RPM. Next, these values were transformed by log2 logarithms to obtain Gaussian distribution.
Identification of differentially expressed genes and functional enrichment analysis
Significant analysis microarray (SAM) (2020. Larsson O, Wahlestedt C, Timmons JA. Considerations when using the significance analysis of microarrays (SAM) algorithm. BMC Bioinformatics 2005; 6: 129, doi: 10.1186/1471-2105-6-129.
https://doi.org/10.1186/1471-2105-6-129...
) is an algorithm used to screen the differentially expressed genes (DEGs). When differential gene expression was simply checked by t-test or variance analysis, high rates of false positive results were produced under repeated tests. However, SAM is effective at correcting the false positive rates through controlling the false discovery rate (FDR) in multiple tests to filter DEGs with significant differences. An absolute value of log2 ratio ≥1.5 and FDR ≤0.001 were set as the threshold for determining the significance of gene expression difference.
The identified DEGs were analyzed in terms of gene ontology (GO) function and pathway enrichment analysis using the DAVID (Database for Annotation, Visualization, and Integrated Discovery) hypergeometric test (2121. Dennis G Jr, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, et al. DAVID: database for annotation, visualization, and integrated discovery. Genome Biol 2003; 4: P3, doi: 10.1186/gb-2003-4-5-p3.
https://doi.org/10.1186/gb-2003-4-5-p3...
).
Construction of co-expression networks and excavation of network modules
Based on information of disease-related gene expression profiles under GC, co-expression networks were constructed by calculating adjacency matrix A of a gene pair using the WGCNA package (<https://cran.r-project.org/web/packages/WGCNA/index.html>) (2222. Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics 2008; 9: 1, doi: 10.1186/1471-2105-9-559.
https://doi.org/10.1186/1471-2105-9-559...
). To calculate the adjacency matrix, an intermediate quantity called the co-expression similarity sij was first defined. The default method defines the co-expression similarity sij as the absolute value of the correlation coefficient between the profiles of nodes i and j:
where xi and xj are the vector expression values of gene i and j, respectively, and Cor is used to evaluate Pearson correlation coefficient of the two vector values.
A weighted network adjacency was defined by raising the co-expression similarity to a power:
With β≥1, the adjacency function calculates the adjacency matrix from expression data. The adjacency in Equation 2 implies that the weighted adjacency aij between two genes is proportional to their similarity on a logarithmic scale, log(aij) = β × log(sij). Pearson correlation coefficient Sij is exponentially transformed into connection coefficient aij, to achieve a reliable network.
Network topology property is taken into account to excavate modules of co-expression networks in WGCNA. This algorithm analyses not only the relationship of two conjoint node genes, but other genes correlated to the two nodes. Connection coefficient aij in co-expression network was turned into weight coefficient Wij by the following formula:
Wij considers overlapping between neighbors of the two conjoint genes i and j. Network modules were excavated after hierarchical cluster analysis of W gene weighted matrix.
Survival analysis
Gastric carcinoma related genes were sorted out on the background of expression profile data, and then were constructed into the co-expression network according to their expression levels, which was next divided into various modules. Samples were hierarchically clustered based on module genes and the significant differences of survival time between samples were analyzed through K-M simple clustering. Finally, analysis of COX single variable regression (2323. Jiang D, Wang H, Li J, Wu Y, Fang M, Yang R. Cox regression model for dissecting genetic architecture of survival time. Genomics 2014; 104: 472–476, doi: 10.1016/j.ygeno.2014.10.002.
https://doi.org/10.1016/j.ygeno.2014.10....
) was carried out between survival time and factors in modules.
Results
Preprocessing of expression profile data
There were 560 miRNAs and 12,474 mRNAs left after preprocessing. Totally 1,165 lncRNAs were obtained from transcripts of mRNA expression profiles matched to HGNC database (http://www.genenames.org/). The expression data of mRNA, lncRNA and miRNA after preprocessing are displayed in Figure 1. After standardization, the distribution values among samples were relatively uniform.
Box plot distribution of TCGA data related to gastric carcinoma samples. A, messenger RNA (mRNA) expression; B, long non-coding RNA (lncRNA) expression; C, microRNA (miRNA) expression. Horizontal axis indicates cancer samples; vertical axis indicates the expression value distribution of the mRNAs, lncRNAs and miRNAs.
Screening of differentially expressed factors and functional enrichment analysis
The differentially expressed mRNA, miRNA and lncRNA were screened by SAM algorithm. The result showed that 908 differentially expressed factors with variance >0.5 were screened, including 324 differentially expressed lncRNAs, 42 differentially expressed miRNAs and 542 differential expressed genes (see Supplementary Table S1). These differentially expressed factors divided the samples into two classes as shown in Figure 2.
Heatmap of sample clustering based on differentially expressed factors. Horizontal axis indicates samples; vertical axis indicates differentially expressed factors.
In all, 542 DEGs were functionally enriched by DAVID analysis. There were 15 major enriched GO terms assigned into the biological process categories, including immune response (GO:0006955), response to stress (GO:0006950), defense response (GO:0006952), cell surface receptor signaling pathway (GO:0007166), cell proliferation (GO:0008283), regulation of immune system process (GO:0002682), cell migration (GO:0016477), leukocyte activation (GO:0045321), inflammatory response (GO:0006954), positive regulation of response to stimulus (GO:0048584), regulation of cell proliferation (GO:0042127), cell differentiation (GO:0030154), cellular developmental process (GO:0048869), cell adhesion (GO:0007155), and biological adhesion (GO:0022610) (Table 1, also see Supplementary Table S2).
Construction of co-expression networks and mining of modules
Based on the above mentioned 908 differential factors (including differentially expressed genes, miRNA and lncRNA), co-expression networks were constructed by using WGCNA package in R language. The gene degrees in the co-expression network were submitted to the power distribution law (2424. Mossa S, Barthelemy M, Stanley HE, Amaral LAN. Truncation of power law behavior in “scale-free” network models due to information filtering. Phys Rev Letters 2002; 88: 138701, doi: 10.1103/PhysRevLett.88.138701.
https://doi.org/10.1103/PhysRevLett.88.1...
), which showed in line with the free scale characteristic in biological networks (Figure 3).
Gene degrees in co-expression network and distribution of corresponding numbers. Horizontal axis indicates values of gene degrees (k); vertical axis indicates proportion of gene with k degree, that is, p (k).
Furthermore, we carried out the module mining of genes in the co-expression networks. As shown in Figure 4, these differentially expressed factors in the co-expression networks were divided into four module groups presented by blue, brown, turquoise and yellow color, which included 73, 39, 89, and 35 differentially expressed factors, respectively. In addition, there were other 672 factors that could not be modular, and were presented in grey.
Clustering results of gene modules in co-expression network. Top, hierarchical clustering of 908 genes in co-expression network. Bottom, divisional module classes. Gray zones indicate genes that do not belong to any modules.
Classification validation of factors in four modules
The following analysis was the classification of 356 samples based on the module factors. It was found that yellow and turquoise modules could separate samples efficiently (group 3 was excluded when classified by yellow module because there were rare death samples in the group), and there were significant differences between the survival curves of isolated samples (P<0.01; Figure 5). By contrast, the genes in blue and brown modules could not separate samples effectively, and therefore, these genes were given up in the following survival time analysis.
Sample clustering based on yellow (A) and turquoise (B) module factors and differences in survival curves of the two groups of samples.
Further validation between factors in yellow and turquoise modules and survival time of GC samples was performed by COX univariate regression analysis. The top 7 factors with significant regression in yellow module are shown in Table 2, including adenosine receptor A3 (ADORA3), toll-like receptor 7 (TLR7), interferon regulatory factor (IRF4), CC chemokine receptor 4 (CCR4), reticulon-1 (RTN1), growth factor receptor-bound protein 2 (GRB2)-binding adaptor protein (GAPT), and GRB2-related adapter protein 2 (GRAP2). There were 10 factors with significant regression in turquoise module, including six genes (guanine nucleotide-binding protein G(o) subunit alpha, GNAO1; isthmin-1, ISM1; cartilage intermediate layer protein, CILP; slit homolog 2 protein, SLIT2; scrapie-responsive protein 1, SCRG1; tumor necrosis factor α-induced protein 8 (TNFAIP8)-like protein 3, TNFAIP8L3), two miRNAs (hsa-mir-183 and hsa-mir-942) and two lncRNA (MIR345 and HCG18) (Table 2).
Genes with top significance of COX univariate regression in the yellow and turquoise modules.
Classification validation of factors in yellow and turquoise modules
GSE62254, the expression profile dataset of gastric cancer, was obtained from GEO database (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE62254), including survival time information of 300 samples, which generated from the GPL570 platform. In the pretreated samples, 20,692 genes were found and their box plot distribution is shown in Figure 6.
Thirty-four genes within yellow module were gained from the GSE62254 database. Using these module genes, samples were well grouped into two classes. Moreover, the survival time in Figure 7A showed significant differences between the two groups of samples (P=0.0144). The similar results were received via 81 genes of turquoise module contained in GSE62254 database: samples were divided into three groups, and their survival times exhibited significant differences (P=0.00128; Figure 7B). These analyses demonstrated that selected genes of yellow and turquoise modules have significant correlation with survival time of GC samples.
Clustering of gastric cancer samples in GSE62254 dataset based on genes in yellow (A) and turquoise (B) modules.
Discussion
Lymph node metastasis and recurrence are the main factors that affect the prognosis of GC, and lymph nodes without metastasis have important immune monitoring functions. Therefore, it is very important to accurately determine the extent and degree of lymph node metastasis, and to carry out rational lymph node dissection. As an early and complicated event in GC metastasis, lymph node metastasis involves a series of functional and regulating genes (44. Fatourou E, Ziogas D, Baltogiannis G. Moving from lymph node metastasis in gastric cancer to biological markers. World J Surg 2010; 34: 1140–1141, doi: 10.1007/s00268-009-0368-9.
https://doi.org/10.1007/s00268-009-0368-...
), and therefore it is more meaningful to mine the co-expression network of various cancer related factors. In the present study, we downloaded large quantities of mRNA-seq and miRNA-seq data from the TCGA database to screen the mRNA, miRNA and lncRNA related to GC lymph node metastasis, and then constructed co-expression networks based on the differential expression of these factors. Compared with previous studies, which only focused on one or two factors of coding genes, miRNAs and lncRNAs (77. Ji R-C. Lymphatic endothelial cells, tumor lymphangiogenesis and metastasis: New insights into intratumoral and peritumoral lymphatics. Cancer Metast Rev 2006; 25: 677–694, doi: 10.1007/s10555-006-9026-y.
https://doi.org/10.1007/s10555-006-9026-...
–1010. Saito T, Kurashige J, Nambara S, Komatsu H, Hirata H, Ueda M, et al. A long non-coding RNA activated by transforming growth factor-β is an independent prognostic marker of gastric cancer. Ann Surg Oncol 2015; 22: 915–922, doi: 10.1245/s10434-015-4554-8.
https://doi.org/10.1245/s10434-015-4554-...
,1515. 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 :1363–1372, doi: 10.1038/onc.2012.156.
https://doi.org/10.1038/onc.2012.156...
,2525. Ying H, Lv J, Ying T, Jin S, Shao J, Wang L, et al. RETRACTED ARTICLE: Gene-gene interaction network analysis of ovarian cancer using TCGA data. J Ovarian Res 2013; 6: 1, doi: 10.1186/1757-2215-6-1.
https://doi.org/10.1186/1757-2215-6-1...
26. Yang Y, Han L, Yuan Y, Li J, Hei N, Liang H. Gene co-expression network analysis reveals common system-level properties of prognostic genes across cancer types. Nat Commun 2014; 5: 321, doi: 10.1038/ncomms4231.
https://doi.org/10.1038/ncomms4231...
–2727. Zhong L, Zhu K, Jin N, Wu D, Zhang J, Guo B, et al. A systematic analysis of miRNA-mRNA paired variations reveals widespread miRNA misregulation in breast cancer. BioMed Res Int 2014; 2014: 291280, doi: 10.1155/2014/291280.
https://doi.org/10.1155/2014/291280...
), the present study performed a systematic analysis of the three factors for the first time. There were 542 DEGs which were functionally enriched into 15 major GO terms in the biological process category, most of which were related to cancer, such as immune response, cell proliferation, cell migration, cell differentiation and cell adhesion.
Thanks to the rapid development in bioinformatics and statistical genomics, gene network analysis has become a powerful approach that can explore the interactions between genes and has been widely applied in gene expression studies of humans and model organisms (1717. Lee I, Date SV, Adai AT, Marcotte EM. A probabilistic functional network of yeast genes. Science 2004; 306: 1555–1558, doi: 10.1126/science.1099511.
https://doi.org/10.1126/science.1099511...
–1919. Ponomarev I, Wang S, Zhang L, Harris RA, Mayfield RD. Gene coexpression networks in human brain identify epigenetic modifications in alcohol dependence. J Neurosc 2012; 32: 1884–1897, doi: 10.1523/JNEUROSCI.3136-11.2012.
https://doi.org/10.1523/JNEUROSCI.3136-1...
). Meanwhile, genes that are highly interconnected within the network are usually involved in the same biological modules or pathways, and therefore, modular analysis also plays an important role in the analysis of gene co-expression network (2828. Stuart JM, Segal E, Koller D, Kim SK. A gene-coexpression network for global discovery of conserved genetic modules. Science 2003; 302: 249–255, doi: 10.1126/science.1087447.
https://doi.org/10.1126/science.1087447...
,2929. Han J-DJ, Bertin N, Hao T, Goldberg DS, Berriz GF, Zhang LV, et al. Evidence for dynamically organized modularity in the yeast protein-protein interaction network. Nature 2004; 430: 88–93, doi: 10.1038/nature02555.
https://doi.org/10.1038/nature02555...
). Several studies have demonstrated the value of analyzing networks based on TCGA database. In the present study, 908 differentially expressed factors including 542 genes, 42 miRNAs and 324 lncRNAs were included in the network analysis, and furthermore, genes in the co-expression networks were used for the modular mining. Four module groups were coded in blue, brown, turquoise and yellow color, which included 73, 39, 89, and 35 differentially expressed factors. Except genes in the group 3 of the yellow module lacking sufficient death number, the other genes in yellow and turquoise modules could separate samples efficiently, and there were significant differences between the survival curves of isolated samples.
Based on the network and modular analysis, seven (ADORA3, TLR7, IRF4, CCR4, RTN1, GAPT, and GRAP2) and ten (GNAO1, ISM1, CILP, SLIT2, SCRG1, TNFAIP8L3, hsa-mir-183, sa-mir-942, MIR345 and HCG18) candidate factors with top significance of COX univariate regression were identified in the yellow and turquoise module, respectively. Within yellow module, ADORA3 (3030. Tsuchiya A, Nishizaki T. Anticancer effect of adenosine on gastric cancer via diverse signaling pathways. World J Gastroenterol 2015; 21: 10931, doi: 10.3748/wjg.v21.i39.10931.
https://doi.org/10.3748/wjg.v21.i39.1093...
), TLR7 (3131. Jiang J, Dong L, Qin B, Guo X, Li H, Shi H, et al. [Expression of Toll-like receptor 7 in gastric cancer cell lines and effects of TLR7 agonist on proliferation and apoptosis of SGC-7901 cells in vitro]. Nan fang yi ke da xue xue bao 2014; 34: 1606–1610.), IR4 (3232. Yamashita M, Toyota M, Suzuki H, Nojima M, Yamamoto E, Kamimae S, et al. DNA methylation of interferon regulatory factors in gastric cancer and noncancerous gastric mucosae. Cancer Sci 2010; 101: 1708–1716, doi: 10.1111/j.1349-7006.2010.01581.x.
https://doi.org/10.1111/j.1349-7006.2010...
), and CCR4 (3333. Cao L, Hu X, Zhang J, Huang G, Zhang Y. The role of the CCL22-CCR4 axis in the metastasis of gastric cancer cells into omental milky spots. J Translat Med 2014; 12: 267, doi: 10.1186/s12967-014-0267-1.
https://doi.org/10.1186/s12967-014-0267-...
34. Chang W-J, Du Y, Zhao X, Ma L-Y, Cao G-W. Inflammation-related factors predicting prognosis of gastric cancer. World J Gastroenterol 2014; 20: 4586–4596, doi: 10.3748/wjg.v20.i16.4586.
https://doi.org/10.3748/wjg.v20.i16.4586...
–3535. Yang Y, Du L, Yang X, Qu A, Zhang X, Zhou C, et al. Aberrant CCR4 expression is involved in tumor invasion of lymph node-negative human gastric cancer. PloS One 2015; 10: e0120059, doi: 10.1371/journal.pone.0120059.
https://doi.org/10.1371/journal.pone.012...
) are genes closely related to lymph node metastasis and survival of GC, and within turquoise module, two miRNAs (miRNA-183 and miRNA-942) (1313. Liu N, Zuo C, Wang X, Chen T, Yang D, Wang J, et al. miR-942 decreases TRAIL-induced apoptosis through ISG12a downregulation and is regulated by AKT. Oncotarget 2014; 5: 4959, doi: 10.18632/oncotarget.2067.
https://doi.org/10.18632/oncotarget.2067...
,3636. Gu W, Gao T, Shen J, Sun Y, Zheng X, Wang J, et al. MicroRNA-183 inhibits apoptosis and promotes proliferation and invasion of gastric cancer cells by targeting PDCD4. Int J Clin Exper Med 2014; 7: 2519–2529.) and three genes (GNAO1, ISM1and SLIT2) are factors involved in the regulation of GC. The above published studies prove that our gene network analysis for screening candidate factors related to GC for further evaluation was reliable.
Therefore, the remaining seven identified factors including three genes in yellow module (RTN1, GAPT, GRAP2) and five factors in turquoise module (CILP, SCRG1, TNFAIP8L3, MIR345 and HCG18) could be new factors related to survival of GC. In fact, there is evidence that these genes may be associated with several diseases and even cancer. For example, RTN1, a neuroendocrine cell specific protein, localized in endoplasmic reticulum, might be involved in the activation of the expression of androgen-responsive genes and related to prostate cancer (3737. Levina E, Ji H, Chen M, Baig M, Oliver D, Ohouo P, et al. Identification of novel genes that regulate androgen receptor signaling and growth of androgen-deprived prostate cancer cells. Oncotarget 2015; 6: 13088, doi: 10.18632/oncotarget.3743.
https://doi.org/10.18632/oncotarget.3743...
). As an adapter protein, GRB2 has been identified as a major mediator in Ras-mitogen-activated protein kinase (MAPK) activation, which is essential for growth factor-induced cell proliferation and differentiation and plays a central role in embryo development and malignant transformation. Therefore, we believe that GAPT and GRAP2 are involved in the activation of MAPK and growth factor-induced cell proliferation and differentiation, which are often associated with the development of cancer (3838. Evan GI, Vousden KH. Proliferation, cell cycle and apoptosis in cancer. Nature 2001; 411: 342–348, doi: 10.1038/35077213.
https://doi.org/10.1038/35077213...
). CILP, an extracellular matrix protein abundant in cartilaginous tissues, is implicated in common musculoskeletal disorders, including osteoarthritis and lumbar disc disease (3939. Mori M, Nakajima M, Mikami Y, Seki S, Takigawa M, Kubo T, et al. Transcriptional regulation of the cartilage intermediate layer protein (CILP) gene. Biochem Biophys Res Commun 2006; 341: 121–127, doi: 10.1016/j.bbrc.2005.12.159.
https://doi.org/10.1016/j.bbrc.2005.12.1...
). It is worth noting that the TNFAIP8 family members are usually associated with immune homeostasis and inflammatory cancer diseases. For example, TNFAIP8 itself usually functions as an oncogenic molecule and it is also associated with enhanced cell survival and inhibition of apoptosis, and TNFAIP8-like 2 (TIPE2) governs immune homeostasis in both the innate and adaptive immune system and prevents hyper-responsiveness (4040. Luan Y, Yao Y, Sheng Z. The tumor necrosis factor-alpha-induced protein 8 family in immune homeostasis and inflammatory cancer diseases. J Biol Regul Homeost Agents 2013; 27: 611–619.). However, the function of TNFAIP8L3 remains unclear. Therefore, our study provides some insight into the emergent properties of prognostic genes, and further investigation of the functional roles of these newly identified factors is urgent for a functional validation system in GC.
In conclusion, our data provides a comprehensive bioinformatics analysis of genes and pathways, which may be involved in the lymph node metastasis of GC. We found a total of 542 genes, 42 miRNAs and 324 lncRNAs, and then constructed the interaction networks of these differentially expressed factors. Furthermore, we conducted functional modules analysis in the network, and found that except for the genes in group 3 of the yellow module, the other genes in the yellow and turquoise modules could separate samples, and therefore these genes could be potential prognostic biomarkers. However, further analyses are still required to reveal their mechanism in the process of tumor genesis and development in GC.
Supplementary Material
Click here to view [pdf].
References
-
1Torre LA, Bray F, Siegel RL, Ferlay J, Lortet-Tieulent J, Jemal A. Global cancer statistics. CA Cancer J Clin 2015; 65: 87–108, doi: 10.3322/caac.21262.
» https://doi.org/10.3322/caac.21262 -
2Xue L, Lü N, Li A, Zou S, Lin D, He Z, et al. A clinicopathological analysis of gastric lymphoma. Chinese J Pathol 2005; 34: 332–336.
-
3Roviello F, Rossi S, Marrelli D, Pedrazzani C, Corso G, Vindigni C, et al. Number of lymph node metastases and its prognostic significance in early gastric cancer: a multicenter Italian study. J Surg Oncol 2006; 94: 275–280, doi: 10.1002/jso.20566.
» https://doi.org/10.1002/jso.20566 -
4Fatourou E, Ziogas D, Baltogiannis G. Moving from lymph node metastasis in gastric cancer to biological markers. World J Surg 2010; 34: 1140–1141, doi: 10.1007/s00268-009-0368-9.
» https://doi.org/10.1007/s00268-009-0368-9 -
5Chan A, Wong B, Lam S-K. Gastric cancer: past, present and future. Can J Gastroenterol 2001; 15: 469–474, doi: 10.1155/2001/850308.
» https://doi.org/10.1155/2001/850308 -
6Nishida T. Treatment of gastric cancer in Asia: the missing link. Lancet Oncol 2009; 10: 1027–1028, doi: 10.1016/S1470-2045(09)70322-5.
» https://doi.org/10.1016/S1470-2045(09)70322-5 -
7Ji R-C. Lymphatic endothelial cells, tumor lymphangiogenesis and metastasis: New insights into intratumoral and peritumoral lymphatics. Cancer Metast Rev 2006; 25: 677–694, doi: 10.1007/s10555-006-9026-y.
» https://doi.org/10.1007/s10555-006-9026-y -
8Shao Y, Chen H, Jiang X, Chen S, Li P, Ye M, et al. Low expression of lncRNA-HMlincRNA717 in human gastric cancer and its clinical significances. Tumor Biol 2014; 35: 9591–9595, doi: 10.1007/s13277-014-2243-z.
» https://doi.org/10.1007/s13277-014-2243-z -
9Wang Y, Zhang D, Wu K, Zhao Q, Nie Y, Fan D. Long noncoding RNA MRUL promotes ABCB1 expression in multidrug-resistant gastric cancer cell sublines. Mol Cell Bio 2014; 34: 3182–3193, doi: 10.1128/MCB.01580-13.
» https://doi.org/10.1128/MCB.01580-13 -
10Saito T, Kurashige J, Nambara S, Komatsu H, Hirata H, Ueda M, et al. A long non-coding RNA activated by transforming growth factor-β is an independent prognostic marker of gastric cancer. Ann Surg Oncol 2015; 22: 915–922, doi: 10.1245/s10434-015-4554-8.
» https://doi.org/10.1245/s10434-015-4554-8 -
11Ueda T, Volinia S, Okumura H, Shimizu M, Taccioli C, Rossi S, et al. Relation between microRNA expression and progression and prognosis of gastric cancer: a microRNA expression analysis. Lancet Oncol 2010; 11: 136–146, doi: 10.1016/S1470-2045(09)70343-2.
» https://doi.org/10.1016/S1470-2045(09)70343-2 -
12Liu X-h, Sun M, Nie F-q, Ge Y-b, Zhang E-b, Yin D-d, et al. Lnc RNA HOTAIR functions as a competing endogenous RNA to regulate HER2 expression by sponging miR-331-3p in gastric cancer. Mol Cancer 2014; 13: 1, doi: 10.1186/1476-4598-13-1.
» https://doi.org/10.1186/1476-4598-13-1 -
13Liu N, Zuo C, Wang X, Chen T, Yang D, Wang J, et al. miR-942 decreases TRAIL-induced apoptosis through ISG12a downregulation and is regulated by AKT. Oncotarget 2014; 5: 4959, doi: 10.18632/oncotarget.2067.
» https://doi.org/10.18632/oncotarget.2067 -
14Tsukamoto Y, Nakada C, Noguchi T, Tanigawa M, Nguyen LT, Uchida T, et al. MicroRNA-375 is downregulated in gastric carcinomas and regulates cell survival by targeting PDK1 and 14-3-3ζ. Cancer Res 2010; 70: 2339–2349, doi: 10.1158/0008-5472.CAN-09-2777.
» https://doi.org/10.1158/0008-5472.CAN-09-2777 -
15Zhao 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 :1363–1372, doi: 10.1038/onc.2012.156.
» https://doi.org/10.1038/onc.2012.156 -
16Sakamoto N, Naito Y, Oue N, Sentani K, Uraoka N, Zarni Oo H, et al. MicroRNA‐148a is downregulated in gastric cancer, targets MMP7, and indicates tumor invasiveness and poor prognosis. Cancer Sci 2014; 105: 236–243, doi: 10.1111/cas.12330.
» https://doi.org/10.1111/cas.12330 -
17Lee I, Date SV, Adai AT, Marcotte EM. A probabilistic functional network of yeast genes. Science 2004; 306: 1555–1558, doi: 10.1126/science.1099511.
» https://doi.org/10.1126/science.1099511 -
18Oldham MC, Konopka G, Iwamoto K, Langfelder P, Kato T, Horvath S, et al. Functional organization of the transcriptome in human brain. Nature Neurosci 2008; 11: 1271–1282, doi: 10.1038/nn.2207.
» https://doi.org/10.1038/nn.2207 -
19Ponomarev I, Wang S, Zhang L, Harris RA, Mayfield RD. Gene coexpression networks in human brain identify epigenetic modifications in alcohol dependence. J Neurosc 2012; 32: 1884–1897, doi: 10.1523/JNEUROSCI.3136-11.2012.
» https://doi.org/10.1523/JNEUROSCI.3136-11.2012 -
20Larsson O, Wahlestedt C, Timmons JA. Considerations when using the significance analysis of microarrays (SAM) algorithm. BMC Bioinformatics 2005; 6: 129, doi: 10.1186/1471-2105-6-129.
» https://doi.org/10.1186/1471-2105-6-129 -
21Dennis G Jr, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, et al. DAVID: database for annotation, visualization, and integrated discovery. Genome Biol 2003; 4: P3, doi: 10.1186/gb-2003-4-5-p3.
» https://doi.org/10.1186/gb-2003-4-5-p3 -
22Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics 2008; 9: 1, doi: 10.1186/1471-2105-9-559.
» https://doi.org/10.1186/1471-2105-9-559 -
23Jiang D, Wang H, Li J, Wu Y, Fang M, Yang R. Cox regression model for dissecting genetic architecture of survival time. Genomics 2014; 104: 472–476, doi: 10.1016/j.ygeno.2014.10.002.
» https://doi.org/10.1016/j.ygeno.2014.10.002 -
24Mossa S, Barthelemy M, Stanley HE, Amaral LAN. Truncation of power law behavior in “scale-free” network models due to information filtering. Phys Rev Letters 2002; 88: 138701, doi: 10.1103/PhysRevLett.88.138701.
» https://doi.org/10.1103/PhysRevLett.88.138701 -
25Ying H, Lv J, Ying T, Jin S, Shao J, Wang L, et al. RETRACTED ARTICLE: Gene-gene interaction network analysis of ovarian cancer using TCGA data. J Ovarian Res 2013; 6: 1, doi: 10.1186/1757-2215-6-1.
» https://doi.org/10.1186/1757-2215-6-1 -
26Yang Y, Han L, Yuan Y, Li J, Hei N, Liang H. Gene co-expression network analysis reveals common system-level properties of prognostic genes across cancer types. Nat Commun 2014; 5: 321, doi: 10.1038/ncomms4231.
» https://doi.org/10.1038/ncomms4231 -
27Zhong L, Zhu K, Jin N, Wu D, Zhang J, Guo B, et al. A systematic analysis of miRNA-mRNA paired variations reveals widespread miRNA misregulation in breast cancer. BioMed Res Int 2014; 2014: 291280, doi: 10.1155/2014/291280.
» https://doi.org/10.1155/2014/291280 -
28Stuart JM, Segal E, Koller D, Kim SK. A gene-coexpression network for global discovery of conserved genetic modules. Science 2003; 302: 249–255, doi: 10.1126/science.1087447.
» https://doi.org/10.1126/science.1087447 -
29Han J-DJ, Bertin N, Hao T, Goldberg DS, Berriz GF, Zhang LV, et al. Evidence for dynamically organized modularity in the yeast protein-protein interaction network. Nature 2004; 430: 88–93, doi: 10.1038/nature02555.
» https://doi.org/10.1038/nature02555 -
30Tsuchiya A, Nishizaki T. Anticancer effect of adenosine on gastric cancer via diverse signaling pathways. World J Gastroenterol 2015; 21: 10931, doi: 10.3748/wjg.v21.i39.10931.
» https://doi.org/10.3748/wjg.v21.i39.10931 -
31Jiang J, Dong L, Qin B, Guo X, Li H, Shi H, et al. [Expression of Toll-like receptor 7 in gastric cancer cell lines and effects of TLR7 agonist on proliferation and apoptosis of SGC-7901 cells in vitro]. Nan fang yi ke da xue xue bao 2014; 34: 1606–1610.
-
32Yamashita M, Toyota M, Suzuki H, Nojima M, Yamamoto E, Kamimae S, et al. DNA methylation of interferon regulatory factors in gastric cancer and noncancerous gastric mucosae. Cancer Sci 2010; 101: 1708–1716, doi: 10.1111/j.1349-7006.2010.01581.x.
» https://doi.org/10.1111/j.1349-7006.2010.01581.x -
33Cao L, Hu X, Zhang J, Huang G, Zhang Y. The role of the CCL22-CCR4 axis in the metastasis of gastric cancer cells into omental milky spots. J Translat Med 2014; 12: 267, doi: 10.1186/s12967-014-0267-1.
» https://doi.org/10.1186/s12967-014-0267-1 -
34Chang W-J, Du Y, Zhao X, Ma L-Y, Cao G-W. Inflammation-related factors predicting prognosis of gastric cancer. World J Gastroenterol 2014; 20: 4586–4596, doi: 10.3748/wjg.v20.i16.4586.
» https://doi.org/10.3748/wjg.v20.i16.4586 -
35Yang Y, Du L, Yang X, Qu A, Zhang X, Zhou C, et al. Aberrant CCR4 expression is involved in tumor invasion of lymph node-negative human gastric cancer. PloS One 2015; 10: e0120059, doi: 10.1371/journal.pone.0120059.
» https://doi.org/10.1371/journal.pone.0120059 -
36Gu W, Gao T, Shen J, Sun Y, Zheng X, Wang J, et al. MicroRNA-183 inhibits apoptosis and promotes proliferation and invasion of gastric cancer cells by targeting PDCD4. Int J Clin Exper Med 2014; 7: 2519–2529.
-
37Levina E, Ji H, Chen M, Baig M, Oliver D, Ohouo P, et al. Identification of novel genes that regulate androgen receptor signaling and growth of androgen-deprived prostate cancer cells. Oncotarget 2015; 6: 13088, doi: 10.18632/oncotarget.3743.
» https://doi.org/10.18632/oncotarget.3743 -
38Evan GI, Vousden KH. Proliferation, cell cycle and apoptosis in cancer. Nature 2001; 411: 342–348, doi: 10.1038/35077213.
» https://doi.org/10.1038/35077213 -
39Mori M, Nakajima M, Mikami Y, Seki S, Takigawa M, Kubo T, et al. Transcriptional regulation of the cartilage intermediate layer protein (CILP) gene. Biochem Biophys Res Commun 2006; 341: 121–127, doi: 10.1016/j.bbrc.2005.12.159.
» https://doi.org/10.1016/j.bbrc.2005.12.159 -
40Luan Y, Yao Y, Sheng Z. The tumor necrosis factor-alpha-induced protein 8 family in immune homeostasis and inflammatory cancer diseases. J Biol Regul Homeost Agents 2013; 27: 611–619.
Publication Dates
-
Publication in this collection
2018
History
-
Received
15 July 2017 -
Accepted
17 Nov 2017