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.

Gene Coexpression Network Characterizing Microenvironmental Heterogeneity and Intercellular Communication in Pancreatic Ductal Adenocarcinoma: Implications of Prognostic Significance and Therapeutic Target

Analysis of Dynamic Molecular Networks for Pancreatic Ductal Adenocarcinoma Progression

Bioinformatics-Based Identification of Tumor Microenvironment-Related Prognostic Genes in Pancreatic Cancer

Single-cell Transcriptome Reveals the Heterogeneity of Malignant Ductal Cells and the Prognostic Value of REG4 and SPINK1 in Primary Pancreatic Ductal Adenocarcinoma

Identification of Key Prognostic Biomarker and Its Correlation with Immune Infiltrates in Pancreatic Ductal Adenocarcinoma

Identification of prognosis-related genes and construction of multi-regulatory networks in pancreatic cancer microenvironment by bioinformatics analysis

Gene Expression Characteristics of Tumor and Adjacent Non-Tumor Tissues of Pancreatic Ductal Adenocarcinoma (PDAC) In-Silico

Leveraging Tumor Microenvironment Infiltration in Pancreatic Cancer to Identify Gene Signatures Related to Prognosis and Immunotherapy Response

Tumor Microenvironment Characteristics of Pancreatic Cancer to Determine Prognosis and Immune-Related Gene Signatures

Deciphering The Prognostic Implications Of The Components And Signatures In The Immune Microenvironment Of Pancreatic Ductal Adenocarcinoma

Gene Coexpression Network Approach to Develop an Immune Prognostic Model for Pancreatic Adenocarcinoma

Genetic Landscape of Prognostic Value in Pancreas Ductal Adenocarcinoma Microenvironment-Reply.

Clinical Functional Proteomics of Intercellular Signalling in Pancreatic Cancer

Identifying Drug Targets in Pancreatic Ductal Adenocarcinoma Through Machine Learning, Analyzing Biomolecular Networks, and Structural Modeling.

Immune Profiling and Prognostic Model of Pancreatic Cancer Using Quantitative Pathology and Single-Cell RNA Sequencing.

Integrated Expression Profiles Analysis Reveals Novel Predictive Biomarker in Pancreatic Ductal Adenocarcinoma

Prognostic Utility of TME-associated Genes in Pancreatic Cancer.

Comprehensive analyses and experimental verification of NETs and an EMT gene signature for prognostic prediction, immunotherapy, and chemotherapy in pancreatic adenocarcinoma

Elevated Expression of Cell Adhesion, Metabolic, and Mucus Secretion Gene Clusters Associated with Tumorigenesis, Metastasis, and Poor Survival in Pancreatic Ductal Adenocarcinoma

Integrating Tertiary Lymphoid Structure-Associated Genes into Computational Models to Evaluate Prognostication and Immune Infiltration in Pancreatic Cancer

Analysis of Immune-Related Signatures Related to CD4+ T Cell Infiltration with Gene Co-Expression Network in Pancreatic Adenocarcinoma.