scholarly journals Identification and Interaction Analysis of Molecular Markers in Pancreatic Ductal Adenocarcinoma by Integrated Bioinformatics Analysis and Molecular Docking Experiments

Author(s):  
Basavaraj Vastrad ◽  
Chanabasayya Vastrad ◽  
Anandkumar Tengli

AbstractThe current investigation aimed to mine therapeutic molecular targets that play an key part in the advancement of pancreatic ductal adenocarcinoma (PDAC). The expression profiling by high throughput sequencing dataset profile GSE133684 dataset was downloaded from the Gene Expression Omnibus (GEO) database. Limma package of R was used to identify differentially expressed genes (DEGs). Functional enrichment analysis of DEGs were performed. Protein-protein interaction (PPI) networks of the DEGs were constructed. An integrated gene regulatory network was built including DEGs, microRNAs (miRNAs), and transcription factors. Furthermore, consistent hub genes were further validated. Molecular docking experiment was conducted. A total of 463 DEGs (232 up regulated and 231 down regulated genes) were identified between very early PDAC and normal control samples. The results of Functional enrichment analysis revealed that the DEGs were significantly enriched in vesicle organization, secretory vesicle, protein dimerization activity, lymphocyte activation, cell surface, transferase activity, transferring phosphorus-containing groups, hemostasis and adaptive immune system. The PPI network and gene regulatory network of up regulated genes and down regulated genes were established, and hub genes were identified. The expression of hub genes (CCNB1, FHL2, HLA-DPA1 and TUBB1) were also validated to be differentially expressed among PDAC and normal control samples. Molecular docking experiment predicted the novel inhibitory molecules for CCNB1 and FHL2. The identification of hub genes in PDAC enhances our understanding of the molecular mechanisms underlying the progression of this disease. These genes may be potential diagnostic biomarkers and/or therapeutic molecular targets in patients with PDAC.

2021 ◽  
Author(s):  
Varun Alur ◽  
Varshita Raju ◽  
Basavaraj Mallikarjunayya Vastrad ◽  
Anandkumar Revanasiddappa Tengli ◽  
Chanabasayya Vastrad ◽  
...  

Gestational diabetes mellitus (GDM) is the metabolic disorder appears during pregnancy. The current investigation aimed to identify central differentially expressed genes (DEGs) in GDM. The transcription profiling by array data (E-MTAB-6418) was obtained from the ArrayExpress database. The DEGs between GDM samples and non GDM samples were analyzed. Functional enrichment analysis were performed using ToppGene. Then we constructed the protein-protein interaction (PPI) network of DEGs by the Search Tool for the Retrieval of Interacting Genes database (STRING) and module analysis was performed. Subsequently, we constructed the miRNA-hub gene network and TF-hub gene regulatory network. The validation of hub genes was performed through receiver operating characteristic curve (ROC). Finally, the candidate small molecules as potential drugs to treat GDM were predicted by using molecular docking. Through transcription profiling by array data, a total of 869 DEGs were detected including 439 up regulated and 430 down regulated genes. Functional enrichment analysis showed these DEGs were mainly enriched in reproduction, cell adhesion, cell surface interactions at the vascular wall and extracellular matrix organization. Ten genes, HSP90AA1, EGFR, RPS13, RBX1, PAK1, FYN, ABL1, SMAD3, STAT3, and PRKCA were associated with GDM, according to ROC analysis. Finally, the most significant small molecules were predicted based on molecular docking. This investigation identified hub genes, signal pathways and therapeutic agents, which might help us, enhance our understanding of the mechanisms of GDM and find some novel therapeutic agents for GDM.


2020 ◽  
Vol 20 (1) ◽  
pp. 5-14 ◽  
Author(s):  
Sichao Chen ◽  
Zeming Liu ◽  
Man Li ◽  
Yihui Huang ◽  
Min Wang ◽  
...  

Aims and Objectives: Among skin cancers, malignant skin melanoma is the leading cause of death. Identification of gene markers of malignant skin melanoma associated with survival may provide new clues for prognosis prediction and treatment. This research aimed to screen out potential prognostic predictors and molecular targets for malignant skin melanoma. Introduction: Information regarding gene expression in skin melanoma and patients’ clinical traits was obtained from the Gene Expression Omnibus database. Weighted gene co-expression network analysis (WGCNA) was applied to build co-expression modules and investigate the association between the modules and clinical traits. Moreover, functional enrichment analysis was performed for clinically significant co-expression modules. Hub genes of these modules were validated via Gene Expression Profiling Interactive Analysis (GEPIA) and the Human Protein Atlas (http:// www.proteinatlas.org). Methods: First, using WGCNA, 9 co-expression modules were constructed by the top 25% differentially expressed genes (4406 genes) from 77 human melanoma samples. Two co-expression modules (magenta and blue modules) were significantly correlated with survival months (r = -0.27, p = 0.02; r = 0.27, p = 0.02, respectively). The results of functional enrichment analysis demonstrated that the magenta module was mainly enriched in the cell cycle process and the blue module was mainly enriched in the immune response process. Additionally, the GEPIA and Human Protein Atlas results suggested that the hub genes CCNB2, ARHGAP30, and SEMA4D were associated with relapse-free survival and overall survival (all p-values < 0.05) and were differentially expressed in melanoma tumors and normal skin. Results and Conclusion: The results provided the framework of co-expression gene modules of skin melanoma and screened out CCNB2, ARHGAP30, and SEMA4D associated with survival as potential prognostic predictors and molecular targets of treatment.


2020 ◽  
Vol 2020 ◽  
pp. 1-15 ◽  
Author(s):  
YunXia Liu ◽  
YeFeng Xu ◽  
Feng Xiao ◽  
JianFeng Zhang ◽  
YiQing Wang ◽  
...  

Gastric cancer (GC) is the most common malignancy of the stomach. This study was aimed at elucidating the regulatory network of circRNA-miRNA-mRNA and identifying the precise inflammation-related targets in GC. The expression profiles of GSE83521, GSE78091, and GSE33651 were obtained from the GEO database. Interactions between miRNAs and circRNAs were investigated by the Circular RNA Interactome, and targets of miRNAs were predicted with miRTarBase. Then, a circRNA/miRNA/mRNA regulatory network was constructed. Also, functional enrichment analysis of selected differentially expressed genes (DEGs) was performed. The inflammation-/GC-related targets were collected in the GeneCards and GenLiP3 database, respectively. And a protein-protein interaction (PPI) network of DE mRNAs was constructed with STRING and Cytoscape to identify hub genes. The genetic alterations, neighboring gene networks, expression levels, and the poor prognosis of hub genes were investigated in cBioPortal, Oncomine, and Human Protein Atlas databases and Kaplan-Meier plotter, respectively. A total of 10 DE miRNAs and 33 DEGs were identified. The regulatory network contained 26 circRNAs, 10 miRNAs, and 1459 mRNAs. Functional enrichment analysis revealed that the selected 33 DEGs were involved in negative regulation of fat cell differentiation, response to wounding, extracellular matrix- (ECM-) receptor interaction, and regulation of cell growth pathways. THBS1, FN1, CALM1, COL4A1, CTGF, and IGFBP5 were selected as inflammation-related hub genes of GC in the PPI network. The genetic alterations in these hub genes were related to amplification and missense mutations. Furthermore, the genes RYR2, ERBB2, PI3KCA, and HELZ2 were connected to hub genes in this study. The hub gene levels in clinical specimens were markedly upregulated in GC tissues and correlated with poor overall survival (OS). Our results suggest that THBS1, FN1, CALM1, COL4A1, CTGF, and IGFBP5 were associated with the pathogenesis of gastric carcinogenesis and may serve as biomarkers and inflammation-related targets for GC.


2020 ◽  
Author(s):  
Yiyuan Zhang ◽  
Rongguo Yu ◽  
Jiayu Zhang ◽  
Eryou Feng ◽  
Haiyang Wang ◽  
...  

Abstract BackgroundOsteoarthritis (OA) is a common chronic disease worldwide. Subchondral bone is an important pathological change in OA and responds more rapidly to adverse loading and events compared to cartilage. However, the pathogenic genes and pathways of subchondral bone are largely unclear.ObjectiveThis study aimed to identify signature differences in genes involved in knee lateral tibial (LT) and medial tibial (MT) plateaus of subchondral bone tissue while exploring their potential molecular mechanisms via bioinformatics analysis.MethodsFirst, the gene expression data of GSE51588 was downloaded from the GEO database. Differentially expressed genes (DEGs) between knee LT and MT were identified, and functional enrichment analyses were performed. Then, a protein-protein interactive network was constructed in order to acquire the hub genes, and modules analysis was conducted using STRING and Cytoscape for further analysis. The enriched hub genes were queried in DGIdb database to find suitable drug candidates in OA.ResultsA total of 202 DEGs (112 upregulated genes and 84 downregulated genes) were determined. In the PPI network, ten hub genes were identified. Five significant modules were identified using the MCODE plugin unit. Functional enrichment analysis revealed the most important signaling pathways. Six of the ten hub genes were targetable by a total of 35 drugs, suggesting their possible therapeutic use for OA .ConclusionsThe identified hub genes and functional enrichment pathways were implicated in the development and progression of subchondral bone in OA, thus improving our understanding of OA and offering molecular targets for future therapeutic modalities.


PeerJ ◽  
2021 ◽  
Vol 9 ◽  
pp. e11321
Author(s):  
Di Zhang ◽  
Pengguang Yan ◽  
Taotao Han ◽  
Xiaoyun Cheng ◽  
Jingnan Li

Background Ulcerative colitis-associated colorectal cancer (UC-CRC) is a life-threatening complication of ulcerative colitis (UC). The mechanisms underlying UC-CRC remain to be elucidated. The purpose of this study was to explore the key genes and biological processes contributing to colitis-associated dysplasia (CAD) or carcinogenesis in UC via database mining, thus offering opportunities for early prediction and intervention of UC-CRC. Methods Microarray datasets (GSE47908 and GSE87466) were downloaded from Gene Expression Omnibus (GEO). Differentially expressed genes (DEGs) between groups of GSE47908 were identified using the “limma” R package. Weighted gene co-expression network analysis (WGCNA) based on DEGs between the CAD and control groups was conducted subsequently. Functional enrichment analysis was performed, and hub genes of selected modules were identified using the “clusterProfiler” R package. Single-gene gene set enrichment analysis (GSEA) was conducted to predict significant biological processes and pathways associated with the specified gene. Results Six functional modules were identified based on 4929 DEGs. Green and blue modules were selected because of their consistent correlation with UC and CAD, and the highest correlation coefficient with the progress of UC-associated carcinogenesis. Functional enrichment analysis revealed that genes of these two modules were significantly enriched in biological processes, including mitochondrial dysfunction, cell-cell junction, and immune responses. However, GSEA based on differential expression analysis between sporadic colorectal cancer (CRC) and normal controls from The Cancer Genome Atlas (TCGA) indicated that mitochondrial dysfunction may not be the major carcinogenic mechanism underlying sporadic CRC. Thirteen hub genes (SLC25A3, ACO2, AIFM1, ATP5A1, DLD, TFE3, UQCRC1, ADIPOR2, SLC35D1, TOR1AIP1, PRR5L, ATOX1, and DTX3) were identified. Their expression trends were validated in UC patients of GSE87466, and their potential carcinogenic effects in UC were supported by their known functions and other relevant studies reported in the literature. Single-gene GSEA indicated that biological processes and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways related to angiogenesis and immune response were positively correlated with the upregulation of TFE3, whereas those related to mitochondrial function and energy metabolism were negatively correlated with the upregulation of TFE3. Conclusions Using WGCNA, this study found two gene modules that were significantly correlated with CAD, of which 13 hub genes were identified as the potential key genes. The critical biological processes in which the genes of these two modules were significantly enriched include mitochondrial dysfunction, cell-cell junction, and immune responses. TFE3, a transcription factor related to mitochondrial function and cancers, may play a central role in UC-associated carcinogenesis.


2020 ◽  
Vol 2020 ◽  
pp. 1-17
Author(s):  
Xiang Qian ◽  
Zhuo Chen ◽  
Sha Sha Chen ◽  
Lu Ming Liu ◽  
Ai Qin Zhang

The study aimed to clarify the potential immune-related targets and mechanisms of Qingyihuaji Formula (QYHJ) against pancreatic cancer (PC) through network pharmacology and weighted gene co-expression network analysis (WGCNA). Active ingredients of herbs in QYHJ were identified by the TCMSP database. Then, the putative targets of active ingredients were predicted with SwissTargetPrediction and the STITCH databases. The expression profiles of GSE32676 were downloaded from the GEO database. WGCNA was used to identify the co-expression modules. Besides, the putative targets, immune-related targets, and the critical module genes were mapped with the specific disease to select the overlapped genes (OGEs). Functional enrichment analysis of putative targets and OGEs was conducted. The overall survival (OS) analysis of OGEs was investigated using the Kaplan-Meier plotter. The relative expression and methylation levels of OGEs were detected in UALCAN, human protein atlas (HPA), Oncomine, DiseaseMeth version 2.0 and, MEXPRESS database, respectively. Gene set enrichment analysis (GSEA) was conducted to elucidate the key pathways of highly-expressed OGEs further. OS analyses found that 12 up-regulated OGEs, including CDK1, PLD1, MET, F2RL1, XDH, NEK2, TOP2A, NQO1, CCND1, PTK6, CTSE, and ERBB2 that could be utilized as potential diagnostic indicators for PC. Further, methylation analyses suggested that the abnormal up-regulation of these OGEs probably resulted from hypomethylation, and GSEA revealed the genes markedly related to cell cycle and proliferation of PC. This study identified CDK1, PLD1, MET, F2RL1, XDH, NEK2, TOP2A, NQO1, CCND1, PTK6, CTSE, and ERBB2 might be used as reliable immune-related biomarkers for prognosis of PC, which may be essential immunotherapies targets of QYHJ.


BMC Cancer ◽  
2020 ◽  
Vol 20 (1) ◽  
Author(s):  
Fuqiang Zu ◽  
Peng Liu ◽  
Huaitao Wang ◽  
Ting Zhu ◽  
Jian Sun ◽  
...  

Abstract Background It is well acknowledged that cancer-related pathways play pivotal roles in the progression of pancreatic cancer (PC). Employing Integrated analysis, we aim to identify the pathway-related ceRNA network associated with PC progression. Methods We divided eight GEO datasets into three groups according to their platform, and combined TCGA and GTEx databases as a group. Additionally, we screened out the differentially expressed genes (DEGs) and performed functional enrichment analysis in each group, and recognized the top hub genes in the most enriched pathway. Furthermore, the upstream of miRNAs and lncRNAs were predicted and validated according to their expression and prognostic roles. Finally, the co-expression analysis was applied to identify a pathway-related ceRNA network in the progression of PC. Results A total of 51 significant pathways that common enriched in all groups were spotted. Enrichment analysis indicated that pathway in cancer was greatly linked with tumor formation and progression. Next, the top 20 hug genes in this pathway were recognized, and stepwise prediction and validation from mRNA to lncRNA, including 11 hub genes, 4 key miRNAs, and 2 key lncRNAs, were applied to identify a meaningful ceRNA network according to ceRNA rules. Ultimately, we identified the PVT1/miR-20b/CCND1 axis as a promising pathway-related ceRNA axis in the progression of PC. Conclusion Overall, we elucidate the pathway-related ceRNA regulatory network of PVT1/miR-20b/CCND1 in the progression of PC, which can be considered as therapeutic targets and encouraging prognostic biomarkers for PC.


2021 ◽  
Vol 2021 ◽  
pp. 1-19
Author(s):  
Qingqing Ye ◽  
Hongbo Chen ◽  
Hongzhen Ma ◽  
Xiaojun Xiang ◽  
Shouci Hu ◽  
...  

Acute kidney injury (AKI) is responsible for significant mortality among hospitalized patients that is especially troubling aged people. An effective self-made Chinese medicine formula, Xiaoyu Xiezhuo Drink (XXD), displayed therapeutic effects on AKI. However, the compositions and underlying mechanisms of XXD remain to be elucidated. In this study, we used the ultra-high-performance liquid chromatography method coupled with hybrid triple quadrupole time-of-flight mass spectrometry (UHPLC-Q-TOF-MS) to investigate the chemical components in XXD. Then, the absorbable components of XXD were identified based on the five principles and inputted into the SwissTargetPrediction and STITCH databases to identify the drug targets. AKI-related targets were collected from the GenCLiP 3, GeneCards, and DisGeNET databases. The crossover genes of XXD and AKI were identified for functional enrichment analysis. The protein-protein interaction (PPI) network of crossover genes was constructed, followed by the identification of hub genes. Subsequently, the effects and potential mechanisms of XXD on AKI predicted by the network pharmacology and bioinformatics analyses were experimentally validated in ischemia-reperfusion (I/R) injury-induced AKI aged mouse models. A total of 122 components in XXD were obtained; among them, 58 components were found that could be absorbed in the blood. There were 800 potential drug targets predicted from the 58 absorbable components in AKI which shared 36 crossover genes with AKI-related targets. The results of functional enrichment analysis indicated that crossover genes mostly associated with the response to oxidative stress and the HIF1 signaling pathway. In the PPI network analysis, 12 hub genes were identified, including ALB, IL-6, TNF, TP53, VEGFA, PTGS2, TLR4, NOS3, EGFR, PPARG, HIF1A, and HMOX1. In AKI aged mice, XXD prominently alleviated I/R injury-induced renal dysfunction, abnormal renal pathological changes, and cellular senescence, inflammation, and oxidative damage with a reduction in the expression level of the inflammatory mediator, α-SMA, collagen-1, F4/80, TP53, VEGFA, PTGS2, TLR4, NOS3, EGFR, PPARG, HIF1A, ICAM-1, TGF-β1, Smad3, and p-Smad3 and an increase of nephridial tissue p-H3, Ki67, HMOX1, MMP-9, and Smad7 levels. In summary, our findings suggest that XXD has renoprotective effects against AKI in aged mice via inhibiting the TGF-β1/Smad3 and HIF1 signaling pathways.


2021 ◽  
Vol 2021 ◽  
pp. 1-13
Author(s):  
Wei Xu ◽  
Jian Xu ◽  
Zhiqiang Wang ◽  
Yuequan Jiang

Objective. Esophageal cancer (ESCA) is one of the most aggressive malignancies globally with an undesirable five-year survival rate. Here, this study was conducted for determining specific functional genes linked with ESCA initiation and progression. Methods. Gene expression profiling of ESCA was curated from TCGA (containing 160 ESCA and 11 nontumor specimens) and GSE38129 (30 paired ESCA and nontumor tissues) datasets. Differential expression analysis was conducted between ESCA and nontumor tissues with adjusted p value <0.05 and |log2fold-change|>1. Weighted gene coexpression network analysis (WGCNA) was conducted for determining the ESCA-specific coexpression modules and genes. Thereafter, ESCA-specific differentially expressed genes (DEGs) were intersected. Functional enrichment analysis was then presented with clusterProfiler package. Protein-protein interaction was conducted, and hub genes were determined. Association of hub genes with pathological staging was evaluated, and survival analysis was presented among ESCA patients. Results. This study determined 91 ESCA-specific DEGs following intersection of DEGs and ESCA-specific genes in TCGA and GSE38129 datasets. They were remarkably linked to cell cycle progression and carcinogenic pathways like the p53 signaling pathway, cellular senescence, and apoptosis. Ten ESCA-specific hub genes were determined, containing ASPM, BUB1B, CCNA2, CDC20, CDK1, DLGAP5, KIF11, KIF20 A, TOP2A, and TPX2. They were prominently associated with pathological staging. Among them, KIF11 upregulation was in relation to undesirable prognosis of ESCA patients. Conclusion. Collectively, we determined ESCA-specific coexpression modules and hub genes, which offered the foundation for future research concerning the mechanistic basis of ESCA.


2021 ◽  
Author(s):  
Basavaraj Mallikarjunayya Vastrad ◽  
Chanabasayya Mallikarjunayya Vastrad

Heart failure (HF) is a complex cardiovascular diseases associated with high mortality. To discover key molecular changes in HF, we analyzed next-generation sequencing (NGS) data of HF. In this investigation, differentially expressed genes (DEGs) were analyzed using limma in R package from GSE161472 of the Gene Expression Omnibus (GEO). Then, gene enrichment analysis, protein-protein interaction (PPI) network, miRNA-hub gene regulatory network and TF-hub gene regulatory network construction, and topological analysis were performed on the DEGs by the Gene Ontology (GO), REACTOME pathway, STRING, HiPPIE, miRNet, NetworkAnalyst and Cytoscape. Finally, we performed receiver operating characteristic curve (ROC) analysis of hub genes. A total of 930 DEGs 9464 up regulated genes and 466 down regulated genes) were identified in HF. GO and REACTOME pathway enrichment results showed that DEGs mainly enriched in localization, small molecule metabolic process, SARS-CoV infections and the citric acid (TCA) cycle and respiratory electron transport. Subsequently, the PPI network, miRNA-hub gene regulatory network and TF-hub gene regulatory network were constructed, and 10 hub genes in these network were focused on by centrality analysis and module analysis. Furthermore, data showed that HSP90AA1, ARRB2, MYH9, HSP90AB1, FLNA, EGFR, PIK3R1, CUL4A, YEATS4 and KAT2B were good diagnostic values. In summary, this study suggests that HSP90AA1, ARRB2, MYH9, HSP90AB1, FLNA, EGFR, PIK3R1, CUL4A, YEATS4 and KAT2B may act as the key genes in HF.


Sign in / Sign up

Export Citation Format

Share Document