Identification of Key Candidate Genes and Pathways in Colorectal Cancer by Integrated Bioinformatical Analysis

Colorectal cancer (CRC) is one of the most common malignant diseases worldwide, but the involved signaling pathways and driven-genes are largely unclear. This study integrated four cohorts profile datasets to elucidate the potential key candidate genes and pathways in CRC. Expression profiles GSE28000, GSE21815, GSE44076 and GSE75970, including 319 CRC and 103 normal mucosa, were integrated and deeply analyzed. Differentially expressed genes (DEGs) were sorted and candidate genes and pathways enrichment were analyzed. DEGs-associated protein–protein interaction network (PPI) was performed. Firstly, 292 shared DEGs (165 up-regulated and 127 down-regulated) were identified from the four GSE datasets. Secondly, the DEGs were clustered based on functions and signaling pathways with significant enrichment analysis. Thirdly, 180 nodes/DEGs were identified from DEGs PPI network complex. Lastly, the most significant 2 modules were filtered from PPI, 31 central node genes were identified and most of the corresponding genes are involved in cell cycle process, chemokines and G protein-coupled receptor signaling pathways. Taken above, using integrated bioinformatical analysis, we have identified DEGs candidate genes and pathways in CRC, which could improve our understanding of the cause and underlying molecular events, and these candidate genes and pathways could be therapeutic targets for CRC.

Identification of core genes and outcomes in hepatocellular carcinoma by bioinformatics analysis

Hepatocellular carcinoma (HCC) is the most common malignant liver disease in the world. However, the mechanistic relationships among various genes and signaling pathways are still largely unclear. In this study, we aimed to elucidate potential core candidate genes and pathways in HCC. The expression profiles GSE14520, GSE25097, GSE29721, and GSE62232, which cover 606 tumor and 550 nontumour samples, were downloaded from the Gene Expression Omnibus (GEO) database. Furthermore, HCC RNA-seq datasets were also downloaded from the Cancer Genome Atlas (TCGA) database. The differentially expressed genes (DEGs) were filtered using R software, and we performed gene ontology (GO) and Kyoto Encyclopedia of Gene and Genome (KEGG) pathway analysis using the online databases DAVID 6.8 and KOBAS 3.0. Furthermore, the protein-protein interaction (PPI) network complex of these DEGs was constructed by Cytoscape software, the molecular complex detection (MCODE) plug-in and the online database STRING. First, a total of 173 DEGs (41 upregulated and 132 downregulated) were identified that were aberrantly expressed in both the GEO and TCGA datasets. Second, GO analysis revealed that most of the DEGs were significantly enriched in extracellular exosomes, cytosol, extracellular region, and extracellular space. Signaling pathway analysis indicated that the DEGs had common pathways in metabolism-related pathways, cell cycle, and biological oxidations. Third, 146 nodes were identified from the DEG PPI network complex, and two important modules with a high degree were detected using the MCODE plug-in. In addition, 10 core genes were identified, TOP2A, NDC80, FOXM1, HMMR, KNTC1, PTTG1, FEN1, RFC4, SMC4, and PRC1. Finally, Kaplan-Meier analysis of overall survival and correlation analysis were applied to these genes. The abovementioned findings indicate that the identified core genes and pathways in this bioinformatics analysis could significantly enrich our understanding of the development and recurrence of HCC; furthermore, these candidate genes and pathways could be therapeutic targets for HCC treatment.

Targeting glypican-4 overcomes 5-FU resistance and attenuates stem cell-like properties via suppression of Wnt/β-catenin pathway in pancreatic cancer cells

The existences of cancer stem cells in patients with pancreatic cancer are considered as pivotal factors contributing to chemoresistance and disease relapse. Glypican-4 (GPC4) is one of the members of the glypicans family, which underlies human congenital malformations and multiple diseases. However, its potential biological function in pancreatic cancer still remains elusive. In this study, we are the first to demonstrate that GPC4 was involved in 5-fluorouracil (5-FU) resistance and pancreatic cancer stemness through comprehensive bioinformatical analysis. Functional experiments showed that knockdown of GPC4 sensitized pancreatic cancer cells to 5-FU and attenuated stem cell-like properties. In terms of mechanism research, knockdown of GPC4 suppressed the activation of Wnt/β-catenin pathway and its downstream targets. Furthermore, the expression of GPC4 was significantly upregulated in pancreatic cancer tissues compared with normal tissues and remarkably correlated with patients' overall survival according to the data derived from the Cancer Genome Atlas database. Taken together, our results suggest that GPC4 is a key regulator in chemoresistance and pancreatic cancer stemness. Thus, targeting GPC4 may serve as a promising strategy for pancreatic cancer therapy.

Identification of potential crucial genes associated with carcinogenesis of clear cell renal cell carcinoma

Clear cell renal cell carcinoma (ccRCC) is a common genitourinary malignancy with high mortality. However, the molecular pathogenesis of ccRCC remains unclear and effective biomarkers for daily practice are still limited. Thus, we aimed to identify the potential crucial genes and pathways associated with carcinogenesis of ccRCC and further analyze the molecular mechanisms implicated in tumorigenesis. In the present study, expression profiles GSE 66270, GSE 53757, GSE 36895, and GSE 76351 were downloaded from GEO database, including 244 matched primary and adjacent normal tissues, furthermore, the level 3 RNAseq dataset (RNAseqV2 RSEM) of KIRC was also downloaded from The Cancer Genome Atlas (TCGA), which consist of 529 ccRCC tumors and 72 normal tissues. Then, differentially expressed genes (DEGs) and pathway enrichment were analyzed by using R software. A total of 129 up- and 123 down-regulated genes were identified, which were aberrantly expressed both in GEO and TCGA data. Second, Gene ontology (GO) analyses revealed that most of the DEGs were significantly enriched in integral component of membrane, extracellular exosome, plasma membrane, cell adhesion, and receptor binding. Signaling pathway analyses indicated that DEGs had common pathways in signal transduction, metabolism, and immune system. Third, hub genes were identified with protein-protein interaction (PPI) network, including PTPRC, TGFB1, EGF, MYC, ITGB2, CTSS, FN1, CCL5, KNG1, and CD86. Additionally, sub-networks analyse was also performed by using MCODE plugin. In conclusion, the novel DEGs and pathways in ccRCC identified in this study may provide new insight into the underlying molecular mechanisms that facilitates RCC carcinogenesis.

ASPM promotes glioblastoma growth by regulating G1 restriction point progression and Wnt-β-catenin signaling

Increasing evidence has indicated that the disorganized expression of certain genes promotes tumour progression. In this study, we elucidate the potential key differentially expressed genes (DEGs) between glioblastoma (GBM) and normal brain tissue by analysing three different mRNA expression profiles downloaded from the Gene Expression Omnibus (GEO) database. DEGs were sorted, and key candidate genes and signalling pathway enrichments were analysed. In our analysis, the highest fold change DEG was found to be abnormal spindle-like microcephaly associated (ASPM). The ASPM expression pattern from the database showed that it is highly expressed in GBM tissue, and patients with high expression of ASPM have a poor prognosis. Moreover, ASPM showed aberrantly high expression in GBM cell lines. Loss-of-function assay indicated that ASPM enhances tumorigenesis in GBM cells in vitro. Xenograft growth verified the oncogenic activity of ASPM in vivo. Furthermore, downregulation of ASPM could arrest the cell cycle of GBM cells at the G0/G1 phase and attenuate the Wnt/β-catenin signalling activity in GBM. These data suggest that ASPM may serve as a new target for the therapeutic treatment of GBM.

Comprehensive Analysis of the Expression and Prognosis for Laminin Genes in Ovarian Cancer

Survival is low in ovarian cancer (OC). Most OC patients demonstrate advanced metastases, and recurrence is common. Dysregulation of laminin interactions is associated with cancer development. However, it is unknown whether laminin subunits can be considered as biomarkers for OC diagnosis, prognosis, and treatment. We used cBioPortal, GEO, ONCOMINE, GEPIA, Human Protein Atlas, Kaplan-Meier Plotter, TIMER, and Metascape to determine the associations among laminin expression, prognosis, and immune cell infiltration in OC. LAMA5, LAMB3, and LAMC2 mRNAs and LAMA3, LAMB1/B2/B3, and LAMC1/C2 proteins were overexpressed in OC tissues compared with normal ovaries. LAMA4, LAMB1, and LAMC1 mRNA upregulation was positively correlated with worse overall survival (OS) and progression-free survival (PFS) in OC. Elevated LAMA2 and LAMC2 mRNA expression levels were related to better PFS or OS, respectively. The results speculated that LAMA5 could potentially be a good prognostic factor in OC. Its expression proves valuable for predicting OS in patients diagnosed with stage Ⅳ and grade 3 OC and PFS in patients diagnosed with all OC stages or grades. LAMB3 and LAMC2 expression was correlated with platinum resistance development. ROC analysis of laminins in OC sets revealed that LAMA2/A4/A5, LAMB1/B2/B3, and LAMC2 could be used to differentiate between malignant tumors and non-neoplastic tissues. LAMA1/A5 and LAMC1 were significantly and negatively correlated with various tumor immune infiltrates (TILs), especially with dendritic cells, CD8⁺ T cells or neutrophil. LAMA4 and LAMB1 might be associated with tumor purity in OC. Overall, LAMA5 and LAMC1 could help predict OC survival and diagnosis and might be deemed important OC oncogenes.

Identification of differentially expressed genes in synovial tissue of rheumatoid arthritis and osteoarthritis in patients

Rheumatoid arthritis (RA) and osteoarthritis (OA) are the common joints disorder in the world. Although they have showed the analogous clinical manifestation and overlapping cellular and molecular foundation, the pathogenesis of RA and OA were different. The pathophysiologic mechanisms of arthritis in RA and OA have not been investigated thoroughly. Thus, the aim of study is to identify the potential crucial genes and pathways associated with RA and OA and further analyze the molecular mechanisms implicated in genesis. First, we compared gene expression profiles in synovial tissue between RA and OA from the National Center of Biotechnology Information (NCBI) Gene Expression Omnibus (GEO) database. Gene Expression Series (GSE) 1919, GSE55235, and GSE36700 were downloaded from the GEO database, including 20 patients of OA and 21 patients of RA. Differentially expressed genes (DEGs) including "CXCL13," "CD247," "CCL5," "GZMB," "IGKC," "IL7R," "UBD///GABBR1," "ADAMDEC1," "BTC," "AIM2," "SHANK2," "CCL18," "LAMP3," "CR1," and "IL32." Second, Gene Ontology analyses revealed that DEGs were significantly enriched in integral component of extracellular space, extracellular region, and plasma membrane in the molecular function group. Signaling pathway analyses indicated that DEGs had common pathways in chemokine signaling pathway, cytokine-cytokine receptor interaction, and cytosolic DNA-sensing pathway. Third, DEGs showed the complex DEGs protein-protein interaction network with the Coexpression of 83.22%, Shared protein domains of 8.40%, Colocalization of 4.76%, Predicted of 2.87%, and Genetic interactions of 0.75%. In conclusion, the novel DEGs and pathways between RA and OA identified in this study may provide new insight into the underlying molecular mechanisms of RA.

Identification of Potential Biomarkers in Association With Progression and Prognosis in Epithelial Ovarian Cancer by Integrated Bioinformatics Analysis

Epithelial ovarian cancer (EOC) is one of the malignancies in women, which has the highest mortality. However, the microlevel mechanism has not been discussed in detail. The expression profiles GSE27651, GSE38666, GSE40595, and GSE66957 including 188 tumor and 52 nontumor samples were downloaded from the Gene Expression Omnibus database. The differentially expressed genes (DEGs) were filtered using R software, and we performed functional analysis using the clusterProfiler. Cytoscape software, the molecular complex detection plugin and database STRING analyzed DEGs to construct protein-protein interaction network. We identified 116 DEGs including 81 upregulated and 35 downregulated DEGs. Functional analysis revealed that they were significantly enriched in the extracellular region and biosynthesis of amino acids. We next identified four bioactive compounds (vorinostat, LY-294002,trichostatin A, and tanespimycin) based on ConnectivityMap. Then 114 nodes were obtained in protein-protein interaction. The three most relevant modules were detected. In addition, according to degree ≥ 10, 14 core genes including FOXM1, CXCR4, KPNA2, NANOG, UBE2C, KIF11, ZWINT, CDCA5, DLGAP5, KIF15, MCM2, MELK, SPP1, and TRIP13 were identified. Kaplan-Meier analysis, Oncomine, and Gene Expression Profiling Interactive Analysis showed that overexpression of FOXM1, SPP1, UBE2C, KIF11, ZWINT, CDCA5, UBE2C, and KIF15 was related to bad prognosis of EOC patients. CDCA5, FOXM1, KIF15, MCM2, and ZWINT were associated with stage. Receiver operating characteristic (ROC) curve showed that messenger RNA levels of these five genes exhibited better diagnostic efficiency for normal and tumor tissues. The Human Protein Atlas database was performed. The protein levels of these five genes were significantly higher in tumor tissues compared with normal tissues. Functional enrichment analysis suggested that all the hub genes played crucial roles in citrate cycle tricarboxylic acid cycle. Furthermore, the univariate and multivariate Cox proportional hazards regression showed that ZWINT was independent prognostic indictor among EOC patients. The genes and pathways discovered in the above studies may open a new direction for EOC treatment.

Identification of Key Biomarkers and Potential Molecular Mechanisms in Renal Cell Carcinoma by Bioinformatics Analysis

Renal cell carcinoma (RCC) is the most common form of kidney cancer, caused by renal epithelial cells. RCC remains to be a challenging public health problem worldwide. Metastases that are resistant to radiotherapy and chemotherapy are the major cause of death from cancer. However, the underlying molecular mechanism regulating the metastasis of RCC is poorly known. Publicly available databases of RCC were obtained from Gene Expression Omnibus (GEO) database. The differentially expressed genes (DEGs) were identified using GEO2R analysis, whereas the Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were performed by Gene Set Enrichment Analysis (GSEA) and Metascape. Protein-protein interaction (PPI) network of DEGs was analyzed by STRING online database, and Cytoscape software was used for visualizing PPI network. Survival analysis of hub genes was conducted using GEPIA online database. The expression levels of hub genes were investigated from The Human Protein Atlas online database and GEPIA online database. Finally, the comparative toxicogenomics database (CTD; http://ctdbase.org) was used to identify hub genes associated with tumor or metastasis. We identified 229 DEGs comprising 135 downregulated genes and 94 upregulated genes. Functional analysis revealed that these DEGs were associates with cell recognition, regulation of immune, negative regulation of adaptive immune response, and other functions. And these DEGs mainly related to P53 signaling pathway, cytokine-cytokine receptor interaction, Natural killer cell mediated cytotoxicity, and other pathways are involved. Ten genes were identified as hub genes through module analyses in the PPI network. Finally, survival analysis of 10 hub genes was conducted, which showed that the MMP2 (matrix metallo peptidase 2), DCN, COL4A1, CASR (calcium sensing receptor), GPR4 (G protein-coupled receptor 4), UTS2 (urotensin 2), and LDLR (low density lipoprotein receptor) genes were significant for survival. In this study, the DEGs between RCC and metastatic RCC were analyzed, which assist us in systematically understanding the pathogeny underlying metastasis of RCC. The MMP2, DCN, COL4A1, CASR, GPR4, UTS2, and LDLR genes might be used as potential targets to improve diagnosis and immunotherapy biomarkers for RCC.

GENE-IS: Time-Efficient and Accurate Analysis of Viral Integration Events in Large-Scale Gene Therapy Data

Integration site profiling and clonality analysis of viral vector distribution in gene therapy is a key factor to monitor the fate of gene-corrected cells, assess the risk of malignant transformation, and establish vector biosafety. We developed the Genome Integration Site Analysis Pipeline (GENE-IS) for highly time-efficient and accurate detection of next-generation sequencing (NGS)-based viral vector integration sites (ISs) in gene therapy data. It is the first available tool with dual analysis mode that allows IS analysis both in data generated by PCR-based methods, such as linear amplification method PCR (LAM-PCR), and by rapidly evolving targeted sequencing (e.g., Agilent SureSelect) technologies. GENE-IS makes use of trimming strategies, customized reference genome, and soft-clipped information with sequential filtering steps to provide annotated IS with clonality information. It is a scalable, robust, precise, and reliable tool for large-scale pre-clinical and clinical data analysis that provides users complete flexibility and control over analysis with a broad range of configurable parameters. GENE-IS is available at https://github.com/G100DKFZ/gene-is.

Identification of core genes and clinical roles in pregnancy-associated breast cancer based on integrated analysis of different microarray profile datasets

More women are delaying child-birth. Thus, the diagnosis of pregnancy-associated breast cancer (PABC) will continue to increase. The aim of this study was to identify core candidate genes of PABC, and the relevance of the genes on the prognosis of PABC. GSE31192 and GSE53031 microarray profile datasets were downloaded from the Gene Expression Omnibus database and differentially expressed genes were analyzed using the R package and GEO2R tool. Then, Gene Ontology and Kyoto Encyclopedia of Gene and Genome pathway enrichment analyses were performed using the Database for Annotation, Visualization, and Integrated Discovery. Moreover, the Search Tool for the Retrieval of Interacting Genes and the Molecular Complex Detection Cytoscape software plug-in were utilized to visualize protein–protein interactions and to screen candidate genes. A total of 239 DEGs were identified in PABC, including 101 up-regulated genes mainly enriched in fatty acid activation and the fibroblast growth factor signaling pathway, while 138 down-regulated genes particularly involved in activation of DNA fragmentation factor and apoptosis-induced DNA fragmentation. Fourteen hub genes with a high degree of connectivity were selected, including CREB1, ARF3, UBA5, SIAH1, KLHL3, HECTD1, MMP9, TRIM69, MEX3C, ASB6, UBE2Q2, FBXO22, EIF4A3, and PXN. Overall survival (OS) analysis of core candidate genes was performed using the Gene Expression Profiling Interactive Analysis and UALCAN websites. High ASB6 expression was associated with worse OS of PABC patients. Molecular subtypes and menopause status were also associated with worse OS for PABC patients. In conclusion, ASB6 could be a potential predictor and therapeutic target in patient with PABC.

Identification of Key Candidate Genes and Pathways in Endometrial Cancer by Integrated Bioinformatical Analysis

Endometrial Cancer is the most common female genital tract malignancy, its pathogenesis is complex, not yet fully described. To identify key genes of Endometrial Cancer we downloaded the gene chip GSE17025 from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) were identified through the GEO2R analysis tool. Functional and pathway enrichment analysis were performed for DEGs using DAVID database. The network of protein–protein-interaction (PPI) was established by STRING website and visualized by Cytoscape. Then, functional and pathway enrichment analysis of DEGS were performed by DAVID database. A total of 1000 significant differences genes were obtained, contain 362 up-regulated genes and 638 down-regulated genes. PCDH10, SLC6A2, OGN, SFRP4, TRH, ANGPTL, FOSB are down-regulated genes. The gene of IGH, CCL20, ELF5, LTF, ASPM expression level in tumor patients are up-regulated. Biological function of enrichment include metabolism of xenobiotics by cytochrome P450, MAPK signaling pathway, Serotonergic synapse, Protein digestion and absorption, IL-17 signaling pathway, Chemokine signaling pathway, HIF-1 signaling pathway, p53 signaling pathway. All in all, the current study to determine endometrial differentially expressed genes and biological function, comprehensive analysis of intrauterine membrane carcinoma pathogenesis mechanism, and might be used as molecular targets and diagnostic biomarkers for the treatment of endometrial cancer.

Identification key genes, key miRNAs and key transcription factors of lung adenocarcinoma

Background

Lung adenocarcinoma (LUAD) is one of the most common cancers worldwide. The etiology and pathophysiology of LUAD remain unclear. The aim of the present study was to identify the key genes, miRNAs and transcription factors (TFs) associated with the pathogenesis and prognosis of LUAD.

Methods

Three gene expression profiles (GSE43458, GSE32863, GSE74706) of LUAD were obtained from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) were identified by GEO2R.The Gene Ontology (GO) terms, pathways, and protein-protein interactions (PPIs) of these DEGs were analyzed. Bases on DEGs, the miRNAs and TFs were predicted. Furthermore, TF-gene-miRNA co-expression network was constructed to identify key genes, miRNAs and TFs by bioinformatic methods. The expressions and prognostic values of key genes, miRNAs and TFs were carried out through The Cancer Genome Atlas (TCGA) database and Kaplan Meier-plotter (KM) online dataset.

Results

A total of 337 overlapped DEGs (75 upregulated and 262 downregulated) of LUAD were identified from the three GSE datasets. Moreover, 851 miRNAs and 29 TFs were identified to be associated with these DEGs. In total, 10 hub genes, 10 key miRNAs and 10 key TFs were located in the central hub of the TF-gene-miRNA co-expression network, and validated using The Cancer Genome Atlas (TCGA) database. Specifically, seven genes (PHACTR2, MSRB3, GHR, PLSCR4, EPB41L2, NPNT, FBXO32), two miRNAs (hsa-let-7e-5p, hsa-miR-17-5p) and four TFs (STAT6, E2F1, ETS1, JUN) were identified to be associated with prognosis of LUAD, which have significantly different expressions between LUAD and normal lung tissue. Additionally, the miRNA/gene co-expression analysis also revealed that hsa-miR-17-5p and PLSCR4 have a significant negative co-expression relationship (r=−0.33, P=1.67e-14) in LUAD.

Conclusions

Our study constructed a regulatory network of TF-gene-miRNA in LUAD, which may provide new insights about the interaction between genes, miRNAs and TFs in the pathogenesis of LUAD, and identify potential biomarkers or therapeutic targets for LUAD.

Identification of a Common Different Gene Expression Signature in Ischemic Cardiomyopathy

The molecular mechanisms underlying the development of ischemic cardiomyopathy (ICM) remain poorly understood. Gene expression profiling is helpful to discover the molecular changes taking place in ICM. The aim of this study was to identify the genes that are significantly changed during the development of heart failure caused by ICM. The differentially expressed genes (DEGs) were identified from 162 control samples and 227 ICM patients. PANTHER was used to perform gene ontology (GO), and Reactome for pathway enrichment analysis. A protein–protein interaction network was established using STRING and Cytoscape. A further validation was performed by real-time polymerase chain reaction (RT-PCR). A total of 255 common DEGs was found. Gene ontology, pathway enrichment, and protein–protein interaction analysis showed that nucleic acid-binding proteins, enzymes, and transcription factors accounted for a great part of the DEGs, while immune system signaling and cytokine signaling displayed the most significant changes. Furthermore, seven hub genes and nine transcription factors were identified. Interestingly, the top five upregulated DEGs were located on chromosome Y, and four of the top five downregulated DEGs were involved in immune and inflammation signaling. Further, the top DEGs were validated by RT-PCR in human samples. Our study explored the possible molecular mechanisms of heart failure caused by ischemic heart disease.

Identification of candidate aberrantly methylated and differentially expressed genes in thyroid cancer

Aberrant methylation of DNA sequences plays a criticle role in finding novel aberrantly methylated genes and pathways in thyroid cancer (THCA). This study aimed to integrate three cohorts profile datasets to find novel aberrantly methylated genes and pathways in THCA. Data of gene expression profiling microarrays (GSE33630 and GSE65144) and gene methylation profiling microarrays (GSE51090) were downloaded from the Gene Expression Omnibus database. Aberrantly methylated and differentially expressed genes were sorted and pathways were analyzed. Functional and enrichment analyses of selected genes were performed using the String database. A protein‐protein interaction network was constructed using the Cytoscape software, and module analysis was performed using Molecular Complex detection. In total, we identified 12 hypomethylation/high‐expression genes and 30 hypermethylation/low‐expression genes at the screening step and, finally, found 6 mostly changed hub genes including PPARGC1A, CREBBP, EP300, CD44, SPP1, and MMP9. Pathway analysis showed that aberrantly methylated differentially expressed genes were mainly associated with the thyroid hormone signaling pathway, AMP‐activated protein kinase (AMPK) signaling pathway, and cell cycle process in THCA. After validation in the Cancer Genome Atlas database, the methylation and expression status of hub genes was significantly altered and the same with our results. Taken together, we identified novel aberrantly methylated genes and pathways in THCA, which could improve our understanding of the cause and underlying molecular events, and these candidate genes could serve as aberrant methylation‐based biomarkers for precise diagnosis and treatment of THCA.

Identification of key candidate genes and pathways in hepatitis B virus-associated acute liver failure by bioinformatical analysis

Hepatitis B virus-associated acute liver failure (HBV-ALF) is a rare but life-threatening syndrome that carried a high morbidity and mortality. Our study aimed to explore the possible molecular mechanisms of HBV-ALF by means of bioinformatics analysis. In this study, genes expression microarray datasets of HBV-ALF from Gene Expression Omnibus were collected, and then we identified differentially expressed genes (DEGs) by the limma package in R. After functional enrichment analysis, we constructed the protein-protein interaction (PPI) network by the Search Tool for the Retrieval of Interacting Genes online database and weighted genes coexpression network by the WGCNA package in R. Subsequently, we picked out the hub genes among the DEGs. A total of 423 DEGs with 198 upregulated genes and 225 downregulated genes were identified between HBV-ALF and normal samples. The upregulated genes were mainly enriched in immune response, and the downregulated genes were mainly enriched in complement and coagulation cascades. Orosomucoid 1 (ORM1), orosomucoid 2 (ORM2), plasminogen (PLG), and aldehyde oxidase 1 (AOX1) were picked out as the hub genes that with a high degree in both PPI network and weighted genes coexpression network. The weighted genes coexpression network analysis found out 3 of the 5 modules that upregulated genes enriched in were closely related to immune system. The downregulated genes enriched in only one module, and the genes in this module majorly enriched in the complement and coagulation cascades pathway. In conclusion, 4 genes (ORM1, ORM2, PLG, and AOX1) with immune response and the complement and coagulation cascades pathway may take part in the pathogenesis of HBV-ALF, and these candidate genes and pathways could be therapeutic targets for HBV-ALF.

Identification of Common Differentially Expressed Genes and Potential Therapeutic Targets in Ulcerative Colitis and Rheumatoid Arthritis

Ulcerative colitis (UC) and rheumatoid arthritis (RA) are immune-mediated inflammatory diseases (IMIDs) with similar symptoms and common genomics. However, the relationship between UC and RA has not been investigated thoroughly. Therefore, this study aimed to establish the differentially expressed genes (DEGs) and potential therapeutic targets in UC and RA. Three microarray datasets (GSE38713, GSE1919, and GSE12251) were selected from the Gene Expression Omnibus (GEO) database for analysis. We used R software to identify the DEGs and performed enrichment analyses. Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) and Cytoscape software were used to construct the protein-protein interaction (PPI) network and identify the hub genes. A regulatory network based on the constructed PPI was generated using StarBase and PROMO databases. We identified a total of 1542 and 261 DEGs in UC and RA. There were 169 common DEGs identified in both UC and RA, including 63 upregulated genes (DEGs1) and nine downregulated genes (DEGs2). The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses of DEGs1 and DEGs2 in the PPI network revealed that the genes enriched were involved in immunity. A total of 45 hub genes were selected based on high scores of correlation; three hub genes (SRGN, PLEK, and FCGR3B) were found to be upregulated in UC and RA, and downregulated in UC patients with response to infliximab treatment. The identification of novel DEGs and hub genes in the current study contributes to a novel perception for latent functional mechanisms and presents potential prognostic indicators and therapeutic targets in UC and RA.

Identification of key candidate genes and pathways in hepatocellular carcinoma by integrated bioinformatical analysis

Hepatocellular carcinoma (HCC) is one of the most common malignant neoplasms worldwide, however the underlying mechanisms and gene signatures of HCC are unknown. In the present study the profile datasets of four cohorts were integrated to elucidate the pathways and candidate genes of HCC. The expression profiles GSE25097, GSE45267, GSE57957 and GSE62232 were downloaded from the Gene Expression Omnibus database, including 436 HCC and 94 normal liver tissues. A total of 185 differentially expressed genes (DEGs) were identified in HCC, including 92 upregulated genes and 92 downregulated genes. Gene ontology (GO) was performed, which revealed that the upregulated DEGs were primarily enriched in cell division, mitotic nuclear division, mitotic cytokinesis and G1/S transition of the mitotic cell cycle. Pathway enrichment was analyzed based on the Kyoto Encyclopedia of Genes and Genomes database to assess the functional relevance of DEGs. The most significant module was selected from protein-protein interactions and 15 important hub genes were identified. The sub-networks of hub genes were involved in cell division, p53 signaling, and T lymphotropic virus type I infection signaling pathways. In conclusion, the present study revealed that the identified DEG candidate genes may promote the understanding of the cause and molecular mechanisms underlying the development of HCC and that these candidates and signal pathways may be potential targets of clinical therapy for HCC.

Integrative bioinformatics analysis of miRNA and mRNA expression profiles and identification of associated miRNA-mRNA network in aortic dissection

BACKGROUND

Aortic dissection (AD) is one of the most lethal cardiovascular diseases. The aim of this study was to identify core genes and pathways revealing pathogenesis in AD.

METHODS

We screened differentially expressed mRNAs and miRNAs using mRNA and miRNA expression profile data of AD from Gene Expression Omnibus. Then functional and pathway enrichment analyses of differential expression genes (DEGs) was performed utilizing the database for annotation, visualization, and integrated discovery (DAVID). Target genes with differential expression miRNAs (DEMIs) were predicted using the miRWalk database, and the intersection between these predictions and DEGs was selected as differentially expressed miRNA-target genes. In addition, a protein-protein interaction (PPI) network and miRNA-mRNA regulatory network were constructed.

RESULTS

In total, 130 DEGs and 47 DEMIs were identified from mRNA and miRNA microarray, respectively, and 45 DEGs were DEMI-target genes. The PPI and miRNA-mRNA network included 79 node genes and 74 node genes, respectively, while 23 hub genes and 2 hub miRNAs were identified. The DEGs, PPI and modules differential expression miRNA-target genes were all mainly enriched in cell cycle, cell proliferation and cell apoptosis signaling pathways.

CONCLUSION

Taken above, the study reveals some candidate genes and pathways potentially involving molecular mechanisms of AD. These findings provide a new insight for research and treatment of AD.

Identification of key candidate genes and pathways in oral squamous cell carcinoma by integrated Bioinformatics analysis

Oral squamous cell carcinoma (OSCC) is one of the most common types of malignant head and neck tumor, which poses a serious threat to human health. In recent years, the incidence of OSCC has been increasing, while the prognosis has not significantly improved. Elucidation of the molecular mechanisms underlying the development of OSCC may provide novel therapeutic strategies. In the present study, the gene expression profiles from 4 datasets, including 244 OSCC and 95 normal oral mucosa samples, were subjected to statistical and Bioinformatics analysis. A total of 34 differentially expressed genes (DEGs) were identified, among which 14 were upregulated and 20 were downregulated in OSCC compared with normal oral mucosa tissues. Gene Ontology enrichment analysis indicated that the DEGs were mainly involved in regulation of the immune response, cell adhesion and cell proliferative processes. The Kyoto Encyclopedia of Genes and Genomes pathway analysis revealed that the DEGs were mainly associated with the phosphoinositide-3 kinase Akt and Toll-like receptor signaling pathway. The key candidate DEGs were identified from the complex protein-protein interaction network, and secreted phosphoprotein 1 (SPP1), integrin subunit α 3 and plasminogen activator, urokinase (PLAU) were confirmed to be significantly associated with the survival rate. Cell Counting Kit-8 and Transwell assays demonstrated that SPP1 and PLAU regulate cell proliferation, migration and invasion. The candidate genes/pathways identified in the present study may include promising diagnostic biomarkers or therapeutic targets for OSCC.

Bioinformatics analysis and identification of potential genes related to pathogenesis of cervical intraepithelial neoplasia

The aim of this study was to explore and identify the key genes and signal pathways contributing to cervical intraepithelial neoplasia (CIN). The gene expression profiles of GSE63514 were downloaded from Gene Expression Omnibus database. Differentially expressed genes (DEGs) were screened performing with packages in software R. After Gene ontology terms, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyzing, and Gene set enrichment analysis (GSEA), weighted gene co-expression network analysis (WGCNA) was used to analyze these genes. Then sub-modules were subsequently analyzed base CIN grade, and protein-protein interaction (PPI) network of DEGs were constructed. 537 DEGs were screened in total, consisting 331 up-regulated genes and 206 down-regulated genes in CIN samples compared to normal samples. The most DEGs were enriched in chromosomal region in cellular component (CC), organelle fission inbiological process (BP) and ATPase activity in molecular function (MF). KEGG pathway enrichment analyzing found the DEGs were mainly concentrated in 10 pathways. The results of GSEA mainly enriched in 4 functional sets: E2F-Targets, G2M-Checkpoint, Mitotic-Spindle and Spermatogenesis. A total of 6 modules were identified by WCGNA. Subsequently, grey module was the highest correlation (Cor=0.78, P=5e-22) and 31 genes were taken as candidate hub genes for CIN high grade risk (weighted correlation coefficients >0.80). Finally, diagnostic analysis showed that in addition to CCDC7, the expression levels of the remaining 13 DEGs have a high diagnostic value (AUC>0.8 and P<0.05). These findings provided a new sight into the understanding of molecular functions for CIN.

Identification of Potential Core Genes Associated With the Progression of Stomach Adenocarcinoma Using Bioinformatic Analysis

Purpose

Stomach adenocarcinoma (STAD) is one of the most frequently diagnosed cancer in the world with both high mortality and high metastatic capacity. Therefore, the present study aimed to investigate novel therapeutic targets and prognostic biomarkers that can be used for STAD treatment.

Materials and Methods

We acquired four original gene chip profiles, namely GSE13911, GSE19826, GSE54129, and GSE65801 from the Gene Expression Omnibus (GEO). The datasets included a total of 114 STAD tissues and 110 adjacent normal tissues. The GEO2R online tool and Venn diagram software were used to discriminate differentially expressed genes (DEGs). Gene ontology (GO) and Kyoto Encyclopedia of Gene and Genome (KEGG) enriched pathways were also performed for annotation and visualization with DEGs. The STRING online database was used to identify the functional interactions of DEGs. Subsequently, we selected the most significant DEGs to construct the protein-protein interaction (PPI) network and to reveal the core genes involved. Finally, the Kaplan-Meier Plotter online database and Gene Expression Profiling Interactive Analysis (GEPIA) were used to analyze the prognostic information of the core DEGs.

Results

A total of 114 DEGs (35 upregulated and 79 downregulated) were identified, which were abnormally expressed in the GEO datasets. GO analysis demonstrated that the majority of the upregulated DEGs were significantly enriched in collagen trimer, cell adhesion, and identical protein binding. The downregulated DEGs were involved in extracellular space, digestion, and inward rectifier potassium channel activity. Signaling pathway analysis indicated that upregulated DEGs were mainly enriched in receptor interaction, whereas downregulated DEGs were involved in gastric acid secretion. A total of 80 DEGs were screened into the PPI network complex, and one of the most important modules with a high degree was detected. Furthermore, 10 core genes were identified, namely COL1A1, COL1A2, FN1, COL5A2, BGN, COL6A3, COL12A1, THBS2, CDH11, and SERPINH1. Finally, the results of the prognostic information further demonstrated that all 10 core genes exhibited significantly higher expression in STAD tissues compared with that noted in normal tissues.

Conclusion

The multiple molecular mechanisms of these novel core genes in STAD are worthy of further investigation and may reveal novel therapeutic targets and biomarkers for STAD treatment.

Do testosterone and sex hormone-binding globulin affect cancer risk? A Mendelian randomization and bioinformatics study

Using Mendelian Randomization (MR) and large-scale Genome-Wide Association Study (GWAS) data, this study aimed to investigate the potential causative relationship between testosterone and sex hormone-binding globulin (SHBG) levels and the onset of several cancers, including pathway enrichment analyses of single nucleotide polymorphisms (SNPs) associated with cancer allowed for a comprehensive bioinformatics approach, which offered a deeper biological understanding of these relationships. The results indicated that increased testosterone levels in women were associated with a higher risk of breast and cervical cancers but a lower risk of ovarian cancer. Conversely, increased testosterone was linked to lower stomach cancer risk for men, whereas high SHBG levels were related to decreased risks of breast and prostate cancers. The corresponding genes of the identified SNPs, as revealed by pathway enrichment analysis, were involved in significant metabolic and proliferative pathways. These findings emphasize the need for further research into the biological mechanisms behind these associations, paving the way for potential targeted interventions in preventing and treating these cancers.

Identification of key candidate genes and pathways in glioblastoma by integrated bioinformatical analysis

Glioblastoma (GBM), characterized by high morbidity and mortality, is one of the most common lethal diseases worldwide. To identify the molecular mechanisms that contribute to the development of GBM, three cohort profile datasets (GSE50161, GSE90598 and GSE104291) were integrated and thoroughly analyzed; these datasets included 57 GBM cases and 22 cases of normal brain tissue. The current study identified differentially expressed genes (DEGs), and analyzed potential candidate genes and pathways. Additionally, a DEGs-associated protein-protein interaction (PPI) network was established for further investigation. Then, the hub genes associated with prognosis were identified using a Kaplan-Meier analysis based on The Cancer Genome Atlas database. Firstly, the current study identified 378 consistent DEGs (240 upregulated and 138 downregulated). Secondly, a cluster analysis of the DEGs was performed based on functions of the DEGs and signaling pathways were analyzed using the enrichment analysis tool on DAVID. Thirdly, 245 DEGs were identified using PPI network analysis. Among them, two co-expression modules comprising of 30 and 27 genes, respectively, and 35 hub genes were identified using Cytoscape MCODE. Finally, Kaplan-Meier analysis of the hub genes revealed that the increased expression of calcium-binding protein 1 (CABP1) was negatively associated with relapse-free survival. To summarize, all enriched Gene Ontology terms and Kyoto Encyclopedia of Genes and Genomes pathways may participate in mechanisms underlying GBM occurrence and progression, however further studies are required. CABP1 may be a key gene associated with the biological process of GBM development and may be involved in a crucial mechanism of GBM progression.

The Dysregulation and Prognostic Analysis of STRIPAK Complex Across Cancers

The striatin-interacting phosphatase and kinase (STRIPAK) is the highly conserved complex, which gains increased attention in physiology and pathology process recently. However, limited studies reported the details of STRIPAK complex in cancers while some results strongly suggested it plays a vital role in tumorigenesis. Hence, we systematically analyzed the molecular and survival profiles of 18 STRIPAK genes to assess the value of STRIPAK complex across cancers. Our findings revealed the low frequencies of DNA aberrances and incomparable expression difference of STRIPAK genes between normal and tumor tissues, but they showed strong prognostic value in cancers, especially the liver hepatocellular carcinoma (LIHC) and kidney renal clear cell carcinoma (KIRC). Interestingly, STRIPAK genes were observed the opposite pattern of survival and expression in the above two cancer types. PPP2R1A and TRAF3IP3 were proposed as the oncogenic genes in LIHC and KIRC, respectively. The STRIPAK genes serve as oncogenes may due to the methylation heterogeneity. Taken together, our comprehensive molecular analysis of STRIPAK complex provides resource to facilitate the understanding of mechanism and utilize the potential therapies to tumors.