Identifying functional evolution processes according to the pathological stages of colorectal cancer Colorectal cancer (CRC) is one of the malignant tumors with high morbidity and mortality. A prevalent method for studying colorectal cancer is to identify differentially expressed genes (DEGs) between control and patient samples, followed by the pathway enrichment analyses. However, many of those studies ignore the fact that different pathological stages of the cancer are often highly different from each other. The mixture of those heterogeneous samples may lack the efficiency of identifying the real DEGs, and loss the opportunity to analyze the dynamic evolution process of cancer. In this study, we develop a feasible framework to identify function evolution processes of cancers according to their pathological stages. Firstly, the limma package was used to identify DEG sets between control and CRC stage I, II, III, and IV samples, separately. Secondly, a pathway interaction network was constructed by taking a comprehensive analysis of pathways at individual stages, and a functional module interaction network was also generated independently by clustering genes into modules in a PPI network. The relationship between pathways and modules in adjacent stages was analyzed for all stages of CRC. A total of 479, 313, 349, and 383 DEGs were identified and they were enriched in 17, 16, 20, and 24 pathways, respectively. A functional evolution network was constructed by using those modules, and two significant evolution processes a2b1c2d1 (Mod1) and a1b2c1d2 (Mod2) were identified which may play critical roles in the development of CRC. The framework proposed in this study can be used to explore molecular mechanisms and evolution processes of CRC.

Identifying Structural Changes in Correlation Networks Models of Cancer Gene Expression by Stage Gene expression analysis using correlation network modeling can help to identify systemslevel cellular changes and cooperation among genes. Network modeling is a relatively novel method for comparing changes across stages of cancer, or between primary tumor tissue and the metastatic tissue. In this study, we develop a pipeline to identify the dynamic changes of cancer in gene expression level through timedependent network analysis of cancer data from The Cancer Genome Atlas (TCGA). A total of 16 correlation networks were built from four (4) stages in four (4) different types of cancers: Thyroid Carcinoma, Colon Adenocarcinoma and Rectum Adenocarcinoma, Stomach Adenocarcinoma, and Kidney Renal Clear Cell Carcinoma. To identify the basic changes in network structure, we performed Jaccard similarity comparison of structurally relevant nodes. We employed mutation analysis to measure and present the timebased changes in mutation rate of genes that are specific to each cancer type. Finally, we present a case study to identify the gene expression changes among primary tumor tissue and metastatic tissue in skin cutaneous melanoma (SKCM).