Identifying key stage-specific genes and transcription factors for gastric cancer based on RNA-sequencing data

To identify gastric cancer (GC)-associated genes and transcription factors (TFs) using RNA-sequencing (RNA-seq) data of Asians.

The RNA-seq data (GSE36968) were downloaded from Gene Expression Omnibus database, including 6 noncancerous gastric tissue samples, 5 stage I GC samples, 5 stage II GC samples, 8 stage III GC samples, and 6 stage IV GC samples. The gene expression values in each sample were calculated using Cuffdiff. Following, stage-specific genes were identified by 1-way analysis of variance and hierarchical clustering analysis. Upstream TFs were identified using Seqpos. Besides, functional enrichment analysis of stage-specific genes was performed by DAVID. In addition, the underlying protein-protein interactions (PPIs) information among stage IV-specific genes were extracted from STRING database and PPI network was constructed using Cytoscape software.

A total of 3576 stage-specific genes were identified, including 813 specifically up-regulated genes in the normal gastric tissues, 2224 stage I and II-specific genes, and 539 stage IV-specific genes. Also, a total of 9 and 11 up-regulated TFs were identified for the stage I and II-specific genes and stage IV-specific genes, respectively. Functional enrichment showed SPARC, MMP17, and COL6A3 were related to extracellular matrix. Notably, 2 regulatory pathways HOXA4-GLI3-RUNX2-FGF2 and HMGA2-PRKCA were obtained from the PPI network for stage IV-specific genes. In the PPI network, TFs HOXA4 and HMGA2 might function via mediating other genes.

These stage-specific genes and TFs might act in the pathogenesis of GC in Asians.