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.     

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.     

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.     

Bioinformatics Analysis of Microarray Data to Reveal Novel Genes Related to Cold-Resistance of Maize

Low temperature has become a major abiotic stress factor that can reduce maize yield and cause a number of economic loss. This study was designed to identify key genes and pathways associated with coldresistance of maize. The gene expression profile GSE46704, including 4 control temperature treated plants and 4 low temperature treated plants, was downloaded from the Gene Expression Omnibus database. Differentially-expressed genes (DEGs) were identified by limma package. Then, protein-protein interaction (PPI) network and module selection were constructed using Cytoscape. Moreover, the DEGs were re-matched based on the Zea mays L. gene ID and symbol data from PlantRegMap. Finally, the re-matched DEGs were performed functional and pathway enrichment analyses by the DAVID online tool. A total of 750 DEGs were screened (including 387 up-regulated and 363 down-regulated genes) In the PPI network, GRMZM2G070837_P01 and GRMZM2G114578_P01 had higher degrees. Besides, carbohydrate metabolic process, starch and sucrose metabolism and biosynthesis of secondary metabolites were significantly enriched in functional and pathway enrichment analysis. GRMZM2G070837_P01 and GRMZM2G114578_P01 might play a critical role in cold-resistance of maize. Meanwhile, carbohydrate metabolic process, starch and sucrose metabolism and biosynthesis of secondary metabolites might function in cold-resistance of maize.     

基于转录组的肝豆状核变性调控网络的构建和分析

本研究旨在利用生物信息学方法构建经铜诱导的ATP7B基因敲除HepG2细胞系的转录调控网络。探讨关键转录因子在肝豆状核变性发生、发展中的潜在作用机制。收集公共基因表达数据库(gene expression omnibus, GEO)中包含野生型、ATP7B基因敲除型、铜诱导的野生型和铜诱导的ATP7B基因敲除型HepG2细胞系数据。筛选由铜诱导产生的差异表达基因(differentially expressed genes,DEGs)后进行基因本体论(gene ontology,GO)、京都基因和基因组百科全书(Kyoto encyclopedia of genes and genomes, KEGG)富集分析。基于蛋白相互作用网络,识别疾病关键基因和功能模块,并对关键功能模块中的基因进行富集分析。最后,构建转录调控网络,筛选核心转录因子。共筛选出1 034个差异表达基因,其中上调525个,下调509个。上、下调关键功能模块分别包括了3785个和3931个基因。关键功能模块中的基因主要定位于细胞-基质连接、染色体、剪接复合体、核糖体等区域,共同参与了mRNA加工、组蛋白修饰、RNA剪切...     

Interaction network analysis revealed biomarkers in myocardial infarction

Myocardial infarction (MI) is a serious heart disease. The cardiac cells of patients with MI will die due to lack of blood for a long time. In this study, we aimed to find new targets for MI diagnosis and therapy. We downloaded GSE22229 including 12 blood samples from healthy persons and GSE29111 from Gene Expression Omnibus including 36 blood samples from MI patients. Then we identified differentially expressed genes (DEGs) in patients with MI compared to normal controls with p value < 0.05 and |logFC| > 1. Furthermore, interaction network and sub-network of these of these DEGs were constructed by NetBox. Linker genes were screened in the Global Network database. The degree of linker genes were calculated by igraph package in R language. Gene ontology and kyoto encyclopedia of genes and genomes pathway analysis were performed for DEGs and network modules. A total of 246 DEGs were identified in MI, which were enriched in the immune response. In the interaction network, LCK, CD247, CD3D, FYN, HLA-DRA, IL2, CD8A CD3E, CD4, CD3G had high degree, among which CD3E, CD4, CD3G were DEGs while others were linker genes screened from Global Network database. Genes in the sub-network were also enriched in the immune response pathway. The genes with high degree may be biomarkers for MI diagnosis and therapy.     

Discovery of core gene families associated with liver metastasis in colorectal cancer and regulatory roles in tumor cell immune infiltration

In this study, we aimed to uncover genes that drive the pathogenesis of liver metastasis in colorectal cancer (CRC), and identify effective genes that could serve as potential therapeutic targets for treating with colorectal liver metastasis patients based on two GEO datasets. Several bioinformatics approaches were implemented. First, differential expression analysis screened out key differentially expressed genes (DEGs) across the two GEO datasets. Based on gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses, we identified the enrichment functions and pathways of the DEGs that were associated with liver metastasis in CRC. Second, immune infiltration analysis identified key immune signature gene sets associated with CRC liver metastasis, among which two key immune gene families (CD and CCL) identified as key DEGs were filtered by protein-protein interaction (PPI) network. Some of the members in these gene families were associated with disease free survival (DFS) or overall survival (OS) in two subtypes of CRC, namely COAD and READ. Finally, functional enrichment analysis of the two gene families and their neighboring genes revealed that they were closely associated with cytokine, leukocyte proliferation and chemotaxis. These results are valuable in comprehending the pathogenesis of liver metastasis in CRC, and are of seminal importance in understanding the role of immune tumor infiltration in CRC. Our study also identified potentially effective therapeutic targets for liver metastasis in CRC including CCL20, CCL24 and CD70.     

IgA 肾病患者与膜性肾病患者外周血单核细胞mRNA分析

为了寻找诊断、鉴别IgA肾病(IgAN)和膜性肾病(MN)的血液特异性标记物,利用公共数据库中的IgAN和MN患者的外周血单核细胞(PBMCs)的转录组表达谱数据集识别特异性生物标记物,为诊断和鉴别提供简便、可靠的依据补充。从公共基因表达数据库(GEO)下载IgAN患者组(n=15)和MN患者组(n=8)芯片数据集,筛选前250个差异表达基因(DEGs)。通过分析筛选关键基因和途径,进行基因本体(GO)富集分析、京都基因与基因组百科全书(KEGG)通路分析和蛋白质与蛋白质相互作用关系(PPI)分析等进一步了解DEGs。通过分析共发现75个显著DEGs,其中73个上调基因,2个下调基因。GO富集分析的生物学过程(BP)主要包括蛋白质转运、内溶酶体到溶酶体转运、趋化因子介导的信号通路作用等。显著富集差异表达基因KEGG通路分析包括Endocytosis和Hepatitis B的相关信号通路。PPI筛选出EPS15、STAT4、CCL2、SUN2、SEC24C、SEC31A、GOLGB1、F2R,RAB12和PTK2B等关键基因。成功筛选出核心差异表达基因,为IgAN和MN的诊断和鉴别提供简便、可靠的依据补充,甚至提供治疗的新靶点。     

Identification of novel biomarkers and candidate small molecule drugs in rheumatoid arthritis and osteoarthritis based on bioinformatics analysis of high‐throughput data

Background: Rheumatoid arthritis (RA) and osteoarthritis (OA) are two major types of joint diseases. The present study aimed to identify hub genes involved in the pathogenesis and further explore the potential treatment targets of RA and OA.Methods: The gene expression profile of {"type":"entrez-geo","attrs":{"text":"GSE12021","term_id":"12021"}}GSE12021 was downloaded from Gene Expression Omnibus (GEO). Total 31 samples (12 RA, 10 OA and 9 NC samples) were used. The differentially expressed genes (DEGs) in RA versus NC, OA versus NC and RA versus OA groups were screened using limma package. We also verified the DEGs in {"type":"entrez-geo","attrs":{"text":"GSE55235","term_id":"55235"}}GSE55235 and {"type":"entrez-geo","attrs":{"text":"GSE100786","term_id":"100786"}}GSE100786. Functional annotation and protein–protein interaction (PPI) network construction of OA‐ and RA‐specific DEGs were performed. Finally, the candidate small molecules as potential drugs to treat RA and OA were predicted in CMap database.Results: 165 up-regulated and 163 down-regulated DEGs between RA and NC samples, 73 up-regulated and 293 down-regulated DEGs between OA and NC samples, 92 up-regulated and 98 down-regulated DEGs between RA and OA samples were identified. Immune response and TNF signaling pathway were significantly enriched pathways for RA‐ and OA‐specific DEGs, respectively. The hub genes were mainly associated with ‘Primary immunodeficiency’ (RA vs. NC group), ‘Ribosome’ (OA vs. NC group), and ‘Chemokine signaling pathway’ (RA vs. OA group). Arecoline and Cefamandole were the most promising small molecule to reverse the RA and OA gene expression.Conclusion: Our findings suggest new insights into the underlying pathogenesis of RA and OA, which may improve the diagnosis and treatment of these intractable chronic diseases.     

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.     

Prognostic value of tumour microenvironment-related genes by TCGA database in rectal cancer

Rectal cancer is a common malignant tumour and the progression is highly affected by the tumour microenvironment (TME). This study intended to assess the relationship between TME and prognosis, and explore prognostic genes of rectal cancer. The gene expression profile of rectal cancer was obtained from TCGA and immune/stromal scores were calculated by Estimation of Stromal and Immune cells in Malignant Tumors using Expression data (ESTIMATE) algorithm. The correlation between immune/stromal scores and survival time as well as clinical characteristics were evaluated. Differentially expressed genes (DEGs) were identified according to the stromal/immune scores, and the functional enrichment analyses were conducted to explore functions and pathways of DEGs. The survival analyses were conducted to clarify the DEGs with prognostic value, and the protein-protein interaction (PPI) network was performed to explore the interrelation of prognostic DEGs. Finally, we validated prognostic DEGs using data from the Gene Expression Omnibus (GEO) database by PrognoScan, and we verified these genes at the protein levels using the Human Protein Atlas (HPA) databases. We downloaded gene expression profiles of 83 rectal cancer patients from The Cancer Genome Atlas (TCGA) database. The Kaplan-Meier plot demonstrated that low-immune score was associated with worse clinical outcome (P = .034), metastasis (M1 vs. M0, P = .031) and lymphatic invasion (+ vs. -, P < .001). A total of 540 genes were screened as DEGs with 539 up-regulated genes and 1 down-regulated gene. In addition, 60 DEGs were identified associated with overall survival. Functional enrichment analyses and PPI networks showed that the DEGs are mainly participated in immune process, and cytokine-cytokine receptor interaction. Finally, 19 prognostic genes were verified by GSE17536 and GSE17537 from GEO, and five genes (ADAM23, ARHGAP20, ICOS, IRF4, MMRN1) were significantly different in tumour tissues compared with normal tissues at the protein level. In summary, our study demonstrated the associations between TME and prognosis as well as clinical characteristics of rectal cancer. Moreover, we explored and verified microenvironment-related genes, which may be the potential key prognostic genes of rectal cancer. Further clinical samples and functional studies are needed to validate this finding.