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.

A study on the role of Taxifolin in inducing apoptosis of pancreatic cancer cells: screening results using weighted gene co-expression network analysis

Integrated Bioinformatics Analysis Identified COL11A1 As an Immune Infiltrates Correlated Prognosticator in Pancreatic Adenocarcinoma

Analysis of Dynamic Molecular Networks for Pancreatic Ductal Adenocarcinoma Progression

Down-expression of CD36 in Pancreatic Adenocarcinoma and Its Correlation with Clinicopathological Features and Prognosis.

Promotion of apoptosis in melanoma cells by taxifolin through the PI3K/AKT signaling pathway: Screening of natural products using WGCNA and CMAP platforms

Identification of Survival-Associated Hub Genes in Pancreatic Adenocarcinoma Based on WGCNA

Comprehensive analysis and experimental verification of the mechanism of anoikis related genes in pancreatic cancer

Integrating Genetic and Transcriptomic Data to Reveal Pathogenesis and Prognostic Markers of Pancreatic Adenocarcinoma

Construction of a Competitive Endogenous RNA Network for Pancreatic Adenocarcinoma Based on Weighted Gene Co-expression Network Analysis and a Prognosis Model

Molecular characterization analysis of PANoptosis-related genes in colorectal cancer based on bioinformatic analysis

A feasibility study on gene therapy of pancreatic carcinoma with Ad-PUMA.

Identification of characteristic markers correlated with Th2 cell infiltration and metabolism molecular subtype in pancreatic adenocarcinoma

Cuproptosis-Related Gene DLAT as a Novel Biomarker Correlated with Prognosis, Chemoresistance, and Immune Infiltration in Pancreatic Adenocarcinoma: A Preliminary Study Based on Bioinformatics Analysis

Identification of genes highly downregulated in pancreatic cancer through a meta-analysis of microarray datasets: implications for discovery of novel tumor-suppressor genes and therapeutic targets

Integrative Analysis Reveals Clinically Relevant Molecular Fingerprints in Pancreatic Cancer

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

Identification of fibroblast-related genes based on single-cell and machine learning to predict the prognosis and endocrine metabolism of pancreatic cancer

Integrated analysis identifies cuproptosis-related gene DLAT and its competing endogenous RNAs network to predict the prognosis of pancreatic adenocarcinoma patients

Identification of potential crucial genes associated with the pathogenesis and prognosis of pancreatic adenocarcinoma

A Novel Cuproptosis-Associated Gene Signature to Predict Prognosis in Patients with Pancreatic Cancer

Molecular Landscape and Prognostic Biomarker Analysis of Advanced Pancreatic Cancer and Predictors of Treatment Efficacy of AG Chemotherapy