Abstract:BackgroundPancreatic ductal adenocarcinoma (PDAC) is one of the deadliest solid tumors. The rapid progression of PDAC results in an advanced stage of patients when diagnosed. However, the dynamic molecular mechanism underlying PDAC progression remains far from clear.MethodsThe microarray GSE62165 containing PDAC staging samples was obtained from Gene Expression Omnibus and the differentially expressed genes (DEGs) between normal tissue and PDAC of different stages were profiled using R software, respectively. The software program Short Time-series Expression Miner was applied to cluster, compare, and visualize gene expression differences between PDAC stages. Then, function annotation and pathway enrichment of DEGs were conducted by Database for Annotation Visualization and Integrated Discovery. Further, the Cytoscape plugin DyNetViewer was applied to construct the dynamic protein-protein interaction networks and to analyze different topological variation of nodes and clusters over time. The phosphosite markers of stage-specific protein kinases were predicted by PhosphoSitePlus database. Moreover, survival analysis of candidate genes and pathways was performed by Kaplan-Meier plotter. Finally, candidate genes were validated by immunohistochemistry in PDAC tissues.ResultsCompared with normal tissues, the total DEGs number for each PDAC stage were 994 (stage I), 967 (stage IIa), 965 (stage IIb), 1027 (stage III), 925 (stage IV), respectively. The stage-course gene expression analysis showed that 30 distinct expressional models were clustered. Kyoto Encyclopedia of Genes and Genomes analysis indicated that the up-regulated DEGs were commonly enriched in five fundamental pathways throughout five stages, including pathways in cancer, small cell lung cancer, ECM-receptor interaction, amoebiasis, focal adhesion. Except for amoebiasis, these pathways were associated with poor PDAC overall survival. Meanwhile, LAMA3, LAMB3, LAMC2, COL4A1 and FN1 were commonly shared by these five pathways and were unfavorable factors for prognosis. Furthermore, by constructing the stage-course dynamic protein interaction network, 45 functional molecular modules and 19 nodes were identified as featured regulators for all PDAC stages, among which the collagen family and integrins were considered as two main regulators for facilitating aggressive progression. Additionally, the clinical relevance analysis suggested that the stage IV featured nodes MLF1IP and ITGB4 were significantly correlated with shorter overall survival. Moreover, 15 stage-specific protein kinases were identified from the dynamic network and CHEK1 was particularly activated at stage IV. Experimental validation showed that MLF1IP, LAMA3 and LAMB3 were progressively increased from tumor initiation to progression.ConclusionsOur study provided a view for a better understanding of the dynamic landscape of molecular interaction networks during PDAC progression and offered potential targets for therapeutic intervention.

Predicting Stage-Specific Cancer Related Genes and Their Dynamic Modules by Integrating Multiple Datasets

Integrating Multiple Datasets to Discover Stage-Specific Cancer Related Genes and Stage-Specific Pathways.

Extracting Stage-Specific and Dynamic Modules Through Analyzing Multiple Networks Associated with Cancer Progression

Prediction of Disease Genes Based on Stage-Specific Gene Regulatory Networks in Breast Cancer

Identifying Disease Genes and Module Biomarkers by Differential Interactions

Identifying Disease Feature Genes Based on Cellular Localized Gene Functional Modules and Regulation Networks

Peeling Off the Hidden Genetic Heterogeneities of Cancers Based on Disease-Relevant Functional Modules

Analysis of Dynamic Molecular Networks for Pancreatic Ductal Adenocarcinoma Progression

Identification of breast cancer risk modules via an integrated strategy

An Integrative Framework of Heterogeneous Genomic Data for Cancer Dynamic Modules Based on Matrix Decomposition.

The Detection of Gene Modules with Overlapping Characteristic via Integrating Multi-omics Data in Six Cancers

Identifying Cancer Prognostic Modules by Module Network Analysis

Multi-omics with Dynamic Network Biomarker Algorithm Prefigures Organ-Specific Metastasis of Lung Adenocarcinoma

Identifying Breast Cancer-Related Genes Based on a Novel Computational Framework Involving KEGG Pathways and PPI Network Modularity

A Cancer Gene Module Mining Method Based on Bio-Network of Multi-Omics Gene Groups

Identification of cancer prognosis-associated functional modules using differential co-expression networks.

An integrated approach to identify causal network modules of complex diseases with application to colorectal cancer.

Predicting Stage-Specific Recurrent Aberrations from Somatic Copy Number Dataset

A Network-Based Approach to Identify Disease-Associated Gene Modules Through Integrating DNA Methylation and Gene Expression.

Integrative Analysis Based on Survival Associated Co-Expression Gene Modules for Predicting Neuroblastoma Patients' Survival Time

Functional modules with disease discrimination abilities for various cancers