Screening of hub genes and pathways in colorectal cancer with microarray technology

Here we intend to identify key genes and pathways in the pathogenesis of colorectal cancer (CRC) through analyzing microarray data with bioinformatic tools. The gene expression profile dataset GSE23878 was downloaded from Gene Expression Omnibus and differentially expressed genes (DEGs) were screened out using Student's t-test. GO function and KEGG pathway enrichment analyses were performed for these DEGs with the DAVID online tool. Interaction network was constructed among the over-represented pathways based on the protein-protein interactions within the pathways. Besides, the protein interaction information obtained from HPRD database were applied to constructed protein-protein interaction networks among the DEGs and hub genes and function module were screened out. A total of 2,296 DEGs were obtained and they were enriched in 34 pathways. An interaction network was constructed among 32 pathways, in which p53 signaling pathway acted as the hub pathway as it showed the highest node degree. The protein-protein interaction network comprised 1,481 interaction relationships among 332 genes which included 40 DEGs. Further analysis revealed that theses DEGs formed 7 function modules and many genes, such as PDGFRB, MET, FZD2, CCND1, PRKCB, ARHGEF6, JUP, WNT2, WNT5A and WNT11 were key genes in the networks. The DEGs and disturbed biological functions uncovered in present study may play important roles in the development of CRC and can contribute to the understanding on molecular mechanisms of CRC. Further these DEGs we obtained can be acted as potential biomarkers for diagnosis and therapy of CRC.