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.

Dynamic Network of Transcription and Pathway Crosstalk to Reveal Molecular Mechanism of MGd-treated Human Lung Cancer Cells.

Identification of Key Genes and Evaluation of Clinical Outcomes in Lung Squamous Cell Carcinoma Using Integrated Bioinformatics Analysis

Molecular Mechanism of Gd@C82(OH)22 Increasing Collagen Expression: Implication for Encaging Tumor

Analysis of Dynamic Molecular Networks for Pancreatic Ductal Adenocarcinoma Progression

Integrative Omics Analysis Reveals Metabolic Features of Ground-Glass Opacity-Associated Lung Cancer

Identification Of Novel Biomarkers And Candidate Small Molecule Drugs In Non-Small-Cell Lung Cancer By Integrated Microarray Analysis

Comprehensive Analyses of M1a Regulator-Mediated Modification Patterns Determining Prognosis in Lower-Grade Glioma (running Title: M1a in LGG)

Identification of Key Genes with Differential Correlations in Lung Adenocarcinoma

Identification Of A Cerna Network In Lung Adenocarcinoma Based On Integration Analysis Of Tumor-Associated Macrophage Signature Genes

Identification of Target Gene and Prognostic Evaluation for Lung Adenocarcinoma Using Gene Expression Meta-Analysis, Network Analysis and Neural Network Algorithms

Characterizing dynamic regulatory programs in mouse lung development and their potential association with tumourigenesis via miRNA-TF-mRNA circuits

Identification of Biomarker and Co-Regulatory Motifs in Lung Adenocarcinoma Based on Differential Interactions

Discovery of a Set of Biomarkers of Human Lung Adenocarcinoma Through Cell-Map Proteomics and Bioinformatics

Identification of Key Biomarkers and Candidate Molecules in Non-Small-Cell Lung Cancer by Integrated Bioinformatics Analysis

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

Functional Analysis of Microrna and Transcription Factor Synergistic Regulatory Network Based on Identifying Regulatory Motifs in Non-Small Cell Lung Cancer

Exploring the Mechanism of Periplogenin in the Treatment of Lung Cancer Based on Network Pharmacology and Molecular Biology

Varied Pathways of Stage Ia Lung Adenocarcinomas Discovered by Integrated Gene Expression Analysis

Uncovering the dynamic effects of DEX treatment on lung cancer by integrating bioinformatic inference and multiscale modeling of scRNA-seq and proteomics data

A Transcriptional Mirna-Gene Network Associated with Lung Adenocarcinoma Metastasis Based on the TCGA Database.

Multidirectional characterization of cellular composition and spatial architecture in human multiple primary lung cancers