Identifying candidate genes for discrimination of ulcerative colitis and Crohn’s disease

In this study we aimed to screen effective biomarkers for differential diagnosis of ulcerative colitis (UC) and Crohn’s disease (CD). By using the gene expression profile dataset GSE24287 including 47 ileal CD, 27 UC and 25 non-inflammatory bowel diseases control downloaded from Gene Expression Omnibus database, we identified the differentially expressed genes (DEGs) between UC patients and controls as well as between CD patients and controls (|log₂FC(fold change)| > 1 and p < 0.05). Then Gene Ontology (GO) functional enrichment analyses were performed for these DEGs in two groups, followed by the construction of weight PPI (protein–protein interaction) networks. Subnets enriched for the PPIs and differentially expressed genes were constructed based on the weight PPI networks. The overlapping genes between the genes in the top 10 subnets with smallest p value and the DEGs were selected as the candidate genes of disease. A total of 75 DEGs were identified in UC group and 87 ones in CD group. There were 69 and 57 specific DEGs in CD group and UC group, respectively. The DEGs in CD group were mainly enriched in “inflammatory response” and “defense response”, while the most significantly enriched GO terms in UC group were “anion transport” and “chemotaxis”. FOS and SOCS3 were identified as candidate genes for CD and other three genes HELB, ZBTB16 and FAM107A were candidate genes for UC. In conclusion, there were distinct genetic alterations between UC and CD. The candidate genes identified in current study may be used as biomarkers for differential diagnosis of CD and UC.     

Bioinformatics Analyses of the Role of Vascular Endothelial Growth Factor in Patients with Non-Small Cell Lung Cancer

Purpose

This study was aimed to identify the expression pattern of vascular endothelial growth factor (VEGF) in non-small cell lung cancer (NSCLC) and to explore its potential correlation with the progression of NSCLC.

Methods

Gene expression profile GSE39345 was downloaded from the Gene Expression Omnibus database. Twenty healthy controls and 32 NSCLC samples before chemotherapy were analyzed to identify the differentially expressed genes (DEGs). Then pathway enrichment analysis of the DEGs was performed and protein-protein interaction networks were constructed. Particularly, VEGF genes and the VEGF signaling pathway were analyzed. The sub-network was constructed followed by functional enrichment analysis.

Results

Total 1666 up-regulated and 1542 down-regulated DEGs were identified. The down-regulated DEGs were mainly enriched in the pathways associated with cancer. VEGFA and VEGFB were found to be the initiating factor of VEGF signaling pathway. In addition, in the epidermal growth factor receptor (EGFR), VEGFA and VEGFB associated sub-network, kinase insert domain receptor (KDR), fibronectin 1 (FN1), transforming growth factor beta induced (TGFBI) and proliferating cell nuclear antigen (PCNA) were found to interact with at least two of the three hub genes. The DEGs in this sub-network were mainly enriched in Gene Ontology terms related to cell proliferation.

Conclusion

EGFR, KDR, FN1, TGFBI and PCNA may interact with VEGFA to play important roles in NSCLC tumorigenesis. These genes and corresponding proteins may have the potential to be used as the targets for either diagnosis or treatment of patients with NSCLC.     

Analysis of the potential pathways and target genes in spinal cord injury using bioinformatics methods

Spinal cord injury (SCI) remains to be the most devastating type of trauma for patients because of long lasting disability and limited response to the acute drug administration and efforts at rehabilitation. With the purpose to identify potential targets for SCI treatment and to gain more insights into the mechanisms of SCI, the microarray data of GSE2270, including 119 raphe magnus (RM) samples and 125 sensorimotor cortex (SMTC) samples, was downloaded from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) were screened in RM group and SMTC group compared with their corresponding controls, respectively. A protein–protein interaction (PPI) network was constructed based on the common DEGs identified in both RM group and SMTC group. Gene ontology (GO) and pathway enrichment analyses of the overlapping DEGs were performed. Furthermore, the common DEGs enriched in each pathway were analyzed to identify significant regulatory elements. Totally, 173 overlapping DEGs (130 up-regulated and 43 down-regulated) were identified in both RM and SMTC samples. These overlapping DEGs were enriched in different GO terms. Pathway enrichment analysis revealed that DEGs were mainly related to inflammation and immunity. CD68 molecule (CD68) was a hub protein in the PPI network. Moreover, the regulatory network showed that ras-related C3 botulinum toxin substrate 2 (RAC2), CD44 molecule (CD44), and actin related protein 2/3 complex (ARPC1B) were hub genes. RAC2, CD44, and ARPC1B may be significantly involved in the pathogenesis of SCI by participating significant pathways such as extracellular matrix-receptor signaling pathway and Toll-like receptor signaling pathway.     

Identification of differentially expressed genes and biological characteristics of colorectal cancer by integrated bioinformatics analysis

Colorectal cancer (CRC) ranks as one of the most common malignant tumors worldwide. Its mortality rate has remained high in recent years. Therefore, the aim of this study was to identify significant differentially expressed genes (DEGs) involved in its pathogenesis, which may be used as novel biomarkers or potential therapeutic targets for CRC. The gene expression profiles of GSE21510, GSE32323, GSE89076, and GSE113513 were downloaded from the Gene Expression Omnibus (GEO) database. After screening DEGs in each GEO data set, we further used the robust rank aggregation method to identify 494 significant DEGs including 212 upregulated and 282 downregulated genes. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed by DAVID and the KOBAS online database, respectively. These DEGs were shown to be significantly enriched in different cancer-related functions and pathways. Then, the STRING database was used to construct the protein–protein interaction network. The module analysis was performed by the MCODE plug-in of Cytoscape based on the whole network. We finally filtered out seven hub genes by the cytoHubba plug-in, including PPBP, CCL28, CXCL12, INSL5, CXCL3, CXCL10, and CXCL11. The expression validation and survival analysis of these hub genes were analyzed based on The Cancer Genome Atlas database. In conclusion, the robust DEGs associated with the carcinogenesis of CRC were screened through the GEO database, and integrated bioinformatics analysis was conducted. Our study provides reliable molecular biomarkers for screening and diagnosis, prognosis as well as novel therapeutic targets for CRC.     

Integrated microRNA–mRNA analysis of coronary artery disease

Although patients with coronary artery disease (CAD) have a high mortality rate, the pathogenesis of CAD is still poorly understood. The purpose of this study was to explore the underlying molecular mechanisms and potential target molecules for CAD. The platelet miRNA (GSE28858) and blood mRNA (GSE42148) expression profiles of patients with CAD and healthy controls were downloaded from Gene Expression Omnibus. Differentially expressed miRNAs and genes (DEGs) were identified by significant analysis of microarray algorithm after data preprocessing. Furthermore, the miRNA-target gene regulatory network was constructed based on miRecords database. The spearman correlation coefficients (ρ) between miRNAs and their target genes were calculated. Six up- (miR-340, miR-545, miR-451, miR454-5p, miR-624 and miR-585) and four down-regulated (miR-199a, miR-17-3p, miR-154 and miR-339) miRNAs were screened. Total 295 target genes of miR-545, miR-451, miR-585 and miR-154 were predicted. Among these 295 target genes, 7 genes were DEGs. Further analysis showed miR-545-TFEC and miR-585-SPOCK1 were highly positively correlated (ρ = 0.808091264; ρ = 0.874680776) in CAD samples. Therefore, differentially expressed miRNAs might participate in the pathogenesis of CAD by regulating their target genes.     

Classifying Integrated Signature Molecules in Macrophages of Rheumatoid Arthritis,Osteoarthritis, and Periodontal Disease: An Omics-Based Study

Rheumatoid arthritis (RA), osteoarthritis (OA), and periodontal disease (PD) are chronic inflammatory diseases that are globally prevalent, and pose a public health concern. The search for a potential mechanism linking PD to RA and OA continues, as it could play a significant role in disease prevention and treatment. Recent studies have linked RA, OA, and PD to Porphyromonas gingivalis (PG), a periodontal bacterium, through a similar dysregulation in an inflammatory mechanism. This study aimed to identify potential gene signatures that could assist in early diagnosis as well as gain insight into the molecular mechanisms of these diseases. The expression data sets with the series IDs {"type":"entrez-geo","attrs":{"text":"GSE97779","term_id":"97779"}}GSE97779, {"type":"entrez-geo","attrs":{"text":"GSE123492","term_id":"123492"}}GSE123492, and {"type":"entrez-geo","attrs":{"text":"GSE24897","term_id":"24897"}}GSE24897 for macrophages of RA, OA synovium, and PG stimulated macrophages (PG-SM), respectively, were retrieved and screened for differentially expressed genes (DEGs). The 72 common DEGs among RA, OA, and PG-SM were further subjected to gene–gene correlation analysis. A GeneMANIA interaction network of the 47 highly correlated DEGs comprises 53 nodes and 271 edges. Network centrality analysis identified 15 hub genes, 6 of which are DEGs (API5, ATE1, CCNG1, EHD1, RIN2, and STK39). Additionally, two significantly up-regulated non-hub genes (IER3 and RGS16) showed interactions with hub genes. Functional enrichment analysis of the genes showed that “apoptotic regulation” and “inflammasomes” were among the major pathways. These eight genes can serve as important signatures/targets, and provide new insights into the molecular mechanism of PG-induced RA, OA, and PD.     

Uncovering the pathogenesis and identifying novel targets of pancreatic cancer using bioinformatics approach

Pancreatic cancer is a uniformly lethal disease that can be difficult to diagnose at its early stage. Thus, our present study aimed to explore the underlying mechanism and identify new targets for this disease. The data GSE16515, including 36 tumor and 16 normal samples were available from Gene Expression Omnibus. Differentially expressed genes (DEGs) were screened out using Robust Multichip Averaging and LIMMA package. Moreover, gene ontology and pathway enrichment analyses were performed to DEGs. Followed with protein–protein interaction (PPI) network construction by STRING and Cytoscape, module analysis was conducted using ClusterONE. Finally, based on PubMed, text mining about these DEGs was carried out. Total 274 up-regulated and 93 down-regulated genes were identified as the common DEGs and these genes were discovered significantly enriched in cell adhesion and extracellular region terms, as well as ECM-receptor interaction pathway. In addition, five modules were screened out from the up-regulated PPI network with none in down-regulated network. Finally, the up-regulated genes, including MIA, MET and CEACAMS, and down-regulated genes, such as FGF, INS and LAPP, had the most references in text mining analysis. Our findings demonstrate that the up- and down-regulated genes play important roles in pancreatic cancer development and might be new targets for the therapy.     

Bioinformatics identification of the methylerythritol phosphate pathway associated genes in <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis> with <Emphasis Type="Italic">ceh1</Emphasis> mutant

The methylerythritol phosphate (MEP) pathway for the production of isoprenoids is recently discovered. The current study aimed to identify MEP pathway disorder-related molecular mechanisms and potential genes in Arabidopsis thaliana. Microarray data (GSE61675) obtained from ceh1 mutant plants and corresponding parental lines were retrieved from Gene Expression Omnibus (GEO) database and were applied for differentially expressed genes (DEGs) screening. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of DEGs were performed. Protein-protein interaction (PPI) network was then constructed and displayed by Cytoscape software. Total 762 DEGs including 620 up-regulated and 142 down-regulated genes were screened. In addition, a great many of DEGs were mainly involved in biosynthesis and metabolism-related pathways, such as stilbenoid, diarylheptanoid, and gingerol biosynthesis, and biosynthesis of terpenoids and steroids. Moreover, a PPI network contained 90 down-regulated genes and 497 up-regulated genes were obtained. Up-regulated DEGs including glutaredoxin (GRX480, cytochrome BC1 synthase (BCS1, syntaxin of plants 121 (SYP121) and A. thaliana MAP kinase 11 (ATMPK11) with higher degree in this network were hub nodes. Pathways including stilbenoid, diarylheptanoid, and gingerol biosynthesis obtained in our study were consistent with previous studies. Importantly, GRX480, BCS1 and ATMPK11 could have close interactions with the MEP pathway and may play important roles in the biosynthesis of isoprenoids.     

Identification of key genes and pathways of thyroid cancer by integrated bioinformatics analysis

Thyroid cancer is a common endocrine malignancy with a rapidly increasing incidence worldwide. Although its mortality is steady or declining because of earlier diagnoses, its survival rate varies because of different tumour types. Thus, the aim of this study was to identify key biomarkers and novel therapeutic targets in thyroid cancer. The expression profiles of GSE3467, GSE5364, GSE29265 and GSE53157 were downloaded from the Gene Expression Omnibus database, which included a total of 97 thyroid cancer and 48 normal samples. After screening significant differentially expressed genes (DEGs) in each data set, we used the robust rank aggregation method to identify 358 robust DEGs, including 135 upregulated and 224 downregulated genes, in four datasets. Gene Ontology and Kyoto Encyclopaedia of Genes and Genomes pathway enrichment analyses of DEGs were performed by DAVID and the KOBAS online database, respectively. The results showed that these DEGs were significantly enriched in various cancer-related functions and pathways. Then, the STRING database was used to construct the protein–protein interaction network, and modules analysis was performed. Finally, we filtered out five hub genes, including LPAR5, NMU, FN1, NPY1R, and CXCL12, from the whole network. Expression validation and survival analysis of these hub genes based on the The Cancer Genome Atlas database suggested the robustness of the above results. In conclusion, these results provided novel and reliable biomarkers for thyroid cancer, which will be useful for further clinical applications in thyroid cancer diagnosis, prognosis and targeted therapy.     

Comprehensive gene expression analysis reveals multiple signal pathways associated with prostate cancer

Prostate cancer (PC) depends on androgenic signaling for growth and survival. To data, the exact molecular mechanism of hormone controlling proliferation and tumorigenesis in the PC remains unclear. Therefore, in this study, we explored the differentially expressed genes (DEGs) and identified featured genes related to hormone stimulus from PC. Two sets of gene expression data, including PC and normal control sample, were downloaded from Gene Expression Omnibus (GEO) database. The t-test was used to identify DEGs between PC and controls. Gene ontology (GO) functional annotation was applied to analyze the function of DEGs and screen hormone-related DEGs. Then these hormone-related DEGs were further analyzed in constructed cancer network and Human Protein Reference Database to screen important signaling pathways they participated in. A total of 912 DEGs were obtained which included 326 up-regulated genes and 586 down-regulated genes. GO functional enrichment analysis identified 50 hormone-related DEGs associated with PC. After pathway and PPI network analysis, we found these hormone-related DEGs participated in several important signaling pathways including TGF-β (TGFB2, TGFB3 and TGFBR2), MAPK (TGFB2, TGFB3 and TGFBR2), insulin (PIK3R3, SHC1 and EIF4EBP1), and p53 signaling pathways (CCND2 and CDKN1A). In addition, a total of five hormone-related DEGs (SHC1, CAV1, RXRA, CDKN1A and SRF) were located in the center of PPI network and 12 hormone-related DEGs formed six protein modules. These important signal pathways and hormone-related DEGs may provide potential therapeutic targets for PC.     

Exploration of the molecular targets and mechanisms of suxiao xintong dropping pills for myocardial infarction by network pharmacology method

Background: Suxiao Xintong dropping pills (SXXTDP), a traditional Chinese medicine, is widely applied for treating myocardial infarction (MI). However, its therapy mechanisms are still unclear. Therefore, this research is designed to explore the molecular mechanisms of SXXTDP in treating MI.Methods: The active ingredients of SXXTDP and their corresponding genes of the active ingredients were retrieved from the Traditional Chinese Medicine Systems Pharmacology (TCMSP) database. MI-related genes were identified via analyzing the expression profiling data (accession number: {"type":"entrez-geo","attrs":{"text":"GSE97320","term_id":"97320"}}GSE97320). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were performed to study the shared genes of drug and disease. Through protein–protein interaction (PPI) network and the Cytoscape plugin cytoHubba, the hub genes were screened out. The compounds and hub targets binding were simulated through molecular docking method.Results: We obtained 21 active compounds and 253 corresponding target genes from TCMSP database. 1833 MI-related genes were identified according to P<0.05 and |log₂FC| ≥ 0.5. 27 overlapping genes between drug and disease were acquired. GO analysis indicated that overlapping genes were mainly enriched in MAP kinase activity and antioxidant activity. KEGG analysis indicated that overlapping genes were mainly enriched in IL-17 signaling pathway and TNF signaling pathway. We obtained 10 hub genes via cytoHubba plugin. Six of the 10 hub genes, including PTGS2, MAPK14, MMP9, MAPK1, NFKBIA, and CASP8, were acted on molecular docking verification with their corresponding compounds of SXXTDP.Conclusion: SXXTDP may exert cardioprotection effect through regulating multiple targets and multiple pathways in MI.     

Identification of hub genes and pathways in adrenocortical carcinoma by integrated bioinformatic analysis

Adrenocortical carcinoma (ACC), a rare malignant neoplasm originating from adrenal cortical cells, has high malignancy and few treatments. Therefore, it is necessary to explore the molecular mechanism of tumorigenesis, screen and verify potential biomarkers, which will provide new clues for the treatment and diagnosis of ACC. In this paper, three gene expression profiles (GSE10927, GSE12368 and GSE90713) were downloaded from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) were obtained using the Limma package. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were enriched by DAVID. Protein‐protein interaction (PPI) network was evaluated by STRING database, and PPI network was constructed by Cytoscape. Finally, GEPIA was used to validate hub genes’ expression. Compared with normal adrenal tissues, 74 up‐regulated DEGs and 126 down‐regulated DEGs were found in ACC samples; GO analysis showed that up‐regulated DEGs were enriched in organelle fission, nuclear division, spindle, et al, while down‐regulated DEGs were enriched in angiogenesis, proteinaceous extracellular matrix and growth factor activity; KEGG pathway analysis showed that up‐regulated DEGs were significantly enriched in cell cycle, cellular senescence and progesterone‐mediated oocyte maturation; Nine hub genes (CCNB1, CDK1, TOP2A, CCNA2, CDKN3, MAD2L1, RACGAP1, BUB1 and CCNB2) were identified by PPI network; ACC patients with high expression of 9 hub genes were all associated with worse overall survival (OS). These hub genes and pathways might be involved in the tumorigenesis, which will offer the opportunities to develop the new therapeutic targets of ACC.     

Identification and validation of key genes associated with non-small-cell lung cancer

Non-small-cell lung cancer (NSCLC) is one of the main causes of death induced by cancer globally. However, the molecular aberrations in NSCLC patients remain unclearly. In the present study, four messenger RNA microarray datasets (GSE18842, GSE40275, GSE43458, and GSE102287) were downloaded from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) between NSCLC tissues and adjacent lung tissues were obtained from GEO2R and the overlapping DEGs were identified. Moreover, functional and pathway enrichment were performed by Funrich, while the protein–protein interaction (PPI) network construction were obtained from STRING and hub genes were visualized and identified by Cytoscape software. Furthermore, validation, overall survival (OS) and tumor staging analysis of selected hub genes were performed by GEPIA. A total of 367 DEGs (95 upregulated and 272 downregulated) were obtained through gene integration analysis. The PPI network consisted of 94 nodes and 1036 edges in the upregulated DEGs and 272 nodes and 464 edges in the downregulated DEGs, respectively. The PPI network identified 46 upregulated and 27 downregulated hub genes among the DEGs, and six (such as CENPE, NCAPH, MYH11, LRRK2, HSD17B6, and A2M) of that have not been identified to be associated with NSCLC so far. Moreover, the expression differences of the mentioned hub genes were consistent with that in lung adenocarcinoma and lung squamous cell carcinoma in the TCGA database. Further analysis showed that all the six hub genes were associated with tumor staging except MYH11, while only the upregulated DEG CENPE was associated with the worse OS of patients with NSCLC. In conclusion, the current study showed that CENPE, NCAPH, MYH11, LRRK2, HSD17B6, and A2M might be the key genes contributed to tumorigenesis or tumor progression in NSCLC, further functional study is needed to explore the involved mechanisms.     

Exploring the Novel Computational Drug Target and Associated Key Pathways of Oral Cancer

Oral cancer (OC) is a serious health concern that has a high fatality rate. The oral cavity has seven kinds of OC, including the lip, tongue, and floor of the mouth, as well as the buccal, hard palate, alveolar, retromolar trigone, and soft palate. The goal of this study is to look into new biomarkers and important pathways that might be used as diagnostic biomarkers and therapeutic candidates in OC. The publicly available repository the Gene Expression Omnibus (GEO) was to the source for the collection of OC-related datasets. {"type":"entrez-geo","attrs":{"text":"GSE74530","term_id":"74530"}}GSE74530, {"type":"entrez-geo","attrs":{"text":"GSE23558","term_id":"23558"}}GSE23558, and {"type":"entrez-geo","attrs":{"text":"GSE3524","term_id":"3524"}}GSE3524 microarray datasets were collected for analysis. Minimum cut-off criteria of |log fold-change (FC)| > 1 and adjusted p < 0.05 were applied to calculate the upregulated and downregulated differential expression genes (DEGs) from the three datasets. After that only common DEGs in all three datasets were collected to apply further analysis. Gene ontology (GO) and pathway analysis were implemented to explore the functional behaviors of DEGs. Then protein–protein interaction (PPI) networks were built to identify the most active genes, and a clustering algorithm was also implemented to identify complex parts of PPI. TF-miRNA networks were also constructed to study OC-associated DEGs in-depth. Finally, top gene performers from PPI networks were used to apply drug signature analysis. After applying filtration and cut-off criteria, 2508, 3377, and 670 DEGs were found for {"type":"entrez-geo","attrs":{"text":"GSE74530","term_id":"74530"}}GSE74530, {"type":"entrez-geo","attrs":{"text":"GSE23558","term_id":"23558"}}GSE23558, and {"type":"entrez-geo","attrs":{"text":"GSE3524","term_id":"3524"}}GSE3524 respectively, and 166 common DEGs were found in every dataset. The GO annotation remarks that most of the DEGs were associated with the terms of type I interferon signaling pathway. The pathways of KEGG reported that the common DEGs are related to the cell cycle and influenza A. The PPI network holds 88 nodes and 492 edges, and CDC6 had the highest number of connections. Four clusters were identified from the PPI. Drug signatures doxorubicin and resveratrol showed high significance according to the hub genes. We anticipate that our bioinformatics research will aid in the definition of OC pathophysiology and the development of new therapies for OC.     

Identification of potential genes related to breast cancer brain metastasis in breast cancer patients

Brain metastases (BMs) usually develop in breast cancer (BC) patients. Thus, the molecular mechanisms of breast cancer brain metastasis (BCBM) are of great importance in designing therapeutic strategies to treat or prevent BCBM. The present study attempted to identify novel diagnostic and prognostic biomarkers of BCBM. Two datasets ({"type":"entrez-geo","attrs":{"text":"GSE125989","term_id":"125989"}}GSE125989 and {"type":"entrez-geo","attrs":{"text":"GSE100534","term_id":"100534"}}GSE100534) were obtained from the Gene Expression Omnibus (GEO) database to find differentially expressed genes (DEGs) in cases of BC with and without brain metastasis (BM). A total of 146 overlapping DEGs, including 103 up-regulated and 43 down-regulated genes, were identified. Functional enrichment analysis showed that these DEGs were mainly enriched for functions including extracellular matrix (ECM) organization and collagen catabolic fibril organization. Using protein–protein interaction (PPI) and principal component analysis (PCA) analysis, we identified ten key genes, including LAMA4, COL1A1, COL5A2, COL3A1, COL4A1, COL5A1, COL5A3, COL6A3, COL6A2, and COL6A1. Additionally, COL5A1, COL4A1, COL1A1, COL6A1, COL6A2, and COL6A3 were significantly associated with the overall survival of BC patients. Furthermore, COL6A3, COL5A1, and COL4A1 were potentially correlated with BCBM in human epidermal growth factor 2 (HER2) expression. Additionally, the miR-29 family might participate in the process of metastasis by modulating the cancer microenvironment. Based on datasets in the GEO database, several DEGs have been identified as playing potentially important roles in BCBM in BC patients.